乱雑位相近似のソースを表示

[[File:Random phase approximation ring diagrams.png|thumb|400px|ファインマンダイアグラムにおいて、乱雑位相近似(RPA)はリングダイアグラムの和として表される。上の太い線は相互作用するGreen関数、細い線は相互作用のないGreen関数、破線は2体の相互作用を表す。]]

'''乱雑位相近似'''（らんざついそうきんじ、{{lang-en|Random Phase Approximation}}, '''RPA'''）とは、元々[[デヴィッド・ボーム]]と{{仮リンク|デヴィッド・パインズ|en|David Pines}}によって展開された[[多体問題 (量子論)|多体系]]における基底状態の[[量子揺らぎ]]及び励起振動状態([[フォノン]])を記述するための近似手法。[[線形応答理論]]における[[摂動論]]的な近似法の一つである。

粒子系（[[電子ガス]]など）が高密度の場合は、乱雑位相近似が妥当な近似であることが分かっている。

同等な近似手法が、多方面（例：[[GW近似]]）で利用、応用されている。

== 第一量子化でのRPA ==
N粒子系（N電子系）における[[密度演算子]]は次のように与えられる。

:<math>\rho(\vec{q})=\sum_{i=1}^{N} \exp(i\vec{q}\cdot\vec{r}_i)</math>

ここで位置座標ベクトル<math>\vec{r}_i</math>が無秩序であれば、逆格子ベクトルと位置座標ベクトルとの積、<math>\vec{q}\cdot\vec{r}_i</math>も無秩序（[[乱雑]]）なので、<math>\rho(\vec{q} \neq 0)</math>からの寄与が<math>\rho(\vec{q} = 0)</math>よりずっと小さいとして無視できる。これを乱雑位相近似(RPA)という。

<math>\vec{q} \neq 0</math>においては、<math>\vec{q}\cdot\vec{r}_i</math>が乱雑なことにより各項の位相も乱雑となり、和の各成分が相殺し合って全体としての寄与が無視できるほど小さくなることによる。勿論、この近似が適用できない場合も多々ある。

== 概要 ==
RPAは1952年と1953年にボームとパインズによって初めて導入された<ref>D. Bohm and D. Pines: ''A Collective Description of Electron Interactions. I. Magnetic Interactions'', Phys. Rev. '''82''', 625–634 (1951) ([http://prola.aps.org/abstract/PR/v82/i5/p625_1 abstract])</ref><ref>D. Pines and D. Bohm: ''A Collective Description of Electron Interactions: II. Collective vs Individual Particle Aspects of the Interactions'', Phys. Rev. '''85''', 338–353 (1952) ([http://prola.aps.org/abstract/PR/v85/i2/p338_1 abstract])</ref><ref>D. Bohm and D. Pines: ''A Collective Description of Electron Interactions: III. Coulomb Interactions in a Degenerate Electron Gas'', Phys. Rev. '''92''', 609–625 (1953) ([http://prola.aps.org/abstract/PR/v92/i3/p609_1 abstract])</ref>。それまで何十年もの間、電子間のミクロな量子力学的相互作用の効果を物質の理論に取り入れようとする試みがあった。ボームとパインズのRPAは、弱く遮蔽された[[クーロン相互作用]]を説明し、電子系における電子の動的な線形応答を記述するために用いられる。

RPAでは、電子は全電位''V''('''r''')（つまり外部[[摂動]]ポテンシャル ''V''<sub>ext</sub>('''''r''''') と遮蔽ポテンシャル ''V''<sub>sc</sub>('''''r''''') の和）にのみ応答すると仮定される。外部摂動ポテンシャルは単一の周波数 ''ω'' で振動すると仮定されるので、このモデルに[[自己無撞着場]]（SCF）法<ref>H. Ehrenreich and M. H. Cohen, [https://doi.org/10.1103/PhysRev.115.786 Phys. Rev. '''115''', 786 (1959)]</ref>を適用すると動的誘電関数 ''ε''<sub>RPA</sub>('''''k''''', ''ω'')が得られる。

誘電関数への全電位の寄与は平均化される仮定するため、波数ベクトル '''''k''''' における電位のみが寄与する。これが乱雑位相近似が意味するものである。結果として生じる誘電関数は[[リンドハード理論|リンドハード誘電関数]]とも呼ばれ<ref>J. Lindhard, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. '''28''', 8 (1954)</ref><ref>N. W. Ashcroft and N. D. Mermin, ''Solid State Physics'' (Thomson Learning, Toronto, 1976)</ref>、電子ガスの多くの性質（[[プラズモン]]など）を正確に予測している<ref>G. D. Mahan, ''Many-Particle Physics'', 2nd ed. (Plenum Press, New York, 1990)</ref>。

RPAは自由度を過大に評価していると50年代後半に批判され、その正当化には理論物理学者の多くの労力が費やされた。[[マレー・ゲルマン]]と[[キース・ブルックナー]]は、高密度の電子ガスのファインマンダイアグラムにおける最低次のチェーンの和からRPAが導かれることを示した<ref>M. Gell-Mann, K.A. Brueckner,  ''Correlation energy of an electron gas at high density'', Phys. Rev. '''106''', 364 (1957)</ref>。

これらの結果の一貫性はRPAの正当化には重要であり、50年代後半と60年代の理論物理学は大きく発展した。

== RPA方程式 ==
まず第0近似として[[ハートリー-フォック]]近似を考える。ハートリー-フォック近似で得られた基底状態には量子揺らぎ効果は含まれてはいない。
そこで、量子揺らぎ効果を含んだ量子状態が一体演算子<math>\hat{F}</math>を用いて次のように与えられると仮定する。

:<math>|\Psi\rangle=e^{i\lambda\hat{F}}|\Phi _{\mathrm{HF}}\rangle</math>

そして、次にこのように与えられた状態を用いて計算されるハミルトニアンの期待値を<math>\lambda</math>に関してテイラー展開すると
次のようになる。

:<math>
\langle\Psi|\hat{H}|\Psi\rangle
=
\langle\Phi _{\mathrm{HF}}|\hat{H}-i\lambda[\hat{F},\hat{H}]+\frac{\lambda^2}{2}[\hat{F},[\hat{H},\hat{F}]]+\dotsb|\Phi _{\mathrm{HF}}\rangle
</math>

<math>[\hat{F},H]</math>の期待値がゼロになるように求めるのがハートリー-フォック近似であるので右辺第2項はゼロとなる。
従って、

:<math>
\begin{align}
\langle\Psi|\hat{H}|\Psi\rangle
&=
\langle\Phi _{\mathrm{HF}}|\hat{H}|\Phi _{\mathrm{HF}}\rangle+\frac{\lambda^2}{2}\langle\Phi _{\mathrm{HF}}|[\hat{F},[\hat{H},\hat{F}]]|\Phi _{\mathrm{HF}}\rangle+\dotsb \\
&=
E _{\mathrm{HF}}
+
\frac{\lambda^2}{2}
\sum_{minj}
\begin{pmatrix}f^*_{mi} & -f_{im}\end{pmatrix}
\begin{pmatrix}
A & B \\ B^* & A^*
\end{pmatrix}_{minj}
\begin{pmatrix}f_{nj} \\ -f^*_{jn}\end{pmatrix}
+\dotsb
\end{align}
</math>

と表されることがわかる。ここで<math>A</math>と<math>B</math>は二重交換関係<math>[X,Y,Z]=\frac{1}{2}[X,[Y,Z]]+\frac{1}{2}[[X,Y],Z]</math>を用いて
:<math>A_{minj}=\langle\Phi _{\mathrm{HF}}|[a_i^\dagger a_m,\hat{H},a_n^\dagger a_j]|\Phi _{\mathrm{HF}}\rangle</math>
:<math>B_{minj}=-\langle\Phi _{\mathrm{HF}}|[a_i^\dagger a_m,\hat{H},a_j^\dagger a_n]|\Phi _{\mathrm{HF}}\rangle</math>
と定義されている。乱雑位相近似は、'''これまでの計算で現れた行列'''
<math>
\begin{pmatrix}
A & B \\ B^* & A^*
\end{pmatrix}
</math>
'''を対角化するための固有値方程式を考え、その固有値と固有ベクトルを求めること'''、という言い方ができる。固有値及び固有ベクトルを求める方程式は'''RPA方程式'''と呼ばれ、次のような形で与えられる。
:<math>
\begin{pmatrix}
A & B \\ B^* & A^*
\end{pmatrix}
\begin{pmatrix}
X^\nu \\ Y^\nu
\end{pmatrix}
=
\hbar\omega_\nu
\begin{pmatrix}
1 & 0 \\ 0 & -1
\end{pmatrix}
\begin{pmatrix}
X^\nu \\ Y^\nu
\end{pmatrix}
</math>

ここで
<math>
\begin{pmatrix}
X^\nu \\ Y^\nu
\end{pmatrix}
</math>
は固有ベクトルであり、<math>\hbar\omega_\nu</math>は固有値であり励起状態を表す。

また、RPA方程式から得られる固有値が正の値をとる時、ハートレーフォック基底状態はエネルギーの極小値であることから
系のエネルギーは安定であることがわかる。しかし、固有値のなかに一つでも負の値のものが含まれる場合、
もはや安定ではなく異なる基底状態(真空)が存在する可能性、つまり[[相転移]]の可能性を示唆している。

固有ベクトルと固有値の存在は量子状態<math>|\nu\rangle</math>が、状態
<math>|mi\rangle=a^\dagger_m a_i|HF\rangle</math>の線形結合を用いて、

:<math>
|\nu\rangle (=O^\dagger_\nu |\Phi_{RPA}\rangle)
= \sum_{mi}(X_{mi}^\nu |mi\rangle-Y_{mi}^\nu |im\rangle)
</math>

と表せることを示している。この時、量子状態<math>|\nu\rangle</math>はその異なるもの同士は直交する、すなわち<math>\langle\nu|\nu'\rangle=\delta_{\nu\nu'}</math>と仮定する。

更に<math>|mi\rangle=a^\dagger_m a_i|HF\rangle</math>の線形結合で定義される状態<math>|\nu\rangle</math>
の最もエネルギーの低い状態(基底状態)<math>|\Phi_{RPA}\rangle</math>を
<math>O_\nu |\Phi_{RPA}\rangle = 0</math>
と定義する。

以上の条件のもとで上述のRPA固有値方程式は

:<math>
\langle \Phi_{RPA} |[O_{\nu'},[H,O^\dagger_\nu]]|\Phi_{RPA} \rangle= \hbar\omega_\nu\delta_{\nu\nu'}
</math>

と等価である。

==応用: 相互作用するボース粒子系のRPA基底状態 ==
ボソン系のRPA真空<math>\left|\mathbf{RPA}\right\rangle</math>は、相関のないボソン真空<math>\left|\mathbf{MFT}\right\rangle</math>と元々のボソン励起<math>\mathbf{a}_{i}^{\dagger}</math>で表すことができる。
:<math>\left|\mathrm{RPA}\right\rangle=\mathcal{N}\mathbf{e}^{Z_{ij}\mathbf{a}_{i}^{\dagger}\mathbf{a}_{j}^{\dagger}/2}\left|\mathrm{MFT}\right\rangle</math>
:<math>\mathcal{N}= \frac{\left\langle  \mathrm{MFT}\right|\left.\mathrm{RPA}\right\rangle}{\left\langle  \mathrm{MFT}\right|\left.\mathrm{MFT}\right\rangle}</math>

ここで''Z''は<math>|Z|\leq 1</math>を満たす対称行列である。

この正規化は次のように計算される。
:<math>\langle 
\mathrm{RPA}|\mathrm{RPA}\rangle=
\mathcal{N}^2 \langle \mathrm{MFT}|
\mathbf{e}^{z_{i}(\tilde{\mathbf{q}}_{i})^2/2}
\mathbf{e}^{z_{j}(\tilde{\mathbf{q}}^{\dagger}_{j})^2/2}
| \mathrm{MFT}\rangle=1
</math>
ここで<math>Z_{ij}=(X^{\mathrm{t}})_{i}^{k} z_{k} X^{k}_{j}</math>
は<math>Z_{ij}</math>の[[特異値分解]]である。
:<math>\tilde{\mathbf{q}}^{i}=(X^{\dagger})^{i}_{j}\mathbf{a}^{j}</math>
:<math>\begin{align}
\mathcal{N}^{-2}&=
\sum_{m_{i}}\sum_{n_{j}}  \frac{(z_{i}/2)^{m_{i}}(z_{j}/2)^{n_{j}}}{m!n!}
\langle \mathrm{MFT}|
\prod_{i\,j}
(\tilde{\mathbf{q}}_{i})^{2 m_{i}}
(\tilde{\mathbf{q}}^{\dagger}_{j})^{2 n_{j}}
| \mathrm{MFT}\rangle \\
&=\prod_{i}
\sum_{m_{i}} (z_{i}/2)^{2 m_{i}} \frac{(2 m_{i})!}{m_{i}!^2} \\
&=\prod_{i}\sum_{m_{i}} (z_{i})^{2 m_{i}} {1/2 \choose m_{i}}=\sqrt{\det(1-|Z|^2)}
\end{align}
</math>

元々の励起と新たな励起の結合は、次のように与えられる。
:<math>\tilde{\mathbf{a}}_{i}=\left(\frac{1}{\sqrt{1-Z^2}}\right)_{ij}\mathbf{a}_{j}+
\left(\frac{1}{\sqrt{1-Z^2}}Z\right)_{ij}\mathbf{a}^{\dagger}_{j}</math>.

==関連記事==
*[[物理学]]
*[[物性物理学]]
*[[原子核物理学]]

== 参考文献 ==
<references />

{{DEFAULTSORT:らんさついそうきんじ}}
[[Category:計算物理学]]
[[Category:量子化学]]