MCSCF

[sanya1024]

Коллеги! поскольку мне надоело многократно объяснять, как делать MCSCF расчеты, а искать по форуму новым пользователям лень, я написала практическое руководство тут: http://photon190573.livejournal.com/tag/mcscf, и буду продолжать эту тему по мере возможности.

[OrganikNV]

Спасибо Вам огромное!!!!

[amge]

Присоединяюсь к благодарностям! Для меня это - самое то что доктор прописал .

[Yu/2]

Да, спасибо. Я там уже пару анонимных (звиняйте, аккаунта нема) комментов оставил.

[sanya1024]

А, это Вы

[iAlex]

del

[sanya1024]

Во-первых, не надо группы Oh. Когда имеете дело с возб. состояниями, от симметрии лучше отказаться, иначе она Вас заведет... куда-то, куда Вам совсем не надо. Берите группу C1 и задавайте все атомы в явном виде.
Дальше, чтобы Вам что-то подсказать (например, что (15,15) -- это неправильно ), я должна видеть Ваши стартовые орбитали. Тогда мы сможем решить, что оставить неактивным (NCORE), а что включить в акт. пространство (NACT и NELS). У Вас ион с замкнутой оболочкой -- но это октаэдр, и как там разлеглись орбитали центрального атома и лигандов -- навскидку не скажу, надо видеть хоть какую-нибудь выдачу (ес-но, с тем ECP, с к-рым Вы считаете, чтобы упрятать в него нужное кол-во электронов).

[iAlex]

del

[sanya1024]

С инпутом в чем проблема: я не могу посмотреть орбитали в $VEC в инпуте, скажем, Кемкрафтом. Я добавляю туда строчки exetyp=check и $GUESS GUESS=moread norb=<полное кол-во орбиталей> prtmo=.t. asis=.t. $end, получаю выдачу за 1 минуту на 1 ядре, заменяю в выдаче
----------------------
INITIAL GUESS ORBITALS
----------------------
на
------------
EIGENVECTORS
------------
и открываю Кемкрафтом. Правда, чтобы программа не ругалась на неправильно заданное акт. пространство, заменяю scftyp=mcscf на rhf -- нам же только стартовые орбитали посмотреть.
Вас интересует основное состояние или возбуждения? пытаться в CASSCF описывать основное состояние с замкнутой оболочкой -- это стрелять из пушки по воробьям. Метод не для этого. Улучшить описание осн. состояния с замкнутной оболочкой можно в MP2 и более высоких порядках ТВ, в CCSD. Поэтому будем считать, что нас интересуют возбужденные состояния.
На орбитальные энергии внимания не обращаем. Для начала можем включить в акт. пространство группы 1t1g и 3a1g. Поскольку симметрию мы отключили, надо просто не забывать, что разбивать пачки вырожденных орбиталей (как Вы сделали при первом удачном запуске с nact=2 nels=2) нельзя -- решение будет перекошенным. Поэтому включаем всю пачку вырожденных состояний целиком. Если решили для начала остановиться на 1t1g и 3a1g, то $det будет выглядеть так:
$DET NACT=4 NCORE=60 NELS=6 NSTATE=7 ITERMX=200 WSTATE(1)=1.,1.,1.,1. cvgtol=1d-7 $END
Что это означает? NACT=4, т.е. 3 штуки из 1t1g и 1 штука 3a1g. На 3-х 1t1g сидят 6 электронов, поэтому NELS=6. NCORE=60 -- это 63 занятых минус 3 штуки 1t1g активных. NSTATE берем с запасом. WSTATE(1)=1.,1.,1.,1. -- это основное (2220) и три варианта однократных возбуждений(2211, 2121, 1221). Это при последующих запусках можно будет и скорректировать (например, окажется, что однократные возбуждения тут по симметрии не проходят). Кстати, если нам не нужны триплеты, то добавим pures=.t. Наоборот, если нужны и триплеты, и синглеты, то ставим pures=.f., заказываем NSTATE=10 (с запасом) и WSTATE(1)=1.,1.,1.,1.,1.,1.,1. -- к синглетам конфигураций 2211, 2121, 1221 добавятся триплеты таких же конфигураций.
Обязательно добавьте $ciinp castrf=.t. $end, чтобы не запилить диск интегралами. Это и будет первым удачным запуском CASSCF расчета.

Дальше можно будет добавлять по одной оболочке (оболочки не разбиваем!). Например, оставляем в занятых 1t1g, а к незанятым добавляем 4t1u. Тогда NELS останется прежним, а NACT увеличится на 3. Если добавляете занятые -- соответственно уменьшайте NCORE.

[sanya1024]

По исходной ссылке -- апдейт. Последняя часть главы методички.