曹轲
,
任凤章
,
王天军
,
刘治军
,
田保红
,
李武会
稀有金属材料与工程
四阶Runge-Kutta法是工程计算中常用的一种求解微分方程的数值计算方法,具有精度高,易收敛等优点.本文在Feynman等人计算方法的基础上,用经典的四阶Runge-Kutta法来求解Thomas-Fermi-Dirac (TFD)方程,进一步提高原计算方法的计算精度.利用该方法求出了元素Cu的TFD方程数值解并计算出一些常见元素在Wigner-Seitz半径处的电子密度.
关键词:
Thomas-Fermi-Dirac方程
,
Runge-Kutta法
,
电子密度
,
Wigner-Seitz半径
