计算化学公社

标题: Szabo书中分子积分与gaussian09算出的不符 [打印本页]

作者
Author: Zork 时间: 2019-5-11 10:33
标题: Szabo书中分子积分与gaussian09算出的不符
Szabo的《Modern Quantum Chemistry》中有给出氢分子的 rhf/sto-3g 的计算结果。我在写程序能实现书一样的结果。但是我用gaussian09 做了相同的计算，并输出了分子积分的结果，发现和Szabo书上给出的结果不一样。（结果如下图）

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

作者
Author: St_Maxwell 时间: 2019-5-11 10:43
本帖最后由 St_Maxwell 于 2019-5-11 15:32 编辑

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

作者
Author: liyuanhe211 时间: 2019-5-11 15:34
看到ID好害怕

作者
Author: 柒月小鱼 时间: 2019-5-11 16:18
看到ID惊了

作者
Author: Zork 时间: 2019-5-11 18:57

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

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

作者
Author: Zork 时间: 2019-5-11 18:58

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

忽略这些细节好吗。

作者
Author: ene 时间: 2019-5-11 19:01
ID瞩目

作者
Author: Zork 时间: 2019-5-11 19:05

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

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

作者
Author: ene 时间: 2019-5-11 19:07

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

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

作者
Author: granvia 时间: 2019-5-11 20:17
看来这里老年人还是蛮多的

作者
Author: Warm_Cloud 时间: 2019-5-11 20:58
高斯加上nosymm，检查一下你的双电子积分是否做对了。

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

作者
Author: beefly 时间: 2019-5-12 18:01
本帖最后由 beefly 于 2019-5-12 18:07 编辑

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)