S η = 1 β ∑ i ω n g ( i ω n ) {\displaystyle S_{\eta }={\frac {1}{\beta }}\sum _{i\omega _{n}}g(i\omega _{n})} Figure 1. Figure 2. 松原振動数の和を評価する上手なやり方は、 z = i ω {\displaystyle z=i\omega } に極 を持つ松原重み関数h η (z )を使う方法である。 ボソンの場合η = +1とフェルミオンの場合η = −1で重み関数は異なる。 重み関数の選択について後述する。 和は、重み関数を使って複素平面での閉曲線積分に置き換えることができる。
S η = 1 β ∑ i ω g ( i ω ) = 1 2 π i β ∮ g ( z ) h η ( z ) d z {\displaystyle S_{\eta }={\frac {1}{\beta }}\sum _{i\omega }g(i\omega )={\frac {1}{2\pi i\beta }}\oint g(z)h_{\eta }(z)\,dz} Fig. 1において、重み関数は虚数軸上に極(赤バツ印)を作る。 閉曲線積分はこれらの極の留数 をピックアップし、これは和に等しい。 閉曲線をg (z )の極(Fig. 2、緑バツ印)を囲むように変形すると、和はg (z )h η (z )の留数の全ての極g (z )についての和によって形式的に遂行される。
S η = − 1 β ∑ z 0 ∈ g ( z ) poles Res g ( z 0 ) h η ( z 0 ) {\displaystyle S_{\eta }=-{\frac {1}{\beta }}\sum _{z_{0}\in g(z){\text{ poles}}}\operatorname {Res} g(z_{0})h_{\eta }(z_{0})} ここで閉曲線が極を時計回りの方向で囲むように変形し、負の留数を生むため、マイナスがつくことに注意。
松原重み関数の選択 ボソン振動数 z = i ω n {\displaystyle z=i\omega _{n}} の極を作るために、どちらの半平面で収束がコントロールされるかに依存して次の2つのタイプの松原重み関数を選ぶことができる。
h B ( 1 ) ( z ) = β 1 − e − β z = − β n B ( − z ) = β ( 1 + n B ( z ) ) {\displaystyle h_{B}^{(1)}(z)={\frac {\beta }{1-e^{-\beta z}}}=-\beta n_{B}(-z)=\beta (1+n_{B}(z))} h B ( 2 ) ( z ) = − β 1 − e β z = β n B ( z ) {\displaystyle h_{B}^{(2)}(z)={\frac {-\beta }{1-e^{\beta z}}}=\beta n_{B}(z)} h B ( 1 ) ( z ) {\displaystyle h_{B}^{(1)}(z)} は左半平面(Re z < 0)での収束をコントロールし、 h B ( 2 ) ( z ) {\displaystyle h_{B}^{(2)}(z)} は右半平面(Re z > 0)での収束をコントロールする。 ここで n B ( z ) = ( e β z − 1 ) − 1 {\displaystyle n_{B}(z)=(e^{\beta z}-1)^{-1}} はボース分布関数 である。
フェルミオン振動数の場合も同様である。 2つのタイプの松原重み関数があり、 z = i ω m {\displaystyle z=i\omega _{m}} に極を作る。
h F ( 1 ) ( z ) = β 1 + e − β z = β n F ( − z ) = β ( 1 − n F ( z ) ) {\displaystyle h_{F}^{(1)}(z)={\frac {\beta }{1+e^{-\beta z}}}=\beta n_{F}(-z)=\beta (1-n_{F}(z))} h F ( 2 ) ( z ) = − β 1 + e β z = − β n F ( z ) {\displaystyle h_{F}^{(2)}(z)={\frac {-\beta }{1+e^{\beta z}}}=-\beta n_{F}(z)} h F ( 1 ) ( z ) {\displaystyle h_{F}^{(1)}(z)} は左半平面(Re z < 0)での収束をコントロールし、 h F ( 1 ) ( z ) {\displaystyle h_{F}^{(1)}(z)} は右半平面(Re z > 0)での収束をコントロールする。 ここで n F ( z ) = ( e β z + 1 ) − 1 {\displaystyle n_{F}(z)=(e^{\beta z}+1)^{-1}} はフェルミ分布関数 である。
グリーン関数への応用では、g (z )は常に次の構造を持つ。
g ( z ) = G ( z ) e − z τ {\displaystyle g(z)=G(z)e^{-z\tau }} これは0 < τ < β で与えられる左半平面で発散する。 収束をコントロールするために、第一のタイプの重み関数は常に h η ( z ) = h η ( 1 ) ( z ) {\displaystyle h_{\eta }(z)=h_{\eta }^{(1)}(z)} と選ぶ。 しかし松原振動数の和が発散しないときは収束をコントロールする必要はない。 そのような場合、どんな松原重み関数を選んでも同じ結果が得られる。
松原振動数の和の表 以下の表に、 いくつかの簡単な有理関数 g (z )での松原振動数の和をまとめる。
S η = 1 β ∑ i ω g ( i ω ) . {\displaystyle S_{\eta }={\frac {1}{\beta }}\sum _{i\omega }g(i\omega ).} η = ±1 は統計的記号である。
g ( i ω ) {\displaystyle g(i\omega )} S η {\displaystyle S_{\eta }} ( i ω − ξ ) − 1 {\displaystyle (i\omega -\xi )^{-1}} − η n η ( ξ ) {\displaystyle -\eta n_{\eta }(\xi )} [1] ( i ω − ξ ) − 2 {\displaystyle (i\omega -\xi )^{-2}} − η n η ′ ( ξ ) = β n η ( ξ ) ( η + n η ( ξ ) ) {\displaystyle -\eta n_{\eta }^{\prime }(\xi )=\beta n_{\eta }(\xi )(\eta +n_{\eta }(\xi ))} ( i ω − ξ ) − n {\displaystyle (i\omega -\xi )^{-n}} − η ( n − 1 ) ! ∂ ξ n − 1 n η ( ξ ) {\displaystyle -{\frac {\eta }{(n-1)!}}\partial _{\xi }^{n-1}n_{\eta }(\xi )} 1 ( i ω − ξ 1 ) ( i ω − ξ 2 ) {\displaystyle {\frac {1}{(i\omega -\xi _{1})(i\omega -\xi _{2})}}} − η ( n η ( ξ 1 ) − n η ( ξ 2 ) ) ξ 1 − ξ 2 {\displaystyle -{\frac {\eta (n_{\eta }(\xi _{1})-n_{\eta }(\xi _{2}))}{\xi _{1}-\xi _{2}}}} 1 ( i ω − ξ 1 ) 2 ( i ω − ξ 2 ) 2 {\displaystyle {\frac {1}{(i\omega -\xi _{1})^{2}(i\omega -\xi _{2})^{2}}}} η ( ξ 1 − ξ 2 ) 2 ( 2 ( n η ( ξ 1 ) − n η ( ξ 2 ) ) ξ 1 − ξ 2 − ( n η ′ ( ξ 1 ) + n η ′ ( ξ 2 ) ) ) {\displaystyle {\frac {\eta }{(\xi _{1}-\xi _{2})^{2}}}\left({\frac {2(n_{\eta }(\xi _{1})-n_{\eta }(\xi _{2}))}{\xi _{1}-\xi _{2}}}-(n_{\eta }^{\prime }(\xi _{1})+n_{\eta }^{\prime }(\xi _{2}))\right)} 1 ( i ω − ξ 1 ) 2 − ξ 2 2 {\displaystyle {\frac {1}{(i\omega -\xi _{1})^{2}-\xi _{2}^{2}}}} η c η ( ξ 1 , ξ 2 ) {\displaystyle \eta c_{\eta }(\xi _{1},\xi _{2})} 1 ( i ω ) 2 − ξ 2 {\displaystyle {\frac {1}{(i\omega )^{2}-\xi ^{2}}}} η c η ( 0 , ξ ) = − 1 2 ξ ( 1 + 2 η n η ( ξ ) ) {\displaystyle \eta c_{\eta }(0,\xi )=-{\frac {1}{2\xi }}(1+2\eta n_{\eta }(\xi ))} ( i ω ) 2 ( i ω ) 2 − ξ 2 {\displaystyle {\frac {(i\omega )^{2}}{(i\omega )^{2}-\xi ^{2}}}} − ξ 2 ( 1 + 2 η n η ( ξ ) ) {\displaystyle -{\frac {\xi }{2}}(1+2\eta n_{\eta }(\xi ))} 1 ( ( i ω ) 2 − ξ 2 ) 2 {\displaystyle {\frac {1}{((i\omega )^{2}-\xi ^{2})^{2}}}} − η 2 ξ 2 ( c η ( 0 , ξ ) + n η ′ ( ξ ) ) {\displaystyle -{\frac {\eta }{2\xi ^{2}}}(c_{\eta }(0,\xi )+n_{\eta }^{\prime }(\xi ))} ( i ω ) 2 ( ( i ω ) 2 − ξ 2 ) 2 {\displaystyle {\frac {(i\omega )^{2}}{((i\omega )^{2}-\xi ^{2})^{2}}}} η 2 ( c η ( 0 , ξ ) − n η ′ ( ξ ) ) {\displaystyle {\frac {\eta }{2}}(c_{\eta }(0,\xi )-n_{\eta }^{\prime }(\xi ))} ( i ω ) 2 + ξ 2 ( ( i ω ) 2 − ξ 2 ) 2 {\displaystyle {\frac {(i\omega )^{2}+\xi ^{2}}{((i\omega )^{2}-\xi ^{2})^{2}}}} − η n η ′ ( ξ ) = β n η ( ξ ) ( η + n η ( ξ ) ) {\displaystyle -\eta n_{\eta }^{\prime }(\xi )=\beta n_{\eta }(\xi )(\eta +n_{\eta }(\xi ))} 1 ( ( i ω ) 2 − ξ 1 2 ) ( ( i ω ) 2 − ξ 2 2 ) {\displaystyle {\frac {1}{((i\omega )^{2}-\xi _{1}^{2})((i\omega )^{2}-\xi _{2}^{2})}}} η ( c η ( 0 , ξ 1 ) − c η ( 0 , ξ 2 ) ) ξ 1 2 − ξ 2 2 {\displaystyle {\frac {\eta (c_{\eta }(0,\xi _{1})-c_{\eta }(0,\xi _{2}))}{\xi _{1}^{2}-\xi _{2}^{2}}}} ( 1 ( i ω ) 2 − ξ 1 2 + 1 ( i ω ) 2 − ξ 2 2 ) 2 {\displaystyle \left({\frac {1}{(i\omega )^{2}-\xi _{1}^{2}}}+{\frac {1}{(i\omega )^{2}-\xi _{2}^{2}}}\right)^{2}} η ( 3 ξ 1 2 + ξ 2 2 2 ξ 1 2 ( ξ 1 2 − ξ 2 2 ) c η ( 0 , ξ 1 ) − n η ′ ( ξ 1 ) 2 ξ 1 2 ) + ( 1 ↔ 2 ) {\displaystyle \eta \left({\frac {3\xi _{1}^{2}+\xi _{2}^{2}}{2\xi _{1}^{2}(\xi _{1}^{2}-\xi _{2}^{2})}}c_{\eta }(0,\xi _{1})-{\frac {n_{\eta }^{\prime }(\xi _{1})}{2\xi _{1}^{2}}}\right)+(1\leftrightarrow 2)} [2] ( 1 ( i ω ) 2 − ξ 1 2 − 1 ( i ω ) 2 − ξ 2 2 ) 2 {\displaystyle \left({\frac {1}{(i\omega )^{2}-\xi _{1}^{2}}}-{\frac {1}{(i\omega )^{2}-\xi _{2}^{2}}}\right)^{2}} η ( − 5 ξ 1 2 − ξ 2 2 2 ξ 1 2 ( ξ 1 2 − ξ 2 2 ) c η ( 0 , ξ 1 ) − n η ′ ( ξ 1 ) 2 ξ 1 2 ) + ( 1 ↔ 2 ) {\displaystyle \eta \left(-{\frac {5\xi _{1}^{2}-\xi _{2}^{2}}{2\xi _{1}^{2}(\xi _{1}^{2}-\xi _{2}^{2})}}c_{\eta }(0,\xi _{1})-{\frac {n_{\eta }^{\prime }(\xi _{1})}{2\xi _{1}^{2}}}\right)+(1\leftrightarrow 2)} [2]
[1] 和は収束しないため、松原重み関数の選択が異なると結果は定数分だけ異なる。
[2] (1 ↔ 2)は、手前の項のインデックス1と2を置き換えたものを表す。
物理学での応用 温度ゼロの極限 極限 β → ∞ {\displaystyle \beta \rightarrow \infty } での松原振動数の和は、虚数振動数の虚軸についての積分に等しい
1 β ∑ i ω = ∫ − i ∞ i ∞ d ( i ω ) 2 π {\displaystyle {\frac {1}{\beta }}\sum _{i\omega }=\int _{-i\infty }^{i\infty }{\frac {\mathrm {d} (i\omega )}{2\pi }}} いくつかの積分は収束しない。 それらは振動数カットオフ Ω {\displaystyle \Omega } を導入し、 Ω → ∞ {\displaystyle \Omega \rightarrow \infty } の極限をとる前に積分から発散部分( Ω {\displaystyle \Omega } -依存)を差し引いて繰り込む必要がある。 例えば自由エネルギーは、対数の積分によって得られる。
η lim Ω → ∞ [ ∫ − i Ω i Ω d ( i ω ) 2 π ( ln ( − i ω + ξ ) − π ξ 2 Ω ) − Ω π ( ln Ω − 1 ) ] = { 0 ξ ≥ 0 , − η ξ ξ < 0 , {\displaystyle \eta \lim _{\Omega \rightarrow \infty }\left[\int _{-i\Omega }^{i\Omega }{\frac {\mathrm {d} (i\omega )}{2\pi }}\left(\ln(-i\omega +\xi )-{\frac {\pi \xi }{2\Omega }}\right)-{\frac {\Omega }{\pi }}(\ln \Omega -1)\right]=\left\{{\begin{array}{cc}0&\xi \geq 0,\\-\eta \xi &\xi <0,\end{array}}\right.} これは、温度ゼロにおいて自由エネルギーは化学ポテンシャルに満たない内部エネルギーと単純な関係にあることを意味している。 また分布関数は次の積分によって得られる。
η lim Ω → ∞ ∫ − i Ω i Ω d ( i ω ) 2 π ( 1 − i ω + ξ − π 2 Ω ) = { 0 ξ ≥ 0 , − η ξ < 0 , {\displaystyle \eta \lim _{\Omega \rightarrow \infty }\int _{-i\Omega }^{i\Omega }{\frac {\mathrm {d} (i\omega )}{2\pi }}\left({\frac {1}{-i\omega +\xi }}-{\frac {\pi }{2\Omega }}\right)=\left\{{\begin{array}{cc}0&\xi \geq 0,\\-\eta &\xi <0,\end{array}}\right.} これは温度ゼロでの階段関数のふるまいを示している。
グリーン関数との関連 時間領域 虚時間区間(0,β )で定義される関数G (τ )を考える。 これはフーリエ級数の観点で与えられる。
G ( τ ) = 1 β ∑ i ω G ( i ω ) e − i ω τ , {\displaystyle G(\tau )={\frac {1}{\beta }}\sum _{i\omega }G(i\omega )e^{-i\omega \tau },} ここで振動数は 2π /β 間隔の離散的な値のみとる。
振動数の選択は、関数G (τ )の境界条件に依存している。 物理学ではG (τ )はグリーン関数の虚時間表現を表す。
G ( τ ) = − ⟨ T τ ψ ( τ ) ψ ∗ ( 0 ) ⟩ . {\displaystyle G(\tau )=-\langle {\mathcal {T}}_{\tau }\psi (\tau )\psi ^{*}(0)\rangle .} これはボソン場の周期的境界条件G (τ +β )=G (τ )を満たす。 一方フェルミオン場では、境界条件は反周期的G (τ + β ) = −G (τ )である。
振動数領域でのグリーン関数G (iω )が与えられたとき、その虚時間表現G (τ )は松原振動数の和によって評価できる。 その和がボソン振動数かフェルミオン振動数のどちらでとるかに依存して、得られるG (τ )は異なる。 これらを区別するため、次を定義する。
G η ( τ ) = { G B ( τ ) , if η = + 1 , G F ( τ ) , if η = − 1 , {\displaystyle G_{\eta }(\tau )={\begin{cases}G_{B}(\tau ),&{\text{if }}\eta =+1,\\G_{F}(\tau ),&{\text{if }}\eta =-1,\end{cases}}} G B ( τ ) = 1 β ∑ i ω n G ( i ω n ) e − i ω n τ , {\displaystyle G_{B}(\tau )={\frac {1}{\beta }}\sum _{i\omega _{n}}G(i\omega _{n})e^{-i\omega _{n}\tau },} G F ( τ ) = 1 β ∑ i ω m G ( i ω m ) e − i ω m τ . {\displaystyle G_{F}(\tau )={\frac {1}{\beta }}\sum _{i\omega _{m}}G(i\omega _{m})e^{-i\omega _{m}\tau }.} ここでτ は区間(0,β )に制限されていることに注意。 境界条件は区間の外にG (τ )を拡張するために用いることができる。 よく用いられる結果を以下の表にまとめる。
G ( i ω ) {\displaystyle G(i\omega )} G η ( τ ) {\displaystyle G_{\eta }(\tau )} ( i ω − ξ ) − 1 {\displaystyle (i\omega -\xi )^{-1}} − e ξ ( β − τ ) n η ( ξ ) {\displaystyle -e^{\xi (\beta -\tau )}n_{\eta }(\xi )} ( i ω − ξ ) − 2 {\displaystyle (i\omega -\xi )^{-2}} e ξ ( β − τ ) n η ( ξ ) ( τ + η β n η ( ξ ) ) {\displaystyle e^{\xi (\beta -\tau )}n_{\eta }(\xi )\left(\tau +\eta \beta n_{\eta }(\xi )\right)} ( i ω − ξ ) − 3 {\displaystyle (i\omega -\xi )^{-3}} − 1 2 e ξ ( β − τ ) n η ( ξ ) ( τ 2 + η β ( β + 2 τ ) n η ( ξ ) + 2 β 2 n η 2 ( ξ ) ) {\displaystyle -{\frac {1}{2}}e^{\xi (\beta -\tau )}n_{\eta }(\xi )\left(\tau ^{2}+\eta \beta (\beta +2\tau )n_{\eta }(\xi )+2\beta ^{2}n_{\eta }^{2}(\xi )\right)} ( i ω − ξ 1 ) − 1 ( i ω − ξ 2 ) − 1 {\displaystyle (i\omega -\xi _{1})^{-1}(i\omega -\xi _{2})^{-1}} − e ξ 1 ( β − τ ) n η ( ξ 1 ) − e ξ 2 ( β − τ ) n η ( ξ 2 ) ξ 1 − ξ 2 {\displaystyle -{\frac {e^{\xi _{1}(\beta -\tau )}n_{\eta }(\xi _{1})-e^{\xi _{2}(\beta -\tau )}n_{\eta }(\xi _{2})}{\xi _{1}-\xi _{2}}}} ( ω 2 + m 2 ) − 1 {\displaystyle (\omega ^{2}+m^{2})^{-1}} e − m τ 2 m + η m cosh m τ n η ( m ) {\displaystyle {\frac {e^{-m\tau }}{2m}}+{\frac {\eta }{m}}\cosh {m\tau }\;n_{\eta }(m)} i ω ( ω 2 + m 2 ) − 1 {\displaystyle i\omega (\omega ^{2}+m^{2})^{-1}} e − m τ 2 − η sinh m τ n η ( m ) {\displaystyle {\frac {e^{-m\tau }}{2}}-\eta \,\sinh {m\tau }\;n_{\eta }(m)}
演算子スイッチング効果 ここでは小さな虚時間が決定的な役割を果たす。 小さな虚時間の負号が変わると、演算子の順番が変わる
⟨ ψ ψ ∗ ⟩ = ⟨ T τ ψ ( τ = 0 + ) ψ ∗ ( 0 ) ⟩ = − G η ( τ = 0 + ) = − 1 β ∑ i ω G ( i ω ) e − i ω 0 + {\displaystyle \langle \psi \psi ^{*}\rangle =\langle {\mathcal {T}}_{\tau }\psi (\tau =0^{+})\psi ^{*}(0)\rangle =-G_{\eta }(\tau =0^{+})=-{\frac {1}{\beta }}\sum _{i\omega }G(i\omega )e^{-i\omega 0^{+}}} .\ ⟨ ψ ∗ ψ ⟩ = η ⟨ T τ ψ ( τ = 0 − ) ψ ∗ ( 0 ) ⟩ = − η G η ( τ = 0 − ) = − η β ∑ i ω G ( i ω ) e i ω 0 + {\displaystyle \langle \psi ^{*}\psi \rangle =\eta \langle {\mathcal {T}}_{\tau }\psi (\tau =0^{-})\psi ^{*}(0)\rangle =-\eta G_{\eta }(\tau =0^{-})=-{\frac {\eta }{\beta }}\sum _{i\omega }G(i\omega )e^{i\omega 0^{+}}} 分布関数 τ = 0でのグリーン関数G (τ )の不連続性のため、分布関数の評価は難しくなる。 次の和を評価するためにボソンとフェルミオンどちらの重み関数を選択することもできるが、結果は異なる。
G ( 0 ) = ∑ i ω ( i ω − ξ ) − 1 {\displaystyle G(0)=\sum _{i\omega }(i\omega -\xi )^{-1}} これは、G (τ )をτ = 0からわずかに遠ざけたとき、収束をコントロールするため、 G ( τ = 0 + ) {\displaystyle G(\tau =0^{+})} での重み関数として h η ( 1 ) ( z ) {\displaystyle h_{\eta }^{(1)}(z)} をとらなければならず、 G ( τ = 0 − ) {\displaystyle G(\tau =0^{-})} では h η ( 2 ) ( z ) {\displaystyle h_{\eta }^{(2)}(z)} をとらなければならないと理解できる。
ボソン
G B ( τ = 0 − ) = 1 β ∑ i ω n e i ω n 0 + i ω n − ξ = − n B ( ξ ) , {\displaystyle G_{B}(\tau =0^{-})={\frac {1}{\beta }}\sum _{i\omega _{n}}{\frac {e^{i\omega _{n}0^{+}}}{i\omega _{n}-\xi }}=-n_{B}(\xi ),} G B ( τ = 0 + ) = 1 β ∑ i ω n e − i ω n 0 + i ω n − ξ = − ( n B ( ξ ) + 1 ) . {\displaystyle G_{B}(\tau =0^{+})={\frac {1}{\beta }}\sum _{i\omega _{n}}{\frac {e^{-i\omega _{n}0^{+}}}{i\omega _{n}-\xi }}=-(n_{B}(\xi )+1).} フェルミオン
G F ( τ = 0 − ) = 1 β ∑ i ω m e i ω m 0 + i ω m − ξ = n F ( ξ ) , {\displaystyle G_{F}(\tau =0^{-})={\frac {1}{\beta }}\sum _{i\omega _{m}}{\frac {e^{i\omega _{m}0^{+}}}{i\omega _{m}-\xi }}=n_{F}(\xi ),} G F ( τ = 0 + ) = 1 β ∑ i ω m e − i ω m 0 + i ω m − ξ = − ( 1 − n F ( ξ ) ) . {\displaystyle G_{F}(\tau =0^{+})={\frac {1}{\beta }}\sum _{i\omega _{m}}{\frac {e^{-i\omega _{m}0^{+}}}{i\omega _{m}-\xi }}=-(1-n_{F}(\xi )).} 自由エネルギー ボソン
1 β ∑ i ω n ln ( β ( − i ω n + ξ ) ) = 1 β ln ( 1 − e − β ξ ) , {\displaystyle {\frac {1}{\beta }}\sum _{i\omega _{n}}\ln(\beta (-i\omega _{n}+\xi ))={\frac {1}{\beta }}\ln(1-e^{-\beta \xi }),} フェルミオン
− 1 β ∑ i ω m ln ( β ( − i ω m + ξ ) ) = − 1 β ln ( 1 + e − β ξ ) . {\displaystyle -{\frac {1}{\beta }}\sum _{i\omega _{m}}\ln(\beta (-i\omega _{m}+\xi ))=-{\frac {1}{\beta }}\ln(1+e^{-\beta \xi }).} ダイアグラムの評価 ここでは、よく使われる単一モードのダイアグラムを評価する。 多重モード問題は、スペクトル関数積分によってアプローチできる。
フェルミオン自己エネルギー Σ ( i ω m ) = − 1 β ∑ i ω n 1 i ω m + i ω n − ε 1 i ω n − Ω = n F ( ε ) − n F ( Ω ) i ω m − ε + Ω . {\displaystyle \Sigma (i\omega _{m})=-{\frac {1}{\beta }}\sum _{i\omega _{n}}{\frac {1}{i\omega _{m}+i\omega _{n}-\varepsilon }}{\frac {1}{i\omega _{n}-\Omega }}={\frac {n_{F}(\varepsilon )-n_{F}(\Omega )}{i\omega _{m}-\varepsilon +\Omega }}.} 粒子-空孔バブル Π ( i ω n ) = 1 β ∑ i ω m 1 i ω m + i ω n − ε 1 i ω m − ε ′ = − n F ( ε ) − n F ( ε ′ ) i ω n − ε + ε ′ . {\displaystyle \Pi (i\omega _{n})={\frac {1}{\beta }}\sum _{i\omega _{m}}{\frac {1}{i\omega _{m}+i\omega _{n}-\varepsilon }}{\frac {1}{i\omega _{m}-\varepsilon '}}=-{\frac {n_{F}(\varepsilon )-n_{F}\left(\varepsilon '\right)}{i\omega _{n}-\varepsilon +\varepsilon '}}.} 粒子-粒子バブル Π ( i ω n ) = − 1 β ∑ i ω m 1 i ω m + i ω n − ϵ 1 − i ω m − ϵ ′ = 1 − n F ( ϵ ) − n F ( ϵ ′ ) i ω n − ϵ − ϵ ′ . {\displaystyle \Pi (i\omega _{n})=-{\frac {1}{\beta }}\sum _{i\omega _{m}}{\frac {1}{i\omega _{m}+i\omega _{n}-\epsilon }}{\frac {1}{-i\omega _{m}-\epsilon '}}={\frac {1-n_{F}(\epsilon )-n_{F}\left(\epsilon '\right)}{i\omega _{n}-\epsilon -\epsilon '}}.} 付録: 分布関数の性質 分布関数 一般的表記 n η {\displaystyle n_{\eta }} はボース分布関数(η = +1)かフェルミ分布関数(η = −1)のどちらかを表す。
n η ( ξ ) = 1 e β ξ − η {\displaystyle n_{\eta }(\xi )={\frac {1}{e^{\beta \xi }-\eta }}} 必要であればボース分布関数とフェルミ分布関数を区別するために、それぞれ記号n B とn F を用いる
n η ( ξ ) = { n B ( ξ ) , if η = + 1 , n F ( ξ ) , if η = − 1. {\displaystyle n_{\eta }(\xi )={\begin{cases}n_{B}(\xi ),&{\text{if }}\eta =+1,\\n_{F}(\xi ),&{\text{if }}\eta =-1.\end{cases}}} 双曲線関数との関係 ボース分布関数は双曲線コタンジェント関数と次の関係にある。
n B ( ξ ) = 1 2 ( coth β ξ 2 − 1 ) . {\displaystyle n_{B}(\xi )={\frac {1}{2}}\left(\operatorname {coth} {\frac {\beta \xi }{2}}-1\right).} フェルミ分布関数は双曲線タンジェント関数と次の関係にある。
n F ( ξ ) = 1 2 ( 1 − tanh β ξ 2 ) . {\displaystyle n_{F}(\xi )={\frac {1}{2}}\left(1-\operatorname {tanh} {\frac {\beta \xi }{2}}\right).} パリティ どちらの分布関数も決まったパリティを持っていない。
n η ( − ξ ) = − η − n η ( ξ ) {\displaystyle n_{\eta }(-\xi )=-\eta -n_{\eta }(\xi )} これは関数 c η {\displaystyle c_{\eta }} を用いて次のようにも書ける。
n η ( − ξ ) = n η ( ξ ) + 2 ξ c η ( 0 , ξ ) {\displaystyle n_{\eta }(-\xi )=n_{\eta }(\xi )+2\xi c_{\eta }(0,\xi )} しかしこれらの導関数は決まったパリティを持つ。
ボーズ・フェルミ変質 ボース分布関数とフェルミ分布関数は、フェルミオン振動数による変数のシフトの下で変質する。
n η ( i ω m + ξ ) = − n − η ( ξ ) {\displaystyle n_{\eta }(i\omega _{m}+\xi )=-n_{-\eta }(\xi )} しかしボソン振動数によるシフトでは違いは生じない。
微分 一次 n B ′ ( ξ ) = − β 4 c s c h 2 β ξ 2 , {\displaystyle n_{B}^{\prime }(\xi )=-{\frac {\beta }{4}}\mathrm {csch} ^{2}{\frac {\beta \xi }{2}},} n F ′ ( ξ ) = − β 4 s e c h 2 β ξ 2 . {\displaystyle n_{F}^{\prime }(\xi )=-{\frac {\beta }{4}}\mathrm {sech} ^{2}{\frac {\beta \xi }{2}}.} 積で表すと、
n η ′ ( ξ ) = − β n η ( ξ ) ( 1 + η n η ( ξ ) ) . {\displaystyle n_{\eta }^{\prime }(\xi )=-\beta n_{\eta }(\xi )(1+\eta n_{\eta }(\xi )).} 温度ゼロの極限では、
n η ′ ( ξ ) = η δ ( ξ ) as β → ∞ . {\displaystyle n_{\eta }^{\prime }(\xi )=\eta \delta (\xi ){\text{ as }}\beta \rightarrow \infty .} 二次 n B ′ ′ ( ξ ) = β 2 4 csch 2 β ξ 2 coth β ξ 2 , {\displaystyle n_{B}^{\prime \prime }(\xi )={\frac {\beta ^{2}}{4}}\operatorname {csch} ^{2}{\frac {\beta \xi }{2}}\operatorname {coth} {\frac {\beta \xi }{2}},} n F ′ ′ ( ξ ) = β 2 4 sech 2 β ξ 2 tanh β ξ 2 . {\displaystyle n_{F}^{\prime \prime }(\xi )={\frac {\beta ^{2}}{4}}\operatorname {sech} ^{2}{\frac {\beta \xi }{2}}\operatorname {tanh} {\frac {\beta \xi }{2}}.} Formula of difference n η ( a + b ) − n η ( a − b ) = − s i n h β b c o s h β a − η c o s h β b . {\displaystyle n_{\eta }(a+b)-n_{\eta }(a-b)=-{\frac {\mathrm {sinh} \beta b}{\mathrm {cosh} \beta a-\eta \,\mathrm {cosh} \beta b}}.} a = 0の場合 n B ( b ) − n B ( − b ) = c o t h β b 2 , {\displaystyle n_{B}(b)-n_{B}(-b)=\mathrm {coth} {\frac {\beta b}{2}},} n F ( b ) − n F ( − b ) = − t a n h β b 2 . {\displaystyle n_{F}(b)-n_{F}(-b)=-\mathrm {tanh} {\frac {\beta b}{2}}.} a → 0の場合 n B ( a + b ) − n B ( a − b ) = coth β b 2 + n B ′ ′ ( b ) a 2 + ⋯ , {\displaystyle n_{B}(a+b)-n_{B}(a-b)=\operatorname {coth} {\frac {\beta b}{2}}+n_{B}^{\prime \prime }(b)a^{2}+\cdots ,} n F ( a + b ) − n F ( a − b ) = − tanh β b 2 + n F ′ ′ ( b ) a 2 + ⋯ . {\displaystyle n_{F}(a+b)-n_{F}(a-b)=-\operatorname {tanh} {\frac {\beta b}{2}}+n_{F}^{\prime \prime }(b)a^{2}+\cdots .} b → 0の場合 n B ( a + b ) − n B ( a − b ) = 2 n B ′ ( a ) b + ⋯ , {\displaystyle n_{B}(a+b)-n_{B}(a-b)=2n_{B}^{\prime }(a)b+\cdots ,} n F ( a + b ) − n F ( a − b ) = 2 n F ′ ( a ) b + ⋯ . {\displaystyle n_{F}(a+b)-n_{F}(a-b)=2n_{F}^{\prime }(a)b+\cdots .} 関数c η 定義:
c η ( a , b ) ≡ − n η ( a + b ) − n η ( a − b ) 2 b . {\displaystyle c_{\eta }(a,b)\equiv -{\frac {n_{\eta }(a+b)-n_{\eta }(a-b)}{2b}}.} ボソンとフェルミオンでは:
c B ( a , b ) ≡ c + ( a , b ) , {\displaystyle c_{B}(a,b)\equiv c_{+}(a,b),} c F ( a , b ) ≡ c − ( a , b ) . {\displaystyle c_{F}(a,b)\equiv c_{-}(a,b).} 双曲線関数との関係 c η ( a , b ) = sinh β b 2 b ( cosh β a − η cosh β b ) . {\displaystyle c_{\eta }(a,b)={\frac {\sinh \beta b}{2b(\cosh \beta a-\eta \cosh \beta b)}}.} c F ( a , b ) {\displaystyle c_{F}(a,b)} は正の定関数であることは明らかである。 数値計算でのオーバーフローを避けるため、tanh関数やcoth関数が用いられる。
c B ( a , b ) = 1 4 b ( coth β ( a − b ) 2 − coth β ( a + b ) 2 ) , {\displaystyle c_{B}(a,b)={\frac {1}{4b}}\left(\operatorname {coth} {\frac {\beta (a-b)}{2}}-\operatorname {coth} {\frac {\beta (a+b)}{2}}\right),} c F ( a , b ) = 1 4 b ( tanh β ( a + b ) 2 − tanh β ( a − b ) 2 ) . {\displaystyle c_{F}(a,b)={\frac {1}{4b}}\left(\operatorname {tanh} {\frac {\beta (a+b)}{2}}-\operatorname {tanh} {\frac {\beta (a-b)}{2}}\right).} a = 0の場合 c B ( 0 , b ) = − 1 2 b coth β b 2 , {\displaystyle c_{B}(0,b)=-{\frac {1}{2b}}\operatorname {coth} {\frac {\beta b}{2}},} c F ( 0 , b ) = 1 2 b tanh β b 2 . {\displaystyle c_{F}(0,b)={\frac {1}{2b}}\operatorname {tanh} {\frac {\beta b}{2}}.} b = 0の場合 c B ( a , 0 ) = β 4 csch 2 β a 2 , {\displaystyle c_{B}(a,0)={\frac {\beta }{4}}\operatorname {csch} ^{2}{\frac {\beta a}{2}},} c F ( a , 0 ) = β 4 sech 2 β a 2 . {\displaystyle c_{F}(a,0)={\frac {\beta }{4}}\operatorname {sech} ^{2}{\frac {\beta a}{2}}.} 低温極限 a = 0において: c F ( 0 , b ) = 1 2 | b | . {\displaystyle c_{F}(0,b)={\frac {1}{2|b|}}.}
b = 0において: c F ( a , 0 ) = δ ( a ) . {\displaystyle c_{F}(a,0)=\delta (a).}
一般的に、
c F ( a , b ) = { 1 2 | b | , if | a | < | b | 0 , if | a | > | b | {\displaystyle c_{F}(a,b)={\begin{cases}{\frac {1}{2|b|}},&{\text{if }}|a|<|b|\\0,&{\text{if }}|a|>|b|\end{cases}}} 関連項目 外部リンク Agustin Nieto: Evaluating Sums over the Matsubara Frequencies . arXiv:hep-ph/9311210 Github repository: MatsubaraSum A Mathematica package for Matsubara frequency summation. 参考文献 ^ A. Abrikosov , L. Gor'kov , I. Dzyaloshinskii : Methods of Quantum Field Theory in Statistical Physics. , New York, Dover Publ., 1975, ISBN 0-486-63228-8