カノニカル分布

提供: miniwiki
移動先:案内検索


正準分布(せいじゅんぶんぷ)、カノニカル分布: canonical distribution)は、統計力学において微視的状態を表現する統計集団の一つである正準集団カノニカルアンサンブル: canonical ensemble)が従う統計分布である。正準集団とは、熱浴English版italiano版との間でのみエネルギーを自由にやりとりできる系を表現する統計集団である。正準分布は、小正準分布大正準分布とは体積が十分に大きい極限(すなわちエネルギーや粒子の出入りが無視できる極限)において熱力学的に等価である。

確率分布

熱浴に接している系が微視的状態 ω をとる確率分布 p(ω) は次式で定義される。

[math]\begin{align} p(\omega) &= \frac{1}{Z(\beta)} \mathrm{e}^{-\beta E(\omega)} \\ Z(\beta) &= \sum_\omega \mathrm{e}^{-\beta E(\omega)} \end{align}[/math]

ここで、E(ω) は系が微視的状態 ω をとるときのエネルギーβ逆温度と呼ばれ、 β = 1/kT(T絶対温度 )で定義される。β は熱浴を特徴付けるパラメータである。kボルツマン定数である。確率分布 p(ω) の分子の exp[−βE(ω)]ボルツマン因子と呼ばれる。 p(ω) の分母に現れた規格化定数 Z(β) は正準分布の分配関数であり、すべての確率p(ω)を足し合わせると1となるよう定義されている。

熱力学との関係

系が微視的状態 ω をとるとき、微視的な物理量O(ω) で与えられるとき、対応する熱力学的な状態量期待値

[math]O(\beta) = \langle O(\omega) \rangle = \sum_\omega O(\omega) p(\omega) = \frac{1}{Z(\beta)}\sum_\omega O(\omega) \mathrm{e}^{-\beta E(\omega)}[/math]

として再現される。特にエネルギーは

[math]E(\beta) = \frac{1}{Z(\beta)}\sum_\omega E(\omega) \mathrm{e}^{-\beta E(\omega)} = -\frac{\partial}{\partial\beta} \ln Z(\beta)[/math]

となる。

自由エネルギー F はエネルギーと

[math]E(\beta) = \frac{\partial}{\partial\beta}\{ \beta F(\beta) \}[/math]

の関係にあり、

[math]F(\beta) = -\frac{1}{\beta} \ln Z(\beta)[/math]

となる。これは、微視的な確率分布に基づく分配関数と熱力学的な状態量の自由エネルギーが関連付けられており、統計力学による熱力学の再現の一例である。

温度を引数にもつ自由エネルギーは完全な熱力学関数であり、ここから様々な状態量が計算される。例えばエントロピー

[math]S(\beta) = k\beta^2 \frac{\partial F(\beta)}{\partial\beta} = k\ln Z(\beta) -k\beta \frac{\partial}{\partial\beta} \ln Z(\beta)[/math]

となり、熱容量

[math]C(\beta) = -\beta\frac{\partial S(\beta)}{\partial\beta} = k\beta^2 \frac{\partial^2}{\partial\beta^2} \ln Z(\beta)[/math]

となる。

体積 V粒子数 N を考慮した系を考えると、圧力 P化学ポテンシャル μ なども統計力学的に表記できる。

[math]\begin{align} P(\beta,V,N) &= \frac{1}{\beta} \frac{\partial}{\partial V} \ln Z(\beta,V,N) \\ \mu(\beta,V,N) &= -\frac{1}{\beta} \frac{\partial}{\partial N} \ln Z(\beta,V,N) \end{align}[/math]


ボルツマンの原理

エントロピーは

[math]S(\beta) = k\beta \{ E(\beta) -F(\beta) \} = k\langle \beta E(\omega) +\ln Z(\beta) \rangle = -k \left\langle \ln \frac{\mathrm{e}^{-\beta E(\omega)}}{Z(\beta)} \right\rangle[/math]

となり、ボルツマンの公式

[math]S = -k\langle \ln p(\omega) \rangle[/math]

をみたす。

量子力学的な表記

量子力学的な系では、微視的状態はヒルベルト空間上の点で表される。特にエネルギー固有状態で代表することが多く、確率分布は

[math]p_i = \frac{1}{Z(\beta)} \mathrm{e}^{-\beta E_i}[/math]

となり、分配関数は

[math]Z(\beta) = \sum_i \mathrm{e}^{-\beta E_i}[/math]

となる。i はエネルギー固有状態を指定する量子数で、 Ei は対応するエネルギー固有値である。

トレースを用いると、分配関数はハミルトニアン テンプレート:Hat により、

[math]Z(\beta) = \mathrm{Tr}\, [\mathrm{e}^{-\beta \hat{H}}][/math]

と表せる。