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

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

概要

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

$H_0=\sum _{\alpha }{\epsilon }_{\alpha }{a}_{\alpha }^{†}{a}_{\alpha }$

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

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

$⟨A_1{A}_2\cdots A_n{⟩}_0=\{\begin{array}{cc}0& (n:\text{odd})\\ \sum _{m=2}^n(\pm 1{)}^m⟨A_1{A}_m{⟩}_0⟨A_2\cdots A_{m+1}{A}_{m-1}\cdots A_n{⟩}_0& (n:\text{even})\end{array}$

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

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

$⟨A_1{A}_2\cdots A_n{⟩}_0=\sum _{P^{\prime }}(\pm 1{)}^P⟨A_{i_1}{A}_{i_2}{⟩}_0⟨A_{i_3}{A}_{i_4}{⟩}_0\cdots ⟨A_{i_{n-1}}{A}_{i_n}{⟩}_0$

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

縮約について、記法

$A^{\bullet }{B}^{\bullet }:=⟨AB{⟩}_0$

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

$⟨A_1{A}_2\cdots A_n{⟩}_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$

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

$\begin{array}{cc}{A}^{\bullet }{B}^{\bullet \bullet }{C}^{\bullet }{D}^{\bullet \bullet }& =\pm (A^{\bullet }{C}^{\bullet })(B^{\bullet \bullet }{D}^{\bullet \bullet })=\pm ⟨AC{⟩}_0⟨BD{⟩}_0\\ A^{\bullet }{B}^{\bullet \bullet }{C}^{\bullet \bullet }{D}^{\bullet }& =(A^{\bullet }{D}^{\bullet })(B^{\bullet \bullet }{C}^{\bullet \bullet })=⟨AD{⟩}_0⟨BC{⟩}_0\end{array}$

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

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

$⟨T_{\tau }[A_1({\tau }_1)A_2({\tau }_2)\cdots A_n({\tau }_n)]{⟩}_0=\sum _{P^{\prime }}(\pm 1{)}^P⟨T_{\tau }[A_{i_1}({\tau }_{i_1})A_{i_2}({\tau }_{i_2})]{⟩}_0⟨T_{\tau }[A_{i_3}({\tau }_{i_3})A_{i_4}({\tau }_{i_4})]{⟩}_0\cdots ⟨T_{\tau }[A_{i_{n-1}}({\tau }_{i_{n-1}})A_{i_n}({\tau }_{i_n})]{⟩}_0$

が成り立つ。

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

$A({\tau }_a{)}^{\circ }B({\tau }_b{)}^{\circ }:=⟨T_{\tau }[A({\tau }_a)B({\tau }_b)]{⟩}_0$

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

$\begin{array}{cc}⟨T_{\tau }[A_1({\tau }_1)A_2({\tau }_2)\cdots A_n({\tau }_n)]{⟩}_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{array}$

と表すことができる。

具体例

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

$⟨ABC{⟩}_0=0$

$\begin{array}{cc}⟨ABCD{⟩}_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 }\\ & =⟨AB{⟩}_0⟨CD{⟩}_0\pm ⟨AC{⟩}_0⟨BD{⟩}_0+⟨AD{⟩}_0⟨BC{⟩}_0\end{array}$

が成り立つ。

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

$⟨T_{\tau }[A({\tau }_a)B({\tau }_b)C({\tau }_c)]{⟩}_0=0$

$\begin{array}{cc}⟨T_{\tau }[A({\tau }_a)B({\tau }_b)C({\tau }_c)D({\tau }_d)]{⟩}_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 }\\ =& ⟨T_{\tau }[A({\tau }_a)B({\tau }_b)]{⟩}_0⟨T_{\tau }[C({\tau }_c)D({\tau }_d)]{⟩}_0\\ & \pm ⟨T_{\tau }[A({\tau }_a)C({\tau }_c)]{⟩}_0⟨T_{\tau }[B({\tau }_b)D({\tau }_d)]{⟩}_0+⟨T_{\tau }[A({\tau }_a)D({\tau }_d)]{⟩}_0⟨T_{\tau }[B({\tau }_b)C({\tau }_c)]{⟩}_0\end{array}$

が成り立つ。

理論の背景

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

$G_{AB}(\tau ,{\tau }^{\prime }):=-⟨T_{\tau }[A(\tau )B({\tau }^{\prime })]⟩$

で定義される2点相関関数である。但し、記号$⟨\cdots ⟩$は

$⟨\cdots ⟩=\frac{\mathrm{T}\mathrm{r}\{e^{-\beta (H-\mu N)}\cdots \}}{\mathrm{T}\mathrm{r}\{e^{-\beta (H-\mu N)}\}}=\frac{\mathrm{T}\mathrm{r}\{e^{-\beta K}\cdots \}}{\mathrm{T}\mathrm{r}\{e^{-\beta K}\}}(K:=H-\mu N)$

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

$\begin{array}{cc}A(\tau )& =e^{+K\tau }A{e}^{-K\tau }\\ A(\tau {)}^{†}& =e^{+K\tau }{A}^{†}{e}^{-K\tau }\end{array}$

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

$T_{\tau }\{A_1({\tau }_1)A_2({\tau }_2)\}=\{\begin{array}{cc}{A}_1({\tau }_1)A_2({\tau }_2)& ({\tau }_1>{\tau }_2)\\ \pm A_2({\tau }_2)A_1({\tau }_1)& ({\tau }_2>{\tau }_1)\end{array}$

を意味する。

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

$H=H_0+V$

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

$\begin{array}{cc}{A}_I(\tau )& =e^{+K_0\tau }A{e}^{-K_0\tau }\\ A_I(\tau {)}^{†}& =e^{+K_0\tau }{A}^{†}{e}^{-K_0\tau }(K_0:=H_0-\mu N)\end{array}$

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

$⟨A_{1I}({\tau }_1)\cdots A_{nI}({\tau }_n){⟩}_0=\frac{\mathrm{T}\mathrm{r}\{e^{-\beta K_0}{A}_{1I}({\tau }_1)\cdots A_{nI}({\tau }_n)\}}{\mathrm{T}\mathrm{r}\{e^{-\beta K_0}\}}$

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

$\begin{array}{cc}⟨a_{\alpha }^{†}{a}_{{\alpha }^{\prime }}{⟩}_0& =\frac{{\delta }_{\alpha {\alpha }^{\prime }}}{1\mp e^{\beta {\epsilon }_{\alpha }}}\\ ⟨a_{\alpha }{a}_{{\alpha }^{\prime }}^{†}{⟩}_0& =\frac{{\delta }_{\alpha {\alpha }^{\prime }}}{1\mp e^{-\beta {\epsilon }_{\alpha }}}\\ ⟨a_{\alpha }^{†}{a}_{{\alpha }^{\prime }}^{†}{⟩}_0& =0\\ ⟨a_{\alpha }{a}_{{\alpha }^{\prime }}{⟩}_0& =0\end{array}:$

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

ガウス過程との関係

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

$P_n(x_1,\cdots ,x_n)=\frac{1}{\sqrt{detA}}\exp (-\sum _{i,j}{A}_{ij}(x_i-⟨X_i⟩)(x_j-⟨X_j⟩))$

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

$⟨X_1{X}_2\cdots X_n⟩=\sum _{\mathrm{a}\mathrm{l}\mathrm{l}\mathrm{p}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{i}\mathrm{t}\mathrm{i}\mathrm{o}\mathrm{n}}\prod ⟨X_{i_1}\cdots X_{i_m}{⟩}_c$

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

$\begin{array}{cc}⟨X_1{X}_2⟩=& ⟨X_1{X}_2{⟩}_c+⟨X_1{⟩}_c⟨X_2{⟩}_c\\ ⟨X_1{X}_2{X}_3⟩=& ⟨X_1{X}_2{X}_3{⟩}_c+⟨X_1{⟩}_c⟨X_2{X}_3{⟩}_c+⟨X_2{⟩}_c⟨X_1{X}_3{⟩}_c+⟨X_3{⟩}_c⟨X_1{X}_2{⟩}_c+⟨X_1{⟩}_c⟨X_2{⟩}_c⟨X_3{⟩}_c\end{array}$

である。

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

$\begin{array}{cc}⟨X_1{X}_2\cdots X_n⟩& =\sum _{\mathrm{a}\mathrm{l}\mathrm{l}\mathrm{p}\mathrm{a}\mathrm{i}\mathrm{r}}\prod ⟨X_i{X}_j⟩\\ & =\sum _{P^{\prime }}⟨X_{i_1}{X}_{i_2}⟩⟨X_{i_3}{X}_{i_4}⟩\cdots ⟨X_{i_{n-1}}{X}_{i_n}⟩\end{array}$

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

参考文献

論文

• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal

書籍

• テンプレート:Cite book
• テンプレート:Cite book
• テンプレート:Cite book
• テンプレート:Cite book
• テンプレート:Cite book

関連項目

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