ボルツマン分布
ボルツマン分布（ボルツマンぶんぷ、英語: Boltzmann distribution）は、一つのエネルギー準位にある粒子の数（占有数）の分布を与える理論式の一つである。ギブス分布とも呼ばれる。気体分子の速度の分布を与えるマクスウェル分布をより一般化したものに相当する。
量子統計力学においては、占有数の分布がフェルミ分布に従うフェルミ粒子と、ボース分布に従うボース粒子の二種類の粒子に大別できる。ボルツマン分布はこの二種類の粒子の違いが現れないような条件におけるフェルミ分布とボーズ分布の近似形（古典近似）である。ボルツマン分布に従う粒子は古典的粒子とも呼ばれる。
核磁気共鳴および電子スピン共鳴などにおいても、磁場の中で分裂した2つの準位の占有率はボルツマン分布に従う。
概要
ボルツマン分布に従う系において、エネルギーが ε に等しい一つの準位にある粒子の数は
[math]f(\epsilon) =\lambda\, \mathrm{e}^{-\beta\epsilon} =\mathrm{e}^{-\beta(\epsilon-\mu)}[/math]
で与えられる。分布関数を特徴付けるパラメータ β は系の温度と解釈され、熱力学温度 T と β=1/kT で関係付けられ、逆温度と呼ばれる。比例係数 λ は活量で、μ は化学ポテンシャルである。比例係数を除いた e−βε=e−ε/kT の項は、エネルギー ε をもつ粒子の割合を表し、ボルツマン因子と呼ばれる。エネルギーが ε の準位の占有数と ε+Δε の準位の占有数の比は
[math]\frac{f(\epsilon+\Delta\epsilon)}{f(\epsilon)} =\mathrm{e}^{-\beta \Delta\epsilon}[/math]
となる[1]。同じ温度では、高いエネルギー（大きな ε）の準位の方が一つの準位あたりの粒子数が小さくなる。また、同じエネルギーの準位でも、高い温度（小さな β、大きな T）の条件では一つの準位あたりの粒子数が大きくなる。
複雑な粒子間相互作用がなく、エネルギー準位の分布が占有数によって変化しないことを仮定する。エネルギーが ε と ε+dε の範囲にある準位の数を g(ε)dε とすれば、この範囲にある粒子の数は f(ε)g(ε)dε で与えられる。系の全粒子数は、全てのエネルギーの範囲で積分して
[math]N =\int_{-\infty}^\infty f(\epsilon)\, g(\epsilon)\, d\epsilon =\lambda \int_{-\infty}^\infty \mathrm{e}^{-\beta\epsilon} g(\epsilon)\, d\epsilon[/math]
で与えられる。また、系の全エネルギーは
[math]E =\int_{-\infty}^\infty \epsilon\, f(\epsilon)\, g(\epsilon)\, d\epsilon =\lambda \int_{-\infty}^\infty \epsilon\, \mathrm{e}^{-\beta\epsilon} g(\epsilon)\, d\epsilon[/math]
で与えられる。
エネルギー準位の分布が離散的な場合は、エネルギーが εi に等しい準位の数を gi として、エネルギーが εi である粒子の数 ni は
[math]n_i =f_i g_i =\lambda\, \mathrm{e}^{-\beta\epsilon_i} g_i[/math]
となり、系の全粒子数と全エネルギーは
[math]N =\sum_i n_i =\lambda \sum_i \mathrm{e}^{-\beta\epsilon_i} g_i ,\quad E =\sum_i \epsilon_i\, n_i =\lambda \sum_i \epsilon_i\, \mathrm{e}^{-\beta\epsilon_i} g_i[/math]
で与えられる。
ボルツマン分布は気体の温度が充分に高く、密度が充分に低く、かつ量子効果が無視されるような系において適用される。βε=ε/kT が大きな値を取るような場合、もしくは状態密度が小さい場合のように古典的粒子として扱うには限界が生じ、かつ粒子の波動関数が実質的に重複していない場合は、ボース=アインシュタイン分布およびフェルミ=ディラック分布の両方がボルツマン分布になる。
例
理想気体
分子のエネルギーは単純に粒子の運動エネルギーで与えられる。
- [math]E_i = {\begin{matrix} \frac{1}{2} \end{matrix}} mv^{2}[/math]
- [math]E_i = {\begin{matrix} \frac{1}{2} \end{matrix}} mv^{2} + mgh[/math]
この場合の気体分子の垂直分布は以下の式で表される。
- [math]{{N_i}\over{N}} = {{g_i e^{-mgh/(k_BT)}}\over{Z(T)}}[/math]
脚注
参考文献
- Gordon M. Barrow 『物理化学』 大門寛、堂免一成 訳、東京化学同人、1999年。
- Tai L. Chow 『科学技術者のための数学ハンドブック』 朝倉書店、2002年。ISBN 4-254-11090-1。