FireFly (PC GAMESS) - ошибка SCF DID NOT CONVERGE...NO MPLE

[BIOHAZARD]

Собственно есть инпут:

Код: Выделить всё

! Minimize (Energy/Geometry) MP2/6-311G(d,p)
 $CONTRL COORD=CART MAXIT=500 MPLEVL=2
MULT=1 RUNTYP=OPTIMIZE SCFTYP=RHF UNITS=ANGS
 $END
 $BASIS GBASIS=N311 NDFUNC=1 NGAUSS=6 NPFUNC=1 POLAR=POPLE
 $END
 $SCF DAMP=.false. DEM=.false. DIIS=.false. DIRSCF=.true. 
EXTRAP=.true. RSTRCT=.false. SHIFT=.false. SOSCF=.true.
 $END
 $GUESS GUESS=HUCKEL $END
 $SYSTEM MWORDS=75 $END
 $DATA
...
 $END

в итоге расчёт заканчивается с ошибкой

...... END OF RHF CALCULATION ......

CPU TIME: STEP = 35745.11 , TOTAL = 35747.5 SECONDS ( 595.8 MIN)
WALL CLOCK TIME: STEP = 35840.90 , TOTAL = 35843.4 SECONDS ( 597.4 MIN)
CPU UTILIZATION: STEP = 99.73%, TOTAL = 99.73%
SCF DID NOT CONVERGE...NO MPLEVL=2 CALCULATION

CPU TIME: STEP = 0.00 , TOTAL = 35747.5 SECONDS ( 595.8 MIN)
WALL CLOCK TIME: STEP = 0.01 , TOTAL = 35843.4 SECONDS ( 597.4 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 99.73%

2705058 WORDS OF DYNAMIC MEMORY USED
EXECUTION OF FIREFLY TERMINATED ABNORMALLY AT 22:42:34 LT 8-SEP-2010

в квантово-химических расчётах я полный профан, так что не понимаю как с этим бороться
просто увеличить MAXIT или и при этом сходимости не будет? может поменять другие настройки?
на этом форуме нашёл пару топиков отдалённо похожих по смыслу, но ответа там нет

то же самое и касаемо pm3

Код: Выделить всё

! Minimize (Energy/Geometry) PM3
 $CONTRL COORD=CART ICHARG=0 MAXIT=5000
MULT=1 RUNTYP=OPTIMIZE SCFTYP=RHF UNITS=ANGS
 $END
 $BASIS GBASIS=PM3 $END
 $SCF DAMP=.false. DEM=.false. DIIS=.false. DIRSCF=.true.
EXTRAP=.true. RSTRCT=.false. SHIFT=.false. SOSCF=.true.
 $END
 $STATPT METHOD=QA NSTEP=50 OPTTOL=0.001 $END
 $FORCE TEMP=298.15  $END
 $GUESS GUESS=HUCKEL $END
 $SYSTEM MWORDS=50 $END
 $DATA
...
 $END

на выходе

...... END OF PROPERTY EVALUATION ......

CPU TIME: STEP = 0.00 , TOTAL = 77.4 SECONDS ( 1.3 MIN)
WALL CLOCK TIME: STEP = 0.00 , TOTAL = 88.2 SECONDS ( 1.5 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 87.76%
1 ***** FAILURE TO LOCATE STATIONARY POINT, SCF HAS NOT CONVERGED
UPDATED HESSIAN, GEOMETRY, AND VECTORS WILL BE PUNCHED FOR RESTART

CPU TIME: STEP = 0.00 , TOTAL = 77.4 SECONDS ( 1.3 MIN)
WALL CLOCK TIME: STEP = 0.00 , TOTAL = 88.2 SECONDS ( 1.5 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 87.76%

197280 WORDS OF DYNAMIC MEMORY USED
EXECUTION OF FIREFLY TERMINATED ABNORMALLY AT 5:29:59 LT 9-SEP-2010

во вложении полные инпут и выходной файлы для первого расчёта (mp2)

[o-oxhem]

$SCF DIIS=.T. DIRSCF=.true.
SHIFT=.T.
$END

И нужно добавить (модифицировать) $SYSTEM kdiag=0 TIMLIM=9999999 mwords=80 $END

Кстати, ошибку Вы диагностировали неверно. Правильно:
SCF IS UNCONVERGED, TOO LITTLE TIME LEFT
TIME TO FORM FOCK OPERATORS= 35670.7 SECONDS ( 205.0 SEC/ITER)
FOCK TIME ON FIRST ITERATION= 183.3, LAST ITERATION= 205.0
TIME TO SOLVE SCF EQUATIONS= 73.1 SECONDS ( 0.4 SEC/ITER)

FINAL ENERGY IS 0.0000000000 AFTER 174 ITERATIONS

[VTur]

Можно еще отсечку интегралов уменьшить icut=11

[BIOHAZARD]

$SCF DIIS=.T. DIRSCF=.true.
SHIFT=.T.
$END

И нужно добавить (модифицировать) $SYSTEM kdiag=0 TIMLIM=9999999 mwords=80 $END

Кстати, ошибку Вы диагностировали неверно. Правильно:
SCF IS UNCONVERGED, TOO LITTLE TIME LEFT
TIME TO FORM FOCK OPERATORS= 35670.7 SECONDS ( 205.0 SEC/ITER)
FOCK TIME ON FIRST ITERATION= 183.3, LAST ITERATION= 205.0
TIME TO SOLVE SCF EQUATIONS= 73.1 SECONDS ( 0.4 SEC/ITER)

FINAL ENERGY IS 0.0000000000 AFTER 174 ITERATIONS

нет, это не помогло
при увеличении числа циклов расчёт scf начинает просто блуждать по кругу, на выходе имеем out файл размером в 20 мб всё с той же ошибкой
в итоге сделал предварительную оптимизацию строения молекулы и дальше пошло считать нормально
видимо слишком уж кривая модель была, хотя gamess us её схавал на ура и пережевал при помощи pm3, результат этого расчёта я и загнал в зс gamess (он же firefly)

[o-oxhem]

Для достижения сходимости часто требуется перепробовать несколько разных вариантов ключевых слов.
И, может быть, стуктур. Последнее часто жутко экономит время. Намёнтанный глаз химика очень помогает в квантовой химии.

[BIOHAZARD]

Намёнтанный глаз химика очень помогает в квантовой химии.

в том то и дело, что на глаз там отличий никаких, но, видимо, на пару долей ангстрема атомы подвинулись и всё сразу встало на свои места