[论文解读] Analytical solution of the Thomas-Fermi equation for atoms

M. Oulne|arXiv (Cornell University)|Nov 2, 2005

Quantum Mechanics and Non-Hermitian Physics被引用 3

一句话总结

本文提出了一种新的三参数解析试函数，用于处理Thomas-Fermi方程，显著提高了原子势能和电离能建模的精度。通过Ritz变分法，该解在 0 ≤ x ≤ 10 范围内的平均误差仅为 0.28%，且在 x=0 处的导数误差低于 2%，其电离能预测结果比以往任何解析近似方法更接近Hartree-Fock结果。

ABSTRACT

An approximate analytical solution of the Thomas-Fermi equation for neutral atoms is obtained, using the Ritz variational method, which reproduces accurately the numerical solution, in the range $0\leq x\leq50$, and its derivative at $x=0$. The proposed solution is used to calculate the total ionization energies of heavy atoms. The obtained results are in good agreement with the Hartree-Fock ones and better than those obtained from previously proposed trial functions by other authors.

研究动机与目标

• 开发一种更精确的Thomas-Fermi方程解析解，使其更符合数值解和物理边界条件。
• 相比现有试函数，改进对重原子总电离能的预测精度。
• 确保解能准确再现原点（x=0）处势能的导数，这是物理一致性的关键基准。
• 提供一种变分试函数，其在势能形状和电离能精度方面均优于以往的一、二、三参数形式。
• 以Hartree-Fock数据为参考标准，验证该方法在重原子（Z ≥ 92）中的有效性。

提出的方法

• 提出一种新的三参数试函数：φ(x) = (1 + α√x + βx e^{-γ√x})² e^{-2α√x}，旨在准确描述x=0处的行为并实现无穷远处的衰减。
• 应用Ritz变分法最小化拉格朗日泛函 L(φ) = ∫₀^∞ F(φ, φ', x) dx，其中 F 定义为 (1/2)(dφ/dx)² + (2/5)(φ^{5/2}/√x)。
• 通过拉格朗日乘子在最小化过程中施加粒子数约束 ∫ρ dv = Z。
• 利用Maple Release 9 数值优化变分参数 α、β 和 γ，以最小化能量泛函并满足归一化条件。
• 通过测试解在再现Thomas-Fermi方程数值解及其在 x=0 处导数的准确性，验证其性能。
• 通过公式 E = (12/7)(2/9π²)^{1/3} (dφ/dx)_{x=0} Z^{7/3} 计算电离能，基于在 x=0 处导数的推导结果。

实验结果

研究问题

• RQ1 能否构造一种新的解析试函数，使其在 0 ≤ x ≤ 50 范围内更精确地再现Thomas-Fermi方程的数值解？
• RQ2 所提出的解在 x=0 处的导数是否比以往试函数更精确，尤其考虑到其在电离能计算中的关键作用？
• RQ3 该新解在重原子（Z ≥ 92）总电离能预测中的精度与Hartree-Fock结果相比如何？
• RQ4 随着原子序数的增加，该新试函数能否保持或提升精度，而不同于以往近似方法的性能下降？
• RQ5 现有试函数为何无法以足够精度再现解及其导数？新形式如何克服这些局限性？

主要发现

• 在 0 ≤ x ≤ 10 范围内，47个点的平均误差仅为 0.28%，显著低于两参数函数（公式2）的1.13%误差，以及其它先前函数超过2.5%的误差。
• 解在 x=0 处的导数为 -1.61623647，与数值结果 -1.58807102 的相对误差低于2%。
• 新函数（公式21）的电离能预测在 Z=100 时误差仅为4.0%，在 Z=106 时为2.0%，且随Z增加而持续改善。
• 新解优于所有先前试函数（公式1、2、3、4、5），其中公式5因在 x=0 处导数无穷大被排除。
• 在 Z=109 时，新方法预测值为44,403哈特里，仅比Hartree-Fock值44,042高出0.8%，表明其在超重原子中具有极高精度。
• 该解在势能形状和原点处导数的精度均表现出色，使其在重元素原子性质计算中更具优势。