Julia写的Obara Saika计算S矩阵结果和gaussian有部分出入

RTChou

思路是利用Obara-Saika recurrence relations降低Primitive GTF的角动量值迭代。代码见附件。解压直接julia HF.jl即可输出S矩阵。与Gaussian16的对比如下，部分矩阵元相等，但是部分矩阵元不等。请问问题出在哪里？我核对过gbs文件一致，也反复检查过矩阵元计算，但是并没有发现具体出错原因

RTChou

hebrewsnabla 发表于 2025-9-22 23:07
建议先排查有的对角元是 0 的问题……这显然不是有出入，而是你的结果本身就有问题了。

您好，您指出了很重要的问题。我重写了一遍Sij，解决了对角项不是1的问题。对于我帖子中的问题，是源于GTF乘积后中心移动我并没有考虑到。目前已经能够和Gaussian结果对上

hebrewsnabla

建议先排查有的对角元是 0 的问题……这显然不是有出入，而是你的结果本身就有问题了。