4.1 相互作用绘景与S矩阵：换导演不换舞台

4.1 相互作用绘景与S矩阵：换导演不换舞台

本节摘要：相互作用绘景让场算符按自由场演化、态矢量承担全部相互作用，使微扰展开可行；S 矩阵定义为时间演化算符在无穷远前后极限下的值，记录入射名册到出射名册的转换幅。本节用二能级玩具模型数值演示时序乘积与朴素指数的稳定差异（0.032 与 0.064），为时序算符 T 的必要性提供物理认证。

承接第三幕的演员阵容：三类场与传播子全部就位，本节第一次让他们同台；通往 4.2 的戴森级数——那里这套时序演化将被逐阶展开。

为什么还要第三种绘景

量子力学的"绘景"是分工协议。薛定谔绘景里态矢量干全部活、算符休息；海森堡绘景反过来。相互作用绘景是折中方案，专为微扰论定制：

分解：      H = H_0 + H_int          （自由部分 + 相互作用部分）
绘景约定：  算符按 H_0 演化（第二、三幕的全部成果继续可用！）
            态矢量按 H_int 演化（只有这部分需要展开处理）

态矢量方程：  i d/dt |ψ_I⟩ = H_int(t) |ψ_I⟩
             其中 H_int(t) 已穿上"自由演化戏服"：
             H_int(t) = e^{iH_0 t} H_int e^{-iH_0 t}

S 矩阵：    S = lim(t→±∞) U_I(t, t') = T exp[-i ∫ H_int(t') dt']
            入射名册 |in⟩ 经 S 得到出射名册：|out⟩ = S|in⟩
​

这份协议的聪明之处在于人力资源的复用：对易关系、模式展开、传播子（第二、三幕的全部机器）都属于自由场，在相互作用绘景里原封不动地继续工作；而相互作用被压缩进态矢量的演化因子里，那里恰好适合做级数展开。实验语言与 S 矩阵也无缝对接：对撞机把探测器读到的出射粒子按名册登记，截面就是 |⟨out|S|in⟩|² 的统计汇总（第四幕开头细算）。

时序算符 T 在 S 的定义式里出场（把晚发生的算符移到左侧），它不是排版偏好。用玩具模型做对照实验：

# 演算 1：时序排序的物理后果（二能级玩具，H(t) = t·sx + t²·sz）
import math
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))]
sx = [[0, 1], [1, 0]]; sz = [[1, 0], [0, -1]]
def expm2(A):
    U = [[1 if i == j else 0 for j in range(2)] for i in range(2)]
    P = [[1 if i == j else 0 for j in range(2)] for i in range(2)]
    for k in range(1, 14):
        P = mmul(P, A)
        for i in range(2):
            for j in range(2):
                U[i][j] += P[i][j]/math.factorial(k)
    return U
def run(N, T=1.0):
    dt = T/N
    Ut = [[1, 0], [0, 1]]            # 时序：晚的左乘
    Ua = [[1, 0], [0, 1]]            # 反时序
    for n in range(N):
        t = (n+0.5)*dt
        H = [[t*sx[i][j] + t*t*sz[i][j] for j in range(2)] for i in range(2)]
        S = expm2([[-1j*H[i][j]*dt for j in range(2)] for i in range(2)])
        Ua = mmul(Ua, S)
        Ut = mmul(S, Ut)
    Havg = [[0.5*sx[i][j] + (1/3)*sz[i][j] for j in range(2)] for i in range(2)]
    Un = expm2([[-1j*Havg[i][j] for j in range(2)] for i in range(2)], )
    return Ut, Ua, Un
for N in (5, 20, 100, 500):
    Ut, Ua, Un = run(N)
    d = lambda A, B: max(abs(A[i][j]-B[i][j]) for i in range(2) for j in range(2))
    print("N=%4d 步: |U_时序 - U_朴素指数| = %.6f ; |U_时序 - U_反时序| = %.6f"
          % (N, d(Ut, Un), d(Ut, Ua)))
# 输出:
# N=   5 步: |U_时序 - U_朴素指数| = 0.030531 ; |U_时序 - U_反时序| = 0.060969
# N=  20 步: |U_时序 - U_朴素指数| = 0.031995 ; |U_时序 - U_反时序| = 0.063890
# N= 100 步: |U_时序 - U_朴素指数| = 0.032089 ; |U_时序 - U_反时序| = 0.064077
# N= 500 步: |U_时序 - U_朴素指数| = 0.032092 ; |U_时序 - U_反时序| = 0.064085
​

三个信息值得逐条读。第一，时序与反时序的差稳定在 0.064，不随步数加密而消失——它是连续极限下的真实差异（马格努斯展开的交换子修正），不是离散化误差。第二，时序乘积与"先积分再取指数"的朴素做法差 0.032，同样稳定——当 H 在不同时刻不对易时，exp(−i∫H dt) 根本不是正确的演化算符。第三，加密步数只让结果更干净（N=20 与 N=500 一致到第三位），数值流程可信。物理源流很清晰：不同时刻的相互作用不对易（先后做两次扰动，顺序影响结果），所以"按时间顺序逐段乘起来"才是唯一合法的记账。

绝热开关：让 S 有定义

S 矩阵取无穷远时刻的极限，但真实的相互作用从不关闭。微扰论的正规处理是"绝热开关"：想象相互作用的耦合常数在遥远过去与遥远未来缓慢趋于零，中间恢复原值。这个数学道具带来两点好处：入射与出射名册由自由场定义清楚（第二、三幕的占据数态）；中间过程的全部修正由 S 吸收。物理学界的共识是：只要不问开关本身的问题（不做涉及无穷长时间累积的极端精密测量），这个手续不改变任何物理预言。第六幕会看到它留下的一处著名的"需要小心"的红外问题，但那是另一幕的剧情。

S 矩阵元的操作含义：
  ⟨f| S |i⟩ = ⟨出射名册 f | 入射名册 i 经全部相互作用后的转换幅⟩
  实验可观测量：|⟨f|S|i⟩|² 配上相空间与流强 -> 截面（5.1 正式定义）

S 的分解：S = 1 + iT    （单位算符 + 真正散射的部分 T）
         1 对应"演员直接穿过舞台没演对手戏"的份额
​

⚠️ 易错点：把 S 矩阵理解为"描述某个瞬间的散射"。它是连接两个无穷远时刻的整本账，中间发生的所有事（包括第六幕的全部圈修正）都已包含在内。初学者常在画图时把 S 元素与某个顶点混同，那是把总账与单笔交易搞混了。

💡 关键直觉：相互作用绘景的本质是"复用自由场的机器，集中火力打相互作用"。工程上等价于把一个无法直接求解的 H 分成"会解的部分"与"小的部分"，然后对小的部分做级数——这套思路在物理以外的领域同样随处可见。

本节要点回顾

• 三绘景分工：相互作用绘景让算符走自由演化、态矢量扛相互作用，微扰展开可行
• S 矩阵：无穷远极限下的演化算符，入射到出射名册的整本转换账
• 时序认证：非对易 H 下时序与反时序差 0.064、与朴素指数差 0.032，稳定不随步数消失
• 绝热开关：数学道具，保证名册良定义，不改变中间过程的物理
• S = 1 + iT：单位算符是"没演对手戏"的直通份额，T 承载全部散射

演化账本有了，下一节把它展开：戴森级数逐阶写下每一笔交易，威克定理把乱账理成传播子的乘积——费曼图就在那个整理台上诞生。