TS+TD-DFT in gaussian

[surius]

чтото на больную голову накатило посчитать TS в первом возбужденном. И тут "внезапно" встал такой вопрос: чтобы посчитать TS 
необходимо подсунуть в расчет силовые константы. И вот что интересно откуда их взять? Т.е. я так понимаю они должны были бы быть посчитаны 
на стадии оптимизации...
но у меня есть только оптимизированныя геометрия для минимума в S1. Или же из тащить за собой с самого S0. Я совсем запутался



заранее спасибо за помощь

[alxyppv]

Ну, если ОЧЕНЬ надо, то ведь можно и численно гессиан посчитать - в том смысле, что градиены аналитические у Вас есть (раз оптимизация есть), значит надо дергать поочередно атомы и считать градиенты... Если система не очень большая, о особых проблемм быть не должно. Если Вы считаете, что TS в возбужденном состоянии близко по строению к TS в основном (поскольку неясно что за процессс, то сложно сказать что-л. определенное) - то можно скормить и гессиан основного сотояния, но можно и сильно провраться...

[sanya1024]

Численные градиенты в TDDFT (да и просто в DFT) в Гамессе -- дело долгое, но в своем роде благодарное. Потому что если расчет оборвался, то остается *.irc файл, в к-рый пишется $VIB. Если этот $VIB приаттачить к инпуту (не забыв подписать в конец $END), то расчет автоматически рестартуется с последнего смещения. Так что рано или поздно доведете дело до конца.