Firefly and PC GAMESS-related discussion club



Learn how to ask questions correctly


Re^2: Convergence problems: Almost converged but oscillating

Miro Moman
edelmiro.moman@gmail.com


Now I am trying to MOREAD the orbitals from the converged Def2-TZVP run to be used with a larger basis set (aug-cc-pVTZ).

I receive the following error:

*** ERROR READING THE MO COEFFICIENTS AT   1  77
    THE INPUT VALUES WERE   2   1

Indeed,  1 77 does not exist :

1 75-2.45175653E-05 4.67451306E-07 8.18963875E-06-2.25179089E-05 1.05807521E-05
1 76 1.49291949E-05
2  1 7.09245510E-05 9.40896541E-06 5.16541291E-05-1.12328743E-04 1.24681158E-04

Do I need to pad with zeroes the $VEC section?


On Sat Jun 26 '10 11:35am, Miro Moman wrote
-------------------------------------------
>Ok. It has readily achieved convergence by switching to SOSCF. SOSCF appears to be more robust than DIIS in this particular case.
>
>
>On Sat Jun 26 '10 10:34am, Miro Moman wrote
>-------------------------------------------
>>Hello,

>>I am trying to optimize the geometry of a [Fe4S4]+ cluster.

>>In fact, the entire structure is [(Fe4S4)(C2H2)(SCH3)3]-2.

>>Having had problems with the full complex, I am now focusing on [Fe4S4]+

>>This is my current input file:

>> $CONTRL SCFTYP=UHF MULT=2 MAXIT=200 dfttyp=pbe0 d5=1 icharg=1 $END
>>!$CONTRL RUNTYP=OPTIMIZE $END
>> $CONTRL INTTYP=HONDO $END
>> $SYSTEM TIMLIM=99999 MWORDS=99 $END
>> $SCF DIIS=.T. $END
>> $statpt nstep=1000 $end
>> $scf dirscf=.t. fdiff=.f. $end
>> $guess guess=moread norb=376 $end
>> $smp httfix=.f. smppar=.t. $end

>>I have tried many different options in a more or less combinatorial way, such as : ICUT=11, ITOL=30, rstrct=.t., fdiff=.t., etc.

>>I have tried with larger and smaller basis sets (currently Def2-TZVP) and also I have used HF instead of DFT. In addition, thinking that the initial geometry could be the culprit, I have tried to optimize it.

>>Often it seems that the density has nearly converged, but it oscillates quite a lot, for instance:

>>  90 89 -6646.252407632     0.000030038   0.012915252   0.000146340     1661471091   9105267
>>  91 90 -6646.252436008    -0.000028377   0.004901678   0.000065611     1661477904   9105225
>>  92 91 -6646.252446605    -0.000010597   0.001334835   0.000051306     1661481388   9105176
>>  93 92 -6646.252449524    -0.000002919   0.002216061   0.000043654     1661486660   9105169
>>  94 93 -6646.252454263    -0.000004739   0.000922602   0.000045303     1661496334   9105149
>>  95 94 -6646.252456240    -0.000001977   0.002316391   0.000044423     1661501828   9105141
>>  96 95 -6646.252461196    -0.000004956   0.001103882   0.000043711     1661505726   9105136
>>  97 96 -6646.252463602    -0.000002406   0.000039273   0.000042713     1661508744   9105140
>>  98 97 -6646.252463639    -0.000000037   0.000755261   0.000041815     1661508666   9105145
>>  99 98 -6646.252462045     0.000001594   0.000914736   0.000042353     1661506603   9105146
>> 100 99 -6646.252463837    -0.000001792   0.004938129   0.000042491     1661508480   9105138

>>Any suggestions?


[ Previous ] [ Next ] [ Index ]           Sat Jun 26 '10 12:06pm
[ Reply ] [ Edit ] [ Delete ]           This message read 823 times