この記事は現在、プロジェクトメンバーによる査読中のものです。草稿段階ですので、内容・表現の正確さについて責任を負いかねます。 リンクを正しく張れていないところが存在しますのでご注意ください。 正式公開まで、いましばらくお待ちください。
グリーン関数を理解しよう(ウィックの定理)

これからいくつかの記事を通して, 物性物理で扱われる絶対零度におけるグリーン関数の理解を目指します. いくつかの定理などの証明は省略して,要点の俯瞰をする方針で行きます. 参考文献として,下に書くMahan先生の本を挙げて おきます.このシリーズでは \hbar=1 とします. 前の記事は 電子とフォノンのグリーン関数 です. 次の記事は ファインマンダイアグラム です.( 目次 )

ちょっと休憩

相互作用 \hat{V}(t) にはどんなものがあるのか, ここで気晴らしもかねて,三つほど挙げてみようと思います. 図において時間は左から右に流れます.

1.電子-電子相互作用

\hat{V}(t) = \dfrac{1}{2}\sum_{\bm{k}\bm{k^\prime}\bm{q}}\sum_{ss^\prime} \dfrac{4 \pi e^2}{q^2} C^\dagger_{\bm{k}+\bm{q},s} C^\dagger_{\bm{k^\prime}-\bm{q},s^\prime} C_{\bm{k^\prime},s^\prime} C_{\bm{k},s} \exp(it(\xi_{\bm{k}+\bm{q}}+\xi_{\bm{k^\prime}-\bm{q}}-\xi_{\bm{k^\prime}}-\xi_{\bm{k}}))\tag{1}

注意として,時間依存性は全て指数関数にまとめられ, C には時間依存性がありません.

chromel-studyGreen04-01.png

2.電子-フォノン相互作用

\hat{V}(t) = \sum_{\bm{q}\bm{k}s} M_{\bm{q}} A_{\bm{q}} C^\dagger_{\bm{k}+\bm{q},s} C_{\bm{k},s}\tag{2}

ここで, M_{\bm{q}}\bm{p} に依存する係数です.

chromel-studyGreen04-02.png

3.電子-フォトン相互作用

これについては,全ハミルトニアンを書きます.

H = \sum_{i} \dfrac{1}{2m} \left[ \bm{p}_i - \dfrac{e}{c}\bm{A}_i \right]^2 + \sum_{i \neq j} \dfrac{e_i e_j}{2r_{ij}} + \sum_{\bm{k} \lambda} \omega_{\bm{k} \lambda} a^\dagger_{\bm{k} \lambda} a_{\bm{k} \lambda} \tag{3}

ここで, \bm{A}_i はベクトルポテンシャルで, a_{\bm{k} \lambda} は光子の消滅演算子でありボゾン演算子です. 実はクーロン反発である1は,この3の相互作用として含まれているそうです.

電子のようなフェルミオンは相互作用で粒子数が保存(陽電子はマイナス一個と数える)のに対し, フォノンの様なボゾンはフェルミオンへ吸収・放出されます.よって,式 (2) の様な 相互作用になるわけです.フォトンもボゾンであるのでフォノンと同じはずですが, 私は式 (3) をどう見ていいか分かりません.

グリーン関数の展開

グリーン関数を徐々に変形していきます. それには 相関関数の計算 の式 (15) で導出した式を使います. 再掲しておきます.

G(\bm{p},t-t^\prime) &= -i \dfrac{_0 \langle | T \hat{C}_{\bm{p}}(t) \hat{C}^\dagger_{\bm{p}}(t^\prime) S(\infty,-\infty) | \rangle_0}{_0 \langle |T  S(\infty,-\infty) | \rangle_0} \tag{4}

これともう一つ, S 行列の展開式( 相互作用表示とS行列 の式 (22) )が必要です.

S(t,t^\prime) = T \exp \left[ -i \int_{t^\prime}^t dt_1 \hat{V}(t_1) \right]\tag{5}

これを

S(\infty, -\infty) = T \exp \left[ -i \int_{-\infty}^\infty dt_1 \hat{V}(t_1) \right]\tag{6}

とします.式 (6) を式 (4) に代入すると,

G(\bm{p},t-t^\prime) &= -i \dfrac{_0 \langle | T \hat{C}_{\bm{p}}(t) \hat{C}^\dagger_{\bm{p}}(t^\prime) \exp \left[ -i \int_{-\infty}^\infty dt_1 \hat{V}(t_1) \right] | \rangle_0}{_0 \langle | S(\infty,-\infty) | \rangle_0} \\&= -i \dfrac{_0 \langle | T \hat{C}_{\bm{p}}(t) \hat{C}^\dagger_{\bm{p}}(t^\prime) \sum_{n=0}^\infty \dfrac{(-i)^{n}}{n!} \int_{-\infty}^\infty dt_1 \cdots \int_{-\infty}^\infty dt_n \hat{V}(t_1) \hat{V}(t_2) \cdots \hat{V}(t_n) | \rangle_0}{_0 \langle | S(\infty,-\infty) | \rangle_0} \\&= \sum_{n=0}^\infty \dfrac{(-i)^{n+1}}{n!} \int_{-\infty}^\infty dt_1 \cdots \int_{-\infty}^\infty dt_n \dfrac{_0 \langle | T \hat{C}_{\bm{p}}(t) \hat{C}^\dagger_{\bm{p}}(t^\prime) \hat{V}(t_1) \hat{V}(t_2) \cdots \hat{V}(t_n) | \rangle_0}{_0 \langle | S(\infty,-\infty) | \rangle_0} \\&= \sum_{n=0}^\infty \dfrac{(-i)^{n+1}}{n!} \int_{-\infty}^\infty dt_1 \cdots \int_{-\infty}^\infty dt_n \dfrac{_0 \langle | T \hat{C}_{\bm{p}}(t) \hat{V}(t_1) \hat{V}(t_2) \cdots \hat{V}(t_n) \hat{C}^\dagger_{\bm{p}}(t^\prime) | \rangle_0}{_0 \langle | S(\infty,-\infty) | \rangle_0}\tag{7}

となります.何度か言っていますが,こうやってシンプルに書けるのは時間順序積の強みです.ですが,書けるだけで全く計算の手間は減っていないのです.なお,最後に変形で時間順序積内でもフェルミオンの演算子の交換の際には符号が変化しなければなりませんが, \hat{C}^\dagger_{\bm{p}}(t^\prime) を符号変化なしで最後に持って行けたのは, \hat{V} が偶数個の演算子を含んでいるからだと思います.

ここで \hat{V}(t_1) をこの記事の最初に出した電子-電子相互作用にします.

\hat{V}(t_1) &= \dfrac{1}{2}\sum_{\bm{k}\bm{k^\prime}\bm{q}}\sum_{ss^\prime} \dfrac{4 \pi e^2}{q^2} C^\dagger_{\bm{k}+\bm{q},s} C^\dagger_{\bm{k^\prime}-\bm{q},s^\prime} C_{\bm{k^\prime},s^\prime} C_{\bm{k},s} \exp(i t_1 (\xi_{\bm{k}+\bm{q}}+\xi_{\bm{k^\prime}-\bm{q}}-\xi_{\bm{k^\prime}}-\xi_{\bm{k}}))\tag{8}

すると例えば,

_0 \langle | T \hat{C}_{\bm{p}}(t) \hat{V}(t_1) \hat{V}(t_2) \hat{C}^\dagger_{\bm{p}}(t^\prime) | \rangle_0\tag{9}

等の計算が必要になります.この場合,一つの相互作用の \hat{V}(t) には二つずつ生成演算子と消滅演算子 が入っているので,式 (9) には五個ずつの生成演算子と消滅演算子が入っていることになります.

ここで振り返りをしておきます.自由な電子グリーン関数は,

_0 \langle | T \hat{C}_\lambda(t) \hat{C}^\dagger_\lambda(t^\prime) | \rangle_0 = i G^{(0)}(\lambda,t-t^\prime)  \tag{10}

でした.( 電子とフォノンのグリーン関数 の式 (4)

一方,時間並進で系が不変と考えるのは自然な事だと思いますが,

それを考えると,同時刻の

_0 \langle | \hat{C}^\dagger_\lambda(t) \hat{C}_\lambda(t) | \rangle_0 &= n_F(\xi_{\bm{k}}) \\_0 \langle | \hat{C}_\lambda(t) \hat{C}^\dagger_\lambda(t) | \rangle_0 &= n_F(-\xi_{\bm{k}}) \tag{11}

となります.( 電子とフォノンのグリーン関数 の式 (14) )

つまり,生成消滅演算子のペアは同時刻か異時刻で値の種類が変わるのです. また重要な注意として, H_0 の基底状態 | \rangle_0 で挟まれる演算子は 生成と消滅がそれぞれ同数でないと期待値がゼロになってしまいます. これは粒子数の異なる波動関数が直交する為です.

ウィックの定理

さあ,予告通り相互作用のあるグリーン関数を自由なグリーン関数で展開しましょう. それをウィックの定理と言います.証明は飛ばします.

【ルール1】

まず,量子数 \lambda が異なる演算子のペアの期待値はゼロになります. つまり, \delta_{\alpha \beta} が必ず付いて

_0 \langle |T \hat{C}_\alpha(t) \hat{C}^\dagger_\beta(t^\prime) | \rangle_0 = \begin{cases}  i \delta_{\alpha \beta} G^{(0)}(\alpha, t-t^\prime) \\  \delta_{\alpha \beta} n_F(\xi_{\bm{k}}) \ \ \ \ \ (t=t^\prime+0) \\  \delta_{\alpha \beta} n_F(-\xi_{\bm{k}}) \ \ \ \ \ (t=t^\prime-0) \end{cases}\tag{12}

となります.生成消滅演算子の非ゼロ値を出す組み合わせはトータルで偶数個ですから,次は四つの展開です.

_0 \langle |T \hat{C}_\alpha(t) \hat{C}^\dagger_\beta(t_1) \hat{C}_\gamma(t_2) \hat{C}^\dagger_\delta(t^\prime) | \rangle_0\tag{13}

これは \alpha = \beta, \gamma = \delta または \alpha = \delta, \beta = \gamma 以外の時はゼロになります. 更に,もう一歩踏み込んで四次の期待値はそれぞれに時間順序積が付いた二次の期待値の多項式に分解できます. その際,添え字 \alpha \sim \delta の交換はフェルミオン演算子の交換の際の符号変化の約束を守ります. つまり,

&_0 \langle |T \hat{C}_\alpha(t) \hat{C}^\dagger_\beta(t_1) \hat{C}_\gamma(t_2) \hat{C}^\dagger_\delta(t^\prime) | \rangle_0 \\&= \delta_{\alpha \beta} \delta_{\gamma \delta} \ _0 \langle |T \hat{C}_\alpha(t) \hat{C}^\dagger_\beta(t_1)| \rangle_0 \ _0\langle |T \hat{C}_\gamma(t_2) \hat{C}^\dagger_\delta(t^\prime) | \rangle_0 \\&- \delta_{\alpha \delta} \delta_{\gamma \beta} \ _0 \langle |T \hat{C}_\alpha(t) \hat{C}^\dagger_\alpha(t^\prime)| \rangle_0 \ _0\langle |T \hat{C}_\gamma(t_2) \hat{C}^\dagger_\gamma(t_1) | \rangle_0 \tag{14}

このペアの数はまず n 個の消滅演算子を固定しましょう.そして,n 個の生成演算子を全ての並び替えで対応させます.つまり, n! 個の項が出てくるわけです.

【ルール2】

電子とフォノン等の異種の演算子が混在しているとき,同種の交換はその粒子の性質に応じ,異種は常に符号変化なく交換する.そして異種の演算子が対になることは無いです.

&_0 \langle |T \hat{C}_{\bm{p}}(t) \hat{C}^\dagger_{\bm{p}_1}(t_1) \hat{A}_{\bm{q}_1}(t_1) \hat{C}_{\bm{p}_2}(t_2) \hat{C}^\dagger_{\bm{p}_3}(t_3) \hat{A}_{\bm{q}_2}(t_2) | \rangle_0 \\&= _0 \langle |T \hat{C}_{\bm{p}}(t) \hat{C}^\dagger_{\bm{p}_1}(t_1) \hat{C}_{\bm{p}_2}(t_2) \hat{C}^\dagger_{\bm{p}_3}(t_3) | \rangle_0 \ _0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_2}(t_2) | \rangle_0\tag{15}

【ルール3】

フォノン演算子もペアリングできます.交換で符号は変わりません.

&_0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_2}(t_2) \hat{A}_{\bm{q}_3}(t_3) \hat{A}_{\bm{q}_4}(t_4) | \rangle_0 \\&= _0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_2}(t_2) | \rangle_0 \ _0 \langle |T \hat{A}_{\bm{q}_3}(t_3) \hat{A}_{\bm{q}_4}(t_4) | \rangle_0 \\&+ _0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_3}(t_3) | \rangle_0 \ _0 \langle |T \hat{A}_{\bm{q}_2}(t_2) \hat{A}_{\bm{q}_4}(t_4) | \rangle_0 \\&+ _0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_4}(t_4) | \rangle_0 \ _0 \langle |T \hat{A}_{\bm{q}_2}(t_2) \hat{A}_{\bm{q}_3}(t_3) | \rangle_0 \tag{16}

この組み合わせの数は 2n 個の演算子があるとして,

\left( \prod_{i=1}^n \ _{2n-2(i-1)} C _{2} \right)/n!\tag{17}

個の項が出てくると思います. 2n=2_2 C _2/1!=12n=4_4C_2 \ _2C_2/2!=3 これは合っていそうです.

そして,物理的に可能なのは \bm{q}_1 のフォノンを作って \bm{q}_2 のフォノンを消すと言う事は,フォノンには生成と消滅の区別が無いため \bm{q}_2 = -\bm{q}_1 でなければなりません.結局は,

&_0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_2}(t_2) \hat{A}_{\bm{q}_3}(t_3) \hat{A}_{\bm{q}_4}(t_4) | \rangle_0 \\&= \delta_{\bm{q}_1 + \bm{q}_2} \delta_{\bm{q}_3 + \bm{q}_4} \ _0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{-\bm{q}_1}(t_2) | \rangle_0 \ _0 \langle |T \hat{A}_{\bm{q}_3}(t_3) \hat{A}_{-\bm{q}_3}(t_4) | \rangle_0 \\&+ \delta_{\bm{q}_1 + \bm{q}_3} \delta_{\bm{q}_2 + \bm{q}_4} \ _0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{-\bm{q}_1}(t_3) | \rangle_0 \ _0 \langle |T \hat{A}_{\bm{q}_2}(t_2) \hat{A}_{-\bm{q}_2}(t_4) | \rangle_0 \\&+ \delta_{\bm{q}_1 + \bm{q}_4} \delta_{\bm{q}_2 + \bm{q}_3} \ _0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{-\bm{q}_1}(t_4) | \rangle_0 \ _0 \langle |T \hat{A}_{\bm{q}_2}(t_2) \hat{A}_{-\bm{q}_2}(t_3) | \rangle_0 \tag{18}

となります.

【ルール4】

同時刻の演算子の期待値は粒子数演算子になります.

_0 \langle | &T \hat{C}^\dagger_{\bm{k}_1}(t_1) \hat{C}_{\bm{k}_2}(t_1)| \rangle_0 \\&= \delta_{\bm{k}_1 = \bm{k}_2} \ _0 \langle | \hat{C}^\dagger_{\bm{k}_1}(t_1) \hat{C}_{\bm{k}_1}(t_1)| \rangle_0 \\&= \delta_{\bm{k}_1 = \bm{k}_2} n_F(\xi_{\bm{k}_1})\tag{19}

ボゾンの場合は,

n_B(\omega_{\bm{k}}) = \dfrac{1}{e^{\beta \omega_k}-1} \tag{20}

【ルール5】

異なる時刻の演算子の期待値はグリーン関数になります. 慣習として,生成演算子が右側に来るようにします. その際,符号変化はお忘れず.

_0 \langle | &T \hat{C}^\dagger_{\bm{k}_1}(t_1) \hat{C}_{\bm{k}_2}(t_2)| \rangle_0 \\&= - \delta_{\bm{k}_1 = \bm{k}_2} \ _0 \langle | T \hat{C}_{\bm{k}_2}(t_2) \hat{C}^\dagger_{\bm{k}_1}(t_1) | \rangle_0\tag{21}

具体例を出しておきましょう.先ほどの例で言えば,

&_0 \langle |T \hat{C}_\alpha(t) \hat{C}^\dagger_\beta(t_1) \hat{C}_\gamma(t_2) \hat{C}^\dagger_\delta(t^\prime) | \rangle_0 \\&= i^2 \delta_{\alpha \beta} \delta_{\gamma \delta} G^{(0)}(\alpha,t-t_1) G^{(0)}(\alpha,t_2-t^\prime) \\&- i^2 \delta_{\alpha \delta} \delta_{\gamma \beta} G^{(0)}(\alpha,t-t^\prime) G^{(0)}(\alpha,t_2-t_1)\tag{22}

であり,

&_0 \langle |T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_2}(t_2) \hat{A}_{\bm{q}_3}(t_3) \hat{A}_{\bm{q}_4}(t_4) | \rangle_0 \\&= \delta_{\bm{q}_1 + \bm{q}_2} \delta_{\bm{q}_3 + \bm{q}_4} D^{(0)}(\bm{q}_1,t_1 - t_2) D^{(0)}(\bm{q}_3,t_3 - t_4) \\&+ \delta_{\bm{q}_1 + \bm{q}_3} \delta_{\bm{q}_2 + \bm{q}_4} D^{(0)}(\bm{q}_1,t_1 - t_3) D^{(0)}(\bm{q}_2,t_2 - t_4) \\&+ \delta_{\bm{q}_1 + \bm{q}_4} \delta_{\bm{q}_2 + \bm{q}_3} D^{(0)}(\bm{q}_1,t_1 - t_4) D^{(0)}(\bm{q}_2,t_2 - t_3) \tag{23}

となります.この様にウィックの定理は時間順序された高次の期待値が全ての可能な時間順序されたペアリングによって表せると言う事を言っています.このウィックの定理が有効なのは H_0 に含まれる生成消滅演算子が一個ずつの時だけです.それにはご注意を.

ウィックの定理の実践例(電子フォノン相互作用)

電子-フォノン相互作用の場合にこれを用いて近似値を出してみましょう. ただし,グリーン関数の分母にある位相因子の _0\langle | S(\infty, -\infty) | \rangle_0 は無視します. 式 (2) を使います.再掲すると,

\hat{V}(t) = \sum_{\bm{q}\bm{k}s} M_{\bm{q}} A_{\bm{q}} C^\dagger_{\bm{k}+\bm{q},s} C_{\bm{k},s}\tag{24}

これを使います.すると,2次までの近似で

G(\bm{p},t-t^\prime) &= G^{(0)}(\bm{p},t-t^\prime) + (-i)^2 \int_{-\infty}^\infty dt_1 \ _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{V}(t_1) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \\&+ \dfrac{(-i)^3}{2!} \int_{-\infty}^\infty dt_1 \int_{-\infty}^\infty dt_2 \ _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{V}(t_1) \hat{V}(t_2) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0\tag{25}

ここで,右辺第二項のフォノンの1次は期待値がゼロになるので

_0 \langle |a_{\bm{q}}| \rangle_0 = _0 \langle |a^\dagger_{\bm{q}}| \rangle_0 = 0\tag{26}

同様の理由で奇数次のボゾンの期待値はゼロです. 計算を進めると,

G(\bm{p},t-t^\prime) &= G^{(0)}(\bm{p},t-t^\prime) + \dfrac{(-i)^3}{2!} \int_{-\infty}^\infty dt_1 \int_{-\infty}^\infty dt_2 \sum_{\bm{q}_1 \bm{q}_2} M_{\bm{q}_1} M_{\bm{q}_2} \ _0 \langle | T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_2}(t_2) | \rangle_0 \\&\times \sum_{\bm{k}_1 \bm{k}_2 s s^\prime} \ _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) \hat{C}_{\bm{k}_1,s}(t_1) \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) \hat{C}_{\bm{k}_2,s^\prime}(t_2) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0\tag{27}

二次近似項のフォノン部分は,

&_0 \langle | T \hat{A}_{\bm{q}_1}(t_1) \hat{A}_{\bm{q}_2}(t_2) | \rangle_0 \\&= i \delta_{\bm{q}_1+\bm{q}_2} D^{(0)}(\bm{q}_1,t_1-t_2)\tag{28}

二次近似項の電子部分は,

&_0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) \hat{C}_{\bm{k}_1,s}(t_1) \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) \hat{C}_{\bm{k}_2,s^\prime}(t_2) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \\ \\&= _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) | \rangle_0 \ _0 \langle | T \hat{C}_{\bm{k}_1,s}(t_1) \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) | \rangle_0 \ _0 \langle | T \hat{C}_{\bm{k}_2,s^\prime}(t_2) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \\ \\&+ _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s^\prime}(t_2) | \rangle_0 \ _0 \langle | T \hat{C}_{\bm{k}_1,s}(t_1) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \ _0 \langle | T \hat{C}_{\bm{k}_2,s^\prime}(t_2) \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) | \rangle_0  \\ \\&+ _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) | \rangle_0 \ _0 \langle | T \hat{C}_{\bm{k}_1,s}(t_1) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \ _0 \langle | T \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) \hat{C}_{\bm{k}_2,s^\prime}(t_2) | \rangle_0  \\ \\&+ _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) | \rangle_0 \ _0 \langle | T \hat{C}_{\bm{k}_2,s^\prime}(t_2) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \ _0 \langle | T \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) \hat{C}_{\bm{k}_1,s}(t_1) | \rangle_0  \\ \\&+ _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \ _0 \langle | T \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) \hat{C}_{\bm{k}_1,s}(t_1) | \rangle_0 \ _0 \langle | T \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) \hat{C}_{\bm{k}_2,s^\prime}(t_2) | \rangle_0  \\ \\&- _0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \ _0 \langle | T  \hat{C}_{\bm{k}_1,s}(t_1) \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) | \rangle_0 \ _0 \langle | T \hat{C}_{\bm{k}_2,s^\prime}(t_2) \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) | \rangle_0\tag{29}

長いので一旦切ります.式はさらに変形でき,

&_0 \langle | T \hat{C}_{\bm{p} \sigma}(t) \hat{C}^\dagger_{\bm{k}_1+\bm{q}_1, s}(t_1) \hat{C}_{\bm{k}_1,s}(t_1) \hat{C}^\dagger_{\bm{k}_2+\bm{q}_2, s}(t_2) \hat{C}_{\bm{k}_2,s^\prime}(t_2) \hat{C}^\dagger_{\bm{p} \sigma}(t^\prime) | \rangle_0 \\ \\&= i^3 \delta_{\bm{p} = \bm{k}_2 = \bm{k}_1 + \bm{q}_1} \delta_{s=s^\prime=\sigma} G^{(0)}(\bm{p},t - t_1) G^{(0)}(\bm{p} - \bm{q}_1,t_1 - t_2) G^{(0)}(\bm{p},t_2 - t^\prime) \\ \\&+ i^3 \delta_{\bm{p} = \bm{k}_1 = \bm{k}_2 - \bm{q}_1} \delta_{s=s^\prime=\sigma} G^{(0)}(\bm{p},t - t_2) G^{(0)}(\bm{p} - \bm{q}_1,t_2 - t_1) G^{(0)}(\bm{p},t_1 - t^\prime) \\ \\&+ i^2 \delta_{\bm{q}_1 = 0} \delta_{\bm{p} = \bm{k}_1} \delta_{s=\sigma} n_F(\xi_{\bm{k_2}}) G^{(0)}(\bm{p},t - t_1) G^{(0)}(\bm{p},t_1 - t^\prime) \\ \\&+ i^2 \delta_{\bm{q}_1 = 0} \delta_{\bm{p} = \bm{k}_2} \delta_{s^\prime=\sigma} n_F(\xi_{\bm{k_1}}) G^{(0)}(\bm{p},t - t_2) G^{(0)}(\bm{p},t_2 - t^\prime) \\ \\&+ i \delta_{\bm{q}_1 = 0} \delta_{\bm{q}_2 = 0} n_F(\xi_{\bm{k_1}})n_F(\xi_{\bm{k_2}}) G^{(0)}(\bm{p},t - t^\prime) \\ \\&- i^3 \delta_{\bm{k}_1 = \bm{k}_2 - \bm{q}_1} \delta_{s^\prime = s} G^{(0)}(\bm{p},t-t^\prime) G^{(0)}(\bm{k}_1,t_1-t_2) G^{(0)}(\bm{k}_1 + \bm{q}_1,t_2-t_1)\tag{30}

となります.この一連の変形は項の順番は入れ替えておらず,順に対応しています. 次回はこれをファインマンダイアグラムという図と対応させて理解します.

今日はここまで,お疲れ様でした.

次の記事は ファインマンダイアグラム です.