BSSE расчет

обсуждение вопросов физической химии и химической физики
physical chemistry and chemical physics: discussions for professionals
Ответить
bar
Сообщения: 2
Зарегистрирован: Чт май 13, 2010 8:50 am

BSSE расчет

Сообщение bar » Чт май 13, 2010 9:32 am

Уважаемые форумчане!

Решил провести BSSE коррекцию энергии связи С-H в метане с помощью программы Gaussian 03. Просмотрев немногочисленные примеры создания gjf-файлов для BSSE-расчетов, никак не мог понять, как правильно составить gjf-файл, используя следующую строку:

# opt b3lyp/dgdzvp geom=connectivity counterpoise=2

атомы с декартовыми координатами

Не определенной ясности, как использовать для этой цели dummy-атомы...

Пробовал провести для первичной практики BSSE-расчет H-H связи на простой молекуле водорода, составив два разных gjf-файлы:

1. # opt b3lyp/dgdzvp geom=connectivity Counterpoise=2

0 1 0 2 0 2
H -2.19584760 0.11954154 0.04746888 1
H -0.72390474 0.44611010 0.02773632 2

2. opt hf/3-21g geom=connectivity counterpoise=2

0 1 0 2 0 2
X -2.15080920 -2.47572812 0.00000000 1
H -3.45080920 -2.47572812 0.00000000 2
H -1.33956084 -1.43344174 -0.24240647 1

Но неизменно возникало сообщение об ошибке: "Inconsistent input #3 in SPerMO", сведения о смысле которой отсутствуют.

Просьба посодействовать в этом вопросе

Ushi
Сообщения: 6
Зарегистрирован: Пн янв 26, 2009 9:28 pm

Re: BSSE расчет

Сообщение Ushi » Чт май 13, 2010 6:25 pm

BSSE -the correction for the error caused by superposition of the basis sets of subsystems или учетом ошибки наложения
базисных наборов подсистем, как я знаю, применяется в случае, если у вас есть комплекс состоящий из двух или более подсистем, там хвост базиса одной молекулы влияет на другую молекулу системы (ошибка уменьшается при более полном базисе)
пишется так, например, для системы H2 Ar
# b3lyp/6-31G counterpoise=2

01 01 01
H x.xxx y.yyy z.zzz 1
H x.xxx y.yyy z.zzz 1
Ar x.xxx y.yyy z.zzz 2

YuraM
Сообщения: 441
Зарегистрирован: Пн ноя 26, 2007 11:07 pm

Re: BSSE расчет

Сообщение YuraM » Пт май 14, 2010 2:12 am

можешь посчитать ее вручную. Для комплекса AB потребуется 5 расчетов.
1. Энрегия A-B
2. Энергия A-(dummyB), где при описании сиситемы на место B кладутся мнимые атомы (те атомов там нет), а базисные ф-ции есть.
3. Энергия (Adummy)-B
4. Энергия A
5. Энергия B

Тогда BSSE-corrected энергия будет - E(BBSE_corrected) = E(A-B)-(E(A-dummyB) - E(A))-(E(Adummy-B) - E(B))
Кто смел тот и съел

Ответить

Вернуться в «физическая химия / physical chemistry»

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

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