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

谢谢Sob老师

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

已经验证，确实是朝向问题，我加上了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

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