松原振動数

熱場の量子論において**松原振動数**の和とは、離散的な虚数振動数についての和のこと。松原武生に因んで名付けられた。松原振動数の和は次の形をとる。
: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*)を使う方法である。
ボソンの場合*η*&nbsp;=&nbsp;+1とフェルミオンの場合*η*&nbsp;=&nbsp;−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&nbsp;*z*&nbsp;h_B^{(2)}(z)は右半平面(Re&nbsp;*z*&nbsp;>&nbsp;0)での収束をコントロールする。
ここでn_B(z)=(e^{\beta z}-1)^{-1}はボース分布関数である。

フェルミオン振動数の場合も同様である。
２つのタイプの松原重み関数があり、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&nbsp;*z*&nbsp;h_F^{(1)}(z)は右半平面(Re&nbsp;*z*&nbsp;>&nbsp;0)での収束をコントロールする。
ここでn_F(z)=(e^{\beta z}+1)^{-1}はフェルミ分布関数である。

グリーン関数への応用では、*g*(*z*)は常に次の構造を持つ。
:g(z)=G(z)e^{-z\tau}
これは0&nbsp;h_\eta(z)=h_\eta^{(1)}(z)と選ぶ。
しかし松原振動数の和が発散しないときは収束をコントロールする必要はない。
そのような場合、どんな松原重み関数を選んでも同じ結果が得られる。

### 松原振動数の和の表
以下の表に、
いくつかの簡単な有理関数*g*(*z*)での松原振動数の和をまとめる。
:S_\eta=\frac{1}{\beta}\sum_{i\omega}g(i\omega).
*η*&nbsp;=&nbsp;±1 は統計的記号である。

[1] 和は収束しないため、松原重み関数の選択が異なると結果は定数分だけ異なる。

[2] (1&nbsp;↔&nbsp;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
これは、温度ゼロにおいて自由エネルギーは化学ポテンシャルに満たない内部エネルギーと単純な関係にあることを意味している。
また分布関数は次の積分によって得られる。
:\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,*β*)で定義される関数*G*(*τ*)を考える。
これはフーリエ級数の観点で与えられる。
: G(\tau)=\frac{1}{\beta}\sum_{i\omega} G(i\omega) e^{-i\omega\tau},
ここで振動数は&nbsp;2/*β*間隔の離散的な値のみとる。

振動数の選択は、関数*G*(*τ*)の境界条件に依存している。
物理学では*G*(*τ*)はグリーン関数の虚時間表現を表す。
: G(\tau)=-\langle \mathcal{T}_\tau \psi(\tau)\psi^*(0) \rangle. 
これはボソン場の周期的境界条件*G*(*τ*+*β*)=*G*(*τ*)を満たす。
一方フェルミオン場では、境界条件は反周期的*G*(*τ*&nbsp;+&nbsp;*β*)&nbsp;=&nbsp;−*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^+}

#### 分布関数
*τ*&nbsp;=&nbsp;0でのグリーン関数*G*(*τ*)の不連続性のため、分布関数の評価は難しくなる。
次の和を評価するためにボソンとフェルミオンどちらの重み関数を選択することもできるが、結果は異なる。
: G(0) = \sum_{i\omega}(i\omega-\xi)^{-1}
これは、*G*(*τ*)を*τ*&nbsp;=&nbsp;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はボース分布関数(*η*&nbsp;=&nbsp;+1)かフェルミ分布関数(*η*&nbsp;=&nbsp;−1)のどちらかを表す。
:n_\eta(\xi)=\frac{1}{e^{\beta\xi}-\eta}
必要であればボース分布関数とフェルミ分布関数を区別するために、それぞれ記号*n**B*と*n**F*を用いる

: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}.

### Formula of difference
: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)=\mathrm{coth}\frac{\beta b}{2},
:n_F(b)-n_F(-b)=-\mathrm{tanh}\frac{\beta b}{2}.

#### ''a'' → 0の場合
: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.

#### ''b'' → 0の場合
: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''<sub>''η''</sub>
定義:
: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).

#### ''a'' = 0の場合
: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}.

#### ''b'' = 0の場合
: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|
\end{cases}

## 関連項目
* 虚時間
* 熱場の量子論

## 外部リンク
:Agustin Nieto: ''Evaluating Sums over the Matsubara Frequencies''. arXiv:hep-ph/9311210
:Github repository: MatsubaraSum A Mathematica package for Matsubara frequency summation.

## 参考文献