[論文レビュー] Semiclassics: The hidden theory behind the success of DFT

Pavel Okun, Kieron Burke|arXiv (Cornell University)|May 10, 2021

Advanced Chemical Physics Studies被引用数 5

ひとこと要約

本論文は、非相対論的電子系における高Z（大きな原子核電荷）極限——特に半古典的極限——が、密度汎関数理論（DFT）の成功の根本的要因であると特定する。この極限における全エネルギー、運動エネルギー、交換エネルギー、相関エネルギーの挙動を分析することで、著者らは局所密度近似（LDA）と一般化勾配近似（GGA）が、WKBおよびオイラー＝マクローリン和公式の技法を用いて導かれる主要補正項とともに、次第に正確になることを示している。主たる貢献は、DFT近似と半古典的展開を結びつける統一的理論枠組みを提供することであり、原子、分子、固体にわたり、その経験的成功を説明するものである。

ABSTRACT

We argue that the success of DFT can be understood in terms of a semiclassical expansion around a very specific limit. This limit was identified long ago by Lieb and Simon for the total electronic energy of a system. This is a universal limit of all electronic structure: atoms, molecules, and solids. For the total energy, Thomas-Fermi theory becomes relatively exact in the limit. The limit can also be studied for much simpler model systems, including non-interacting fermions in a one-dimensional well, where the WKB approximation applies for individual eigenvalues and eigenfunctions. Summation techniques lead to energies and densities that are functionals of the potential. We consider several examples in one dimension (fermions in a box, in a harmonic well, in a linear half-well, and in the Pöschl-Teller well. The effects of higher dimension are also illustrated with the three-dimensional harmonic well and the Bohr atom, non-interacting fermions in a Coulomb well. Modern density functional calculations use the Kohn-Sham scheme almost exclusively. The same semiclassical limit can be studied for the Kohn-Sham kinetic energy, for the exchange energy, and for the correlation energy. For all three, the local density approximation appears to become relatively exact in this limit. Recent work, both analytic and numerical, explores how this limit is approached, in an effort to deduce the leading corrections to the local approximation. A simple scheme, using the Euler-Maclaurin summation formula, is the result of many different attempts at this problem. In very simple cases, the correction formulas are much more accurate than standard density functionals. Several functionals are already in widespread use in both chemistry and materials that incorporate these limits, and prospects for the future are discussed.

研究の動機と目的

• 多数体問題の複雑さにもかかわらず、DFTが多様な系にわたりなぜこれほどうまく機能するのかを理解すること。
• DFT近似が次第に正確になる普遍的な極限——特に高い原子核電荷（Z）極限——を特定すること。
• 半古典的技法を用いて、運動エネルギー、交換エネルギー、相関エネルギーに対する局所密度近似（LDA）の主要補正項を導出すること。
• 単純な1次元モデル系（例：井戸型ポテンシャル内のフェルミ粒子、調和井戸）の挙動を、現実の3次元原子や分子の挙動と結びつけること。
• 現代のDFT汎関数の理論的基盤を提供し、それらの正確性が普遍的な半古典的極限への収束に起因することを示すこと。

提案手法

• Lieb-Simon（LS）極限を分析し、原子核電荷Zと電子数が特定の方法でスケーリングされる状況で、全エネルギーに関してトマス＝フェルミ理論が正確になる普遍的領域を同定すること。
• 1次元モデル系における1粒子固有値および固有関数にWKB近似を適用し、ポテンシャルのエネルギーおよび密度汎関数を導出すること。
• オイラー＝マクローリン和公式を用いて、局所密度近似に対する補正項を体系的に導出し、標準的な汎関数を上回る精度を向上させること。
• 3次元系への拡張——特に3次元調和井戸およびボーア原子——を実施し、1次元の結果が高次元でも有効であることを検証すること。
• 勾配展開および端効果補正（境界や転送点由来）を導出し、半古典的展開における非一様性を説明すること。
• 境界条件および尖鋭特異点（cusp singularities）の役割を理解するために、マスロフ指数による系の分類を実施すること。

Figure 1 : Accurate radial density of Xe (blue) and its TF approximation (red), in atomic units.

実験結果

リサーチクエスチョン

• RQ1 Kohn-Sham DFTにおける局所密度近似（LDA）および一般化勾配近似（GGA）が、多様な系にわたりなぜこれほど高い精度を示すのか。
• RQ2 高Z極限（Z → ∞）が、DFT近似が正確になる普遍的な漸近的領域としてどのように機能するのか。
• RQ3 この極限における運動エネルギー、交換エネルギー、相関エネルギーに対するLDAの主要補正項は何か。
• RQ4 WKBおよびオイラー＝マクローリン和公式といった半古典的技法を用いて、1次元モデル系におけるポテンシャルの正確な汎関数をどのように導出できるか。
• RQ5 なぜ交換汎関数は一般的に均一電子ガス極限を満たすが、相関汎関数はそうではないのか。この違いが重い原子における性能にどのように影響を与えるか。

主な発見

• Z → ∞ かつ電子数を適切にスケーリングするLieb-Simon極限では、トマス＝フェルミ理論が全電子エネルギーに関して相対的に正確になり、主要補正はスコット補正である。
• 中性原子では、全エネルギーの展開は E(Z) = -c₀Z⁷ᐟ³ + ½Z² - c₂Z⁵ᐟ³ + ... の形を取り、c₀ ≈ 0.768745 および c₂ ≈ 0.269900 は半古典的解析から導かれる基本定数である。
• フェルミ粒子が井戸型ポテンシャルや調和井戸に閉じ込められた1次元系では、オイラー＝マクローリン公式により得られる補正式は、標準的なDFT汎関数よりも顕著に高い精度を示す。
• Kohn-Sham形式における運動エネルギー密度は、Zの増加に従い、局所極限に急速に近づくが、相関エネルギー補正はZ ln Zに比例するのみであり、支配的になるには物理的に不自然に大きなZが必要である。
• 勾配展開は有限系における境界効果や転送点効果を捉えられず、標準的なGGA汎関数に含まれない追加の端効果補正が必要となる。
• マスロフ指数による分類により、硬い壁（例：3次元径方向問題における原点）や尖鋭特異点（例：クーロンポテンシャル）を有する系は、半古典的取り扱いが異なり、汎関数の精度に影響を与えることが明らかになった。

Figure 2 : The exact noble gas energies (black circles) compared with the expansion in Eq. ( 1 ): TF (blue), with first order correction (red), second order (orange).