Moreau-Yosida-based Kohn-Sham Inversion for Periodic Systems

技術要約：周期系におけるMoreau–Yosidaに基づくKohn–Sham反転法

問題提起
密度汎関数理論（DFT）は電子構造計算の礎石であるが、正確な交換相関（xc）汎関数の開発は依然として大きな課題である。主要な困難は、厳密な普遍的密度汎関数と広く用いられている近似との関係に関する数学的な洞察が限られていることにある。これらの近似を改善するための重要なステップは、「逆Kohn–Sham（KS）問題」である。すなわち、既知の基底状態密度から有効なKSポテンシャル（特にxcポテンシャル）を再構成することである。いくつかの反転スキームが存在するものの、周期系（固体材料）における堅牢かつ効率的な数値的手法は依然として未解決の問題である。既存のアプローチは、厳密な関数解析的基礎を欠いていることが多く、周期的な設定における安定性や収束性に課題を抱えている。

手法
著者らは、Moreau–Yosida（MY）正則化に基づいた、周期系における密度・ポテンシャル反転のための数学的に厳密なフレームワークを開発している。この手法は以下のステップで進行する：

1. 関数解析的フレームワーク： 理論は、密度のための周期的な同次ソボレフ空間 $H^{-1}_{\text{per,hom}}(\Omega)$ と、その双対空間であるポテンシャルのための $H^{1}_{\text{per,hom}}(\Omega)$ の中で定式化される。この空間の選択は極めて重要であり、平均場Hartreeエネルギーを自然に扱い、ポテンシャルのゲージ自由度（定数を除いて定義されるもの）を処理することを可能にする。
2. Moreau–Yosida 正則化： 普遍的密度汎関数（具体的には運動エネルギー制約付き探索汎関数）は、MYエンベロープを用いて正則化される。これは、$F(\sigma) + \frac{1}{2\varepsilon}\|\sigma - \rho_{\text{gs}}\|^2_X$ という形式の汎関数の最小化を伴う。ここで $\rho_{\text{gs}}$ はターゲットとなる基底状態密度であり、$\varepsilon > 0$ は正則化パラメータである。
3. 近接写像（Proximal Mapping）： 正則化された汎関数の最小化因子は「近接密度」$\rho_\varepsilon$ である。有効なKSポテンシャルは、以下の極限操作によって回収される：
   $$v_{\text{KS}} = \lim_{\varepsilon \to 0^+} \frac{1}{\varepsilon} J(\rho_\varepsilon - \rho_{\text{gs}})$$
   ここで $J$ は密度空間とポテンシャル空間の間の双対写像である。
4. 理論的保証： 主要な解析的貢献は、運動エネルギー制約付き探索汎関数が選択されたトポロジーにおいて下半連続（l.s.c.）であることを証明した点にある。これにより、近接点の存在と反転スキームの妥当性が保証される。
5. 数値実装： 近接密度は、以下の3つの異なるアプローチを用いて計算される：
   • 直接最小化： オービタルに対するエネルギー汎関数の直接的な最小化（BFGSまたはニュートン法を使用）。
   • 自己無撞着場（SCF）： KS的なEuler-Lagrange方程式を反復的に解く。
   • 近接点アルゴリズム： 非平滑凸最適化のために特別に設計された反復スキーム。

主な貢献
本論文には主に3つの貢献がある：

1. 一般的な反転スキーム： 著者らは、反射的バナッハ空間における密度に対する一般的なMYベースの反転アルゴリズムを定式化した。彼らは、このスキームが汎関数の下半連続性と近接写像の非拡大性にのみ依存し、入力密度の誤差伝播を制限することを証明している。
2. 周期系への適用： フレームワークは周期的な同次ソボレフ空間へと特化されている。著者らは、このトポロジーにおける運動エネルギー汎関数の下半連続性を証明し、これによりZhao-Morrison-Parr (ZMP) 法を周期的な設定において初めて確固たる関数解析的基礎の上に置いた。これにより、劣微分における最小ノルムの要素として厳密なxcポテンシャルを決定することが可能となる。
3. 数値解析： 本論文は、周期的な1次元Gross–Pitaevskii方程式、および3次元バルク材料（シリコン、ガリウム砒素、塩化ナトリウム、塩化カリウム）を用いた反転スキームの包括的な数値調査を提供している。

結果

• 収束性： 数値実験により、近接密度 $\rho_\varepsilon$ が $\varepsilon \to 0$ において基底状態密度 $\rho_{\text{gs}}$ に強収束することが示された。再構成されたポテンシャルは参照となるxcポテンシャルに収束する。
• アルゴリズムの性能： Gross–Pitaevier方程式において、完全なニュートン法を用いた直接最小化（DFTKに実装）は、特に小さな $\varepsilon$ において、SCFおよび準ニュートン法よりも優れた収束性を示した。近接点アルゴリズムはPDEソルバーと同様の収束特性を示したが、より高い計算オーバーヘッドを伴った。
• バルク材料： 本手法は、様々なバルク材料のxcポテンシャルを正常に回収した。結果は、収束率がターゲットポテンシャルの滑らかさに依存することを示しており、鋭い特徴（例：原子核付近）を解像するには、より小さな $\varepsilon$ の値が必要となる。
• ガイディング汎関数： ガイディング汎関数の構成（Hartreeエネルギー、外部ポテンシャル）の影響を調査した。外部ポテンシャルをガイディング汎関数に含めることで、再構成されたポテンシャルの精度が大幅に向上した。
• 安定性： 反転スキームは、摂動を加えた入力密度（フーリエ切断によってシミュレート）を用いてテストされた。結果は、再構成されたポテンシャルの精度が入力の誤差の大きさに制限されることを示しており、正則化アプローチの安定性を裏付けている。

意義
本論文の主要な意義は、周期系におけるKS反転のための数学的に厳密で数値的に安定したフレームワークを提供することにある。周期的な同次ソボレフ空間における運動エネルギー汎関数の下半連続性を確立することにより、著者らは厳密な交換相関ポテンシャルを回収するためのMoreau–Yosida正則化の妥当性を検証している。本研究は、抽象的な凸解析と実用的なDFTアプリケーションの間の溝を埋め、密度汎関数の分析および改善のための堅牢なツールを提供している。さらに、数値結果は、MYベースのアプローチが、特に数値ノイズや周期境界条件に対する堅牢性が不可欠な固体システムにおいて、既存の反転スキームに代わる実行可能な選択肢であることを示唆している。