Hellmann-Feynman力的ab initio计算及在化学键研究中的应用 Ⅱ.AB INITIO FORTRAN IV通用程序ABHF

本文根据我们推导出的Gauss函数的力积分的公式和力分析的公式,在ab initio FORTRAN IV程序MQAB-80的基础上编制了可计算H-F力的ab initio FORTRAN IV通用程序ABHF,本程序可使用多中心GTO基组,可计算包含20个原子42个原子轨道的开壳层和闭壳层的分子中作用于指定原子指定方向的 Hellmann-Feynman力,依据LCAO-MO理论提供力分析的数据,使用程序ABHF和多中心GTO基组计算H-F力,在HF分子的例子中得到作用于F核的净力是—0.0225753a.u.,作用于H核的净力是—0.0570231 a.u.,超过了用4—31G和6—31G基组的计算结果,也超过了用Clementi的HF较精确的波函数(~(16)STO)的计算结果,在能量、偶极矩、维里值的计算上取得了全面改善。