gamess-us ошибка при оптимизации.

[rebel]

При оптимизации молекулы с 78 атомами в симметрии С1, Gamess выдает следующую ошибку:
INITIAL GUESS ORBITALS GENERATED BY HUCKEL ROUTINE.
HUCKEL GUESS REQUIRES 9690948 WORDS.
EQLRAT HAS FAILED TO CONVERGE FOR ROOT 863
0*** EVVRSP PARAMETERS ***
*** N = 931 ***
*** NVECT = 931 ***
*** LENA = 433846 ***
*** NV = 1140 ***
*** IORDER = 0 ***
*** IERR = 863 ***
UNABLE TO DIAGONALIZE IN COPROJ
----------------------------------------------------------
Как ее решить?)

[avelon]

либо измените метод расчета начальных орбиталей например на HCORE или попробуйте изменить метод диагонализации например KDIAG=2

[rebel]

А можно получить ключ с опциями для него?)(где это надо поставить)

[avelon]

$SYSTEM KDIAG=2 $END
$GUESS GUESS=HCORE $END
кстати правку по программе никто не отменял))

[rebel]

Cпасибо, вроде идет оптимизация на кластере) правда уже 5 часов на 16 узлах)

[avelon]

получилось посчитать?

[rebel]

К сожалению не оптимизировал в геометрии C1, поставил симметрию C2v (тетрафенилпорфин). Надеюсь пройдет.

[avelon]

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

[rebel]

Результат таков. В ADF оптимизация выполнена,в Gaussian тоже, но в Gamess-Us упорно она не проходит. Прикрепляю файлы с геометрией в декартовых координатах. В АДф был использован слетеровский базис - QZ4P и метод OPBE, в Gaussian B3LYP и базис 6-31G(d). Я понимаю, что геометрия получилась разной, так как программы разные и методы, базисы разные. Но не может же она так сильно отличаться? Молекула Тетрафенилпорфирин.
--------------------------------------------------------------------------------------------------------------------------------
Ключи использовал для Gamess такие как подсказали здесь. (метод диагонализации и Guess). Может в Gamess-Us нужно использовать другой метод оптимизации? я слышал, что по умолчанию задается метод Berny, но есть и другие. Или причина в другом?