4.2 戴森级数与威克定理：剧本的展开与理账

4.2 戴森级数与威克定理：剧本的展开与理账

本节摘要：戴森级数把时序演化算符按耦合常数展开为幂级数，第 n 阶对应对相互作用顶点的 n 重时空积分；威克定理把每阶中乱序的场算符乘积整理为"全收缩乘正常序"之和，而收缩恰是费曼传播子。本节用数值枚举验证收缩计数 (n−1)!!（两场 1、四场 3、六场 15、八场 105），并讨论组合爆炸对微扰论适用范围的界定。

承接 4.1 的时序演化：本节把它逐阶展开并整理成可计算的形式；通往 4.3 的费曼图——图形化正是这套理账术的可视化界面。

逐阶展开的账本

戴森级数的形式骨架很朴素，就是"指数函数的级数"加上时序记号：

U_I(t,0) = T exp[-i ∫₀ᵗ H_int(t')dt']
         = 1 + (-i)∫ H_int(t₁)dt₁ + (-i)²/2! ∫∫ T H(t₁)H(t₂) + ...

以 phi⁴ 理论为例（L_int = -λ/4! φ⁴，H_int = ∫ d³x λ/4! φ⁴）：
  一阶项：单个顶点，四条 φ 线从顶点长出
  二阶项：两个顶点、八条线，时序积分遍及全部先后排法
  每多一阶：多一个顶点、多一个耦合常数、多一层积分
​

级数的每一项都是"若干时空点上各做一次相互作用"的历史分账。λ 小的时候，阶数越高贡献越小，级数快速收敛——这正是"微扰"二字的含义。QED 的耦合 α ≈ 1/137，收敛极快；但第六幕会看到，高阶项里藏着圈积分的发散，收敛性要与重整化联手才能兑现成与实验的对账。

威克定理：把乱账理成乘积

级数的麻烦在于场算符不会自动排好队：第 n 阶出现 4n 个 φ 的乘积（φ⁴ 理论），它们分属不同时空点、互不对易。威克 1950 年的定理解决了理账问题：

威克定理：任意多个场的乘积
  = 正常序（全部湮灭算符右移）
  + 单收缩 × 正常序
  + 双收缩 × 正常序 + ...
  + 全部场两两收缩配对（全收缩项）

收缩的定义：两点的收缩 = 时序两点函数 = 传播子
  contract(phi(x) phi(y)) = ⟨0|T phi(x)phi(y)|0⟩ = Delta_F(x-y)
​

定理的威力在于把"算符的乱序多项式"换算成"传播子的乘积"——前者无法直接作用于名册，后者是普通的复数函数。当外线数目固定（比如 2→2 散射，四个外场），只有把其余场全部收缩掉的项才对振幅有贡献，于是计算化为数图与积分。收缩的配对数有精确公式，数值枚举一次：

# 演算 1：n 个场的全收缩配对数枚举（威克定理的组合学）
def pairings(items):
    if not items:
        yield []
        return
    first = items[0]
    for i in range(1, len(items)):
        rest = items[1:i] + items[i+1:]
        for p in pairings(rest):
            yield [(first, items[i])] + p
import math
for n in (2, 4, 6, 8):
    cnt = sum(1 for _ in pairings(list(range(n))))
    print("n=%d 个场 -> %3d 种全收缩   （公式 (n-1)!! = %d：%s）"
          % (n, cnt, math.prod(range(1, n, 2)), cnt == math.prod(range(1, n, 2))))
# 输出:
# n=2 个场 ->   1 种全收缩   （公式 (n-1)!! = 1：True）
# n=4 个场 ->   3 种全收缩   （公式 (n-1)!! = 3：True）
# n=6 个场 ->  15 种全收缩   （公式 (n-1)!! = 15：True）
# n=8 个场 -> 105 种全收缩   （公式 (n-1)!! = 105：True）
​

1、3、15、105——双阶乘序列。注意增长速度：从八场到十二场（三个顶点的真空图），配对数跳到 10395。组合爆炸是微扰论的双刃剑：每一阶的图数比耦合常数的衰减涨得更快，级数实际上是渐近的（对有限耦合最终发散）；但在 α ~ 1/137 的领域，前三四阶已经给出天文精度，第六幕的反常磁矩就是这个幸福的极致案例。

图 威克定理的理账流程

正常序与真空图

威克定理里的"正常序"（湮灭算符全部右移）有一个重要的副产品：如果一项里留下未收缩的场对，它作用在真空上会给出真空期待值——这就是真空图（不连任何外线的闭合圈）。它们贡献的是 S 矩阵的整体相位（真空到真空的账），对一切散射截面的干涉结构无影响，常规计算里直接丢弃；但它们并非无意义——第六幕会看到，真空图的求和正是"真空能量"的重整化对象，2.4 的零点账单在圈图层面的化身。

# 演算 2：正常序的真空期望消失（截断矩阵验证）
import math
N = 6
a = [[0j]*N for _ in range(N)]
for n in range(1, N):
    a[n-1][n] = math.sqrt(n)
ad = [[a[j][i] for j in range(N)] for i in range(N)]
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))]
e0 = [1, 0, 0, 0, 0, 0]                     # |0>
def act_v(M, v):
    return [sum(M[i][k]*v[k] for k in range(len(v))) for i in range(len(M))]
print("⟨0|phi_正常序|0⟩ 型： a|0> = ", [round(v.real, 6) for v in act_v(a, e0)], "（零向量）")
print("⟨0| a+ a |0⟩ = ", round(sum(e0[i]*act_v(mmul(ad, a), e0)[i] for i in range(N)).real, 6), "（真空无粒子可数）")
print("⟨0| a a+ |0⟩ = ", round(sum(e0[i]*act_v(mmul(a, ad), e0)[i] for i in range(N)).real, 6), "（对易子的真空账 = 1）")
# 输出:
# ⟨0|phi_正常序|0⟩ 型： a|0> =  [0, 0, 0, 0, 0, 0] （零向量）
# ⟨0| a+ a |0⟩ =  0.0 （真空无粒子可数）
# ⟨0| a a+ |0⟩ =  1.0 （对易子的真空账 = 1）
​

三行矩阵元把正常序的账目纪律说尽：正常序的算符串在真空上取零；只有"违规"留下未正常序的对（如 a a†）才贡献 1——这一份正是收缩，回到传播子。真空图的系统性出现与丢弃规则由此有据可依。

⚠️ 易错点：威克定理只对自由场（高斯理论）成立。相互作用一旦打开，场不再自由，定理不能直接套用——微扰论的逻辑恰恰是"相互作用只出现在顶点上，传播仍是自由场的"，4.1 的绘景分工保证了这个用法合法。

💡 关键直觉：戴森级数加威克定理的组合，本质是把一个不可解的多体问题换成"数所有可能历史分账，逐项求和"。路径积分（7.1）将提供同一件事的另一种记账法——不拆顶点、直接对所有场构型积分——两种记法在自由场与微扰层面给出完全相同的结果。

本节要点回顾

• 戴森级数：时序指数展开为耦合常数的幂级数，每阶多一个顶点与一层积分
• 威克定理：乱序乘积化为"收缩 × 正常序"之和，收缩即传播子
• 组合计数：全收缩数 (n−1)!!（1、3、15、105 逐一验证），组合爆炸划定渐近边界
• 正常序纪律：正常序算符真空期待为零，未配对的部分即真空图账目
• 真空图去向：贡献 S 的整体相位，截面计算丢弃，但成为第六幕真空能量重整化的对象

理账术完备，下一节启用图形界面：把配对画成线、顶点画成交点，费曼规则表把每张图机械翻译成公式——第一张全流程对账单（φ⁴ 树图截面）随之出炉。