[論文レビュー] Thomas-Fermi Theory

Élliott H. Lieb|ArXiv.org|Mar 27, 2000

History and advancements in chemistry参考文献 3被引用数 5

ひとこと要約

この論文は、大きな原子および分子の基底状態エネルギーおよび電子密度を、粒子数制約の下で最小化される凸エネルギー汎関数として、Thomas-Fermi (TF)理論を変分近似として定式化する。主な結果は、核電荷 $ Z \to \infty $ の極限において、TF理論が量子力学的エネルギーおよび密度を漸近的に再現することであり、その有効性と限界について、Tellerの定理による分子不安定性に関する厳密な根拠が与えられる。

ABSTRACT

This is a brief review of the Thomas-Fermi ("statistical" or "density functional") theory of atoms and molecules.

研究の動機と目的

• Thomas-Fermi理論を、大きな原子および分子における平均場近似として、厳密な数学的基盤を確立すること。
• 核電荷 $ Z \to \infty $ の極限におけるTF理論の漸近的妥当性を調査すること。
• 特に大 $ Z $ の領域において、TFエネルギー汎関数と真の量子力学的基底状態エネルギー $ E^Q(N) $ の関係を明確にすること。
• 特にTellerの定理による分子不安定性を経由して、TF理論における分子の安定性を分析すること。
• 大 $ Z $ 極限における量子力学的電子密度がTF密度にどのように収束するかを検討すること。

提案手法

• TFエネルギー汎関数 $ \mathcal{E}(\rho) $ を、運動エネルギー、電子-核、電子-電子、核反発の項の和として定式化し、$ \rho \in L^{5/3} $、$ \rho \geq 0 $、$ \int \rho = N $ を満たすものとする。
• $ \mathcal{E}(\rho) $ を最小化してTF基底状態エネルギー $ E^{\rm{TF}}(N) $ を定義し、汎関数の凸性を用いて $ N \leq Z $ の場合に最小解が存在かつ一意であることを保証する。
• Euler-Lagrange方程式としてのThomas-Fermi方程式を導出：$ \gamma \rho^{2/3} = [\Phi(x) - \mu]_+ $、ここで $ \Phi(x) = V(x) - \int |x-y|^{-1} \rho(y) dy $。
• 化学ポテンシャル $ \mu $ の物理的解釈を確立：$ \mu(N) = -\partial E^{\rm{TF}}(N)/\partial N $ であり、$ \mu(Z) = 0 $ であることを示す。
• Lieb-Simonの極限 $ Z \to \infty $、$ N = \lambda Z $、スケーリングされた核位置 $ R_j = Z^{-1/3} R_j^0 $ を用いて、$ \lim_{Z \to \infty} E^{\rm{TF}}(\lambda Z)/E^Q(\lambda Z) = 1 $ を証明する。
• Lieb-ThirringおよびLieb-Oxfordの不等式を適用し、修正されたTF汎関数が真の量子エネルギーの下界を提供することを示す。

実験結果

リサーチクエスチョン

• RQ1 核電荷 $ Z \to \infty $ の極限において、Thomas-Fermiエネルギー汎関数 $ E^{\rm{TF}}(N) $ は、真の量子力学的基底状態エネルギー $ E^Q(N) $ に収束するか？
• RQ2 最小解 $ \rho^{\rm{TF}}_N $ が存在する粒子数 $ N $ の範囲は何か？ これは全核電荷 $ Z $ とどのように関係するか？
• RQ3 分子がTF理論で不安定であると述べるTellerの定理は、大 $ Z $ 極限における真の量子力学的挙動をどのように反映しているか？
• RQ4 大原子におけるTF電子密度の漸近的挙動は何か？ 量子力学的密度と比較するとどうなるか？
• RQ5 von Weizsäcker項およびDirac交換項といった補正が、TF汎関数およびその物理的予測にどのように影響を与えるか？

主な発見

• $ N \leq Z $ の場合、TFエネルギー汎関数 $ \mathcal{E}(\rho) $ は一意な最小解 $ \rho^{\rm{TF}}_N $ を持ち、$ E^{\rm{TF}}(N) $ は厳密に凸かつ単調減少である。$ N > Z $ の場合、$ E^{\rm{TF}}(N) = E^{\rm{TF}}(Z) $ となる。
• 中性原子のTF電子密度は、$ \rho^{\rm{TF}}_{\rm{atom}}(x; N=\lambda Z, Z) = Z^2 \rho^{\rm{TF}}_{\rm{atom}}(Z^{1/3}x; N=\lambda, Z=1) $ とスケーリングされ、サイズが $ Z^{-1/3} $ に縮小することが示される。
• 中性TF密度の $ |x| $ が大きい領域における漸近的挙動は、$ \rho^{\rm{TF}}_{\rm{atom}}(x) \sim \gamma^3 (3/\pi)^3 |x|^{-6} $ であり、$ Z $ に依存せず、$ Z^{-1/3} $ スケールで量子力学的挙動と一致する。
• $ Z \to \infty $ かつ $ N = \lambda Z $ の極限において、比 $ E^{\rm{TF}}(\lambda Z)/E^Q(\lambda Z) \to 1 $ となる。これは、大原子においてTF理論が漸近的に正確であることを証明する。
• Tellerの定理は量子力学的領域でも成立する：$ N = \sum N_j \leq Z $ に対して $ E^{\rm{TF}}(N) > \sum_j E^{\rm{TF}}_{\rm{atom}}(N_j, Z_j) $ であり、核反発 $ U $ のため分子不安定性が生じることを示す。
• 修正された $ \gamma $ を用いたTF汎関数は、真の量子エネルギー $ E^Q(N) $ の下界を提供する。von Weizsäcker項およびDirac交換項といった補正は、大 $ Z $ 極限においても主要な特徴を保つ。