松原振動数
2017年8月26日 (土) 13:03 (UTC)
熱場の量子論において松原振動数の和とは、離散的な虚数振動数についての和のこと。松原武生に因んで名付けられた。松原振動数の和は次の形をとる。
S_\eta = \frac{1}{\beta}\sum_{i\omega_n} g(i\omega_n)
ここで\beta = \hbar / k_B Tは逆温度で、振動数\omega_nは次の2種類のどちらかである(ただしn\in\mathbb{Z})。
ボソン振動数: \omega_n=\frac{2n\pi}{\beta}
フェルミオン振動数: \omega_n=\frac{(2n+1)\pi}{\beta}
g(z=i\omega)がz\to\inftyの極限でz^{-1}よりも速く0に収束するとき、この和は収束する。
ボソン振動数についての和は S_B (with \eta=+1)と表され、フェルミオン振動数についての和はS_F (with \eta=-1)と表される。ここで\etaは統計的な記号である。
熱場の量子論に加えて、松原振動数の和は固体物理学における有限温度でのファインマン・ダイアグラムを考える上で重要な役割を果たす。
一般的にファインマン・ダイアグラムは、T=0\,\text{K}では積分\int_{T=0}  \mathrm{ d}\omega \ g(\omega )で表されるが、有限温度では和S_\etaで与えられる。
松原振動数の和を評価する上手なやり方は、 z=i\omegaに極を持つ松原重み関数hη(z)を使う方法である。
ボソンの場合η = +1とフェルミオンの場合η = −1で重み関数は異なる。
重み関数の選択について後述する。
和は、重み関数を使って複素平面での閉曲線積分に置き換えることができる。
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_\eta=-\frac 1 \beta \sum_{z_0\in g(z)\text{ poles} \operatorname{Res} g(z_0) h_\eta(z_0)
ここで閉曲線が極を時計回りの方向で囲むように変形し、負の留数を生むため、マイナスがつくことに注意。
ボソン振動数z=i\omega_nの極を作るために、どちらの半平面で収束がコントロールされるかに依存して次の2つのタイプの松原重み関数を選ぶことができる。
h_B^{(1)}(z)=\frac{\beta}{1-e^{-\beta z}=-\beta n_B(-z)=\beta(1+n_B(z))
h_B^{(2)}(z)=\frac{-\beta}{1-e^{\beta z}=\beta n_B(z)
h_B^{(1)}(z)は左半平面(Re z < 0)での収束をコントロールし、h_B^{(2)}(z)は右半平面(Re z > 0)での収束をコントロールする。
ここでn_B(z)=(e^{\beta z}-1)^{-1}はボース分布関数である。
フェルミオン振動数の場合も同様である。
2つのタイプの松原重み関数があり、z=i\omega_mに極を作る。
h_F^{(1)}(z)=\frac{\beta}{1+e^{-\beta z}=\beta n_F(-z)=\beta(1-n_F(z))
h_F^{(2)}(z)=\frac{-\beta}{1+e^{\beta z}=-\beta n_F(z)
h_F^{(1)}(z)は左半平面(Re z < 0)での収束をコントロールし、h_F^{(1)}(z)は右半平面(Re z > 0)での収束をコントロールする。
ここでn_F(z)=(e^{\beta z}+1)^{-1}はフェルミ分布関数である。
グリーン関数への応用では、g(z)は常に次の構造を持つ。
g(z)=G(z)e^{-z\tau}
これは0 < τ < βで与えられる左半平面で発散する。
収束をコントロールするために、第一のタイプの重み関数は常にh_\eta(z)=h_\eta^{(1)}(z)と選ぶ。
しかし松原振動数の和が発散しないときは収束をコントロールする必要はない。
そのような場合、どんな松原重み関数を選んでも同じ結果が得られる。
以下の表に、
いくつかの簡単な有理関数g(z)での松原振動数の和をまとめる。
S_\eta=\frac{1}{\beta}\sum_{i\omega}g(i\omega).
η = ±1 は統計的記号である。
[1] 和は収束しないため、松原重み関数の選択が異なると結果は定数分だけ異なる。
[2] (1 ↔ 2)は、手前の項のインデックス1と2を置き換えたものを表す。
極限\beta\rightarrow\inftyでの松原振動数の和は、虚数振動数の虚軸についての積分に等しい
\frac{1}{\beta}\sum_{i\omega}=\int_{-i\infty}^{i\infty}\frac{\mathrm{d}(i\omega)}{2\pi}
いくつかの積分は収束しない。
それらは振動数カットオフ\Omegaを導入し、\Omega\rightarrow\inftyの極限をとる前に積分から発散部分(\Omega-依存)を差し引いて繰り込む必要がある。
例えば自由エネルギーは、対数の積分によって得られる。
\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\geq0, \\
-\eta\xi & \xi<0,
\end{array}
\right.
これは、温度ゼロにおいて自由エネルギーは化学ポテンシャルに満たない内部エネルギーと単純な関係にあることを意味している。
また分布関数は次の積分によって得られる。
\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\geq0, \\
-\eta & \xi<0,
\end{array}
\right.
これは温度ゼロでの階段関数のふるまいを示している。
虚時間区間(0,β)で定義される関数G(τ)を考える。
これはフーリエ級数の観点で与えられる。
G(\tau)=\frac{1}{\beta}\sum_{i\omega} G(i\omega) e^{-i\omega\tau},
ここで振動数は 2pi/β間隔の離散的な値のみとる。
振動数の選択は、関数G(τ)の境界条件に依存している。
物理学ではG(τ)はグリーン関数の虚時間表現を表す。
G(\tau)=-\langle \mathcal{T}_\tau \psi(\tau)\psi^*(0) \rangle.
これはボソン場の周期的境界条件G(τ+β)=G(τ)を満たす。
一方フェルミオン場では、境界条件は反周期的G(τ + β) = −G(τ)である。
振動数領域でのグリーン関数G(iω)が与えられたとき、その虚時間表現G(τ)は松原振動数の和によって評価できる。
その和がボソン振動数かフェルミオン振動数のどちらでとるかに依存して、得られるG(τ)は異なる。
これらを区別するため、次を定義する。
G_\eta(\tau)= \begin{cases}
G_B(\tau), & \text{if } \eta = +1, \\
G_F(\tau), & \text{if } \eta = -1,
\end{cases}
G_B(\tau)=\frac{1}{\beta}\sum_{i\omega_n}G(i\omega_n)e^{-i\omega_n\tau},
G_F(\tau)=\frac{1}{\beta}\sum_{i\omega_m}G(i\omega_m)e^{-i\omega_m\tau}.
ここでτは区間(0,β)に制限されていることに注意。
境界条件は区間の外にG(τ)を拡張するために用いることができる。
よく用いられる結果を以下の表にまとめる。
ここでは小さな虚時間が決定的な役割を果たす。
小さな虚時間の負号が変わると、演算子の順番が変わる
\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^+}.\
\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) = \sum_{i\omega}(i\omega-\xi)^{-1}
これは、G(τ)をτ = 0からわずかに遠ざけたとき、収束をコントロールするため、G(\tau=0^+)での重み関数としてh_\eta^{(1)}(z)をとらなければならず、G(\tau=0^-)ではh_\eta^{(2)}(z)をとらなければならないと理解できる。
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(\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(\tau=0^-)=\frac{1}{\beta}\sum_{i\omega_m}\frac{e^{i\omega_m 0^+}{i\omega_m-\xi}=n_F(\xi),
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)).
\frac{1}{\beta}\sum_{i\omega_n} \ln(\beta(-i\omega_n+\xi))=\frac{1}{\beta}\ln(1-e^{-\beta\xi}),
-\frac{1}{\beta}\sum_{i\omega_m} \ln(\beta(-i\omega_m+\xi))=-\frac{1}{\beta}\ln(1+e^{-\beta\xi}).
ここでは、よく使われる単一モードのダイアグラムを評価する。
多重モード問題は、スペクトル関数積分によってアプローチできる。
\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 }.
\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'}.
\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_\etaはボース分布関数(η = +1)かフェルミ分布関数(η = −1)のどちらかを表す。
n_\eta(\xi)=\frac{1}{e^{\beta\xi}-\eta}
必要であればボース分布関数とフェルミ分布関数を区別するために、それぞれ記号nBとnFを用いる
n_\eta(\xi)= \begin{cases}
n_B(\xi), & \text{if } \eta = +1, \\
n_F(\xi), & \text{if } \eta = -1.
\end{cases}
ボース分布関数は双曲線コタンジェント関数と次の関係にある。
n_B(\xi)=\frac{1}{2}\left(\operatorname{coth}\frac{\beta\xi}{2}-1\right).
フェルミ分布関数は双曲線タンジェント関数と次の関係にある。
n_F(\xi)=\frac{1}{2}\left(1-\operatorname{tanh}\frac{\beta\xi}{2}\right).
どちらの分布関数も決まったパリティを持っていない。
n_\eta(-\xi)=-\eta-n_\eta(\xi)
これは関数c_\etaを用いて次のようにも書ける。
n_\eta(-\xi)=n_\eta(\xi)+2\xi c_\eta(0,\xi)
しかしこれらの導関数は決まったパリティを持つ。
ボース分布関数とフェルミ分布関数は、フェルミオン振動数による変数のシフトの下で変質する。
n_\eta(i\omega_m+\xi)=-n_{-\eta}(\xi)
しかしボソン振動数によるシフトでは違いは生じない。
n_B^\prime(\xi)=-\frac{\beta}{4}\mathrm{csch}^2\frac{\beta \xi}{2},
n_F^\prime(\xi)=-\frac{\beta}{4}\mathrm{sech}^2\frac{\beta \xi}{2}.
積で表すと、
n_\eta^\prime(\xi)= -\beta n_\eta(\xi)(1+\eta n_\eta(\xi)).
温度ゼロの極限では、
n_\eta^\prime(\xi)=\eta\delta(\xi) \text{ as } \beta\rightarrow\infty.
n_B^{\prime\prime}(\xi)=\frac{\beta^2}{4}\operatorname{csch}^2\frac{\beta \xi}{2}\operatorname{coth}\frac{\beta \xi}{2},
n_F^{\prime\prime}(\xi)=\frac{\beta^2}{4}\operatorname{sech}^2\frac{\beta \xi}{2}\operatorname{tanh}\frac{\beta \xi}{2}.
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)=\mathrm{coth}\frac{\beta b}{2},
n_F(b)-n_F(-b)=-\mathrm{tanh}\frac{\beta b}{2}.
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)=-\operatorname{tanh}\frac{\beta b}{2}+n_F^{\prime\prime}(b)a^2+\cdots.
n_B(a+b)-n_B(a-b)=2n_B^\prime(a)b+\cdots,
n_F(a+b)-n_F(a-b)=2n_F^\prime(a)b+\cdots.
定義:
c_\eta(a,b)\equiv-\frac{n_\eta(a+b)-n_\eta(a-b)}{2b}.
ボソンとフェルミオンでは:
c_B(a,b)\equiv c_+(a,b),
c_F(a,b)\equiv c_-(a,b).
c_\eta(a,b)=\frac{\sinh\beta b}{2b(\cosh\beta a-\eta\cosh\beta b)}.
c_F(a,b)は正の定関数であることは明らかである。
数値計算でのオーバーフローを避けるため、tanh関数やcoth関数が用いられる。
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)=\frac{1}{4b}\left(\operatorname{tanh}\frac{\beta(a+b)}{2} - \operatorname{tanh}\frac{\beta(a-b)}{2}\right).
c_B(0,b)=-\frac{1}{2b}\operatorname{coth}\frac{\beta b}{2},
c_F(0,b)=\frac{1}{2b}\operatorname{tanh}\frac{\beta b}{2}.
c_B(a,0)=\frac{\beta}{4}\operatorname{csch}^2\frac{\beta a}{2},
c_F(a,0)=\frac{\beta}{4}\operatorname{sech}^2\frac{\beta a}{2}.
a = 0において: c_F(0,b)=\frac{1}{2|b〕.
b = 0において: c_F(a,0)=\delta(a).
一般的に、
c_F(a,b)=\begin{cases}
\frac{1}{2|b〕, & \text{if } |a|<|b| \\
0,  & \text{if } |a|>|b|
\end{cases}