Kohn-Sham反演の数学的保証付き方法

基本情報

• 論文ID: 2409.04372
• タイトル: Kohn-Sham inversion with mathematical guarantees
• 著者: Michael F. Herbst (EPFL)、Vebjørn H. Bakkestuen (オスロメトロポリタン大学)、Andre Laestadius (オスロメトロポリタン大学・オスロ大学)
• 分類: physics.chem-ph、math-ph、math.MP、quant-ph
• 発表時期: 2024年9月（arXiv v3: 2025年5月5日）
• 論文リンク: https://arxiv.org/abs/2409.04372

摘要

本論文は、精密なMoreau-Yosida正則化法を用いて周期系の交換相関ポテンシャル(exchange-correlation potential)を取得する。本研究は、厳密な数学原理と効率的な数値実装の間の深い関連性を明らかにし、Moreau-Yosidaベースの反演法を実際の物理系に初めて適用したものである。著者らは数学的に厳密な反演アルゴリズムを開発し、代表的なバルク材料(体シリコン、砒化ガリウム、塩化カリウム)上で検証した。本アルゴリズムは厳密な誤差界を構築し数値検証を行うことができ、Kohn-Sham反演法の分析に新しい道を開き、近似汎関数の数学的手法の発展を促進する可能性がある。

研究背景と動機

核心問題

密度汎関数理論(DFT)は化学、材料科学および固体物理学における不可欠なツールである。Kohn-Sham(KS)形式では、DFTのすべての未知量は交換相関(xc)汎関数に集中しており、近似が必要である。DFTは原理的には正確であるが、KS-DFTは特定の物理シナリオにおいて依然として課題に直面している。特に:

1. 分数電荷問題: 分数電子電荷を伴うプロセス(解離または電荷移動励起など)を正確に記述することが困難
2. バンドギャップ問題: 半導体バンドギャップが系統的に過小評価される

研究の重要性

1. 理論的欠陥: 精密な普遍密度汎関数と一般的に使用される近似の間に数学的理解が欠けており、新しいより良い汎関数を厳密に構築することが困難
2. 反演問題: KS反演(与えられた基底状態密度から精密なxcポテンシャルを求める)は正方向のKS-DFT問題ほど十分に研究されていない
3. 数学的関連性: 初期の研究は分数電荷およびバンドギャップ問題とKS-DFTの非微分可能性の間の深い関連性を確立している

既存方法の限界

• 既存のKS反演法(van Leeuwen-Baerends、Zhao-Morrison-Parr、Wu-Yangなど)は厳密な数学的保証を欠いている
• 堅牢で効率的な数値スキームが欠けている
• 厳密な誤差界と収束性分析がない

研究動機

最近の理論的結果49は、Moreau-Yosida(MY)正則化形式の数学的極限を通じてxcポテンシャルを得ることができることを示唆している。MY正則化は精密な普遍汎関数の非微分可能性問題に対処し、これは上述の物理的問題と密接に関連している。本論文は初めてこの理論的枠組みを実際の物理系に適用する。

核心的貢献

1. 初回実装: 実際の物理系上でMoreau-Yosida枠組みに基づくKS反演を初めて実装
2. 数学的厳密性: 厳密な数学的保証を持つ反演アルゴリズムを開発し、精密な反演公式(方程式10)を確立
3. 誤差界: 逆KS問題に対する初めての厳密な誤差界(方程式14、16、17)を確立し、数値検証
4. 非拡張性証明: 近接写像の(確実な)非拡張性質(方程式12)を証明
5. 実用的応用: 3つの代表的なバルク材料(Si、GaAs、KCl)上でアルゴリズムを成功裏に検証
6. 理論-実践の橋: 厳密な数学理論と実際の数値実装の間の関連性を確立

方法の詳細説明

タスク定義

入力: 精密な基底状態密度 $\rho_{gs}$(実験データ、全配置相互作用、結合クラスタまたは量子モンテカルロ計算から得られる可能性がある)

出力: 対応する交換相関ポテンシャル $v_{xc}$。補助的な非相互作用系において当該密度を再現できるようにする

制約: 密度が非相互作用v-表現可能であると仮定(すなわち、ある種のポテンシャルが存在して $\rho_{gs}$ も非相互作用基底状態密度である)

数学的枠組み

1. 関数空間の選択

• 密度空間: $D = H^{-1}_{per}(\Omega, \mathbb{C})$(周期Sobolev空間)
• ポテンシャル空間: $V = H^1_{per}(\Omega, \mathbb{C})$(Dの双対空間)
• ノルム定義: $\|u\|^2_{H^s_{per}} = \sum_G (1 + |G|^2)^s |\hat{u}_G|^2$ ここでGは逆格子ベクトル

2. 双対写像

双対写像 $J: D \to V$ は以下のように定義される: $J(\rho) = \{v \in V : \|v\|^2_V = \|\rho\|^2_D = \langle v, \rho \rangle\}$

選択された関数空間では、双対写像は明示的な形式を持つ: $J[\rho](r) = (\Phi * \rho)(r) = \int_{\mathbb{R}^3} \frac{\rho(r')}{4\pi|r-r'|} e^{-|r-r'|} d^3r'$ ここで $\Phi(r) = e^{-|r|}/(4\pi|r|)$ はYukawaポテンシャルであり、この形式は数値的に処理しやすい。

3. Moreau-Yosida正則化

導出密度汎関数は以下のように定義される: $\mathcal{F}(\rho) = T(\rho) + E_H(\rho) + \int_\Omega v_{ext}\rho$ ここで $T(\rho)$ は運動エネルギー汎関数、$E_H(\rho)$ はHartree寄与である。

重要な最適化問題: $\mathcal{E}(\rho; \rho_{gs}) = \mathcal{F}(\rho) + \frac{1}{2\varepsilon}\|\rho - \rho_{gs}\|^2_D$

この汎関数を最小化して近接密度 $\rho^\varepsilon_{gs} = \arg\min_\rho \mathcal{E}(\rho; \rho_{gs})$ を得る

4. 反演公式

交換相関ポテンシャルは以下の極限を通じて得られる: $v_{xc}(r) = \lim_{\varepsilon \to 0^+} \frac{1}{\varepsilon} \int_{\mathbb{R}^3} \frac{\rho^\varepsilon_{gs}(r') - \rho_{gs}(r')}{4\pi|r-r'|} e^{-|r-r'|} d^3r'$

これは本論文の核心的な理論的結果であり、近接密度からxcポテンシャルへの明示的な計算公式を提供する。

モデルアーキテクチャ

数値実装戦略

1. 軌道パラメータ化: 正規直交軌道 $\Phi = (\psi_1, \ldots, \psi_{N_b})$ を使用して密度をパラメータ化: $\rho_\Phi(r) = 2\sum_{i=1}^{N_b} |\psi_i(r)|^2$
2. エネルギー表現(方程式15): $\mathcal{E}(\Phi, \rho_{gs}) = \sum_{i=1}^{N_b} \int_\Omega |\nabla\psi_i|^2 + E_H(\rho_\Phi) + \int_\Omega v_{ext}\rho_\Phi + \frac{1}{2\varepsilon}\|\rho_\Phi - \rho_{gs}\|^2_D$
3. 最適化方法:
   • BFGSベースの準ニュートン法を使用
   • Stiefel多様体の幾何学的構造に適応(軌道の正規直交性を保持)
   • 停止基準: オプティマイザが機械精度に達するか、$\rho^\varepsilon_{gs}$ の反復変化が0.01ε未満
4. ε列: 1から約$10^{-7}$の範囲の指数関数的に減少する列を採用

技術的革新点

1. 非拡張性証明

近接写像 $\rho \mapsto \rho^\varepsilon$ が(確実な)非拡張作用素であることを証明: $\|\rho^\varepsilon - \tilde{\rho}^\varepsilon\|_D \leq \|\rho - \tilde{\rho}\|_D$

証明の考え方:

• $-\frac{1}{\varepsilon}J(\rho^\varepsilon - \rho) \in \partial\mathcal{F}(\rho^\varepsilon)$ を利用
• 部分微分の極大単調性を通じて
• Hölder不等式を適用

2. 誤差界理論

比率 $Q_\varepsilon(\Delta\rho) := \frac{\|\rho^\varepsilon_{gs} - \tilde{\rho}^\varepsilon_{gs}\|_D}{\|\Delta\rho\|_D} \leq 1$ を定義

主要な誤差界(方程式14): $\|v^\varepsilon_{xc} - \tilde{v}^\varepsilon_{xc}\|_V \leq \frac{1 + Q_\varepsilon(\Delta\rho)}{\varepsilon}\|\Delta\rho\|_D$

精密な界(方程式16): $\left\|v^\varepsilon_{xc} - \tilde{v}^\varepsilon_{xc} - \frac{1}{\varepsilon}J(\Delta\rho)\right\|_V \leq \frac{Q_\varepsilon(\Delta\rho)}{\varepsilon}\|\Delta\rho\|_D$

比率 $R_\varepsilon$ と $S_\varepsilon$ を導入し、以下を証明(方程式17): $0 \leq 1 - Q_\varepsilon(\Delta\rho) \leq R_\varepsilon(\Delta\rho) \leq 1 + Q_\varepsilon(\Delta\rho) \leq 2$

3. 従来の方法との違い

• 従来の方法: 厳密な数学的保証を欠き、通常は発見的最適化に基づく
• 本論文の方法:
  • 凸分析とBanach空間理論に基づく
  • 収束性保証を提供($\rho^\varepsilon_{gs} \to \rho_{gs}$ when $\varepsilon \to 0^+$)
  • 計算可能な誤差界
  • 汎関数の非微分可能性に対処

実験設定

材料系

3つの代表的なバルク材料を研究:

1. 体シリコン(Si): 典型的な半導体
2. 砒化ガリウム(GaAs): 化合物半導体
3. 塩化カリウム(KCl): イオン結晶

計算の詳細

正方向計算(参照)

• xc汎関数: PBE汎関数
• 疑似ポテンシャル: PBE pseudodojo標準疑似ポテンシャル(非線形核修正を含む)
• k点間隔: 最大0.12 Å$^{-1}$
• 運動エネルギーカットオフ: 推奨値の約2倍(高精度を確保)
• ソフトウェア: Density-Functional ToolKit (DFTK)

反演計算

• 同じ疑似ポテンシャル近似を使用(Kleiman-Bylander非局所ポテンシャル項を含む)
• ε列: 1から約$10^{-7}$の指数関数的に減少する列
• 最適化停止基準: 機械精度または $\Delta\rho^\varepsilon_{gs} < 0.01\varepsilon$

誤差導入方法

誤差界をテストするため、Fourier基底切断を通じて制御された摂動 $\Delta\rho$ を導入:

• 異なるカットオフエネルギー $E_{cut}$(15、25、35、45 Ha)
• $E_{cut} = 45$ Haが無摂動参照
• 対応する $\|\Delta\rho\|_D$ を計算

実験設定の説明

著者らは「逆犯罪」設定(正方向と反演が同じモデルと離散化基を使用)を採用したことを認めているが、以下の理由から強調している:

1. 数学理論の厳密性を検証
2. 反演で得られた密度とポテンシャルを参照値と直接比較
3. 将来の研究では他の高精度方法からの参照密度を使用

実験結果

主要な結果: 精密な反演(ノイズなし)

1. 体シリコン(図2)

• ポテンシャルの回復: 結晶の高対称経路(O → (001) → O' → (110) → O'' → (111) → O)に沿ってポテンシャルをプロット
• 収束性能:
  • $\varepsilon \sim 10^{-6}$: 相対誤差は10%未満
  • εが1桁減少: 誤差がさらに1桁減少
• 空間特性: ポテンシャルの最も鋭い特性の近くで、点ごとの収束がより遅く、相対誤差がより大きい

2. 砒化ガリウム(図3)

• 類似の経路に沿ってプロット(Ga-Ga結合の間から開始)
• シリコンと比較して、同じε値で絶対相対誤差がやや大きい
• 全体的に参照ポテンシャルが正確に回復される

3. 塩化カリウム(図4)

• カリウム(K)原子から開始する経路
• GaAsと類似の誤差特性
• 3つの材料すべてで参照ポテンシャルが成功裏に回復される

重要な発見: 追加のノイズがない場合($\Delta\rho = 0$)、アルゴリズムは3つの材料すべてのxcポテンシャルを正確に回復でき、方法の有効性を検証している。

ノイズを含む反演と誤差推定

1. 収束挙動(図5)

異なる基底切断によって導入された摂動の収束への影響を研究:

• 重要な観察: $\varepsilon > \|\Delta\rho\|_{L^2_{per}}$ である限り、ポテンシャルの収束特性は変わらない
• より小さいεの場合、ポテンシャルは参照(V-ノルム)から逸脱し始める
• 異なるカットオフエネルギー(15、25、35 Ha)は異なる $\|\Delta\rho\|_D$ に対応

2. 非拡張性の検証(図6)

比率 $Q_\varepsilon(\Delta\rho) = \|\rho^\varepsilon_{gs} - \tilde{\rho}^\varepsilon_{gs}\|_D / \|\Delta\rho\|_D$ を計算:

• 理論的界: $0 \leq Q_\varepsilon \leq 1$(近接写像の非拡張性により保証)
• 数値結果:
  • 大きいε値: $Q_\varepsilon \ll 1$
  • $\varepsilon \to 0^+$: $Q_\varepsilon \to 1^-$
  • 理論的予測と完全に一致

3. 誤差界の検証(図7)

比率 $S_\varepsilon$(図7上):

• 定義: $S_\varepsilon(\Delta\rho) := \varepsilon \|v^\varepsilon_{xc} - \tilde{v}^\varepsilon_{xc} - \frac{1}{\varepsilon}J(\Delta\rho)\|_V / \|\Delta\rho\|_D$
• 理論的界: $0 \leq S_\varepsilon \leq Q_\varepsilon \leq 1$
• 数値的性能:
  • 大きいε: $S_\varepsilon$ は小さい(3つのほぼゼロの量の差異により、界からわずかに逸脱)
  • 小さいε: $Q_\varepsilon$ で設定された界と完全に一致
  • $\varepsilon \to 0^+$: $S_\varepsilon \to 1$

比率 $R_\varepsilon$(図7下):

• 定義: $R_\varepsilon(\Delta\rho) := \varepsilon \|v^\varepsilon_{xc} - \tilde{v}^\varepsilon_{xc}\|_V / \|\Delta\rho\|_D$
• 理論的界: $1 - Q_\varepsilon \leq R_\varepsilon \leq 1 + Q_\varepsilon$
• 数値的性能:
  • $Q_\varepsilon$ で設定された界を厳密に遵守
  • 下界 $R_\varepsilon \geq 1 - Q_\varepsilon$ に密接に追従
  • 大きいε値: $R_\varepsilon \approx 1$
  • 小さいε: 下界に近い

数値的課題: $\varepsilon \lesssim 5 \times 10^{-6}$ の場合、問題は数値的に課題となり、2つの量の傾向に小さな振動が現れる。

実験的発見の要約

1. 理論検証: 数値計算は理論的予測の誤差界と非拡張性質と完全に一致
2. 堅牢性: アルゴリズムは密度摂動に対して良好な堅牢性を示す($\varepsilon > \|\Delta\rho\|$ の範囲内)
3. 普遍性の指示: $Q_\varepsilon$ は $\Delta\rho$ に独立した定数で推定される可能性がある(ただしパラメータはεと導出汎関数に依存)
4. 実用性: 3つの異なる種類の材料系への成功した応用
5. 数値精度: $\varepsilon \sim 10^{-7}$ でも安定した計算を維持

関連研究

KS反演法の発展の歴史

初期の研究(1988-1994)

• Aryasetiawan & Stott (1988): 有効ポテンシャル法
• Knorr & Godby (1992): 量子モンテカルロ研究モデル半導体
• Görling (1992): 電子密度からKSポテンシャルと波動関数の決定
• van Leeuwen & Baerends (1994): 正しい漸近挙動を持つxcポテンシャル

最適化方法(2002-2007)

• Wu-Yang (2002、2003): 直接最適化法
• Zhao-Morrison-Parr (1994): ZMP法
• Bulat et al. (2007): 有限基底組における最適化有効ポテンシャル

最近の進展(2013-2024)

• ソフトウェアパッケージ開発:
  • n2v (Shi、Chávez、Wasserman、2022)
  • KS-pies (Nam et al.、2021)
• 固体系への拡張:
  • Aouina et al. (2023): 固体基底状態密度の精密なKS補助系
  • Ravindran et al. (2024): 固体における密度反演局所xcポテンシャル
• 理論分析:
  • Burke組: 密度駆動誤差分析
  • Gould (2023): 「Lieb-response」法

数学理論の基礎

DFTにおけるMoreau-Yosida正則化の応用

• Kvaal et al. (2014): 微分可能だが精密なDFT形式
• Laestadius et al. (2018、2019): Banach空間上の一般化KS反復
• Penz et al. (2019): 有限次元正則化KS反復の保証収束
• Penz、Csirik、Laestadius (2023): MY正則化からの密度-ポテンシャル反演(本論文の直接的な理論的基礎)

本論文の相対的な利点

1. 初回実装: MY枠組みを実際の物理系に初めて適用
2. 厳密な保証: 数学的に厳密な誤差界を提供(前例がない)
3. 理論-実践の結合: 抽象的な数学理論を計算可能な数値スキームに変換
4. 一般的な枠組み: 周期系に適用可能で、より複雑な系への拡張が可能
5. 誤差分析: 既存の方法の発見的な誤差推定を超える

結論と議論

主要な結論

1. 方法の有効性: MY正則化に基づく厳密なKS反演アルゴリズムを成功裏に開発・検証
2. 理論的貢献:
   • 明示的な反演公式(方程式10)を確立
   • 近接写像の非拡張性(方程式12)を証明
   • 初めての厳密な誤差界(方程式14、16、17)を導出
3. 数値検証: 3つの代表的なバルク材料上で理論的予測を検証
4. 橋渡し作用: 数学分析、数値スキーム、物理的近似の間の関連性を確立

限界

1. 理論的側面

• 非局所ポテンシャル: 現在の理論的枠組みはまだ疑似ポテンシャルの非局所効果を含んでいない(数値実装では使用されているが)
• 関数空間の選択: $H^{-1}_{per}$ と $H^1_{per}$ の選択は合理的だが、他の選択がより最適である可能性がある

2. 数値的側面

• 逆犯罪: 正方向と反演が同じモデルを使用し、将来は独立した参照密度の供給源を使用する必要がある
• ε列の最適化: 現在は単純な指数関数的減少列を使用し、より効率的な選択が存在する可能性がある
• 停止基準: 0.01εの発見的基準はさらに最適化される可能性がある
• 計算コスト: 各ε値に対して最適化問題を解く必要がある

3. 適用範囲

• 現在は周期的な絶縁体系に限定
• 3つの比較的単純な材料でのみ検証
• システムサイズは数百の電子に限定

将来の方向

短期的な目標

1. 参照密度の供給源: 平均場理論(半局所DFTを超える)からの密度への応用
2. 非局所ポテンシャル理論: 非局所効果を含むように理論的枠組みを完善
3. 関数空間の最適化: 他の関数空間選択の効果を探索
4. 近似誤差界: $Q_\varepsilon$ が定数である可能性の観察に基づいて、より実用的な誤差推定を開発

長期的な目標

1. 汎関数開発: 厳密な反演スキームを利用して新しい近似汎関数を支援
2. Hohenberg-Kohn写像: 密度-ポテンシャル写像の深い理解
3. 量子埋め込み: 量子埋め込み技術への応用
4. 最適化有効ポテンシャル: 最適化有効ポテンシャル法の改善
5. 複雑な系: より大きく複雑な材料系への拡張

深い評価

利点

1. 方法の革新性(★★★★★)

• 理論的突破: MY正則化理論を実際の物理系に初めて成功裏に適用し、理論と実践の隔たりを埋める
• 数学的厳密性: KS反演分野で前例のない数学的保証を提供
• 誤差界: 計算可能で検証可能な初めての厳密な誤差界を確立
• 非拡張性の利用: 凸分析の非拡張性質を巧妙に利用して誤差理論を構築

2. 実験の十分性(★★★★☆)

• 複数材料の検証: 3つの異なる種類の材料(半導体、イオン結晶)上で検証
• 系統的なテスト:
  • 精密な反演(ノイズなし)
  • ノイズを含む反演
  • 誤差界の検証
  • 収束性分析
• 明確な可視化: 実空間ポテンシャル図と誤差図を通じて結果を直感的に表示
• 定量的分析: 詳細な数値データと比率分析を提供

改善の余地:

• より多くの材料種(金属、強相関系など)を追加できる
• 他の反演法との定量的比較を行える
• 計算効率分析を提供できる

3. 結果の説得力(★★★★★)

• 理論-実験の一致: 数値結果は理論的予測と完全に一致
• 誤差界の厳密性: すべての比率が理論的界内にある
• 収束の明確性: $\varepsilon \to 0$ の収束挙動を明確に示す
• 堅牢性の検証: 密度摂動に対する方法の安定性を証明

4. 執筆の明確性(★★★★★)

• 構造の合理性: 理論→数値実装→実験検証の論理が明確
• 数学的表現: 厳密だが読みやすく、物理的直感を適切に使用
• 図表の質: 高品質なポテンシャル図と誤差分析図
• 再現性: 完全なオープンソースコード(GitHub)とデータ(Zenodo)を提供

不足

1. 方法の限界

• 計算コスト: ε値の列に対して最適化問題を解く必要があり、従来の方法より高くなる可能性がある
• ε選択: ε列を適応的に選択するための理論的指導が欠けている
• 関数空間への依存: 結果は特定の関数空間選択に依存し、その最適性は十分に探索されていない

2. 実験設定

• 逆犯罪: 著者が認める限界で、将来の研究で解決する必要がある
• 材料の多様性: 比較的単純な3つの材料のみでテスト
• ベンチマーク比較: 他の反演法(Wu-Yang、ZMP)との直接的な定量的比較が欠けている

3. 理論の完全性

• 非局所ポテンシャル: 理論的枠組みはまだ実際に使用される非局所疑似ポテンシャルを含んでいない
• 近似界: 現在の誤差界は $Q_\varepsilon$ の計算を必要とし、実際の応用では実行不可能である可能性がある
• 最適性: 提案された方法がある意味で最適であることは証明されていない

影響力

1. 分野への貢献(★★★★★)

• パラダイムシフト: KS反演に厳密な数学的枠組みを導入し、分野の研究方法を変える可能性がある
• 理論的基礎: より信頼性の高い反演スキームの将来の開発のための堅実な基礎を提供
• 誤差分析の新しい方向: 密度-ポテンシャル反演問題の厳密な誤差推定の新しい方向を開く
• 学際的な橋: 関数解析、凸最適化、量子化学を結びつける

2. 実用的価値(★★★☆☆)

• 現在の段階: 主に概念実証であり、直接的な実用性は限定的
• 将来の可能性:
  • 汎関数開発プロセスを改善する可能性
  • 量子埋め込めに対するツールを提供
  • DFTの基礎的問題(バンドギャップ、分数電荷)の理解を支援
• 計算コスト: 日常的な計算に使用するにはさらに最適化が必要

3. 再現性(★★★★★)

• オープンソースコード: Julia実装の完全版(DFTKベース)
• データの公開: 原始データはZenodo上で公開(DOI: 10.5281/zenodo.14894064)
• 詳細なドキュメント: 方法とパラメータの説明が詳細
• ソフトウェアエコシステム: 成熟したDFTKプラットフォームに基づき、拡張が容易

適用シナリオ

理想的な応用シナリオ

1. 汎関数開発: 厳密な誤差制御が必要な近似汎関数の構築
2. ベンチマークテスト: 他の反演法の厳密な参照標準を提供
3. 理論研究: DFTの基礎的問題(非微分可能性、v-表現可能性)の探索
4. 方法学研究: 新しい数値反演技術の開発

現在の制限シナリオ

1. 大規模系: 計算コストが応用を制限する可能性がある
2. 金属系: 現在の実装は絶縁体に限定
3. 強相関系: これらの系でテストされていない
4. リアルタイム応用: 迅速な反演が必要なシナリオには不適切

将来の拡張シナリオ

1. 量子埋め込め: 埋め込め法の核心成分として
2. 機械学習: ML汎関数の高品質な訓練データを提供
3. 不確実性の定量化: 誤差界を利用した不確実性分析
4. マルチスケールシミュレーション: 異なる精度レベル間での情報伝達

参考文献(主要文献)

理論的基礎

1. Penz、Csirik、Laestadius (2023): "Density-potential inversion from Moreau–Yosida regularization"、Electron. Struct. 5、014009 - 本論文の直接的な理論的基礎
2. Penz et al. (2019): "Guaranteed convergence of a regularized Kohn-Sham iteration in finite dimensions"、Phys. Rev. Lett. 123、037401
3. Laestadius et al. (2018): "Generalized Kohn–Sham iteration on Banach spaces"、J. Chem. Phys. 149、164103

DFT基礎

4. Hohenberg & Kohn (1964): "Inhomogeneous electron gas"、Phys. Rev. 136、B864
5. Kohn & Sham (1965): "Self-consistent equations including exchange and correlation effects"、Phys. Rev. 140、A1133
6. Levy (1979): "Universal variational functionals of electron densities"、Proc. Natl. Acad. Sci. USA 76、6062

反演法

7. van Leeuwen & Baerends (1994): "Exchange-correlation potential with correct asymptotic behavior"、Phys. Rev. A 49、2421
8. Wu & Yang (2003): "A direct optimization method for calculating density functionals"、J. Chem. Phys. 118、2498
9. Shi & Wasserman (2021): "Inverse Kohn–Sham Density Functional Theory: Progress and challenges"、J. Phys. Chem. Lett. 12、5308

数値ツール

10. Herbst、Levitt、Cancès (2021): "DFTK: A Julian approach for simulating electrons in solids"、Proceedings of the JuliaCon Conference 3、69

総合評価

推奨読者: DFT理論研究者、量子化学方法論者、計算材料学者、数値解析と凸最適化に関心のある研究者。