乱雑位相近似

乱雑位相近似（らんざついそうきんじ、, RPA）とは、元々デヴィッド・ボームとによって展開された多体系における基底状態の量子揺らぎ及び励起振動状態(フォノン)を記述するための近似手法。線形応答理論における摂動論的な近似法の一つである。
粒子系（電子ガスなど）が高密度の場合は、乱雑位相近似が妥当な近似であることが分かっている。
同等な近似手法が、多方面（例：GW近似）で利用、応用されている。
第一量子化でのRPA.
N粒子系（N電子系）における密度演算子は次のように与えられる。
ここで位置座標ベクトルformula_2が無秩序であれば、逆格子ベクトルと位置座標ベクトルとの積、formula_3も無秩序（乱雑）なので、formula_4からの寄与がformula_5よりずっと小さいとして無視できる。これを乱雑位相近似(RPA)という。
formula_6においては、formula_3が乱雑なことにより各項の位相も乱雑となり、和の各成分が相殺し合って全体としての寄与が無視できるほど小さくなることによる。勿論、この近似が適用できない場合も多々ある。
概要.
RPAは1952年と1953年にボームとパインズによって初めて導入された。それまで何十年もの間、電子間のミクロな量子力学的相互作用の効果を物質の理論に取り入れようとする試みがあった。ボームとパインズのRPAは、弱く遮蔽されたクーロン相互作用を説明し、電子系における電子の動的な線形応答を記述するために用いられる。
RPAでは、電子は全電位"V"(r)（つまり外部摂動ポテンシャル "V"ext(r) と遮蔽ポテンシャル "V"sc(r) の和）にのみ応答すると仮定される。外部摂動ポテンシャルは単一の周波数 "ω" で振動すると仮定されるので、このモデルに自己無撞着場（SCF）法を適用すると動的誘電関数 "ε"RPA(k, "ω")が得られる。
誘電関数への全電位の寄与は平均化される仮定するため、波数ベクトル k における電位のみが寄与する。これが乱雑位相近似が意味するものである。結果として生じる誘電関数はリンドハード誘電関数とも呼ばれ、電子ガスの多くの性質（プラズモンなど）を正確に予測している。
RPAは自由度を過大に評価していると50年代後半に批判され、その正当化には理論物理学者の多くの労力が費やされた。マレー・ゲルマンとキース・ブルックナーは、高密度の電子ガスのファインマンダイアグラムにおける最低次のチェーンの和からRPAが導かれることを示した。
これらの結果の一貫性はRPAの正当化には重要であり、50年代後半と60年代の理論物理学は大きく発展した。
RPA方程式.
まず第0近似としてハートリー-フォック近似を考える。ハートリー-フォック近似で得られた基底状態には量子揺らぎ効果は含まれてはいない。
そこで、量子揺らぎ効果を含んだ量子状態が一体演算子formula_8を用いて次のように与えられると仮定する。
そして、次にこのように与えられた状態を用いて計算されるハミルトニアンの期待値をformula_10に関してテイラー展開すると
次のようになる。
formula_12の期待値がゼロになるように求めるのがハートリー-フォック近似であるので右辺第2項はゼロとなる。
従って、
と表されることがわかる。ここでformula_14とformula_15は二重交換関係formula_16
と等価である。
応用: 相互作用するボース粒子系のRPA基底状態.
ボソン系のRPA真空formula_17は、相関のないボソン真空formula_18と元々のボソン励起formula_19で表すことができる。
ここで"Z"はformula_22を満たす対称行列である。
この正規化は次のように計算される。
ここでformula_24
はformula_25の特異値分解である。
元々の励起と新たな励起の結合は、次のように与えられる。