一种TFD方程的数值解法及一些常见元素的电子密度求解.pdfVIP

稀有金属材料与工程 Vo 1.41 , No.3 第41 卷第3 期 RARE METAL MATERlALSAND ENGINEERING March2012 2012 年 3 月 一种 TFD 方程的数值解法及一些常见元素 的电子密度求解 曹 柯 1 ，任凤章 1 气王天军 1 ，刘治军 1 ，田保红 1 ，李式会 1 (1.河南科技大学，河南洛阳 471003) (2. 洞南省有色金属材料科学与加工技术重点实验室，河南洛阳 471003) 摘 要z 四阶 Runge-Kutta 法是工程计算中常用的一种求解微分方程的数值计算方法，具有精度高，易收敛等优点.本 文在 Feynman 等人计算方法的基础上，用经典的四阶 Runge-Kutta 法来求解Thomas-Fermi-Dirac (TFD) 方程，进一步 提高原计算方法的计算精度。利用该方法求出了元素 Cu 的 TFD 方程数值解并计算出一些常见元素在 Wigner-Seitz 半径 处的电子密度。 关键词: Thomas-Fermi-Dirac 方程; Runge-Kutta 法:电子密度; Wigner-Seitz 半径 中图法分类号: TGll l.1 文献标识码:A 文章编号: 1002-185X(2012)03-0506-04 亏 1926 年和 1928 年， Thomas[l]和 Fermi[2]各自独立 均要用到元素的电子密度数据，因而 TFD 方程的数值 解研究一直是计算材料科学中研究的热点[HAle 地提出了一个描述电子在原子中分布的统计模型，称 为 Thomas-Fermi (简称 TF) 模型。 1930 年， Dirac[3] 在数值分析中， Runge-Kutta 法是用于模拟常微分 将电子的交换作用引入到 TF 模型中，将其发展为 方程解的一类重要的显式或隐式迭代法，具有精度高， Thomas-Fermi-Dirac (简称 TFD) 模型，被广泛地称 收敛稳定(在一定条件下〉等优点，在工程计算中得 到了广泛的应用[14] 。本文利用经典的四阶 Runge-Kutta 为电子的统计理论。 1984 年，程开甲等人问根据经典 TFD 模型，在 Latter 等人口]的计算基础上对 TFD 方程 法对 TFD 方程进行了数值求解，重新计算了一些常见 的数值求解方法作了新的改进，并于 1996 年求解出了 元素在 Wigner-Seitz 半径处的电子密度，这将为 TFDC 一些常见元素的电子密度，对 TFD 模型进行了新的挖 电子理论在材料研究和设计中的应用提供方便。 掘，提出了后来被广泛称为 Thomas-Fermi-Dirac­ 1 TFD 方程 Cheng (简称 TFDC) 的电子理论模型闷。在 TFDC 电 Thomas[l]和 Fermi[2]建立的 TF 模型，其电子密度 子理论电子密度连续的边界条件下，程开甲等人已在 凝聚态材料方面成功地开展了许多应用研究z 如计算薄