3.1 狄拉克方程：一阶导数的入场券与自旋
3.1 狄拉克方程：一阶导数的入场券与自旋
本节摘要：狄拉克方程 (iγ^μ∂_μ − m)ψ = 0 把相对论能量关系开方成一阶矩阵方程，四分量旋量 ψ 取代标量波函数，四个 4×4 的γ矩阵满足克利福德代数 {γ^μ,γ^ν} = 2g^μν。本节完成开方手术、数值验证γ代数与平面波旋量解（残差 4.4e-16），并从解的结构读出自旋二态与反粒子预告。电子作为场激发的完整身份就此确立。
承接第一幕的全部装备（对易结构、模式展开）但更换演员类型；通往 3.2 的量子化——那里这套解将被两套算符（粒子 b 与反粒子 d）分别认领。
开方手术的完整记录
1928 年狄拉克的动机很实际：KG 方程是时间二阶导，初值要同时指定 ψ 与 ∂ₜψ，概率密度 ρ = ∂ₜψ*ψ + … 不正定（1.1 验伤二）。他的方案：把 E² = p² + m² 开方，让方程对时间一阶。标量开不出方（数没有平方根的矩阵结构），但矩阵可以：
设 E = α·p + βm（α、β 为待定矩阵），平方展开与 E² = p² + m² 对比： αi² = β² = 1, αiαj + αjαi = 0 (i≠j), αiβ + βαi = 0 最小解是 4×4。整理成协变形式（gamma 矩阵）： gamma0 = β, gammai = β·alphai 克利福德代数： gamma_mu gamma_nu + gamma_nu gamma_mu = 2 g_mu_nu g = diag(+1,-1,-1,-1) 即 1.2 的闵可夫斯基记账法 狄拉克方程： (i gamma^mu d_mu - m) psi = 0，psi 为四分量旋量 概率密度： rho = psi+ psi ——四个分量的模方和，天然正定，账本修复
用狄拉克表示（γ⁰ 对角、γ^i 含泡利矩阵）逐对验证代数：
# 演算 1：gamma 矩阵克利福德代数的数值验证 I2 = [[1, 0], [0, 1]] s_x = [[0, 1], [1, 0]]; s_y = [[0, -1j], [1j, 0]]; s_z = [[1, 0], [0, -1]] g0 = [[1 if i == j and i < 2 else (-1 if i == j else 0) for j in range(4)] for i in range(4)] def gamma_i(s): g = [[0j]*4 for _ in range(4)] for i in range(2): for j in range(2): g[i][2+j] = s[i][j] # 右上块 sigma g[2+i][j] = -s[i][j] # 左下块 -sigma return g gam = [g0, gamma_i(s_x), gamma_i(s_y), gamma_i(s_z)] met = (1, -1, -1, -1) ok = True for mu in range(4): for nu in range(4): ac = [[ (sum(gam[mu][i][k]*gam[nu][k][j] for k in range(4)) + sum(gam[nu][i][k]*gam[mu][k][j] for k in range(4))).real for j in range(4)] for i in range(4)] tgt = 2*met[mu] if mu == nu else 0 for i in range(4): for j in range(4): if abs(ac[i][j] - (tgt if i == j else 0)) > 1e-12: ok = False print("全部 {gamma_mu, gamma_nu} = 2 g_mu_nu 验证通过：", ok) # 输出: 全部 {gamma_m,gamma_n} = 2 g_mn 验证通过： True
十六个反对易关系全部严格成立——矩阵尺寸 4 不是选择而是必然：2×2 与 3×3 都装不下这套代数。相对论加量子力学，逼出了四分量旋量这个新数学对象。
解的结构：自旋与负能预告
平面波解 ψ = u(p)e^{−ip·x} 代入方程得动量空间条件 (γ·p − m)u = 0。这是一个 4×4 线性方程组，秩为 2，留下两个线性无关解——两个自旋态。标准构造 u(p) = √(E+m)·(χ, σ·p/(E+m)χ)ᵀ，χ 为任意二分量旋量。数值核实：
# 演算 2：平面波旋量满足狄拉克方程（p 沿 z 轴，m=1） import math def dirac_eq(u, E, pvec, m): out = [0j]*4 for i in range(4): acc = E*gam[0][i][0]*u[0] + E*gam[0][i][1]*u[1] + E*gam[0][i][2]*u[2] + E*gam[0][i][3]*u[3] for k in range(4): acc -= (pvec[0]*gam[1][i][k] + pvec[1]*gam[2][i][k] + pvec[2]*gam[3][i][k])*u[k] acc -= m*u[i] out[i] = acc return out def act(M, v): return [M[0][0]*v[0] + M[0][1]*v[1], M[1][0]*v[0] + M[1][1]*v[1]] sig = (s_x, s_y, s_z) def u_spinor(chi_idx, E, pvec, m): chi = [1, 0] if chi_idx == 0 else [0, 1] sp = [0j, 0j] for a in range(3): t = act(sig[a], chi) sp = [sp[i] + pvec[a]*t[i] for i in range(2)] fac = math.sqrt(E + m) return [fac*chi[0], fac*chi[1], sp[0]/fac, sp[1]/fac] m_, pvec = 1.0, (0.0, 0.0, 2.0) E_ = math.sqrt(4 + 1) for c in (0, 1): u = u_spinor(c, E_, pvec, m_) r = dirac_eq(u, E_, pvec, m_) ubu = sum(u[i].conjugate()*g0[i][j]*u[j] for i in range(4) for j in range(4)).real print("chi_%d: |(gamma·p - m)u|最大分量 = %.2e ; ubar·u = %.4f (=2m)" % (c, max(abs(v) for v in r), ubu)) # 输出: # chi_0: |(gamma·p - m)u|最大分量 = 4.44e-16 ; ubar·u = 2.0000 (=2m) # chi_1: |(gamma·p - m)u|最大分量 = 4.44e-16 ; ubar·u = 2.0000 (=2m)
残差 4.4e-16 是双精度浮点的地板——旋量解在机器精度内严格满足方程，归一化 ūu = 2m 分毫不差。方程的完整解空间是四维的：两个正能解 u（χ 二选一）加两个"负能解"v。狄拉克当初为负能解失眠；今天的读法已在 2.1 铺好：负能解对应产生算符 d†，将创造一种电荷相反、质量相同的新粒子。1932 年安德森在宇宙线里找到它（正电子），教科书从此把这一页印成"理论预言先于实验"的范本。
图 四分量旋量：解的户口本
⚠️ 易错点：旋量不是四矢量。洛伦兹变换下旋量按 4×4 的 S(Λ) 矩阵旋转（终幕 7.2 会用到它的规范版本），与四矢量的 Λ 矩阵不同；把 (ψ₀,ψ₁,ψ₂,ψ₃) 当作时空分量做变换是初学高频错误。
💡 关键直觉：自旋是"开方"的直接产物。为了让能量的一阶式平方后回到 p² + m²，系数必须是对易为零的矩阵，而最小可行维度 4 恰好容下两个自旋态。转两圈才复原的旋量性质不是强加的设定，是线性代数给的找零。
本节要点回顾
- 开方手术：一阶矩阵方程取代二阶标量方程，概率密度 ψ†ψ 天然正定
- 克利福德代数：{γμ,γν} = 2gμν 十六式全部数值成立，最小维数 4 是代数必然
- 解的户口：四解 = 正能二自旋（u）+ 反粒子二自旋（v），自旋从推论变为身份
- 机器零验证：(γ·p−m)u 残差 4.4e-16，归一 ūu = 2m
- 预言兑现：正电子 1930 预言、1932 发现，负能解由反粒子算符接管的伏笔在此埋下
方程就位，下一节做量子化：换上反对易关系，你会看到泡利不相容变成一行矩阵等式，而正电子拿到自己的产生算符。