为什么普朗克公司中，用能量离散化就可以解决黑体辐射问题，物理机制是什么？或者统计规律是什么，为什么连续的就不行 普朗克用能量离散化解决黑体辐射问题，核心机制在于：离散化改变了统计物理中“一份能量”与“热运动典型能量”的比较方式，从而在高频端天然地“冻结”了那些本该被激发的振动模式。 经典连续理论之所以失败，是因为它不加区分地给所有模式分配相同的平均能量，导致频率越高、模式越多、总能量越发散。 经典连续理论为什么“不行” 经典物理用瑞利—金斯公式描述黑体辐射。它基于两个看似无懈可击的假设： 腔内有无限多的电磁振动模式，频率越高，模式数量越多（按频率的平方增长）。 能量均分定理：在热平衡下，每个模式的平均能量都等于 \(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\) 很大时，分母中的指数项急剧增大，平均能量被指数式压制到近乎为零。高频模式数量再多，乘以一个趋近于零的平均能量，总贡献就变得有限。紫外灾难由此被“截断”。 所以，离散化并不直接改变电磁场的模式数量，它改变的是能量分配的统计规则：让高频模式因为“单价太高”而无法参与热平衡的能量分配。普朗克本人起初只把这当作一个数学上的“绝望之举”，但后来的发展表明，这个假设触及了能量交换过程本身的不连续性。 平均能量是怎么算出来的 平均能量 \(\bar{E}=\frac{h\nu}{e^{h\nu/kT}-1}\) 的推导，核心是把一个振子看成与温度为 \(T\) 的热库平衡的系统，用离散能级的玻尔兹曼分布求统计平均。下面给出完整计算。 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}\)。 ...

不使用球坐标的三维谐振子方程 如果不使用球坐标，改用直角坐标处理三维谐振子中的狄拉克方程，思路是：把三维问题分解成三个一维狄拉克谐振子，再用升降算符或多项式方法求解。 下面取最常见的等量标量-矢量耦合： \[ S(\mathbf r)=V(\mathbf r)=\frac12 m\omega^2 r^2 \] 此时狄拉克方程可写为： \[ \left[\boldsymbol{\alpha}\cdot\mathbf p+\beta\big(m+S\big)\right]\psi=(E-V)\psi \] 在自然单位 \(\hbar=c=1\) 下，令 \(S=V=\frac12 m\omega^2 r^2\)，则： \[ \left[\boldsymbol{\alpha}\cdot\mathbf p+\beta m+\frac12 m\omega^2 r^2(\beta+1)\right]\psi=E\psi \] 关键简化 取狄拉克矩阵的表示： \[ \beta= \begin{pmatrix} I&0\ 0&-I \end{pmatrix},\qquad \boldsymbol{\alpha}= \begin{pmatrix} 0&\boldsymbol{\sigma}\ \boldsymbol{\sigma}&0 \end{pmatrix} \] 把四分量旋量写成： \[ \psi= \begin{pmatrix} \phi\ \chi \end{pmatrix} \] 其中 \(\phi,\chi\) 是二分量旋量。代入后得到耦合方程： \[ \boldsymbol{\sigma}\cdot\mathbf p,\chi + m\phi + \frac12 m\omega^2 r^2(1+1)\phi = E\phi \] \[ \boldsymbol{\sigma}\cdot\mathbf p,\phi - m\chi + \frac12 m\omega^2 r^2(-1+1)\chi = E\chi \] ...

狄拉克所担忧的“无穷大问题”，最终并不是被“消除”了，而是通过重整化（Renormalization）技术被“驯服”和“回避”了。物理学家们发展出一套严格的数学程序，将计算中出现的无穷大吸收到对电子质量和电荷的重新定义中。 🤔 无穷大从何而来 在量子电动力学中，电子会不断发射和吸收虚光子，这些虚光子又可能产生短暂的虚电子-正电子对。计算这些过程时，会涉及到对虚粒子所有可能动量的积分，结果就是发散的。这种无穷大主要出现在两个地方：电子的自能（导致质量修正）和真空极化（导致电荷修正）。 🧮 解决方案：重整化 20世纪40年代末，朝永振一郎、朱利安·施温格、理查德·费曼和弗里曼·戴森等人完善了重整化方法。其核心思想是： 参数重定义：我们实验上测量的电子质量和电荷（观测值），实际上已经包含了这些虚粒子相互作用的贡献。理论公式中原始的“裸”质量和“裸”电荷本身是无穷大的，但经过重新定义后，理论计算出的可观测值与实验相符。 吸收发散：通过引入正规化（如维度正规化）暂时让积分变得有限，然后通过添加抵消项，将无穷大精确地抵消掉，最终留下有限的、与正规化参数无关的物理结果。 ⚡ 实验验证与遗留争议 重整化后的量子电动力学（QED）取得了惊人成功。它极其精确地预测了兰姆位移（Lamb shift）和电子的反常磁矩（g-2），理论值与实验值在小数点后多位都高度吻合。 然而，狄拉克和费曼本人至死都对这种“把无穷大从后门扫出去”的做法不满。狄拉克曾批评说，这就像“忽视一个无穷大的量，仅仅因为你不想要它”，在数学上并不合理。费曼也认为重整化是“一种戏法”。尽管如此，重整化至今仍是量子场论的标准工具，只是它仍属于一种基于微扰论的系统性“回避”，而非对理论基础的彻底重构。

下面从四分量狄拉克方程出发，在非相对论极限下推出泡利方程。全程使用四分量形式，不预先化成二分量方程。 1. 有电磁场的四分量狄拉克方程 取 Dirac 表象： \[ \psi= \begin{pmatrix} \psi_1\ \psi_2\ \psi_3\ \psi_4 \end{pmatrix} = \begin{pmatrix} \phi\ \chi \end{pmatrix}, \qquad \phi= \begin{pmatrix} \psi_1\ \psi_2 \end{pmatrix}, \quad \chi= \begin{pmatrix} \psi_3\ \psi_4 \end{pmatrix} \] 最小耦合： \[ \boldsymbol{\pi}=-i\hbar\nabla-\frac{e}{c}\mathbf{A}, \qquad i\hbar\partial_t\to i\hbar\partial_t-eA_0 \] 四分量方程写成二分块形式： \[ i\hbar\partial_t\phi = eA_0\phi+mc^2\phi+c,\boldsymbol{\sigma}\cdot\boldsymbol{\pi},\chi \tag{1} \] \[ i\hbar\partial_t\chi = eA_0\chi-mc^2\chi+c,\boldsymbol{\sigma}\cdot\boldsymbol{\pi},\phi \tag{2} \] 2. 提取非相对论极限 非相对论极限下，能量主要来自静止质量： \[ \psi\sim e^{-imc^2t/\hbar} \begin{pmatrix} \tilde\phi\ \tilde\chi \end{pmatrix} \] 令 \[ \phi=e^{-imc^2t/\hbar}\tilde\phi, \qquad \chi=e^{-imc^2t/\hbar}\tilde\chi \] 代入 (1)、(2)： ...

使用上面的4分量方程来证明g=2，不要使用二分量方程 下面用你前面列出的四分量狄拉克方程，在非相对论极限下推导电子与磁场的相互作用，并证明朗德 \(g\) 因子为 \[ \boxed{g=2} \] 全程只用四分量方程，不用预先化成二分量方程。 1. 有电磁场时的四分量方程 取 Dirac 表象，\(\psi=(\psi_1,\psi_2,\psi_3,\psi_4)^T\)。有电磁场时，最小耦合为 \[ i\hbar\partial_t \to i\hbar\partial_t-eA_0, \qquad -i\hbar\nabla \to -i\hbar\nabla-\frac{e}{c}\mathbf{A} \] 令 \[ \boldsymbol{\pi}=-i\hbar\nabla-\frac{e}{c}\mathbf{A} \] 则四分量方程可写成二分块形式： \[ i\hbar\frac{\partial}{\partial t} \begin{pmatrix}\phi \\ \chi\end{pmatrix} = eA_0 \begin{pmatrix}\phi \\ \chi\end{pmatrix} + mc^2 \begin{pmatrix}\phi \\ -\chi\end{pmatrix} + c,\boldsymbol{\sigma}\cdot\boldsymbol{\pi} \begin{pmatrix}\chi \\ \phi\end{pmatrix} \] 其中 \[ \phi=\begin{pmatrix}\psi_1 \\ \psi_2\end{pmatrix}, \qquad \chi=\begin{pmatrix}\psi_3 \\ \psi_4\end{pmatrix} \] 写成两个方程： \[ i\hbar\partial_t\phi = eA_0\phi+mc^2\phi+c,\boldsymbol{\sigma}\cdot\boldsymbol{\pi},\chi \tag{1} \] \[ i\hbar\partial_t\chi = eA_0\chi-mc^2\chi+c,\boldsymbol{\sigma}\cdot\boldsymbol{\pi},\phi \tag{2} \] ...

狄拉克方程 \[ (i\gamma^\mu \partial_\mu - m)\psi = 0 \] （自然单位制 \(\hbar=c=1\)）在 Dirac 表示下，取 \[ \psi=\begin{pmatrix}\psi_1\ \psi_2\ \psi_3\ \psi_4\end{pmatrix}, \qquad \gamma^0=\begin{pmatrix}I&0\0&-I\end{pmatrix},\quad \gamma^i=\begin{pmatrix}0&\sigma^i\-\sigma^i&0\end{pmatrix} \] 则四个分量方程为： \[ i\frac{\partial \psi_1}{\partial t} = m\psi_1 + \left(-i\frac{\partial}{\partial x}+\frac{\partial}{\partial y}\right)\psi_4 + i\frac{\partial \psi_3}{\partial z} \] \[ i\frac{\partial \psi_2}{\partial t} = m\psi_2 + \left(-i\frac{\partial}{\partial x}-\frac{\partial}{\partial y}\right)\psi_3 - i\frac{\partial \psi_4}{\partial z} \] \[ i\frac{\partial \psi_3}{\partial t} = -m\psi_3 + \left(-i\frac{\partial}{\partial x}+\frac{\partial}{\partial y}\right)\psi_2 + i\frac{\partial \psi_1}{\partial z} \] \[ i\frac{\partial \psi_4}{\partial t} = -m\psi_4 + \left(-i\frac{\partial}{\partial x}-\frac{\partial}{\partial y}\right)\psi_1 - i\frac{\partial \psi_2}{\partial z} \] ...

好，我们来解氢原子能级。这是狄拉克方程最经典的应用之一。 不过要先说清楚：完整的求解过程比较长，因为要处理库仑势、角动量分解、径向方程。我先把关键步骤和结果写出来，让你看到能级是怎么出来的。 一、从自由方程到氢原子 自由狄拉克方程： \[ (\gamma^\mu p_\mu - mc)\psi=0 \] 在电磁场中，做最小耦合： \[ p_\mu \to p_\mu - \frac{e}{c}A_\mu \] 氢原子只有库仑势： \[ A_0=\frac{Ze}{r},\qquad \vec{A}=0 \] 所以： \[ \left[\gamma^0\left(p_0-\frac{Ze^2}{r}\right)+\vec{\gamma}\cdot\vec{p}-mc\right]\psi=0 \] 注意这里 \(p_0=E/c\)，\(Ze^2/r\) 是势能除以 \(c\) 后的量。 二、分块形式 用 \(\psi_L,\psi_R\)： \[ \left(p_0-\frac{Ze^2}{r}-mc\right)\psi_L+(\vec{\sigma}\cdot\vec{p})\psi_R=0 \] \[ (\vec{\sigma}\cdot\vec{p})\psi_L+\left(-p_0-\frac{Ze^2}{r}-mc\right)\psi_R=0 \] 这里 \(Ze^2/r\) 是势能项，量纲是动量。 三、角动量分解 在球坐标下，角动量算符： \[ \vec{J}=\vec{L}+\vec{S} \] 其中： \[ \vec{S}=\frac12\vec{\Sigma} \] 狄拉克方程的解可以按 \(j\) 分类。 定义： \[ K=\gamma^0(\vec{\Sigma}\cdot\vec{L}+1) \] 它的本征值： \[ \kappa=\pm(j+\tfrac12) \] 对于 \(j=l+\frac12\)： \[ \kappa=-(j+\tfrac12) \] 对于 \(j=l-\frac12\)： ...

分块矩阵形式 好，我们直接写出这两个方程和 \(\vec{\sigma}\cdot\vec{p}\) 的具体矩阵。 一、两个二分量方程 前面已经得到： \[ (p_0-m)\psi_L+(\vec{\sigma}\cdot\vec{p})\psi_R=0 \] \[ (\vec{\sigma}\cdot\vec{p})\psi_L+(p_0+m)\psi_R=0 \] 其中： \[ \psi_L=\begin{pmatrix}\psi_1\\psi_2\end{pmatrix},\qquad \psi_R=\begin{pmatrix}\psi_3\\psi_4\end{pmatrix} \] 二、\(\vec{\sigma}\cdot\vec{p}\) 的具体值 泡利矩阵： \[ \sigma_1=\begin{pmatrix}0&1\1&0\end{pmatrix} \] \[ \sigma_2=\begin{pmatrix}0&-i\i&0\end{pmatrix} \] \[ \sigma_3=\begin{pmatrix}1&0\0&-1\end{pmatrix} \] 所以： \[ \vec{\sigma}\cdot\vec{p} \sigma_1 p_1+\sigma_2 p_2+\sigma_3 p_3 \] 逐项写： \[ \sigma_1 p_1= \begin{pmatrix}0&p_1\p_1&0\end{pmatrix} \] \[ \sigma_2 p_2= \begin{pmatrix}0&-ip_2\ip_2&0\end{pmatrix} \] \[ \sigma_3 p_3= \begin{pmatrix}p_3&0\0&-p_3\end{pmatrix} \] 相加： \[ \vec{\sigma}\cdot\vec{p} \begin{pmatrix} p_3 & p_1-ip_2\ p_1+ip_2 & -p_3 \end{pmatrix} \] 三、把两个方程完全展开成普通方程 第一个方程 \[ (p_0-m)\psi_L+(\vec{\sigma}\cdot\vec{p})\psi_R=0 \] ...

两个分量方程：能否直接把第二个的X代入第一个，从而消去X 可以，但要注意两点：这不是简单地“把第二个方程解出 \(\chi\) 再代入第一个”就完事，因为第二个方程里同时含有 \(\chi\) 和 \(\phi\)，而且含有时间导数与空间导数。要真正消去 \(\chi\)，需要对第二个方程再作用一次相应的算子，利用一阶方程把 \(\chi\) 消掉，最后得到 \(\phi\) 的二阶方程。 1. 先把方程写清楚 你写的两个方程是狄拉克方程在 Weyl/手征基下的形式。通常写成 \[ i(\partial_t+ie\phi)\phi+i\sigma^i(\partial_i+ieA_i)\chi-m\phi=0, \] \[ -i(\partial_t+ie\phi)\chi-i\sigma^i(\partial_i+ieA_i)\phi-m\chi=0. \] 这里 \(\phi,\chi\) 是二分量旋量。为了方便，定义协变导数 \[ D_t=\partial_t+ie\phi,\qquad D_i=\partial_i+ieA_i. \] 于是方程变为 \[ iD_t\phi+i\sigma^iD_i\chi-m\phi=0, \tag{1} \] \[ -iD_t\chi-i\sigma^iD_i\phi-m\chi=0. \tag{2} \] 2. 为什么不能直接“解出 \(\chi\) 代入”？ 从 (2) 看，\(\chi\) 满足 \[ (-iD_t-m)\chi-i\sigma^iD_i\phi=0. \] 形式上可以写成 \[ (iD_t+m)\chi=-i\sigma^iD_i\phi. \] 如果 \(iD_t+m\) 有逆，可以形式解出 \[ \chi=-(iD_t+m)^{-1}i\sigma^iD_i\phi. \] 但这里 \((iD_t+m)\) 是一个微分算子，不是普通数。把它代入 (1) 会得到 \[ iD_t\phi+i\sigma^iD_i\left[-(iD_t+m)^{-1}i\sigma^jD_j\phi\right]-m\phi=0, \] 这虽然形式上消去了 \(\chi\)，但出现了非局域的逆算子，而且时间导数在分母上，实际使用很不方便。物理上更标准的方法是消去 \(\chi\)，得到一个只含 \(\phi\) 的二阶方程。 ...