黑体辐射与紫外灾难解析
了解黑体光谱的形状、古典物理学为何预测紫外灾难，以及普朗克如何透过能量量子化解决这个问题。
加热一块金属时，随著温度升高，它发出的光会改变：先是暗红色，接著变成橙色，最后变成黄白色。根据简单的古典物理预测，似乎很容易得出一个结论。如果温度越高，物质中的带电粒子就振动得越剧烈，那么波长越短的电磁波或许就会携带越来越多的能量。若将这种推论推展到极端，紫外线的能量就应该占据主导地位。
但事实并非如此。真实的热辐射光谱会先上升至一个峰值，然后在短波长区域急剧下降。古典物理学却预测了相反的结果–发射的能量会无限制地增加。这个失败后来被称为紫外灾难，而要解决这个问题，物理学家就必须改变他们对能量特性的看法。
01黑体光谱实际上是什么样子？
想像一个中空的腔体，其中一面墙上有一个小孔。进入小孔的辐射会在腔体内部不断反射，几乎不可能再次逃出。因此，这个小孔的行为几乎就像一个理想的吸收体。
如果将腔体维持在固定温度，内部的辐射最终会达到热平衡。从小孔射出的辐射，其光谱主要由温度决定，而不是由制成腔体的材料决定。
这就是黑体背后有用的物理模型：一种理想物体，能吸收所有入射的电磁辐射，并发射与其温度相对应的特征热辐射光谱。
现在来作出预测。如果温度升高，光谱会有什么变化？
我们会观察到两件事：
- 物体整体发射的能量增加；以及
- 发射最强时所对应的波长会移向较短的波长。
光谱并不是一个狭窄的尖峰。高温物体会在广泛的波长范围内发射辐射。其峰值波长 \(\lambda_{\text{max}}\)，就是曲线达到最大值时所对应的波长。
这个峰值的移动可用维恩位移定律描述：
\[ \lambda_{\text{max}}T=b \]
其中：
- \(\lambda_{\text{max}}\) 是以公尺为单位的峰值波长；
- \(T\) 是以克耳文为单位的绝对温度；以及
- \(b=2.898\times10^{-3}\text{ m K}\) 是维恩位移常数。
因此，当 \(T\) 增加时，\(\lambda_{\text{max}}\) 就会减小。
这解释了我们熟悉的金属从红热变成白热的现象。它还告诉我们一件更重要的事：自然界存在某种机制，能防止短波长辐射无限制地增加。
计算示例：高温灯丝在哪个波长发射最强？
一条灯丝的温度为 \(3000\text{ K}\)。计算其黑体光谱达到峰值时的波长，并说明这个波长是否位于可见光范围内。
步骤 1
\[ \lambda_{\text{max}}T=b \]
步骤 2
\[ \lambda_{\text{max}}=\frac{b}{T} \]
步骤 3
\[ \lambda_{\text{max}} =\frac{2.898\times10^{-3}}{3000} =9.66\times10^{-7}\text{ m} \]
步骤 4
\[ 9.66\times10^{-7}\text{ m}=966\text{ nm} \]
峰值波长约为 \(966\text{ nm}\)，位于红外线范围内。
这并不代表灯丝完全不会发出可见光。由于光谱涵盖广泛的波长范围，仍有部分辐射落在可见光范围内。只是光谱峰值恰好位于红光波长之外。
02为什么古典物理学预测紫外线辐射会过多？
要了解古典物理学的预测在哪里出错，最好先不要再从发光物体的表面思考，而是回到腔体本身。
电磁波可以在腔体内形成驻波图样。只有某些特定的图样能够容纳在腔体内，这有点像吉他弦上允许形成的驻波，只不过腔体是三维的，而且波是电磁波。
每一种允许存在的驻波图样都称为一个模式。
在低频率下，可能存在的模式相对较少。在较高频率下，则有更多不同的驻波图样能够容纳在腔体内。
事实上，在一个很小的频率区间内，可用模式的数量大致随 \(f^2\) 增加，其中 \(f\) 是频率。
因此，高频率会带来一种危险的组合：
- 可用的模式数量更多；以及
- 古典物理学认为，每个模式平均应获得大致相同的热能。
第二个观念来自古典的能量均分定理。
对于处于热平衡的电磁模式，古典统计力学预测其平均能量为
\[ \overline{E}=kT \]
其中 \(k=1.381\times10^{-23}\text{ J K}^{-1}\) 是波兹曼常数。
想像学校餐厅有一条荒谬的规定：每一位到场的学生都会得到一份三明治。只要学生人数有限，这个规定就行得通。但现在假设，当你沿著队伍走向紫外线区域时，学生人数会不断增加，没有上限。那么三明治的总需求量就会变成无限大。
在这个比喻中，学生代表电磁模式，而每份三明治代表古典物理学所预测的能量 \(kT\)。
这个比喻的限制在于，能量并不是由某个人实际分发出去的，而电磁模式也不是一群互不相关、饥肠辘辘的青少年。这个比喻真正有用的地方在于数量的计算：古典物理学为每个可用模式分配非零的平均能量，而高频模式的数量却持续增加。
03瑞利－金斯定律的预测
将古典模式计数与能量均分定理结合，就能得到瑞利－金斯定律。以频率表示时，每单位频率的能量密度为
\[ u(f,T)=\frac{8\pi f^2kT}{c^3} \]
其中：
- \(u(f,T)\) 描述辐射能量如何随频率分布；
- \(f\) 是以赫兹为单位的频率；
- \(T\) 是以克耳文为单位的温度；
- \(k\) 是波兹曼常数；以及
- \(c\) 是光速。
请注意其中最重要的部分：
\[ u(f,T)\propto f^2 \]
在固定温度下，频率加倍，古典物理学预测的能量密度就会变成原来的四倍。
如果频率变成三倍，能量密度就会变成九倍。
当 \(f\) 趋近无限大时，预测的能量密度也会趋近无限大。
由于紫外线的频率比可见光高，这种在高频端失控增加的预测便被称为紫外灾难。
这场「灾难」是理论上的灾难，而不是实际物体遭遇的灾难。高温烤箱并没有偷偷向你的厨房发射无限量的紫外线。
04为什么这个预测看起来合情合理？
古典模型并不愚蠢。在频率足够低时，它的表现相当良好。
这一点很重要。一个好的替代理论不应该直接抛弃已经成功的古典物理学，而应该解释旧理论为什么在某个范围内有效，却在另一个范围内失效。
从实验结果来看，瑞利－金斯定律在长波长范围内，对黑体辐射的预测相当合理。当波长逐渐缩短时，才会出现严重的偏差。
因此，问题可以更精确地表述为：
为什么高频模式获得的热能会比古典能量均分定理所预测的更少？
05普朗克的改变：能量以与频率相关的离散份量交换
马克斯・普朗克引入了一项新的限制。与频率为 \(f\) 的模式相关的能量交换，是以离散的份量进行，其能量与下式有关：
\[ E=hf \]
其中 \(h=6.626\times10^{-34}\text{ J s}\) 是普朗克常数。
更精确地说，在普朗克的模型中，允许的振荡子能量可以写成
\[ E_n=nhf \]
其中 \(n=0,1,2,\ldots\)。
关键在于每一个能量阶距的大小：
\[ \Delta E=hf \]
低频模式的能量阶距较小。高频模式的能量阶距则较大。
想像一间游乐场，里面的不同游戏需要不同数量的代币才能开始玩。古典能量均分定理实际上假设每种游戏都可以吸收任意微小的能量。普朗克的规则则表示，高频游戏所需的最低代币数量要大得多。
在固定温度下，可用的热能有限。因此，高频模式更难被激发。
这个比喻的限制在于，模式并不是真的在累积热能「代币」，但它抓住了重要的竞争关系：量子能量尺度 \(hf\) 正在与热能尺度 \(kT\) 比较。
这种比较是理解黑体辐射最有用的方法之一：
\[ \frac{hf}{kT} \]
如果 \(hf\ll kT\)，能量阶距相对于可用热能而言很小，因此古典物理学的行为是一个良好的近似。
如果 \(hf\) 变得与 \(kT\) 相当，或大于 \(kT\)，离散的能量阶距就会产生显著影响。
06量子化如何消除紫外灾难
普朗克的模型给出腔体模式的平均能量：
\[ \overline{E} = \frac{hf}{e^{hf/(kT)}-1} \]
这个公式看起来比古典结果 \(kT\) 复杂，但它所表达的物理意义，其实比公式乍看之下更清楚。
在低频率下，\(hf\) 相对于 \(kT\) 很小，因此这个表示式会趋近古典值 \(kT\)。
在非常高的频率下，\(hf/(kT)\) 会变得很大。分母中的指数项会极快地增长，因此每个高频模式的平均能量会逐渐趋近零。
高频模式的数量仍然很多。
关键的改变在于，这些模式不再各自被分配 \(kT\) 的能量。
将普朗克的模式平均能量与相同的腔体模式计数结合，就会得到
\[ u(f,T) = \frac{8\pi hf^3}{c^3} \frac{1}{e^{hf/(kT)}-1} \]
起初，随著 \(f\) 增加，可用模式的数量也增加，因此光谱会上升。
然而，当 \(hf\) 相对于 \(kT\) 变得太大时，激发模式所需的能量就会过高。指数抑制效应开始占上风，光谱便会下降。
如此一来，光谱就会形成有限的峰值，而不是在紫外线区域无限制地暴增。
计算示例：古典近似在什么时候失效？
考虑温度为 \(5800\text{ K}\) 的热辐射，这大致是太阳可见表面的温度量级。比较波长为 \(10\,\mu\text{m}\) 和 \(200\text{ nm}\) 时的量子修正。
在相同波长下，普朗克预测与相应的古典瑞利－金斯预测之间的比值可写成
\[ \frac{\text{Planck}}{\text{classical}} = \frac{x}{e^x-1} \]
其中
\[ x=\frac{hc}{\lambda kT} \]
步骤 1
这里，
\[ \lambda=10\,\mu\text{m}=1.0\times10^{-5}\text{ m} \]
因此
\[ x= \frac{(6.626\times10^{-34})(3.00\times10^8)} {(1.0\times10^{-5})(1.381\times10^{-23})(5800)} \approx0.248 \]
步骤 2
\[ \frac{x}{e^x-1} = \frac{0.248}{e^{0.248}-1} \approx0.88 \]
普朗克的结果约为古典预测的 \(88\%\)。在这个较长的波长下，两种模型的结果相当接近。
步骤 3
\[ \lambda=200\text{ nm}=2.00\times10^{-7}\text{ m} \]
\[ x= \frac{(6.626\times10^{-34})(3.00\times10^8)} {(2.00\times10^{-7})(1.381\times10^{-23})(5800)} \approx12.4 \]
步骤 4
\[ \frac{x}{e^x-1} = \frac{12.4}{e^{12.4}-1} \approx5.1\times10^{-5} \]
在 \(200\text{ nm}\) 时，普朗克的预测只有古典结果的约 \(0.005\%\)。
这就是以数值呈现的紫外灾难。古典预测在短波长下并不是只有轻微的误差，而是最终会变得极度不准确。
07最容易令人误解的观念
一种常见的解读是：普朗克透过设定一个最高允许频率来解决问题。
事实并非如此。
不存在一道硬性的界线，使辐射在某个频率突然变得不可能。黑体光谱在高频端有一条平滑的尾部。
量子模型的做法，是让高频激发变得越来越不可能，因为所需的能量阶距 \(hf\) 相对于 \(kT\) 会越来越大。
这个差异很重要。任意设定一个紫外线截止频率，或许能避免总能量变成无限大，但这种做法无法自然地重现观察到的平滑光谱，以及光谱随温度系统性变化的现象。
08温度升高时，实际上改变了什么？
假设温度升高。
很容易让人说：「光谱会向左移。」
但这种说法并不完整。
温度较高的黑体会在广泛的波长范围内发出更强的辐射。同时，光谱峰值会移向较短的波长。
这个原因同样可以从 \(hf\) 与 \(kT\) 之间的竞争看出。提高 \(T\) 会增加热能尺度 \(kT\)。原本相对于可用热能而言代价太高的频率，现在就变得比较容易被激发。
因此，光谱中的高频部分会有更多模式被激发，而特征峰值则会移向较高的频率，亦即较短的波长。
维恩位移定律
\[ \lambda_{\text{max}}T=b \]
就是这个现象可测量的结果。
09古典预测与量子预测的并列比较
| 特征 | 古典瑞利－金斯模型 | 普朗克模型 | 观察结果 |
|---|---|---|---|
| 长波长 | 表现良好 | 趋近古典结果 | 两者都符合观察 |
| 高频率 | 能量持续增加 | 高频模式受到强烈抑制 | 光谱下降 |
| 总发射能量 | 若包含所有频率，则发散至无限大 | 有限 | 有限 |
| 光谱是否有峰值 | 高频端不存在有限的峰值 | 有 | 有 |
| 能量交换 | 实际上是连续的 | 以与 \(hf\) 相关的离散份量进行量子化 | 量子模型符合光谱 |
| 温度升高的影响 | 每个模式的古典能量增加，但紫外灾难仍然存在 | 更高的频率变得较容易被激发 | 峰值移向较短波长 |
重要的教训并不只是「古典物理学错了」。古典物理学其实是更完整的量子描述在低频率下的极限情况。
10问题与解答
问题 1
一个黑体的温度为 \(4000\text{ K}\)。
计算它的峰值波长。第二个黑体的温度为 \(8000\text{ K}\)。不必重复完整的计算，直接写出第二个物体的峰值波长，并解释这个变化。
使用 \(b=2.898\times10^{-3}\text{ m K}\)。
解答 1
温度为 \(4000\text{ K}\) 的物体，其峰值波长约为 \(725\text{ nm}\)；而温度为 \(8000\text{ K}\) 的物体，其峰值波长约为 \(362\text{ nm}\)。
对第一个物体，
\[ \lambda_{\text{max}} = \frac{b}{T} = \frac{2.898\times10^{-3}}{4000} = 7.245\times10^{-7}\text{ m} \]
因此，
\[ \lambda_{\text{max}}\approx725\text{ nm} \]
对第二个物体，温度加倍。由于
\[ \lambda_{\text{max}}\propto\frac{1}{T} \]
当 \(T\) 加倍时，峰值波长就会减半：
\[ \lambda_{\text{max}} = \frac{725}{2} \approx362\text{ nm} \]
容易犯的错误是说，较热物体的整个光谱只是移到了 \(362\text{ nm}\)。事实并非如此。光谱仍然涵盖广泛的波长范围，其强度会增加，而峰值则移向较短的波长。
问题 2
根据瑞利－金斯定律，
\[ u(f,T)=\frac{8\pi f^2kT}{c^3} \]
同一个腔体内的两个频率分别为 \(3.0\times10^{14}\text{ Hz}\) 和 \(9.0\times10^{14}\text{ Hz}\)。
在相同温度下：
- 计算这两个频率所对应的古典光谱能量密度之比。
- 解释为什么将这种推论延伸至任意高的频率，会导致紫外灾难。
解答 2
瑞利－金斯模型预测，在 \(9.0\times10^{14}\text{ Hz}\) 时的光谱能量密度，是在 \(3.0\times10^{14}\text{ Hz}\) 时的九倍；而且，随著频率持续按 \(f^2\) 增长，将所有频率都纳入计算时，总能量就会变成无限大。
在温度固定时，除了 \(f^2\) 以外的所有因子都会相消：
\[ \frac{u_2}{u_1} = \frac{f_2^2}{f_1^2} \]
代入数值，
\[ \frac{u_2}{u_1} = \left( \frac{9.0\times10^{14}} {3.0\times10^{14}} \right)^2 = 3^2 = 9 \]
因此，
\[ u_2=9u_1 \]
这个计算本身还不会造成灾难。问题在于，如果不设上限，就一直套用同一条规则。
当 \(f\) 增加时，
\[ u(f,T)\propto f^2 \]
因此，预测的能量密度会持续增加，而不是先达到峰值，再逐渐下降。
一个容易让人走错的思路是：「每个个别频率仍然有有限的数值，所以没有问题。」真正改变结论的关键条件是，黑体光谱包含连续的频率范围。将越来越高频率所对应的预测贡献全部加起来，就无法得到有限的总能量。
然而，真实光谱在高频端会急剧下降。
问题 3
一个腔体的温度为 \(6000\text{ K}\)。使用下式比较 \(500\text{ nm}\) 和 \(200\text{ nm}\) 的辐射：
\[ x=\frac{hf}{kT}=\frac{hc}{\lambda kT} \]
取 \(h=6.626\times10^{-34}\text{ J s}\)、\(c=3.00\times10^8\text{ m s}^{-1}\)，以及 \(k=1.381\times10^{-23}\text{ J K}^{-1}\)。
- 计算每个波长所对应的 \(x\)。
- 指出哪个波长受到量子效应的抑制较强。
- 解释为什么在这种情况下，古典能量均分定理会变得特别不准确。
解答 3
计算结果约为：在 \(500\text{ nm}\) 时，\(x=4.80\)；在 \(200\text{ nm}\) 时，\(x=12.0\)。因此，\(200\text{ nm}\) 的辐射受到的抑制强得多。
对于 \(500\text{ nm}\)，
\[ \lambda=5.00\times10^{-7}\text{ m} \]
而且
\[ x= \frac{(6.626\times10^{-34})(3.00\times10^8)} {(5.00\times10^{-7})(1.381\times10^{-23})(6000)} \approx4.80 \]
对于 \(200\text{ nm}\)，
\[ \lambda=2.00\times10^{-7}\text{ m} \]
因此
\[ x= \frac{(6.626\times10^{-34})(3.00\times10^8)} {(2.00\times10^{-7})(1.381\times10^{-23})(6000)} \approx12.0 \]
在 \(200\text{ nm}\) 时，量子能量尺度 \(hf\) 约为热能尺度 \(kT\) 的十二倍。因此，激发这种模式需要相对较大的能量阶距。
容易让人误信的古典论证是：不论频率是多少，每个模式都应该获得平均能量 \(kT\)。这个假设忽略了最低量子能量阶距会随频率增加的事实。
由于在 \(200\text{ nm}\) 时 \(hf\gg kT\)，高频激发不太可能发生。因此，普朗克的表示式所分配给这些模式的平均能量，会远低于古典能量均分定理的预测。
问题 4
一名学生将黑体光谱绘制成以波长为横轴的图，发现图形在 \(500\text{ nm}\) 处达到峰值。
学生立即声称，同一辐射以频率为横轴绘图时，峰值必定出现在
\[ f=\frac{c}{500\text{ nm}} \]
因为波长与频率满足 \(c=f\lambda\)。
这个结论一定正确吗？
使用以下事实：对应的小区间内所包含的能量必须满足
\[ B_\nu\,d\nu=B_\lambda\,|d\lambda| \]
而且
\[ \left|\frac{d\lambda}{d\nu}\right|=\frac{c}{\nu^2} \]
其中 \(B_\lambda\) 是每单位波长的光谱辐亮度，而 \(B_\nu\) 是每单位频率的光谱辐亮度。
解答 4
不一定。频率光谱的峰值通常不会出现在波长光谱峰值的简单倒数所对应的位置。
关系式
\[ c=f\lambda \]
确实可以将某个特定波长对应到某个特定频率。问题在于，假设改变横轴变数后，纵轴所代表的量仍然保持不变。
光谱描述的是每个区间内的能量。相同大小的波长区间，并不对应相同大小的频率区间。
根据
\[ B_\nu\,d\nu=B_\lambda\,|d\lambda| \]
可得
\[ B_\nu = B_\lambda \left|\frac{d\lambda}{d\nu}\right| \]
因此
\[ B_\nu = B_\lambda\frac{c}{\nu^2} \]
等价地，由于 \(\lambda=c/\nu\)，
\[ B_\nu=B_\lambda\frac{\lambda^2}{c} \]
这个乘数并不是常数，而是会随波长或频率改变。
这会改变图形的形状，因而也会改变其最大值的位置。
容易让人走错的思路，是将光谱当成一组只需重新标示横轴就能直接转换的点。真正的问题在于，每个绘图高度所代表的是「每单位横轴区间的能量」。从波长区间改成频率区间时，区间大小也随之改变。
这并不会改变实际的辐射，只会改变表示同一辐射的方式。
问题 5
假设你不知道普朗克的精确公式。你只知道，一个成功的理论必须在频率足够高时，使每个腔体模式的平均能量低于古典值。
一位理论学家提出以下一般形式：
\[ \overline{E} = kT\, g\left(\frac{hf}{kT}\right) \]
其中 \(g(x)\) 具有以下性质：
- 当 \(x\ll1\) 时，\(g(x)\rightarrow1\)；
- 当 \(x\gg1\) 时，\(g(x)\) 会迅速下降并趋近零。
每单位频率的电磁模式数量仍然与 \(f^2\) 成正比。
即使不知道 \(g\) 的确切数学形式：
- 解释为什么这个模型在低频率下会恢复瑞利－金斯定律的结果。
- 解释这个模型如何避免紫外灾难。
- 证明光谱的特征频率必须与 \(T\) 成正比。
- 因此解释为什么特征波长必须与 \(1/T\) 成正比。
- 不使用详细积分，而是根据光谱的形状，预测总辐射能量密度如何随温度变化。
解答 5
这个模型会自动在低频率下呈现古典行为，抑制足够高的频率，预测 \(f_{\text{characteristic}}\propto T\)，预测 \(\lambda_{\text{characteristic}}\propto1/T\)，并得出总辐射能量与 \(T^4\) 成正比的结果。
从提出的模式平均能量开始：
\[ \overline{E} = kT\, g\left(\frac{hf}{kT}\right) \]
在低频率下，
\[ \frac{hf}{kT}\ll1 \]
因此，根据题目给出的规则，
\[ g\left(\frac{hf}{kT}\right)\rightarrow1 \]
所以
\[ \overline{E}\rightarrow kT \]
这正是古典能量均分定理的结果。由于模式密度仍然与 \(f^2\) 成正比，
\[ u(f,T)\propto f^2kT \]
因此可以恢复瑞利－金斯定律的行为。
在高频率下，
\[ \frac{hf}{kT}\gg1 \]
而 \(g\) 会迅速下降并趋近零。尽管模式数量仍然随 \(f^2\) 增加，但每个模式的平均能量可以下降得足够快，使高频尾部逐渐减小，而不是发散至无限大。
这就是关键的竞争关系。即使模式数量更多，只要每个模式的能量下降得足够快，就不一定会出现紫外灾难。
现在定义
\[ x=\frac{hf}{kT} \]
光谱中的某个特征，例如峰值，只要函数 \(g\) 本身不随温度改变，就会出现在某个特定的 \(x\) 值，假设这个值为 \(x=x_0\)。
那么
\[ x_0=\frac{hf_{\text{peak}}}{kT} \]
因此
\[ f_{\text{peak}} = \frac{x_0k}{h}T \]
乘在 \(T\) 前面的所有量都是常数。因此，
\[ f_{\text{peak}}\propto T \]
由于波长与频率的关系为
\[ \lambda=\frac{c}{f} \]
可得
\[ \lambda_{\text{peak}}\propto\frac{1}{T} \]
这就是维恩位移定律背后的基本结构。
同样的推理还隐藏著另一个结果。光谱的形式为
\[ u(f,T)\propto f^2kT\, g\left(\frac{hf}{kT}\right) \]
对于缩放变数 \(x\) 的固定值，频率与 \(T\) 成正比。因此，\(f^2kT\) 这一部分所决定的纵向尺度，会与下式成正比：
\[ T^2\times T=T^3 \]
同时，由于特征频率与 \(T\) 成正比，在频率图上，光谱的水平宽度也会随 \(T\) 成正比。
因此，光谱下方的面积会与下式成正比：
\[ T^3\times T=T^4 \]
所以，总辐射能量密度与
\[ T^4 \]
成正比。
一个容易令人误解的答案是，普朗克的解决方法只是「截断紫外线」。这个一般性论证揭示了更深层的道理。成功的量子修正会透过 \(hf/(kT)\) 这个比值，将频率与温度联系起来。单是这个比较，就能解释为什么古典物理学会在低频率下重新成立、为什么紫外线发散会消失、为什么光谱峰值会随温度移动，以及为什么物体越热，总发射能量就会增加得越快。
同样的概念–将 \(hf\) 这类量子能量尺度与另一个物理上可用的能量尺度进行比较–在研究光子与光电效应时也会再次派上用场。