光子气体与黑体辐射

在统计力学的框架下，光子气体（Photon Gas）是一个极其重要的理想模型。它描述了处于热平衡状态下的电磁辐射场，即在封闭容器（黑体腔）内，光子与容器壁不断进行吸收与发射过程，最终达到的一种宏观稳定的状态。

光子气体的研究不仅是理解黑体辐射（Blackbody Radiation）的核心，也是连接量子力学与热力学的重要桥梁。与传统的理想气体不同，光子气体具有独特的统计特性，例如其粒子数是不守恒的，且其化学势始终为零。本文将从统计分布、态密度、能量密度以及热力学性质四个维度，深入探讨光子气体与黑体辐射的理论构建。

光子的统计特性与化学势

光子是自旋为 1 的玻色子（Boson），因此其遵循玻色-爱因斯坦统计（Bose-Einstein Statistics）。在统计力学中，描述系统粒子分布的函数为：

$$\langle n_i \rangle = \frac{1}{e^{\beta(\epsilon_i - \mu)} - 1}$$

其中，$\beta = 1/k_B T$，$\epsilon_i$ 是能级能量，$\mu$ 是化学势。

对于光子气体，有一个至关重要的物理事实：光子的数量不是守恒的。在黑体腔内，光子可以被容器壁吸收，也可以由容器壁发射。为了使系统的自由能达到极小值，在热力学平衡条件下，必须满足粒子数 $N$ 对自由能 $F$ 的偏导数为零，即：

$$\left( \frac{\partial F}{\partial N} \right)_{T, V} = \mu = 0$$

因此，光子气体的分布函数简化为：

$$\langle n(\omega) \rangle = \frac{1}{e^{\hbar \omega / k_B T} - 1}$$

这一特性使得光子气体的热力学行为完全由温度 $T$ 决定，而不需要考虑粒子数的变化。

态密度（Density of States）的推导

要计算光子气体的宏观物理量，首先需要确定在给定的频率范围内，系统允许存在的量子态数量，即态密度 $g(\omega)$。

考虑一个边长为 $L$ 的立方体黑体腔，根据电磁波的边界条件，其波矢 $\mathbf{k}$ 是量子化的。在 $k$ 空间中，每个量子态占据的体积约为 $(\pi/L)^3$。由于电磁波具有两个独立的偏振方向，因此在 $k$ 空间中，单位体积内的态数为：

$$g(k) dk = 2 \times \left( \frac{L}{\pi} \right)^3 \times \frac{1}{8} \times 4\pi k^2 dk = \frac{V k^2}{\pi^2} dk$$

利用光子的色散关系 $\omega = ck$（其中 $c$ 为光速），我们可以将关于 $k$ 的表达式转化为关于频率 $\omega$ 的表达式：

$$g(\omega) d\omega = \frac{V \omega^2}{\pi^2 c^3} d\omega$$

这里的 $g(\omega)$ 即为单位体积内的态密度。可以看到，态密度随频率的平方增长，这意味着高频光子拥有更多的可用状态。

普朗克辐射定律

结合玻色-爱因斯坦分布和态密度，我们可以推导出黑体辐射的能量谱密度。

首先，单个频率 $\omega$ 对应的能量为 $\hbar \omega$。单位体积内频率在 $\omega$ 到 $\omega + d\omega$ 之间的能量密度 $u(\omega, T)$ 为：

$$u(\omega, T) d\omega = \frac{1}{V} \cdot \hbar \omega \cdot g(\omega) \cdot \langle n(\omega) \rangle d\omega$$

代入之前的表达式：

$$u(\omega, T) = \frac{\hbar \omega^3}{\pi^2 c^3} \frac{1}{e^{\hbar \omega / k_B T} - 1}$$

这就是著名的普朗克定律（Planck's Law）。它成功地解决了经典物理学中的“紫外灾难”问题：在经典理论中，能量密度随频率趋于无穷大，而普朗克定律通过引入量子化能级 $\hbar \omega$，使得在高频区域（$\hbar \omega \gg k_B T$）能量密度呈指数级衰减，从而保证了物理上的自洽性。

热力学性质与宏观定律

通过对普朗克分布进行积分，我们可以导出光子气体的宏观热力学规律。

1. 斯特藩-玻尔兹曼定律 (Stefan-Boltzmann Law)

计算单位体积内的总能量密度 $u(T)$：

$$u(T) = \int_{0}^{\infty} u(\omega, T) d\omega = \frac{\hbar}{\pi^2 c^3} \int_{0}^{\infty} \frac{\omega^3}{e^{\hbar \omega / k_B T} - 1} d\omega$$

通过变量代换 $x = \hbar \omega / k_B T$，可以得到：

$$u(T) = \left( \frac{k_B T}{\hbar c} \right)^4 \frac{\hbar}{\pi^2 c^3} \int_{0}^{\infty} \frac{x^3}{e^x - 1} dx$$

积分项 $\int_{0}^{\infty} \frac{x^3}{e^x - 1} dx = \frac{\pi^4}{15}$ 是一个常数。因此：

$$u(T) = a T^4$$

其中 $a = \frac{\pi^2 k_B^4}{15 \hbar^3 c^3}$。这表明黑体辐射的总能量密度与温度的四次方成正比，即斯特藩-玻尔兹曼定律。

2. 维恩位移定律 (Wien's Displacement Law)

通过寻找 $u(\omega, T)$ 函数的最大值点（即令 $\frac{\partial u}{\partial \omega} = 0$），可以得到波峰频率 $\omega_{max}$ 与温度的关系：

$$\frac{\hbar \omega_{max}}{k_B T} = \text{constant} \implies \lambda_{max} T = \text{constant}$$

这解释了为什么物体温度越高，发出的光颜色越偏向蓝紫色（短波长）。

3. 压强与能量密度的关系

对于光子气体，其压强 $P$ 与能量密度 $u$ 之间存在一个简单的比例关系。通过热力学基本方程或统计力学中的配分函数可以证明：

$$P = \frac{1}{3} u(T)$$

这一结论在相对论性流体动力学中非常关键，它反映了光子作为无质量粒子，其动量与能量之间存在 $p = E/c$ 的线性关系。

总结

光子气体模型是统计力学中处理非守恒粒子系统的典范。通过引入量子化概念，我们不仅完美解释了黑体辐射的谱分布，还推导出了斯特藩-玻尔兹曼定律和维恩位移定律等宏观热力学规律。理解光子气体的本质——即化学势为零的玻色气体——是深入研究量子统计、宇宙学背景辐射以及凝聚态物理中激子（Exciton）等准粒子行为的基础。