Ringo——我的第一个量化程序

目前 Modern Quantum Chemistry 上氢分子和氦氢离子(HeH+)的两个例子已经学完了。

然后我就尝试写出了自己的第一个量化程序,用于 HeH+ 的 RHF/STO-3G 的计算。

程序是用 Fortran 写的。

代码和文档上传到了Github:https://github.com/St-Maxwell/Ringo

因为基本上所有内容都在文档中,所以这里不再详细介绍了。

如果有人对量化编程有兴趣,希望我的代码和文档能帮助到你。

随着我未来对理论化学的学习,会逐渐为这个代码的增加更多的内容。


运行结果

与Gaussian的结果的对比

2 thoughts on “Ringo——我的第一个量化程序”

  1. 个人建议哈,或许你在电子积分这一部分花费了很多时间,所以你可以尝试一下pyscf这个package.(http://sunqm.github.io/pyscf/) 然后通过调用numpy里面的方程,写起来更容易一些。这里面的处理好在你可以直接调用不同形式的电子积分,并以矩阵的形式呈现出来,对你写算法或许帮助更大。而且,有一个小建议,你可以用矩阵乘法和数组之间的操作来避免连续的loop,感觉这个在运算更大的体系的时候效率没那么高吧。不过讲真你看的挺快的,这么早开始写RHF真不容易。个人觉得MQC这本书最值得码出来的部分就是RHF,CIS 和MP2 吧,加油。

    1. 谢谢评论。现在写看书写代码都算是课余的事情,所以纯粹是写给自己玩的。这种情况下我会倾向于自己造轮子。当然,也要学习别人的代码,我会去看pyscf是如何写的,也还找了一个大佬写的电子积分程序。(假如以后搞科研要写程序的话,感觉老老实实去调用别人的代码比较好23333333)

发表评论

电子邮件地址不会被公开。 必填项已用*标注