ボーズ・アインシュタイン分布関数の導出


 ボーズ・アインシュタイン分布関数に従う粒子は、スピンが整数であるボーズ粒子である。 一方、スピンが半整数(例えばスピン1/2の電子など)のフェルミ粒子の場合はフェルミ分布関数に従い、導出にはパウリの排他原理を用いた。 ボーズ粒子とフェルミ粒子の統計性の違いが分布関数の違いを生んでいる。

 ここでは、以下のボーズ・アインシュタイン分布関数を導出する。

    \begin{eqnarray*} f(E)=\frac{1}{\exp{\left( \frac{E-\mu}{k_{\rm B}T} \right)}-1} \end{eqnarray*}

【参考】フェルミ分布関数の導出


ボーズ・アインシュタイン分布関数の導出

あるエネルギー状態数

 異なる状態でも同一のエネルギーになる。そのことは、自由電子のエネルギー E=\frac{\hbar^2}{2m}(k_x^2 + k_y^2 +k_z^2) について、{\bf k}=(1\, 0\, 0){\bf k}=(0\, 1\, 0) が同じエネルギーになることからわかる。 そのことはボーズ粒子の場合も同様である。

 はじめに、エネルギー E_i にある状態数 w_i を求めるため、 とりうる状態を Z_i として、ボーズ粒子が N_i 個ある場合を考える。フェルミ粒子はパウリの排他原理を満たしたが、ボーズ粒子はそのような制限がないため、複数のボーズ粒子が同一の状態をとってもよい


 かんたんな例として \textcolor{red}{6} (=N_i) 個のボールを \textcolor{blue}{3}(=Z_i) 個の袋にいれる問題を考える。 このとき、1個の袋に2個以上のボールが入っていても良い。 ボールを袋に詰める詰め方は図のように x^{\textcolor{red}{6}} の係数を求める問題に帰着する。

つまり、

    \begin{eqnarray*} (x^0+x^1+x^2+x^3+\cdots)^{\textcolor{blue}{3}} \end{eqnarray*}

を展開したときの x^{\textcolor{red}{6}} の係数を求めることと同じである(重複組み合わせによる別の考え方はページ一番下の補足に書いた)。

 一般化すると、

    \begin{eqnarray*} (x^0+x^1+x^2+x^3+\cdots)^{\textcolor{blue}{Z_i}}\quad\cdots\quad(1) \end{eqnarray*}

における x^{\textcolor{red}{N_i}} の係数が求める w_i になる。

式(1) のカッコ内は無限等比級数(初項1, 等比 x)であるため

    \begin{eqnarray*} 1+x+x^2+x^3+\cdots = \frac{1}{1-x} \end{eqnarray*}

と書ける。これより

    \begin{eqnarray*} f(x)&\equiv& (x^0+x^1+x^2+x^3+\cdots)^{\textcolor{blue}{Z_i}}=\left( \frac{1}{1-x} \right)^{\textcolor{blue}{Z_i}}\\ &=&(1-x)^{-\textcolor{blue}{Z_i}} \end{eqnarray*}

となる。


 さらに f(x)x=0 まわりでテイラー展開^{[*]} して x^{\textcolor{red}{N_i}} の係数を求める。そのために、f(x)N_i階微分を求めておく。

    \begin{eqnarray*} f'(x)&=&Z_i \left( 1-x \right)^{-Z_i-1}\\ f''(x)&=&Z_i(Z_i+1) \left( 1-x \right)^{-Z_i-2}\\ \cdots\\ f^{N_i}&=&Z_i(Z_i+1)\cdots(Z_i+N_i-1) \left( 1-x \right)^{-Z_i-N_i} \end{eqnarray*}

 これより、f(x) のテイラー展開は

    \begin{eqnarray*} f(x)&=&f(0)+xf'(0)+\frac{1}{2!}f''(0)+\cdots+\frac{1}{N_i!}f^{(N_i)}(0)+\cdots\\ &=&1+Z_i x + \frac{Z_i(Z_i+1)}{2!}+\cdots \\ &&\quad + \textcolor{red}{\frac{Z_i(Z_i+1)+\cdots (Z_i+N_i-1)}{N_i!} x^{N_i}}+\cdots \end{eqnarray*}

となり、x^{\textcolor{red}{N_i}} の係数が

    \begin{eqnarray*} \frac{Z_i(Z_i+1)+\cdots (Z_i+N_i-1)}{N_i!}=\frac{(Z_i+N_i-1)!}{N_i! (Z_i-1)!} \end{eqnarray*}

となる。(これは「考え方2」でみるように組み合わせの総数 _{Z_i+N_i-1}C_{N_i} に対応する。)


* x=0 におけるテイラー展開(マクローリン展開)

    \begin{eqnarray*} f(x)&=&f(0)+xf'(0)+\frac{1}{2!}f''(0)+\cdots+\frac{1}{N!}f^{(N)}(0)+\cdots\\ &=&\sum_{k=0}^{\infty} \frac{1}{k!}f^{(k)}(0)\\ &&\left(f^{(k)}(0)=\frac{d^k f(x)}{dx^k}\Big|_{x=0}\right) \end{eqnarray*}


全状態の数を求める

 これが i 番目のエネルギー E_i にある状態数である。 i についての総乗が、全てのエネルギーを考えたときの「占有のしかたの総数 W」 である。

    \begin{eqnarray*} W=\Pi_i w_i=\Pi_i\, \frac{(Z_i+N_i-1)!}{N_i! (Z_i-1)!} \end{eqnarray*}


エントロピーを最大にする

 総数 W がわかるとボルツマンの原理よりエントロピー S がわかる。

    \begin{eqnarray*} S=k_{\rm B} \ln W \end{eqnarray*}

 ここで、k_{\rm B} はボルツマン定数である。 熱力学によると、系が熱平衡状態にあるときエントロピー S最大になる(非平衡であれば S は増大する)。 \ln は単調増加関数であるため、S を最大化するためには \ln W\quad (W) を最大化すれば良い。 つまり、次に示すような束縛条件のもと、\ln W を最大化していく。


ラグランジュの未定乗数法を用いる

 全粒子数 N 一定、全エネルギー E 一定の条件

    \begin{eqnarray*} \begin{cases} N=\sum_i N_i\\ E=\sum_i N_i E_i \end{cases} \end{eqnarray*}

を課して、\ln W を最大化する。 そのためにラグランジュの未定乗数法で求めていく。

 2つのラグランジュの未定乗数を \alpha,\beta と置いて、

    \begin{eqnarray*} F(\alpha,\beta,N_i)=\ln W +\alpha(N-\sum_i N_i )+\beta(E-\sum_i N_i E_i) \end{eqnarray*}

を最大にする。

    \begin{eqnarray*} \frac{\partial F(\alpha,\beta,N_i)}{\partial N_i} =0 \quad \cdots (2) \end{eqnarray*}

N_i に関して上に凸の関数であるため、この条件で最大値を求めることができる。

ここで、スターリングの公式を用いて \ln W を近似しておく。


ポイント

スターリングの公式は以下の通りである。Nが十分大きい時、

    \begin{eqnarray*} \ln N! \simeq N\ln N - N \end{eqnarray*}

と近似できる。


 これを用いて \ln W を簡単にする。

    \begin{eqnarray*} \ln W&=&\ln \left[\Pi_i \frac{(Z_i+N_i-1)!}{N_i!(Z_i-1)!}\right]\\ &=&\sum_i \ln\left[\frac{Z_i+N_i-1!}{N_i!(Z_i-1)!}\right]\\ &\textcolor{red}{\simeq}&\sum_i\Bigg[(Z_i+N_i-1)\ln (Z_i+N_i-1) \cancel{-(Z_i+N_i-1)}\\ &&\quad - N_i\ln N_i + \cancel{N_i} -(Z_i-1) \ln (Z_i-1) + \cancel{Z_i-1}\Bigg]\\ &=&\sum_i\Bigg[(Z_i+N_i-1)\ln (Z_i +N_i-1)\\ &&\quad - N_i\ln N_i -(Z_i-1) \ln (Z_i-1)\Bigg] \end{eqnarray*}

これを N_i で偏微分すると、

    \begin{eqnarray*} \frac{\partial \ln W}{\partial N_i}&=&\ln(Z_i+N_i-1)+1-\ln N_i -1 \\ &=&\ln \frac{Z_i+N_i-1}{N_i} \end{eqnarray*}

 したがって、式(2)の N_i についての偏微分は

    \begin{eqnarray*} \frac{\partial F(\alpha,\beta,N_i)}{\partial N_i} &=&\ln \frac{Z_i+N_i-1}{N_i} - \alpha -\beta E_i \end{eqnarray*}

(\alpha,\beta の項の偏微分の計算については補足に記載^{[1]})


 よって、

    \begin{eqnarray*} &&\ln\frac{Z_i+N_i-1}{N_i}=\alpha+\beta E_i\\ \Leftrightarrow \quad && \frac{Z_i}{N_i}=\exp\left(\alpha+\beta E_i\right)-1+\cancel{\frac{1}{N_i}}\\ \Leftrightarrow \quad && \frac{N_i}{Z_i}= \frac{1}{\exp\left(\alpha+\beta E_i\right)-1} \end{eqnarray*}

となる。途中で状態数、粒子数が多いという条件 1 <  <  Z_i,N_i より \frac{1}{N_i} を無視した。 得られた結果は考えているエネルギー E_i をとる状態数 Z_i に対して、ボーズ粒子が粒子数 N_i だけ占有している場合の確率を表す


E_i連続変数として扱えば、

    \begin{eqnarray*} f(E)=\frac{1}{\exp\left(\alpha+\beta E\right)-1}\quad \cdots (3) \end{eqnarray*}

がボーズ・アインシュタイン分布関数となる。

\alpha,\beta については熱力学の関係式から

    \begin{eqnarray*} \begin{cases} \alpha = \frac{-\mu}{k_{\rm B}T}\\ \beta = \frac{1}{k_{\rm B}T} \end{cases} \end{eqnarray*}

となる(参考:フェルミ分布関数の導出[熱力学の式から定数を求める])。\mu は系の化学ポテンシャルである。

 この\alpha,\beta を用いて、式(3)からボーズ・アインシュタイン分布関数がわかる。

    \begin{eqnarray*} f(E)=\frac{1}{\exp{\left( \frac{E-\mu}{k_{\rm B}T} \right)}-1} \end{eqnarray*}



補足[1]: N_i に関する偏微分のところを簡単に書いておく。

    \begin{eqnarray*} \frac{\partial \sum_i N_i E_i}{\partial N_i} \end{eqnarray*}

について、\sum_i における i はダミーであるので、本当は

    \begin{eqnarray*} \frac{\partial \sum_j N_j E_j}{\partial N_i} \end{eqnarray*}

とすべきであった。そうすると、

    \begin{eqnarray*} \sum_j N_j E_j &=& N_0 E_0 + N_1 E_1 + \cdots + N_i E_i + \cdots \end{eqnarray*}

であるため、

    \begin{eqnarray*} \frac{\partial \sum_j N_j E_j}{\partial N_i} = E_i \end{eqnarray*}

になる。


プランクの分布関数

 ラグランジュの未定乗数 \alpha は粒子数一定の条件からあらわれた。 粒子数一定の条件を課さない場合は分布関数中に \alpha はあらわれない。したがって、 ボーズ・アインシュタイン分布関数は

    \begin{eqnarray*} f(E)=\frac{1}{\exp\left(\frac{E}{k_{\rm B}T}\right)-1} \end{eqnarray*}

プランクの分布関数になる。


 フェルミ分布関数の形

    \begin{eqnarray*} f(E)=\frac{1}{\exp\left(\frac{E-\mu}{k_{\rm B}T}\right)+1} \end{eqnarray*}

から、フェルミ分布関数は 0\leq f(E)\leq 1 であった。 これはパウリの排他原理より、ある状態にフェルミ粒子が存在する確率が1より大きくならないことをあらわす(パウリの排他原理を反映)。

 一方で、ボーズ・アインシュタイン分布関数あるいはプランクの分布関数では \exp(x)-1 の分母をもつため、f(E) が1より大きくなることがある。 特に高温で k_{\rm B}T が大きい場合は分母が0に近く(\expが1に近く)ため、粒子数が増大する。


状態数を求める考え方2

\textcolor{red}{6} (=N_i) 個のボールを \textcolor{blue}{3}(=Z_i) 個の袋にいれる問題について別の考え方を書いておく。 これは6個のボールを3つの袋に分ける重複組合せである(この名前より意味が大切である)。

 とにかく6つを3つに分けるため、6個のボールと2 (3-1) 個の仕切り棒を用意すればよい(図)。

 これを並び替えて

    \begin{eqnarray*} \frac{8!}{2!6!} \end{eqnarray*}

が組み合わせの総数となる。

 一般化すると、N_i 個のボールと Z_i-1 個の棒の並びの総数を求めれば良い。つまり

    \begin{eqnarray*} w_i=\frac{(Z_i+N_i-1)}{N_i! (Z_i-1)!} \end{eqnarray*}

となる。




コメントを残す

メールアドレスが公開されることはありません。