氣體動力論與分子速率分布

從「一群碰撞的小球」到統計分布

一個盛滿氣體的容器看似靜止，內部卻是約 $10^{23}$ 個分子以每秒數百公尺的速率不斷碰撞、飛行的混沌劇場。氣體動力論（kinetic theory of gases）的核心主張是：宏觀的壓力、溫度等量，其實只是這群分子微觀運動的統計平均。本篇要做的，不是停在「分子撞牆產生壓力」的圖像，而是嚴謹地從牛頓力學與機率論導出壓力公式、能量均分，並完整推導 Maxwell–Boltzmann 速率分布。

從動量傳遞導出壓力

考慮邊長 $L$ 的立方容器，內有 $N$ 個質量 $m$ 的分子。先看單一分子沿 $x$ 方向以速度分量 $v_x$ 撞擊右壁。彈性碰撞使動量改變 $\Delta p_x = 2 m v_x$。兩次撞同一面壁之間，分子需往返一趟，時間 $\Delta t = 2L / v_x$。因此該分子對右壁的平均施力為

$$ f_x = \frac{\Delta p_x}{\Delta t} = \frac{2 m v_x}{2L / v_x} = \frac{m v_x^2}{L}. $$

對 $N$ 個分子求和，總力除以面積 $A = L^2$ 即得壓力。引入體積 $V = L^3$：

$$ P = \frac{\sum_i m v_{x,i}^2}{L^3} = \frac{N m}{V} \langle v_x^2 \rangle, $$

其中 $\langle v_x^2 \rangle$ 為 $v_x^2$ 的系綜平均。由於空間各向同性，$\langle v_x^2 \rangle = \langle v_y^2 \rangle = \langle v_z^2 \rangle$，而 $\langle v^2 \rangle = \langle v_x^2 \rangle + \langle v_y^2 \rangle + \langle v_z^2 \rangle = 3 \langle v_x^2 \rangle$。於是

$$ PV = \frac{1}{3} N m \langle v^2 \rangle = \frac{2}{3} N \left( \frac{1}{2} m \langle v^2 \rangle \right). $$

對照理想氣體狀態方程 $PV = N k_B T$，立刻讀出平均平動動能與溫度的對應：

$$ \frac{1}{2} m \langle v^2 \rangle = \frac{3}{2} k_B T. $$

這正是能量均分定理（equipartition theorem）在三個平動自由度上的展現：每個平方項自由度貢獻 $\tfrac{1}{2} k_B T$。

Maxwell–Boltzmann 速度分布的推導

上面只用到二階矩 $\langle v^2 \rangle$，但分子速率究竟如何「分布」？Maxwell 的精妙論證如下。設速度分布函數可分離為三個方向的乘積，且僅依賴速率大小（各向同性）：

$$ f(\vec{v}) = g(v_x)\, g(v_y)\, g(v_z) = \phi(v_x^2 + v_y^2 + v_z^2). $$

要同時滿足「乘積形式」與「只依賴 $v^2$」，唯一的函數形式是指數函數。取對數後 $\ln g(v_x) + \ln g(v_y) + \ln g(v_z)$ 須為 $v_x^2+v_y^2+v_z^2$ 的函數，故 $\ln g(v_x) = a + b v_x^2$，即

$$ g(v_x) = C\, e^{-\alpha v_x^2}, $$

其中 $\alpha > 0$（否則無法歸一化）。由能量考量與配分函數可定出 $\alpha = m / (2 k_B T)$。歸一化 $\int_{-\infty}^{\infty} g(v_x)\, dv_x = 1$ 配合高斯積分 $\int e^{-\alpha x^2} dx = \sqrt{\pi/\alpha}$，得 $C = \sqrt{m/(2\pi k_B T)}$。

要得到速率（而非速度向量）的分布，需把三維速度空間轉為球座標。速率介於 $v$ 與 $v + dv$ 的殼層體積為 $4\pi v^2\, dv$，因此

$$ \boxed{\,F(v) = 4\pi \left( \frac{m}{2\pi k_B T} \right)^{3/2} v^2\, e^{-\frac{m v^2}{2 k_B T}}\,} $$

這就是 Maxwell–Boltzmann 速率分布。$v^2$ 因子使 $F(0)=0$，指數因子壓抑高速尾端，兩者競爭造就一個非對稱的峰。

由 $F(v)$ 可導出三個常用特徵速率。令 $\dfrac{dF}{dv}=0$ 得最可能速率

$$ v_p = \sqrt{\frac{2 k_B T}{m}}. $$

對 $v$、$v^2$ 取期望值（利用 $\int_0^\infty v^n e^{-\alpha v^2} dv$ 的標準結果）分別得平均速率與方均根速率：

$$ \bar{v} = \sqrt{\frac{8 k_B T}{\pi m}}, \qquad v_{\text{rms}} = \sqrt{\frac{3 k_B T}{m}}. $$

三者比值固定：$v_p : \bar{v} : v_{\text{rms}} = \sqrt{2} : \sqrt{8/\pi} : \sqrt{3} \approx 1 : 1.128 : 1.225$，恆有 $v_p < \bar{v} < v_{\text{rms}}$，正反映分布的右偏特性。

定量小範例：室溫氮氣的方均根速率

求 $T = 300\,\text{K}$ 時氮氣 $\mathrm{N_2}$ 的 $v_{\text{rms}}$。氮分子莫耳質量 $M = 28.0\,\text{g/mol} = 28.0\times10^{-3}\,\text{kg/mol}$，單分子質量

$$ m = \frac{M}{N_A} = \frac{28.0\times10^{-3}}{6.022\times10^{23}} \approx 4.65\times10^{-26}\,\text{kg}. $$

代入 $v_{\text{rms}} = \sqrt{3 k_B T / m}$，其中 $k_B = 1.38\times10^{-23}\,\text{J/K}$：

$$ v_{\text{rms}} = \sqrt{\frac{3 \times (1.38\times10^{-23}) \times 300}{4.65\times10^{-26}}} = \sqrt{\frac{1.242\times10^{-20}}{4.65\times10^{-26}}}. $$

括號內 $\approx 2.67\times10^{5}\,\text{m}^2/\text{s}^2$，開根號得

$$ v_{\text{rms}} \approx 517\,\text{m/s}. $$

這比室溫聲速（約 $346\,\text{m/s}$）略大，符合物理直覺——聲速本就源於分子熱運動傳遞擾動，且約為 $v_{\text{rms}}$ 的 $\sqrt{\gamma/3}$ 倍。對應的平均速率 $\bar v = v_{\text{rms}}\sqrt{8/(3\pi)} \approx 476\,\text{m/s}$。

平均自由徑與輸運的橋樑

速率分布還不足以描述「碰撞頻率」。把分子視為直徑 $d$ 的硬球，一個分子掃過的碰撞截面為 $\sigma = \pi d^2$。考慮相對運動（其他分子也在動），引入 $\sqrt{2}$ 修正因子，平均自由徑（mean free path）為

$$ \lambda = \frac{1}{\sqrt{2}\, n \sigma} = \frac{k_B T}{\sqrt{2}\, \pi d^2 P}, $$

其中 $n = N/V$ 為數密度。$\lambda$ 與壓力成反比：常壓下空氣分子 $\lambda \sim 70\,\text{nm}$，高真空下可達公尺量級。$\lambda$ 與 $\bar v$ 進一步決定黏滯係數、熱導率與擴散係數，是連結微觀分布與宏觀輸運現象的關鍵橋樑。

深入探討（研究所視角）

上述推導本質上是平衡態統計，更深的框架來自 Ludwig Boltzmann 的輸運方程。Boltzmann 方程描述單粒子分布函數 $f(\vec{r}, \vec{v}, t)$ 的演化：

$$ \frac{\partial f}{\partial t} + \vec{v}\cdot\nabla_{\vec r} f + \frac{\vec{F}}{m}\cdot\nabla_{\vec v} f = \left( \frac{\partial f}{\partial t} \right)_{\text{coll}}, $$

左側為相空間中沿哈密頓流的對流（Liouville 定理的單粒子約化），右側碰撞積分項以二體碰撞核耦合不同速度態。Maxwell–Boltzmann 分布正是使碰撞項恆為零的局部平衡解，這也是 Boltzmann 的 $H$ 定理所保證的終態：定義 $H = \int f \ln f\, d^3v$，可證 $dH/dt \le 0$，等號僅在分布為 Maxwell 形式時成立。$H$ 與熵 $S = -k_B H$ 的對應，使動力論成為熱力學第二定律的微觀奠基，也引出時間箭頭與可逆微觀動力學之間著名的 Loschmidt 佯謬。

在哈密頓表述下，氣體可視為 $6N$ 維相空間 $(\{q_i\}, \{p_i\})$ 中的一點，其演化遵循正則方程 $\dot q_i = \partial \mathcal{H}/\partial p_i$、$\dot p_i = -\partial \mathcal{H}/\partial q_i$。Liouville 定理斷言相空間機率密度沿軌跡守恆（不可壓縮流），這是整個統計力學系綜理論的支柱。Maxwell–Boltzmann 分布在此框架下對應正則系綜（canonical ensemble）的速度邊際分布，其權重 $e^{-\beta \mathcal{H}}$ 中的 $\beta = 1/(k_B T)$ 由配分函數 $Z = \int e^{-\beta \mathcal{H}}\, d\Gamma$ 統攝全部熱力學量——例如自由能 $F = -k_B T \ln Z$。

古典結果在低溫或高密度下失效，須讓位給量子統計。當熱波長 $\lambda_{\text{dB}} = h/\sqrt{2\pi m k_B T}$ 接近粒子間距時，全同粒子的不可分辨性主導行為：玻色子服從 Bose–Einstein 分布、費米子服從 Fermi–Dirac 分布，兩者在 $\lambda_{\text{dB}} \ll n^{-1/3}$ 的稀薄高溫極限收斂回 Maxwell–Boltzmann。這解釋了金屬中電子氣的簡併壓、白矮星的穩定，以及 1995 年實現的玻色–愛因斯坦凝聚。

理論連結上，動力論可由更基本的 BBGKY 階層（描述 $s$ 體分布函數的耦合方程組）系統地截斷導出，Boltzmann 方程即是在「分子混沌假設」（Stosszahlansatz，假設碰撞前兩粒子無關聯）下的二體截斷。當代發展更把這套框架推廣到電漿（Vlasov 方程，含長程庫侖場與集體模）、稀薄氣動力學（太空船重返大氣的 DSMC 數值模擬）乃至相對論流體與重離子碰撞的夸克–膠子電漿輸運。從一群碰撞小球出發，動力論最終成為連接微觀力學與宏觀連續體、古典與量子的核心理論支點。