计算化学公社

标题: Julia写的Obara Saika计算S矩阵结果和gaussian有部分出入 [打印本页]

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

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

作者
Author: RTChou 时间: 2025-9-23 01:02

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

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

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