ブロッホ＝ドミニシスの定理

統計物理学において、ブロッホ＝ドミニシスの定理（英: Bloch-De Dominicis theorem）とは、量子多体系における熱平均で定義された多点相関関数を、2点相関関数の組み合わせ和に分解する定理。場の量子論の真空期待値に関するウィックの定理に対し、有限温度の系での類似版に相当しており、ウィックの定理とも呼ばれる。物理学者 松原武生によって、温度グリーン関数の理論展開ともに導入された[1]。定理の名は、最初に完全な証明を与えた物理学者C. ブロッホとC. T. ドミニシスに因む[2]。

概要

A を生成演算子a_α^†、または消滅演算子a_α、もしくはそれらを虚時間で相互作用表示したものとする。ここで、相互作用表示において、非摂動系の自由ハミルトニアンは

${\displaystyle H_{0}=\sum _{\alpha }\varepsilon _{\alpha }a_{\alpha }^{\,\dagger }a_{\alpha }}$

のように2次形式で表されているとする。また平均値${\displaystyle \langle \cdot \rangle _{0}}$は${\displaystyle e^{-\beta (H_{0}-\mu N)}}$による非摂動系でのグランドカノニカル分布の熱平均を表すものとする。

このとき、この熱平均で定義されるn点相関関数は、nが偶数である場合のみゼロにならず、

${\displaystyle \langle A_{1}A_{2}\cdots A_{n}\rangle _{0}=\left\{{\begin{matrix}0&\quad (n:{\mbox{odd}})\\\sum _{m=2}^{n}(\pm 1)^{m}\langle A_{1}A_{m}\rangle _{0}\langle A_{2}\cdots A_{m+1}A_{m-1}\cdots A_{n}\rangle _{0}&\quad (n:{\mbox{even}})\end{matrix}}\right.}$

が成り立つ。ここで現れる2点相関関数は縮約（contraction）と呼ばれる。また、(±1)^mの項はフェルミ粒子での演算子の順番の並べ替えにおいて、隣合う演算子同士を置き換える際に生じる符号の反転を表しており、符号は正がボーズ粒子、負がフェルミ粒子に対応するものとする。以降、本項に現れる複合の符号は全て、上がボーズ粒子、下がフェルミ粒子に対応するものとする。

さらに、nが偶数であるときは、この結果を繰り返し適用することで、

${\displaystyle \langle A_{1}A_{2}\cdots A_{n}\rangle _{0}=\sum _{P}'(\pm 1)^{P}\langle A_{i_{1}}A_{i_{2}}\rangle _{0}\langle A_{i_{3}}A_{i_{4}}\rangle _{0}\cdots \langle A_{i_{n-1}}A_{i_{n}}\rangle _{0}}$

と全ての演算子の縮約の組み合わせ和に分解できる。ここで、P は(1, 2, …, n)→(i₁, i₂, …, i_n)なる置換を表し、和Σ'において、ここの縮約で対となる演算子は、i_k-1 < i_kを満たし、全体としてはi₁ < i₃ < …< i_nの順序が満たされる項について和をとるものする。また、(±1)^Pは、フェルミ粒子の隣合う演算子同士の並べ替えの際に付与する正負の符号の変化を表し、ボーズ粒子については+1、フェルミ粒子については演算子の並べ替えの回数に応じた符号（置換の符号）を与えるものとする。これをブロッホ＝ドミニシスの定理またはウィックの定理と呼ぶ[3][4]。

縮約について、記法

${\displaystyle A^{\bullet }B^{\bullet }:=\langle AB\rangle _{0}}$

を導入すれば、全ての可能な縮約の組み合わせをとるというブロッホ＝ドミニシスの定理は

${\displaystyle \langle A_{1}A_{2}\cdots A_{n}\rangle _{0}=A_{1}^{\bullet }A_{2}^{\bullet }A_{3}^{\bullet \bullet }\cdots A_{n}^{\bullet \bullet \bullet }+A_{1}^{\bullet }A_{2}^{\bullet \bullet }A_{3}^{\bullet }\cdots A_{n}^{\bullet \bullet \bullet }+A_{1}^{\bullet }A_{2}^{\bullet \bullet }A_{3}^{\bullet \bullet }\cdots A_{n}^{\bullet \bullet \bullet }+\cdots }$

とも表すことができる。但し、複数個の演算子の縮約については、

${\displaystyle {\begin{aligned}A^{\bullet }B^{\bullet \bullet }C^{\bullet }D^{\bullet \bullet }&=\pm (A^{\bullet }C^{\bullet })(B^{\bullet \bullet }D^{\bullet \bullet })=\pm \langle AC\rangle _{0}\langle BD\rangle _{0}\\A^{\bullet }B^{\bullet \bullet }C^{\bullet \bullet }D^{\bullet }&=(A^{\bullet }D^{\bullet })(B^{\bullet \bullet }C^{\bullet \bullet })=\langle AD\rangle _{0}\langle BC\rangle _{0}\end{aligned}}}$

のように、同じ右付き添え字の演算子同士の縮約を行い、さらにフェルミ粒子の場合には演算子の入れ替えの回数に応じた符号を与えるものとする。

A =A (τ)が虚時間での相互作用表示の演算子であるとし、虚時間に対する時間順序積をとる場合にも、同様にブロッホ＝ドミニシスの定理

${\displaystyle \langle T_{\tau }[A_{1}(\tau _{1})A_{2}(\tau _{2})\cdots A_{n}(\tau _{n})]\rangle _{0}=\sum _{P}'(\pm 1)^{P}\langle T_{\tau }[A_{i_{1}}(\tau _{i_{1}})A_{i_{2}}(\tau _{i_{2}})]\rangle _{0}\langle T_{\tau }[A_{i_{3}}(\tau _{i_{3}})A_{i_{4}}(\tau _{i_{4}})]\rangle _{0}\cdots \langle T_{\tau }[A_{i_{n-1}}(\tau _{i_{n-1}})A_{i_{n}}(\tau _{i_{n}})]\rangle _{0}}$

が成り立つ。

この場合も縮約について、

${\displaystyle A(\tau _{a})^{\circ }B(\tau _{b})^{\circ }:=\langle T_{\tau }[A(\tau _{a})B(\tau _{b})]\rangle _{0}}$

となる記法を導入すれば、

${\displaystyle {\begin{aligned}\langle T_{\tau }[A_{1}(\tau _{1})A_{2}(\tau _{2})\cdots A_{n}(\tau _{n})]\rangle _{0}=&A_{1}(\tau _{1})^{\circ }A_{2}(\tau _{2})^{\circ }A_{3}(\tau _{3})^{\circ \circ }\cdots A_{n}(\tau _{n})^{\circ \circ \circ }\\&+A_{1}(\tau _{1})^{\circ }A_{2}(\tau _{2})^{\circ \circ }A_{3}(\tau _{3})^{\circ }\cdots A_{n}(\tau _{n})^{\circ \circ \circ }+A_{1}(\tau _{1})^{\circ }A_{2}(\tau _{2})^{\circ \circ }A_{3}(\tau _{3})^{\circ \circ }\cdots A_{n}(\tau _{n})^{\circ \circ \circ }+\cdots \end{aligned}}}$

と表すことができる。

具体例

ブロッホ＝ドミニシスの定理により、3点相関関数、4点相関関数について

${\displaystyle \langle ABC\rangle _{0}=0}$

${\displaystyle {\begin{aligned}\langle ABCD\rangle _{0}&=A^{\bullet }B^{\bullet }C^{\bullet \bullet }D^{\bullet \bullet }+A^{\bullet }B^{\bullet \bullet }C^{\bullet }D^{\bullet \bullet }+A^{\bullet }B^{\bullet \bullet }C^{\bullet \bullet }D^{\bullet }\\&=\langle AB\rangle _{0}\langle CD\rangle _{0}\pm \langle AC\rangle _{0}\langle BD\rangle _{0}+\langle AD\rangle _{0}\langle BC\rangle _{0}\end{aligned}}}$

が成り立つ。

時間順序積をとる場合にも同様に

${\displaystyle \langle T_{\tau }[A(\tau _{a})B(\tau _{b})C(\tau _{c})]\rangle _{0}=0}$

${\displaystyle {\begin{aligned}\langle T_{\tau }[A(\tau _{a})B(\tau _{b})C(\tau _{c})D(\tau _{d})]\rangle _{0}=&A(\tau _{a})^{\circ }B(\tau _{b})^{\circ }C(\tau _{c})^{\circ \circ }D(\tau _{d})^{\circ \circ }\\&+A(\tau _{a})^{\circ }B(\tau _{b})^{\circ \circ }C(\tau _{c})^{\circ }D(\tau _{d})^{\circ \circ }+A(\tau _{a})^{\circ }B(\tau _{b})^{\circ \circ }C(\tau _{c})^{\circ \circ }D(\tau _{d})^{\circ }\\=&\langle T_{\tau }[A(\tau _{a})B(\tau _{b})]\rangle _{0}\langle T_{\tau }[C(\tau _{c})D(\tau _{d})]\rangle _{0}\\&\pm \langle T_{\tau }[A(\tau _{a})C(\tau _{c})]\rangle _{0}\langle T_{\tau }[B(\tau _{b})D(\tau _{d})]\rangle _{0}+\langle T_{\tau }[A(\tau _{a})D(\tau _{d})]\rangle _{0}\langle T_{\tau }[B(\tau _{b})C(\tau _{c})]\rangle _{0}\end{aligned}}}$

が成り立つ。

理論の背景

量子多体系における温度グリーン関数の理論では、温度グリーン関数によって系の様々な物理量を求めることができるともに、摂動計算を系統的に行うことができる。ここで演算子A, B の温度グリーン関数は

${\displaystyle G_{AB}(\tau ,\tau '):=-\langle T_{\tau }[A(\tau )B(\tau ')]\rangle }$

で定義される2点相関関数である。但し、記号${\displaystyle \langle \cdots \rangle }$は

${\displaystyle \langle \cdots \rangle ={\frac {\mathrm {Tr} \{e^{-\beta (H-\mu N)}\cdots \}}{\mathrm {Tr} \{e^{-\beta (H-\mu N)}\}}}={\frac {\mathrm {Tr} \{e^{-\beta K}\cdots \}}{\mathrm {Tr} \{e^{-\beta K}\}}}\quad (K:=H-\mu N)}$

で定義されるグランドカノニカル分布での熱平均であり、H はハミルトニアン、N は数演算子、βは逆温度、μは化学ポテンシャルを表す。またA (τ)は

${\displaystyle {\begin{aligned}A(\tau )&=e^{+K\tau }Ae^{-K\tau }\\A(\tau )^{\dagger }&=e^{+K\tau }A^{\dagger }e^{-K\tau }\end{aligned}}}$

で定義される虚時間τ=it についてのハイゼンベルク表示の演算子である。T_τは虚時間についての時間順序積であり、

${\displaystyle T_{\tau }\{A_{1}(\tau _{1})A_{2}(\tau _{2})\}=\left\{{\begin{matrix}A_{1}(\tau _{1})A_{2}(\tau _{2})&\quad (\tau _{1}>\tau _{2})\\\pm A_{2}(\tau _{2})A_{1}(\tau _{1})&\quad (\tau _{2}>\tau _{1})\end{matrix}}\right.}$

を意味する。

一般に温度グリーン関数の計算において、相互作用にある系では、

${\displaystyle H=H_{0}+V\,}$

とハミルトニアンを可解な非摂動項と相互作用を含む摂動項に分け、相互作用表示の演算子

${\displaystyle {\begin{aligned}A_{I}(\tau )&=e^{+K_{0}\tau }Ae^{-K_{0}\tau }\\A_{I}(\tau )^{\dagger }&=e^{+K_{0}\tau }A^{\dagger }e^{-K_{0}\tau }\quad (K_{0}:=H_{0}-\mu N)\end{aligned}}}$

に対して、摂動計算を行うことが必要となる。このとき、摂動計算において、

${\displaystyle \langle A_{1I}(\tau _{1})\cdots A_{nI}(\tau _{n})\rangle _{0}={\frac {\mathrm {Tr} \{e^{-\beta K_{0}}A_{1I}(\tau _{1})\cdots A_{nI}(\tau _{n})\}}{\mathrm {Tr} \{e^{-\beta K_{0}}\}}}}$

という高次の相関関数が現れる。ブロッホ＝ドミニシスの定理は、こうした多点相関関数を縮約

${\displaystyle {\begin{aligned}\langle a_{\alpha }^{\,\dagger }a_{\alpha '}\rangle _{0}&={\frac {\delta _{\alpha \alpha '}}{1\mp e^{\beta \varepsilon _{\alpha }}}}\\\langle a_{\alpha }a_{\alpha '}^{\,\dagger }\rangle _{0}&={\frac {\delta _{\alpha \alpha '}}{1\mp e^{-\beta \varepsilon _{\alpha }}}}\\\langle a_{\alpha }^{\,\dagger }a_{\alpha '}^{\,\dagger }\rangle _{0}&=0\\\langle a_{\alpha }a_{\alpha '}\rangle _{0}&=0\end{aligned}}:}$

に分解し、実際の計算を可能にする。

ガウス過程との関係

ブロッホ＝ドミニシスの定理は、古典系におけるガウス過程の持つ性質を量子系に拡張したものに相当する [5]。実際、分布が

${\displaystyle P_{n}(x_{1},\cdots ,x_{n})={\frac {1}{\sqrt {\det {A}}}}\exp {{\bigl (}-\sum _{i,j}A_{ij}(x_{i}-\langle X_{i}\rangle )(x_{j}-\langle X_{j}\rangle ){\bigr )}}}$

で与えられるガウス過程[6]を考えると

${\displaystyle \langle X_{1}X_{2}\cdots X_{n}\rangle =\sum _{\operatorname {all\,\,partition} }\prod \langle X_{i_{1}}\cdots X_{i_{m}}\rangle _{c}}$

が成り立つ。ここで<…>はこの分布に対する期待値、<…>_cはキュムラントを表すものとする。 また、右辺の和は、X₁,…,X_nをいくつかの集まりに分割する全ての組み合わせにわたってとるものである。例えば、3点相関関数、4点相関関数については、

${\displaystyle {\begin{aligned}\langle X_{1}X_{2}\rangle =&\langle X_{1}X_{2}\rangle _{c}+\langle X_{1}\rangle _{c}\langle X_{2}\rangle _{c}\\\langle X_{1}X_{2}X_{3}\rangle =&\langle X_{1}X_{2}X_{3}\rangle _{c}+\langle X_{1}\rangle _{c}\langle X_{2}X_{3}\rangle _{c}+\langle X_{2}\rangle _{c}\langle X_{1}X_{3}\rangle _{c}+\langle X_{3}\rangle _{c}\langle X_{1}X_{2}\rangle _{c}+\langle X_{1}\rangle _{c}\langle X_{2}\rangle _{c}\langle X_{3}\rangle _{c}\end{aligned}}}$

である。

さらに全てのX_iについて、⟨X_i⟩=0であるとするならば、ブロッホ＝ドミニシスの定理と同様にn点相関関数は、nが偶数である場合のみゼロにならず、

${\displaystyle {\begin{aligned}\langle X_{1}X_{2}\cdots X_{n}\rangle &=\sum _{\operatorname {all\,\,pair} }\prod \langle X_{i}X_{j}\rangle \\&=\sum _{P}'\langle X_{i_{1}}X_{i_{2}}\rangle \langle X_{i_{3}}X_{i_{4}}\rangle \cdots \langle X_{i_{n-1}}X_{i_{n}}\rangle \end{aligned}}}$

と2点相関関数の組み合わせ和に分解される。

脚注

1. T. Matsubara, Prog. Theor. Phys., 14, p.351 (1955)
2. C. Bloch and C. T. de Dominicis, Nucl. Phys., 7, p.459 (1958)
3. A. L. Fetter and J. D. Walecka (2003)
4. 阿部龍蔵 (1992)
5. 今田正俊 (2004)
6. A =(A_ij)は正定値行列

参考文献

論文

• G. C. Wick (1950). “The Evaluation of the Collision Matrix”. Phys. Rev. 80: 268. doi:10.1103/PhysRev.80.268. 
• T. Matsubara (1955). “A New Approach to Quantum-Statistical Mechanics”. Prog. Theor. Phys. 14: 351. doi:10.1143/PTP.14.351. 
• C. Bloch; C. T. de Dominicis (1958). “Un développement du potentiel de gibbs d'un système quantique composé d'un grand nombre de particules”. Nucl. Phys. 7: 459. doi:10.1016/0029-5582(58)90285-2. 
• M. Gaudin (1960). “Une démonstration simplifiée du théorème de wick en mécanique statistique”. Nucl. Phys. 15: 89. doi:10.1016/0029-5582(60)90285-6. 

書籍

• Alexander L. Fetter, John Dirk Walecka (2003). Quantum Theory of Many-Particle Systems. Dover Publications. ISBN 978-0486428277 
• 阿部龍蔵『統計力学』（2版）東京大学出版会、1992年。ISBN 978-4130621342。 
• 西川恭治、森弘之『統計物理学』朝倉書店〈朝倉物理学大系〉、2000年。ISBN 978-4254136807。 
• 今田正俊『統計物理学』丸善、2004年。ISBN 978-4621074831。 
• 高田康民『多体問題』朝倉書店〈朝倉物理学大系〉、1999年。ISBN 978-4-254-13679-1。 

関連項目

• 場の量子論 - ウィックの定理
• 量子多体系 - 温度グリーン関数