计算化学公社

标题: 一个关于重叠积分计算的疑惑 [打印本页]

作者
Author: nagami 时间: 2015-5-3 14:12
标题: 一个关于重叠积分计算的疑惑
本帖最后由 nagami 于 2015-5-3 14:37 编辑

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

作者
Author: sobereva 时间: 2015-5-3 20:53
注意分子朝向，会影响积分值。我不清楚你的NH3的朝向是什么样的。
高斯输出的重叠矩阵是没有错的。只要能与高斯的输出对上结果就肯定没问题。

作者
Author: 卡开发发 时间: 2015-5-4 01:55
已经验证，确实是朝向问题，我加上了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

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

作者
Author: nagami 时间: 2015-5-4 09:11

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

谢谢Sob老师

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