Разные basis sets для разных элементов?

[Леонид]

В одной из статей видел что используют в вычислениях разные basis sets для различных атомов в молекуле, к примеру метал вычисляется одним, азот другим а углерод третим basis set.
Дайте пожалуста пример файла где использован такой метод вычисления.

[Shiver]

Каким пакетом пользуетесь-то?

[Леонид]

Прошу простить
Я действительно не указал софт. Пользуюсь Gaussian98.

[EvgeniX]

Вот пример со внешним ECP базисом для Zr:
------------------------
%mem=512Mb
#p opt=z-matrix pbepbe/genecp

Cp2ZrCl2 (PBE/ECP28MHF)

0 1
C
C 1 1.4
C 2 1.4 1 108.0
C 1 1.4 2 108.0 3 -1.0
C 4 1.4 1 108.0 2 1.0
H 1 1.1 2 126.0 4 178.0
H 2 1.1 3 126.0 1 177.0
H 3 1.1 2 126.0 1 179.0
H 4 1.1 5 126.0 1 -177.0
H 5 1.1 4 126.0 1 -179.0
Zr 1 2.6 2 74.0 3 65.0
C 11 2.6 1 173.0 6 1.0
C 12 1.4 11 74.0 1 -125.0
C 13 1.4 12 108.0 11 -65.0
C 12 1.4 13 108.0 14 1.0
C 15 1.4 12 108.0 13 -1.0
H 12 1.1 13 126.0 15 -178.0
H 13 1.1 14 126.0 12 -177.0
H 14 1.1 13 126.0 12 -179.0
H 15 1.1 16 126.0 12 177.0
H 16 1.1 15 126.0 12 179.0
Cl 11 2.5 1 87.0 6 49.0
Cl 11 2.5 1 88.0 6 -51.0

C 0
aug-cc-pvdz
****
Zr 0
S 3 1.0
5.8737890 -0.8394311
4.2872700 1.4322944
1.4641370 0.3960967
S 1 1.0
0.8312450 1.0
S 1 1.0
0.3623350 1.0
S 1 1.0
0.0793570 1.0
S 1 1.0
0.0339200 1.0
S 1 1.0
0.0110000 1.0
P 2 1.0
2.8742240 -4.1490377
2.1199010 3.8441952
P 2 1.0
0.8513640 0.5739450
0.4372920 0.4590635
P 1 1.0
0.2029040 1.0
P 1 1.0
0.0632590 1.0
P 1 1.0
0.0225230 1.0
D 4 1.0
2.5801370 -0.0915734
1.8553960 0.2237318
0.6807540 0.4769343
0.2566910 0.5279852
D 1 1.0
0.0894000 1.0
D 1 1.0
0.0300000 1.0
****
Cl 0
aug-cc-pvdz
****
H 0
aug-cc-pvdz
****

Zr 0
Zr-ECP 4 28
G POTENTIAL
1
2 1.000000 0.000000
S-G POTENTIAL
2
2 7.820000 146.202289
2 3.910000 18.323348
P-G POTENTIAL
2
2 6.990000 99.859830
2 3.500000 13.655859
D-G POTENTIAL
2
2 5.440000 32.237066
2 2.660000 7.660252
F-G POTENTIAL
2
2 7.510000 -21.069954
2 3.770000 -3.113384
------------------------

[Леонид]

Даный ввод логичен?
Ради примера использовал просто молекулу этанола:
%mem=6MW
%nproc=1
# opt=(tight) hf/gen scf=(qc,tight,direct,maxcycle=1000)

ethanol example

0 1
C
H 1 B1
H 1 B2 2 A1
C 1 B3 2 A2 3 D1
H 4 B4 1 A3 2 D2
H 4 B5 1 A4 2 D3
H 4 B6 1 A5 2 D4
O 1 B7 4 A6 5 D5
H 8 B8 1 A7 4 D6

B1 1.08524010
B2 1.08524010
B3 1.52454066
B4 1.08205877
B5 1.08205877
B6 1.08362390
B7 1.44434552
B8 0.96606627
A1 108.14929326
A2 110.18470376
A3 109.49555185
A4 109.49555185
A5 110.87721903
A6 106.19779937
A7 110.85780472
D1 120.48957464
D2 61.00929257
D3 179.73288693
D4 -59.62891025
D5 -59.36179718
D6 180.00000000

O 0
6-311++g(2p)
****

C H 0
3-21g
****