Быстрый тест метода

вопросы строения молекул и квантовой химии
kbob
Сообщения: 1705
Зарегистрирован: Сб окт 17, 2009 9:48 am

Re: Быстрый тест метода

Сообщение kbob » Ср окт 25, 2017 7:18 am

Droog_Andrey писал(а):
Вс окт 22, 2017 8:08 pm
Droog_Andrey писал(а):
Чт окт 12, 2017 3:53 pm
3) вертикальная энергия возбуждения C32+ линейной геометрии (R = 1.27Å) из 1Σg в 3Σu без учёта ZPE (эВ)
CASSCF(10,15) / cc-pVQZ даёт +0.33 эВ.
Это максимальный вариант? В молекуле 10 электронов, на каждый атом (1S + 2S + 2Px + 2Py + 2Pz) орбиталей. Причем как я заметил aug базис улучшения не дал.
If you are not part of the solution, you are part of the precipitate.

Аватара пользователя
Droog_Andrey
Сообщения: 2669
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Быстрый тест метода

Сообщение Droog_Andrey » Ср окт 25, 2017 7:23 am

В частице 16 электронов. Ну и для дикатиона диффузные функции дают мало толку.
2^74207281-1 is prime!

kbob
Сообщения: 1705
Зарегистрирован: Сб окт 17, 2009 9:48 am

Re: Быстрый тест метода

Сообщение kbob » Ср окт 25, 2017 7:53 am

Droog_Andrey писал(а):
Ср окт 25, 2017 7:23 am
В частице 16 электронов. Ну и для дикатиона диффузные функции дают мало толку.
Понятно - все электроны кроме 1S в активном пространстве. Для дикатиона диффузные функции не дают толку. а вот для возбужденного состояния дадут - его энергия понизится, но незначительно как я заметил.
If you are not part of the solution, you are part of the precipitate.

Аватара пользователя
Droog_Andrey
Сообщения: 2669
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Быстрый тест метода

Сообщение Droog_Andrey » Ср окт 25, 2017 11:36 pm

По-хорошему надо CASPT3 считать.
2^74207281-1 is prime!

kbob
Сообщения: 1705
Зарегистрирован: Сб окт 17, 2009 9:48 am

Re: Быстрый тест метода

Сообщение kbob » Чт окт 26, 2017 6:06 am

Droog_Andrey писал(а):
Ср окт 25, 2017 11:36 pm
По-хорошему надо CASPT3 считать.
Я в ORCA MRCI посчитал на CASSCF(6,9) на большее времени не хватило.
If you are not part of the solution, you are part of the precipitate.

Аватара пользователя
Droog_Andrey
Сообщения: 2669
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Быстрый тест метода

Сообщение Droog_Andrey » Сб янв 06, 2018 7:18 pm

Посчитал C32+ в CASSCF(10,23,RAS(0,0,2,12))/maug-cc-pVQZ.
То есть замораживаем 1s, берём следующие 23 орбитали, но на верхние 12 пускаем не более двух электронов.
Получилось:
Синглет -112.61868
Триплет -112.60655
Энергия +0.33 эВ.
2^74207281-1 is prime!

Аватара пользователя
Droog_Andrey
Сообщения: 2669
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Быстрый тест метода

Сообщение Droog_Andrey » Вс июл 01, 2018 10:11 pm

Droog_Andrey писал(а):
Чт окт 12, 2017 3:53 pm
Придумал три простые задачи, расположенные по возрастанию вычислительной сложности:
1) сродство к электрону атома водорода (эВ)
2) межъядерное расстояние в потенциальном минимуме димера натрия (Å)
3) вертикальная энергия возбуждения C32+ линейной геометрии (R = 1.27Å) из 1Σg в 3Σu без учёта ZPE (эВ)
Пункт 2) лучше заменить на межъядерное расстояние в триплетной (3Πu) молекуле фтора (около 2.2Å). А чтобы не упускать при этом семейство металлов, пункт 1) заменяем на сродство к электрону атома лития (0.618 эВ).

В итоге, правда, оказались в рамках второго периода. Но, тем не менее, набор задач весьма капризен к методу. Ни один известный мне DFT-функционал и близко не справился со всеми тремя.
2^74207281-1 is prime!

Аватара пользователя
Droog_Andrey
Сообщения: 2669
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: Быстрый тест метода

Сообщение Droog_Andrey » Пт ноя 26, 2021 5:30 am

Чтобы не оставаться в рамках второго периода, можно добавить высоту барьера H2 + H = [H---H---H] = H + H2. А чтобы не было засилия энергетических характеристик, стоит добавить межъядерное расстояние в потенциальном минимуме катиона HeBe+ в основном состоянии (2Σ+). Итого получим (в порядке возрастания числа электронов в системе):

1) высота барьера H2 + H = [H---H---H] = H + H2 без учёта ZPE, референсное значение 0,421 эВ
2) сродство к электрону атома лития, референсное значение 0,618 эВ
3) межъядерное расстояние в потенциальном минимуме катиона HeBe+ в основном состоянии 2Σ+, референсное значение 2,94Å
4) вертикальная энергия возбуждения C32+ линейной геометрии (R = 1.27Å) из 1Σg в 3Σu без учёта ZPE, референсное значение 0,33 эВ
5) междъядерное расстояние в потенциальном минимуме в молекуле фтора в низшем триплетном состоянии 3Πu, референсное значение 2,2Å
2^74207281-1 is prime!

Ответить

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

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

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