[論文レビュー] 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方程式の新しい3パラメータ解析的試行関数を提示する。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）におけるポテンシャルの微分係数を物理的整合性の重要な基準として正確に再現することを目的とする。
• ポテンシャルの形状とイオン化エネルギーの両方において、1パラメータ、2パラメータ、および従来の3パラメータ形式を凌駕する変分試行関数の提供を目的とする。
• 重い原子（Z ≥ 92）のハートリー・フォックデータを基準基準として用いて、手法の妥当性を検証することを目的とする。

提案手法

• 原点での正しい挙動と無限大における減衰を捉えることを目的として、新しい3パラメータ試行関数 φ(x) = (1 + α√x + βx e^{-γ√x})² e^{-2α√x} を提案する。
• ラグランジアン汎関数 L(φ) = ∫₀^∞ F(φ, φ', x) dx を最小化するRitz変分法を適用し、F は (1/2)(dφ/dx)² + (2/5)(φ^{5/2}/√x) として定義される。
• 粒子数制約 ∫ρ dv = Z は、最小化の過程でラグランジュ乗数を用いて強制する。
• 変分パラメータ α, β, γ は、正規化条件を満たしながらエネルギー汎関数を最小化するように、Maple Release 9 を用いて数値的に最適化する。
• 得られた解の正確性は、Thomas-Fermi方程式の数値解および x=0 における微分係数の再現性をテストすることで検証される。
• イオン化エネルギーは、導関数 (dφ/dx)_{x=0} を用いた式 E = (12/7)(2/9π²)^{1/3} (dφ/dx)_{x=0} Z^{7/3} により算出される。

実験結果

リサーチクエスチョン

• RQ1 0 ≤ x ≤ 50 の範囲で、数値解をより正確に再現できる新しい解析的試行関数を構築可能か？
• RQ2 提案された解は、イオン化エネルギー計算において極めて重要な x=0 における微分係数を、従来の試行関数よりもより正確に再現できるか？
• RQ3 重い原子（Z ≥ 92）の全イオン化エネルギーについて、新しい解の精度はハートリー・フォック結果と比較してどの程度か？
• RQ4 原子番号が増加するに従い、従来の近似とは異なり、新しい試行関数は精度を維持または向上させられるか？
• RQ5 既存の試行関数が、解およびその微分を十分な精度で再現できない理由は何か？また、新しい形式はこれらの制限をどのように克服するか？

主な発見

• 提案された解は、0 ≤ x ≤ 10 の範囲で47点において平均誤差0.28%を達成し、2パラメータ関数（式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ハートリーを予測し、ハートリー・フォック値の44,042ハートリーからわずか0.8%の上回りにとどまり、超重元素に対しても高い精度を示す。
• ポテンシャルの形状と原点における微分係数の両方において高い正確性を示すため、重い元素の原子的性質の計算において優れた性能を発揮する。