我是一名计算机科学专业的毕业生,我正在编写一个电子结构计算程序,我一直坚持使用 6-31G 基组为具有更高原子序数的原子(p、d、f 壳)形成基函数。
例如,考虑一个高斯 g 为
g(alpha) = (2*alpha/pi)^(3/4) * exp(-alpha*r^2).
查看 6-31G 基组的 Carbon (1s2 2s2 2p2) 基组数据
C 0
S 6 1.00
0.3047524880D+04 0.1834737132D-02
0.4573695180D+03 0.1403732281D-01
0.1039486850D+03 0.6884262226D-01
0.2921015530D+02 0.2321844432D+00
0.9286662960D+01 0.4679413484D+00
0.3163926960D+01 0.3623119853D+00
SP 3 1.00
0.7868272350D+01 -0.1193324198D+00 0.6899906659D-01
0.1881288540D+01 -0.1608541517D+00 0.3164239610D+00
0.5442492580D+00 0.1143456438D+01 0.7443082909D+00
SP 1 1.00
0.1687144782D+00 0.1000000000D+01 0.1000000000D+01
在这种情况下,碳有多少个基函数?我一直在关注一个系列讲座,该系列对一个简单的 H2 分子进行 HF SCF 计算,其中每个氢原子都有两个基函数,根据 6-31G 基组。
H 0
S 3 1.00
0.1873113696D+02 0.3349460434D-01
0.2825394365D+01 0.2347269535D+00
0.6401216923D+00 0.8137573261D+00
S 1 1.00
0.1612777588D+00 1.0000000
我可以将两个基函数写成形式
b1= 0.0334*g(18.73) + 0.2347*g(2.825) + 0.8137*g(0.64);
b2= g(0.161);
以类似的方式,我将如何编写碳原子的基函数。任何建议都会非常有帮助,因为我被困在这几天并且无法为此找到任何实质性资源。
我确实偶然发现了这篇上一篇文章,但很难弄清楚 Carbon 的 9 个基函数是如何产生的。