一个关于重叠积分计算的疑惑

nagami

大家好，请教个问题，在徐光宪的从头算书里有一个氨分子的HF计算例子，给出了一系列分子积分的数据，我用自己写的分子积分跑了下做个检验调试，也是用他上面的基组形式做的，用sto-3g跑计算的结果一致，他的基组和sto-3g略微差异，这没多大关系。我请一位朋友用gaussian跑了对应NH3分子的积分数据以及他自己写的量化程序计算，都做过对比，发现我们两者不一致。我和徐的一致，他和gaussian的一致。发现（p,s）重叠积分不一样，依照氨分子的对称性，ps积分应该保持对称，为什么gaussian计算的结果却不是呢？数据如下
是否构造基组的形式不同，但归一化数据又正确。

sobereva

注意分子朝向，会影响积分值。我不清楚你的NH3的朝向是什么样的。
高斯输出的重叠矩阵是没有错的。只要能与高斯的输出对上结果就肯定没问题。

卡开发发

已经验证，确实是朝向问题，我加上了nosymm就能够Gaussian重复徐光宪书上的结果，最近也在筹备写SCF的程序，大家共勉！

试试看下面的结果拿到Gaussian算，结果与徐老先生结果完全一致：

%chk=NH3.chk
# hf/sto-3g nosymm iop(3/33=1) scf=conver=6

Title Card Required

0 1
N                  0.16359854    0.16359854    0.16359854   
H                  1.15062790    0.00000000    0.00000000
H                  0.00000000    1.15062790    0.00000000
H                  0.00000000    0.00000000    1.15062790

输出结果：
*** Overlap ***
                1             2             3             4             5
      1  0.100000D+01
      2  0.235038D+00  0.100000D+01
      3  0.000000D+00  0.000000D+00  0.100000D+01
      4  0.000000D+00  0.000000D+00  0.000000D+00  0.100000D+01
      5  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.100000D+01
      6  0.581692D-01  0.493788D+00  0.427508D+00 -0.708587D-01 -0.708587D-01
      7  0.581692D-01  0.493788D+00 -0.708587D-01  0.427508D+00 -0.708587D-01
      8  0.581692D-01  0.493788D+00 -0.708587D-01 -0.708587D-01  0.427508D+00
                6             7             8
      6  0.100000D+01
      7  0.213322D+00  0.100000D+01
      8  0.213322D+00  0.213322D+00  0.100000D+01

Shannon

确实是对称性问题。  刚去检查了一下那个代码，结果发现昨天的代码 的输入是直接从高斯的输出中提取出来的， 坐标已经被高斯处理过一遍了，因此和高斯一致却和 徐老先生的不一样。
从对称性来看，如果按照徐先生书上的 几何坐标， N的 py轨道应该和三个氢原子有完全一样的重叠。

nagami

sobereva 发表于 2015-5-3 20:53
注意分子朝向，会影响积分值。我不清楚你的NH3的朝向是什么样的。
高斯输出的重叠矩阵是没有错的。只要能 ...

谢谢Sob老师