粒子物理中的散射截面计算
散射截面是粒子物理中连接量子场论计算与实验测量的核心量。它描述入射粒子与靶粒子发生特定相互作用并产生末态的概率,量纲为面积。常用单位包括 barn、mb、pb、fb,其中 (1,\text{barn}=10^{-28},\text{m}^2),且 (1,\text{GeV}^{-2}\approx 0.389,\text{mb})。在电磁相互作用、弱相互作用和强相互作用过程中,截面计算都依赖相同的洛伦兹不变相空间与跃迁振幅框架。
按观测方式,截面可分为:
- 总截面 (\sigma_{\text{tot}}):所有末态之和。
- 微分截面 (d\sigma/d\Omega):给定立体角内的分布。
- 双微分截面 (d^2\sigma/(d\Omega dE)):同时测量角度与能量。
- 排他截面与包含截面:前者指定完整末态,后者只指定部分粒子。
实验上,事件率 (R) 与截面关系为
[
R = \mathcal{L},\sigma,
]
其中 (\mathcal{L}) 是瞬时亮度。理论计算的目标通常是给出 (\sigma) 或 (d\sigma/d\Omega) 的预言。
从跃迁振幅到截面公式
在量子场论中,(S) 矩阵写为
[
S_{fi}=\delta_{fi}+i(2\pi)^4\delta^{(4)}(P_f-P_i)\mathcal{M}{fi}.
]
(\mathcal{M}{fi}) 是费曼规则给出的不变振幅。对 (2\to n) 过程,微分截面为
[
d\sigma=\frac{1}{4E_A E_B |v_A-v_B|}
|\mathcal{M}|^2
(2\pi)^4\delta^{(4)}(P_f-P_i)
\prod_f \frac{d^3p_f}{(2\pi)^3 2E_f}.
]
其中通量因子可写成洛伦兹不变形式:
[
4E_A E_B |v_A-v_B|
=4\sqrt{(p_A\cdot p_B)^2-m_A^2m_B^2}.
]
因此,计算截面的核心步骤是:写出费曼图、得到 (\mathcal{M})、求模方并自旋求和/平均,最后对末态相空间积分。
两体散射与曼德斯坦变量
对 (2\to2) 过程 (A+B\to C+D),常用曼德斯坦变量:
[
s=(p_A+p_B)^2,\quad
t=(p_A-p_C)^2,\quad
u=(p_A-p_D)^2.
]
它们满足
[
s+t+u=m_A^2+m_B^2+m_C^2+m_D^2.
]
在质心系中,若初末态粒子质量可忽略,则
[
\frac{d\sigma}{d\Omega}
=\frac{|\mathcal{M}|^2}{64\pi^2 s}.
]
若末态动量大小不同,还需乘以 (|p_f|/|p_i|)。曼德斯坦变量直接对应能量与动量转移,是电磁场能量动量守恒在散射运动学中的体现。
示例:电子—正电子湮灭为 μ 子对
考虑树图 QED 过程
[
e^+e^-\to \mu^+\mu^-.
]
在质心系高能极限下,忽略电子与 μ 子质量,自旋平均后的振幅模方为
[
\overline{|\mathcal{M}|^2}=e^4(1+\cos^2\theta).
]
代入两体微分截面公式:
[
\frac{d\sigma}{d\Omega}
=\frac{\alpha^2}{4s}(1+\cos^2\theta),
]
其中 (\alpha=e^2/(4\pi))。对立体角积分:
[
\sigma=\int \frac{d\sigma}{d\Omega}d\Omega
=\frac{4\pi\alpha^2}{3s}.
]
若取 (\sqrt{s}=10,\text{GeV}),(s=100,\text{GeV}^2),(\alpha\approx1/137),则
[
\sigma\approx 2.23\times10^{-6},\text{GeV}^{-2}
\approx 0.87,\text{nb}.
]
该结果与实验在树图精度内一致,是检验 QED 的经典过程。
电磁场能量动量的角色
在 QED 中,电磁场四维势 (A_\mu) 与场强张量
[
F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu
]
构成拉氏量
[
\mathcal{L}=-\frac14 F_{\mu\nu}F^{\mu\nu}
+\bar{\psi}(i\gamma^\mu D_\mu-m)\psi.
]
电磁场能量动量张量
[
T^{\mu\nu}=-F^{\mu\rho}F^\nu_{\ \rho}
+\frac14 g^{\mu\nu}F_{\rho\sigma}F^{\rho\sigma}
]
决定了光子传播子、能量守恒与规范不变性。费曼规则中的光子传播子 (-ig_{\mu\nu}/q^2) 正来自电磁场动能项,而 (\delta^{(4)}(P_f-P_i)) 则来自时空平移不变性。因此,散射截面计算本质上是电磁场能量动量在量子层面的重新分配。
数值计算与工具链
现代截面计算通常借助自动化工具完成:
- 用 FeynArts、MadGraph5_aMC@NLO 或 CalcHEP 生成费曼图与振幅。
- 用 FormCalc、Helicity 方法或张量约化计算 (|\mathcal{M}|^2)。
- 用 Vegas、Cuba 等算法进行相空间积分。
- 处理红外与紫外发散,引入重整化与因子化。
- 对强子过程,还需卷积部分子分布函数并接入 Pythia、Sherpa 进行部分子簇射。
掌握从振幅到截面的完整链条,是理解粒子物理实验、电磁场能量动量传递以及新物理搜索的基础。