統計力学 において、分配関数 (ぶんぱいかんすう、英 : Partition function )または状態和 (じょうたいわ、英 : state sum, sum over states )は、ある系 の物理量 の統計集団 的平均を計算する際に用いられる規格化定数を指す。単に分配関数と呼ぶときはカノニカル分布 における分配関数を指し、ドイツ語で状態和を表す語Zustandssummeに由来する記号Z で表す[ 1] 。一方、グランドカノニカル分布 において同様の役割を担う関数を大分配関数 (だいぶんぱいかんすう、英 : Grand partition function )と呼び、
Ξ
{\displaystyle \Xi \,}
あるいは
Z
{\displaystyle {\mathcal {Z}}}
で表す。
系の取りうる全ての状態の集合を Ω とし、系が状態 ω ∈Ω にあるときのエネルギーを
E
(
ω
)
{\displaystyle {\mathcal {E}}(\omega )}
とするとき、分配関数 Z (β ) は
Z
(
β
)
=
∑
ω
∈
Ω
exp
{
−
β
E
(
ω
)
}
{\displaystyle Z(\beta )=\sum _{\omega \in \Omega }\exp\{-\beta {\mathcal {E}}(\omega )\}}
によって定義される。和の中の
exp
{
−
β
E
(
ω
)
}
{\displaystyle \exp\{-\beta {\mathcal {E}}(\omega )\}}
はボルツマン因子 と呼ばれる。カノニカルアンサンブル は熱浴と接触する閉鎖系を表現するアンサンブルである。パラメータ β は熱浴を特徴づける量で、熱浴の温度と解釈される。熱力学温度 T とは β =1/kT の関係にあり、逆温度 と呼ばれる。k はボルツマン定数 である。分配関数に定数を乗じることはエネルギーの基準値をずらすことに等しい。分配関数の大きさそのものには意味がない。
系がエネルギーE の状態を取る確率は
p
(
E
)
=
g
E
e
−
β
E
Z
{\displaystyle p(E)={\frac {g_{E}e^{-\beta E}}{Z}}}
で与えられる。ここでgE はエネルギーE の状態の縮退度 である。系が取りうるエネルギーE にわたる分子の和は
∑
E
g
E
e
−
β
E
=
∑
ω
∈
Ω
e
−
β
E
(
ω
)
=
Z
{\displaystyle \sum _{E}g_{E}e^{-\beta E}=\sum _{\omega \in \Omega }e^{-\beta E(\omega )}=Z}
であり、確率p (E ) の和は分配関数によって1に規格化される。物理量A の期待値は
⟨
A
⟩
=
∑
E
A
(
E
)
p
(
E
)
=
1
Z
∑
E
g
E
A
(
E
)
e
−
β
E
{\displaystyle \langle A\rangle =\sum _{E}A(E)p(E)={\frac {1}{Z}}\sum _{E}g_{E}A(E)e^{-\beta E}}
となる。特にE については、
⟨
E
⟩
=
1
Z
∑
E
g
E
E
e
−
β
E
=
−
1
Z
∂
Z
∂
β
=
−
∂
ln
Z
∂
β
{\displaystyle \langle E\rangle ={\frac {1}{Z}}\sum _{E}g_{E}Ee^{-\beta E}=-{\frac {1}{Z}}{\frac {\partial Z}{\partial \beta }}=-{\frac {\partial \ln {Z}}{\partial \beta }}}
と分配関数の対数微分 で表される。
量子系 においては、系の状態はヒルベルト空間 上の状態ベクトル
|
ψ
⟩
{\displaystyle \vert \psi \rangle }
で表される。ある状態における物理量は量子論的 な演算子 で与えられ、特にエネルギーはハミルトン演算子
H
^
{\displaystyle {\hat {\mathcal {H}}}}
で与えられる。したがって、分配関数は
Z
(
β
)
=
∑
ψ
⟨
ψ
|
exp
{
−
β
H
^
}
|
ψ
⟩
{\displaystyle Z(\beta )=\sum _{\psi }\langle \psi \vert \exp\{-\beta {\hat {\mathcal {H}}}\}\vert \psi \rangle }
となる。
状態ベクトルはパラメータ n で指定される正規直交完全系
|
n
⟩
{\displaystyle \vert n\rangle }
により
|
ψ
⟩
=
∑
n
c
n
|
n
⟩
,
⟨
ψ
|
=
∑
n
c
¯
n
⟨
n
|
{\displaystyle \vert \psi \rangle =\sum _{n}c_{n}\vert n\rangle ,~\langle \psi \vert =\sum _{n}{\bar {c}}_{n}\langle n\vert }
と展開される。状態ベクトルに対する和は展開係数に関する積分に置き換えられるので
Z
(
β
)
=
∏
l
∫
d
c
l
d
c
¯
l
∑
m
,
n
c
n
c
¯
m
⟨
m
|
exp
{
−
β
H
^
}
|
n
⟩
=
∑
m
,
n
∏
l
∫
d
c
l
d
c
¯
l
c
n
c
¯
m
⟨
m
|
exp
{
−
β
H
^
}
|
n
⟩
=
C
∑
m
,
n
δ
m
,
n
⟨
m
|
exp
{
−
β
H
^
}
|
n
⟩
=
C
∑
n
⟨
n
|
exp
{
−
β
H
^
}
|
n
⟩
{\displaystyle {\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
(
β
)
=
∑
n
⟨
n
|
exp
{
−
β
H
^
}
|
n
⟩
{\displaystyle Z(\beta )=\sum _{n}\langle n\vert \exp\{-\beta {\hat {\mathcal {H}}}\}\vert n\rangle }
となる。トレース を用いれば
Z
(
β
)
=
t
r
[
exp
{
−
β
H
^
}
]
{\displaystyle Z(\beta )=\mathrm {tr} [\exp\{-\beta {\hat {\mathcal {H}}}\}]}
と表現できる。
量子系では通常はハミルトン演算子を対角化 するエネルギー固有状態 を用いて表現される。エネルギー量子数 i と対応するエネルギー固有値 Ei により
Z
(
β
)
=
∑
i
e
−
β
E
i
{\displaystyle Z(\beta )=\sum _{i}\mathrm {e} ^{-\beta E_{i}}}
となる。
ここで ∑i は全てのエネルギー固有状態についての和であり、縮退 などがある場合には注意を要する。
古典系では、状態変数は連続的に変化するので、状態毎の和をとることが出来ない。そこで、粗視化 を行い、位置と運動量が「あまり変わらない」状態を同一の状態と考える。
例えば、1次元空間内の1粒子からなる系では、量子状態が位相空間 において「面積」
2
π
ℏ
{\displaystyle 2\pi \hbar }
に1つの割合で分布すると考え、ボルツマン因子e -βE の位相空間上の積分を
2
π
ℏ
{\displaystyle 2\pi \hbar }
で割ったものを分配関数と定義する。
Z
(
β
)
=
1
2
π
ℏ
∬
d
p
d
q
e
−
β
H
(
p
,
q
)
{\displaystyle 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
(
β
,
N
)
=
1
N
!
(
2
π
ℏ
)
N
d
∬
⋯
∫
d
d
p
1
⋯
d
d
p
N
d
d
q
1
⋯
d
d
q
N
e
−
β
H
(
p
1
,
⋯
,
p
N
,
q
1
,
⋯
,
q
N
)
{\displaystyle 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!は、粒子が区別出来ないことによる状態の数え過ぎを補正するための項である。
系の取りうる全ての状態の集合を Ω とし、系が状態 ω ∈Ω にあるときのエネルギーを
E
(
ω
)
{\displaystyle {\mathcal {E}}(\omega )}
、粒子数を
N
(
ω
)
{\displaystyle {\mathcal {N}}(\omega )}
とするとき、大分配関数 Ξ (β ,μ ) は
Ξ
(
β
,
μ
)
=
∑
ω
∈
Ω
exp
{
−
β
E
(
ω
)
+
β
μ
N
(
ω
)
}
{\displaystyle \Xi (\beta ,\mu )=\sum _{\omega \in \Omega }\exp\{-\beta {\mathcal {E}}(\omega )+\beta \mu {\mathcal {N}}(\omega )\}}
によって定義される。グランドカノニカルアンサンブル は熱浴、粒子浴と接触する解放系を表現するアンサンブルである。パラメータ μ は粒子浴の化学ポテンシャル である。
集合 Ω を粒子数 N によって
Ω
(
N
)
=
{
ω
∈
Ω
;
N
(
ω
)
=
N
}
{\displaystyle \Omega (N)=\{\omega \in \Omega ;{\mathcal {N}}(\omega )=N\}}
Ω
=
∐
N
Ω
(
N
)
{\displaystyle \Omega =\coprod _{N}\Omega (N)}
の非交和 に分解する。これを用いて大分配関数を変形すれば
Ξ
(
β
,
μ
)
=
∑
N
∑
ω
∈
Ω
(
N
)
exp
{
−
β
E
(
ω
)
+
β
μ
N
(
ω
)
}
=
∑
N
e
β
μ
N
∑
ω
∈
Ω
(
N
)
exp
{
−
β
E
(
ω
)
}
=
∑
N
e
β
μ
N
Z
(
β
,
N
)
=
∑
N
λ
N
Z
(
β
,
N
)
{\displaystyle {\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}}}
となる。ここで λ =eβμ は活量 である。大分配関数は粒子数 N の分配関数の母関数 と見ることができる。
分配関数は統計力学を熱力学 に関係付ける上で重要な関数である。
系のヘルムホルツエネルギー F (β ) は
F
(
β
)
=
−
1
β
ln
Z
(
β
)
{\displaystyle F(\beta )=-{\frac {1}{\beta }}\ln Z(\beta )}
で定義される。
温度の関数として表されたヘルムホルツエネルギーは完全な熱力学関数 であり、系の熱力学的な性質の全てを導くことが可能である。
この式はカノニカルアンサンブルにおいて、マクロな熱力学関数をミクロな統計力学に基づいて導く式である。
大分配関数を用いて定義される
J
(
β
,
μ
)
=
−
1
β
ln
Ξ
(
β
,
μ
)
{\displaystyle J(\beta ,\mu )=-{\frac {1}{\beta }}\ln \Xi (\beta ,\mu )}
はグランドポテンシャル と呼ばれる。温度と化学ポテンシャルの関数としてのグランドポテンシャルも完全な熱力学関数であり、グランドカノニカルアンサンブルにおいて、統計力学に基づいて熱力学関数を導く式である。
別の表現として、逆温度 β の関数として表された以下の関数も完全な熱力学関数になっている。
Ψ
(
β
)
=
−
β
F
(
β
)
=
ln
Z
(
β
)
{\displaystyle \Psi (\beta )=-\beta F(\beta )=\ln Z(\beta )}
q
(
β
,
α
)
=
−
β
J
(
β
,
α
/
β
)
=
ln
Ξ
(
β
,
α
/
β
)
{\displaystyle q(\beta ,\alpha )=-\beta J(\beta ,\alpha /\beta )=\ln \Xi (\beta ,\alpha /\beta )}
Ψ をマシュー関数 (英語版 ) 、 q をクラマース関数 (英語版 ) という。
熱力学関数どうしがルジャンドル変換 で関係づけられていることに対応して、分配関数はラプラス変換 を通じて結びついている[ 2] 。状態密度 Ω(E , V , N ) 、分配関数 Z (β , V , N ) および大分配関数 Ξ(β , V , μ ) の間には
Z
(
β
,
V
,
N
)
=
∫
0
∞
e
−
β
E
Ω
(
E
,
V
,
N
)
d
E
,
Ξ
(
β
,
V
,
μ
)
=
∑
N
=
0
∞
e
β
μ
N
Z
(
T
,
V
,
N
)
{\displaystyle {\begin{aligned}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{aligned}}}
の関係がある。
また、等温定圧集団 については分配関数 Z (β , V , N ) から
Z
(
T
,
P
,
N
)
:=
∫
0
∞
Z
(
T
,
V
,
N
)
e
−
P
V
/
k
T
d
V
{\displaystyle {\mathcal {Z}}(T,P,N):=\int _{0}^{\infty }Z(T,V,N)e^{-PV/kT}\mathrm {d} V}
で与えられるT -P 分配関数を用いて、
G
(
T
,
P
,
N
)
=
−
k
T
ln
Z
(
T
,
P
,
N
)
{\displaystyle G(T,P,N)=-kT\ln {\mathcal {Z}}(T,P,N)}
でギブス自由エネルギー を表すことができる。