松原振動数

熱場の量子論において松原振動数の和とは、離散的な虚数振動数についての和のこと。松原武生に因んで名付けられた。松原振動数の和は次の形をとる。 S η = 1 β ∑ i ω n g ( i ω n ) {\displaystyle S_{\eta }={\frac {1}{\beta }}\sum _{i\omega _{n}}g(i\omega _{n})} ここで β = ℏ / k B T {\displaystyle \beta =\hbar /k_{B}T} は逆温度で、振動数 ω n {\displaystyle \omega _{n}} は次の2種類のどちらかである(ただし n ∈ Z {\displaystyle n\in \mathbb {Z} } )。 ボソン振動数: ω n = 2 n π β {\displaystyle \omega _{n}={\frac {2n\pi }{\beta }}} フェルミオン振動数: ω n = ( 2 n + 1 ) π β {\displaystyle \omega _{n}={\frac {(2n+1)\pi }{\beta }}} g ( z = i ω ) {\displaystyle g(z=i\omega )} が z → ∞ {\displaystyle z\to \infty } の極限で z − 1 {\displaystyle z^{-1}} よりも速く0に収束するとき、この和は収束する。 ボソン振動数についての和は S B {\displaystyle S_{B}} (with η = + 1 {\displaystyle \eta =+1} )と表され、フェルミオン振動数についての和は S F {\displaystyle S_{F}} (with η = − 1 {\displaystyle \eta =-1} )と表される。ここで η {\displaystyle \eta } は統計的な記号である。 熱場の量子論に加えて、松原振動数の和は固体物理学における有限温度でのファインマン・ダイアグラムを考える上で重要な役割を果たす。 一般的にファインマン・ダイアグラムは、 T = 0 K {\displaystyle T=0\,{\text{K}}} では積分 ∫ T = 0 d ω   g ( ω ) {\displaystyle \int _{T=0}\mathrm {d} \omega \ g(\omega )} で表されるが、有限温度では和 S η {\displaystyle S_{\eta }} で与えられる。 松原振動数の和を評価する上手なやり方は、 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}} はボース分布関数である。 フェルミオン振動数の場合も同様である。 ２つのタイプの松原重み関数があり、 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 は統計的記号である。 [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(τ)を拡張するために用いることができる。 よく用いられる結果を以下の表にまとめる。 ここでは小さな虚時間が決定的な役割を果たす。 小さな虚時間の負号が変わると、演算子の順番が変わる ⟨ ψ ψ ∗ ⟩ = ⟨ 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 }}} 必要であればボース分布関数とフェルミ分布関数を区別するために、それぞれ記号nBとnFを用いる 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}}.} 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}}.} 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}}.} 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 .} 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 η ( 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).} 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}}.} 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.