RPA as a Hessian Closure: Effective Functionals and Source-Variable Duality Across DFT, LR-TDDFT, 1RDMFT, and MBPT

本論文は、ランダム位相近似（RPA）を共通のソース変数階層内におけるヘッセ行列閉包近似として定義する統一的な変分フレームワークを提案し、それによって密度汎関数理論、線形応答時間依存密度汎関数理論、一粒子既約密度行列汎関数理論、および多体摂動論の間の首尾一貫した理論的関連性を確立するものである。

要約

大きなアイデア：RPAは「簡略化された地図」である

想像してみてください。あなたは巨大で複雑な都市（量子物理学の世界）をナビゲートしようとしています。手元には、舗装のひび割れ一つ、木の一本、さらには人々の動きまでをも完璧に描き出した、1:1スケールの精密な地図があります。これが**「厳密な理論（Exact Theory）」**です。正確ではありますが、あまりにも詳細すぎるため、素早い計算を行ったり、全体像を把握したりするには不向きです。

この論文は、**RPA（ランダム位相近似）とは、特定のツールや特定の公式、あるいは特定の種類の地図のことではないと主張しています。むしろ、RPAとは「簡略化の手法」**なのです。それは、その完璧で圧倒的な地図から、主要な道路だけを残して細部を無視することで、使い勝手の良い簡略版を作成するための「ルール」なのです。

著者であるNan Sheng氏は、この簡略化のルールは、都市を上空から眺める（密度）、時間の経過とともに変化する様子を見る（時間依存）、3Dモデルを見る（簡約密度行列）、あるいは都市の交通の全履歴を見る（グリーン関数）といった、どの視点においても同様に機能すると主張しています。

コアとなる概念：「ヘシアン」は硬さの測定器

この簡略化がどのように機能するかを理解するために、論文では**「ヘシアン（Hessian）」**と呼ばれる数学的概念を導入しています。

• 比喩： 都市が巨大で柔軟なトランポリンでできていると想像してください。ヘシアンとは、あらゆる地点におけるトランポリンの「硬さ」や「弾力性」を測る指標です。
  • もしあなたがトランポリンを押し下げた（力を加えた）とき、ヘシアンはそれがどれだけ跳ね返るか（応答）を正確に教えてくれます。
  • **厳密なヘシアン（Exact Hessian）**には、布地、スプリング、風、そして飛び跳ねる人々の重さなど、あらゆる微細な相互作用が含まれています。これは完璧な硬さの測定器です。

論文によれば、RPAとは、「どの部分の硬さを残し、どの部分を捨てるか」を決定する行為なのです。

都市を見る4つの方法（4つのレベル）

論文では、この「簡略化のルール」を、システムを記述する4つの異なる方法に適用できることを示しています。これらは、同じ物理問題を見ている4つの異なるカメラ、あるいはレンズだと考えてください。

1. 静的密度（「スナップショット」）：

   • 何を見るか： ある特定の瞬間における群衆の密度。今、人々はどこに立っているのか？
   • 簡略化： 主要な群衆の圧力（「ハートリー（Hartree）」項）は残しますが、人々が互いに囁き合うような複雑な仕組み（「交換相関（exchange-correlation）」項）は無視します。
   • 結果： 群衆密度のシンプルな地図。

2. 動的密度（「ビデオ」）：

   • 何を見るか： 時間とともに変化する群衆の密度。突然の出来事に対して、群衆はどう動き、反応するのか？
   • 簡略化： 主要な群衆の圧力は残しますが、複雑で時間差のある「囁き声」は無視します。
   • 結果： 実物よりも計算しやすい、群衆の動きのビデオ。

3. 等時間二局所的（「3Dモデル」）：

   • 何を見るか： 単に人々がどこにいるかだけでなく、同じ瞬間に隣人とどのように繋がっているか。空間的に詳細なモデルです。
   • 簡略化： 主要な圧力と、隣人との直接的な「手繋ぎ（交換）」は残しますが、複雑で間接的なソーシャルネットワークは無視します。
   • 結果： 管理可能な範囲内に収まった、詳細な3Dモデル。

4. 時空二局所的（「フル・シミュレーション」）：

   • 何を見るか： 最も完全な視点。空間と時間の両方において、あらゆる個人、その繋がり、そして動きを同時に追跡します。これは「グリーン関数」のレベルです。
   • 簡略化： 主要な圧力と直接的な相互作用は残しますが、複雑で既約な背景ノイズは切り捨てます。
   • 結果： 実行可能な程度にちょうど簡略化された、最も強力なシミュレーション。

決定的な発見：地図は必ずしも一致しない

ここが、この論文の主張の中で最も重要な部分です。

通常、科学者はこう考えるかもしれません。「もし『スナップショット（レベル1）』を簡略化してから『ビデオ（レベル2）』に変換したなら、それは『フル・シミュレーション（レベル4）』を簡略化してから『ビデオ』に変換したものと同じ結果になるはずだ」と。

しかし、論文はこう言っています。「いいえ、それは正しくありません」

• 比喩： あなたが都市の高解像度写真を持っていると想像してください。
  • 経路A： 写真をぼかしてシンプルにしてから、それをアニメーション化しようとする。
  • 経路B： 高解像度の写真をまずアニメーション化してから、そのビデオをぼかす。
  • 結果： 最終的な「ぼやけたビデオ」は、どちらの手順を踏んだかによって、見た目が変わってしまうのです！

論文は、「RPAによる簡略化」は**「どのカメラ（変数）から始めたか」**に依存することを証明しています。

• 「静的密度」のカメラから得られる「RPA」は、たとえ同じ物理現象を記述しようとしていても、「フル・シミュレーション」のカメラから得られる「RPA」とは、数学的に異なるオブジェクトです。
• これらは同じ概念の「並行した実現形態」ではありますが、互換性はありません。単に入れ替えることはできず、特定の作業を行うために、適切な方を選択しなければならないのです。

論文の主張のまとめ

1. RPAは「ヘシアン・クロージャ（Hessian Closure）」である： それは、主要な相互作用を残し、複雑で既約な残りカスを捨てることで、システムの「硬さ（応答）」を簡略化する特定の手法です。
2. それはあらゆる場所で機能する： この論理は、単純な密度、時間依存密度、あるいは複雑な量子シミュレーションなど、どのような視点においても適用可能です。
3. 文脈が重要である： 得られる具体的な結果は、「どのようにシステムを見ているか」によって決まります。「密度」の計算から得られる「RPA」は、構造的に「グリーン関数」の計算から得られる「RPA」とは異なります。彼らは双子ではなく、従兄弟のような関係なのです。

この論文は、新しい応用例や臨床的な用途を提示するものではありません。既存の理論の理解の仕方を再構成し、それらがすべて共通の「簡略化エンジン（ヘシアン・クロージャ）」を共有しながらも、出発点に応じて異なる結果を生み出すことを示しているのです。

技術概要

技術要約：ヘッシアン・クロージャとしてのRPA

問題提起
ランダム位相近似（RPA）は、密度汎関数理論（DFT）、線形応答時間依存密度汎関数理論（LR-TDDFT）、一粒子簡約密度行列汎関数理論（1RDMFT）、およびグリーン関数多体摂動論（MBPT）など、多くの多体物理学において遍在する概念である。しかし、「RPA」という用語は、環状図の再総和（ring-diagram resummations）、密度応答近似、あるいは小振幅平均場理論といった、一連の異なる構成を記述するために用いられることが多い。これらには統一された形式的な定義が存在しない。この単一の形式的言語の欠如により、3つの重要な構成要素間の構造的な区別が不明瞭になっている。すなわち、基本変数の選択、その変数に関連する厳密な応答理論、そして応答理論を閉じる特定の近似である。その結果、異なる分野におけるRPAの定式化が、同一の基礎となる変分対象を表しているのか、あるいは単に関連はあるが異なる近似を表しているのかが不明確となっている。

手法
本論文は、ソース―変数双対性（source–variable duality）とヘッシアン解析に基づく統一的な変分フレームワークを提案する。その手法は以下のステップで進行する。

1. ソース―変数双対性： 著者らは、物理系が共役な一対の変数、すなわち基本変数 $q$（例：密度、グリーン関数）とソース $J$（例：ポテンシャル）によって記述される一般的な枠組みを確立する。ソース側にはエネルギーに類する汎関数 $E[J]$ が定義され、変数側にはルジャンドル変換（双対変分構成）を通じて対応する有効汎関数 $\Gamma[q]$ が定義される。
2. 逆ヘッシアンとしての厳密な線形応答： 厳密な線形応答は、有効汎関数 $\Gamma[q]$ のヘッシアンによって厳密に制御されることを本論文は示している。具体的には、厳密な応答演算子 $\chi$ は、ヘッシアンの逆（符号を除いて）である：$\chi = -(\delta^2 \Gamma / \delta q^2)^{-1}$。
3. ヘッシアンの分解： 厳密なヘッシアンは、参照項（$\Gamma_{ref}$）、明示的に保持される二次の相互作用カーネル（$K_{int}$）、および既約な残差（$\Phi$）の3つの部分に分解される。
   $$\frac{\delta^2 \Gamma}{\delta q^2} = \frac{\delta^2 \Gamma_{ref}}{\delta q^2} + K_{int} + \frac{\delta^2 \Phi}{\delta q^2}$$
4. RPAの定義： RPAは、ヘッシアンのレベルにおいて既約な残差（$\Phi$）を破棄し、参照項と明示的な相互作用カーネルのみを保持することによって得られる**クロージャ近似（closure approximation）**として形式的に定義される。
   $$\frac{\delta^2 \Gamma_{RPA}}{\delta q^2} := \frac{\delta^2 \Gamma_{ref}}{\delta q^2} + K_{int}$$

著者らは、密度の記述に対する2つの独立した拡張によって整理された、4つの具体的な変分レベルに対してこの抽象的な定義を適用している。

• 静的密度レベル (DFT): ソースは局所的なスカラーポテンシャルであり、変数は静的な密度である。
• 動的密度レベル (LR-TDDFT): ソースは時間依存の局所ポテンシャルであり、変数は時間依存の密度である。
• 等時刻二局所レベル (1RDMFT): ソースは非局所的な一粒子ポテンシャルであり、変数は等時刻の一粒子簡約密度行列（1RDM）である。
• 時空二局域レベル (MBPT): ソースは空間および時間において二局所的であり、変数は一粒子グリーン関数である。

主要な貢献と結果

• 統一的定義： 本論文は、RPAを特定の図式的再総和や方程式の運動としてではなく、有効汎関数の厳密なヘッシアンの特定の切断（truncation）として、厳密な定義を提供している。これにより、異なる理論の背後にある共通の変分構造を孤立させている。
• 階層構造： 本研究は、DFT、LR-TDDFT、1RDMFT、およびMBPTを接続する二方向の階層をマッピングしている。
  • 静的DFTからLR-TDDFTへの移動は、動的な拡張（時間依存性の追加）を表す。
  • 静的DFTから1RDMFTへの移動は、空間的な非局所的拡張（局所密度から二局所的な1RDMへの移行）を表す。
  • MBPTはこれら両方の拡張を組み合わせたものである。
• 射影の非可換性： 重要な結果として、異なるレベルにおけるRPAクロージャは、射影の下で可換ではないことが示された。カーネルを射影によって縮小することは、その逆（ヘッシアン）を射影することとは異なる。したがって、「DFT-RPA」、「LR-TDDFT-RPA」、「1RDMFT-RPA」、および「MBPT-RPA」は、異なる変数とチャネルに適用された同一のヘッシアン原理の並行した実現であり、単に異なる表記法で書かれた同一の近似ではない。
• 具体的な実現：
  • DFTにおいては、このクロージャは交換相関カーネル（$f_{xc} \approx 0$）を破棄し、非相互作用応答とクーロンカーネルのみを保持することに対応する。
  • LR-TDDFTにおいては、このクロージャは周波数依存の交換相関カーネルを破棄し、断熱接続揺動散逸（ACFD）公式で使用される標準的な動的RPA応答を与える。
  • 1RDMFTにおいては、直接RPAクロージャは、密度チャネルに投影されたハートリー項のみを保持する。また、本論文は、全二局所空間に作用するハートリー・フォック交換ヘッシアンを保持する「xRPA」（交換を含む）クロージャも定義している。
  • MBPTにおいては、このクロージャは、二粒子ヘッシアンの粒子・ホールチャネルにおける裸のクーロン相互作用を保持し、既約な二粒子カーネル（$K_{rem} \approx 0$）を破棄することに対応する。これにより、標準的なベテ・サルピーター方程式の形式が回収される。

意義と主張
本論文の主な貢献は、新しい計算式や応用を導入することではなく、概念的な明晰さを提供することであると著者らは主張している。RPAを「ヘッシアン・クロージャ」として枠付けることで、著者らは以下のことを達成している。

1. 近似の分解： 変数の選択、厳密な応答理論、およびクロージャ近似を分離し、「RPA」がヘッシアンの特定の切断を指しており、理論全体を指すのではないことを明確にした。
2. 関係の明確化： DFT、LR-TDDFT、1RDMFT、およびMBPTの間の構造的な関係を明らかにし、それらをソース―変数双対性の異なる点として示した。
3. 関連手法の再文脈化： GW、$G_0W_0$、QRPAなどの関連する近似をこの階層内に位置づけている。例えば、GWはRPAの定義そのものではなく、RPA型のヘッシアン・クロージャから導出された遮蔽相互作用に基づいた、ダウンストリームの自己エネルギー近似として記述されている。
4. 限定的な範囲： 著者らは、特定の分野における標準的な定式化を置き換えたり、その応用をレビューしたりすることを目的としていないことを明示している。代わりに、それらの背後にある共通の変分内容を述べることを目的としている。また、非微分的な汎関数の凸解析的取り扱いや、1RDMおよびグリーン関数理論における有用な許容領域の決定に関する未解決の問題についても認めている。

要約すると、本論文は、RPAとは根本的に有効汎関数の厳密なヘッシアンのクロージャであり、その様々な形態は、基本変数と応答チャネルの選択によって決定される、この原理の異なる実現であると論じている。