[統計力学] マクスウェル・ボルツマン分布

February 16, 2020      

エネルギーが位置と運動量によって表されるような古典的な系では,相空間の微小部分に関するカノニカル分布は以下のように表される.

dwn(p,q)=AeEn(p,q)/T dpdq\begin{equation} dw_n(p, q) = A e^{-E_n(p, q) / T}~dpdq \end{equation}

エネルギーE(p,q)E(p, q)は運動エネルギーK(p)K(p)と位置エネルギーU(q)U(q)の和として表され,それぞれの分布が正規化されているはずだ.

dwp=aeK(p)/T dpdwq=beU(q)/T dq\begin{gather} dw_p = ae^{-K(p)/T}~dp \\ dw_q = be^{-U(q)/T}~dq \end{gather}

ある質量mmの分子,運動量に関する確率分布は次のように表される.

dwp=(12πmkBT)32exp[px2+py2+pz22mkBT]dpxdpydpz\begin{equation} dw_p = \left( \frac{1}{2\pi m k_\mathrm{B} T} \right)^{\frac{3}{2}} \exp \left[ -\frac{p_x^2+p_y^2+p_z^2}{2mk_\mathrm{B}T} \right] dp_x dp_y dp_z \end{equation}

運動量を速度に書き換えれば,マクスウェル・ボルツマン分布が得られる.

dwv=(m2πkBT)32exp[m(vx2+vy2+vz2)2kBT]dvxdvydvz\begin{equation} % \label{eq:M-B} dw_v = \left( \frac{m}{2\pi k_\mathrm{B} T} \right)^{\frac{3}{2}} \exp \left[ -\frac{m(v_x^2+v_y^2+v_z^2)}{2k_\mathrm{B}T} \right] dv_x dv_y dv_z \end{equation}

単位体積内の分子の個数をかければ,個数分布に書き換えることが出来る.

dNv=NV(m2πkBT)32exp[m(vx2+vy2+vz2)2kBT]dvxdvydvz\begin{equation} dN_v = \frac{N}{V} \left( \frac{m}{2\pi k_\mathrm{B} T} \right)^{\frac{3}{2}} \exp \left[ -\frac{m(v_x^2+v_y^2+v_z^2)}{2k_\mathrm{B}T} \right] dv_x dv_y dv_z \end{equation}