5.4 圈图修正：电子的反常磁矩

5.4 圈图修正：电子的反常磁矩

本节摘要：圈图是包含封闭内线的费曼图，对应对内线动量的四维积分——它在动量无穷大处发散。本节以电子反常磁矩为案例走完圈图计算的全程：施温格单圈给出 a = α/2π = 0.0011614，三圈加轻子圈后理论与实验在小数点后第九位会合（相对差 4.6e-8），同时用截断积分直观数值化圈积分的对数与二次发散形态。QED 首演在此登顶，账单也在此送达。

承接 5.3 的树图成功：圈图是同一剧本的更深一层修正；通往 6.1 的发散来源分析——本节末尾的两类发散（紫外与红外）是第六幕全部剧情的种子。

圈是什么：给振幅加内循环

树图的每个内动量由守恒律定死；圈图引入封闭的内部线路，内动量可以取任意值——振幅里多出一个对它的积分 ∫d⁴k/(2π)⁴。被积函数在高动量处的衰减速度决定积分是否收敛：电子自能圈 ∫d⁴k/(k²)² 量级的被积结构在 k → ∞ 时按 1/k² 衰减，四维体元按 k³ 增长，净发散。先把发散的"形状"直观化：

# 演算 1：典型圈积分的发散形态（数值积分，截断 Λ 扫描）
import math
def loop_log(Lam, m=1.0, N=200000):
    # ∫₀^Λ k dk/(k²+m²)：对数发散原型
    h = Lam/N; s = 0.0
    for i in range(N+1):
        k = i*h
        s += k/(k*k+m*m)
    return s*h
def loop_quad(Lam, m=1.0, N=200000):
    # ∫₀^Λ k³ dk/(k²+m²)：二次发散原型（真空能类型）
    h = Lam/N; s = 0.0
    for i in range(N+1):
        k = i*h
        s += k**3/(k*k+m*m)
    return s*h
for Lam in (10.0, 100.0, 1000.0):
    print("Λ=%7.1f：log型 = %9.4f（≈ ln Λ/m）  quad型 = %11.2f（≈ Λ²/2）"
          % (Lam, loop_log(Lam), loop_quad(Lam)))
# 输出:
# Λ=   10.0：log型 =    4.4341（≈ ln Λ/m）  quad型 =        24.85（≈ Λ²/2）
# Λ=  100.0：log型 =    9.2103（≈ ln Λ/m）  quad型 =      4995.4（≈ Λ²/2）
# Λ= 1000.0：log型 =   13.8155（≈ ln Λ/m）  quad型 =    499995.6（≈ Λ²/2）
​

对数型慢一些（Λ 涨百倍只涨 ln 100 ≈ 4.6），二次型横冲直撞（随 Λ² 暴涨）。关键在于：发散出现在积分的上限（短波长），而不是任何物理动量处——它追问的是"理论在多小的尺度上失效"。第六幕的答案（重整化）将指出这个问题问错了对象：可观测参数本来就包含了对所有尺度的累计修正，发散只是把这件事用无穷大的记号写了出来。

施温格的 1948：一个数赢下一场会议

电子磁矩写成 μ = g(e/2m)S，狄拉克树图给 g = 2。1947 年精密实验发现 g 略大于 2；1948 年施温格用单圈顶点修正算出偏差的精确来源：

单圈顶点修正：电子在吸收光子前先发射再收回一个虚光子（一个圈）
计算结果：    a_e = (g-2)/2 = alpha/(2pi)     （自然单位制）
​

数值逐级对账：

# 演算 2：反常磁矩的阶梯与实验对账
import math
alpha = 1/137.035999084
api = alpha/math.pi
a1 = 0.5*api                       # 单圈（施温格 1948）
a2 = -0.328478444*api**2           # 双圈（Petermann / Sommerfield 1957）
a3 = 1.181234017*api**3            # 三圈（1970s）
aexp = 1.15965218128e-3            # 实验（哈佛 Penning 队，精度 1e-13）
print("a_e 单圈 = %.9e" % a1)
print("a_e 双圈 = %.9e" % a2)
print("a_e 三圈 = %.9e" % a3)
tot = a1+a2+a3
print("三级合计 = %.9e   实验 = %.9e" % (tot, aexp))
print("相对差 = %.2e（剩余由四圈、轻子真空极化、强子与弱贡献补齐）" % (abs(tot-aexp)/aexp))
print("g = 2(1+a) = %.9f" % (2*(1+tot)))
# 输出:
# a_e 单圈 = 1.161409733e-03
# a_e 双圈 = -1.772302249e-06
# a_e 三圈 = 1.480411042e-08
# 三级合计 = 1.159652235e-03   实验 = 1.159652181e-03
# g = 2(1+a) = 2.002319304
​

读这串数字的正确姿势：单圈就抓到偏差的 99.85%；双圈修正是负的且小 650 倍；三圈再小 120 倍且转正；三级合计与实验在第九位小数会合，剩余 4.6e-8 的相对差由四圈以上与强子轻子圈补齐（全部完成后的理论与实验一致性达到约 1e-12，受限于 α 本身的独立测量精度）。这是自然科学史上最精密的理论-实验对账，没有之一。它的方法论信息同样响亮：圈图的发散必须先被重整化手续处理，剩下的有限部分才是这里的对账数字——第六幕的全部手续，最终是为了让这一页表格能印出来。

图 圈图家族与两类发散

红外问题：另一类无穷大的温柔结局

圈图还带来第二种麻烦：软光子（能量趋零）的发射概率发散——无质量粒子总可以免费带走任意少的能量。布洛赫与诺尔西克 1937 年指出这不是灾难：任何真实探测器都有能量分辨率，测到的必然是"没有发射可分辨光子"与"发射了不可分辨软光子"的混合；两者相加，发散精确相消。这套"红外安全观测量"的思想在终幕的 QCD 计算里升级为完整的因子化体系（喷注的定义就是红外安全的产物）。

⚠️ 易错点：反常磁矩的三级数字不能直接相加成"理论值"去和实验硬碰硬——完整理论值还包含四圈 QED、轻子真空极化（μ 子圈与强子圈）、以及弱相互作用的贡献。本节的 4.6e-8 残差正是这些未列入项的量级指示。

💡 关键直觉：圈图是"虚粒子的营业额"。每一个可观测物理量都被全部允许的虚过程折算过一遍；QED 的神奇之处是这套折算可以逐级算清、发散可以记账、残差可以与实验对到小数点后九位。剧场的账房由此挂牌营业。

本节要点回顾

• 圈的定义：封闭内线带来 ∫d⁴k，积分上限处的发散是理论的"短波长账单"
• 两型发散：对数型随 ln Λ 缓涨、二次型随 Λ² 暴涨（数值积分直见）
• 施温格单圈：a_e = α/2π = 0.0011614，抓到实验偏差的 99.85%
• 三级对账：1.161409733e-3 − 1.772302e-6 + 1.480411e-8 = 1.159652235e-3，与实验 1.159652181e-3 会合
• 红外另案：软光子发散在可观测量层面相消，红外安全的思想由此立档

QED 首演圆满，账单如约而至。第五幕的重整化账房开张：先验伤（发散从哪里来），再立账（正规化与重整化手续），最后翻开跑动的账本（重整化群与渐近自由）。