спин-орбитальное взаимодействие

[bezceller]

добрый день
мне нужно расчитать спин-орб. взаим. для дуплета в firefly, почитав мануал, я пришел к выводу, что это невозможно сделать в этой программе... так ли это и какие есть алтренативы (беслпатные)?

[sanya1024]

Альтернатива -- GAMESS-US. Читайте мануал (ищите по слову spin-orbit).

[bezceller]

извиняюсь, неточно задал вопрос... хотел посчитать неограниченным хф, так что теперь пересчитываю rohf

[sanya1024]

ХФ не советую -- энергии невозмущенных состояний получаются никуда не годными. Делают либо в DFT (имхо, тоже энергии не всегда получаются адекватными), либо в CASSCF с уточнением энергий по QDPT или даже (там, где ресурсы позволяют) SOCI (second order CI).

Гамессу в качестве входных данных требуются энергии взаимодействующих состояний (соответственно, они уже должны быть посчитаны более-менее прилично) и орбитали, тоже достаточно прилично описывающие взаимодействующие состояния. Очевидно, что ХФ энергии и орбитали тут не годятся.

[bezceller]

спасибо за ответ
что то я плохо в последнее время выражаю свои мысли (приболел, что ли... или другое что)

вобщем, я хотел посчитать uhf в dft при помощи кв (CITYP = GUGA)...
ведь, насколько я понимаю, различие между rhf и uhf в учёте спаривания электронов (rhf - только спаренные электроны с одинаковой попарно энергией, uhf - неспаренные электроны, а rohf - rhf "ядро" с 1м неспаренным элеткроном) ?

ладно, надо больше умных книжек читать...

[sanya1024]

вобщем, я хотел посчитать uhf в dft при помощи кв (CITYP = GUGA)...

????
Взять DFT-шные орбитали и сделать на них GUGA CI? имхо, сомнительный метод. Вдобавок это никак не будет аналогом UHF -- хотя бы потому, что в GUGA CI строятся конфигурации, соответствующие определенному спину. Если хотите получить все возможные значения спина (при данном S_z), тогда уж делайте ALDET CI. Кстати, в описанном выше способе CASSCF+QDPT тоже орбитали получают с CISTEP=ALDET.

ладно, надо больше умных книжек читать...

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

[sanya1024]

Кстати, полезная статья: http://wires.wiley.com/WileyCDA/WiresAr ... CMS83.html -- pdf доступен.

[bezceller]

спасибо