S1 оптимизация

[ruderal]

Ой, а что это Вы сделали? PURES=.f. в сочетании с WSTATE(1)=1,1 привело к тому, что Вы усредняли по основному состоянию S0 и первому -- нет, не возбужденному S1 -- а первому триплету, T1. И то, что Вы делали -- было оптимизацией состояния T1 (RUNTYP=optimize, ISTATE=2, а ключ method=conic в $statpt и группа $mcaver с ключами отсутствуют). Короче, Вы сделали совсем не то, что собирались. А где Вы образец инпута брали? поищите на сайте и на форуме FireFly, там точно есть нужные файлы.

Попытался найти коническое пересечение, но получается какой-то глюк(
Может быть, найдете ошибку в инпуте?

[ruderal]

поменял guga на aldet - сейчас считает

[sanya1024]

По идее, GUGA и ALDET должны приводить к одному и тому же результату. А что за глюк был?

[ruderal]

По идее, GUGA и ALDET должны приводить к одному и тому же результату. А что за глюк был?

DELTA ENERGY IS 574.118032405676
NSERCH= 0 ENERGY= -2009.4131134
...
MAXIMUM GRADIENT = 0.0000000 RMS GRADIENT = 0.0000000
FORCE CONSTANT MATRIX NOT UPDATED --- TAKING FIRST STEP
не знаю, что за глюк был...

с aldet геометрия в итоге развалилась, и расчет гессиана вылетел. аутпут прикрепил

[sanya1024]

***** MINIMUM ENERGY CROSSING POINT LOCATED *****
Структура выглядит нормальной для конического пересечения. "Разрыв" связи N1-C2 на самом деле никакой не разрыв, просто она удлинилась сильнее, чем ChemCraft считает допустимым для связи.
А частоты зачем считали? частоты имеют смысл для минимума или седла, т.е. точки, где ППЭ гладкая. А в конических пересечениях она запросто может не быть гладкой (оставаясь непрерывной) -- откуда там вторые производные? Вот они и не посчитались, и программа обругалась, отсюда и абнормальное завершение.

[ruderal]

***** MINIMUM ENERGY CROSSING POINT LOCATED *****
Структура выглядит нормальной для конического пересечения. "Разрыв" связи N1-C2 на самом деле никакой не разрыв, просто она удлинилась сильнее, чем ChemCraft считает допустимым для связи.

Большое спасибо за ответ.
Теперь было бы логично просканировать по N1-C2 связи S1 и S0, да?

[sanya1024]

Ну... можно сделать что-то вроде линейного транзита между минимумом S0 и найденным коническим пересечением -- просто ряд точек в SA-CASSCF с последующим XMCQDPT. Если делать релаксированный скан по длине связи N1-C2, то релаксировать надо S1, а за энергией S0 просто следить, какая она получается. Понятно, что в точности структуру основного состояния при таком скане не получите, но по крайней мере точка будет близкой (не вижу причин для обратного). Если же наоборот, релаксировать S0, то скорее всего от точки конического пересечения окажетесь далеко.