3.2 费米子量子化与反粒子：一模式一座位

3.2 费米子量子化与反粒子：一模式一座位

本节摘要：旋量场的量子化把对易关系换成反对易关系 {b,b†} = 1，占据数只能取 0 或 1——泡利不相容从公理降级为代数推论。本节用约当-维格纳矩阵在 8×8 空间实测反对易代数与占据数的严格幂等性（n² = n），登记粒子与反粒子两套算符，并把 1.3 的诺特电荷升级为算符 Q = Σ(b†b − d†d)。物质的稳定性从此有了代数担保。

承接 3.1 的四解结构：u 解与 v 解各配一套算符；通往 3.3 的离散对称性——正反粒子之间的镜像关系 C 将在那里正式定名。

规矩要换：从对易到反对易

第一幕的标量场用对易关系 [a, a†] = 1，得到整数占据 0,1,2,…。如果旋量场照抄，会立刻撞上两堵墙：其一，自旋统计定理（与 2.5 的因果性证明同源）禁止自旋半奇数的粒子用对易关系量子化——因果律与正定能量不能同时保住；其二，实验事实直白：原子壳层结构要求每个量子态至多一个电子。出路是把换位运算整个换掉：

费米量子化：  {b_i, b_j+} = δ_ij     （反对易子 AB+BA）
              {b_i, b_j} = 0 = {b_i+, b_j+}
推论一：      b_i+ b_i+ = 0          （同一个态上产生两个粒子 = 零向量）
推论二：      占据数 n_i = b_i+ b_i 只取 0 或 1（幂等：n² = n）
推论三：      交换两个费米子的态自动差一个负号（反对称波函数）
​

泡利不相容在这里不是被假设的规矩，而是 {b, b†}=1 的第二行推论。用约当-维格纳构造在三格座位上实测这套代数：

# 演算 1：约当-维格纳链上的反对易关系与幂等占据（3 个座位 -> 8 维）
def kp(A, B):
    r, s = len(B), len(B[0])
    return [[A[i//r][j//s]*B[i % r][j % s] for j in range(len(A[0])*s)] for i in range(len(A)*r)]
def mmul(A, B):
    return [[sum(A[i][k]*B[k][j] for k in range(len(B))) for j in range(len(B[0]))] for i in range(len(A))]
Zm = [[1, 0], [0, -1]]; Pm = [[0, 0], [1, 0]]; Id = [[1, 0], [0, 1]]
def c_lower(site, L=3):
    mats = [Zm if s < site else (Pm if s == site else Id) for s in range(L)]
    out = mats[0]
    for M in mats[1:]:
        out = kp(out, M)
    return out
def dag(M):
    return [[M[j][i].conjugate() for j in range(len(M))] for i in range(len(M[0]))]
cs = [c_lower(s) for s in range(3)]
csd = [dag(c) for c in cs]
def antic(A, B):
    n = len(A)
    return [[(sum(A[i][k]*B[k][j] for k in range(n)) + sum(B[i][k]*A[k][j] for k in range(n))).real
             for j in range(n)] for i in range(n)]
ok = True
for i in range(3):
    for j in range(3):
        v = antic(cs[i], csd[j])
        tgt = 1 if i == j else 0
        for r in range(8):
            for c in range(8):
                if abs(v[r][c] - (tgt if r == c else 0)) > 1e-12: ok = False
print("全部 {c_i, c_j+} = delta_ij 验证：", ok)
n1 = mmul(csd[1], cs[1])
nn = mmul(n1, n1)
print("max |n^2 - n| = %.2e  （严格幂等：占据只能是 0 或 1）" %
      max(abs(nn[i][j]-n1[i][j]) for i in range(8) for j in range(8)))
# 输出:
# 全部 {c_i, c_j+} = delta_ij 验证： True
# max |n^2 - n| = 0.00e+00  （严格幂等：占据只能是 0 或 1）
​

n² − n 精确为零——不是浮点零，是矩阵多项式的恒等式。对比玻色一方（2.2 的截断矩阵对角 0,1,2,3,4），两种统计的差别一行代码就能晒出来。值得多想一步这行代数买到的东西：原子里的电子必须填不同的座位，座位的排布方式就是元素周期表；白矮星与中子星抵抗引力坍缩的压强（简并压）同样是它。你身体里每一个原子不塌缩，都在向 {b,b†}=1 交税。

两套算符与电荷账本

狄拉克场的完整模式展开需要两套产生湮灭算符：b、b† 认领正能解 u（电子），d、d† 认领 v 解（正电子）。历史版本（狄拉克海）把反粒子读成"负能海里的空穴"，直觉生动但账目笨重：真空要先被填成无限满的海。现代记账（福克空间直和）干净得多：

场展开：  psi(x) = Σ_s ∫ d^3p [ b_ps u_s(p) e^{-ipx} + d_ps+ v_s(p) e^{+ipx} ]
哈密顿量：H = Σ ∫ E_p ( b+ b + d+ d )        粒子反粒子平权，全正
电荷算符：Q = Σ ∫ ( b+ b - d+ d )            电子 +1，正电子 -1
真空：    b|0> = d|0> = 0                     无需先填海
​

注意 e^{+ipx} 与 d† 配对——负能解的"时间倒流"感在算符层面被翻译成"产生一个正能量反粒子"，这正是 2.1 埋下伏笔的兑现。电荷算符的形式继承 1.3 的诺特流：Q 的被积表达式正是 j⁰ = ψ†ψ 的正规序版本。粒子与反粒子电荷相反、质量严格相同——CPT 定理（3.3）担保这个"相同"精确到小数点后所有位。

一模式一座位的剧场后果

⚠️ 易错点：反对易关系不等于"负号的对易关系"。{A,B} = AB+BA 与 [A,B] = AB−BA 是两种不同的双线性；把费米算符当"带负号的玻色算符"做代数，符号错误会一路传染到第五幕的圈图计算（费米子圈的负号就藏在这些细节里）。

💡 关键直觉：自旋统计定理把"自旋 1/2"与"反对易"锁死在一起，不是约定而是因果律的要求。用错统计的假想世界要么违反微观因果，要么能量无下界——剧场的账本与法律必须同时成立。

本节要点回顾

• 换规矩：费米量子化用 {b,b†} = δ，泡利不相容是它的代数推论
• 幂等实测：约当-维格纳链上 n² − n 严格为零，占据只取 0 或 1
• 双套算符：b 认领电子、d 认领正电子，负能解翻转为反粒子的产生支
• 电荷账本：Q = Σ(b†b − d†d)，继承诺特流，正反电荷相加为零
• 稳定性来源：物质不坍缩、周期表结构、简并压，全部押在反对易关系上

正反粒子已是对等演员，下一节给他们照三面镜子：C、P、T 三种离散变换，外加费曼图计算的核心零件——旋量传播子与自旋求和。