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

查看数: 275 | 评论数: 2 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2025-9-22 16:47

思路是利用Obara-Saika recurrence relations降低Primitive GTF的角动量值迭代。代码见附件。解压直接julia HF.jl即可输出S矩阵。与Gaussian16的对比如下，部分矩阵元相等，但是部分矩阵元不等。请问问题出在哪里？ ...

回复 Reply

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

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

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

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

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

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