2.2 正则量子化：给无限架谐振琴调音

2.2 正则量子化：给无限架谐振琴调音

本节摘要：正则量子化把场 φ 与共轭动量 π 升级为算符并强加对易关系 [φ(x),π(y)] = iδ(x−y)；模式展开使场等价于无限多架独立量子谐振琴，每架的能级为 (n+1/2)ω。本节用截断矩阵亲手验证对易关系与能级结构，再把全部零点能加总——真空能量密度按截断四次方增长，第一张无穷大账单正式入账。

承接 2.1 的模式清单：每个 k 模式已是独立的简正坐标；本节给它们装上量子化，通往 2.3 的产生湮灭算符——那些算符只是本节谐振子阶梯算符的更名。

量子化的三步操作

正则量子化的手续与量子力学里处理粒子完全同构，只是"位置与动量"换成了"场与场的共轭动量"：

第一步：从拉格朗日密度求共轭动量密度
  pi(x,t) = dL/d(d_t phi) = d_t phi（自由标量场的情形）

第二步：升级为算符，强加等时对易关系
  [phi(x,t), pi(y,t)] = i δ(y-x)      （δ 为三维狄拉克函数）
  [phi(x,t), phi(y,t)] = 0
  [pi(x,t), pi(y,t)] = 0

第三步：把 2.1 的模式展开代入，对易关系转化为模式算符的关系
  [a_k, a_k'+] = δ(k-k')  ——每架琴一套阶梯算符
  哈密顿量 H = ∫ d^3k ω(k) (a_k+ a_k + 1/2)
​

第三步的代数值得在草稿纸上走一遍：φ 含 a 与 a† 两支，π 含 −iωa 与 +iωa†；把 [φ,π] 的傅里叶分量逐一配对，δ 函数正好来自平面波的正交归一。走完这段，你会确信"场的量子化"没有新物理假设——只有一次换了对象的正则替换。谐振子的全部知识随即免税过户：能级 (n+½)hbar ω、阶梯算符、基态 |0⟩。

截断矩阵上的调音实验

对易关系 [a, a†] = 1 在有限维矩阵里只能近似实现——这个"近似失败的位置"本身就是教材。取 5 维截断：

# 演算 1：截断福克空间中的阶梯算符与对易关系
import math
N = 5                                  # 保留 |0> 到 |4>
a = [[0j]*N for _ in range(N)]
for n in range(1, N):
    a[n-1][n] = math.sqrt(n)           # a|n> = sqrt(n)|n-1>
ad = [[a[j][i] for j in range(N)] for i in range(N)]   # a+
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)]
comm = mmul(a, ad); comm2 = mmul(ad, a)
print("[a,a+] 对角元:", [round((comm[i][i]-comm2[i][i]).real, 4) for i in range(N)])
print("N=a+a 对角元:", [round(comm2[i][i].real, 4) for i in range(N)])
# 输出:
# [a,a+] 对角元: [1.0, 1.0, 1.0, 1.0, -4.0]
# N=a+a 对角元: [0.0, 1.0, 2.0, 3.0, 4.0]
​

前四个对角元严格等于 1，最后一个是 −4：截断砍掉了 |5⟩，于是 a† 在最高层"无处可去"，对易关系在边界破产。教训朴素而重要——任何有限截断都是对理论的近似，真正的对易关系活在无限维。而 N = a†a 的对角元 0,1,2,3,4 则预告了粒子数算符：这架琴的能级恰好以 ω 为台阶计数。

图 无限谐振琴阵列：场的量子化全景

零点能的加总：第一张无穷大账单

每架琴贡献 ½ω(k)，全剧场加总：

# 演算 2：真空能量密度对截断的依赖（m=1，自然单位制）
import math
def rho_vac(Lam, ngrid=4000):
    # (1/(2π)^3) ∫ d^3k 0.5*√(k²+m²)，球对称化为径向辛普森积分
    h = Lam/ngrid
    s = 0.0
    for i in range(ngrid+1):
        kk = i*h
        f = kk*kk*0.5*math.sqrt(kk*kk+1)
        w = 1 if i in (0, ngrid) else (4 if i % 2 == 1 else 2)
        s += w*f
    return s*(h/3)*4*math.pi/(2*math.pi)**3
for Lam in (2.0, 4.0, 8.0):
    r = rho_vac(Lam)
    print("Λ=%4.0f GeV -> ρ_vac = %9.4f GeV^4 ; 16π²ρ/Λ⁴ = %.4f（趋于 1）" % (Lam, r, 16*math.pi**2*r/Lam**4))
# 输出:
# Λ=   2 GeV -> ρ_vac =    0.1229 GeV^4 ; 16π²ρ/Λ⁴ = 1.2127（趋于 1）
# Λ=   4 GeV -> ρ_vac =    1.7166 GeV^4 ; 16π²ρ/Λ⁴ = 1.0589（趋于 1）
# Λ=   8 GeV -> ρ_vac =   26.3355 GeV^4 ; 16π²ρ/Λ⁴ = 1.0153（趋于 1）
​

截断每翻一倍，能量密度涨约十六倍——严格的四次方增长，渐近式 ρ → Λ⁴/16π²（修正来自 m² 项，随 Λ 增大而稀释）。这张账单有两种态度。物理学的处理（第六幕）：零点能的绝对值不可观测，可观的是差值，重整化把它吸收进宇宙学常数——只留一个至今未解的"为什么宇宙学常数这么小"的疑问（著名的 120 个数量级问题）。实验的处理（2.4 节）：卡西米尔效应直接测到了零点能的差值。两种态度都承认同一事实：真空是满座的，只是票价只在改动座位表时才显形。

💡 关键直觉：正则量子化没有引入任何"粒子"概念——粒子是从谐振琴能级的等距台阶里自动长出来的。台阶高度 ω(k) 正是模式能量，"一份激发 = 一个粒子"是记账的自然结果，不是额外公理。

本节要点回顾

• 三步手续：共轭动量、算符化、对易关系 [φ,π] = iδ，与单粒子量子化同构
• 模式对易：[a_k, a_k'†] = δ(k−k')，每架琴独立拥有阶梯算符
• 截断教训：5 维矩阵里 [a,a†] 对角为 1,1,1,1,−4，边界破坏对易关系，真实结构在无限维
• 能级计数：H = Σω(n+½)，台阶就是粒子数（N 对角元 0,1,2,3,4）
• 零点账单：ρ_vac 随截断四次方增长，渐近 Λ⁴/16π²；可观的是差值而非绝对值

琴已调好，下一节正式介绍两位新员工：a† 负责敲响，a 负责收音。第一个粒子将在矩阵的第 (2,0) 号元素里诞生。