Firefly and PC GAMESS-related discussion club


 
Learn how to ask questions correctly  
 
 
We are NATO-free zone
 



Re: Problem with SCF after couple steps of geometry optimization

Solntsev Pasha
solntsev@univ.kiev.ua


OK, i played around and used D5 symmetry for the molecule (i didn't want). It looks like it works. I am trying to use opttol=1d-5 but have some fluctuation.

Pavel.



On Mon Apr 1 '13 6:11am, Solntsev Pasha wrote
---------------------------------------------
>Dear friends.

>I have a problem. I am trying to optimize geometry of ferrocene^(+1) molecule. I use cc-pVDZ basis set for all my calculations. Initially, i optimized geometry of ferrocene+hessian and then used equilibrium geometry to optimize the geometry of ferrocinium cation using UHF wave function MULT=2 icharg=1. I also used $statpt ststep=0.15  $end.

>Then i started new calculation with DFT level of theory BPW91 functional I used converged geometry and MO vectors from the HF job to start the DFT job. I was able to converge the first SCF only with SOSCF=.t. ,  DIIS=.t. didn't work. So, using soscf=.t. i started geometry optimization. However, after two steps of geometry optimization the SCF can't be converged. So far, i used default GDIIS method. But i am not sure what parameters i can use to tune gdiis algorithm.

>I used C1 symmetry for all my calculations. I also wasn't sure do we need to use d5=.t. for the spherical basis set and C1 point group?

>Relevant part of the input file presented below. Any advice?

> $CONTRL MAXIT=300 NZVAR=1
> inttyp=hondo icut=11 itol=30
> exetyp=run d5=.t. dfttyp=bpw91
> runtyp=optimize scftyp=uhf  
> icharg=1 mult=2 $END
> $SYSTEM TIMLIM=2736 mwords=20  $END
> $zmat auto=.t. dlc=.t. $end
> $SCF FDIFF=.FALSE. dirscf=.t. soscf=.t. diis=.f. nconv=7
> rstrct=.t.  $END
> $guess guess=moread norb=249 $end
> $basis gbasis=cc-pVDZ extfil=.t.  $end
> $statpt opttol=1d-6 nstep=200 hssend=.t. ststep=0.1 $end
> $DATA

>Pavel.
>


[ Previous ] [ Next ] [ Index ]           Tue Apr 2 '13 1:43am
[ Reply ] [ Edit ] [ Delete ]           This message read 694 times