分配関数

「分配係数」とは異なります。

この項目では、統計力学について説明しています。数学における一般化された分配函数については「分配函数 (数学)」を、場の理論の分配函数については「分配函数 (場の量子論)」をご覧ください。

統計力学において、分配関数（ぶんぱいかんすう、英: Partition function）または状態和（じょうたいわ、英: state sum, sum over states）は、ある系の物理量の統計集団的平均を計算する際に用いられる規格化定数を指す。単に分配関数と呼ぶときはカノニカル分布における分配関数を指し、ドイツ語で状態和を表す語Zustandssummeに由来する記号Zで表す^[1]。一方、グランドカノニカル分布において同様の役割を担う関数を大分配関数（だいぶんぱいかんすう、英: Grand partition function）と呼び、 $\Xi \,$ あるいは ${\mathcal {Z}}$ で表す。

分配関数

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

$Z(\beta )=\sum _{\omega \in \Omega }\exp\{-\beta {\mathcal {E}}(\omega )\}$

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

量子系

量子系においては、系の状態はヒルベルト空間上の状態ベクトル $\vert \psi \rangle$ で表される。ある状態における物理量は量子論的な演算子で与えられ、特にエネルギーはハミルトン演算子 ${\hat {\mathcal {H}}}$ で与えられる。したがって、分配関数は

$Z(\beta )=\sum _{\psi }\langle \psi \vert \exp\{-\beta {\hat {\mathcal {H}}}\}\vert \psi \rangle$

となる。状態ベクトルはパラメータ n で指定される正規直交完全系 $\vert n\rangle$ により

$\vert \psi \rangle =\sum _{n}c_{n}\vert n\rangle ,~\langle \psi \vert =\sum _{n}{\bar {c}}_{n}\langle n\vert$

と展開される。状態ベクトルに対する和は展開係数に関する積分に置き換えられるので

${\begin{aligned}Z(\beta )&=\prod _{l}\int dc_{l}d{\bar {c}}_{l}\sum _{m,n}c_{n}{\bar {c}}_{m}\langle m\vert \exp\{-\beta {\hat {\mathcal {H}}}\}\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{aligned}}$

となる。分配関数の大きさそのものには意味がないので係数 C を除くことができて、最終的には

$Z(\beta )=\sum _{n}\langle n\vert \exp\{-\beta {\hat {\mathcal {H}}}\}\vert n\rangle$

となる。トレースを用いれば

$Z(\beta )=\mathrm {tr} [\exp\{-\beta {\hat {\mathcal {H}}}\}]$

と表現できる。

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

$Z(\beta )=\sum _{i}\mathrm {e} ^{-\beta E_{i}}$

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

古典系

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

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

$Z(\beta )={\frac {1}{2\pi \hbar }}\iint {\mathrm {d} }p\,{\mathrm {d} }q\,e^{-\beta H(p,q)}$

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

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

$Z(\beta ,N)={\frac {1}{N!\,(2\pi \hbar )^{Nd}}}\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})}$

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

大分配関数

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

$\Xi (\beta ,\mu )=\sum _{\omega \in \Omega }\exp\{-\beta {\mathcal {E}}(\omega )+\beta \mu {\mathcal {N}}(\omega )\}$

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

分配関数との関係

集合 Ω を粒子数 N によって

$\Omega (N)=\{\omega \in \Omega ;{\mathcal {N}}(\omega )=N\}$

$\Omega =\coprod _{N}\Omega (N)$

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

${\begin{aligned}\Xi (\beta ,\mu )&=\sum _{N}\sum _{\omega \in \Omega (N)}\exp\{-\beta {\mathcal {E}}(\omega )+\beta \mu {\mathcal {N}}(\omega )\}\\&=\sum _{N}\mathrm {e} ^{\beta \mu N}\sum _{\omega \in \Omega (N)}\exp\{-\beta {\mathcal {E}}(\omega )\}\\&=\sum _{N}\mathrm {e} ^{\beta \mu N}Z(\beta ,N)\\&=\sum _{N}\lambda ^{N}Z(\beta ,N)\\\end{aligned}}$