5.1

Problem 5.1

schwarzChapter 5

习题 5.1

来源: 第5章, PDF第67页


5.1 Show that the differential cross section for 222 \rightarrow 2 scattering with piμ+pAμpfμ+pBμp_i^\mu + p_A^\mu \rightarrow p_f^\mu + p_B^\mu in the rest frame of particle AA can be written as

dσdΩ=164π2mA[EB+Ef(1pipfcosθ)]1pfpiM2,(5.54)\frac{d\sigma}{d\Omega} = \frac{1}{64\pi^2 m_A} \left[ E_B + E_f \left( 1 - \frac{|\vec{p}_i|}{|\vec{p}_f|} \cos\theta \right) \right]^{-1} \frac{|\vec{p}_f|}{|\vec{p}_i|} |\mathcal{M}|^2, \tag{5.54}

where θ\theta is the angle between pi\vec{p}_i and pf\vec{p}_f, EB=(pfpi)2+mB2E_B = \sqrt{(\vec{p}_f - \vec{p}_i)^2 + m_B^2} and Ef=pf2+mf2E_f = \sqrt{\vec{p}_f^2 + m_f^2}.

习题 5.1 - 解答


对于 222 \rightarrow 2 散射过程 piμ+pAμpfμ+pBμp_i^\mu + p_A^\mu \rightarrow p_f^\mu + p_B^\mu,其微分截面的一般公式由相空间积分与不变矩阵元平方 M2|\mathcal{M}|^2 以及入射通量因子构成:

dσ=14(pipA)2mi2mA2M2dΠLIPSd\sigma = \frac{1}{4 \sqrt{(p_i \cdot p_A)^2 - m_i^2 m_A^2}} |\mathcal{M}|^2 d\Pi_{\text{LIPS}}

其中,洛伦兹不变相空间 (Lorentz Invariant Phase Space, LIPS) 为:

dΠLIPS=(2π)4δ(4)(pi+pApfpB)d3pf(2π)32Efd3pB(2π)32EBd\Pi_{\text{LIPS}} = (2\pi)^4 \delta^{(4)}(p_i + p_A - p_f - p_B) \frac{d^3p_f}{(2\pi)^3 2E_f} \frac{d^3p_B}{(2\pi)^3 2E_B}

1. 运动学与通量因子分析 在靶粒子 AA 的静止参考系(实验室系)中,各粒子的四维动量可表示为: pAμ=(mA,0)p_A^\mu = (m_A, \vec{0}) piμ=(Ei,pi)p_i^\mu = (E_i, \vec{p}_i) pfμ=(Ef,pf)p_f^\mu = (E_f, \vec{p}_f) pBμ=(EB,pB)p_B^\mu = (E_B, \vec{p}_B) 首先计算入射通量因子。由于 pipA=EimAp_i \cdot p_A = E_i m_A,通量因子可以化简为:

4(pipA)2mi2mA2=4Ei2mA2mi2mA2=4mAEi2mi2=4mApi4 \sqrt{(p_i \cdot p_A)^2 - m_i^2 m_A^2} = 4 \sqrt{E_i^2 m_A^2 - m_i^2 m_A^2} = 4 m_A \sqrt{E_i^2 - m_i^2} = 4 m_A |\vec{p}_i|

2. 相空间积分化简 将四维动量守恒的 δ\delta 函数分离为能量部分和三维动量部分: δ(4)(pi+pApfpB)=δ(Ei+mAEfEB)δ(3)(pipfpB)\delta^{(4)}(p_i + p_A - p_f - p_B) = \delta(E_i + m_A - E_f - E_B) \delta^{(3)}(\vec{p}_i - \vec{p}_f - \vec{p}_B) 利用 δ(3)(pipfpB)\delta^{(3)}(\vec{p}_i - \vec{p}_f - \vec{p}_B)d3pBd^3p_B 进行积分,这要求末态粒子 BB 的动量满足 pB=pipf\vec{p}_B = \vec{p}_i - \vec{p}_f。此时粒子 BB 的能量变为: EB=pB2+mB2=(pipf)2+mB2E_B = \sqrt{\vec{p}_B^2 + m_B^2} = \sqrt{(\vec{p}_i - \vec{p}_f)^2 + m_B^2} 积分后的相空间变为:

dΠLIPS=1(2π)2d3pf4EfEBδ(Ei+mAEfEB)\int d\Pi_{\text{LIPS}} = \frac{1}{(2\pi)^2} \int \frac{d^3p_f}{4 E_f E_B} \delta(E_i + m_A - E_f - E_B)

在球坐标系下展开 d3pf=pf2dpfdΩd^3p_f = |\vec{p}_f|^2 d|\vec{p}_f| d\Omega,其中 dΩd\Omega 是出射粒子 ff 的立体角微元。此时微分截面可以写为:

dσdΩ=14mApi116π2pf2dpfEfEBM2δ(Ei+mAEfEB)\frac{d\sigma}{d\Omega} = \frac{1}{4 m_A |\vec{p}_i|} \frac{1}{16\pi^2} \int \frac{|\vec{p}_f|^2 d|\vec{p}_f|}{E_f E_B} |\mathcal{M}|^2 \delta(E_i + m_A - E_f - E_B)

3. 能量守恒 δ\delta 函数的积分与雅可比行列式 为了完成对 pf|\vec{p}_f| 的积分,需要利用一维 δ\delta 函数的性质: f(x)δ(g(x))dx=[f(x)g(x)]g(x)=0\int f(x) \delta(g(x)) dx = \left[ \frac{f(x)}{|g'(x)|} \right]_{g(x)=0} 令总能量函数为 Etot(pf)=Ef+EBE_{\text{tot}}(|\vec{p}_f|) = E_f + E_B,我们需要计算其对 pf|\vec{p}_f| 的导数。 已知 Ef2=pf2+mf2E_f^2 = |\vec{p}_f|^2 + m_f^2,对其两边求导得: 2EfEfpf=2pf    Efpf=pfEf2 E_f \frac{\partial E_f}{\partial |\vec{p}_f|} = 2 |\vec{p}_f| \implies \frac{\partial E_f}{\partial |\vec{p}_f|} = \frac{|\vec{p}_f|}{E_f} 对于 EBE_B,展开其动量平方项: EB2=pi2+pf22pipfcosθ+mB2E_B^2 = |\vec{p}_i|^2 + |\vec{p}_f|^2 - 2 |\vec{p}_i| |\vec{p}_f| \cos\theta + m_B^2 其中 θ\thetapi\vec{p}_ipf\vec{p}_f 之间的夹角。在给定散射角 θ\theta 的情况下对 pf|\vec{p}_f| 求导: 2EBEBpf=2pf2picosθ    EBpf=pfpicosθEB2 E_B \frac{\partial E_B}{\partial |\vec{p}_f|} = 2 |\vec{p}_f| - 2 |\vec{p}_i| \cos\theta \implies \frac{\partial E_B}{\partial |\vec{p}_f|} = \frac{|\vec{p}_f| - |\vec{p}_i| \cos\theta}{E_B} 因此,总能量对 pf|\vec{p}_f| 的导数为:

(Ef+EB)pf=pfEf+pfpicosθEB=pfEB+Ef(pfpicosθ)EfEB\frac{\partial (E_f + E_B)}{\partial |\vec{p}_f|} = \frac{|\vec{p}_f|}{E_f} + \frac{|\vec{p}_f| - |\vec{p}_i| \cos\theta}{E_B} = \frac{|\vec{p}_f| E_B + E_f (|\vec{p}_f| - |\vec{p}_i| \cos\theta)}{E_f E_B}

提取公因子 pfEfEB\frac{|\vec{p}_f|}{E_f E_B},得到:

(Ef+EB)pf=pfEfEB[EB+Ef(1pipfcosθ)]\frac{\partial (E_f + E_B)}{\partial |\vec{p}_f|} = \frac{|\vec{p}_f|}{E_f E_B} \left[ E_B + E_f \left( 1 - \frac{|\vec{p}_i|}{|\vec{p}_f|} \cos\theta \right) \right]

4. 综合结果 将上述导数代入 δ\delta 函数的积分公式中,完成对 pf|\vec{p}_f| 的积分:

pf2EfEBδ(Ei+mAEfEB)dpf=pf2EfEB×EfEBpf[EB+Ef(1pipfcosθ)]\int \frac{|\vec{p}_f|^2}{E_f E_B} \delta(E_i + m_A - E_f - E_B) d|\vec{p}_f| = \frac{|\vec{p}_f|^2}{E_f E_B} \times \frac{E_f E_B}{|\vec{p}_f| \left[ E_B + E_f \left( 1 - \frac{|\vec{p}_i|}{|\vec{p}_f|} \cos\theta \right) \right]}

消去 EfEBE_f E_B 并约分 pf|\vec{p}_f|,得到:

pfEB+Ef(1pipfcosθ)\frac{|\vec{p}_f|}{E_B + E_f \left( 1 - \frac{|\vec{p}_i|}{|\vec{p}_f|} \cos\theta \right)}

最后,将此结果代回微分截面的表达式中,并合并常数因子 14mApi×116π2=164π2mApi\frac{1}{4 m_A |\vec{p}_i|} \times \frac{1}{16\pi^2} = \frac{1}{64\pi^2 m_A |\vec{p}_i|}

dσdΩ=164π2mApipfEB+Ef(1pipfcosθ)M2\frac{d\sigma}{d\Omega} = \frac{1}{64\pi^2 m_A |\vec{p}_i|} \frac{|\vec{p}_f|}{E_B + E_f \left( 1 - \frac{|\vec{p}_i|}{|\vec{p}_f|} \cos\theta \right)} |\mathcal{M}|^2

pi|\vec{p}_i| 移至分子中 pf|\vec{p}_f| 的分母位置,即可得到最终的微分截面公式:

dσdΩ=164π2mA[EB+Ef(1pipfcosθ)]1pfpiM2\boxed{ \frac{d\sigma}{d\Omega} = \frac{1}{64\pi^2 m_A} \left[ E_B + E_f \left( 1 - \frac{|\vec{p}_i|}{|\vec{p}_f|} \cos\theta \right) \right]^{-1} \frac{|\vec{p}_f|}{|\vec{p}_i|} |\mathcal{M}|^2 }
5.2

Problem 5.2

schwarzChapter 5

习题 5.2

来源: 第5章, PDF第67页


5.2 Show that dΠLIPSd\Pi_{\text{LIPS}} is Lorentz invariant and verify Eq. (5.21).

Referenced Equations:

Equation (5.21):

dΠLIPS(2π)4δ4(Σp)×final states jd3pj(2π)312Epj(5.21)d\Pi_{\text{LIPS}} \equiv (2\pi)^4\delta^4(\Sigma p) \times \prod_{\text{final states } j} \frac{d^3p_j}{(2\pi)^3} \frac{1}{2E_{p_j}} \tag{5.21}

习题 5.2 - 解答


洛伦兹不变相空间(Lorentz Invariant Phase Space, LIPS)测度是量子场论中计算散射截面和衰变宽度的核心部分。要证明 dΠLIPSd\Pi_{\text{LIPS}} 的洛伦兹不变性并验证公式 (5.21),需要从明显具有洛伦兹不变性的四维动量空间积分出发,并引入物理上的在壳(on-shell)条件。

对于任意单个末态粒子 jj,其四维动量为 pjμ=(pj0,pj)p_j^\mu = (p_j^0, \mathbf{p}_j)。在没有任何约束的情况下,动量空间的积分测度为 d4pj(2π)4\frac{d^4p_j}{(2\pi)^4}。由于 d4pjd^4p_j 是四维体积元,其在洛伦兹变换下是不变的(洛伦兹变换矩阵的雅可比行列式的绝对值为 1)。

物理上的末态粒子必须满足质量壳条件 pj2=mj2p_j^2 = m_j^2 且能量为正 pj0>0p_j^0 > 0。我们可以通过引入狄拉克 δ\delta 函数和阶跃函数 θ(pj0)\theta(p_j^0) 来施加这些物理约束。构造如下的单粒子积分测度: dμ(pj)=d4pj(2π)4(2π)δ(pj2mj2)θ(pj0)d\mu(p_j) = \frac{d^4p_j}{(2\pi)^4} (2\pi) \delta(p_j^2 - m_j^2) \theta(p_j^0)

分析该表达式的洛伦兹不变性:

  1. d4pj(2π)4\frac{d^4p_j}{(2\pi)^4} 是洛伦兹不变的。
  2. pj2mj2p_j^2 - m_j^2 是洛伦兹标量,因此 δ(pj2mj2)\delta(p_j^2 - m_j^2) 是洛伦兹不变的。
  3. 对于类时或类光矢量(由 δ(pj2mj2)\delta(p_j^2 - m_j^2) 保证),其时间分量 pj0p_j^0 的符号在正规正时洛伦兹变换(proper orthochronous Lorentz transformations)下保持不变,因此 θ(pj0)\theta(p_j^0) 也是洛伦兹不变的。

综上所述,单粒子测度 dμ(pj)d\mu(p_j) 是明显洛伦兹不变的。

接下来,我们对 pj0p_j^0 进行积分以化简该测度。利用 δ\delta 函数的复合函数性质: δ(f(x))=iδ(xxi)f(xi)\delta(f(x)) = \sum_i \frac{\delta(x - x_i)}{|f'(x_i)|} 其中 xix_if(x)=0f(x)=0 的单根。 在这里,变量为 pj0p_j^0,函数为 f(pj0)=pj2mj2=(pj0)2pj2mj2f(p_j^0) = p_j^2 - m_j^2 = (p_j^0)^2 - |\mathbf{p}_j|^2 - m_j^2。 令 Epj=pj2+mj2E_{p_j} = \sqrt{|\mathbf{p}_j|^2 + m_j^2},则方程 f(pj0)=0f(p_j^0) = 0 的根为 pj0=±Epjp_j^0 = \pm E_{p_j}。 对 pj0p_j^0 求导得到 fpj0=2pj0\left| \frac{\partial f}{\partial p_j^0} \right| = 2|p_j^0|

因此,δ\delta 函数可以展开为: δ(pj2mj2)=δ(pj0Epj)2Epj+δ(pj0+Epj)2Epj\delta(p_j^2 - m_j^2) = \frac{\delta(p_j^0 - E_{p_j})}{2E_{p_j}} + \frac{\delta(p_j^0 + E_{p_j})}{2E_{p_j}}

代入单粒子测度表达式中,并利用阶跃函数 θ(pj0)\theta(p_j^0) 滤除负能量解(即第二项): dμ(pj)=d3pjdpj0(2π)3[δ(pj0Epj)2Epj+δ(pj0+Epj)2Epj]θ(pj0)d\mu(p_j) = \frac{d^3p_j dp_j^0}{(2\pi)^3} \left[ \frac{\delta(p_j^0 - E_{p_j})}{2E_{p_j}} + \frac{\delta(p_j^0 + E_{p_j})}{2E_{p_j}} \right] \theta(p_j^0) dμ(pj)=d3pjdpj0(2π)3δ(pj0Epj)2Epjd\mu(p_j) = \frac{d^3p_j dp_j^0}{(2\pi)^3} \frac{\delta(p_j^0 - E_{p_j})}{2E_{p_j}}

pj0p_j^0 积分后,我们得到三维动量空间中的单粒子相空间测度: dμ(pj)=d3pj(2π)312Epj\int d\mu(p_j) = \frac{d^3p_j}{(2\pi)^3} \frac{1}{2E_{p_j}} 由于 dμ(pj)d\mu(p_j) 是洛伦兹不变的,推导出的三维形式 d3pj(2π)32Epj\frac{d^3p_j}{(2\pi)^3 2E_{p_j}} 必然也是洛伦兹不变的。

对于包含多个末态粒子的散射或衰变过程,总的相空间测度是所有单粒子相空间测度的乘积,并附加一个全局的能量-动量守恒约束。假设初态总四维动量为 PinitialP_{\text{initial}},末态总四维动量为 jpj\sum_j p_j,能量-动量守恒由四维 δ\delta 函数表示: (2π)4δ4(Pinitialjpj)(2\pi)^4 \delta^4\left(P_{\text{initial}} - \sum_j p_j\right) 由于四维 δ\delta 函数的自变量是四维动量之差(洛伦兹矢量),该 δ\delta 函数在洛伦兹变换下也是不变的。

将全局守恒律与所有末态粒子的洛伦兹不变测度相乘,即可得到总的洛伦兹不变相空间(LIPS)测度: dΠLIPS=(2π)4δ4(Pinitialjpj)final states j(d3pj(2π)312Epj)d\Pi_{\text{LIPS}} = (2\pi)^4 \delta^4\left(P_{\text{initial}} - \sum_j p_j\right) \prod_{\text{final states } j} \left( \frac{d^3p_j}{(2\pi)^3} \frac{1}{2E_{p_j}} \right)

题目公式 (5.21) 中的 δ4(Σp)\delta^4(\Sigma p)δ4(Pinitialjpj)\delta^4(P_{\text{initial}} - \sum_j p_j) 的简写形式,表示初末态总动量守恒。因为该表达式完全由明显洛伦兹不变的四维 δ\delta 函数和洛伦兹不变的单粒子测度乘积构成,所以 dΠLIPSd\Pi_{\text{LIPS}} 整体是洛伦兹不变的。由此验证了公式 (5.21):

dΠLIPS(2π)4δ4(Σp)×final states jd3pj(2π)312Epj\boxed{ d\Pi_{\text{LIPS}} \equiv (2\pi)^4\delta^4(\Sigma p) \times \prod_{\text{final states } j} \frac{d^3p_j}{(2\pi)^3} \frac{1}{2E_{p_j}} }
5.3

Problem 5.3

schwarzChapter 5

习题 5.3

来源: 第5章, PDF第67页


5.3 A muon decays to an electron, an electron anti-neutrino and a muon neutrino, μeνμνˉe\mu^- \rightarrow e^- \nu_\mu \bar{\nu}_e. The matrix element for this process, ignoring the electron and neutrino masses, is given by M2=32GF2(m22mE)mE|\mathcal{M}|^2 = 32G_F^2 (m^2 - 2mE)mE, where mm is the mass of the muon and EE is the energy of the outgoing νe\nu_e. GF=1.166×105 GeV2G_F = 1.166 \times 10^{-5} \text{ GeV}^{-2} is the Fermi constant.

(a) Perform the integral over dΠLIPSd\Pi_{\text{LIPS}} to show that the decay rate is

Γ=GF2m5192π3.(5.55)\Gamma = \frac{G_F^2 m^5}{192\pi^3}. \tag{5.55}

(b) Compare your result to the observed values m=106 MeVm = 106 \text{ MeV} and τ=Γ1=2.20 μs\tau = \Gamma^{-1} = 2.20 \text{ } \mu\text{s}. How big is the discrepancy as a percentage? What might account for the discrepancy?

习题 5.3 - 解答


习题 5.3 分析与解答

(a) 衰变率 Γ\Gamma 的推导

在静止系中,μ\mu 子的四维动量为 P=(m,0)P = (m, \vec{0})。设出射的 ee^-νμ\nu_\muνˉe\bar{\nu}_e 的四维动量分别为 p1p_1p2p_2p3p_3,对应的能量为 E1E_1E2E_2E3E_3。题目中给定的 νˉe\bar{\nu}_e 能量为 E=E3E = E_3

根据衰变率的黄金定则,三体衰变的微分衰变率为: dΓ=12mM2dΠLIPSd\Gamma = \frac{1}{2m} |\mathcal{M}|^2 d\Pi_{\text{LIPS}} 其中,洛伦兹不变相空间 (LIPS) 测度为: dΠLIPS=(2π)4δ(4)(Pp1p2p3)d3p1(2π)32E1d3p2(2π)32E2d3p3(2π)32E3d\Pi_{\text{LIPS}} = (2\pi)^4 \delta^{(4)}(P - p_1 - p_2 - p_3) \frac{d^3p_1}{(2\pi)^3 2E_1} \frac{d^3p_2}{(2\pi)^3 2E_2} \frac{d^3p_3}{(2\pi)^3 2E_3}

首先对 p2\vec{p}_2 积分,利用三维动量 δ\delta 函数消去 d3p2d^3p_2,此时 p2=(p1+p3)\vec{p}_2 = -(\vec{p}_1 + \vec{p}_3),且 E2=p2E_2 = |\vec{p}_2|(忽略中微子质量): dΠLIPS=18(2π)5d3p1d3p3E1E2E3δ(mE1E2E3)d\Pi_{\text{LIPS}} = \frac{1}{8(2\pi)^5} \frac{d^3p_1 d^3p_3}{E_1 E_2 E_3} \delta(m - E_1 - E_2 - E_3)

d3p1d^3p_1d3p3d^3p_3 展开为球坐标。由于系统具有旋转对称性,可以先对 p1p_1 的方向积分得到 4π4\pi,再设 p1\vec{p}_1p3\vec{p}_3 的夹角为 θ\thetad3p1d3p3=(4πE12dE1)(2πE32dE3dcosθ)=8π2E12E32dE1dE3dcosθd^3p_1 d^3p_3 = (4\pi E_1^2 dE_1) (2\pi E_3^2 dE_3 d\cos\theta) = 8\pi^2 E_1^2 E_3^2 dE_1 dE_3 d\cos\theta 由余弦定理,E22=p1+p32=E12+E32+2E1E3cosθE_2^2 = |\vec{p}_1 + \vec{p}_3|^2 = E_1^2 + E_3^2 + 2E_1 E_3 \cos\theta。在给定 E1E_1E3E_3 时,对 cosθ\cos\theta 的微分可以转化为对 E2E_2 的微分: 2E2dE2=2E1E3dcosθ    dcosθ=E2E1E3dE22E_2 dE_2 = 2E_1 E_3 d\cos\theta \implies d\cos\theta = \frac{E_2}{E_1 E_3} dE_2 代入相空间表达式: dΠLIPS=8π2E12E328(2π)5E1E2E3E2E1E3dE1dE2dE3δ(mE1E2E3)=132π3dE1dE2dE3δ(mE1E2E3)d\Pi_{\text{LIPS}} = \frac{8\pi^2 E_1^2 E_3^2}{8(2\pi)^5 E_1 E_2 E_3} \frac{E_2}{E_1 E_3} dE_1 dE_2 dE_3 \delta(m - E_1 - E_2 - E_3) = \frac{1}{32\pi^3} dE_1 dE_2 dE_3 \delta(m - E_1 - E_2 - E_3)

利用能量 δ\delta 函数对 E2E_2 积分,得到 E2=mE1E3E_2 = m - E_1 - E_3。此时相空间简化为: dΠLIPS=132π3dE1dE3d\Pi_{\text{LIPS}} = \frac{1}{32\pi^3} dE_1 dE_3 接下来确定积分上下限。由物理条件 1cosθ1-1 \le \cos\theta \le 1 可得: (E1E3)2E22(E1+E3)2(E_1 - E_3)^2 \le E_2^2 \le (E_1 + E_3)^2 代入 E2=mE1E3E_2 = m - E_1 - E_3,解得运动学边界条件: E1m2,E3m2,E1+E3m2E_1 \le \frac{m}{2}, \quad E_3 \le \frac{m}{2}, \quad E_1 + E_3 \ge \frac{m}{2}

题目给定的矩阵元平方 M2=32GF2(m22mE)mE|\mathcal{M}|^2 = 32G_F^2 (m^2 - 2mE)mE 仅依赖于 E3=EE_3 = E。因此,我们可以先对 E1E_1 积分。对于固定的 EEE1E_1 的积分范围是 [m2E,m2][\frac{m}{2} - E, \frac{m}{2}]m2Em2dE1=m2(m2E)=E\int_{\frac{m}{2} - E}^{\frac{m}{2}} dE_1 = \frac{m}{2} - \left(\frac{m}{2} - E\right) = E 将积分后的相空间代入衰变率公式: dΓ=12mM2E32π3dE=12m[32GF2(m22mE)mE]E32π3dE=GF22π3(m22mE)E2dEd\Gamma = \frac{1}{2m} |\mathcal{M}|^2 \frac{E}{32\pi^3} dE = \frac{1}{2m} \left[ 32G_F^2 (m^2 - 2mE)mE \right] \frac{E}{32\pi^3} dE = \frac{G_F^2}{2\pi^3} (m^2 - 2mE) E^2 dE

最后,对 EE00m/2m/2 积分得到总衰变率: Γ=GF22π30m/2(m2E22mE3)dE=GF22π3[m2E332mE44]0m/2\Gamma = \frac{G_F^2}{2\pi^3} \int_0^{m/2} (m^2 E^2 - 2m E^3) dE = \frac{G_F^2}{2\pi^3} \left[ m^2 \frac{E^3}{3} - 2m \frac{E^4}{4} \right]_0^{m/2} Γ=GF22π3(m524m532)=GF22π3(m596)\Gamma = \frac{G_F^2}{2\pi^3} \left( \frac{m^5}{24} - \frac{m^5}{32} \right) = \frac{G_F^2}{2\pi^3} \left( \frac{m^5}{96} \right) Γ=GF2m5192π3\boxed{ \Gamma = \frac{G_F^2 m^5}{192\pi^3} }


(b) 理论值与观测值的对比及误差分析

已知参数:

  • m=106 MeV=0.106 GeVm = 106 \text{ MeV} = 0.106 \text{ GeV}
  • GF=1.166×105 GeV2G_F = 1.166 \times 10^{-5} \text{ GeV}^{-2}
  • 自然单位制转换常数 6.582×1025 GeVs\hbar \approx 6.582 \times 10^{-25} \text{ GeV}\cdot\text{s}

将数值代入 (a) 中得到的公式计算理论衰变率 Γcalc\Gamma_{\text{calc}}Γcalc=(1.166×105)2×(0.106)5192π31.3596×1010×1.3382×1055953.23.056×1019 GeV\Gamma_{\text{calc}} = \frac{(1.166 \times 10^{-5})^2 \times (0.106)^5}{192\pi^3} \approx \frac{1.3596 \times 10^{-10} \times 1.3382 \times 10^{-5}}{5953.2} \approx 3.056 \times 10^{-19} \text{ GeV} 对应的理论寿命 τcalc\tau_{\text{calc}} 为: τcalc=Γcalc=6.582×1025 GeVs3.056×1019 GeV2.154×106 s=2.154 μs\tau_{\text{calc}} = \frac{\hbar}{\Gamma_{\text{calc}}} = \frac{6.582 \times 10^{-25} \text{ GeV}\cdot\text{s}}{3.056 \times 10^{-19} \text{ GeV}} \approx 2.154 \times 10^{-6} \text{ s} = 2.154 \text{ } \mu\text{s}

与观测值 τobs=2.20 μs\tau_{\text{obs}} = 2.20 \text{ } \mu\text{s} 相比,百分比误差(Discrepancy)为: Discrepancy=τobsτcalcτobs×100%=2.202.1542.20×100%2.1%\text{Discrepancy} = \frac{| \tau_{\text{obs}} - \tau_{\text{calc}} |}{\tau_{\text{obs}}} \times 100\% = \frac{2.20 - 2.154}{2.20} \times 100\% \approx 2.1\% Discrepancy2.1%\boxed{ \text{Discrepancy} \approx 2.1\% }

导致该误差的物理与数值原因分析(What might account for the discrepancy):

  1. μ\mu 子质量的近似(主要原因): 衰变率对 μ\mu 子质量极其敏感(Γm5\Gamma \propto m^5)。真实的 μ\mu 子质量约为 105.66 MeV105.66 \text{ MeV},而题目中近似取为 106 MeV106 \text{ MeV}。质量上约 0.32%0.32\% 的高估,会导致衰变率被高估约 5×0.32%1.6%5 \times 0.32\% \approx 1.6\%,从而使计算出的寿命偏短。这解释了绝大部分的误差。
  2. 忽略了电子质量 mem_e 题目假设电子无质量。实际上 me0.511 MeVm_e \approx 0.511 \text{ MeV},非零的电子质量会略微减小可用的相空间体积。引入电子质量会产生一个相空间修正因子 f(me2/mμ2)18(me/mμ)20.9998f(m_e^2/m_\mu^2) \approx 1 - 8(m_e/m_\mu)^2 \approx 0.9998,这会使理论衰变率降低约 0.02%0.02\%,使寿命略微变长。
  3. 忽略了辐射修正(Radiative Corrections): 树图级别的计算没有包含量子电动力学(QED)的高阶修正(如虚光子圈图修正和真实光子的轫致辐射)。QED 修正会给衰变率乘上一个因子 1α2π(π2254)0.996\approx 1 - \frac{\alpha}{2\pi}(\pi^2 - \frac{25}{4}) \approx 0.996,这会使衰变率降低约 0.4%0.4\%,进一步弥补理论寿命与实验寿命之间的差距。
5.4

Problem 5.4

schwarzChapter 5

习题 5.4

来源: 第5章, PDF第67页


5.4 Repeat the e+eμ+μe^+ e^- \rightarrow \mu^+ \mu^- calculation in Section 5.3 using circular polarizations.

习题 5.4 - 解答


在质心系中考察过程 e(p)e+(p)μ(k)μ+(k)e^-(p) e^+(p') \rightarrow \mu^-(k) \mu^+(k')。由于电子质量 meEm_e \ll EEE 为质心系单粒子能量),我们可以采用高能近似 me0m_e \approx 0,而保留 μ\mu 子的质量 mμm_\mu。在此近似下,电子的螺旋度守恒,只有初始状态为 eReL+e^-_R e^+_LeLeR+e^-_L e^+_R 时,跃迁矩阵元才不为零。

1. 运动学与电子流 设电子束沿 zz 轴入射,μ\mu 子出射角为 θ\theta。四动量分别为: p=(E,0,0,E),p=(E,0,0,E)p = (E, 0, 0, E), \quad p' = (E, 0, 0, -E) k=(E,k),k=(E,k)k = (E, \vec{k}), \quad k' = (E, -\vec{k}) 其中 k=k(sinθ,0,cosθ)\vec{k} = k(\sin\theta, 0, \cos\theta),且 k=k=E2mμ2k = |\vec{k}| = \sqrt{E^2 - m_\mu^2}

在手征表示下,无质量电子的非零流矩阵元 jμ=vˉ(p)γμu(p)j^\mu = \bar{v}(p') \gamma^\mu u(p) 为: jRLμ=vˉL(p)γμuR(p)=2E(0,1,i,0)j_{RL}^\mu = \bar{v}_L(p') \gamma^\mu u_R(p) = 2E(0, -1, -i, 0) jLRμ=vˉR(p)γμuL(p)=2E(0,1,i,0)j_{LR}^\mu = \bar{v}_R(p') \gamma^\mu u_L(p) = 2E(0, -1, i, 0)

2. μ\mu 子的螺旋度旋量 对于具有质量 mμm_\muμ\mu 子,沿 k\vec{k} 方向的螺旋度基旋量为: ξ+=(cosθ2sinθ2),ξ=(sinθ2cosθ2)\xi_+ = \begin{pmatrix} \cos\frac{\theta}{2} \\ \sin\frac{\theta}{2} \end{pmatrix}, \quad \xi_- = \begin{pmatrix} -\sin\frac{\theta}{2} \\ \cos\frac{\theta}{2} \end{pmatrix} 利用手征表示下的旋量公式 u(p)=(pσξpσˉξ)u(p) = \begin{pmatrix} \sqrt{p \cdot \sigma} \xi \\ \sqrt{p \cdot \bar{\sigma}} \xi \end{pmatrix}v(p)=(pσηpσˉη)v(p) = \begin{pmatrix} \sqrt{p \cdot \sigma} \eta \\ -\sqrt{p \cdot \bar{\sigma}} \eta \end{pmatrix},可得 μ\mu^- 的旋量: u(k,+)=(Ekξ+E+kξ+),u(k,)=(E+kξEkξ)u(k, +) = \begin{pmatrix} \sqrt{E-k} \xi_+ \\ \sqrt{E+k} \xi_+ \end{pmatrix}, \quad u(k, -) = \begin{pmatrix} \sqrt{E+k} \xi_- \\ \sqrt{E-k} \xi_- \end{pmatrix} 对于 μ+\mu^+,动量方向为 k-\vec{k}。右旋正电子(螺旋度 +1+1)对应的物理自旋与动量同向,其 vv 旋量需取相反的自旋基 η=ξ+\eta = \xi_+;左旋正电子(螺旋度 1-1)对应 η=ξ\eta = \xi_-v(k,+)=(E+kξ+Ekξ+),v(k,)=(EkξE+kξ)v(k', +) = \begin{pmatrix} \sqrt{E+k} \xi_+ \\ -\sqrt{E-k} \xi_+ \end{pmatrix}, \quad v(k', -) = \begin{pmatrix} \sqrt{E-k} \xi_- \\ -\sqrt{E+k} \xi_- \end{pmatrix}

3. μ\mu 子流与矩阵元计算 μ\mu 子流定义为 Jμ(λ,λ+)=uˉ(k,λ)γμv(k,λ+)J^\mu(\lambda_-, \lambda_+) = \bar{u}(k, \lambda_-) \gamma^\mu v(k', \lambda_+)。 总散射振幅为 M=e2sjμJμ\mathcal{M} = \frac{e^2}{s} j_\mu J^\mu。由于 jRLj_{RL}jLRj_{LR} 只有横向空间分量,我们只需要计算 J1±iJ2J^1 \pm i J^2。 在手征表示下,uˉγiv=u2σiv2u1σiv1\bar{u} \gamma^i v = u_2^\dagger \sigma^i v_2 - u_1^\dagger \sigma^i v_1。引入升降算符 σ±=σ1±iσ2\sigma^\pm = \sigma^1 \pm i \sigma^2,可得: J1±iJ2=u1σ±v1+u2σ±v2J^1 \pm i J^2 = -u_1^\dagger \sigma^\pm v_1 + u_2^\dagger \sigma^\pm v_2 计算所需的泡利矩阵元为: ξ+σ+ξ+=sinθ,ξ+σ+ξ=1+cosθ,ξσ+ξ+=cosθ1,ξσ+ξ=sinθ\xi_+^\dagger \sigma^+ \xi_+ = \sin\theta, \quad \xi_+^\dagger \sigma^+ \xi_- = 1+\cos\theta, \quad \xi_-^\dagger \sigma^+ \xi_+ = \cos\theta-1, \quad \xi_-^\dagger \sigma^+ \xi_- = -\sin\theta ξ+σξ+=sinθ,ξ+σξ=cosθ1,ξσξ+=1+cosθ,ξσξ=sinθ\xi_+^\dagger \sigma^- \xi_+ = \sin\theta, \quad \xi_+^\dagger \sigma^- \xi_- = \cos\theta-1, \quad \xi_-^\dagger \sigma^- \xi_+ = 1+\cos\theta, \quad \xi_-^\dagger \sigma^- \xi_- = -\sin\theta

对于 eReL+e^-_R e^+_L 初态,MRL=e2s(2E)(J1+iJ2)=e22E(J1+iJ2)\mathcal{M}_{RL} = \frac{e^2}{s} (2E)(J^1 + i J^2) = \frac{e^2}{2E} (J^1 + i J^2)。代入旋量分量计算四种螺旋度组合:

  1. (+,+)(+, +): J1+iJ2=(Ekξ+)σ+(E+kξ+)+(E+kξ+)σ+(Ekξ+)=2E2k2sinθ=2mμsinθJ^1 + i J^2 = -(\sqrt{E-k}\xi_+^\dagger)\sigma^+(\sqrt{E+k}\xi_+) + (\sqrt{E+k}\xi_+^\dagger)\sigma^+(-\sqrt{E-k}\xi_+) = -2\sqrt{E^2-k^2}\sin\theta = -2m_\mu\sin\theta M(RL++)=e2mμEsinθ\mathcal{M}(RL \to ++) = -e^2 \frac{m_\mu}{E} \sin\theta
  2. (+,)(+, -): J1+iJ2=(Ek)(1+cosθ)(E+k)(1+cosθ)=2E(1+cosθ)J^1 + i J^2 = -(E-k)(1+\cos\theta) - (E+k)(1+\cos\theta) = -2E(1+\cos\theta) M(RL+)=e2(1+cosθ)\mathcal{M}(RL \to +-) = -e^2 (1+\cos\theta)
  3. (,+)(-, +): J1+iJ2=(E+k)(cosθ1)(Ek)(cosθ1)=2E(1cosθ)J^1 + i J^2 = -(E+k)(\cos\theta-1) - (E-k)(\cos\theta-1) = 2E(1-\cos\theta) M(RL+)=e2(1cosθ)\mathcal{M}(RL \to -+) = e^2 (1-\cos\theta)
  4. (,)(-, -): J1+iJ2=(E+kξ)σ+(Ekξ)+(Ekξ)σ+(E+kξ)=2mμsinθJ^1 + i J^2 = -(\sqrt{E+k}\xi_-^\dagger)\sigma^+(\sqrt{E-k}\xi_-) + (\sqrt{E-k}\xi_-^\dagger)\sigma^+(-\sqrt{E+k}\xi_-) = 2m_\mu\sin\theta M(RL)=e2mμEsinθ\mathcal{M}(RL \to --) = e^2 \frac{m_\mu}{E} \sin\theta

同理,对于 eLeR+e^-_L e^+_R 初态,MLR=e22E(J1iJ2)\mathcal{M}_{LR} = \frac{e^2}{2E} (J^1 - i J^2),利用 σ\sigma^- 的矩阵元可得: M(LR++)=e2mμEsinθ,M(LR+)=e2(1cosθ)\mathcal{M}(LR \to ++) = -e^2 \frac{m_\mu}{E} \sin\theta, \quad \mathcal{M}(LR \to +-) = e^2 (1-\cos\theta) M(LR+)=e2(1+cosθ),M(LR)=e2mμEsinθ\mathcal{M}(LR \to -+) = -e^2 (1+\cos\theta), \quad \mathcal{M}(LR \to --) = e^2 \frac{m_\mu}{E} \sin\theta

4. 非极化微商截面 对所有初始自旋求平均,对所有末态自旋求和: 14spinsM2=14(λ,λ+MRL2+λ,λ+MLR2)\frac{1}{4} \sum_{\text{spins}} |\mathcal{M}|^2 = \frac{1}{4} \left( \sum_{\lambda_-,\lambda_+} |\mathcal{M}_{RL}|^2 + \sum_{\lambda_-,\lambda_+} |\mathcal{M}_{LR}|^2 \right) 由于 MRL2|\mathcal{M}_{RL}|^2MLR2|\mathcal{M}_{LR}|^2 给出相同的贡献,我们有: 14spinsM2=12e4[2mμ2E2sin2θ+(1+cosθ)2+(1cosθ)2]\frac{1}{4} \sum_{\text{spins}} |\mathcal{M}|^2 = \frac{1}{2} e^4 \left[ 2\frac{m_\mu^2}{E^2}\sin^2\theta + (1+\cos\theta)^2 + (1-\cos\theta)^2 \right] =e4[mμ2E2sin2θ+1+cos2θ]=e4[1+mμ2E2+(1mμ2E2)cos2θ]= e^4 \left[ \frac{m_\mu^2}{E^2}\sin^2\theta + 1 + \cos^2\theta \right] = e^4 \left[ 1 + \frac{m_\mu^2}{E^2} + \left(1 - \frac{m_\mu^2}{E^2}\right)\cos^2\theta \right]

代入质心系两体相空间微商截面公式 dσdΩ=164π2skpM2\frac{d\sigma}{d\Omega} = \frac{1}{64\pi^2 s} \frac{|\vec{k}|}{|\vec{p}|} \overline{|\mathcal{M}|^2},其中 kp=kE=1mμ2E2\frac{|\vec{k}|}{|\vec{p}|} = \frac{k}{E} = \sqrt{1 - \frac{m_\mu^2}{E^2}},且 α=e24π\alpha = \frac{e^2}{4\pi}dσdΩ=e464π2s1mμ2E2[1+mμ2E2+(1mμ2E2)cos2θ]\frac{d\sigma}{d\Omega} = \frac{e^4}{64\pi^2 s} \sqrt{1 - \frac{m_\mu^2}{E^2}} \left[ 1 + \frac{m_\mu^2}{E^2} + \left(1 - \frac{m_\mu^2}{E^2}\right)\cos^2\theta \right]

化简后得到最终结果: dσdΩ=α24s1mμ2E2[1+mμ2E2+(1mμ2E2)cos2θ]\boxed{ \frac{d\sigma}{d\Omega} = \frac{\alpha^2}{4s} \sqrt{1 - \frac{m_\mu^2}{E^2}} \left[ 1 + \frac{m_\mu^2}{E^2} + \left(1 - \frac{m_\mu^2}{E^2}\right) \cos^2\theta \right] } 此结果与 Section 5.3 中使用迹定理(Trace theorems)计算得到的结果完全一致。

5.5

Problem 5.5

schwarzChapter 5

习题 5.5

来源: 第5章, PDF第67,68页


5.5 One of the most important scattering experiments ever was Rutherford's gold foil experiment. Rutherford scattering is αNαN\alpha N \rightarrow \alpha N, where NN is some atomic nucleus and α\alpha is an α\alpha-particle (helium nucleus). It is an almost identical process to Coulomb scattering (ep+ep+e^- p^+ \rightarrow e^- p^+).

(a) Look up or calculate the classical Rutherford scattering cross section. What assumptions go into its derivation?

(b) We showed that the quantum mechanical cross section for Coulomb scattering in Eq. (5.41) follows either from the Born approximation or from quantum field theory. Start from the formula for Coulomb scattering and make the appropriate replacements for αN\alpha N scattering.

(c) Draw the Feynman diagram for Rutherford scattering. What is the momentum of the virtual photon, kμk_{\mu}, in terms of the scattering angle and the energy of the incoming α\alpha-particle? (d) Substitute in for k4k^4 and rewrite the cross section in terms of the kinetic energy of the α\alpha-particle. Show that Rutherford's classical formula is reproduced. (e) Why are the classical and quantum answers the same? Could you have known this ahead of time? (f) Would the cross section for eeeee^- e^- \rightarrow e^- e^- also be given by the Coulomb scattering cross section?

Referenced Equations:

Equation (5.41):

(dσdΩ)CM=164π2mp2M2e4me24π21k4,(5.41)\left( \frac{d\sigma}{d\Omega} \right)_{\text{CM}} = \frac{1}{64\pi^2 m_p^2} |\mathcal{M}|^2 \sim \frac{e^4 m_e^2}{4\pi^2} \frac{1}{\vec{k}^4}, \tag{5.41}

习题 5.5 - 解答


(a) 经典卢瑟福散射截面及其假设

在经典力学中,考虑一个电荷量为 q1=zeq_1 = z e 的入射粒子与一个电荷量为 q2=Zeq_2 = Z e 的靶核发生库仑散射。通过求解经典轨道方程(双曲线轨道),可以得到瞄准距 bb 与散射角 θ\theta 的关系: b=q1q28πEKcotθ2b = \frac{q_1 q_2}{8\pi E_K} \cot\frac{\theta}{2} 其中 EK=12mv2E_K = \frac{1}{2}m v^2 是入射粒子的动能(采用自然单位制或 Heaviside-Lorentz 单位制,库仑力常数为 1/4π1/4\pi)。 微分散射截面由 dσ=2πbdbd\sigma = 2\pi b |db| 和立体角微元 dΩ=2πsinθdθd\Omega = 2\pi \sin\theta d\theta 给出: dσdΩ=bsinθdbdθ=(q1q216πEKsin2(θ/2))2\frac{d\sigma}{d\Omega} = \frac{b}{\sin\theta} \left| \frac{db}{d\theta} \right| = \left( \frac{q_1 q_2}{16\pi E_K \sin^2(\theta/2)} \right)^2 对于 α\alpha 粒子,电荷数 z=2z=2,代入 q1=2eq_1 = 2eq2=Zeq_2 = Ze,得到经典卢瑟福散射截面: (dσdΩ)classical=Z2e464π2EK2sin4(θ/2)\boxed{ \left(\frac{d\sigma}{d\Omega}\right)_{\text{classical}} = \frac{Z^2 e^4}{64 \pi^2 E_K^2 \sin^4(\theta/2)} }

推导该公式所依赖的假设包括:

  1. 靶核无限重:忽略靶核的反冲,质心系与实验室系重合。
  2. 纯经典库仑相互作用:假设 α\alpha 粒子和原子核均为无内部结构的经典点电荷,仅受静电斥力作用。
  3. 非相对论极限α\alpha 粒子的运动速度远小于光速,动能可表示为 EK=p2/2mαE_K = p^2/2m_\alpha
  4. 无辐射损失:忽略粒子在加速偏转过程中产生的轫致辐射。
  5. 无自旋效应:将粒子视为无自旋的经典质点。

(b) 从库仑散射公式推广至 αN\alpha N 散射

题目给出的电子-质子库仑散射截面公式为: (dσdΩ)CMe4me24π2k4\left( \frac{d\sigma}{d\Omega} \right)_{\text{CM}} \sim \frac{e^4 m_e^2}{4\pi^2 \vec{k}^4} 对于 α\alpha 粒子与原子核 NN 的散射,需要对相互作用顶点和运动学质量进行替换:

  1. 电荷替换:电子电荷 ee 替换为 α\alpha 粒子的电荷 2e2e;质子电荷 ee 替换为靶核电荷 ZeZe。因此,散射振幅平方中的耦合常数因子由 e4e^4 变为 (2eZe)2=4Z2e4(2e \cdot Ze)^2 = 4Z^2 e^4
  2. 质量替换:入射电子质量 mem_e 替换为 α\alpha 粒子质量 mαm_\alpha

进行上述替换后,量子力学给出的 αN\alpha N 散射截面公式为: dσdΩ=4Z2e4mα24π2k4=Z2e4mα2π2k4\boxed{ \frac{d\sigma}{d\Omega} = \frac{4 Z^2 e^4 m_\alpha^2}{4\pi^2 \vec{k}^4} = \frac{Z^2 e^4 m_\alpha^2}{\pi^2 \vec{k}^4} }


(c) 费曼图与虚光子动量

费曼图描述: 卢瑟福散射在量子场论中对应于 tt-通道的树图(Tree-level tt-channel diagram)。入射 α\alpha 粒子(动量 p1p_1)与靶核 NN(动量 p2p_2)通过交换一个虚光子(动量 kk)发生相互作用,散射后变为出射 α\alpha 粒子(动量 p3p_3)和出射靶核 NN(动量 p4p_4)。

虚光子动量 kμk_\mu 的计算: 根据动量守恒,虚光子的四维动量为 kμ=(p1p3)μk_\mu = (p_1 - p_3)_\mu。 由于靶核被假设为无限重,它不吸收能量,因此 α\alpha 粒子发生的是弹性散射,能量守恒给出 E1=E3=EE_1 = E_3 = E。 虚光子的能量分量为 k0=E1E3=0k_0 = E_1 - E_3 = 0。 虚光子的三维动量大小平方为: k2=(p1p3)2=p12+p322p1p3cosθ\vec{k}^2 = (\vec{p}_1 - \vec{p}_3)^2 = |\vec{p}_1|^2 + |\vec{p}_3|^2 - 2|\vec{p}_1||\vec{p}_3|\cos\theta 由于弹性散射 p1=p3=p|\vec{p}_1| = |\vec{p}_3| = p,利用半角公式 1cosθ=2sin2(θ/2)1 - \cos\theta = 2\sin^2(\theta/2),得到: k2=2p2(1cosθ)=4p2sin2θ2\vec{k}^2 = 2p^2(1 - \cos\theta) = 4p^2 \sin^2\frac{\theta}{2}α\alpha 粒子的动能 EK=p2/2mαE_K = p^2 / 2m_\alpha 代换 p2=2mαEKp^2 = 2m_\alpha E_K,得到: k2=8mαEKsin2θ2\vec{k}^2 = 8 m_\alpha E_K \sin^2\frac{\theta}{2} 因此,虚光子的四维动量 kμk_\mu 及其平方 k2k^2 为: kμ=(0,p1p3),k2=k2=8mαEKsin2θ2\boxed{ k_\mu = (0, \vec{p}_1 - \vec{p}_3), \quad k^2 = -|\vec{k}|^2 = -8 m_\alpha E_K \sin^2\frac{\theta}{2} }


(d) 验证量子截面与经典公式的一致性

将 (c) 中得到的 k2\vec{k}^2 平方后代入 (b) 中的截面公式: k4=(8mαEKsin2θ2)2=64mα2EK2sin4θ2\vec{k}^4 = \left( 8 m_\alpha E_K \sin^2\frac{\theta}{2} \right)^2 = 64 m_\alpha^2 E_K^2 \sin^4\frac{\theta}{2} 代入量子力学截面公式: dσdΩ=Z2e4mα2π2(64mα2EK2sin4(θ/2))\frac{d\sigma}{d\Omega} = \frac{Z^2 e^4 m_\alpha^2}{\pi^2 \left( 64 m_\alpha^2 E_K^2 \sin^4(\theta/2) \right)} 消去分子分母中的 mα2m_\alpha^2,得到: dσdΩ=Z2e464π2EK2sin4(θ/2)\boxed{ \frac{d\sigma}{d\Omega} = \frac{Z^2 e^4}{64 \pi^2 E_K^2 \sin^4(\theta/2)} } 这与 (a) 中推导的经典卢瑟福散射公式完全一致。


(e) 经典与量子结果相同的原因

原因分析: 经典力学与量子力学(玻恩近似)在卢瑟福散射中给出完全相同的微分散射截面,这是 V(r)1/rV(r) \propto 1/r 库仑势的一个独有的数学奇迹。在量子力学中,库仑势的高阶微扰项(玻恩级数的高阶项)仅对散射振幅贡献一个无法被观测到的纯相位因子(Coulomb phase),因此一阶微扰(树图级)给出的截面就是精确的。这在深层次上与库仑/开普勒问题中隐藏的 SO(4)SO(4) 动力学对称性(龙格-楞次矢量的守恒)密切相关。

能否提前预知: 可以提前预知。从物理学史的角度看,卢瑟福最初完全基于经典力学推导出了该公式,并且该公式与盖革-马尔登的 α\alpha 粒子散射实验数据完美吻合。既然微观粒子的真实散射过程本质上是量子力学的,而经典公式却能完美拟合实验,这就意味着对于该特定势场,经典推导与量子推导必然给出相同的数学结果。


(f) 电子-电子散射 (eeeee^- e^- \rightarrow e^- e^-) 是否适用该公式?

不适用。电子-电子散射(Møller 散射)的截面不能用上述简单的库仑散射公式描述,原因如下:

  1. 全同粒子效应(不可分辨性):两个电子是全同费米子。在量子场论中,除了 tt-通道图外,还必须包含交换两个出射电子的 uu-通道图。这两个振幅会发生量子干涉,且截面必须在交换两个出射粒子时保持对称。
  2. 反冲效应不可忽略:靶电子与入射电子质量相等,靶粒子无限重的假设完全失效,必须严格处理质心系与实验室系的运动学转换及反冲动能。
  3. 自旋与相对论效应:电子是自旋为 1/21/2 的粒子,自旋-自旋相互作用(磁矩相互作用)在相对论能量下变得极其重要,这会显著改变散射截面的角分布。
5.6

Problem 5.6

schwarzChapter 5

习题 5.6

来源: 第5章, PDF第68页


5.6 In Section 5.3 we found that the e+eμ+μe^+ e^- \rightarrow \mu^+ \mu^- cross section had the form dσdΩ=e464π2ECM2(1+cos2θ)\frac{d\sigma}{d\Omega} = \frac{e^4}{64\pi^2 E_{\text{CM}}^2} (1 + \cos^2 \theta) in the center-of-mass frame. (a) Work out the Lorentz-invariant quantities s=(pe++pe)2s = (p_{e^+} + p_{e^-})^2, t=(pμpe)2t = (p_{\mu^-} - p_{e^-})^2 and u=(pμ+pe)2u = (p_{\mu^+} - p_{e^-})^2 in terms of ECME_{\text{CM}} and cosθ\cos \theta (still assuming mμ=me=0m_{\mu} = m_e = 0). (b) Derive a relationship between s,ts, t and uu. (c) Rewrite dσdΩ\frac{d\sigma}{d\Omega} in terms of s,ts, t and uu. (d) Now assume mμm_{\mu} and mem_e are non-zero. Derive a relationship between s,ts, t and uu and the masses.

习题 5.6 - 解答


本题主要考察 222 \to 2 散射过程中的运动学变量(Mandelstam 变量)的定义、计算及其在质心系(Center-of-Mass frame)下的具体形式。

(a) 在质心系中,设正负电子的束流方向为 zz 轴。在零质量极限(me=mμ=0m_e = m_\mu = 0)下,粒子的能量 EE 与动量大小相等。总质心能量为 ECME_{\text{CM}},则每个入射粒子的能量为 E=ECM/2E = E_{\text{CM}}/2。 入射正负电子的四维动量可表示为: pe=(E,0,0,E)p_{e^-} = (E, 0, 0, E) pe+=(E,0,0,E)p_{e^+} = (E, 0, 0, -E) 出射正负缪子的四维动量(散射角为 θ\theta)可表示为: pμ=(E,Esinθ,0,Ecosθ)p_{\mu^-} = (E, E\sin\theta, 0, E\cos\theta) pμ+=(E,Esinθ,0,Ecosθ)p_{\mu^+} = (E, -E\sin\theta, 0, -E\cos\theta)

根据 Mandelstam 变量的定义,利用闵可夫斯基度规 (+,,,)(+, -, -, -) 计算内积: 对于 sss=(pe++pe)2=(2E,0,0,0)2=4E2=ECM2s = (p_{e^+} + p_{e^-})^2 = (2E, 0, 0, 0)^2 = 4E^2 = E_{\text{CM}}^2 对于 ttt=(pμpe)2=pμ2+pe22pμpet = (p_{\mu^-} - p_{e^-})^2 = p_{\mu^-}^2 + p_{e^-}^2 - 2p_{\mu^-} \cdot p_{e^-} 由于零质量极限下 p2=0p^2 = 0,有: t=2pμpe=2(EEEcosθE)=2E2(1cosθ)=ECM22(1cosθ)t = -2p_{\mu^-} \cdot p_{e^-} = -2(E \cdot E - E\cos\theta \cdot E) = -2E^2(1 - \cos\theta) = -\frac{E_{\text{CM}}^2}{2}(1 - \cos\theta) 对于 uuu=(pμ+pe)2=pμ+2+pe22pμ+pe=2pμ+peu = (p_{\mu^+} - p_{e^-})^2 = p_{\mu^+}^2 + p_{e^-}^2 - 2p_{\mu^+} \cdot p_{e^-} = -2p_{\mu^+} \cdot p_{e^-} u=2[EE(Ecosθ)E]=2E2(1+cosθ)=ECM22(1+cosθ)u = -2[E \cdot E - (-E\cos\theta) \cdot E] = -2E^2(1 + \cos\theta) = -\frac{E_{\text{CM}}^2}{2}(1 + \cos\theta)

最终结果为: s=ECM2,t=ECM22(1cosθ),u=ECM22(1+cosθ)\boxed{s = E_{\text{CM}}^2, \quad t = -\frac{E_{\text{CM}}^2}{2}(1 - \cos\theta), \quad u = -\frac{E_{\text{CM}}^2}{2}(1 + \cos\theta)}

(b) 将 (a) 中求得的 s,t,us, t, u 表达式直接相加: s+t+u=ECM2ECM22(1cosθ)ECM22(1+cosθ)s + t + u = E_{\text{CM}}^2 - \frac{E_{\text{CM}}^2}{2}(1 - \cos\theta) - \frac{E_{\text{CM}}^2}{2}(1 + \cos\theta) s+t+u=ECM2ECM22[(1cosθ)+(1+cosθ)]=ECM2ECM22(2)=0s + t + u = E_{\text{CM}}^2 - \frac{E_{\text{CM}}^2}{2} \left[ (1 - \cos\theta) + (1 + \cos\theta) \right] = E_{\text{CM}}^2 - \frac{E_{\text{CM}}^2}{2}(2) = 0 这体现了在无质量极限下,Mandelstam 变量之和为零的运动学特性。 s+t+u=0\boxed{s + t + u = 0}

(c) 需要将微分散射截面 dσdΩ=e464π2ECM2(1+cos2θ)\frac{d\sigma}{d\Omega} = \frac{e^4}{64\pi^2 E_{\text{CM}}^2} (1 + \cos^2 \theta)s,t,us, t, u 表示。 由 (a) 可知 ECM2=sE_{\text{CM}}^2 = s。 为了构造出 1+cos2θ1 + \cos^2\theta,我们考察 t2+u2t^2 + u^2t2+u2=[s2(1cosθ)]2+[s2(1+cosθ)]2t^2 + u^2 = \left[-\frac{s}{2}(1 - \cos\theta)\right]^2 + \left[-\frac{s}{2}(1 + \cos\theta)\right]^2 t2+u2=s24[(1cosθ)2+(1+cosθ)2]=s24(2+2cos2θ)=s22(1+cos2θ)t^2 + u^2 = \frac{s^2}{4} \left[ (1 - \cos\theta)^2 + (1 + \cos\theta)^2 \right] = \frac{s^2}{4} (2 + 2\cos^2\theta) = \frac{s^2}{2}(1 + \cos^2\theta) 因此,角度依赖项可以写为: 1+cos2θ=2(t2+u2)s21 + \cos^2\theta = \frac{2(t^2 + u^2)}{s^2} 将其代入微分散射截面公式中: dσdΩ=e464π2s2(t2+u2)s2=e432π2s3(t2+u2)\frac{d\sigma}{d\Omega} = \frac{e^4}{64\pi^2 s} \cdot \frac{2(t^2 + u^2)}{s^2} = \frac{e^4}{32\pi^2 s^3} (t^2 + u^2) 最终结果为: dσdΩ=e4(t2+u2)32π2s3\boxed{\frac{d\sigma}{d\Omega} = \frac{e^4 (t^2 + u^2)}{32\pi^2 s^3}}

(d) 当质量不为零时,四维动量的平方满足在壳条件(On-shell condition): pe2=me2,pe+2=me2,pμ2=mμ2,pμ+2=mμ2p_{e^-}^2 = m_e^2, \quad p_{e^+}^2 = m_e^2, \quad p_{\mu^-}^2 = m_\mu^2, \quad p_{\mu^+}^2 = m_\mu^2 根据 Mandelstam 变量的定义展开: s=(pe++pe)2=pe+2+pe2+2pe+pe=2me2+2pe+pes = (p_{e^+} + p_{e^-})^2 = p_{e^+}^2 + p_{e^-}^2 + 2p_{e^+} \cdot p_{e^-} = 2m_e^2 + 2p_{e^+} \cdot p_{e^-} t=(pμpe)2=pμ2+pe22pμpe=mμ2+me22pμpet = (p_{\mu^-} - p_{e^-})^2 = p_{\mu^-}^2 + p_{e^-}^2 - 2p_{\mu^-} \cdot p_{e^-} = m_\mu^2 + m_e^2 - 2p_{\mu^-} \cdot p_{e^-} u=(pμ+pe)2=pμ+2+pe22pμ+pe=mμ2+me22pμ+peu = (p_{\mu^+} - p_{e^-})^2 = p_{\mu^+}^2 + p_{e^-}^2 - 2p_{\mu^+} \cdot p_{e^-} = m_\mu^2 + m_e^2 - 2p_{\mu^+} \cdot p_{e^-} 将三者相加: s+t+u=4me2+2mμ2+2pe(pe+pμpμ+)s + t + u = 4m_e^2 + 2m_\mu^2 + 2p_{e^-} \cdot (p_{e^+} - p_{\mu^-} - p_{\mu^+}) 根据四维动量守恒定律,有 pe+pe+=pμ+pμ+p_{e^-} + p_{e^+} = p_{\mu^-} + p_{\mu^+},移项可得: pe+pμpμ+=pep_{e^+} - p_{\mu^-} - p_{\mu^+} = -p_{e^-} 将其代入求和式中: s+t+u=4me2+2mμ2+2pe(pe)=4me2+2mμ22pe2s + t + u = 4m_e^2 + 2m_\mu^2 + 2p_{e^-} \cdot (-p_{e^-}) = 4m_e^2 + 2m_\mu^2 - 2p_{e^-}^2 再次利用在壳条件 pe2=me2p_{e^-}^2 = m_e^2,得到: s+t+u=4me2+2mμ22me2=2me2+2mμ2s + t + u = 4m_e^2 + 2m_\mu^2 - 2m_e^2 = 2m_e^2 + 2m_\mu^2 这符合任意 1+23+41+2 \to 3+4 散射过程的普适运动学恒等式 i=14mi2\sum_{i=1}^4 m_i^2。 最终结果为: s+t+u=2me2+2mμ2\boxed{s + t + u = 2m_e^2 + 2m_\mu^2}