aug(df,pd)-6-311G**

[1177]

Подскажите, пожалуйста, как в Gaussian'е создать следующие базисы:

The aug(d,p)-6-311G** basis set is the 6-311G** basis set augmented with diffuse d functions on carbon atoms [αd(C) = 0.1565] and diffuse p functions on hydrogen atoms [αp(H) = 0.1875]. The aug(df,pd)-6-311G** basis set is the 6-311G** basis set augmented with the diffuse d functions and diffuse f functions on carbon atoms [αf(C) = 0.2] and the diffuse p functions and diffuse d functions on hydrogen atoms [αd(H) = 0.25].

Оригинальная работа в обменнике. Заранее благодарю

[Yu/2]

http://gaussian.com/g_ur/k_gen.htm

[1177]

Чего-то я не понимаю

Написал:

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

# HF/6-311G** ExtraBasis Test
 
[C1 dimer: aug(d,p)-6-311G**]
 
  0  1
C 0   -0.267339    0.282701    0.006189
H 0    0.741460    0.709728    0.021331
H 0   -0.582212    0.128477   -1.031752
H 0   -0.961535    0.970190    0.501497
H 0   -0.264416   -0.677719    0.532974
C 0    2.517715   -2.389123   -0.018555
H 0    3.564342   -2.712326   -0.019444
H 0    2.473461   -1.295099   -0.022909
H 0    2.017656   -2.770613    0.878657
H 0    2.013420   -2.777740   -0.910213

C 0 
D 1 1.00 
 0.1565000000D+00  1.0000000000D+00
****

H 0 
P 1 1.00 
 0.1875000000D+00  1.0000000000D+00
****

Gaussian выдает 94 bf (140 prim. gaus.), а должно быть [J. Phys. Chem. A, Vol. 102, No. 12, 1998, 2092] 118 bf (160 prim. gaus.)

[1177]

[Darth Vasya]

Так, навскидку... Не может это быть связано с дихотомией pure vs. cartesian d-functions?

[1177]

Возможно...

[1177]

Все равно не получается...
Может я изначально иду не тем путем?

[alxyppv]

Gaussian очень капризен в синтаксисе. Уберите пустую строчку между дополнительными функциями для разных атомов, т.е. надо

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

......
H 0    2.013420   -2.777740   -0.910213

C 0 
D 1 1.00 
 0.1565000000D+00  1.0000000000D+00
****
H 0 
P 1 1.00 
 0.1875000000D+00  1.0000000000D+00
****

а не

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

.....
H 0    2.013420   -2.777740   -0.910213

C 0 
D 1 1.00 
 0.1565000000D+00  1.0000000000D+00
****

H 0 
P 1 1.00 
 0.1875000000D+00  1.0000000000D+00
****

Иначе он (скорее всего) просто не считывает функции для водорода.

При этом пустая строка между координатами и доп. базисными функциями должна быть сохранена:)

[1177]

Сработало!!!
Спасибо!