BSSE расчет

[bar]

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

Решил провести 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]

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]

можешь посчитать ее вручную. Для комплекса 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))