Gaussian - архивная секция, параметр PUHF

вопросы строения молекул и квантовой химии
Ответить
Аватара пользователя
EvgeniX
Сообщения: 2780
Зарегистрирован: Пт апр 27, 2007 5:32 am

Gaussian - архивная секция, параметр PUHF

Сообщение EvgeniX » Пн окт 09, 2017 10:15 pm

Доброй ночи!

Кто подскажет, что означают параметры "PUHF=", "PMP2-0=" в архивной секции log-файла Gaussian?

Код: Выделить всё

 1\1\ WCSS.PL-BEM-ATLANTACITY\SP\UMP2-FC\CC-pVTZ\C1H2(3)\EVGENIY\08-Oct
 -2017\0\\#p UMP2/cc-pVTZ sp\\13_CH2rr\\0,3\C,0,0.000000001,0.,0.059683
 536\H,0,-1.003543593,0.,-0.355320828\H,0,1.003543585,0.,-0.355320831\\
 Version=ES64L-G09RevE.01\State=3-B1\HF=-38.9373111\MP2=-39.055181\PUHF
 =-38.9402706\PMP2-0=-39.0571601\S2=2.017076\S2-1=2.005728\S2A=2.000077
 \RMSD=3.090e-09\PG=C02V [C2(C1),SGV(H2)]\\@

Аватара пользователя
Гесс
Сообщения: 13067
Зарегистрирован: Ср фев 15, 2012 11:19 pm

Re: Gaussian - архивная секция, параметр PUHF

Сообщение Гесс » Пн окт 09, 2017 11:44 pm

Посмотрел свои недавние аутпуты, все клоздшеллнутых систем, там такого нет, так что это явно фишка анрестриктнутых расчетов.
Гугл говорит:

Another approach is to run an unrestricted calculation then project out the spin contamination after the wave function has been obtained (PUHF, PMP2).

Projected UHF (PUHF) annihilates all spin contaminants from the self-consistent UHF wave function. The projected energy is evaluated as the expectation of the projected wave function и дана ссылка на http://aip.scitation.org/doi/10.1063/1.450026

Аватара пользователя
EvgeniX
Сообщения: 2780
Зарегистрирован: Пт апр 27, 2007 5:32 am

Re: Gaussian - архивная секция, параметр PUHF

Сообщение EvgeniX » Вт окт 10, 2017 5:51 am

В Gaussian 09 добавил iop(5/14=1), который должен отключать спиновую аннигиляцию.

В итоге для радикала HCC* для MP2 получил следующее:
Со спиновой аннигиляцией:
HF=-76.1760983\MP2=-76.4144725\PUHF=-76.2020311\PMP2-0=-76.4382021\S2=1.192476\S2-1=1.103847\S2A=0.880225\

Без неё:
HF=-76.1760983\MP2=-76.4144725\PUHF=-76.2020311\PMP2-0=-76.4424721\S2=1.192476\S2-1=1.103847\S2A=0.\
Диполные/квадрупольные моменты тоже становятся другими...

Gaussview в качестве энергии выводит значение "MP2=" (выходит, без проецирования спина). Но почему с отключенной спиновой аннигиляцией отличаются значения "HF" и "PUHF" ("MP2" и "PMP2-0")?

И ещё, этот самый iop(5/14=1) не работает, когда задаётся другой алгоритм scf (например, "scf(qc)").

Ответить

Вернуться в «квантовая химия и моделирование»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 7 гостей