DFT Analytic Hessian in GAMESS ?

[timn]

Уважаемые форумчане!

Подскажите, пожалуйста, доступно ли в какой либо реализации GAMESS
аналитическое вычисление Гессиана при использовании DFT ?

На сколько я понял из документации на US GAMESS (GAMESS VERSION = 11 APR 2008 (R1)),
при указании, например, DFTTYP=B3LYP в секции $CONTRL автоматически
выбирается метод METHOD=SEMINUM для вычисления вторых производных от энергии...
Мало того, что это затратно по времени, так еще и дает ошибку в частотах валентных
колебаний порядка 30 см-1... Gaussian, кстати, умеет считать вторые производные
аналитически для всех DFT методов...

Спасибо!

[VTur]

Обещали скоро сделать в ФайерФлай.

если $FORCE PURIFY =.TRUE. VIBANL=.TRUE. NVIB=2 DECOMP =.t. $END , то разность с ГАУССИАНОМ не более +-10 см-1. Это для для гармонических частот.
А для гармонического приближения ошибка может быль любой.

[timn]

Спасибо!
Не подскажите, ФайерФлай для linux бывает?...

[VTur]

Конечно!
http://classic.chem.msu.su/gran/gamess/index.html

[OrganikNV]

Извините за возможно дурацкий вопрос, но я прошу объясните простыми словами, чем конкретно полезна функция NVIB в группе $FORCE при расчете Гессиана?

[sanya1024]

Все очень просто: когда делается численный расчет гессиана (в DFT -- всегда численный), каждый атом смещают из минимума по очереди вдоль осей x, y и z, в смещенной геометрии вычисляют градиент, и из этих градиентов вычисляют матрицу вторых производных энергии по координатам -- гессиан. Теоретически, достаточно смещать атомы только в одном направлении от минимума (скажем, вперед по каждой из xyz координат), но реально нужно делать два таких смещения (NVIB=2), вперед и назад -- для точности. Если этого не сделать, могут вылезти непонятные мнимые частоты.

[OrganikNV]

огромное спасибо за подробный ответ!!!