グランドカノニカル分布

大正準分布（だいせいじゅんぶんぷ）、グランドカノニカル分布（英: grand canonical distribution）は、統計力学において系の微視的状態を表現する統計集団の一つである大正準集団、グランドカノニカルアンサンブル（英: grand canonical ensemble）に従う統計分布である。グランドカノニカル分布は熱浴（heat bath）との間でエネルギーを、粒子浴との間で粒子数を自由にやりとりできる系（開放系）に対応する統計分布で、小正準分布、正準分布とは体積が十分に大きい極限において熱力学的に等価である。

確率分布

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

[math]p(\omega) = \frac{1}{\Xi(\beta,\boldsymbol{\mu})} e^{-\beta E(\omega)+\beta \sum_i \mu_i N_i(\omega)}[/math]

ここで、E(ω) と N_i(ω) はそれぞれ系が微視的状態 ω をとるときのエネルギーと粒子数（i は粒子の種類）で、β μ_i はリザバーを特徴付けるパラメータでそれぞれ温度と化学ポテンシャルである。 β は絶対温度 T と β=1/kT の関係にあり、逆温度と呼ばれる。k はボルツマン定数である。

確率分布 p(ω) の分母に現れた規格化定数 Ξ(β,μ) はグランドカノニカル分布の大分配関数であり、次式で定義される。

[math]\Xi(\beta, \boldsymbol{\mu}) = \sum_\omega e^{-\beta E(\omega)+\beta \sum_i \mu_i N_i(\omega)}[/math]

熱力学との関係

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

[math]O(\beta, \boldsymbol{\mu}) = \langle O(\omega) \rangle = \sum_\omega O(\omega) p(\omega) = \frac{1}{\Xi(\beta, \boldsymbol{\mu})} \sum_\omega O(\omega) e^{-\beta E(\omega)+\beta \sum_i \mu_i N_i(\omega)}[/math]

として再現される。 特に粒子数は

[math]N_i(\beta, \boldsymbol{\mu}) = \frac{1}{\Xi(\beta, \boldsymbol{\mu})} \sum_\omega N_i(\omega) e^{-\beta E(\omega)+\beta \sum_i \mu_i N_i(\omega)} = \frac{1}{\beta}\frac{\partial}{\partial\mu_i} \ln \Xi(\beta,\boldsymbol{\mu})[/math]

となり、エネルギーは

[math]E(\beta, \boldsymbol{\mu}) = \frac{1}{\Xi(\beta, \boldsymbol{\mu})} \sum_\omega E(\omega) e^{-\beta E(\omega)+\beta \sum_i \mu_i N_i(\omega)} = -\frac{\partial}{\partial\beta} \ln \Xi(\beta,\boldsymbol{\mu}) +\sum_i \frac{\mu_i}{\beta}\frac{\partial}{\partial\mu_i} \ln \Xi(\beta,\boldsymbol{\mu})[/math]

となる。

系のグランドポテンシャルを

[math]J(\beta,\boldsymbol{\mu}) = -\frac{1}{\beta} \ln \Xi(\beta, \boldsymbol{\mu})[/math]

として定義すると、グランドポテンシャルは完全な熱力学関数であり、カノニカル分布における自由エネルギーと同様に、他の状態量を計算することができる。

量子理想気体

グランドカノニカル分布は粒子が生成・消滅する系でも使えるため、場の量子論における量子理想気体の平衡状態について記述する際に便利である。

理想気体なので粒子間の相互作用が無く、一粒子のエネルギー固有状態を考えればよい。 一粒子のエネルギー固有状態 j にある粒子数を n_j とし、対応する一粒子のエネルギー固有値を ε_j とすると、微視的状態 ω は粒子数 n_j の組によって指定される。

[math]N(n_1, \ldots, n_j, \ldots ) = \sum_{j=1}^\infty n_j,~ E(n_1,\ldots, n_j, \ldots ) = \sum_{j=1}^\infty \epsilon_j n_j[/math]

これはグランドカノニカル分布においては、全エネルギー及び全粒子数について拘束条件が無い（一定である必要が無い）為に行える操作であり、大分配関数は次のように書き直せる。

[math]\Xi(\beta, \mu) = \sum_{n_1} \cdots \sum_{n_j} \cdots \left( \prod_{j=1}^\infty \exp[-\beta (\epsilon_j-\mu) n_j] \right) = \prod_{j=1}^\infty \left( \sum_{n_j} \exp[-\beta (\epsilon_j-\mu) n_j] \right) = \prod_{j=1}^\infty \Xi^{(j)}(\beta, \mu)[/math]

このように、全体の大分配関数を固有状態 j の大分配関数の各々の積として表せる。 これが、グランドカノニカル分布が他の統計分布と比べて量子理想気体を記述する際に使い勝手の良い理由である。

ボゾン

一粒子のエネルギー固有値 ε_j をもつ固有状態jについて、ボゾンの場合、粒子数 n_j は 0 以上の全ての整数値をとりうるので大分配関数は、

[math]\Xi^{(j)}(\beta, \mu) = \sum_{n_j=0}^\infty e^{-\beta(\epsilon_j -\mu) n_j } = \frac{1}{1 - e^{-\beta(\epsilon_j -\mu) } }[/math]

となる。これから固有状態jの粒子数（占有数）の期待値を計算する。これはマクロな観測量では無いが、期待値を求めておくと量子理想気体などの解析に便利である。結果は、

[math]\langle n_j \rangle = \frac{1}{\beta} \frac{\partial}{\partial\mu} \ln\Xi^{(j)}(\beta, \mu) = -\frac{1}{\beta} \frac{\partial}{\partial\mu} \ln(1 -e^{-\beta(\epsilon_j -\mu) }) = \frac{1}{e^{\beta(\epsilon_j -\mu) } -1}[/math]

となる。これがボース分布関数である。

フェルミオン

一粒子のエネルギー固有値ε_jをもつ固有状態jについて、フェルミオンの場合、粒子数n_jは0もしくは1のみをとるので大分配関数は、

[math]\Xi^{(j)}(\beta, \mu) = \sum_{n_j=0}^1 e^{-\beta(\epsilon_j -\mu) n_j } = 1 + e^{-\beta(\epsilon_j -\mu) }[/math]

となる。これから粒子数の期待値を計算すると、

[math]\langle n_j \rangle = \frac{1}{\beta} \frac{\partial}{\partial\mu} \ln\Xi^{(j)}(\beta, \mu) = \frac{1}{\beta} \frac{\partial}{\partial\mu} \ln(1 +e^{-\beta(\epsilon_j -\mu) }) = \frac{1}{e^{\beta(\epsilon_j -\mu) } +1}[/math]

となる。これがフェルミ分布関数である。

量子力学的な表記

ヒルベルト空間の正規直交基底をe_iとして、任意の演算子[math]\hat{A}[/math]のトレースを、

[math] \mathbf{Tr} [ \hat{A} ] = \sum_i \langle e_i | \hat{A} | e_i \rangle[/math]

と定義する。これを用いると大分配関数はハミルトニアン[math]\hat{H}[/math]と粒子数演算子[math]\hat{N}[/math]を用いて、

[math] \Xi(\beta , \mu) = \mathbf{Tr} [ e^{- \beta ( \hat{H} - \mu \hat{N})} ][/math]

と表せる。

関連項目

• 化学ポテンシャル