Dftatom公司

Dftatom:一个稳健和通用的薛定谔和狄拉克解算器,用于原子结构计算。提出了径向薛定谔方程、狄拉克方程和Kohn-Sham方程的鲁棒通用解。该公式允许一般的势和网格:均匀的,指数的,或其他由节点分布和导数函数定义的。对于给定的网格类型,可以通过增加网格点的数目来系统地控制收敛。径向积分是使用渐近形式,龙格库塔和隐式Adams方法的组合。为了提高鲁棒性和速度,本征函数由二分法和摄动法相结合确定。采用外向泊松积分提高核心区的精度,使重原子(如铀)的总能量的绝对精度达到$10^{-8}$Hartree。详细的收敛性研究,并提供计算参数,以达到通常要求的精度在实践中。与原子序数为$Z=1-92$的分析和当前基准密度函数结果进行了比较,验证并提供了对当前基准的改进。一个高效的、模块化的fortran95实现dfatom作为开放源代码提供,包括示例、测试和与其他语言接口的包装器;其中特别强调单个例程的独立性(无全局变量)、可重用性和通用性。