不使用球坐标的三维谐振子方程 如果不使用球坐标，改用直角坐标处理三维谐振子中的狄拉克方程，思路是：把三维问题分解成三个一维狄拉克谐振子，再用升降算符或多项式方法求解。 下面取最常见的等量标量-矢量耦合： \[ 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\) 的二阶方程。 ...

写出具体推导过程 下面把“从狄拉克方程到非相对论极限下出现 \(-\frac{e\hbar}{2mc}\boldsymbol{\sigma}\cdot\mathbf{B}\)”这个过程完整写一遍。为了不使符号过于臃肿，取 \(\hbar=c=1\)，最后再恢复常数。 1. 狄拉克方程与最小耦合 自由狄拉克方程： \[ (i\gamma^\mu\partial_\mu - m)\psi = 0 \] 引入电磁场，做最小耦合： \[ \partial_\mu \to D_\mu = \partial_\mu + i e A_\mu \] 于是： \[ (i\gamma^\mu D_\mu - m)\psi = 0 \] 即： \[ \left[i\gamma^\mu(\partial_\mu + ieA_\mu) - m\right]\psi = 0 \] 分开时间与空间： \[ i\gamma^0\partial_t + i\gamma^i\partial_i - e\gamma^0 A_0 - e\gamma^i A_i - m = 0 \] 用 \(A_0=\phi\)，\(A_i\) 对应矢量势 \(\mathbf{A}\)。取狄拉克表示： \[ \gamma^0 = \begin{pmatrix} I & 0 \\ 0 & -I \end{pmatrix},\qquad \gamma^i = \begin{pmatrix} 0 & \sigma^i \\ -\sigma^i & 0 \end{pmatrix} \] ...