RHF方法编程笔记-4(完) 一些简单物理量的计算

zyniso

       本次主要介绍：在完成了RHF计算后，如何获得一些基本的物理量，如体系的动能，势能，<V>/<T>比值，Mulliken电荷，偶极矩，四极矩，八极矩，十六极矩。       至此，RHF方法编程笔记系列至此全部介绍完毕，笔者会在不久后上传一个完整版本的RHF计算的C++程序源代码。今后将会以不定期更新的方式介绍其他量化编程的内容：如UHF（其实从RHF过渡到UHF，笔者在编程过程中，涉及的代码改动不超过20行，各位老师和同学可以先行尝试）；SCF计算的加速方法（DIIS方法，一般会使得SCF计算在10步内收敛）；MP2（涉及一步标度为N^5的积分变换）；其他双电子积分的计算方法（O-S，Rys等）；DFT（至少一年后。。。）；解析导数的计算等等。
       最后，希望本系列的笔记能对热爱编程的各位有所帮助，任何涉及笔记内容的问题，欢迎随时交流讨论。

granvia

“体系总的势能”应该还包括核与核直接的静电排斥，最好能指明一下