量子危机的诞生：黑体辐射与能量子概念的历史考察

基尔霍夫与黑体辐射问题的起源

1859年深秋，柏林科学院收到了一篇来自海德堡的论文，作者是年仅26岁的 Gustav Kirchhoff。这篇题为《关于可燃物燃烧所产生火焰的颜色与燃烧产物发光性质之间的关系》的论文，在当时并未引起广泛关注。然而，正是这篇看似普通的论文，奠定了热辐射理论的基石，并在四十一年后催生了物理学史上最伟大的革命之一——量子论的诞生。

Kirchhoff 在这篇论文中提出了一个后来被称为基尔霍夫辐射定律（Kirchhoff's law of thermal radiation）的重要原理：对于任何处于热平衡的物体，其发射能力和吸收能力的比值是一个普适函数，仅取决于辐射的波长和温度，而与物体的具体材质无关。这一定律的数学表述为：

$$\frac{E_1(\lambda, T)}{A_1(\lambda, T)} = \frac{E_2(\lambda, T)}{A_2(\lambda, T)} = f(\lambda, T)$$

其中 $E_i(\lambda, T)$ 表示物体 $i$ 在温度 $T$ 下对波长 $\lambda$ 的发射能力（emissive power），$A_i(\lambda, T)$ 表示其吸收能力（absorptivity）。这个普适函数 $f(\lambda, T)$ 不依赖于任何物质的特殊属性，它反映的是自然界中热辐射的一般规律。

黑体：理想辐射体的概念

为了研究这个普适函数的具体形式，物理学家需要引入一个理想化的参照物——黑体（black body）。黑体的定义源于其完美的吸收特性：在任何波长上，黑体的吸收能力都等于1，即 $A_{\text{black}}(\lambda, T) = 1$。这意味着黑体能够吸收所有照射到其表面的辐射，不反射也不透射任何光线——正因如此，它在肉眼看来呈现为完美的黑色。

然而，黑体的意义远不止于“完美的吸收体”。更具洞察力的观察是：由于处于热平衡的黑体既不反射外界辐射，也不透射外界辐射，它所发出的热辐射必然完全来自其自身内部的热运动。因此，黑体辐射的光谱分布直接反映了热平衡条件下电磁辐射与物质相互作用的最一般规律。基尔霍夫定律的推论因此具有了深刻的物理含义：

$$f(\lambda, T) = E_{\text{black}}(\lambda, T)$$

这意味着基尔霍夫所寻找的普适函数，恰恰就是黑体在热平衡状态下发射的辐射强度分布。找到这个函数的数学形式，就等于揭示了热辐射的终极规律。

黑体模型的实验探索

十九世纪下半叶，随着冶金工业和光谱分析技术的发展，实验物理学家开始系统测量黑体辐射的光谱分布。最具影响力的实验装置是由 Wilhelm Otto von Lummer 和 Ernst Pringsheim 在1895年前后完善的空腔辐射器（cavity radiator）。这一装置的核心思想极为精妙：如果在一个带有小孔的密封腔体内部壁涂上煤烟或其他强吸收材料，那么任何从小孔进入腔体的辐射，都将在腔内经历无数次反射和吸收，几乎不可能从小孔逃逸。换言之，这个小孔就是一个人工制造的完美黑体。

Lummer 和 Pringsheim 在德国帝国物理技术研究所（Physikalisch-Technische Reichsanstalt）进行了长达数年的精密测量，他们测量了从 800K 到 2000K 温度范围内，不同波长处的辐射强度。这些实验数据清晰地表明了黑体辐射的两个基本特征：其一，在任何给定温度下，辐射强度随波长变化的曲线都呈现出一个单一的峰值；其二，随着温度升高，这个峰值的位置向短波方向（即高频方向）移动。这后一个规律后来被称为维恩位移定律（Wien's displacement law）。

维恩公式：经验规律的半理论推导

在实验数据积累的同时，理论物理学家也在尝试从热力学基本原理出发，推导黑体辐射的普适函数。其中最杰出的成就是 Wilhelm Wien 于1893年给出的位移定律的推导。Wien 的论证基于一个简洁而深刻的物理思想：多普勒效应。

Wien 设想，如果一个热力学系统整体以速度 $v$ 运动，那么其内部辐射场的光谱分布将会发生变化。考虑一个沿某个方向以速度 $v$ 运动的空腔，由于多普勒效应，腔内辐射的波长会发生收缩：波长 $\lambda$ 变为 $\lambda' = \lambda(1 - v/c)$（在非相对论近似下）。同时，根据相对论或经典力学的不同框架，系统的温度也会相应改变。Wien 证明，如果温度为 $T$ 的系统以速度 $v$ 运动，其“表观温度” $T'$ 满足 $T' = T(1 - v/c)$。

基于这一分析，Wien 推导出黑体辐射光谱分布必须满足的一般形式：

$$\rho_\lambda(\lambda, T) = \frac{1}{\lambda^5} \phi(\lambda T)$$

其中 $\rho_\lambda(\lambda, T)$ 是单位波长间隔内的辐射能量密度，$\phi$ 是一个仅依赖于乘积 $\lambda T$ 的普适函数。这一结果通常称为维恩位移定律的数学表述。它表明，如果我们知道某一温度下的辐射光谱，原则上可以推知所有温度下的光谱分布——它们之间通过一个简单的标度变换相联系。

维恩公式的提出

基于位移定律的形式以及部分实验数据，Wien 在1896年进一步提出了一个具体的经验公式：

$$\rho_\lambda(\lambda, T) = \frac{C_1}{\lambda^5} e^{-C_2/(\lambda T)}$$

其中 $C_1$ 和 $C_2$ 是两个经验常数。这个公式通常写作维恩公式（Wien's law）或维恩辐射公式。它在短波区域与实验数据吻合良好，特别是能够精确描述紫外线波段的辐射分布。

然而，Wien 公式的提出并非纯粹的经验拟合。Wien 本人曾尝试从更基本的物理原理出发论证这一公式，其中最具影响力的是他对熵与概率之间关系的探讨。Wien 假设黑体辐射的熵 $S$ 与辐射场的某种状态参数之间满足某种微分关系，并据此推导出上述指数形式。这一推导虽然不够严格，但在当时的物理学界产生了重要影响，并启发了后来普朗克的工作。

瑞利-金斯公式：经典物理的极限

正当维恩公式在短波区域大放异彩时，另一批物理学家从经典物理学的框架出发，试图推导出黑体辐射的全频谱分布。其中最著名的成果是瑞利勋爵（Lord Rayleigh，后晋封为瑞利男爵）和 James Jeans 在1900年前后独立完成的工作。

瑞利-金斯公式的推导

瑞利-金斯的推导基于经典统计力学和电磁理论的结合，其逻辑链条可以概括如下：

第一步：辐射场的自由度计数

在经典物理中，电磁场被视为连续的波动系统。一个有限体积 $V$ 中的电磁驻波必须满足边界条件（通常是波节或波腹出现在腔壁）。对于给定波长 $\lambda$ 的辐射，允许存在的振动模式数量为：

$$dN = \frac{8\pi V}{\lambda^4} d\lambda$$

这一结果来自对波矢空间中模式密度的积分。注意 $\lambda \propto 1/k$，所以 $dN \propto k^3 dk \propto d\lambda/\lambda^4$。

第二步：能量均分定理

根据经典统计力学的能量均分定理（equipartition theorem），每个自由度在热平衡时的平均能量为 $k_B T$（这里 $k_B$ 是玻尔兹曼常数）。对于电磁场，每个振动模式对应两个自由度（分别对应两个偏振方向），因此每个模式的平均能量为：

$$\langle \varepsilon \rangle = 2 \times \frac{1}{2}k_B T = k_B T$$

第三步：辐射能量密度

综合以上两步，单位波长间隔内的辐射能量密度为：

$$\rho_\lambda(\lambda, T) = \frac{8\pi V}{\lambda^4} \cdot \frac{k_B T}{V} = \frac{8\pi k_B T}{\lambda^4}$$

写成更常见的频率形式（利用 $\lambda = c/\nu$，$d\lambda = (c/\nu^2)d\nu$）：

$$\rho_\nu(\nu, T) = \frac{8\pi \nu^2 k_B T}{c^3}$$

这就是著名的瑞利-金斯公式（Rayleigh-Jeans law）。

公式的命运：紫外灾难

瑞利-金斯公式在长波区域与实验数据吻合良好，这是经典物理的合理预期。然而，当这一公式被外推到短波区域时，灾难性的后果出现了：按照这一公式，辐射能量密度与 $\nu^2$ 成正比，这意味着当频率趋于无穷大（即波长趋于零）时，辐射能量密度将无限发散。这一荒谬的结果被 Kelvin 勋爵形象地称为“紫外灾难”（ultraviolet catastrophe）。

紫外灾难的物理含义是：在经典物理的框架下，热平衡辐射场包含无限多的高频模式，而每个模式都从热库中分得 $k_B T$ 的能量。如果高频模式的数目没有上限（实际上在经典理论中确实没有上限），那么总能量必然发散。这一定量上的荒谬预示着经典物理的某个基本假设必然存在问题。

普朗克公式：内插法的胜利

1900年10月19日，在德国物理学会的一次会议上，Max Planck 提交了一篇简短的论文，标题为《论维恩辐射定律的改进》。在这篇论文中，Planck 提出了一个全新的黑体辐射公式：

$$\rho_\lambda(\lambda, T) = \frac{C_1}{\lambda^5} \frac{1}{e^{C_2/(\lambda T)} - 1}$$

或等价地，用频率形式表示为：

$$\rho_\nu(\nu, T) = \frac{8\pi h \nu^3}{c^3} \frac{1}{e^{h\nu/(k_B T)} - 1}$$

这就是著名的普朗克公式（Planck's law），其中 $h$ 和 $k_B$ 分别是后来以 Planck 和 Boltzmann 命名的自然常数。

内插法的数学技巧

值得注意的是，Planck 本人最初并未从第一性原理推导出这一公式。他的策略是一种数学内插法：已知维恩公式在短波极限下正确，而瑞利-金斯公式在长波极限下正确，Planck 试图找到一个函数，它在 $\lambda T \to 0$（短波）时趋近于维恩公式，在 $\lambda T \to \infty$（长波）时趋近于瑞利-金斯公式。

维恩公式可以写成：

$$\rho_\lambda^{\text{Wien}} = \frac{C_1}{\lambda^5} e^{-C_2/(\lambda T)}$$

而瑞利-金斯公式则为：

$$\rho_\lambda^{\text{RJ}} = \frac{8\pi k_B T}{\lambda^4}$$

注意到如果将指数函数 $e^{C_2/(\lambda T)}$ 在 $C_2/(\lambda T)$ 很小的条件下展开：

$$e^{C_2/(\lambda T)} \approx 1 + \frac{C_2}{\lambda T} + \cdots$$

那么：

$$e^{C_2/(\lambda T)} - 1 \approx \frac{C_2}{\lambda T} \quad (\text{当 } \lambda T \to \infty)$$

代入普朗克公式的长波极限：

$$\rho_\lambda \approx \frac{C_1}{\lambda^5} \frac{\lambda T}{C_2} = \frac{C_1}{C_2} \frac{T}{\lambda^4}$$

只要令 $\frac{C_1}{C_2} = 8\pi k_B$，即可得到瑞利-金斯公式的形式。

而在短波极限（$\lambda T \to 0$）下，$e^{C_2/(\lambda T)} - 1 \approx e^{C_2/(\lambda T)}$，普朗克公式自动约化为维恩公式。

常数的物理化

然而，Planck 的贡献远不止于找到一个内插公式。他进一步确定了两个经验常数 $C_1$ 和 $C_2$ 与基本自然常数之间的关系。Planck 证明：

$$C_1 = 2\pi h c^2, \quad C_2 = \frac{h c}{k_B}$$

因此普朗克公式可以写成完全由自然常数决定的形式：

$$\boxed{\rho_\nu(\nu, T) = \frac{8\pi h \nu^3}{c^3} \frac{1}{e^{h\nu/(k_B T)} - 1}}$$

这一公式将黑体辐射问题中出现的所有经验常数，归结为三个真正的基本自然常数：光速 $c$（来自相对论和电磁理论）、玻尔兹曼常数 $k_B$（来自热力学和统计力学）以及一个全新的常数 $h$。Planck 将 $h$ 称为作用量子（quantum of action），其量纲为能量乘以时间（即作用量的量纲），现代物理学中称之为普朗克常数。

玻尔兹曼熵公式与普朗克的统计论证

如果说1900年10月19日的论文解决了黑体辐射光谱分布的数学问题，那么仅仅六周之后的1900年12月14日，Planck 在德国物理学会的另一次报告中，给出了这一公式的物理论证。这一次，Planck 引入了后来被视为量子论诞生标志的概念——能量子（energy quantum）。

熵的统计基础

Planck 的论证起点是玻尔兹曼熵公式（Boltzmann entropy formula）：

$$S = k_B \ln W$$

其中 $W$ 是系统给定宏观状态所对应的微观状态数（或称热力学概率），$k_B$ 是玻尔兹曼常数。这一公式如今被视为连接统计力学与热力学的枢纽。然而，在1900年，玻尔兹曼的统计观念尚未被所有物理学家完全接受。Planck 本人虽然最终成为玻尔兹曼统计思想的坚定支持者，但他的这一转变恰恰发生在量子论诞生的关键时刻。

谐振子的能量离散化

Planck 论证的核心是线性谐振子（linear harmonic oscillator）模型。在经典电动力学中，带电粒子（如电子）在平衡位置附近的微小振动可以近似为线性谐振子，而热平衡时空腔内的辐射场正是通过这些振子与电磁场不断交换能量而维持平衡的。

Planck 假设腔壁上的谐振子只能具有分立的能量值，而不是经典物理所允许的连续取值。具体而言，频率为 $\nu$ 的谐振子，其能量只能是某个基本单元的整数倍：

$$E_n = n h \nu, \quad n = 0, 1, 2, 3, \cdots$$

这里 $h\nu$ 便是能量子——频率为 $\nu$ 的振子所能吸收或发射的最小能量单元。

熵的计算

现在考虑 $N$ 个频率相同（均为 $\nu$）、处于温度 $T$ 的热平衡谐振子。设这些振子的总能量为 $E_N$，那么每个振子的平均能量为 $\bar{E} = E_N/N$。

按照玻尔兹曼的统计方法，我们需要计算这 $N$ 个振子在满足总能量约束 $E_N = \sum_i E_i$ 的条件下，所有可能的能量分配方式（即微观状态数）$W$。

如果振子的能量可以连续变化，则 $W$ 依赖于 $E_i$ 的具体数值；但如果能量被限制为离散值 $0, h\nu, 2h\nu, \cdots$，则我们需要计算将总能量 $E_N$ 表示为 $P$ 个能量子（每个大小为 $h\nu$）分配到 $N$ 个振子上的方式数。

这一问题等价于：将 $P$ 个不可区分的能量子分配到 $N$ 个振子上，每个振子可以接受任意多个（包括零个）能量子。数学上，这相当于在 $N+P$ 个位置中插入 $N-1$ 个分隔符的方式数：

$$W = \frac{(N + P - 1)!}{P!(N-1)!}$$

其中 $P = E_N/(h\nu)$ 是能量子的总数。

吉布斯熵与普朗克熵

利用斯特林公式（Stirling's approximation）$n! \approx n^n e^{-n} \sqrt{2\pi n}$，我们可以计算 $S = k_B \ln W$：

$$\ln W \approx (N+P)\ln(N+P) - N\ln N - P\ln P$$

因此：

$$S = k_B \left[ (N+P)\ln(N+P) - N\ln N - P\ln P \right]$$

将 $P = E_N/(h\nu)$ 和 $\bar{E} = E_N/N$ 代入，并进行代数运算，最终得到：

$$S = k_B N \left[ \left(1 + \frac{\bar{E}}{h\nu}\right)\ln\left(1 + \frac{\bar{E}}{h\nu}\right) - \frac{\bar{E}}{h\nu}\ln\left(\frac{\bar{E}}{h\nu}\right) \right]$$

从热力学关系 $\frac{1}{T} = \frac{\partial S}{\partial E}$ 出发，可以推导出每个振子的平均能量：

$$\bar{E} = \frac{h\nu}{e^{h\nu/(k_B T)} - 1}$$

这正是普朗克分布中频率为 $\nu$ 的辐射模式所应具有的平均能量。将这一结果代入辐射场的模式密度 $8\pi \nu^2/c^3$，立即得到普朗克辐射公式。

能量子概念的历史意涵

普朗克在1900年12月14日提出的能量子假说，其革命性意义在于它打破了两个根深蒂固的经典物理信念：

连续性与可分性

经典物理认为，自然界中的一切物理量都可以连续变化。无论是位置、速度、能量还是其他任何可观测量，在经典力学的框架下都取值于连续的实数域。这一信念可以追溯到古希腊哲学家芝诺（Zeno）的悖论，但在牛顿力学建立之后，它获得了现代科学的支撑：经典力学的基本方程（如牛顿第二定律）是关于连续变量的微分方程，这些方程的解在时间上演化是连续变化的。

能量子的引入则断言：能量的变化不是连续的，而是以 $h\nu$ 为基本单元的离散跳跃。这意味着在谐振子与辐射场之间交换能量时，只能以“整个包”的形式进行，每次交换的能量必须是 $h\nu$ 的整数倍。这与经典电动力学中“带电粒子可以以任意小的速率连续辐射能量”的图像截然相