4、普朗克公式的推导（二）：从平均能量到辐射能量密度、普朗克公式的最终形式
好，咱们接着往下走。上一章我们搞定了黑体辐射的「微观模型」——谐振子，也找到了它的平均能量表达式。但说实话，那会儿的平均能量还是个「光杆司令」，跟实验测到的辐射能量密度还差着一层窗户纸。
这一章，我就带你亲手把这层纸捅破。说白了，就是要把单个谐振子的平均能量，乘以「空腔里到底有多少个谐振子」，再换算成单位体积、单位频率里的能量分布。最终，那个让整个物理学界为之颤抖的普朗克公式，就会完整地摆在你面前。
4.1 从「一个振子」到「一群振子」：模式密度的概念
我个人习惯，遇到这种「从微观到宏观」的跳跃，先问自己一个问题：空腔里到底有多少个谐振子在振动？
你想想看，黑体空腔是一个封闭的金属盒子，里面的电磁波是驻波。每一个驻波模式，就对应一个谐振子。所以，我们需要的其实是——单位体积内，频率在 ν 到 ν+dν 之间的驻波模式有多少个。这个量，物理上叫 模式密度 g(ν)。
这个推导，经典电动力学已经帮我们算好了。我记得当年学的时候，老师花了整整两节课推这个，其实核心就三步：
- 驻波条件：在边长为 L 的立方体空腔里，波矢 k 的三个分量必须满足 kx = nxπ/L，nx 是正整数。这就像琴弦两端固定，波长必须刚好「卡」在腔体里。
- k 空间计数：在 k 空间（波矢空间）里，每个模式占据的体积是 (π/L)³。那么，在半径为 k、厚度为 dk 的球壳里，模式数就是球壳体积除以单个模式体积，再乘以 1/8（因为 nx, ny, nz 都取正数）。
- 考虑偏振：电磁波是横波，每个 k 对应两个独立的偏振方向。所以模式数要乘以 2。
把 k = 2πν/c 代进去，一通整理，你就得到：
模式密度公式（单位体积，单位频率间隔）
g(ν) = 8πν² / c³
这个结果，经典物理和量子物理都认。它纯粹是几何和波动学的结论，跟能量是连续还是离散没关系。所以，我们直接拿来用。
一个小经验：我刚开始做计算的时候，总忘记乘那个「2」（偏振）。结果算出来的能量密度正好少一半，跟实验数据对不上，查了半天才发现。嗯，偏振这个因子，你最好在公式旁边标个记号。
4.2 经典物理的「死胡同」：瑞利-金斯公式
有了模式密度，经典物理是怎么走的呢？
经典统计力学告诉我们，每个谐振子的平均能量是 kT（能量均分定理）。那么，辐射能量密度 u(ν,T) 就是：
u(ν,T) = g(ν) × 平均能量 = (8πν² / c³) × kT
这就是著名的 瑞利-金斯公式。
你一看就知道问题在哪——u(ν,T) 随 ν² 增长，频率越高，能量密度越大，到紫外区直接爆炸。这就是物理学史上著名的「紫外灾难」。我当年第一次看到这个结果，第一反应是：「这公式是不是写错了？」后来发现，没写错，经典物理在这个问题上，就是走到了死胡同。
| 公式名称 | 平均能量 | 能量密度 u(ν,T) | 高频行为 | 结论 |
|---|---|---|---|---|
| 瑞利-金斯 | kT | (8πν²/c³) · kT | ν² → ∞ | 紫外灾难 |
| 普朗克（最终） | hν/(ehν/kT-1) | (8πν²/c³) · hν/(ehν/kT-1) | e-hν/kT → 0 | 与实验完美吻合 |
你看，差别就在平均能量上。经典物理给了 kT，普朗克给了 hν/(ehν/kT-1)。一个简单的替换，结果天差地别。
4.3 普朗克的「神来之笔」：量子化平均能量
上一章我们推导了，如果谐振子的能量是离散的（ε = nhν），那么它的平均能量就是：
普朗克平均能量
⟨E⟩ = hν / (ehν/kT - 1)
这个式子，你仔细看看它的极限行为：
- 低频极限（hν ≪ kT）：ehν/kT ≈ 1 + hν/kT，代入得 ⟨E⟩ ≈ kT。跟经典结果一致。这说明在低频区，量子效应不明显，经典物理「蒙对了」。
- 高频极限（hν ≫ kT）：ehν/kT 巨大，⟨E⟩ ≈ hν·e-hν/kT，指数衰减。这就把高频能量「压」住了，紫外灾难自然消失。
我曾经在给学生讲课时打过一个比方：经典平均能量 kT 像是一个「无底洞」，频率越高，能量往里灌得越多；而普朗克平均能量像是一个「有底的水桶」，每个能量子 hν 就是桶的容量，频率越高，桶越小，能装的水（能量）反而越少。这个图像，你想想看，是不是很直观？
4.4 最终形式：普朗克公式诞生
现在，我们把两样东西拼在一起：
- 模式密度：g(ν) = 8πν² / c³
- 平均能量：⟨E⟩ = hν / (ehν/kT - 1)
辐射能量密度 u(ν,T) 就是它们的乘积：
普朗克黑体辐射公式（最终形式）
u(ν,T) = (8πν² / c³) · [hν / (ehν/kT - 1)]
或者写成更紧凑的形式：
u(ν,T) = (8πhν³ / c³) · [1 / (ehν/kT - 1)]
这个公式，就是整个黑体辐射问题的终极答案。它给出了任意温度 T 下，黑体空腔内在频率 ν 处的辐射能量密度。1900 年 12 月 14 日，普朗克在德国物理学会上公布了它，这一天后来被公认为量子力学的诞生日。
注意单位：u(ν,T) 的单位是 J·s/m³（焦耳·秒每立方米），它表示「单位体积内、单位频率间隔里的辐射能量」。如果你看到公式里写的是 u(λ,T)（按波长表示），形式会不一样，因为 dν = -c/λ² dλ，代入后要小心符号和系数。我建议你始终用频率形式作为「基准」，需要波长形式时再转换。
4.5 知识结构总览
下面这张图，把整个推导的逻辑链条串起来了。你可以把它当作本章的「思维导图」：
4.6 几个容易踩的坑
推导过程看似简单，但有几个细节，我建议你特别留意：
- 模式密度里的「8π」是怎么来的？ 很多人只记公式，不记推导。我建议你至少手推一遍：从驻波条件 → k 空间球壳体积 → 乘以 1/8 → 乘以 2（偏振）→ 换成频率。走一遍，你就再也不会忘了那个 8π 是 8 × (π)，而不是 (8π) 的什么神秘系数。
- 平均能量公式的分母是 ehν/kT - 1，不是 ehν/kT + 1。玻色-爱因斯坦统计才是 -1，费米-狄拉克统计是 +1。光子是玻色子，所以这里是减号。我曾经见过有人把符号搞反，结果算出来的能量密度在低温下变成负的——这物理上显然不可能。
- 公式里的 h 是普朗克常数，不是约化普朗克常数 ħ。虽然 ħ = h/2π 在量子力学里更常用，但普朗克公式里用的是原始的 h。别搞混了。
一个实用技巧：如果你需要快速验证自己的推导是否正确，可以用 T = 6000 K（太阳表面温度）代入，看看峰值波长是不是在 500 nm 左右（可见光范围）。如果峰值跑到 X 射线区了，那你的公式肯定有问题。这个「物理直觉检验」比任何数学检查都管用。
好，到这一步，普朗克公式的完整推导就结束了。从模式密度到平均能量，再到最终的乘积形式，每一步都有清晰的物理图像支撑。你拿着这个公式，就能计算任何温度下黑体的辐射光谱——这，就是量子力学送给人类的第一份礼物。