1.2 场的技术图纸与作用量：拉格朗日密度

1.2 场的技术图纸与作用量：拉格朗日密度

本节摘要：场的动力学由拉格朗日密度编码，作用量对场构型取极值给出欧拉-拉格朗日场方程。本节从最小作用量原理推出克莱因-戈登方程，用闵可夫斯基内积验证四动量的不变性，并完成最简场构型的数值求解检查。这是第二幕量子化操作的全部图纸。

承接 1.1 的验伤报告：旧剧本废弃后，新实体"场"需要一个标准的记谱法；本节铺完图纸，通往 1.3 的守恒账本——作用量的对称性将在那里兑换成守恒量。

图纸的第一张：场与作用量

场是给时空每一点配一个物理量的对象。标量场 φ(x) 在每点给一个数（温度场、希格斯场的经典对应），矢量场给一个矢量（电磁势 A），旋量场给一个带自旋指标的旋量（第三幕的电子场）。剧场隐喻在这里很贴切：场是遍布整个舞台的"介质"，而拉格朗日密度就是这部剧的"总谱"——它写明每一处时空、每一种场、以及它们的导数如何组合出被积分的量。

作用量 S 是拉格朗日密度对全时空的积分。最小作用量原理说：真实演出的场构型，使 S 取驻值。对场做小变分 δφ，要求 δS 为零，分部积分把导数从 δφ 上挪走，就得到欧拉-拉格朗日场方程。以自由实标量场为例，总谱只有两项：

拉格朗日密度  L = 1/2 * (d_phi/dx)^mu (d_phi/dx)_mu  -  1/2 * m^2 * phi^2
              动能项：场的时空梯度        质量项：每个位置为非零场值付出的代价

作用量        S = ∫ d^4x  L

变分步骤：
  dS/d(phi) 逐点为零
  -> 对梯度和时间导数分部积分（边界项取零）
  -> (Box + m^2) phi = 0     即克莱因-戈登方程
  其中 Box = d_t^2 - d_x^2 - d_y^2 - d_z^2（达朗贝尔算符）
​

注意质量项的符号约定：动能项与质量项相减，是保证平面波解的频率满足 ω² = k² + m²（自然单位制）的唯一选择。这个色散关系就是第二幕谐振子分解的入口。

图纸的第二张：闵可夫斯基的记账法

相对论的账本用度规 (+,−,−,−)：两个四矢量的内积是时间分量相乘减去空间分量点乘。动能项里的 (dφ)^μ(dφ)_μ 正是用这套记账法压缩的。用数值把不变性落到实处：

# 演算 1：闵可夫斯基内积与质壳条件（度规 +,-,-,-）
E, px, py, pz = 5.0, 4.0, 0.0, 0.0      # 单位 GeV，自然单位制
p2 = E*E - px*px - py*py - pz*pz
print("p·p =", p2, "GeV^2 ；若粒子质量 m=3 GeV，则 p·p=m^2 =", p2 == 9.0)
# 换一个惯性系（沿 x 轴 boost，速度 0.6c，gamma=1.25）：
g, v = 1.25, 0.6
E_p  = g*(E - v*px)      # 洛伦兹变换后的能量
px_p = g*(px - v*E)      # 变换后的动量
p2_p = E_p*E_p - px_p*px_p
print("变换后: E'=%.4f, p_x'=%.4f, p'·p'=%.4f" % (E_p, px_p, p2_p))
# 输出:
# p·p = 9.0 GeV^2 ；若粒子质量 m=3 GeV，则 p·p=m^2 = True
# 变换后: E'=1.2500, p_x'=-2.7500, p'·p'=9.0000
​

同一笔账在两个惯性系里分毫不差——这就是"质壳条件"p·p = m² 的含义。拉格朗日密度里的每一项都按这套记账法书写，作用量自然是洛伦兹不变的，场方程也就自动继承了相对论。1.1 的穿帮二（负能解）在图纸层面依然存在（E = ±√(p²+m²)），但它的最终解法要等量子化把"概率密度"这个角色撤掉。

把方程放到检验台上

图纸是否正确，最便宜的办法是拿一个平面波解代入方程，看残差是否为零。

# 演算 2：克莱因-戈登方程的平面波解数值检查（m=1, k=2, 自然单位制）
import math
m, k = 1.0, 2.0
w = math.sqrt(k*k + m*m)                # 色散关系 w = sqrt(k^2+m^2)
x, t = 0.7, 1.3
phi = math.cos(k*x - w*t)               # 试解 phi = cos(kx - wt)
residual = (-w*w + k*k + m*m)*phi       # (d_t^2 - d_x^2 + m^2)phi 的解析结果
print("w = %.6f,  方程残差 = %.2e" % (w, residual))
w_wrong = 1.01*w                        # 故意用错的频率
residual_wrong = (-w_wrong*w_wrong + k*k + m*m)*math.cos(k*x - w_wrong*t)
print("频率偏差 1%% 时残差 = %.6f（不为零，色散关系是刚性约束）" % residual_wrong)
# 输出:
# w = 2.236068,  方程残差 = -5.67e-17   （浮点噪声级别，试解成立）
# 频率偏差 1% 时残差 = -0.003501（不为零，色散关系是刚性约束）
​

残差在 10 的负 17 次方量级，纯属性能优良的浮点噪声；而频率偏离色散关系百分之一，残差立刻可见。这张图纸是刚性的：任何求解、任何量子化都必须尊重 ω² = k² + m²。2.1 节将专门解剖这条色散曲线，它是"质量是频率的下限"这一剧场规则的全部来源。

⚠️ 易错点：初学者常把"拉格朗日密度"与"拉格朗日量"混用。密度是逐点的函数，拉格朗日量是密度对空间积分；作用量再对时间积分。量子场论的惯用对象是密度，因为它显式洛伦兹协变。

💡 关键直觉：作用量原理把"运动方程"从公理降级为定理。你只需选择总谱（拉格朗日密度），演出规则自动涌现——第二幕加相互作用（第四幕）、换规范场（终幕）时，改的从来只是总谱上几行字。

给总谱加一行字：微扰如何进入图纸

图纸的价值在于"改一行字就能换一部剧"。在自由标量总谱上加相互作用项 −λφ⁴/4!（第四幕 φ⁴ 理论的同款），变分手续原封不动重跑一遍：新项对场方程的贡献是 −λφ³/6，自由方程 (□ + m²)φ = 0 随之变成 (□ + m²)φ = −λφ³/6。用数值看线性解在这行新字下的处境：

# 演算 3：加 -λφ⁴/4! 后，线性平面波解不再精确（m=1, k=2, λ=0.01）
import math
k, m, lam, amp = 2.0, 1.0, 0.01, 1.0
w = math.sqrt(k*k + m*m)
maxR = 0.0
xs = [i*0.05 for i in range(int(2*math.pi/k/0.05)+1)]
for x in xs:
    for t in (0.0, 0.37, 0.83, 1.29):
        phi = amp*math.cos(k*x - w*t)
        R = lam*phi**3/6.0          # 线性解代入修正方程后的残差
        maxR = max(maxR, abs(R))
print("w = %.4f；线性解代入修正方程的最大残差 = %.4e" % (w, maxR))
print("相对幅度 maxR/amp = %.4e = lam*amp^2/6 = %.4e" % (maxR, lam*amp**2/6))
# 输出:
# w = 2.2361；线性解代入修正方程的最大残差 = 1.6667e-03
# 相对幅度 maxR/amp = 1.6667e-03 = lam*amp^2/6 = 1.6667e-03
​

残差处处不为零，但很小：弱耦合（λ = 0.01）时，自由场的平面波仍是"近似合格的戏"，只是场值处处带着 λ·amp²/6 量级的账目误差。这个比例正是微扰论的入场宣言——第四幕以后的所有计算都建立在这句承诺上：耦合越弱，自由剧目离真实演出越近。而量子化手续对总谱的任何一行字一视同仁：改的只是图纸，量的仍是同一套场坐标与共轭动量。这就是为什么场论教材把"写对拉格朗日密度"当作唯一的创造性步骤——剩下的都是机械执行。

本节要点回顾

• 场与总谱：拉格朗日密度 L(φ, dφ) 是场动力学的完整编码
• 作用量驻值：变分加分部积分给出欧拉-拉格朗日场方程，KG 方程是自由标量场的结果
• 色散关系：ω² = k² + m² 由质量项符号决定，是后续一切谐振子分解的锚点
• 闵可夫斯基记账：(+,−,−,−) 度规下 p·p = m² 在惯性系间不变，作用量自动协变
• 图纸的刚性：平面波试解的残差检查显示色散关系不容偏离

图纸有了，下一站立账本：作用量的每一种连续对称性，都能在诺特定理的窗口兑换出一份守恒量。