经典连续理论为什么“不行”
经典物理用瑞利—金斯公式描述黑体辐射。它基于两个看似无懈可击的假设：
- 腔内有无限多的电磁振动模式，频率越高，模式数量越多（按频率的平方增长）。
- 能量均分定理：在热平衡下，每个模式的平均能量都等于 \(kT\)（\(k\) 是玻尔兹曼常数，\(T\) 是温度），与频率无关。
两者相乘，高频区的模式数量无限多，每个又都分到固定的 \(kT\)，总能量密度在频率趋于无穷时发散，这就是“紫外灾难”。经典框架下，能量被视为可以无限细分的连续量，没有任何机制阻止高频模式吸收能量。
离散化后的物理机制与统计规律
普朗克假设，频率为 \(\nu\) 的振子，其能量只能取 \(0, h\nu, 2h\nu, 3h\nu, \dots\) 这些分立值，不能连续变化。这一改动在统计上的后果是决定性的：
一个模式要被激发，至少需要一份能量 \(h\nu\)。 在温度为 \(T\) 的热平衡中，热运动能提供的典型能量尺度约为 \(kT\)。于是存在一个竞争关系：
- 低频模式（\(h\nu \ll kT\)）：一份能量很小，热运动轻易就能“凑出”好几份，模式被正常激发，平均能量接近经典值 \(kT\)。这保证了与实验低频端吻合。
- 高频模式（\(h\nu \gg kT\)）：一份能量远大于热运动能提供的典型值，该模式“买不起”哪怕最便宜的一份能量，因此几乎完全不被激发，对辐射的贡献趋近于零。
用统计力学的语言说，离散化后，一个振子的平均能量不再是固定的 \(kT\)，而是变成： \[ \bar{E} = \frac{h\nu}{e^{h\nu/kT}-1} \] 当 \(h\nu\) 很大时，分母中的指数项急剧增大，平均能量被指数式压制到近乎为零。高频模式数量再多，乘以一个趋近于零的平均能量，总贡献就变得有限。紫外灾难由此被“截断”。
所以，离散化并不直接改变电磁场的模式数量，它改变的是能量分配的统计规则：让高频模式因为“单价太高”而无法参与热平衡的能量分配。普朗克本人起初只把这当作一个数学上的“绝望之举”，但后来的发展表明，这个假设触及了能量交换过程本身的不连续性。
1. 能级与玻尔兹曼因子
频率为 \(\nu\) 的振子，能量只能取 \[ E_n = n h\nu,\qquad n=0,1,2,\dots \] 在温度 \(T\) 下，振子处于能级 \(n\) 的概率正比于玻尔兹曼因子 \(e^{-E_n/kT}\)。定义 \[ x \equiv \frac{h\nu}{kT} \] 则第 \(n\) 个能级的相对概率为 \(e^{-nx}\)。
2. 配分函数
归一化常数（单振子配分函数）为 \[ Z = \sum_{n=0}^{\infty} e^{-nx} = \frac{1}{1-e^{-x}} \] 这里用了等比级数求和，要求 \(x>0\)（即 \(T>0\)）。
3. 平均能量的定义与计算
平均能量按概率加权： \[ \bar{E} = \frac{\sum_{n=0}^{\infty} E_n e^{-E_n/kT}}{\sum_{n=0}^{\infty} e^{-E_n/kT}} = \frac{\sum_{n=0}^{\infty} n h\nu, e^{-nx}}{\sum_{n=0}^{\infty} e^{-nx}} \]
分子中的求和可以用配分函数对 \(x\) 求导得到： \[ \sum_{n=0}^{\infty} n, e^{-nx} = -\frac{d}{dx}\sum_{n=0}^{\infty} e^{-nx} = -\frac{d}{dx}\frac{1}{1-e^{-x}} = \frac{e^{-x}}{(1-e^{-x})^2} \]
于是 \[ \bar{E} = h\nu \cdot \frac{e^{-x}/(1-e^{-x})^2}{1/(1-e^{-x})} = h\nu \cdot \frac{e^{-x}}{1-e^{-x}} = \frac{h\nu}{e^{x}-1} \] 代回 \(x=h\nu/kT\)，即得 \[ \boxed{\bar{E} = \frac{h\nu}{e^{h\nu/kT}-1}} \]
4. 两个极限，看出离散化的作用
低频极限 \(h\nu \ll kT\)，即 \(x\ll 1\)。展开 \(e^x-1\approx x\)，得 \[ \bar{E} \approx \frac{h\nu}{h\nu/kT} = kT \] 回到经典能量均分结果。低频模式“买得起”能量，行为连续。
高频极限 \(h\nu \gg kT\)，即 \(x\gg 1\)。分母指数增长，得 \[ \bar{E} \approx h\nu, e^{-h\nu/kT} \to 0 \] 高频模式被指数式压低，几乎不被激发。
5. 与连续情形对比
如果能量连续，求和变成积分，配分函数为 \[ Z_{\text{连续}} = \int_0^{\infty} e^{-E/kT},dE = kT \] 平均能量为 \[ \bar{E}_{\text{连续}} = \frac{\int_0^{\infty} E e^{-E/kT}dE}{\int_0^{\infty} e^{-E/kT}dE} = kT \] 与频率无关。这正是紫外灾难的统计根源：每个模式不论频率高低都固定分到 \(kT\)。
所以，平均能量公式的推导本身只是标准统计力学，关键物理输入只有一条：能级是离散的、最小间隔为 \(h\nu\)。正是这个最小间隔与 \(kT\) 的比较，决定了该模式能否被热激发。
计算
能量 \(E\) 连续取值 \(0\) 到 \(\infty\)，概率密度正比于玻尔兹曼因子 \(e^{-E/kT}\)。配分函数为 \[ Z_{\text{连续}}=\int_0^{\infty} e^{-E/kT},dE = kT \]
平均能量： \[ \bar{E}_{\text{连续}}=\frac{\int_0^{\infty} E,e^{-E/kT},dE}{\int_0^{\infty} e^{-E/kT},dE} \]
分子用分部积分： \[ \int_0^{\infty} E,e^{-E/kT},dE = (kT)^2 \]
所以 \[ \boxed{\bar{E}_{\text{连续}} = \frac{(kT)^2}{kT} = kT} \]
关键对比
| 情形 | 平均能量 | 与频率的关系 |
|---|---|---|
| 连续能量 | \(kT\) | 无关，每个模式都分到 \(kT\) |
| 离散能量 | \(\dfrac{h\nu}{e^{h\nu/kT}-1}\) | 高频时指数式趋近于零 |
连续情形下，无论频率多高，每个模式都固定获得 \(kT\) 的平均能量。而高频模式数量按 \(\nu^2\) 无限增长，两者相乘导致总能量密度发散——这就是紫外灾难。离散化后，高频模式因 \(h\nu \gg kT\) 而几乎不被激发，平均能量被指数压制，发散才被截断。
\[ x \equiv \frac{h\nu}{kT} \]
则离散平均能量为
\[ \bar{E}_{\text{离散}} = \frac{h\nu}{e^{x}-1} = kT\cdot\frac{x}{e^{x}-1} \]
连续平均能量为 \(kT\)。两者之差为
\[ \Delta \equiv \bar{E}_{\text{连续}} - \bar{E}_{\text{离散}} = kT - \frac{h\nu}{e^{x}-1} = kT\left(1-\frac{x}{e^{x}-1}\right) \]
1. 精确表达式
\[ \boxed{\Delta = kT\left(1-\frac{x}{e^{x}-1}\right) = kT\cdot\frac{e^{x}-1-x}{e^{x}-1}} \]
由于 \(e^{x}-1-x>0\)（\(x>0\)），恒有 \(\Delta>0\)，即离散平均能量总是小于连续情形的 \(kT\)。离散化“扣掉”了一部分能量，扣得多少取决于 \(x=h\nu/kT\)。
2. 小 \(x\) 展开（低频，\(h\nu \ll kT\)）
先展开 \(\dfrac{x}{e^{x}-1}\)。用 \(e^{x}-1=x+\frac{x^2}{2}+\frac{x^3}{6}+\cdots\)，得
\[ \frac{x}{e^{x}-1} = \frac{1}{1+\frac{x}{2}+\frac{x^2}{6}+\cdots} = 1-\frac{x}{2}+\frac{x^2}{12}-\frac{x^4}{720}+\cdots \]
（这是伯努利数展开，注意没有 \(x^3\) 项。）于是
\[ \bar{E}_{\text{离散}} = kT\left(1-\frac{x}{2}+\frac{x^2}{12}-\frac{x^4}{720}+\cdots\right) \]
差值：
\[ \Delta = kT\left(\frac{x}{2}-\frac{x^2}{12}+\frac{x^4}{720}-\cdots\right) \]
代回 \(x=h\nu/kT\)：
\[ \boxed{\Delta \approx \frac{h\nu}{2} - \frac{(h\nu)^2}{12kT} + \frac{(h\nu)^4}{720(kT)^3}-\cdots} \]
一阶修正为 \(\frac{h\nu}{2}\)，与 \(T\) 无关；也就是说低频端离散化带来的能量亏损主要由零点能式的 \(\frac{h\nu}{2}\) 主导（普朗克原始假设没有零点能，这里是数学展开的体现）。相对误差：
\[ \frac{\Delta}{kT}\approx \frac{x}{2}=\frac{h\nu}{2kT}\ll 1 \]
所以低频时离散与连续几乎一致。
3. 大 \(x\) 展开（高频，\(h\nu \gg kT\)）
此时 \(e^{-x}\ll 1\)，
\[ \frac{x}{e^{x}-1} = x e^{-x}\frac{1}{1-e^{-x}} = x\left(e^{-x}+e^{-2x}+e^{-3x}+\cdots\right) \]
所以
\[ \bar{E}_{\text{离散}} = kT, x\left(e^{-x}+e^{-2x}+\cdots\right) = h\nu\left(e^{-x}+e^{-2x}+\cdots\right) \]
差值：
\[ \Delta = kT - h\nu\left(e^{-x}+e^{-2x}+\cdots\right) \approx kT - h\nu, e^{-x} \]
即
\[ \boxed{\Delta \approx kT - h\nu, e^{-h\nu/kT}\quad (h\nu\gg kT)} \]
此时 \(\bar{E}_{\text{离散}}\approx h\nu,e^{-x}\to 0\)，而连续值仍是 \(kT\)，差值趋近于 \(kT\)。相对差 \(\Delta/kT\to 1\)：高频模式在离散情形下几乎完全不被激发，连续情形却仍分到 \(kT\)。
4. 小结
| 区域 | 离散平均能量 | 与 \(kT\) 的差 \(\Delta\) |
|---|---|---|
| \(h\nu\ll kT\) | \(\approx kT-\frac{h\nu}{2}+\frac{(h\nu)^2}{12kT}\) | \(\approx \frac{h\nu}{2}\)，相对差 \(\frac{h\nu}{2kT}\ll1\) |
| \(h\nu\sim kT\) | \(kT\cdot\frac{x}{e^x-1}\) | \(kT\left(1-\frac{x}{e^x-1}\right)\)，量级为 \(kT\) |
| \(h\nu\gg kT\) | \(\approx h\nu,e^{-h\nu/kT}\to0\) | \(\approx kT\)，相对差趋近 1 |
离散与连续的差别本质上是 \(x=h\nu/kT\) 的函数：低频处差别只是 \(O(h\nu)\) 的小修正，高频处差别是量级 \(kT\) 的根本性压制。正是这个高频端的指数压制，使得总能量密度 \(\int \bar{E}(\nu),g(\nu),d\nu\) 收敛，紫外灾难被消除。
黑体辐射的普朗克公司薛定谔方程，从而得到薛定谔方程的概率解释，也就是能量处于多少的概率
1. 玻尔兹曼统计的核心结构
玻尔兹曼统计给出：系统处于能量为 \(E_n\) 的状态的概率为 \[ P_n = \frac{e^{-E_n/kT}}{Z},\qquad Z=\sum_n e^{-E_n/kT} \] 平均能量 \(\bar{E}=\sum_n E_n P_n\)。普朗克推导中，把振子能级取为 \(E_n=nh\nu\)，就得到 \[ \bar{E}=\frac{h\nu}{e^{h\nu/kT}-1} \] 这一步是纯统计力学的，没有用到薛定谔方程。关键结构是：一个指数权重 \(e^{-E/kT}\) 乘以态密度/能级简并，再做归一化。
2. 普朗克公式与“概率”的形式对应
普朗克黑体辐射公式（按频率的能量密度）为 \[ u(\nu),d\nu = \frac{8\pi h\nu^3}{c^3}\frac{1}{e^{h\nu/kT}-1},d\nu \] 它可以读成两部分：
- \(\dfrac{8\pi\nu^2}{c^3}d\nu\)：模式数（态密度）；
- \(\dfrac{h\nu}{e^{h\nu/kT}-1}\)：每个模式的平均能量。
而“某个模式被激发到第 \(n\) 个能级”的概率是 \[ P_n \propto e^{-nh\nu/kT} \] 这正是玻尔兹曼因子。所以黑体辐射谱本质上是能量量子数 \(n\) 的玻尔兹曼分布在连续频率上的表现。这里已经出现了“能量取某个值的概率正比于指数权重”的统计规律。
3. 从玻尔兹曼因子到量子力学的形式桥梁
薛定谔方程本身不是从黑体辐射推出来的。但两者共享一个深刻的数学结构：指数型权重与微分方程的本征值问题。具体联系有三条。
（a）虚时间与热核
把量子力学的时间演化算符 \(e^{-iHt/\hbar}\) 做虚时间替换 \(t\to -i\hbar\beta\)（\(\beta=1/kT\)），得到 \[ e^{-iHt/\hbar};\longrightarrow; e^{-\beta H} \] 这正是玻尔兹曼因子算符化后的形式。配分函数写成 \[ Z=\mathrm{Tr},e^{-\beta H}=\sum_n e^{-\beta E_n} \] 其中 \(E_n\) 是哈密顿量 \(H\) 的本征值，即定态薛定谔方程 \[ H\psi_n=E_n\psi_n \] 的解。于是：玻尔兹曼统计里的能级 \(E_n\)，在量子力学里就是薛定谔方程的本征值；玻尔兹曼权重 \(e^{-\beta E_n}\) 就是本征态上的热权重。 这是严格的形式对应（统计力学与量子力学的“虚时间”联系），也是路径积分量子统计的基础。
（b）本征值问题的统一
普朗克处理的是谐振子，其能级 \(E_n=(n+\frac12)h\nu\) 后来正是薛定谔方程对谐振子的解。薛定谔在1926年建立波动方程后，立刻用它解谐振子，得到的能级间隔正是 \(h\nu\)，与普朗克的离散能级一致。所以“能量离散”这件事，在旧量子论里是假设，在薛定谔方程里是边界条件加算符本征值问题的自然结果。离散化的物理机制从“假设”变成了“本征值谱”。
（c）概率幅与概率
玻尔兹曼统计给的是经典概率 \(P_n\propto e^{-\beta E_n}\)。量子力学给的则是概率幅 \(\psi_n\)，测量到能量 \(E_n\) 的概率是 \(|\langle \psi_n|\Psi\rangle|^2\)（玻恩规则）。两者不是同一个概率，但结构上平行：
- 统计力学：权重 \(e^{-\beta E_n}\)，归一化得概率；
- 量子力学：振幅 \(c_n=\langle \psi_n|\Psi\rangle\)，模方归一化得概率。
“能量处于多少的概率”在量子力学里就是：把态按 \(H\) 的本征态展开 \[ |\Psi\rangle=\sum_n c_n|\psi_n\rangle,\qquad H|\psi_n\rangle=E_n|\psi_n\rangle \] 则测得能量 \(E_n\) 的概率为 \[ P(E_n)=|c_n|^2 \] 而 \(c_n\) 的演化由薛定谔方程决定。热平衡时，若系统与热库耦合，密度矩阵对角元趋于 \(e^{-\beta E_n}/Z\)，又回到玻尔兹曼分布。这是两条线索的汇合点：量子概率的统计平均在热平衡下还原为玻尔兹曼权重。
4. 严格说清楚“哪些能推、哪些不能”
- 不能从普朗克黑体辐射公式“推出”薛定谔方程。黑体辐射只给出能级离散和玻尔兹曼权重，不给出波函数、叠加原理、相位演化。
- 能建立的是形式对应：\(e^{-\beta H}\) 与 \(e^{-iHt/\hbar}\) 同源（解析延拓），\(H\) 的本征值同时是玻尔兹曼能级和薛定谔定态能量。普朗克的 \(E_n=nh\nu\) 是谐振子本征值谱的早期版本。
- 概率解释（玻恩规则）不是从玻尔兹曼统计直接得到的；它是量子力学独立的公设。但热平衡下量子概率的分布确实还原为玻尔兹曼形式，这是两者一致的体现。
5. 一句话总结这条链
玻尔兹曼统计给出“能量取 \(E_n\) 的权重 \(e^{-E_n/kT}\)”；普朗克把 \(E_n\) 离散化为 \(nh\nu\)，解决了黑体辐射；薛定谔方程把 \(E_n\) 变成哈密顿量的本征值，把“权重”升级为“概率幅的模方”；热平衡下概率幅模方的统计分布又回到玻尔兹曼权重。所以这条链是形式对应与历史启发，不是逻辑推导：离散能级从假设变成本征值谱，概率从经典权重变成量子模方，而两者在虚时间替换下严格对应。
如果你愿意，下一步可以把谐振子的薛定谔方程解出来，显式验证 \(E_n=(n+\frac12)h\nu\) 与普朗克能级的对应，以及热平衡密度矩阵如何给出 \(\bar{E}=\frac{h\nu}{e^{h\nu/kT}-1}\)。