计算化学公社

标题: 用Fortran 90改写Modern Quantum Chemistry附录B中的HF自洽场源码 [打印本页]

作者
Author: zhaixiaoyi001 时间: 2026-1-29 08:47
标题: 用Fortran 90改写Modern Quantum Chemistry附录B中的HF自洽场源码
本帖最后由 zhaixiaoyi001 于 2026-1-30 09:20 编辑

作为练习，用Fortran 90改写了HeH+自洽场代码，主要改动如下：
1、去掉了跳转，改写了do循环和if语句。
2、保留了原代码中的变量和电子积分计算方式。
3、用lapack做对角化，有部分变量是冗余的没有删除。

作者
Author: SeptemberMy 时间: 2026-1-29 19:39
如果你愿意的话可以把这个程序添加到一个Modern Quantum Chemistry汉化的项目中，直接提PR即可，https://github.com/nominhanggai/szaboqc

作者
Author: slxc920113 时间: 2026-1-29 23:46
代码的327行有语法错误，编译后运行会报错，改成

write(*,'(/,4X,A,F10.6,/)') 'DELTA(CONVERGENCE OF DENSITY MATRIX) = ', delta

复制代码

我重写了一个python的版本，逻辑和结构上会清晰很多。

(, 下载次数 Times of downloads: 7)

作者
Author: zhaixiaoyi001 时间: 2026-1-30 08:46

slxc920113 发表于 2026-1-29 23:46
代码的327行有语法错误，编译后运行会报错，改成

应该是那个39H字符串报错了，这块的Fortran 77语法我没修改，有的编译器识别不了，之前用ifx和gfortran编译倒是没报错。我再修改一下

作者
Author: zhaixiaoyi001 时间: 2026-1-30 10:06

SeptemberMy 发表于 2026-1-29 19:39
如果你愿意的话可以把这个程序添加到一个Modern Quantum Chemistry汉化的项目中，直接提PR即可，https://gi ...

又修改了一下代码，已提PR。

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