Szabo书中分子积分与gaussian09算出的不符

Zork

Szabo的《Modern Quantum Chemistry》中有给出 氢分子 的 rhf/sto-3g 的计算结果。我在写程序能实现书一样的结果。但是我用gaussian09 做了相同的计算，并输出了分子积分的结果，发现和Szabo书上给出的结果不一样。（结果如下图）

顺便说一下，我遇到的另一个问题。我写的程序在算稍微大一点的分子（比如 H2O，没错，水就是稍微大的分子）的时候，能量会在一个大数和一个小数之间震荡。不是个例，我试了好几个稍微大的分子，都是这样的。显然这两个数值都不是正确的结果。（水的能量大概在-74hartree，我算的结果如下图。）

St_Maxwell

可能是因为基组定义不同，书里的STO-3G的系数和高斯里的不一样。

liyuanhe211

看到ID好害怕

柒月小鱼

看到ID惊了

Zork

St_Maxwell 发表于 2019-5-11 10:43
可能是因为基组定义不同，书里的STO-3G的系数和高斯里的不一样。

是一样的，我用gen参数，在在下面手工输入基组，然后用iop把基组信息输出，和我用sto-3g这个参数直接指定基组，得到的结果是一样的。而且，szabo书上的氢的基组，和emsl库里的数据是一样的。

Zork

liyuanhe211 发表于 2019-5-11 15:34
看到ID好害怕

忽略这些细节好吗。

ene

ID瞩目

Zork

ene 发表于 2019-5-11 19:01
ID瞩目

我错了，这个id似乎不能改了。不然我立马换一个。

ene

Zork 发表于 2019-5-11 19:05
我错了，这个id似乎不能改了。不然我立马换一个。

完全没有这个意思，楼主不要介意，正常用就好

granvia

看来这里老年人还是蛮多的

Warm_Cloud

高斯加上nosymm，检查一下你的双电子积分是否做对了。

Zork

问题解决了，，，原来原因是，g09 默认的单位制中长度是 埃 。 而 szabo书上的长度单位是 a.u. 。高斯加上 units(au) 参数（表明输入的坐标使用原子单位制），结果就一样了。 感谢各位帮助，打扰了。学量化途中一定会遇到其他问题，望各位不吝赐教（ps: 起名时我年少无知。）

beefly

我感觉，在不引起歧义的情况下（特别是涉及不好的事情），任何id都可以用。比如，我用了爱因斯坦的头像，爱大爷从未托梦说不行。