29.1

Problem 29.1

schwarzChapter 29

习题 29.1

来源: 第29章, PDF第614页


29.1 The dominant production mechanism for Higgs bosons at LEP was e+eZHe^+ e^- \rightarrow ZH. Calculate the total cross section for this process at tree-level in the Standard Model. How many 100 GeV Higgs bosons would there have been when LEP ran at 206 GeV?

习题 29.1 - 解答


物理背景与过程分析

在标准模型中,LEP 上的希格斯玻色子主要产生机制为 Higgs-strahlung 过程:e+(p1)+e(p2)Z(q)Z(k1)+H(k2)e^+(p_1) + e^-(p_2) \rightarrow Z^*(q) \rightarrow Z(k_1) + H(k_2)。 在树图阶,该过程通过 ss-通道的虚 ZZ 玻色子交换进行。由于 LEP 运行能量 s=206 GeV\sqrt{s} = 206 \text{ GeV} 远高于 ZZ 玻色子质量(mZ91.2 GeVm_Z \approx 91.2 \text{ GeV}),我们可以忽略 ZZ 玻色子传播子中的宽度 ΓZ\Gamma_Z。同时,电子质量 mem_e 极小,在极高能下可完全忽略(me0m_e \to 0)。

1. 矩阵元与截面推导

顶点与传播子规则:

  1. e+eZe^+e^-Z 顶点:igcosθWγμ(cVcAγ5)-i \frac{g}{\cos\theta_W} \gamma^\mu (c_V - c_A \gamma^5),其中 cV=14+sin2θWc_V = -\frac{1}{4} + \sin^2\theta_WcA=14c_A = -\frac{1}{4}。为书写简便,记 gZ=gcosθWg_Z = \frac{g}{\cos\theta_W}
  2. ZZHZZH 顶点:igZmZgνρi g_Z m_Z g^{\nu\rho}
  3. ZZ 玻色子传播子:i(gμνqμqν/mZ2)smZ2\frac{-i(g_{\mu\nu} - q_\mu q_\nu/m_Z^2)}{s - m_Z^2}。由于电子质量忽略不计,由狄拉克方程可知与电子流缩并的 qμq_\mu 项为零,因此传播子可简化为 igμνsmZ2\frac{-i g_{\mu\nu}}{s - m_Z^2}

写出不变矩阵元 M\mathcal{M} iM=vˉ(p1)[igZγμ(cVcAγ5)]u(p2)[igμνsmZ2][igZmZgνρϵρ(k1)]i\mathcal{M} = \bar{v}(p_1) \left[ -i g_Z \gamma^\mu (c_V - c_A \gamma^5) \right] u(p_2) \left[ \frac{-i g_{\mu\nu}}{s - m_Z^2} \right] \left[ i g_Z m_Z g^{\nu\rho} \epsilon^*_\rho(k_1) \right] 化简得: M=gZ2mZsmZ2vˉ(p1)γμ(cVcAγ5)u(p2)ϵμ(k1)\mathcal{M} = \frac{g_Z^2 m_Z}{s - m_Z^2} \bar{v}(p_1) \gamma^\mu (c_V - c_A \gamma^5) u(p_2) \epsilon^*_\mu(k_1)

计算自旋平均的矩阵元平方: 对初态电子自旋求平均(因子 14\frac{1}{4}),对末态 ZZ 玻色子极化求和 ϵμϵν=gμν+k1μk1νmZ2Pμν\sum \epsilon^*_\mu \epsilon_\nu = -g_{\mu\nu} + \frac{k_{1\mu} k_{1\nu}}{m_Z^2} \equiv P_{\mu\nu}M2=14gZ4mZ2(smZ2)2Tr[̸p1γμ(cVcAγ5)̸p2γν(cVcAγ5)]Pμν\overline{|\mathcal{M}|^2} = \frac{1}{4} \frac{g_Z^4 m_Z^2}{(s - m_Z^2)^2} \text{Tr}\left[ \not{p}_1 \gamma^\mu (c_V - c_A \gamma^5) \not{p}_2 \gamma^\nu (c_V - c_A \gamma^5) \right] P_{\mu\nu} 利用 γ5\gamma^5 的反交换关系,迹可化为: Tr[]=Tr[̸p1γμ̸p2γν(cV2+cA2)2cVcA̸p1γμ̸p2γνγ5]\text{Tr}[\dots] = \text{Tr}\left[ \not{p}_1 \gamma^\mu \not{p}_2 \gamma^\nu (c_V^2 + c_A^2) - 2c_V c_A \not{p}_1 \gamma^\mu \not{p}_2 \gamma^\nu \gamma^5 \right] 由于 PμνP_{\mu\nu} 是对称张量,含 γ5\gamma^5 的反对称项在缩并后为零。计算剩余的迹: Tr[]=4(cV2+cA2)(p1μp2ν+p1νp2μ(p1p2)gμν)\text{Tr}[\dots] = 4(c_V^2 + c_A^2) \left( p_1^\mu p_2^\nu + p_1^\nu p_2^\mu - (p_1 \cdot p_2) g^{\mu\nu} \right) 将其与极化张量 PμνP_{\mu\nu} 缩并。在质心系中,p1p2=s/2p_1 \cdot p_2 = s/2。设末态 ZZ 玻色子的动量大小为 k|\vec{k}|,出射角为 θ\theta,则有 EZ2mZ2=k2E_Z^2 - m_Z^2 = |\vec{k}|^2。经过标准运动学代数运算,缩并结果为: (p1μp2ν+p1νp2μs2gμν)(gμν+k1μk1νmZ2)=s(1+k22mZ2sin2θ)\left( p_1^\mu p_2^\nu + p_1^\nu p_2^\mu - \frac{s}{2} g^{\mu\nu} \right) \left( -g_{\mu\nu} + \frac{k_{1\mu} k_{1\nu}}{m_Z^2} \right) = s \left( 1 + \frac{|\vec{k}|^2}{2 m_Z^2} \sin^2\theta \right) 因此,矩阵元平方为: M2=gZ4mZ2s(smZ2)2(cV2+cA2)(1+k22mZ2sin2θ)\overline{|\mathcal{M}|^2} = \frac{g_Z^4 m_Z^2 s}{(s - m_Z^2)^2} (c_V^2 + c_A^2) \left( 1 + \frac{|\vec{k}|^2}{2 m_Z^2} \sin^2\theta \right)

计算总截面: 两体相空间微分为 dΦ2=k16π2sdΩd\Phi_2 = \frac{|\vec{k}|}{16\pi^2 \sqrt{s}} d\Omega。微分散射截面为: dσdΩ=12sk16π2sM2=k32π2s3/2M2\frac{d\sigma}{d\Omega} = \frac{1}{2s} \frac{|\vec{k}|}{16\pi^2 \sqrt{s}} \overline{|\mathcal{M}|^2} = \frac{|\vec{k}|}{32\pi^2 s^{3/2}} \overline{|\mathcal{M}|^2} 对立体角 dΩ=2πsinθdθd\Omega = 2\pi \sin\theta d\theta 积分,利用 11(1+Asin2θ)d(cosθ)=2+43A=2(1+23A)\int_{-1}^1 (1 + A \sin^2\theta) d(\cos\theta) = 2 + \frac{4}{3}A = 2(1 + \frac{2}{3}A)σ=k8πs3/2gZ4mZ2s(smZ2)2(cV2+cA2)(1+k23mZ2)\sigma = \frac{|\vec{k}|}{8\pi s^{3/2}} \frac{g_Z^4 m_Z^2 s}{(s - m_Z^2)^2} (c_V^2 + c_A^2) \left( 1 + \frac{|\vec{k}|^2}{3 m_Z^2} \right) 利用费米常数关系 GF=2g28mW2=2gZ28mZ2G_F = \frac{\sqrt{2} g^2}{8 m_W^2} = \frac{\sqrt{2} g_Z^2}{8 m_Z^2},可得 gZ4=128GF2mZ4g_Z^4 = 128 G_F^2 m_Z^4。代入上式并引入两体运动学函数 λ=(1mH2/smZ2/s)24mH2mZ2/s2\lambda = (1 - m_H^2/s - m_Z^2/s)^2 - 4 m_H^2 m_Z^2 / s^2,其中 k=s2λ1/2|\vec{k}| = \frac{\sqrt{s}}{2} \lambda^{1/2},最终得到树图阶总截面的解析表达式: σ(e+eZH)=2GF2mZ43πs(cV2+cA2)λ1/2(λ+12mZ2/s)(1mZ2/s)2\boxed{ \sigma(e^+e^- \to ZH) = \frac{2 G_F^2 m_Z^4}{3 \pi s} (c_V^2 + c_A^2) \frac{\lambda^{1/2} \left( \lambda + 12 m_Z^2/s \right)}{(1 - m_Z^2/s)^2} }

2. 数值计算与事例数预估

代入标准模型参数:

  • 质心能:s=206 GeV    s=42436 GeV2\sqrt{s} = 206 \text{ GeV} \implies s = 42436 \text{ GeV}^2
  • 玻色子质量:mZ=91.1876 GeVm_Z = 91.1876 \text{ GeV}mH=100 GeVm_H = 100 \text{ GeV}
  • 弱混合角:sin2θW0.2312\sin^2\theta_W \approx 0.2312
  • 费米常数:GF=1.16637×105 GeV2G_F = 1.16637 \times 10^{-5} \text{ GeV}^{-2}

计算耦合常数因子: cV=0.25+0.2312=0.0188,cA=0.25    cV2+cA20.06285c_V = -0.25 + 0.2312 = -0.0188, \quad c_A = -0.25 \implies c_V^2 + c_A^2 \approx 0.06285 计算运动学因子: mZ2/s0.1959,mH2/s0.2356m_Z^2/s \approx 0.1959, \quad m_H^2/s \approx 0.2356 λ=(10.23560.1959)24(0.2356)(0.1959)0.1387    λ1/20.3724\lambda = (1 - 0.2356 - 0.1959)^2 - 4(0.2356)(0.1959) \approx 0.1387 \implies \lambda^{1/2} \approx 0.3724 λ1/2(λ+12mZ2/s)(1mZ2/s)2=0.3724×(0.1387+2.3508)(10.1959)21.433\frac{\lambda^{1/2} (\lambda + 12 m_Z^2/s)}{(1 - m_Z^2/s)^2} = \frac{0.3724 \times (0.1387 + 2.3508)}{(1 - 0.1959)^2} \approx 1.433 计算截面(利用转换关系 1 GeV2=0.389×109 pb1 \text{ GeV}^{-2} = 0.389 \times 10^9 \text{ pb}): σ=2×(1.166×105)2×(91.19)43π×42436×0.06285×1.433×(0.389×109 pb)\sigma = \frac{2 \times (1.166 \times 10^{-5})^2 \times (91.19)^4}{3 \pi \times 42436} \times 0.06285 \times 1.433 \times (0.389 \times 10^9 \text{ pb}) σ1.65 pb\boxed{ \sigma \approx 1.65 \text{ pb} }

计算产生的希格斯玻色子数量: 在 2000 年 LEP 运行于 s206 GeV\sqrt{s} \approx 206 \text{ GeV} 附近时,四个实验组(ALEPH, DELPHI, L3, OPAL)每个探测器收集到的积分亮度约为 L200 pb1\mathcal{L} \approx 200 \text{ pb}^{-1},总积分亮度约为 Ltot800 pb1\mathcal{L}_{\text{tot}} \approx 800 \text{ pb}^{-1}。 产生的 100 GeV 希格斯玻色子总数为: N=σ×Ltot=1.65 pb×800 pb1N = \sigma \times \mathcal{L}_{\text{tot}} = 1.65 \text{ pb} \times 800 \text{ pb}^{-1} N1320 events\boxed{ N \approx 1320 \text{ events} } (注:若仅按单个探测器计算,则约为 330 个事例。)

29.2

Problem 29.2

schwarzChapter 29

习题 29.2

来源: 第29章, PDF第614页


29.2 e+ehadronse^+ e^- \rightarrow \text{hadrons} in the Standard Model. (a) Calculate the rate for the total cross section σtot(e+ehadrons)\sigma_{\text{tot}}(e^+ e^- \rightarrow \text{hadrons}) in the Standard Model at tree-level including both ZZ-boson and photon contributions and their interference. The contribution using photons alone was calculated in Section 26.3. (b) Calculate σtot\sigma_{\text{tot}} at 1-loop. (c) Plot the total cross section as a function of center-of-mass energy showing separately the photon contribution, the ZZ-boson contribution, and their sum. Plot also the sum ignoring interference between the ZZ-boson and photon contributions. When can interference be ignored?

习题 29.2 - 解答


(a) 树图阶总截面 σtot(e+ehadrons)\sigma_{\text{tot}}(e^+ e^- \rightarrow \text{hadrons})

在标准模型中,正负电子湮灭产生强子(即夸克-反夸克对 e+eqqˉe^+ e^- \rightarrow q \bar{q})的树图阶过程包含 ss-通道的光子 (γ\gamma) 交换和 ZZ 玻色子交换。

设电子和夸克与 ZZ 玻色子的矢量和轴矢量耦合常数分别为 cVfc_V^fcAfc_A^f。在标准约定下,相互作用顶点为 igcosθWγμ(cVfcAfγ5)-i \frac{g}{\cos\theta_W} \gamma^\mu (c_V^f - c_A^f \gamma_5),其中: cVf=12T3fQfsin2θW,cAf=12T3fc_V^f = \frac{1}{2} T_3^f - Q_f \sin^2\theta_W, \quad c_A^f = \frac{1}{2} T_3^f 对于电子,Qe=1,T3e=1/2Q_e = -1, T_3^e = -1/2;对于夸克,QfQ_f 为分数电荷,T3f=±1/2T_3^f = \pm 1/2

总散射振幅为 M=Mγ+MZ\mathcal{M} = \mathcal{M}_\gamma + \mathcal{M}_Z,其中: Mγ=e2QeQfs(vˉeγμue)(uˉfγμvf)\mathcal{M}_\gamma = \frac{e^2 Q_e Q_f}{s} (\bar{v}_e \gamma^\mu u_e) (\bar{u}_f \gamma_\mu v_f) MZ=g2cos2θW1sMZ2+iMZΓZ[vˉeγμ(cVecAeγ5)ue][uˉfγμ(cVfcAfγ5)vf]\mathcal{M}_Z = \frac{g^2}{\cos^2\theta_W} \frac{1}{s - M_Z^2 + i M_Z \Gamma_Z} [\bar{v}_e \gamma^\mu (c_V^e - c_A^e \gamma_5) u_e] [\bar{u}_f \gamma_\mu (c_V^f - c_A^f \gamma_5) v_f]

对初始自旋求平均并对末态自旋和颜色(Nc=3N_c = 3)求和,计算平方项 14M2\frac{1}{4} \sum |\mathcal{M}|^2

  1. 纯光子贡献14Mγ2=Nce4Qe2Qf2s22s2(1+cos2θ)\frac{1}{4} \sum |\mathcal{M}_\gamma|^2 = N_c \frac{e^4 Q_e^2 Q_f^2}{s^2} 2s^2 (1 + \cos^2\theta) 对立体角积分 dΩ(1+cos2θ)=16π3\int d\Omega (1+\cos^2\theta) = \frac{16\pi}{3},得到: σγ=4πα23sNcQf2\sigma_\gamma = \frac{4\pi \alpha^2}{3s} N_c Q_f^2

  2. 干涉项贡献: 计算 2Re(MγMZ)2\text{Re}(\mathcal{M}_\gamma^* \mathcal{M}_Z) 时,迹的计算会产生对称部分(正比于 cVecVf(1+cos2θ)c_V^e c_V^f (1+\cos^2\theta))和反对称部分(正比于 cAecAfcosθc_A^e c_A^f \cos\theta)。由于 cosθ\cos\theta 在全立体角积分下为零,轴矢量干涉项对总截面的贡献严格为零。仅保留矢量部分: σγZ=4πα23sNc(2QeQfcVecVfsin2θWcos2θW)s(sMZ2)(sMZ2)2+MZ2ΓZ2\sigma_{\gamma Z} = \frac{4\pi \alpha^2}{3s} N_c \left( \frac{2 Q_e Q_f c_V^e c_V^f}{\sin^2\theta_W \cos^2\theta_W} \right) \frac{s(s - M_Z^2)}{(s - M_Z^2)^2 + M_Z^2 \Gamma_Z^2}

  3. ZZ 玻色子贡献: 同理,纯 ZZ 交换的平方项中,交叉项(正比于 cVcAcosθc_V c_A \cos\theta)积分后消失,仅留下 (cV2+cA2)(c_V^2 + c_A^2) 项: σZ=4πα23sNc((cVe)2+(cAe)2)((cVf)2+(cAf)2)sin4θWcos4θWs2(sMZ2)2+MZ2ΓZ2\sigma_Z = \frac{4\pi \alpha^2}{3s} N_c \frac{((c_V^e)^2 + (c_A^e)^2)((c_V^f)^2 + (c_A^f)^2)}{\sin^4\theta_W \cos^4\theta_W} \frac{s^2}{(s - M_Z^2)^2 + M_Z^2 \Gamma_Z^2}

将上述三项相加,代入 Qe=1Q_e = -1,并对所有运动学允许的夸克味道 ff(通常为 u,d,s,c,bu, d, s, c, b)求和,得到树图阶总截面: σtottree=4πα23sNcf[Qf22QfcVecVfs(sMZ2)sin2θWcos2θWDZ(s)2+((cVe)2+(cAe)2)((cVf)2+(cAf)2)s2sin4θWcos4θWDZ(s)2]\boxed{ \sigma_{\text{tot}}^{\text{tree}} = \frac{4\pi \alpha^2}{3s} N_c \sum_f \left[ Q_f^2 - \frac{2 Q_f c_V^e c_V^f s (s - M_Z^2)}{\sin^2\theta_W \cos^2\theta_W |D_Z(s)|^2} + \frac{((c_V^e)^2 + (c_A^e)^2)((c_V^f)^2 + (c_A^f)^2) s^2}{\sin^4\theta_W \cos^4\theta_W |D_Z(s)|^2} \right] } 其中 DZ(s)2=(sMZ2)2+MZ2ΓZ2|D_Z(s)|^2 = (s - M_Z^2)^2 + M_Z^2 \Gamma_Z^2


(b) 单圈阶 (1-loop) 总截面

在单圈阶,QCD 修正来源于末态夸克的虚胶子交换(顶点修正)和实胶子辐射(e+eqqˉge^+ e^- \rightarrow q \bar{q} g)。 由于在当前能标下夸克质量可以近似为零,夸克的矢量流和轴矢量流是守恒的(手征极限下无反常维度)。因此,QCD 对矢量流(光子和 ZZ 的矢量部分)和轴矢量流(ZZ 的轴矢量部分)的修正因子是完全相同的。

根据标准的微扰 QCD 计算,单圈阶的修正表现为一个整体的乘性因子 (1+αs(s)π)(1 + \frac{\alpha_s(s)}{\pi})。因此,包含 1-loop QCD 修正的总截面为: σtot1-loop=σtottree(1+αs(s)π)\boxed{ \sigma_{\text{tot}}^{\text{1-loop}} = \sigma_{\text{tot}}^{\text{tree}} \left( 1 + \frac{\alpha_s(s)}{\pi} \right) }


(c) 总截面随质心能量的变化及干涉项分析

以下是使用 LaTeX pgfplots 绘制的各部分截面随质心能量 s\sqrt{s} 变化的示意图代码及物理分析:

\begin{center}
\begin{tikzpicture}
\begin{axis}[
    width=12cm, height=8cm,
    xmode=log, ymode=log,
    xmin=10, xmax=200,
    ymin=10, ymax=100000,
    xlabel={Center-of-mass energy $\sqrt{s}$ [GeV]},
    ylabel={Total cross section $\sigma_{\text{tot}}$ [pb]},
    legend pos=south west,
    domain=10:200,
    samples=200,
    grid=both,
    thick
]
% 纯光子贡献 (正比于 1/s)
\addplot[blue, dashed] {318000/(x^2)};
\addlegendentry{Photon ($\gamma$) contribution}

% 纯 Z 玻色子贡献 (Breit-Wigner 共振峰)
\addplot[red, dashed] {217000*(x^2)/((x^2 - 91.2^2)^2 + (91.2*2.5)^2)};
\addlegendentry{$Z$-boson contribution}

% 包含干涉项的总和
\addplot[black, solid] {318000/(x^2) + 217000*(x^2)/((x^2 - 91.2^2)^2 + (91.2*2.5)^2) + 16100*(x^2 - 91.2^2)/((x^2 - 91.2^2)^2 + (91.2*2.5)^2)};
\addlegendentry{Total sum (with interference)}

% 忽略干涉项的总和
\addplot[green, dotted, very thick] {318000/(x^2) + 217000*(x^2)/((x^2 - 91.2^2)^2 + (91.2*2.5)^2)};
\addlegendentry{Sum ignoring interference}
\end{axis}
\end{tikzpicture}
\end{center}

物理分析与干涉项的忽略条件:

  1. 截面行为:纯光子截面随能量呈 1/s1/s 衰减;纯 ZZ 玻色子截面在 s=MZ91.2 GeV\sqrt{s} = M_Z \approx 91.2 \text{ GeV} 处展现出显著的 Breit-Wigner 共振峰。
  2. 干涉项何时可以被忽略? 干涉项在几乎所有能量下都可以被安全地忽略(即图中黑色实线与绿色虚线几乎完全重合)。原因有两点:
    • ZZ 极点附近 (sMZ\sqrt{s} \approx M_Z):干涉项正比于 sMZ2s - M_Z^2。在共振峰处,传播子的实部穿过零点,导致干涉项严格为零。
    • 在远离 ZZ 极点的区域:干涉项的系数正比于电子的矢量耦合常数 cVe=14+sin2θWc_V^e = -\frac{1}{4} + \sin^2\theta_W。由于实验测得的弱混合角 sin2θW0.2315\sin^2\theta_W \approx 0.2315,这使得 cVe0.01851c_V^e \approx -0.0185 \ll 1。电子的矢量耦合由于偶然的参数相消而极小,导致 γZ\gamma-Z 干涉项受到强烈的参数压低。

因此,在计算总截面时,忽略干涉项是一个极其优良的近似。

29.3

Problem 29.3

schwarzChapter 29

习题 29.3

来源: 第29章, PDF第614页


29.3 Higgs decays. (a) Calculate the rate for HbbˉH \rightarrow b\bar{b} in the Standard Model. (b) Calculate the rate for HggH \rightarrow gg in the Standard Model. The dominant contribution to this comes from a triangle loop diagram involving top quarks. (c) Calculate the rate for HγγH \rightarrow \gamma\gamma in the Standard Model. Include contributions both from top loops and from loops of WW bosons. (d) Plot the branching ratios for HbbˉH \rightarrow b\bar{b}, HggH \rightarrow gg and HγγH \rightarrow \gamma\gamma as a function of Higgs mass.

习题 29.3 - 解答


(a) Calculate the rate for HbbˉH \rightarrow b\bar{b} in the Standard Model

在标准模型中,Higgs 玻色子与费米子的汤川耦合(Yukawa coupling)拉格朗日量为: LYukawambvhbˉb\mathcal{L}_{\text{Yukawa}} \supset -\frac{m_b}{v} h \bar{b} b 其中 v=(2GF)1/2246 GeVv = (\sqrt{2} G_F)^{-1/2} \approx 246 \text{ GeV} 是 Higgs 场的真空期望值。由此可得 HbbˉH \rightarrow b\bar{b} 的顶点 Feynman 规则为 imbv-i \frac{m_b}{v}

衰变过程的矩阵元为: M=imbvuˉ(p1)v(p2)\mathcal{M} = -i \frac{m_b}{v} \bar{u}(p_1) v(p_2) 对末态自旋求和并包含底夸克的颜色自由度 Nc=3N_c = 3,矩阵元模方为: M2=Ncmb2v2Tr[(̸p1+mb)(̸p2mb)]=3mb2v2×4(p1p2mb2)\sum |\mathcal{M}|^2 = N_c \frac{m_b^2}{v^2} \text{Tr}\left[ (\not{p}_1 + m_b)(\not{p}_2 - m_b) \right] = 3 \frac{m_b^2}{v^2} \times 4(p_1 \cdot p_2 - m_b^2) 利用运动学关系 pH=p1+p2p_H = p_1 + p_2,可得 mH2=2mb2+2p1p2m_H^2 = 2m_b^2 + 2p_1 \cdot p_2,即 2p1p2=mH22mb22p_1 \cdot p_2 = m_H^2 - 2m_b^2。代入上式: M2=12mb2v2(mH22mb22mb2)=6mb2v2mH2(14mb2mH2)\sum |\mathcal{M}|^2 = \frac{12 m_b^2}{v^2} \left( \frac{m_H^2 - 2m_b^2}{2} - m_b^2 \right) = \frac{6 m_b^2}{v^2} m_H^2 \left( 1 - \frac{4m_b^2}{m_H^2} \right) 两体衰变率的相空间公式为 Γ=12mH18π14mb2mH2M2\Gamma = \frac{1}{2m_H} \frac{1}{8\pi} \sqrt{1 - \frac{4m_b^2}{m_H^2}} \sum |\mathcal{M}|^2。代入矩阵元模方并使用 v2=1/(2GF)v^2 = 1/(\sqrt{2} G_F),得到: Γ(Hbbˉ)=3mb2mH8πv2(14mb2mH2)3/2=32GFmb2mH8π(14mb2mH2)3/2\Gamma(H \rightarrow b\bar{b}) = \frac{3 m_b^2 m_H}{8\pi v^2} \left( 1 - \frac{4m_b^2}{m_H^2} \right)^{3/2} = \frac{3 \sqrt{2} G_F m_b^2 m_H}{8\pi} \left( 1 - \frac{4m_b^2}{m_H^2} \right)^{3/2} 最终结果为: Γ(Hbbˉ)=32GFmb2mH8π(14mb2mH2)3/2\boxed{ \Gamma(H \rightarrow b\bar{b}) = \frac{3 \sqrt{2} G_F m_b^2 m_H}{8\pi} \left( 1 - \frac{4m_b^2}{m_H^2} \right)^{3/2} }


(b) Calculate the rate for HggH \rightarrow gg in the Standard Model

Higgs 玻色子不带色荷,因此 HggH \rightarrow gg 过程在树图级别为零,主要贡献来自于重夸克(主要是顶夸克 tt)的三角圈图。 在 mH2mtm_H \ll 2m_t 的极限下,可以通过计算顶夸克圈图得到 Higgs 与胶子的有效相互作用拉格朗日量(由 QCD 迹反常给出): Leff=αs12πvhGμνaGaμν\mathcal{L}_{\text{eff}} = \frac{\alpha_s}{12\pi v} h G_{\mu\nu}^a G^{a\mu\nu} 由此导出的矩阵元为: M=αs3πvδab(k1μk2ν(k1k2)gμν)ϵ1μ(k1)ϵ2ν(k2)\mathcal{M} = \frac{\alpha_s}{3\pi v} \delta^{ab} \left( k_1^\mu k_2^\nu - (k_1 \cdot k_2) g^{\mu\nu} \right) \epsilon_{1\mu}^*(k_1) \epsilon_{2\nu}^*(k_2) 对极化和颜色求和(注意胶子是全同粒子,相空间需乘以 1/2!1/2!): M2=(αs3πv)2×8×2(k1k2)2=16αs29π2v2(mH22)2=4αs2mH49π2v2\sum |\mathcal{M}|^2 = \left( \frac{\alpha_s}{3\pi v} \right)^2 \times 8 \times 2 (k_1 \cdot k_2)^2 = \frac{16 \alpha_s^2}{9\pi^2 v^2} \left( \frac{m_H^2}{2} \right)^2 = \frac{4 \alpha_s^2 m_H^4}{9\pi^2 v^2} 衰变率为: Γeff=12mH18π12!M2=αs2mH372π3v2=2GFαs2mH372π3\Gamma_{\text{eff}} = \frac{1}{2m_H} \frac{1}{8\pi} \frac{1}{2!} \sum |\mathcal{M}|^2 = \frac{\alpha_s^2 m_H^3}{72\pi^3 v^2} = \frac{\sqrt{2} G_F \alpha_s^2 m_H^3}{72\pi^3} 为了包含精确的顶夸克质量依赖,引入圈图形状因子 A1/2(τt)A_{1/2}(\tau_t),其中 τt=mH2/4mt2\tau_t = m_H^2 / 4m_t^2。在重顶夸克极限下 A1/2(τt)4/3A_{1/2}(\tau_t) \rightarrow 4/3。将有效衰变率除以极限值 4/32=16/9|4/3|^2 = 16/9 并乘上精确形状因子模方,得到完整的衰变率: Γ(Hgg)=GFαs2mH3642π3A1/2(τt)2\boxed{ \Gamma(H \rightarrow gg) = \frac{G_F \alpha_s^2 m_H^3}{64 \sqrt{2} \pi^3} \left| A_{1/2}(\tau_t) \right|^2 } 其中形状因子定义为: A1/2(τ)=2τ2[τ+(τ1)f(τ)]A_{1/2}(\tau) = \frac{2}{\tau^2} \left[ \tau + (\tau - 1) f(\tau) \right] 对于 τ1\tau \le 1(即 mH2mtm_H \le 2m_t),函数 f(τ)=arcsin2τf(\tau) = \arcsin^2\sqrt{\tau}


(c) Calculate the rate for HγγH \rightarrow \gamma\gamma in the Standard Model

与双胶子衰变类似,HγγH \rightarrow \gamma\gamma 也是单圈诱导过程。不同之处在于,除了带电费米子(主要是顶夸克)的贡献外,还有 WW 玻色子圈图的贡献。 费米子圈和 WW 玻色子圈的贡献可以分别用形状因子 A1/2(τf)A_{1/2}(\tau_f)A1(τW)A_1(\tau_W) 参数化。总振幅是这两部分贡献的相干叠加。

根据 QED 迹反常和圈图计算,衰变率的精确表达式为: Γ(Hγγ)=GFα2mH31282π3fNcfQf2A1/2(τf)+A1(τW)2\boxed{ \Gamma(H \rightarrow \gamma\gamma) = \frac{G_F \alpha^2 m_H^3}{128 \sqrt{2} \pi^3} \left| \sum_f N_{cf} Q_f^2 A_{1/2}(\tau_f) + A_1(\tau_W) \right|^2 } 其中:

  • α\alpha 是精细结构常数。
  • NcfN_{cf} 是费米子的颜色数(夸克为 3,轻子为 1),QfQ_f 是费米子的电荷(顶夸克为 +2/3+2/3)。
  • τi=mH2/4mi2\tau_i = m_H^2 / 4m_i^2
  • 费米子形状因子 A1/2(τ)A_{1/2}(\tau) 与 (b) 问中定义相同。
  • 矢量玻色子形状因子 A1(τ)A_1(\tau) 定义为: A1(τ)=1τ2[2τ2+3τ+3(2τ1)f(τ)]A_1(\tau) = - \frac{1}{\tau^2} \left[ 2\tau^2 + 3\tau + 3(2\tau - 1) f(\tau) \right] 物理说明:对于 mH125 GeVm_H \approx 125 \text{ GeV},顶夸克贡献 A1/2(τt)1.38A_{1/2}(\tau_t) \approx 1.38,而 WW 玻色子贡献 A1(τW)8.3A_1(\tau_W) \approx -8.3。两者符号相反,表明费米子圈和规范玻色子圈之间存在强烈的相消干涉(Destructive interference),且 WW 玻色子圈的贡献占据主导地位。

(d) Plot the branching ratios for HbbˉH \rightarrow b\bar{b}, HggH \rightarrow gg and HγγH \rightarrow \gamma\gamma as a function of Higgs mass

分支比(Branching Ratio, BR)定义为特定通道的衰变率除以 Higgs 的总衰变宽度 Γtot\Gamma_{\text{tot}}。总宽度必须包含所有可能的衰变通道(如 WW,ZZ,τ+τ,ccˉWW^*, ZZ^*, \tau^+\tau^-, c\bar{c} 等)。

以下是这三个通道分支比随 Higgs 质量 mHm_H 变化的解析特征与趋势描述(等效于绘图分析):

  1. HbbˉH \rightarrow b\bar{b} 曲线

    • 在低质量区(mH<135 GeVm_H < 135 \text{ GeV}),由于底夸克是该区域内最重的运动学允许费米子,该通道占据绝对主导。
    • mH=100 GeVm_H = 100 \text{ GeV} 时,BR(bbˉ)80%\text{BR}(b\bar{b}) \approx 80\%
    • mH=125 GeVm_H = 125 \text{ GeV} 时,BR(bbˉ)58%\text{BR}(b\bar{b}) \approx 58\%
    • mH>140 GeVm_H > 140 \text{ GeV} 时,随着 HWWH \rightarrow WW^*ZZZZ^* 阈值的打开,规范玻色子通道迅速主导,导致 bbˉb\bar{b} 的分支比呈指数级断崖式下降。
  2. HggH \rightarrow gg 曲线

    • 这是一个单圈过程,但由于强耦合常数 αs\alpha_s 较大且顶夸克汤川耦合极强,其分支比在低质量区相当可观。
    • 100 GeV<mH<135 GeV100 \text{ GeV} < m_H < 135 \text{ GeV} 范围内,BR(gg)\text{BR}(gg) 呈现一个平缓的峰,数值维持在 5%8%5\% \sim 8\% 之间(125 GeV125 \text{ GeV} 时约为 8%8\%)。
    • 同样在 mH>140 GeVm_H > 140 \text{ GeV} 后,受 WWWW^* 通道压制而迅速衰减。
  3. HγγH \rightarrow \gamma\gamma 曲线

    • 依赖于电磁耦合 α\alpha 且为纯圈图过程,因此分支比极小。
    • 曲线形状与 gggg 类似,在 mH125 GeVm_H \approx 125 \text{ GeV} 附近达到峰值,最大分支比约为 0.22%0.22\%
    • 尽管数值微小,但由于双光子末态在探测器中具有极高的能量分辨率和干净的背景,它是发现 125 GeV125 \text{ GeV} Higgs 玻色子的最重要黄金通道之一。

分支比随 mHm_H 变化的特征数据表(近似值,用于构建曲线轮廓)

| mHm_H (GeV) | BR(Hbbˉ)\text{BR}(H \rightarrow b\bar{b}) | BR(Hgg)\text{BR}(H \rightarrow gg) | BR(Hγγ)\text{BR}(H \rightarrow \gamma\gamma) | 主导物理效应 | | :---: | :---: | :---: | :---: | :--- | | 100 | 80%\sim 80\% | 6%\sim 6\% | 0.15%\sim 0.15\% | bbˉb\bar{b} 绝对主导 | | 115 | 70%\sim 70\% | 7%\sim 7\% | 0.20%\sim 0.20\% | WWWW^* 贡献开始显现 | | 125 | 58%\sim 58\% | 8%\sim 8\% | 0.22%\sim 0.22\% | 物理 Higgs 质量点 | | 140 | 35%\sim 35\% | 5%\sim 5\% | 0.15%\sim 0.15\% | WWWW^* 成为主要通道 | | 160 | 4%\sim 4\% | 1%\sim 1\% | 0.02%\sim 0.02\% | WWWW 真实阈值打开,完全压制其他通道 |

29.4

Problem 29.4

schwarzChapter 29

习题 29.4

来源: 第29章, PDF第614,615页


29.4 Partial wave unitarity. (a) Calculate the matrix element for longitudinal WL+WLWL+WLW_L^+ W_L^- \rightarrow W_L^+ W_L^- scattering in the Standard Model. (b) Show that the high-energy behavior of this matrix element is reproduced using the Goldstone boson equivalence theorem.

(c) Does this give a stronger unitarity constraint than the one using WL+ZLWL+ZLW_L^+ Z_L \rightarrow W_L^+ Z_L scattering?

习题 29.4 - 解答


(a) 纵向 WL+WLWL+WLW_L^+ W_L^- \rightarrow W_L^+ W_L^- 散射的矩阵元

在标准模型中,计算纵向规范玻色子在极高能极限(s,tmW2s, t \gg m_W^2)下的散射矩阵元。在高能下,纵向极化矢量可以近似为 ϵLμ(p)pμ/mW\epsilon_L^\mu(p) \approx p^\mu / m_W

树图级别的散射振幅包含两部分贡献:规范玻色子交换(接触图、ss-道和 tt-道 γ/Z\gamma/Z 交换)以及 Higgs 玻色子交换(ss-道和 tt-道)。

  1. 规范部分的贡献: 将 ϵLμpμ/mW\epsilon_L^\mu \approx p^\mu / m_W 代入规范相互作用图中,会产生随能量增长的 O(E4/mW4)\mathcal{O}(E^4/m_W^4)O(E2/mW2)\mathcal{O}(E^2/m_W^2) 项。由于标准模型非阿贝尔规范群的结构,所有 O(E4/mW4)\mathcal{O}(E^4/m_W^4) 的发散项在接触图和 s,ts, t-道交换图之间精确抵消。规范部分剩余的 O(E2/mW2)\mathcal{O}(E^2/m_W^2) 领头项为: Mgaugeg24mW2(s+t)\mathcal{M}_{\text{gauge}} \approx \frac{g^2}{4 m_W^2} (s+t)

  2. Higgs 交换的贡献: 利用 hWWhWW 顶点 igmWgμν-i g m_W g_{\mu\nu},可以计算 ss-道和 tt-道 Higgs 交换的振幅: Mh,s=(igmWgμν)ϵ1μϵ2νismh2(igmWgρσ)ϵ3ρϵ4σg2mW2(p1p2/mW2)(p3p4/mW2)smh2=g24mW2s2smh2\mathcal{M}_{h,s} = (-i g m_W g_{\mu\nu}) \epsilon_1^\mu \epsilon_2^\nu \frac{i}{s - m_h^2} (-i g m_W g_{\rho\sigma}) \epsilon_3^\rho \epsilon_4^\sigma \approx - g^2 m_W^2 \frac{(p_1 \cdot p_2 / m_W^2)(p_3 \cdot p_4 / m_W^2)}{s - m_h^2} = - \frac{g^2}{4 m_W^2} \frac{s^2}{s - m_h^2} 同理,对于 tt-道有: Mh,tg24mW2t2tmh2\mathcal{M}_{h,t} \approx - \frac{g^2}{4 m_W^2} \frac{t^2}{t - m_h^2} 将两者相加,并利用代数恒等式 x2xmh2=x+mh2+mh4xmh2\frac{x^2}{x - m_h^2} = x + m_h^2 + \frac{m_h^4}{x - m_h^2} 展开: Mhg24mW2(s+t+2mh2+mh4smh2+mh4tmh2)\mathcal{M}_h \approx - \frac{g^2}{4 m_W^2} \left( s + t + 2 m_h^2 + \frac{m_h^4}{s - m_h^2} + \frac{m_h^4}{t - m_h^2} \right)

  3. 总矩阵元: 将规范部分与 Higgs 部分相加,随能量增长的 O(E2)\mathcal{O}(E^2)(s+t)(s+t) 再次精确抵消,得到高能极限下有限的矩阵元: M=Mgauge+Mhg2mh24mW2(2+mh2smh2+mh2tmh2)\mathcal{M} = \mathcal{M}_{\text{gauge}} + \mathcal{M}_h \approx - \frac{g^2 m_h^2}{4 m_W^2} \left( 2 + \frac{m_h^2}{s - m_h^2} + \frac{m_h^2}{t - m_h^2} \right) 代入 WW 玻色子质量关系 mW=gv/2m_W = g v / 2,即 g2mh24mW2=mh2v2\frac{g^2 m_h^2}{4 m_W^2} = \frac{m_h^2}{v^2},最终得到: M(WL+WLWL+WL)mh2v2(ssmh2+ttmh2)\boxed{ \mathcal{M}(W_L^+ W_L^- \rightarrow W_L^+ W_L^-) \approx - \frac{m_h^2}{v^2} \left( \frac{s}{s - m_h^2} + \frac{t}{t - m_h^2} \right) }


(b) 使用 Goldstone 玻色子等效定理 (GBET) 验证

Goldstone 玻色子等效定理指出,在 EmWE \gg m_W 时,纵向规范玻色子的散射振幅等效于其对应的 Goldstone 玻色子的散射振幅:M(WL+WLWL+WL)M(π+ππ+π)\mathcal{M}(W_L^+ W_L^- \rightarrow W_L^+ W_L^-) \approx \mathcal{M}(\pi^+ \pi^- \rightarrow \pi^+ \pi^-)

在标准模型中,Higgs 二重态参数化为 H=12(2π+v+h+iπ0)H = \frac{1}{\sqrt{2}} \begin{pmatrix} \sqrt{2} \pi^+ \\ v + h + i \pi^0 \end{pmatrix}。标量势为 V(H)=λ(HHv2/2)2V(H) = \lambda (H^\dagger H - v^2/2)^2。 展开 HHv22=vh+12h2+π+π+12(π0)2H^\dagger H - \frac{v^2}{2} = v h + \frac{1}{2} h^2 + \pi^+ \pi^- + \frac{1}{2} (\pi^0)^2,提取与 π+ππ+π\pi^+ \pi^- \rightarrow \pi^+ \pi^- 相关的相互作用项: Lλ(π+π)22λvhπ+π\mathcal{L} \supset - \lambda (\pi^+ \pi^-)^2 - 2 \lambda v h \pi^+ \pi^- 由此可得四点接触顶点为 4iλ-4i\lambda,三点顶点 hπ+πh\pi^+\pi^-2iλv-2i\lambda v。计算树图振幅(包含接触图、ss-道和 tt-道 Higgs 交换): M(π+ππ+π)=4λ+(2λv)2smh2+(2λv)2tmh2\mathcal{M}(\pi^+ \pi^- \rightarrow \pi^+ \pi^-) = -4 \lambda + \frac{(-2 \lambda v)^2}{s - m_h^2} + \frac{(-2 \lambda v)^2}{t - m_h^2} 利用 Higgs 质量关系 mh2=2λv2m_h^2 = 2 \lambda v^2,可替换 λ=mh22v2\lambda = \frac{m_h^2}{2 v^2}M(π+ππ+π)=2mh2v2mh4/v2smh2mh4/v2tmh2=mh2v2(2+mh2smh2+mh2tmh2)\mathcal{M}(\pi^+ \pi^- \rightarrow \pi^+ \pi^-) = - \frac{2 m_h^2}{v^2} - \frac{m_h^4 / v^2}{s - m_h^2} - \frac{m_h^4 / v^2}{t - m_h^2} = - \frac{m_h^2}{v^2} \left( 2 + \frac{m_h^2}{s - m_h^2} + \frac{m_h^2}{t - m_h^2} \right) 化简后得到: M(π+ππ+π)=mh2v2(ssmh2+ttmh2)\boxed{ \mathcal{M}(\pi^+ \pi^- \rightarrow \pi^+ \pi^-) = - \frac{m_h^2}{v^2} \left( \frac{s}{s - m_h^2} + \frac{t}{t - m_h^2} \right) } 这与 (a) 中直接计算纵向 WW 玻色子得到的高能行为完全一致。


(c) 与 WL+ZLWL+ZLW_L^+ Z_L \rightarrow W_L^+ Z_L 散射的幺正性约束对比

为了得到幺正性约束,我们需要计算 smh2s \gg m_h^2 极限下的第 00 分波振幅 a0=132π11Mdcosθa_0 = \frac{1}{32\pi} \int_{-1}^1 \mathcal{M} d\cos\theta

  1. 对于 WL+WLWL+WLW_L^+ W_L^- \rightarrow W_L^+ W_L^-: 在 smh2s \gg m_h^2 极限下,矩阵元趋于常数 M2mh2v2\mathcal{M} \approx - \frac{2 m_h^2}{v^2}a0(WL+WL)132π11(2mh2v2)dcosθ=mh28πv2a_0(W_L^+ W_L^-) \approx \frac{1}{32\pi} \int_{-1}^1 \left( - \frac{2 m_h^2}{v^2} \right) d\cos\theta = - \frac{m_h^2}{8 \pi v^2}

  2. 对于 WL+ZLWL+ZLW_L^+ Z_L \rightarrow W_L^+ Z_L: 根据 GBET,这等效于 π+π0π+π0\pi^+ \pi^0 \rightarrow \pi^+ \pi^0 散射。从标量势中提取相关项 Lλπ+π(π0)2\mathcal{L} \supset - \lambda \pi^+ \pi^- (\pi^0)^2,得到四点顶点为 2iλ-2i\lambda。由于 π+π0\pi^+ \pi^0 态的电荷守恒,只有 tt-道可以交换 Higgs 玻色子: M(π+π0π+π0)=2λ+(2λv)2tmh2=mh2v2(1+mh2tmh2)=mh2v2ttmh2\mathcal{M}(\pi^+ \pi^0 \rightarrow \pi^+ \pi^0) = -2 \lambda + \frac{(-2 \lambda v)^2}{t - m_h^2} = - \frac{m_h^2}{v^2} \left( 1 + \frac{m_h^2}{t - m_h^2} \right) = - \frac{m_h^2}{v^2} \frac{t}{t - m_h^2}smh2s \gg m_h^2 极限下,除了极小角区域外 ttmh21\frac{t}{t - m_h^2} \approx 1,矩阵元趋于 Mmh2v2\mathcal{M} \approx - \frac{m_h^2}{v^2}a0(WL+ZL)132π11(mh2v2)dcosθ=mh216πv2a_0(W_L^+ Z_L) \approx \frac{1}{32\pi} \int_{-1}^1 \left( - \frac{m_h^2}{v^2} \right) d\cos\theta = - \frac{m_h^2}{16 \pi v^2}

结论: 对比两者分波振幅的绝对值,有 a0(WL+WL)=2a0(WL+ZL)|a_0(W_L^+ W_L^-)| = 2 |a_0(W_L^+ Z_L)|。根据部分波幺正性条件 Re a01/2|\text{Re } a_0| \le 1/2,振幅越大给出的质量上限越严格。因此: 是 (Yes),WL+WLWL+WL 散射给出了比 WL+ZLWL+ZL 更强的幺正性约束。\boxed{ \text{是 (Yes),} W_L^+ W_L^- \rightarrow W_L^+ W_L^- \text{ 散射给出了比 } W_L^+ Z_L \rightarrow W_L^+ Z_L \text{ 更强的幺正性约束。} }

29.5

Problem 29.5

schwarzChapter 29

习题 29.5

来源: 第29章, PDF第615页


29.5 Figure 29.2 includes a number of experimental constraints on the CKM matrix. (a) The parameter ϵK\epsilon_K is what we were calling ϵ\epsilon in Section 29.5.2. Why do the curves marked ϵK\epsilon_K have the shape they do? That is, what combination of CKM elements is ϵK\epsilon_K sensitive to? (b) What could you measure to produce the curves marked Δmd\Delta m_d or Vub|V_{ub}|?

习题 29.5 - 解答


(a) ϵK\epsilon_K 曲线的形状及其敏感的 CKM 矩阵元组合

参数 ϵK\epsilon_K 描述了中性 K 介子系统 (K0Kˉ0K^0-\bar{K}^0) 中的间接 CP 破坏。在标准模型中,K0Kˉ0K^0-\bar{K}^0 混合主要由包含上型夸克 (u,c,tu, c, t) 和 WW 玻色子的盒图 (box diagrams) 贡献。ϵK\epsilon_K 正比于混合矩阵元 M12M_{12} 的虚部: ϵKIm(M12)\epsilon_K \propto \text{Im}(M_{12})

根据 GIM 机制和 OPE (算符乘积展开),M12M_{12} 的表达式依赖于 CKM 矩阵元的组合 λi=VisVid\lambda_i = V_{is}^* V_{id}(其中 i=c,ti = c, t),具体形式为: M12ηccλc2S0(xc)+ηttλt2S0(xt)+2ηctλcλtS0(xc,xt)M_{12} \propto \eta_{cc} \lambda_c^2 S_0(x_c) + \eta_{tt} \lambda_t^2 S_0(x_t) + 2\eta_{ct} \lambda_c \lambda_t S_0(x_c, x_t) 其中 xi=mi2/mW2x_i = m_i^2/m_W^2S0S_0 为 Inami-Lim 函数,ηij\eta_{ij} 为微扰 QCD 修正因子。

为了在 (ρˉ,ηˉ)(\bar{\rho}, \bar{\eta}) 平面上分析其形状,我们引入 Wolfenstein 参数化: Vcdλ,Vcs1    λcλV_{cd} \approx -\lambda, \quad V_{cs} \approx 1 \implies \lambda_c \approx -\lambda Vtd=Aλ3(1ρˉiηˉ),VtsAλ2    λtA2λ5(1ρˉiηˉ)V_{td} = A\lambda^3(1-\bar{\rho}-i\bar{\eta}), \quad V_{ts} \approx -A\lambda^2 \implies \lambda_t \approx -A^2\lambda^5(1-\bar{\rho}-i\bar{\eta})

由于 λc\lambda_c 在此近似下为实数,Im(λc2)0\text{Im}(\lambda_c^2) \approx 0。CP 破坏(虚部)完全来源于包含顶夸克的项:

  1. 交叉项 λcλt\lambda_c \lambda_t 的虚部: Im(λcλt)Im[(λ)(A2λ5(1ρˉiηˉ))]=A2λ6ηˉ\text{Im}(\lambda_c \lambda_t) \approx \text{Im}\left[ (-\lambda) \cdot (-A^2\lambda^5(1-\bar{\rho}-i\bar{\eta})) \right] = -A^2\lambda^6 \bar{\eta}
  2. 顶夸克项 λt2\lambda_t^2 的虚部: Im(λt2)=2Re(λt)Im(λt)2[A2λ5(1ρˉ)][A2λ5ηˉ]=2A4λ10ηˉ(1ρˉ)\text{Im}(\lambda_t^2) = 2\text{Re}(\lambda_t)\text{Im}(\lambda_t) \approx 2[-A^2\lambda^5(1-\bar{\rho})][A^2\lambda^5\bar{\eta}] = -2A^4\lambda^{10}\bar{\eta}(1-\bar{\rho})

将这些虚部代入 ϵK\epsilon_K 的表达式中,我们得到: ϵKηˉ[C1+C2(1ρˉ)]\epsilon_K \propto \bar{\eta} \left[ C_1 + C_2 (1 - \bar{\rho}) \right] 其中 C1A2λ6ηctS0(xc,xt)C_1 \propto A^2\lambda^6 \eta_{ct} S_0(x_c, x_t)C2A4λ10ηttS0(xt)C_2 \propto A^4\lambda^{10} \eta_{tt} S_0(x_t) 均为正的常数。

在实验测量中 ϵK\epsilon_K 是一个定值,因此在 (ρˉ,ηˉ)(\bar{\rho}, \bar{\eta}) 平面上,上述方程的形式为 ηˉ(K1K2ρˉ)=const\bar{\eta}(K_1 - K_2 \bar{\rho}) = \text{const}。这是一个双曲线方程。这就是为什么图中标记为 ϵK\epsilon_K 的限制曲线呈现出双曲线带的形状。

关于 ϵK\epsilon_K 敏感的 CKM 矩阵元组合,由上述推导可知,它直接依赖于 Im(λcλt)\text{Im}(\lambda_c \lambda_t)Im(λt2)\text{Im}(\lambda_t^2)The parameter ϵK is sensitive to the combinations Im(VcsVcdVtsVtd) and Im((VtsVtd)2).\boxed{\text{The parameter } \epsilon_K \text{ is sensitive to the combinations } \text{Im}(V_{cs}^* V_{cd} V_{ts}^* V_{td}) \text{ and } \text{Im}((V_{ts}^* V_{td})^2).}


(b) 产生 Δmd\Delta m_dVub|V_{ub}| 曲线的物理测量

1. Δmd\Delta m_d 曲线 Δmd\Delta m_d 是中性 B 介子系统 (B0Bˉ0B^0-\bar{B}^0,即 BdB_d 介子) 的质量差。与 K 介子系统类似,它由盒图主导。但由于顶夸克质量极大 (mtmc,mum_t \gg m_c, m_u),顶夸克在环图中的贡献占据绝对主导地位。 ΔmdVtbVtd2S0(xt)\Delta m_d \propto |V_{tb} V_{td}^*|^2 S_0(x_t) 由于 Vtb1V_{tb} \approx 1,我们有 ΔmdVtd2\Delta m_d \propto |V_{td}|^2。在 Wolfenstein 参数化下: Vtd2=Aλ3(1ρˉiηˉ)2=A2λ6[(1ρˉ)2+ηˉ2]|V_{td}|^2 = |A\lambda^3(1-\bar{\rho}-i\bar{\eta})|^2 = A^2\lambda^6 \left[ (1-\bar{\rho})^2 + \bar{\eta}^2 \right] 这在 (ρˉ,ηˉ)(\bar{\rho}, \bar{\eta}) 平面上是一个以 (1,0)(1, 0) 为圆心的。 为了产生这条曲线,需要测量 B0Bˉ0B^0-\bar{B}^0 混合的振荡频率。 To produce the Δmd curves, one must measure the time-dependent oscillation frequency of B0Bˉ0 mixing.\boxed{\text{To produce the } \Delta m_d \text{ curves, one must measure the time-dependent oscillation frequency of } B^0 \leftrightarrow \bar{B}^0 \text{ mixing.}} (实验上通常通过测量 B0B^0 介子衰变到特定味最终态(如半轻子衰变 B0D+νB^0 \to D^- \ell^+ \nu)的时间演化率来实现)。

2. Vub|V_{ub}| 曲线 Vub|V_{ub}| 对应于底夸克到上夸克的带电流弱衰变 (buWb \to u W^-)。在 Wolfenstein 参数化下: Vub=Aλ3(ρˉiηˉ)    Vub2=A2λ6(ρˉ2+ηˉ2)V_{ub} = A\lambda^3(\bar{\rho}-i\bar{\eta}) \implies |V_{ub}|^2 = A^2\lambda^6 (\bar{\rho}^2 + \bar{\eta}^2) 这在 (ρˉ,ηˉ)(\bar{\rho}, \bar{\eta}) 平面上是一个以原点 (0,0)(0, 0) 为圆心的。 为了产生这条曲线,需要测量不含粲夸克的 B 介子半轻子衰变过程的衰变率(分支比)。 To produce the Vub curves, one must measure the branching ratios of charmless semileptonic B decays (e.g., Bπν or inclusive BXuν).\boxed{\text{To produce the } |V_{ub}| \text{ curves, one must measure the branching ratios of charmless semileptonic B decays (e.g., } B \to \pi \ell \nu \text{ or inclusive } B \to X_u \ell \nu \text{).}}

29.6

Problem 29.6

schwarzChapter 29

习题 29.6

来源: 第29章, PDF第615页


29.6 Show that with general Dirac and Majorana mass matrices, there are three phases in the PNMS matrix, while if the mass matrix is purely Dirac, there is only one. How many phases are there if the masses are purely Majorana?

习题 29.6 - 解答


物理背景与参数计数基础

在标准模型及其扩展中,轻子部分的带电荷流(Charged Current, CC)相互作用由 PMNS 矩阵(Pontecorvo-Maki-Nakagawa-Sakata matrix)UU 描述: LCC=g2α=e,μ,τi=1,2,3lˉLαγμUαiνLiWμ+h.c.\mathcal{L}_{\text{CC}} = -\frac{g}{\sqrt{2}} \sum_{\alpha=e,\mu,\tau} \sum_{i=1,2,3} \bar{l}_{L\alpha} \gamma^\mu U_{\alpha i} \nu_{Li} W_\mu^- + \text{h.c.} 其中 lLαl_{L\alpha} 是带电轻子的质量本征态,νLi\nu_{Li} 是中微子的质量本征态。

对于 NN 代费米子,PMNS 矩阵 UU 是一个 N×NN \times N 的幺正矩阵。一个一般的 N×NN \times N 幺正矩阵包含 N2N^2 个独立实参数。这些参数可以分解为:

  1. 旋转角(欧拉角)的个数,对应于 N×NN \times N 正交矩阵的参数量:nθ=12N(N1)n_{\theta} = \frac{1}{2}N(N-1)
  2. 复相位的个数:nϕ=N212N(N1)=12N(N+1)n_{\phi} = N^2 - \frac{1}{2}N(N-1) = \frac{1}{2}N(N+1)

在物理上,并非所有相位都是可观测的。我们可以通过对费米子场进行重新定相(Rephasing)来吸收掉一部分无物理意义的相位。具体能吸收多少个相位,取决于中微子质量项的性质(Dirac 型还是 Majorana 型)。本题中代数 N=3N=3


1. 纯 Dirac 质量矩阵的情况 (Purely Dirac Mass Matrix)

如果中微子是纯 Dirac 费米子,其质量项形式为 Lmass=νˉRMDνL+h.c.\mathcal{L}_{\text{mass}} = -\bar{\nu}_R M_D \nu_L + \text{h.c.}。 在这种情况下,我们可以对带电轻子场和中微子场进行独立的全局相位变换,而不改变其对角化的质量项(因为 ψˉψψˉeiαeiαψ=ψˉψ\bar{\psi}\psi \to \bar{\psi}e^{-i\alpha}e^{i\alpha}\psi = \bar{\psi}\psi): lLαeiϕαlLα(α=e,μ,τ)l_{L\alpha} \to e^{i\phi_\alpha} l_{L\alpha} \quad (\alpha = e, \mu, \tau) νLieiθiνLi(i=1,2,3)\nu_{Li} \to e^{i\theta_i} \nu_{Li} \quad (i = 1, 2, 3) 这提供了 NN 个带电轻子相位和 NN 个中微子相位自由度。

在带电荷流相互作用中,PMNS 矩阵的变换为: UαieiϕαUαieiθiU_{\alpha i} \to e^{-i\phi_\alpha} U_{\alpha i} e^{i\theta_i} 这似乎提供了 2N2N 个可以用来消除 UU 中相位的自由度。然而,如果所有带电轻子和中微子都旋转相同的相位(即 ϕe=ϕμ=ϕτ=θ1=θ2=θ3\phi_e=\phi_\mu=\phi_\tau=\theta_1=\theta_2=\theta_3),UαiU_{\alpha i} 保持不变。这个整体的相位对应于轻子数守恒,不能用来消除 PMNS 矩阵内部的相对相位。 因此,真正独立的、可用于吸收 PMNS 矩阵相位的自由度个数为 2N12N - 1

物理相位的个数为幺正矩阵的总相位数减去可吸收的相位数: nphysical=12N(N+1)(2N1)=12(N1)(N2)n_{\text{physical}} = \frac{1}{2}N(N+1) - (2N - 1) = \frac{1}{2}(N-1)(N-2) 代入 N=3N=3nphysical=12(31)(32)=1n_{\text{physical}} = \frac{1}{2}(3-1)(3-2) = 1 这唯一的一个相位就是 Dirac CP 破坏相位 δCP\delta_{\text{CP}}

If the mass matrix is purely Dirac, there is 1 phase in the PMNS matrix.\boxed{\text{If the mass matrix is purely Dirac, there is } 1 \text{ phase in the PMNS matrix.}}


2. 一般的 Dirac 与 Majorana 质量矩阵 (General Dirac and Majorana Mass Matrices)

在最一般的 Seesaw 机制(如 Type-I Seesaw)中,同时存在 Dirac 质量项 MDM_D 和右手中微子的 Majorana 质量项 MRM_R。低能有效理论下,通过对角化 6×66 \times 6 的质量矩阵,我们会得到 3 个极重的中微子和 3 个极轻的中微子。关键在于,这 3 个轻中微子的有效质量项是 Majorana 型 的: Lmasseff=12νLcMeffνL+h.c.\mathcal{L}_{\text{mass}}^{\text{eff}} = -\frac{1}{2} \overline{\nu_L^c} M_{\text{eff}} \nu_L + \text{h.c.} 其中 νLc=CνˉLT\nu_L^c = C \bar{\nu}_L^T 是电荷共轭场。

对于 Majorana 质量项,如果我们尝试对中微子场进行相位变换 νLieiθiνLi\nu_{Li} \to e^{i\theta_i} \nu_{Li},质量项会发生改变: νLicmiνLie2iθiνLicmiνLi\overline{\nu_{Li}^c} m_i \nu_{Li} \to e^{2i\theta_i} \overline{\nu_{Li}^c} m_i \nu_{Li} 为了保证物理基底中 Majorana 质量 mim_i 为实数且为正,我们失去了对中微子场进行任意重新定相的自由度(除了平凡的符号翻转,但这不能消除连续相位)。

因此,我们只能利用带电轻子场的 NN 个相位自由度 lLαeiϕαlLαl_{L\alpha} \to e^{i\phi_\alpha} l_{L\alpha} 来吸收 PMNS 矩阵中的相位。 物理相位的个数变为: nphysical=12N(N+1)N=12N(N1)n_{\text{physical}} = \frac{1}{2}N(N+1) - N = \frac{1}{2}N(N-1) 代入 N=3N=3nphysical=12(3)(2)=3n_{\text{physical}} = \frac{1}{2}(3)(2) = 3 这 3 个相位通常被参数化为 1 个 Dirac 相位 δ\delta 和 2 个 Majorana 相位 α1,α2\alpha_1, \alpha_2

With general Dirac and Majorana mass matrices, there are 3 phases in the PMNS matrix.\boxed{\text{With general Dirac and Majorana mass matrices, there are } 3 \text{ phases in the PMNS matrix.}}


3. 纯 Majorana 质量矩阵的情况 (Purely Majorana Mass Matrix)

如果质量矩阵是纯 Majorana 的(例如,标准模型中没有右手中微子,轻中微子质量完全由维度为 5 的 Weinberg 算符 cΛ(LH~)(LH~)\frac{c}{\Lambda}(L \tilde{H})(L \tilde{H}) 自发破缺后直接生成),低能下的物理状态与上述一般情况完全相同:轻中微子依然是 Majorana 费米子

由于中微子本质上是 Majorana 粒子,其质量项 12νLcMνL\frac{1}{2} \overline{\nu_L^c} M \nu_L 同样打破了轻子数守恒,禁止了对 νLi\nu_{Li} 场的重新定相。 吸收相位的逻辑与第 2 部分完全一致:只有带电轻子的 N=3N=3 个相位可以被吸收。

计算过程不变: nphysical=12N(N+1)N=12(3)(4)3=3n_{\text{physical}} = \frac{1}{2}N(N+1) - N = \frac{1}{2}(3)(4) - 3 = 3

If the masses are purely Majorana, there are also 3 phases in the PMNS matrix.\boxed{\text{If the masses are purely Majorana, there are also } 3 \text{ phases in the PMNS matrix.}}

29.7

Problem 29.7

schwarzChapter 29

习题 29.7

来源: 第29章, PDF第615页


29.7 Neutrino oscillations. (a) Neutrinos are produced in the Sun predominantly through the reaction p+p+ed+νep + p + e^- \rightarrow d + \nu_e. What is the momentum of the neutrinos produced this way? (b) Consider a two-neutrino flavor system. The mass eigenstates evolve in time as

ν1=eiE1t(cosθνe+sinθνμ),ν2=eiE2t(sinθνe+cosθνμ),(29.117)\begin{aligned} |\nu_1\rangle &= e^{-iE_1t} \Big( \cos \theta |\nu_e\rangle + \sin \theta |\nu_\mu\rangle \Big), \\ |\nu_2\rangle &= e^{-iE_2t} \Big( -\sin \theta |\nu_e\rangle + \cos \theta |\nu_\mu\rangle \Big), \end{aligned} \tag{29.117}
(29.118)\tag{29.118}

where θ\theta is the mixing angle. Show that the probability of finding a solar neutrino as an electron neutrino after a time TT is given by

P=1sin2(2θ)sin2(E2E1)T2.(29.119)P = 1 - \sin^2(2\theta) \sin^2 \frac{(E_2 - E_1)T}{2}. \tag{29.119}

(c) Take the high-energy limit EmνE \gg m_\nu to show that the probability of finding a solar neutrino with energy EE as an electron neutrino at a distance LL is given by

P=1sin2(2θ)sin2Δm2L4E.(29.120)P = 1 - \sin^2(2\theta) \sin^2 \frac{\Delta m^2 L}{4E}. \tag{29.120}

(d) How far should you put your detector from a reactor producing 4MeV\sim 4\text{MeV} neutrinos assuming Δm2=7.5×105 eV2\Delta m^2 = 7.5 \times 10^{-5} \text{ eV}^2 to see the largest effect?

习题 29.7 - 解答


(a) 在太阳中,通过 p+p+ed+νep + p + e^- \rightarrow d + \nu_e(即 pep 反应)产生的电子中微子是单能的。由于太阳核心的温度对应的热能(keV\sim \text{keV})远小于参与反应粒子的静止质量,可以近似认为初始状态的质子和电子处于静止状态。

根据能量守恒和动量守恒定律: Ei=2mp+me=Ed+EνE_i = 2m_p + m_e = E_d + E_\nu pi=0=pd+pν    pd=pν=p\vec{p}_i = 0 = \vec{p}_d + \vec{p}_\nu \implies |\vec{p}_d| = |\vec{p}_\nu| = p 其中 mp938.27 MeV/c2m_p \approx 938.27 \text{ MeV}/c^2 为质子质量,me0.511 MeV/c2m_e \approx 0.511 \text{ MeV}/c^2 为电子质量,md1875.61 MeV/c2m_d \approx 1875.61 \text{ MeV}/c^2 为氘核质量。中微子质量极小,可近似为零,因此其能量 Eν=pcE_\nu = p c。氘核的能量为 Ed=md2c4+p2c2E_d = \sqrt{m_d^2 c^4 + p^2 c^2}

代入能量守恒方程: 2mpc2+mec2=md2c4+p2c2+pc2m_p c^2 + m_e c^2 = \sqrt{m_d^2 c^4 + p^2 c^2} + p c 令初始总能量 Ei=2mpc2+mec21877.05 MeVE_i = 2m_p c^2 + m_e c^2 \approx 1877.05 \text{ MeV},移项并平方: (Eipc)2=md2c4+p2c2(E_i - p c)^2 = m_d^2 c^4 + p^2 c^2 Ei22Eipc+p2c2=md2c4+p2c2E_i^2 - 2E_i p c + p^2 c^2 = m_d^2 c^4 + p^2 c^2 解得中微子的动量 pp 为: p=Ei2md2c42Eicp = \frac{E_i^2 - m_d^2 c^4}{2E_i c} 代入数值计算: p=(1877.05)2(1875.61)22×1877.05 MeV/c5403.83754.1 MeV/c1.44 MeV/cp = \frac{(1877.05)^2 - (1875.61)^2}{2 \times 1877.05} \text{ MeV}/c \approx \frac{5403.8}{3754.1} \text{ MeV}/c \approx 1.44 \text{ MeV}/c p1.44 MeV/c\boxed{p \approx 1.44 \text{ MeV}/c}

(b) 根据题目给出的质量本征态随时间的演化关系,在 t=0t=0 时刻有: ν1(0)=cosθνe+sinθνμ|\nu_1(0)\rangle = \cos \theta |\nu_e\rangle + \sin \theta |\nu_\mu\rangle ν2(0)=sinθνe+cosθνμ|\nu_2(0)\rangle = -\sin \theta |\nu_e\rangle + \cos \theta |\nu_\mu\rangle 将上述关系反转,可以将初始的电子中微子味本征态 νe|\nu_e\rangle 用质量本征态表示: νe=cosθν1(0)sinθν2(0)|\nu_e\rangle = \cos \theta |\nu_1(0)\rangle - \sin \theta |\nu_2(0)\rangle 经过时间 TT 后,质量本征态各自累积了动力学相位,系统状态演化为: ν(T)=cosθeiE1Tν1(0)sinθeiE2Tν2(0)|\nu(T)\rangle = \cos \theta e^{-iE_1T} |\nu_1(0)\rangle - \sin \theta e^{-iE_2T} |\nu_2(0)\rangleν1(0)|\nu_1(0)\rangleν2(0)|\nu_2(0)\rangle 的表达式代入上式,提取 νe|\nu_e\rangle 的系数即可得到在时间 TT 找到电子中微子的概率幅 A\mathcal{A}A=νeν(T)=cos2θeiE1T+sin2θeiE2T\mathcal{A} = \langle \nu_e | \nu(T) \rangle = \cos^2 \theta e^{-iE_1T} + \sin^2 \theta e^{-iE_2T} 对应的概率 PP 为概率幅的模平方: P=A2=(cos2θeiE1T+sin2θeiE2T)(cos2θeiE1T+sin2θeiE2T)=cos4θ+sin4θ+cos2θsin2θ(ei(E2E1)T+ei(E2E1)T)=(cos2θ+sin2θ)22cos2θsin2θ+2cos2θsin2θcos((E2E1)T)=12cos2θsin2θ[1cos((E2E1)T)]\begin{aligned} P &= |\mathcal{A}|^2 = \left( \cos^2 \theta e^{-iE_1T} + \sin^2 \theta e^{-iE_2T} \right) \left( \cos^2 \theta e^{iE_1T} + \sin^2 \theta e^{iE_2T} \right) \\ &= \cos^4 \theta + \sin^4 \theta + \cos^2 \theta \sin^2 \theta \left( e^{i(E_2 - E_1)T} + e^{-i(E_2 - E_1)T} \right) \\ &= (\cos^2 \theta + \sin^2 \theta)^2 - 2\cos^2 \theta \sin^2 \theta + 2\cos^2 \theta \sin^2 \theta \cos((E_2 - E_1)T) \\ &= 1 - 2\cos^2 \theta \sin^2 \theta \left[ 1 - \cos((E_2 - E_1)T) \right] \end{aligned} 利用三角恒等式 2cosθsinθ=sin(2θ)2\cos\theta\sin\theta = \sin(2\theta) 以及 1cosx=2sin2(x/2)1 - \cos x = 2\sin^2(x/2),可得: P=1sin2(2θ)sin2(E2E1)T2\boxed{P = 1 - \sin^2(2\theta) \sin^2 \frac{(E_2 - E_1)T}{2}}

(c) 在极端相对论极限(高能极限)下,EmνE \gg m_\nu,中微子的能量可以进行泰勒展开近似: Ei=p2+mi2p+mi22pE_i = \sqrt{p^2 + m_i^2} \approx p + \frac{m_i^2}{2p} 由于 EpE \approx p(自然单位制 c=1c=1),能量差可以近似为: E2E1m22m122E=Δm22EE_2 - E_1 \approx \frac{m_2^2 - m_1^2}{2E} = \frac{\Delta m^2}{2E} 同时,中微子以接近光速运动,在时间 TT 内传播的距离 LTL \approx T。将这些近似代入 (b) 中得到的相位项: (E2E1)T2Δm22EL2=Δm2L4E\frac{(E_2 - E_1)T}{2} \approx \frac{\Delta m^2}{2E} \frac{L}{2} = \frac{\Delta m^2 L}{4E} 因此,在距离 LL 处探测到电子中微子的概率公式化为: P=1sin2(2θ)sin2Δm2L4E\boxed{P = 1 - \sin^2(2\theta) \sin^2 \frac{\Delta m^2 L}{4E}}

(d) 要观察到最大的中微子振荡效应,即让电子中微子存活概率 PP 达到极小值,需要满足振荡项的相位为 π/2\pi/2 的奇数倍。对于第一次最大效应(即最短距离),有: Δm2L4E=π2    L=2πEΔm2\frac{\Delta m^2 L}{4E} = \frac{\pi}{2} \implies L = \frac{2\pi E}{\Delta m^2} 恢复国际单位制中的 \hbarcc,相位项无量纲化要求: L=2πEcΔm2c4L = \frac{2\pi E \hbar c}{\Delta m^2 c^4} 已知参数为:E=4 MeV=4×106 eVE = 4 \text{ MeV} = 4 \times 10^6 \text{ eV}Δm2c4=7.5×105 eV2\Delta m^2 c^4 = 7.5 \times 10^{-5} \text{ eV}^2,以及普朗克常数与光速的乘积 c197.3 MeVfm=1.973×107 eVm\hbar c \approx 197.3 \text{ MeV}\cdot\text{fm} = 1.973 \times 10^{-7} \text{ eV}\cdot\text{m}。 代入数值计算探测器的最佳放置距离 LLL=2π×(4×106 eV)×(1.973×107 eVm)7.5×105 eV2L = \frac{2\pi \times (4 \times 10^6 \text{ eV}) \times (1.973 \times 10^{-7} \text{ eV}\cdot\text{m})}{7.5 \times 10^{-5} \text{ eV}^2} L4.9597.5×105 m66120 m66.1 kmL \approx \frac{4.959}{7.5 \times 10^{-5}} \text{ m} \approx 66120 \text{ m} \approx 66.1 \text{ km} L66.1 km\boxed{L \approx 66.1 \text{ km}}

29.8

Problem 29.8

schwarzChapter 29

习题 29.8

来源: 第29章, PDF第615页


29.8 Show that when you integrate out the right-handed neutrinos in Eq. (29.63), a dimension-5 operator like that in Eq. (29.65) results. What is the exact relationship between MijM_{ij} and M~ij\widetilde{M}_{ij}?

Referenced Equations:

Equation (29.63):

Lmass=YijeLˉiHeRjYijνLˉiH~νRjiMij(νRi)cνRj+h.c.(29.63)\mathcal{L}_{\text{mass}} = -Y_{ij}^e \bar{L}^i H e_R^j - Y_{ij}^\nu \bar{L}^i \tilde{H} \nu_R^j - i M_{ij} (\nu_R^i)^c \nu_R^j + h.c. \tag{29.63}

Equation (29.65):

Ldim-5=M~ij(LˉiH~)(H~Lj).(29.65)\mathcal{L}_{\text{dim-5}} = -\widetilde{M}_{ij} (\bar{L}^i \tilde{H}) (\tilde{H} L^j)^\dagger . \tag{29.65}

习题 29.8 - 解答


题目分析与物理背景

本题要求通过积分掉(integrate out)重质量的右手中微子 νR\nu_R(即应用第一类跷跷板机制,Type-I Seesaw Mechanism),从包含狄拉克质量项和马约拉纳质量项的拉格朗日量中推导出低能有效理论下的维度为5的温伯格算符(Weinberg operator),并求出有效质量矩阵 M~ij\widetilde{M}_{ij} 与右手中微子马约拉纳质量矩阵 MijM_{ij} 之间的精确关系。

注:题干公式中的 iMij(νRi)cνRj-i M_{ij} (\nu_R^i)^c \nu_R^j(H~Lj)(\tilde{H} L^j)^\dagger 存在常见的排版或OCR识别错误。根据 Schwartz《Quantum Field Theory and the Standard Model》原书第29章 Eq. (29.63) 和 Eq. (29.65),标准的马约拉纳质量项应为 12MijνRciνRj-\frac{1}{2} M_{ij} \overline{\nu_R^{ci}} \nu_R^j,维度为5的算符应为 M~ij(LˉiH~)(H~TLcj)-\widetilde{M}_{ij} (\bar{L}^i \tilde{H}) (\tilde{H}^T L^{cj})。以下推导将基于标准物理形式进行。

推导过程

在远低于右手中微子质量标度(EME \ll M)的低能有效理论中,我们可以忽略 νR\nu_R 的动能项,直接通过其经典运动方程(Euler-Lagrange 方程)将其积分掉。

提取拉格朗日量中与 νR\nu_R 相关的项: LνR=YijνLˉiH~νRj12MijνRciνRj+h.c.\mathcal{L}_{\nu_R} = - Y_{ij}^\nu \bar{L}^i \tilde{H} \nu_R^j - \frac{1}{2} M_{ij} \overline{\nu_R^{ci}} \nu_R^j + h.c.

利用电荷共轭关系 νRc=CνˉRT\nu_R^c = C \bar{\nu}_R^T 以及 νRc=νRTC\overline{\nu_R^c} = \nu_R^T C,可以将马约拉纳质量项写为: LνR=LˉH~YννR12νRTCMνR+h.c.\mathcal{L}_{\nu_R} = - \bar{L} \tilde{H} Y^\nu \nu_R - \frac{1}{2} \nu_R^T C M \nu_R + h.c. 其中 CC 是电荷共轭矩阵,满足 CT=CC^T = -CC1=CC^{-1} = -CMM 是对称的味空间质量矩阵(MT=MM^T = M)。

νR\nu_R 进行变分以求解其运动方程。由于 νR\nu_R 是 Grassmann 数,变分时需注意交换带来的符号变化: δ(νRTCMνR)=δνRTCMνR+νRTCMδνR\delta \left( \nu_R^T C M \nu_R \right) = \delta \nu_R^T C M \nu_R + \nu_R^T C M \delta \nu_R 利用 Grassmann 数的转置性质,第二项可以改写为: (νRTCMδνR)T=δνRTMTCTνR=δνRTMCνR=δνRTCMνR(\nu_R^T C M \delta \nu_R)^T = -\delta \nu_R^T M^T C^T \nu_R = \delta \nu_R^T M C \nu_R = \delta \nu_R^T C M \nu_R (这里用到了 MT=MM^T = MCT=CC^T = -C)。因此: δ(νRTCMνR)=2δνRTCMνR\delta \left( \nu_R^T C M \nu_R \right) = 2 \delta \nu_R^T C M \nu_R

拉格朗日量对 νR\nu_R 的变分为: δLνR=LˉH~YνδνRδνRTCMνR=δνRT[(Yν)T(LˉH~)T+CMνR]=0\delta \mathcal{L}_{\nu_R} = - \bar{L} \tilde{H} Y^\nu \delta \nu_R - \delta \nu_R^T C M \nu_R = - \delta \nu_R^T \left[ (Y^\nu)^T (\bar{L} \tilde{H})^T + C M \nu_R \right] = 0 由此得到 νR\nu_R 的运动方程: CMνR=(Yν)T(LˉH~)TC M \nu_R = - (Y^\nu)^T (\bar{L} \tilde{H})^T 左乘 C1=CC^{-1} = -C 并假设 MM 可逆,解得 νR\nu_RMνR=C(Yν)T(LˉH~)T    νR=M1C(Yν)T(LˉH~)TM \nu_R = C (Y^\nu)^T (\bar{L} \tilde{H})^T \implies \nu_R = M^{-1} C (Y^\nu)^T (\bar{L} \tilde{H})^T

将求得的 νR\nu_R 代回原拉格朗日量中。首先利用运动方程 νRTCM=LˉH~Yν\nu_R^T C M = - \bar{L} \tilde{H} Y^\nu 化简质量项: 12νRTCMνR=12(LˉH~Yν)νR=12LˉH~YννR- \frac{1}{2} \nu_R^T C M \nu_R = - \frac{1}{2} \left( - \bar{L} \tilde{H} Y^\nu \right) \nu_R = \frac{1}{2} \bar{L} \tilde{H} Y^\nu \nu_R 代入后,有效拉格朗日量化简为: Leff=LˉH~YννR+12LˉH~YννR+h.c.=12LˉH~YννR+h.c.\mathcal{L}_{\text{eff}} = - \bar{L} \tilde{H} Y^\nu \nu_R + \frac{1}{2} \bar{L} \tilde{H} Y^\nu \nu_R + h.c. = - \frac{1}{2} \bar{L} \tilde{H} Y^\nu \nu_R + h.c.

接着,将 νR\nu_R 的显式解代入上式: Leff=12LˉH~YνM1C(Yν)T(LˉH~)T+h.c.\mathcal{L}_{\text{eff}} = - \frac{1}{2} \bar{L} \tilde{H} Y^\nu M^{-1} C (Y^\nu)^T (\bar{L} \tilde{H})^T + h.c. 利用电荷共轭旋量的定义 Lc=CLˉTL^c = C \bar{L}^T,我们可以将 C(LˉH~)TC (\bar{L} \tilde{H})^T 改写为: C(LˉH~)T=C(H~TLˉT)=H~T(CLˉT)=H~TLcC (\bar{L} \tilde{H})^T = C (\tilde{H}^T \bar{L}^T) = \tilde{H}^T (C \bar{L}^T) = \tilde{H}^T L^c 因此,有效拉格朗日量化为: Leff=12(LˉH~)YνM1(Yν)T(H~TLc)+h.c.\mathcal{L}_{\text{eff}} = - \frac{1}{2} (\bar{L} \tilde{H}) Y^\nu M^{-1} (Y^\nu)^T (\tilde{H}^T L^c) + h.c. 恢复代数指标 i,ji, jLeff=12(YνM1(Yν)T)ij(LˉiH~)(H~TLcj)+h.c.\mathcal{L}_{\text{eff}} = - \frac{1}{2} (Y^\nu M^{-1} (Y^\nu)^T)_{ij} (\bar{L}^i \tilde{H}) (\tilde{H}^T L^{cj}) + h.c.

结论与关系式

将上述推导得到的 Leff\mathcal{L}_{\text{eff}} 与目标维度为5的算符(Weinberg operator)进行对比: Ldim-5=M~ij(LˉiH~)(H~TLcj)+h.c.\mathcal{L}_{\text{dim-5}} = -\widetilde{M}_{ij} (\bar{L}^i \tilde{H}) (\tilde{H}^T L^{cj}) + h.c. 可以清晰地看出,积分掉右手中微子后确实自然产生了一个维度为5的算符。对比系数,我们得到有效质量矩阵 M~ij\widetilde{M}_{ij} 与右手中微子质量矩阵 MijM_{ij} 之间的精确关系为:

M~ij=12(YνM1(Yν)T)ij\boxed{ \widetilde{M}_{ij} = \frac{1}{2} \left( Y^\nu M^{-1} (Y^\nu)^T \right)_{ij} }
29.9

Problem 29.9

schwarzChapter 29

习题 29.9

来源: 第29章, PDF第615页


29.9 Show that when multiple generations are rotated, then the θ\theta angle shifts by argdet(RL)\arg \det (R^\dagger L).

习题 29.9 - 解答


习题分析与物理背景

在标准模型或一般的规范场论中,多代费米子的手征变换会由于量子反常(即阿德勒-贝尔-杰克威反常,ABJ Anomaly)导致路径积分测度发生改变。这种测度的变化等效于在拉格朗日量中引入了一个拓扑项,从而导致强 CP 破坏项中的 θ\theta 角发生平移。

同时,物理上可观测的强 CP 相位 θˉ\bar{\theta} 是由规范场的 θ\theta 角和费米子质量矩阵 MM 的行列式相位共同决定的,即 θˉ=θ+argdetM\bar{\theta} = \theta + \arg \det M。由于 θˉ\bar{\theta} 是物理可观测的,它必须在手征变换下保持不变。我们可以通过 Fujikawa 路径积分测度方法或物理不变量 θˉ\bar{\theta} 的守恒性来推导 θ\theta 角的平移量。


推导过程

方法一:基于 Fujikawa 路径积分测度方法

设理论中包含多代(多味)狄拉克费米子 ψ\psi。对左手和右手费米子场分别作幺正的味空间旋转(手征变换): ψLψL=LψL\psi_L \to \psi_L' = L \psi_L ψRψR=RψR\psi_R \to \psi_R' = R \psi_R 其中 LLRR 是作用在代空间上的幺正矩阵,且 ψL,R=PL,Rψ=1γ52ψ\psi_{L,R} = P_{L,R} \psi = \frac{1 \mp \gamma_5}{2} \psi

我们可以将该变换写为整体的狄拉克旋量形式: ψUψ=(LPL+RPR)ψ=exp(iαV+iαAγ5)ψ\psi \to U \psi = (L P_L + R P_R) \psi = \exp(i \alpha_V + i \alpha_A \gamma_5) \psi 其中 αV\alpha_VαA\alpha_A 是代空间中的厄米矩阵。比较左右手投影部分,可以得到: L=exp[i(αVαA)]L = \exp[i(\alpha_V - \alpha_A)] R=exp[i(αV+αA)]R = \exp[i(\alpha_V + \alpha_A)]

对上述两式在代空间取行列式: detL=exp[iTr(αVαA)]\det L = \exp[i \text{Tr}(\alpha_V - \alpha_A)] detR=exp[iTr(αV+αA)]\det R = \exp[i \text{Tr}(\alpha_V + \alpha_A)] 由于 LLRR 是幺正矩阵,它们的行列式是纯相位,因此可以提取出轴矢量生成元的迹: Tr(2αA)=Tr(αV+αA)Tr(αVαA)=ilndetR(ilndetL)\text{Tr}(2\alpha_A) = \text{Tr}(\alpha_V + \alpha_A) - \text{Tr}(\alpha_V - \alpha_A) = -i \ln \det R - (-i \ln \det L) Tr(2αA)=argdetRargdetL\text{Tr}(2\alpha_A) = \arg \det R - \arg \det L

根据 Fujikawa 方法,在手征变换下,费米子的路径积分测度 DψDψˉ\mathcal{D}\psi \mathcal{D}\bar{\psi} 并不具有不变性,而是会产生一个雅可比行列式(Jacobian),该雅可比行列式给出了手征反常: DψDψˉDψDψˉexp(id4xg232π2FμνaF~aμνTr(2αA))\mathcal{D}\psi \mathcal{D}\bar{\psi} \to \mathcal{D}\psi \mathcal{D}\bar{\psi} \exp\left( -i \int d^4x \frac{g^2}{32\pi^2} F_{\mu\nu}^a \tilde{F}^{a\mu\nu} \text{Tr}(2\alpha_A) \right) 测度的这一变化等效于在有效拉格朗日量中增加了一项: δL=Tr(2αA)g232π2FμνaF~aμν\delta \mathcal{L} = - \text{Tr}(2\alpha_A) \frac{g^2}{32\pi^2} F_{\mu\nu}^a \tilde{F}^{a\mu\nu}

规范场论中原本的 θ\theta 项为 Lθ=θg232π2FμνaF~aμν\mathcal{L}_\theta = \theta \frac{g^2}{32\pi^2} F_{\mu\nu}^a \tilde{F}^{a\mu\nu}。因此,测度变化等效于 θ\theta 角发生了如下平移: Δθ=Tr(2αA)=(argdetRargdetL)=argdetLargdetR\Delta \theta = - \text{Tr}(2\alpha_A) = - (\arg \det R - \arg \det L) = \arg \det L - \arg \det R 利用复数辐角和幺正矩阵的性质 argdetR1=argdetR=argdetR\arg \det R^{-1} = \arg \det R^\dagger = - \arg \det R,上式可化简为: Δθ=arg(detR)+arg(detL)=arg(detRdetL)=argdet(RL)\Delta \theta = \arg(\det R^\dagger) + \arg(\det L) = \arg(\det R^\dagger \det L) = \arg \det (R^\dagger L)


方法二:基于物理可观测量 θˉ\bar{\theta} 的不变性

考虑多代费米子的质量项: Lm=ψˉLMψRψˉRMψL\mathcal{L}_m = - \bar{\psi}_L M \psi_R - \bar{\psi}_R M^\dagger \psi_L 在手征变换 ψLLψL\psi_L \to L \psi_LψRRψR\psi_R \to R \psi_R 下,质量项变为: ψˉLLMRψR+h.c.- \bar{\psi}_L L^\dagger M R \psi_R + \text{h.c.} 这意味着质量矩阵发生了变换:MM=LMRM \to M' = L^\dagger M R。 变换后质量矩阵行列式的相位变为: argdetM=argdet(LMR)=argdetMargdetL+argdetR\arg \det M' = \arg \det(L^\dagger M R) = \arg \det M - \arg \det L + \arg \det R

物理上唯一可观测的强 CP 破坏参数是 θˉ\bar{\theta},它定义为: θˉ=θ+argdetM\bar{\theta} = \theta + \arg \det M 由于 θˉ\bar{\theta} 是物理不变量,在手征变换前后必须保持相等,即 θˉ=θˉ\bar{\theta}' = \bar{\theta}θ+argdetM=θ+argdetM\theta' + \arg \det M' = \theta + \arg \det M 代入 argdetM\arg \det M' 的表达式: θ+argdetMargdetL+argdetR=θ+argdetM\theta' + \arg \det M - \arg \det L + \arg \det R = \theta + \arg \det M 解得新的 θ\theta 角为: θ=θ+argdetLargdetR\theta' = \theta + \arg \det L - \arg \det R 因此,θ\theta 角的平移量为: Δθ=θθ=argdetLargdetR=argdet(RL)\Delta \theta = \theta' - \theta = \arg \det L - \arg \det R = \arg \det(R^\dagger L)

结论

无论是通过微观的路径积分测度反常计算,还是通过宏观的物理不变量守恒推导,多代费米子在经历手征旋转时,θ\theta 角的平移量均严格为 argdet(RL)\arg \det (R^\dagger L)

Δθ=argdet(RL)\boxed{\Delta \theta = \arg \det (R^\dagger L)}