ECPを番号で指定する方法について
無機固体を研究している大学4年のものです。gaussian16を用いて計算を行っています。
現在、ECPをある特定の原子にのみ割り当てたいのですが、できなくて困っております。
b3lyp/genECPのキーワードで座標の下に基底関数を
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 19 20 23 24 25 0
6-31g
****
17 18 21 22 0
6-31+g(d)
****
Sn 0
S 1 1.00
0.5418000 1.0000000
S 1 1.00
0.3784000 1.0000000
S 1 1.00
0.0926000 1.0000000
P 1 1.00
1.0470000 1.0000000
P 1 1.00
0.1932000 1.0000000
P 1 1.00
0.0630000 1.0000000
P 1 1.00
0.0174 1.0
D 1 1.00
0.1860 1.0
****
SN 0
SN-ECP 3 46
f potential
5
0 0.2169021 -0.0665380
1 33.3097182 -29.7469416
2 9.8846832 -61.8307278
2 2.5019593 -18.9410634
2 0.7833323 -2.2384982
s-f potential
5
0 155.6062293 2.9482777
1 56.5145464 63.7614718
2 22.6583648 230.5276788
2 6.5837080 85.0713314
2 1.5947702 28.8215644
p-f potential
5
0 40.4666705 2.6941946
1 13.6218440 44.8271079
2 4.1092688 57.6223370
2 1.5046137 -44.6853312
2 1.3646762 57.1351815
d-f potential
5
0 60.1915052 3.0259297
1 20.4961142 39.5341071
2 8.4438328 68.4686504
2 2.6103790 29.1123764
2 0.6501317 5.8428577
27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 0
lanl2dz
のように書いてあります。計算結果のファイルを見るとECPを考慮した基底関数のMulliken chargesを確認すると
内殻の電子も計算しているように見えます。何か問題がわかる方いらしたらコメントお願いします。