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

RTChou · 发表于 Post on 2025-9-22 16:47:04

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

hebrewsnabla · 发表于 Post on 2025-9-22 23:07:30

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

RTChou · 发表于 Post on 2025-9-23 01:02:25

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[算法与编程] Julia写的Obara Saika计算S矩阵结果和gaussian有部分出入

浏览过的版块