分配関数

統計力学において、**分配関数**（ぶんぱいかんすう、）または**状態和**（じょうたいわ、）はカノニカル分布の理論で導入される量である。カノニカル分布では、温度、体積、粒子数が指定された熱平衡状態において、系があるエネルギー
の状態をとる確率はボルツマン因子
に比例する。このとき、分配関数 は系がとりうる全ての状態についてのボルツマン因子の和 で定義される。また、はボルツマン因子を分配関数で規格化した
で与えられる。分配関数を表す記号 はドイツ語で状態和を表す語Zustandssummeに由来し、マックス・プランクによって導入されたものである。化学分野では記号  の代わりに記号  を用いることがあり、IUPACは両者の表記を採用している。統計力学の形成において、アンサンブル理論を導入したウィラード・ギブズは相積分、マックス・プランクは状態和と呼び、後にラルフ・ファウラーは分配関数と名付けた。熱力学との対応において、ヘルムホルツの自由エネルギー は分配関数との関係で結びつき、熱平衡状態における系の熱力学的量は分配関数から全て求められる。

一方、グランドカノニカル分布において同様の役割を担う関数を**大分配関数**（だいぶんぱいかんすう、）と呼び、\Xi\,あるいは\mathcal{Z}で表す。

## 分配関数
系の取りうる全ての状態の集合を  とし、系が状態  にあるときのエネルギーを \mathcal{E}(\omega) とするとき、分配関数  は

によって定義される。和の中の \exp\{ -\beta\mathcal{E}(\omega) \} はボルツマン因子と呼ばれる。カノニカルアンサンブルは熱浴と接触する閉鎖系を表現するアンサンブルである。パラメータ  は熱浴を特徴づける量で、熱浴の温度と解釈される。熱力学温度  とは  の関係にあり、逆温度と呼ばれる。 はボルツマン定数である。分配関数に定数を乗じることはエネルギーの基準値をずらすことに等しい。分配関数の大きさそのものには意味がない。

熱平衡状態において、系がエネルギーの状態を取る確率は
: P(E_{i})=\frac{g_i e^{- \beta E_i}}{Z} 
で与えられる。ここではエネルギーの状態の縮退度であり、これはを満たす状態 の数である 。系が取りうるエネルギーにわたる分子の和は
: \sum_{i}g_i e^{- \beta E_{i}}=\sum_{\omega \isin \Omega} e^{- \beta \mathcal{E}(\omega)}=Z 
であり、確率の和は分配関数によって1に規格化される。物理量の期待値は
: \langle A \rangle =\sum_{i}A(E_{i})P(E_{i})=\frac{1}{Z}\sum_{i} g_{i} A(E_{i})e^{- \beta E_{i}} 
となる。特にエネルギーについては、
: 
\langle E \rangle =\frac{1}{Z}\sum_{i}g_i E e^{- \beta E_{i}} =-\frac{1}{Z}\frac{\partial Z}{\partial \beta} 
=-\frac{\partial \ln{Z}}{\partial \beta}

と分配関数の対数微分で表される。

### 量子系
量子系においては、系の状態はヒルベルト空間上の状態ベクトル \vert \psi \rangle で表される。ある状態における物理量は量子論的な演算子で与えられ、特にエネルギーはハミルトン演算子 \hat{\mathcal{H}} で与えられる。したがって、分配関数は
 \} \vert \psi \rangle
}}
となる。
状態ベクトルはパラメータ  で指定される正規直交完全系 \vert n \rangle により

と展開される。状態ベクトルに対する和は展開係数に関する積分に置き換えられるので
 \} \vert n \rangle \\
 &=\sum_{m,n} \prod_l \int dc_l d\bar{c}_l c_n \bar{c}_m
 \langle m \vert \exp\{ -\beta\hat{\mathcal{H}} \} \vert n \rangle \\
 &=C \sum_{m,n} \delta_{m,n} \langle m \vert \exp\{ -\beta\hat{\mathcal{H}} \} \vert n \rangle \\
 &=C \sum_n \langle n \vert \exp\{ -\beta\hat{\mathcal{H}} \} \vert n \rangle \\
\end{align} 
}}
となる。分配関数の大きさそのものには意味がないので係数  を除くことができて、最終的には
 \} \vert n \rangle
}}
となる。トレースを用いれば
 \} ]
}}
と表現できる。

量子系では通常はハミルトン演算子を対角化するエネルギー固有状態を用いて表現される。エネルギー量子数  と対応するエネルギー固有値  により

となる。
ここで  は全てのエネルギー固有状態についての和であり、縮退などがある場合には注意を要する。

### 古典系
古典系では、状態変数は連続的に変化するので、状態毎の和をとることが出来ない。そこで、粗視化を行い、位置と運動量が「あまり変わらない」状態を同一の状態と考える。

例えば、1次元空間内の1粒子からなる系では、量子状態が位相空間において「面積」に1つの割合で分布すると考え、ボルツマン因子の位相空間上の積分をで割ったものを分配関数と定義する。

ここで、は位相空間上の点におけるハミルトニアンである。

これは系が次元空間内の個の同一粒子からなる場合にも簡単に拡張できて、

\iint\!\cdots\!\int {\mathrm d}^d p_1 \cdots {\mathrm d}^d p_N \, {\mathrm d}^d q_1 \cdots {\mathrm d}^d q_N \,
e^{-\beta H({\mathbf p}_1, \cdots, {\mathbf p}_N,{\mathbf q}_1, \cdots, {\mathbf q}_N)} }}

ここで、は、粒子が区別出来ないことによる状態の数え過ぎを補正するための項である。

## 独立な系への分解
系のハミルトニアンが

:H=H_a+H_b+H_c+ \cdots 

と独立な項に分けられ、対応するエネルギーが

:E=E_a+E_b+E_c+ \cdots 

という部分和に表される場合、分配関数は

:Z=Z_a \cdot Z_b \cdot Z_c \cdots 

という積の形で表される。
粒子間の相互作用のない粒子数 の系においては、分配関数は1粒子の分配関数によって、

:Z=(Z_1)^N 

との乗の形に表される。

## 大分配関数
系の取りうる全ての状態の集合を  とし、系が状態  にあるときのエネルギーを \mathcal{E}(\omega)、粒子数を \mathcal{N}(\omega) とするとき、大分配関数  は

によって定義される。グランドカノニカルアンサンブルは熱浴、粒子浴と接触する解放系を表現するアンサンブルである。パラメータ  は粒子浴の化学ポテンシャルである。

### 分配関数との関係
集合  を粒子数  によって

の非交和に分解する。これを用いて大分配関数を変形すれば

となる。ここで  は活量である。大分配関数は粒子数  の分配関数の母関数と見ることができる。

## 熱力学との関係
分配関数は統計力学を熱力学に関係付ける上で重要な関数である。
系のヘルムホルツエネルギー  は

で定義される。
温度の関数として表されたヘルムホルツエネルギーは完全な熱力学関数であり、系の熱力学的な性質の全てを導くことが可能である。
この式はカノニカルアンサンブルにおいて、マクロな熱力学関数をミクロな統計力学に基づいて導く式である。

大分配関数を用いて定義される

はグランドポテンシャルと呼ばれる。温度と化学ポテンシャルの関数としてのグランドポテンシャルも完全な熱力学関数であり、グランドカノニカルアンサンブルにおいて、統計力学に基づいて熱力学関数を導く式である。

別の表現として、逆温度  の関数として表された以下の関数も完全な熱力学関数になっている。

 を、  をという。

## 状態密度 、分配関数  、大分配関数  の関係
熱力学関数どうしがルジャンドル変換で関係づけられていることに対応して、状態密度 、分配関数  および大分配関数  の間はラプラス変換を通じて結びついている。
:\begin{align}Z(\beta,V,N)&=\int_0^\infty e^{-\beta E}\Omega(E,V,N)\mathrm{d}E,\\\Xi(\beta,V,\mu)&=\sum _{N=0}^{\infty }e^{\beta\mu N}Z(T,V,N)\end{align}
の関係がある。

また、等温定圧集団については分配関数  から
:\mathcal{Z}(T,P,N):=\int_0^\infty Z(T,V,N)e^{-PV/kT}\mathrm{d}V
で与えられる*T*-*P*分配関数を用いて、
:G(T,P,N)=-kT\ln\mathcal{Z}(T,P,N)
でギブス自由エネルギーを表すことができる。

## 分配関数の例
### 調和振動子
;古典系
個の独立な調和振動子から構成される古典系を考える。調和振動子の質量をとし、角振動数をとする。系のハミルトニアンは

:
H_N=\sum_{i=1}^{N}{h(q_i, p_i)}=\sum_{i=1}^{N} \biggl ( \frac{p_{i}^{\,2}}{2m}
+\frac{1}{2}m \omega^{2} q_{i}^{\, 2} \biggr )


となる。但し、

:
h(q_i, p_i)=\frac{p_{i}^{\,2}}{2m}+\frac{1}{2} m \omega^{2} q_{i}^{\,2}


は一つの調和振動子のハミルトニアンであり、それぞれの調和振動子は区別できるとする。このとき、分配関数は調和振動子間の相互作用が無いことから

: 
\begin{align}
Z(\beta, N) &= \frac{1}{(2 \pi \hbar)^{N}}
\iint\!\cdots\!\int \mathrm{d} p_1 \cdots \mathrm{d} p_N \, \mathrm{d} q_1 \cdots \mathrm{d} q_N \,
e^{-\beta H_N} \\
&=
\biggl ( \frac{1}{(2 \pi \hbar)} \iint \mathrm{d} p_i\mathrm{d} q_i \, e^{-\beta h(q_i,p_i)}
\biggr )^N \\
&=\bigl ( Z(\beta,1) \bigr ) ^N
\end{align}


との積で表される。ここでのにわたる積分がに、にわたる積分がになることから

:
Z(\beta, 1)=\frac{1}{\hbar \beta \omega}

:
Z(\beta, N)= \biggl ( \frac{1}{\hbar \beta \omega} \biggr )^{N}


となる。

;量子系
個の独立な調和振動子から構成される量子系を考える。調和振動子の質量をとし、角振動数をとする。一つの調和振動子のハミルトニアンに対するエネルギー固有値は

:
E_{n_i}=\hbar \omega \biggl (n_{i}+ \frac{1}{2} \biggr ) \quad (n_i=0, 1, 2, \cdots )


であり、系のハミルトニアンのエネルギー固有値は

:
E=\sum_{i=1}^{N}E_{n_i}=\sum_{i=1}^{N}\hbar \omega \biggl (n_{i}+ \frac{1}{2} \biggr )


になる。このとき、分配関数は

:
\begin{align}
Z ( \beta ,N ) &=\sum_{n_1=0}^{\infty}\cdots \sum_{n_N=0}^{\infty} \exp{\biggl (- \beta \sum_{i=1}^{N} E_{n_i} } \biggr ) \\
&=\sum_{n_1=0}^{\infty}\cdots \sum_{n_N=0}^{\infty} \prod_{i=1}^{N} e^{- \beta  E_{n_i}}\\
&=\prod_{i=1}^{N}\sum_{n_i=0}^{\infty} e^{-\beta E_{n_i}} \\
& =(Z(\beta_,1))^{N}
\end{align}


となる。一つの調和振動子におけるは

:
Z(\beta ,1)=\sum_{n=0}^{\infty}e^{-\beta E_n}
=e^{-\frac{1}{2}\beta \hbar \omega}\sum_{n=0}^{\infty}e^{-\beta \hbar \omega n }


となる。これは等比級数であるから、

:
Z(\beta ,1)=\frac{e^{-\frac{1}{2}\beta \hbar \omega}}{1-e^{- \beta \hbar \omega}}
= \biggl [ 2\sinh{\frac{\beta \hbar \omega}{2}} \biggr ]^{-1}


と求まり、

:
Z(\beta ,N)=\bigl ( Z(\beta,1) \bigr ) ^N
= \biggl [ 2\sinh{\frac{\beta \hbar \omega}{2}} \biggr ]^{-N}


になる。

### 単原子分子の理想気体
個の単原子分子の粒子から構成される理想気体の古典系を考える。単原子分子の質量をとする。番目の粒子の正準座標を、正準運動量をとすると系のハミルトニアンは

:
H_N=\sum_{i=1}^{N}{h(q_i, p_i)}=\sum_{i=1}^{N}\frac{1}{2m}\bigl ( p_{ix}^{\, 2}+p_{iy}^{\, 2}+p_{iz}^{\, 2}\bigr )


となる。但し、

:
h(q_i, p_i)=\frac{1}{2m}\bigl ( p_{ix}^{\, 2}+p_{iy}^{\, 2}+p_{iz}^{\, 2}\bigr )


は1粒子のハミルトニアンである。このとき、分配関数は粒子間の相互作用が無いことから

: 
\begin{align}
Z(\beta, N) &= \frac{1}{N!\,(2 \pi \hbar)^{3N}}
\iint\!\cdots\!\int \mathrm{d}^3 p_1 \cdots \mathrm{d}^3 p_N \, \mathrm{d}^3 q_1 \cdots \mathrm{d}^3 q_N \,
e^{-\beta H_N} \\
&=\frac{1}{N!}
\biggl ( \frac{1}{(2 \pi \hbar)^{3}} \iint \mathrm{d}^3 p_i\mathrm{d}^3 q_i \, e^{-\beta h(q_i,p_i)}
\biggr )^N \\
&=\frac{1}{N!}\bigl ( Z(\beta,1) \bigr ) ^N
\end{align}


との積で表される。ここでのにわたる積分が体積、にわたる積分がになることから

:
Z(\beta, 1)=V \biggl ( \frac{m}{2\pi \hbar^2 \beta} \biggr )^{\frac{3}{2}}
=V\biggl ( \frac{2\pi m}{h^2 \beta} \biggr )^{\frac{3}{2}}

:
Z(\beta, N)=\frac{V^N}{N!} \biggl ( \frac{m}{2\pi \hbar^2 \beta} \biggr )^{\frac{3N}{2}}=\frac{V^N}{N!} \biggl ( \frac{2\pi m}{h^2 \beta} \biggr )^{\frac{3N}{2}}


となる。

## 脚注
### 注

### 出典

## 参考文献
*
*
* 
*
*
* 

## 関連記事
* 統計力学
* 統計集団
** カノニカルアンサンブル - グランドカノニカルアンサンブル
* 状態数
* 熱力学

## 外部リンク
*