四阶Runge-Kutta法是工程计算中常用的一种求解微分方程的数值计算方法,具有精度高,易收敛等优点.本文在Feynman等人计算方法的基础上,用经典的四阶Runge-Kutta法来求解Thomas-Fermi-Dirac (TFD)方程,进一步提高原计算方法的计算精度.利用该方法求出了元素Cu的TFD方程数值解并计算出一些常见元素在Wigner-Seitz半径处的电子密度.
参考文献
| [1] | Thomas L H .[J].Proc Comb Phil Sco,1926,23:542. |
| [2] | Fermi E Z .[J].Physica,1928,48:73. |
| [3] | Dirac P A M .[J].Proceedings of the Cambridge Philosophical Society,1930,26:376. |
| [4] | 程开甲,樊启科,高占鹏.TF（或TFD）模型中原子的边界势及状态方程[J].物理学报,1984(02):176-192. |
| [5] | Latter Richard .[J].Chemical Physics,1956,24(02):280. |
| [6] | 程开甲,程漱玉.论材料科学的理论基础[J].自然科学进展,1996(01):12. |
| [7] | 程开甲,程漱玉.论固体膜形成纳米管的机理[J].稀有金属材料与工程,2003(10):773-776. |
| [8] | 程开甲,程漱玉.论位错的稳定存在[J].稀有金属材料与工程,2002(02):81-83. |
| [9] | 程开甲,程漱玉.论薄膜的超导电性[J].稀有金属材料与工程,2008(10):1693-1695. |
| [10] | 李世春.TFDC相图[J].自然科学进展,2003(11):1154-1159. |
| [11] | 王国灿,丁培柱,郑成德.量子力学中TFD方程边值问题的存在性[J].应用数学和力学,2000(02):215-220. |
| [12] | Feynman R P;Metropolis N;Teller E .[J].Physical Review,1949,75:1561. |
| [13] | Epele LN.;Canal CAG.;Ponciano JA.;Fanchiotti H. .Pade approximant approach to the Thomas-Fermi problem[J].Physical Review, A,1999(1):280-283. |
| [14] | 林成森.数值计算方法[M].北京:科学出版社,2005:212. |
| [15] | Metropolis N;Reitz J R .[J].Journal of Chemical Physics,1951,19(05):555. |
| [16] | 周玉;武高辉.材料分析测试技术[M].哈尔滨:哈尔滨工业大学出版社,2007:275. |
上一张
下一张
上一张
下一张
计量
- 下载量()
- 访问量()
文章评分
- 您的评分：
-
10%
-
20%
-
30%
-
40%
-
50%