2.3 产生与湮灭算符：演员的上下台动作
2.3 产生与湮灭算符：演员的上下台动作
本节摘要：产生算符 a† 与湮灭算符 a 是量子场论的全部动词系统：a† 把模式从第 n 级激发推到第 n+1 级（粒子上台），a 反向（粒子谢幕）。本节验证算符代数与归一化，定义粒子数算符，并把场算符 φ(x) 本身解读为"在 x 点产生一个粒子的加权探针"。至此"粒子是场的激发"完成数学落地。
承接 2.2 的谐振琴阵列：阶梯算符已经装好，本节给它们发工牌并定义"粒子"；通往 2.4 的福克空间——所有占据数状态的正式登记处。
工牌与动作规范
2.2 的矩阵已经写出了动作：a 的第 (n−1, n) 号元素是 √n。把它翻译成剧场语言：
a†|n⟩ = √(n+1) |n+1⟩ 产生一个粒子（上台） a |n⟩ = √n |n−1⟩ 湮灭一个粒子（谢幕） a |0⟩ = 0 真空无法再空（谢幕无门） 粒子数算符 N = a†a, N|n⟩ = n|n⟩ 对易关系 [N, a†] = a+, [N, a] = -a （计数器与动作的标准换位）
√n 因子不是装饰：它保证归一化。在截断矩阵上把这件事连同矩阵元一起验掉：
# 演算 1：占据数的矩阵元与归一化 import math N = 5 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(N)) for j in range(N)] for i in range(N)] aa = mmul(a, a); adad = mmul(ad, ad) print("⟨0| a a+ |0⟩ =", round(aa[0][0].real, 4), " （a|0⟩=0 之外的部分=1，单粒子态归一）") print("⟨0| a a a+ a+ |0⟩ =", round(mmul(aa, adad)[0][0].real, 4), " （(a+)²|0⟩ 的模方，即 √2! 的平方）") # 输出: # ⟨0| a a+ |0⟩ = 1.0 # ⟨0| a a a+ a+ |0⟩ = 2.0
第二个矩阵元等于 2，意味着两粒子态的归一化系数是 1/√2!。推广开：n 个全同玻色子的态是 (a†)ⁿ|0⟩/√(n!)——对称化自动完成，不需要手工排列组合。单粒子剧本里"多粒子波函数要对称化"的繁琐规则，在场论里是算符代数的免费赠品。
多模式与场算符
真实剧场有无限多架琴（每个 k 一架），通用态写成 |n_k1, n_k2, …⟩。产生算符也按模式分工：a†_k 只动第 k 架。把所有琴按位置重新编组，就得到场算符——它不是新对象，而是产生湮灭算符的位置空间打包：
phi(x) = ∫ d^3k/(2π)³ 1/√(2ω_k) [ a_k e^{-ikx} + a_k+ e^{+ikx} ]
解读这行打包式：场算符在 x 点作用到一个态上，会以 1/√(2ω) 的权重产生（或湮灭）各种动量的粒子。它就是"在 x 点探测/制造一次激发"的探针。第四幕的相互作用项（如 φ⁴）全部由这个探针组装，费曼图的每条线都是一个 φ——分镜脚本的原材料在此出厂。
两点函数给出第一份定量账目：
# 演算 2：真空两点函数 ⟨0|phi(x)phi(0)|0⟩（一维模式求和，m=1） import math def G(xx, L=40.0, nk=8000): s = 0.0; dk = 2*L/nk for i in range(nk): k = -L + i*dk s += math.cos(k*xx)/(2*math.sqrt(k*k+1))*dk/(2*math.pi) return s prev = None for x in (0.0, 1.0, 2.0, 3.0): g = G(x) r = "" if prev is None else " 相邻比值=%.4f（大 x 时趋于 e^-1=0.3679）" % (g/prev) print("G(%.0f) = %.6f%s" % (x, g, r)) prev = g # 输出: # G(0) = 0.697446 # G(1) = 0.070034 相邻比值=0.1004（大 x 时趋于 e^-1=0.3679） # G(2) = 0.016154 相邻比值=0.2307（大 x 时趋于 e^-1=0.3679） # G(3) = 0.006290 相邻比值=0.3894（大 x 时趋于 e^-1=0.3679）
这条曲线说三件事。其一，两点函数只依赖间隔 x−y（平移不变性被模式积分自动编码，平移对称性兑换的能量守恒隐在其中）。其二，它随距离指数式衰减，衰减率由质量 m 控制——重粒子是"尺寸小"的激发，轻粒子更弥散，光子（m=0）则长程相关（这就是电磁力平方反比的场论起源）。其三，G(0) 非零且实际是发散的（被数值截断压到 0.697）：真空里任意两点的场都在关联涨落，"空"是关联意义上的空，不是取值为零。
⚠️ 易错点：把 φ(x) 当作"粒子在 x 的波函数"。φ(x) 是算符，作用在任何态上都会改变粒子数；单粒子波函数则是 ⟨0|φ(x)|1⟩ 这样的矩阵元。混淆这两者是把场论读回单粒子剧本的老毛病。
💡 关键直觉：a† 不创造"某个位置上的粒子"，它创造"某个动量模式上的一份激发"。位置说法只是所有动量模式的相干叠加。粒子的位置从来是衍生概念——在相对论性理论里，这是被迫的严谨。
相干态：上台动作的经典极限
产生湮灭算符最常见的合作演出是相干态 |α⟩ = e^{−|α|²/2} Σ αⁿ/√(n!) |n⟩——激光的量子态原型。它满足 a|α⟩ = α|α⟩（湮灭算符的本征态：谢幕一次，名册形状不变，只换系数），占据数服从泊松分布。在截断名册上验证这两条性质：
# 演算 3：相干态（alpha=2）的泊松画像与本征态验证 import math alpha = 2.0 Nmax = 40 cn = [math.exp(-abs(alpha)**2/2)*alpha**n/math.sqrt(math.factorial(n)) for n in range(Nmax+1)] norm = sum(c*c for c in cn) expN = sum(n*c*c for n, c in enumerate(cn)) expN2 = sum(n*n*c*c for n, c in enumerate(cn)) print("截断 40 阶：归一化 = %.10f，<N> = %.6f，ΔN² = %.6f" % (norm, expN, expN2 - expN**2)) for n in (0, 2, 4, 6, 8): print("P(%d) = %.6f" % (n, cn[n]**2)) dev = 0.0 for n in range(1, Nmax+1): dev = max(dev, abs(math.sqrt(n)*cn[n] - alpha*cn[n-1])) print("a|alpha> = alpha|alpha> 的最大系数偏差 = %.2e" % dev) # 输出: # 截断 40 阶：归一化 = 1.0000000000，<N> = 4.000000，ΔN² = 4.000000 # P(0) = 0.018316 # P(2) = 0.146525 # P(4) = 0.195367 # P(6) = 0.104196 # P(8) = 0.029770 # a|alpha> = alpha|alpha> 的最大系数偏差 = 1.11e-16
三个数字都是教材结论的数值替身：均值 4 = |α|²，方差 4 = 均值——泊松分布的指纹，ΔN = √⟨N⟩ 意味着相对涨落 1/√⟨N⟩ 随强度增大而消失，这正是"经典电磁波"在名册上的出处；本征方程的最大系数偏差停留在浮点噪声级（10 的负 16 次方），√n 因子的阶梯被 α 的指数压得严丝合缝。经典的光不是没有量子结构，而是占据数大到上台与谢幕动作频繁到无法分辨——剧场坐满时，观众看到的是灯光的整体明暗，而不是一次次起身与落座。
本节要点回顾
- 动词系统：a† 上台、a 谢幕、a|0⟩=0（真空无戏可下），√n 因子保归一化
- 计数器：N = a†a 本征值取非负整数，[N,a±] 换位关系标准
- 免费对称化：(a†)ⁿ|0⟩/√(n!) 自动给出玻色对称多体态
- 场即探针：φ(x) 是产生湮灭算符的位置打包，第四幕相互作用的全部原料
- 关联尺寸：两点函数指数衰减，衰减长度由质量倒数（康普顿波长）设定
单架琴的账目理清，下一节把所有琴的花名册装订成册——福克空间，并第一次认真打量那个"空"的座位：真空。