Semilocal exchange functionals from the exact-exchange condition for the hydrogen atom: Hydrogenic exactness and recovery of Rydberg-like bound states

技術要約：水素原子の厳密交換条件に基づく半局所交換汎関数

問題提起
標準的な半局所密度汎関数（PBEのようなGGAやSCANのようなmeta-GGA）は、根本的な構造的欠陥を抱えている。すなわち、それらの交換ポテンシャル（$v_x$）は有限系の外領域において指数関数的に減衰し、正確なクーロン的な$-1/r$の裾（テイル）を再現できない。この欠陥は、リدバーグ状態（Rydberg states）の不在、不安定な負イオン、系統的なイオン化エネルギーの過小評価、および自己相互作用誤差（SIE）による過剰な電荷の非局在化といった、よく知られた一連の失敗を引き起こす。モデルポテンシャル（例：LB94、mBJ）や範囲分離ハイブリッドがこれに対処しているが、前者は変分的な一貫性に欠け（親となるエネルギー汎関数が存在しない）、後者は高い計算コストを要する。既存の半局所的な試みであるArmiento-Kümmel（AK13）汎関数は漸近挙動を改善しているが、水素原子に対する厳密な条件から導出されたものではない。

手法
本研究は、1993年にGillとPople（GP93）によって導出された常微分方程式を再検討するものである。この方程式は、GGAが水素原子の正確な交換ポテンシャルを再現するために満たすべき条件を定義している。著者らは、この方程式を逆問題として完全に解くことで、独自の増幅因子、すなわちGP93因子（$F_{GP93}$）を導出する。

手法は以下の3つの層で構成される：

1. 水素的厳密性（Hydrogenic Exactness）： GP93因子は、水素の1s密度に対して厳密となるように構築されている。一様な座標スケーリングにより、これは経験的パラメータなしに、すべての水素様1sイオンに対して厳密な固有値（$\epsilon_{1s} = -Z^2/2$）を与える。
2. 漸近挙動： 大きな勾配領域において、GP93因子は$O[s(\ln s)^{2/3}]$の形の非有界な成長を示す。これはAK13の$s \ln s$の成長とは異なる。この特定の成長により、交換ポテンシャルは水素密度上で正確な$-1/r$の裾を回復する。
3. スイッチング機構： GP93因子の発散的な挙動が多電子系を不安定化させるのを防ぐため、著者らはこれを「スイッチされた」meta-GGA形式の中に組み込んでいる。
   • $\alpha$-スイッチ： 指標 $\alpha = (\tau - \tau^W)/\tau_{unif}$ を用いて、一電子的な領域を検出する。GP93項は $\alpha \approx 0$（一電子軌道領域）においてのみ活性化される。
   • $q$-スイッチ（トリプル・スイッチ）： 原子的な裾（裾が必要な領域）と共有結合の中心（裾が不安定化を招く領域）を区別するために、簡約されたラプラシアン指標 $q = \nabla^2 n / (4(3\pi^2)^{2/3}n^{5/3})$ が導入される。原子的な裾では $q \to +\infty$ となるが、結合中心では $q \lesssim 0$ となる。スイッチ関数 $D(q)$ は、結合中心においてGP93項を抑制しつつ、原子的な裾においてはそれを保持する。
   • 軌道独立性： 運動エネルギー密度 $\tau$ は、Perdew-Constantin（PC07）の構成に基づく軌道自由運動エネルギー密度汎関数（KEDF）に置き換えられ、これにより $\tau^W \le \tau$ が強制される。これにより、全体系が純粋な密度汎関数となり、一般化されたコーン・シャム（gKS）ソルバーを必要としなくなる。

主な結果

• 水素的厳密性： 本汎関数は、水素様イオン（$Z=1$から$5$）に対して、ミリハートリーレベルで正確な固有値$\epsilon_{1s} = -Z^2/2$ を再現する。
• リドバーグ様束縛状態： 一電子的な$s$殻HOMOを持つ系（H、He、および大きなコアを持つ擬ポテンシャルを用いたLi/Na/K）において、本汎関数は一連の束縛された仮想コーン・シャム状態（4〜8準位）を回復する。対照的に、標準的な汎関数（PBE、SCAN、TPPS、M06-L）およびグローバルハイブリッド（B3LYP、PBE0）は、これらの系においてゼロのリドバーグ状態しか束縛しない。
• イオン化エネルギー： テストされた$s$殻系において、最高被占軌道固有値（$-\epsilon_{HOMO}$）は、PBE、SCAN、またはB3LYPと比較して、実験的なイオン化ポテンシャルに著しく近い（平均絶対誤差 $\approx 0.80$ eV）。
• 選択性： 本汎関数は、$\alpha > 0$ となる$p$殻原子（Ne、Ar）に対してGP93項を正しく抑制し、設計通り、ゼロのリドバーグ状態をもたらす。
• 分子の安定性： トリプル・スイッチは、未スイッチのGP93因子が失敗するような共有結合中心（例：$H_2$、$N_2$）における数値的な発散を正常に抑制しつつ、原子系における正しい漸近的裾を維持する。
• 相関との適合性： 本汎関数は、一電子系において消失するLYP相関汎関数と組み合わされており、これにより交換成分の水素的厳密性が保持される。

意義と主張
本論文は、既存のすべての半局所的汎関数（SCANを含む）の定性的失敗である、一電子的な$s$殻HOMOを持つ系における正確な$-1/r$の裾と束縛リドバーグ状態を回復する、非経験的な半局所交換汎関数を構築したと主張している。

著者らは、その主張の範囲を明示的に限定している：

• 成功領域： 本汎関数は、一電子的な性質（リドバーグ状態、漸近挙動、水素的準位、SIE）および一電子的な$s$殻HOMOを持つ系のために設計されている。
• 限界：
  • 熱化学： 本汎関数は原子化エネルギーを対象としていない。水素的厳密性を強制することは、PBEにおいて最適化された誤差補償値からH原子の絶対エネルギーを逸脱させるため、Hを含む熱化学には適さない。
  • $p$殻系： 漸近的な改善は、$\alpha$-スイッチ機構により、$p$殻HOMO（例：Ne、Ar）には拡張されない。
  • 定量的IP： Heについては正確であるが、水素を基準とした較正のため、より重い系におけるイオン化エネルギーには系統的な誤差が含まれる可能性がある。
  • 応答特性： 分極率などの高次微分は、スイッチ境界の感度のため、対象範囲外である。
• 実装： 本論文では、固定密度上での非自己整合的な実装を示している。ガウス基底における完全な自己整合的実装は、ラプラシアン項（$v_{lapl}$）のために不安定になることを注記しているが、ラプラシアンを安定して扱える実空間擬ポテンシャル実装には自然に適していると論じている。

本研究は、単一の非経験的な入力――GP93の水素厳密条件――が、いかにして正確な$-1/r$の裾と水素的準位を同時に生成できるかを示しており、非局所的な厳密交換を必要とせずに、半局所的な計算コストでリドバーグ物理学への道を開くものである。