17.1

Problem 17.1

peskinChapter 17

习题 17.1

来源: 第17章, PDF第595,596页


17.1 Two-loop renormalization group relations.

(a) In higher orders of perturbation theory, the expression for the QCD β\beta function will be a series

β(g)=b0(4π)2g3b1(4π)4g5b2(4π)6g7+.\beta(g) = -\frac{b_0}{(4\pi)^2}g^3 - \frac{b_1}{(4\pi)^4}g^5 - \frac{b_2}{(4\pi)^6}g^7 + \cdots.

Integrate the renormalization group equation and show that the running coupling constant is now given by

αs(Q2)=4πb0[1log(Q2/Λ2)b1b02loglog(Q2/Λ2)(log(Q2/Λ2))2+],\alpha_s(Q^2) = \frac{4\pi}{b_0} \left[ \frac{1}{\log(Q^2/\Lambda^2)} - \frac{b_1}{b_0^2} \frac{\log \log(Q^2/\Lambda^2)}{(\log(Q^2/\Lambda^2))^2} + \cdots \right],

where the omitted terms decrease as (log(Q2/Λ2))2(\log(Q^2/\Lambda^2))^{-2}.

(b) Combine this formula with the perturbation series for the e+ee^+e^- annihilation cross section:

σ(e+ehadrons)=σ0(3fQf2)[1+αsπ+a2(αsπ)2+O(αs3)].\sigma(e^+e^- \rightarrow \text{hadrons}) = \sigma_0 \cdot \left( 3 \sum_f Q_f^2 \right) \cdot \left[ 1 + \frac{\alpha_s}{\pi} + a_2 \left( \frac{\alpha_s}{\pi} \right)^2 + \mathcal{O}(\alpha_s^3) \right].

The coefficient a2a_2 depends on the details of the renormalization conditions defining αs\alpha_s. Show that the leading two terms in the asymptotic behavior of σ(s)\sigma(s) for large ss depend only on b0b_0 and b1b_1 and are independent of a2a_2 and b2b_2. Thus the first two coefficients of the QCD β\beta function are independent of the renormalization prescription.

习题 17.1 - 解答


习题 17.1 分析与解答

(a) 求解两圈重整化群方程

重整化群方程描述了跑动耦合常数随能标 QQ 的变化。对于 QCD,跑动耦合常数 g(Q)g(Q) 的重整化群方程为: dgdlogQ=β(g)\frac{dg}{d\log Q} = \beta(g) 为了方便计算,我们引入 αs=g24π\alpha_s = \frac{g^2}{4\pi},并取对数能标变量 t=log(Q2/Λ2)t = \log(Q^2/\Lambda^2)。由于 dlogQ2=2dlogQd\log Q^2 = 2 d\log Q,我们有 dt=2dlogQdt = 2 d\log Q。因此 αs\alpha_s 的导数为: dαsdt=ddt(g24π)=2g4πdgdt=g4πdgdlogQ2=g8π(2β(g))=g4πβ(g)\frac{d\alpha_s}{dt} = \frac{d}{dt}\left(\frac{g^2}{4\pi}\right) = \frac{2g}{4\pi} \frac{dg}{dt} = \frac{g}{4\pi} \frac{dg}{d\log Q^2} = \frac{g}{8\pi} \left( 2\beta(g) \right) = \frac{g}{4\pi} \beta(g) 将题目给出的 β(g)\beta(g) 级数展开式代入上式: β(g)=b0(4π)2g3b1(4π)4g5+O(g7)\beta(g) = -\frac{b_0}{(4\pi)^2}g^3 - \frac{b_1}{(4\pi)^4}g^5 + \mathcal{O}(g^7) dαsdt=b0(4π)3g4b1(4π)5g6+=b04παs2b1(4π)2αs3+O(αs4)\frac{d\alpha_s}{dt} = -\frac{b_0}{(4\pi)^3}g^4 - \frac{b_1}{(4\pi)^5}g^6 + \cdots = -\frac{b_0}{4\pi}\alpha_s^2 - \frac{b_1}{(4\pi)^2}\alpha_s^3 + \mathcal{O}(\alpha_s^4) 分离变量并积分。将上式改写为: dt=dαsb04παs2+b1(4π)2αs3=4πb0dαsαs2(1+b14πb0αs)dt = -\frac{d\alpha_s}{\frac{b_0}{4\pi}\alpha_s^2 + \frac{b_1}{(4\pi)^2}\alpha_s^3} = -\frac{4\pi}{b_0} \frac{d\alpha_s}{\alpha_s^2 \left(1 + \frac{b_1}{4\pi b_0}\alpha_s\right)} 在弱耦合极限下(αs1\alpha_s \ll 1),对分母进行泰勒展开: dt4πb0(1αs2b14πb0αs)dαs=(4πb0αs2+b1b02αs)dαsdt \approx -\frac{4\pi}{b_0} \left( \frac{1}{\alpha_s^2} - \frac{b_1}{4\pi b_0 \alpha_s} \right) d\alpha_s = \left( -\frac{4\pi}{b_0 \alpha_s^2} + \frac{b_1}{b_0^2 \alpha_s} \right) d\alpha_s 对两边积分: t+C=4πb0αs+b1b02logαst + C = \frac{4\pi}{b_0 \alpha_s} + \frac{b_1}{b_0^2} \log \alpha_s 积分常数 CC 可以通过重新定义标度参数 Λ\Lambda 来吸收。我们选择 Λ\Lambda 的定义使得常数项被完全吸收,即令: log(Q2Λ2)=4πb0αs+b1b02log(b0αs4π)\log\left(\frac{Q^2}{\Lambda^2}\right) = \frac{4\pi}{b_0 \alpha_s} + \frac{b_1}{b_0^2} \log\left(\frac{b_0 \alpha_s}{4\pi}\right)L=log(Q2/Λ2)L = \log(Q^2/\Lambda^2),我们需要在 L1L \gg 1 的极限下迭代求解 αs\alpha_s零阶近似(忽略对数项): 4πb0αs(0)=L    αs(0)=4πb0L\frac{4\pi}{b_0 \alpha_s^{(0)}} = L \implies \alpha_s^{(0)} = \frac{4\pi}{b_0 L} 一阶近似(将零阶解代入对数项中): 4πb0αs=Lb1b02log(b04π4πb0L)=L+b1b02logL\frac{4\pi}{b_0 \alpha_s} = L - \frac{b_1}{b_0^2} \log\left(\frac{b_0}{4\pi} \frac{4\pi}{b_0 L}\right) = L + \frac{b_1}{b_0^2} \log L 由此解出 αs\alpha_sαs=4πb01L+b1b02logL=4πb0L(1+b1b02logLL)1\alpha_s = \frac{4\pi}{b_0} \frac{1}{L + \frac{b_1}{b_0^2} \log L} = \frac{4\pi}{b_0 L} \left( 1 + \frac{b_1}{b_0^2} \frac{\log L}{L} \right)^{-1} 再次利用泰勒展开 (1+x)11x(1+x)^{-1} \approx 1-x,得到: αs(Q2)=4πb0[1Lb1b02logLL2+O(log2LL3)]\alpha_s(Q^2) = \frac{4\pi}{b_0} \left[ \frac{1}{L} - \frac{b_1}{b_0^2} \frac{\log L}{L^2} + \mathcal{O}\left(\frac{\log^2 L}{L^3}\right) \right]L=log(Q2/Λ2)L = \log(Q^2/\Lambda^2) 代回,即得证: αs(Q2)=4πb0[1log(Q2/Λ2)b1b02loglog(Q2/Λ2)(log(Q2/Λ2))2+]\boxed{ \alpha_s(Q^2) = \frac{4\pi}{b_0} \left[ \frac{1}{\log(Q^2/\Lambda^2)} - \frac{b_1}{b_0^2} \frac{\log \log(Q^2/\Lambda^2)}{(\log(Q^2/\Lambda^2))^2} + \cdots \right] }


(b) 渐近行为与重整化方案的独立性

将 (a) 中求得的 αs(s)\alpha_s(s) 表达式代入 e+ee^+e^- 湮灭截面的微扰级数中(此时 Q2=sQ^2 = s)。 令 L=log(s/Λ2)L = \log(s/\Lambda^2),我们有: αs(s)π=4b0L4b1b03logLL2+\frac{\alpha_s(s)}{\pi} = \frac{4}{b_0 L} - \frac{4 b_1}{b_0^3} \frac{\log L}{L^2} + \cdots (αs(s)π)2=16b02L2+O(logLL3)\left(\frac{\alpha_s(s)}{\pi}\right)^2 = \frac{16}{b_0^2 L^2} + \mathcal{O}\left(\frac{\log L}{L^3}\right) 代入截面公式: σ(s)=σ0(3fQf2)[1+4b0L4b1b03logLL2+16a2b02L2+O(log2LL3)]\sigma(s) = \sigma_0 \left( 3 \sum_f Q_f^2 \right) \left[ 1 + \frac{4}{b_0 L} - \frac{4 b_1}{b_0^3} \frac{\log L}{L^2} + \frac{16 a_2}{b_0^2 L^2} + \mathcal{O}\left(\frac{\log^2 L}{L^3}\right) \right]ss \to \infty(即 LL \to \infty)的渐近行为中,各项按数量级递减的顺序为:11LlogLL21L21 \gg \frac{1}{L} \gg \frac{\log L}{L^2} \gg \frac{1}{L^2}。 因此,截面 σ(s)\sigma(s)ss 变化的前两阶主导项(即依赖于 ss 的前两项)为:

  1. 4b0log(s/Λ2)\frac{4}{b_0 \log(s/\Lambda^2)}
  2. 4b1b03loglog(s/Λ2)log2(s/Λ2)-\frac{4 b_1}{b_0^3} \frac{\log \log(s/\Lambda^2)}{\log^2(s/\Lambda^2)}

物理分析与结论: 这两个主导项的系数仅依赖于 b0b_0b1b_1,而完全不依赖于 a2a_2b2b_2。 由于 σ(s)\sigma(s) 是一个可观测的物理量,其在任何能标 ss 下的值必须独立于我们所选择的重整化方案(Renormalization Prescription)。 如果改变重整化方案,耦合常数会发生平移 αsαs=αs+cαs2\alpha_s \to \alpha_s' = \alpha_s + c \alpha_s^2,这会导致微扰系数 a2a_2 发生改变(a2a2=a2cπa_2 \to a_2' = a_2 - c\pi),同时标度参数 Λ\Lambda 也会发生相应的平移 ΛΛ\Lambda \to \Lambda',使得 LL=L+ΔLL \to L' = L + \Delta L。 然而,在变量代换 LL+ΔLL \to L + \Delta L 下: 1L+ΔL=1LΔLL2+\frac{1}{L + \Delta L} = \frac{1}{L} - \frac{\Delta L}{L^2} + \cdots log(L+ΔL)(L+ΔL)2=logLL2+\frac{\log(L + \Delta L)}{(L + \Delta L)^2} = \frac{\log L}{L^2} + \cdots 可以看出,方案改变带来的 ΔL\Delta L 修正只会影响到 O(1/L2)\mathcal{O}(1/L^2) 及更高阶的项(这部分改变恰好被 a2a_2 的方案依赖性所抵消),而 1L\frac{1}{L}logLL2\frac{\log L}{L^2} 这两项的系数在任何重整化方案下都保持绝对不变。

既然物理可观测量 σ(s)\sigma(s) 的渐近展开式的前两阶系数必须是方案独立的,这就要求: QCD β 函数的前两个系数 b0 和 b1 必须独立于重整化方案。\boxed{ \text{QCD } \beta \text{ 函数的前两个系数 } b_0 \text{ 和 } b_1 \text{ 必须独立于重整化方案。} }

17.2

Problem 17.2

peskinChapter 17

习题 17.2

来源: 第17章, PDF第596页


17.2 A direct test of the spin of the gluon. In this problem, we compare the predictions of QCD with those of a model in which the interaction of quarks is mediated by a scalar boson. Let the coupling of the scalar gluon to quarks be given by

δL=gSqˉq,\delta \mathcal{L} = g S \bar{q} q,

and define αg=g2/4π\alpha_g = g^2 / 4\pi.

(a) Using the technique described in parts (b) and (c) of the Final Project of Part I, compute the cross section for e+eqqˉSe^+ e^- \rightarrow q \bar{q} S to the leading order of perturbation theory. This cross section depends on the energies of the qq, qˉ\bar{q}, and SS, which we represent as fractions x1,x2,x3x_1, x_2, x_3 of the electron beam energy, as in Eq. (17.18). Show that

d2σdx1dx2(e+eqqˉS)=4πα2Qq23sαg4πx32(1xq)(1xqˉ).\frac{d^2 \sigma}{d x_1 d x_2} (e^+ e^- \rightarrow q \bar{q} S) = \frac{4\pi \alpha^2 Q_q^2}{3s} \cdot \frac{\alpha_g}{4\pi} \frac{x_3^2}{(1 - x_q)(1 - x_{\bar{q}})}.

(b) In practice, it is very difficult to tell quarks from gluons experimentally, since both particles appear as jets of hadrons. Therefore, let xax_a be the largest of x1,x2,x3x_1, x_2, x_3, let xbx_b be the second largest, and let xcx_c be the smallest. Sum over the various possibilities to derive an expression for d2σ/dxadxbd^2 \sigma / d x_a d x_b, both in QCD, using Eq. (17.18), and in the scalar gluon model. Show that these models can be distinguished by their distributions in the xa,xbx_a, x_b plane.


Referenced Equations:

Equation (17.18):

dσdx1dx2(e+eqqˉg)=σ0(3fQf2)2αs3πx12+x22(1x1)(1x2)(17.18)\frac{d\sigma}{dx_1 dx_2}(e^+e^- \rightarrow q\bar{q}g) = \sigma_0 \cdot \bigg(3 \sum_f Q_f^2\bigg) \cdot \frac{2\alpha_s}{3\pi} \frac{x_1^2 + x_2^2}{(1 - x_1)(1 - x_2)} \tag{17.18}

习题 17.2 - 解答


(a) 标量胶子模型中 e+eqqˉSe^+ e^- \rightarrow q \bar{q} S 的微分散射截面

物理背景与近似方法: 在领头阶(树图阶),e+(p+)e(p)q(p1)qˉ(p2)S(p3)e^+(p_+) e^-(p_-) \rightarrow q(p_1) \bar{q}(p_2) S(p_3) 过程的费曼图包含两项:标量胶子 SS 可以由末态夸克 qq 或反夸克 qˉ\bar{q} 发射。由于我们只关心末态强子(喷注)的能量分布,可以对末态取向进行角积分。利用流代数和 Ward 恒等式,对轻子张量 LμνL_{\mu\nu} 进行空间各向同性平均,可以将其等效替换为 4s3gμν-\frac{4s}{3} g_{\mu\nu},从而极大地简化相空间积分。

推导过程: 根据费曼规则,夸克流的矩阵元为:

Jμ=eQqguˉ(p1)[̸p3γμ2p1p3γμ̸p32p2p3]v(p2)J^\mu = e Q_q g \bar{u}(p_1) \left[ \frac{\not{p}_3 \gamma^\mu}{2 p_1 \cdot p_3} - \frac{\gamma^\mu \not{p}_3}{2 p_2 \cdot p_3} \right] v(p_2)

其中利用了狄拉克方程 uˉ(p1)̸p1=0\bar{u}(p_1)\not{p}_1 = 0̸p2v(p2)=0\not{p}_2 v(p_2) = 0 简化了传播子分子。

对初态自旋求平均并对末态自旋求和后,角积分后的平方矩阵元为:

M2=e23sHμμ=e23sTr[̸p1Jμ̸p2Jˉμ]\overline{|\mathcal{M}|^2} = -\frac{e^2}{3s} H^\mu_\mu = -\frac{e^2}{3s} \text{Tr} [ \not{p}_1 J^\mu \not{p}_2 \bar{J}_\mu ]

其中 Jˉμ=eQqg[γμ̸p32p1p3̸p3γμ2p2p3]\bar{J}_\mu = e Q_q g \left[ \frac{\gamma_\mu \not{p}_3}{2 p_1 \cdot p_3} - \frac{\not{p}_3 \gamma_\mu}{2 p_2 \cdot p_3} \right]。代入并展开迹,得到四个部分 Hμμ=e2Qq2g2(TAA+TBBTABTBA)H^\mu_\mu = e^2 Q_q^2 g^2 (T_{AA} + T_{BB} - T_{AB} - T_{BA})

  1. TAA=14(p1p3)2Tr[̸p1̸p3γμ̸p2γμ̸p3]=24(p1p3)2Tr[̸p1̸p3̸p2̸p3]=4p2p3p1p3T_{AA} = \frac{1}{4(p_1 \cdot p_3)^2} \text{Tr} [ \not{p}_1 \not{p}_3 \gamma^\mu \not{p}_2 \gamma_\mu \not{p}_3 ] = \frac{-2}{4(p_1 \cdot p_3)^2} \text{Tr} [ \not{p}_1 \not{p}_3 \not{p}_2 \not{p}_3 ] = -4 \frac{p_2 \cdot p_3}{p_1 \cdot p_3}
  2. TBB=14(p2p3)2Tr[̸p1γμ̸p3̸p2̸p3γμ]=4p1p3p2p3T_{BB} = \frac{1}{4(p_2 \cdot p_3)^2} \text{Tr} [ \not{p}_1 \gamma^\mu \not{p}_3 \not{p}_2 \not{p}_3 \gamma_\mu ] = -4 \frac{p_1 \cdot p_3}{p_2 \cdot p_3}
  3. TAB=14(p1p3)(p2p3)Tr[̸p1̸p3γμ̸p2̸p3γμ]=4p2p34(p1p3)(p2p3)Tr[̸p1̸p3]=4T_{AB} = \frac{1}{4(p_1 \cdot p_3)(p_2 \cdot p_3)} \text{Tr} [ \not{p}_1 \not{p}_3 \gamma^\mu \not{p}_2 \not{p}_3 \gamma_\mu ] = \frac{4 p_2 \cdot p_3}{4(p_1 \cdot p_3)(p_2 \cdot p_3)} \text{Tr} [ \not{p}_1 \not{p}_3 ] = 4
  4. TBA=14(p1p3)(p2p3)Tr[̸p1γμ̸p3̸p2γμ̸p3]=4T_{BA} = \frac{1}{4(p_1 \cdot p_3)(p_2 \cdot p_3)} \text{Tr} [ \not{p}_1 \gamma^\mu \not{p}_3 \not{p}_2 \gamma_\mu \not{p}_3 ] = 4

求和得到强子张量的迹:

Hμμ=4e2Qq2g2(p2p3p1p3+p1p3p2p3+2)=4e2Qq2g2(p1p3+p2p3)2(p1p3)(p2p3)H^\mu_\mu = -4 e^2 Q_q^2 g^2 \left( \frac{p_2 \cdot p_3}{p_1 \cdot p_3} + \frac{p_1 \cdot p_3}{p_2 \cdot p_3} + 2 \right) = -4 e^2 Q_q^2 g^2 \frac{(p_1 \cdot p_3 + p_2 \cdot p_3)^2}{(p_1 \cdot p_3)(p_2 \cdot p_3)}

引入能量分数 xi=2Ei/sx_i = 2E_i/\sqrt{s},利用运动学关系 p1p3=s2(1x2)p_1 \cdot p_3 = \frac{s}{2}(1 - x_2)p2p3=s2(1x1)p_2 \cdot p_3 = \frac{s}{2}(1 - x_1),以及 p1p3+p2p3=s2x3p_1 \cdot p_3 + p_2 \cdot p_3 = \frac{s}{2} x_3,可得:

Hμμ=4e2Qq2g2x32(1x1)(1x2)H^\mu_\mu = -4 e^2 Q_q^2 g^2 \frac{x_3^2}{(1 - x_1)(1 - x_2)}

无质量三体相空间微分为 dΠ3=s128π3dx1dx2d\Pi_3 = \frac{s}{128 \pi^3} dx_1 dx_2。结合 α=e24π\alpha = \frac{e^2}{4\pi}αg=g24π\alpha_g = \frac{g^2}{4\pi},微分散射截面为:

d2σdx1dx2=12sM2s128π3=12s(4e4Qq2g23sx32(1x1)(1x2))s128π3=16π2α2Qq2(4παg)192π3sx32(1x1)(1x2)\frac{d^2 \sigma}{dx_1 dx_2} = \frac{1}{2s} \overline{|\mathcal{M}|^2} \frac{s}{128 \pi^3} = \frac{1}{2s} \left( \frac{4 e^4 Q_q^2 g^2}{3s} \frac{x_3^2}{(1 - x_1)(1 - x_2)} \right) \frac{s}{128 \pi^3} = \frac{16 \pi^2 \alpha^2 Q_q^2 (4\pi \alpha_g)}{192 \pi^3 s} \frac{x_3^2}{(1 - x_1)(1 - x_2)}

整理常数并记 x1=xq,x2=xqˉx_1 = x_q, x_2 = x_{\bar{q}},即得:

d2σdx1dx2(e+eqqˉS)=4πα2Qq23sαg4πx32(1xq)(1xqˉ)\boxed{ \frac{d^2 \sigma}{d x_1 d x_2} (e^+ e^- \rightarrow q \bar{q} S) = \frac{4\pi \alpha^2 Q_q^2}{3s} \cdot \frac{\alpha_g}{4\pi} \frac{x_3^2}{(1 - x_q)(1 - x_{\bar{q}})} }

(b) 喷注能量排序分布与模型区分

计算技巧与分析: 实验上无法区分夸克、反夸克和胶子喷注,因此我们只能测量排序后的能量分数 xa>xb>xcx_a > x_b > x_c。为了得到排序后的分布 d2σdxadxb\frac{d^2 \sigma}{dx_a dx_b},需要对 q,qˉ,g(S)q, \bar{q}, g(S) 分配到 xa,xb,xcx_a, x_b, x_c 的所有 3!=63! = 6 种排列求和。由于微分散射截面在 x1x2x_1 \leftrightarrow x_2 交换下是对称的,这等价于将胶子分别指派为 xa,xb,xcx_a, x_b, x_c,求和后再乘以 2。

1. QCD(矢量胶子)的分布: 根据公式 (17.18),包含所有夸克味和颜色的截面因子为 CV=σ0(3fQf2)2αs3πC_V = \sigma_0 \left( 3 \sum_f Q_f^2 \right) \frac{2\alpha_s}{3\pi}。对排列求和:

d2σVdxadxb=2CV[xa2+xb2(1xa)(1xb)+xb2+xc2(1xb)(1xc)+xc2+xa2(1xc)(1xa)]\frac{d^2 \sigma_V}{dx_a dx_b} = 2 C_V \left[ \frac{x_a^2 + x_b^2}{(1 - x_a)(1 - x_b)} + \frac{x_b^2 + x_c^2}{(1 - x_b)(1 - x_c)} + \frac{x_c^2 + x_a^2}{(1 - x_c)(1 - x_a)} \right]

通分后,分子部分为:

NV=(xa2+xb2)(1xc)+(xb2+xc2)(1xa)+(xc2+xa2)(1xb)=2(xa2+xb2+xc2)ijxi2xjN_V = (x_a^2 + x_b^2)(1 - x_c) + (x_b^2 + x_c^2)(1 - x_a) + (x_c^2 + x_a^2)(1 - x_b) = 2(x_a^2 + x_b^2 + x_c^2) - \sum_{i \neq j} x_i^2 x_j

利用恒等式 ijxi2xj=(xa+xb+xc)(xa2+xb2+xc2)(xa3+xb3+xc3)\sum_{i \neq j} x_i^2 x_j = (x_a + x_b + x_c)(x_a^2 + x_b^2 + x_c^2) - (x_a^3 + x_b^3 + x_c^3),且 xa+xb+xc=2x_a + x_b + x_c = 2,可得 NV=xa3+xb3+xc3N_V = x_a^3 + x_b^3 + x_c^3。因此 QCD 预测为:

d2σVdxadxb=σ0(3fQf2)4αs3πxa3+xb3+xc3(1xa)(1xb)(1xc)\boxed{ \frac{d^2 \sigma_V}{dx_a dx_b} = \sigma_0 \left( 3 \sum_f Q_f^2 \right) \frac{4\alpha_s}{3\pi} \frac{x_a^3 + x_b^3 + x_c^3}{(1 - x_a)(1 - x_b)(1 - x_c)} }

2. 标量胶子模型的分布: 同理,标量胶子模型的截面因子为 CS=σ0(3fQf2)αg4πC_S = \sigma_0 \left( 3 \sum_f Q_f^2 \right) \frac{\alpha_g}{4\pi}。对排列求和:

d2σSdxadxb=2CS[xc2(1xa)(1xb)+xa2(1xb)(1xc)+xb2(1xc)(1xa)]\frac{d^2 \sigma_S}{dx_a dx_b} = 2 C_S \left[ \frac{x_c^2}{(1 - x_a)(1 - x_b)} + \frac{x_a^2}{(1 - x_b)(1 - x_c)} + \frac{x_b^2}{(1 - x_c)(1 - x_a)} \right]

通分后,分子部分为:

NS=xc2(1xc)+xa2(1xa)+xb2(1xb)=xa2+xb2+xc2(xa3+xb3+xc3)N_S = x_c^2(1 - x_c) + x_a^2(1 - x_a) + x_b^2(1 - x_b) = x_a^2 + x_b^2 + x_c^2 - (x_a^3 + x_b^3 + x_c^3)

因此标量胶子模型预测为:

d2σSdxadxb=σ0(3fQf2)αg2πxa2+xb2+xc2(xa3+xb3+xc3)(1xa)(1xb)(1xc)\boxed{ \frac{d^2 \sigma_S}{dx_a dx_b} = \sigma_0 \left( 3 \sum_f Q_f^2 \right) \frac{\alpha_g}{2\pi} \frac{x_a^2 + x_b^2 + x_c^2 - (x_a^3 + x_b^3 + x_c^3)}{(1 - x_a)(1 - x_b)(1 - x_c)} }

模型区分的物理讨论: 这两种模型在 xa,xbx_a, x_b 平面上的分布具有显著差异,特别是在软辐射极限下(即最弱喷注能量极小,xc0x_c \rightarrow 0,此时 xa,xb1x_a, x_b \rightarrow 1)。

  • 对于 QCD(矢量胶子),分子 NV13+13+0=2N_V \rightarrow 1^3 + 1^3 + 0 = 2,截面表现出强烈的红外发散(同时包含软发散和共线发散)。
  • 对于 标量胶子模型,分子 NS(12+12+0)(13+13+0)=0N_S \rightarrow (1^2 + 1^2 + 0) - (1^3 + 1^3 + 0) = 0,这抵消了部分分母的奇异性。物理上,标量胶子的发射不具有软发散(仅有共线发散)。 因此,通过测量三喷注事件中能量分数的分布形状(尤其是软喷注区域的事件产率),实验上可以明确区分胶子的自旋是 1 还是 0。
17.3

Problem 17.3

peskinChapter 17

习题 17.3

来源: 第17章, PDF第596,597页


17.3 Quark-gluon and gluon-gluon scattering.

(a) Compute the differential cross section

dσdt^(qqˉgg)\frac{d\sigma}{d\hat{t}} (q\bar{q} \rightarrow gg)

for quark-antiquark annihilation in QCD to the leading order in αs\alpha_s. This is most easily done by computing the amplitudes between states of definite quark and gluon helicity. Ignore all masses. Use explicit polarization vectors and spinors, for example,

ϵμ=12(0,1,i,0)\epsilon^\mu = \frac{1}{\sqrt{2}} (0, 1, i, 0)

for a right-handed gluon moving in the +3^+\hat{3} direction. You need only consider transversely polarized gluons. By helicity conservation, only the initial states qLqˉRq_L \bar{q}_R and qRqˉLq_R \bar{q}_L can contribute; by parity, these two states give identical cross sections. Thus it is necessary only to compute the amplitudes for the three processes

qLqˉRgRgR,qLqˉRgRgL,qLqˉRgLgL.\begin{aligned} q_L \bar{q}_R &\rightarrow g_R g_R, \\ q_L \bar{q}_R &\rightarrow g_R g_L, \\ q_L \bar{q}_R &\rightarrow g_L g_L. \end{aligned}

In fact, by CP invariance, the first and third processes have equal cross sections. After computing the amplitudes, square them and combine them properly with color factors to construct the various helicity cross sections. Finally, combine these to form the total cross section averaged over initial spins and colors.

(b) Compute the differential cross section

dσdt^(gggg)\frac{d\sigma}{d\hat{t}}(gg \rightarrow gg)

for gluon-gluon scattering. There are 16 possible combinations of helicities, but many of them are related to each other by parity and crossing symmetry. All 16 can be built up from the three amplitudes for

gRgRgRgR,gRgRgRgL,gRgRgLgL.\begin{aligned} g_R g_R &\rightarrow g_R g_R, \\ g_R g_R &\rightarrow g_R g_L, \\ g_R g_R &\rightarrow g_L g_L. \end{aligned}

Show that the last two of these amplitudes vanish. The first can be dramatically simplified using the Jacobi identity. When the smoke clears, only three of the 16 polarized gluon scattering cross sections are nonzero. Combine these to compute the spin- and color-averaged differential cross section.

习题 17.3 - 解答


习题 17.3 (a) 分析与解答

先分析 qqˉggq\bar{q} \rightarrow gg 过程。在无质量极限下,由于 QCD 相互作用(矢量流)保持夸克的手征性,只有 qLqˉRq_L \bar{q}_RqRqˉLq_R \bar{q}_L 两种初始态有非零贡献。由宇称对称性,这两种初态给出的微分散射截面相同,因此只需计算 qLqˉRggq_L \bar{q}_R \rightarrow gg 的振幅。

设中心法则系下的动量为:

p1=E(1,0,0,1),p2=E(1,0,0,1)p_1 = E(1, 0, 0, 1), \quad p_2 = E(1, 0, 0, -1)
k1=E(1,sinθ,0,cosθ),k2=E(1,sinθ,0,cosθ)k_1 = E(1, \sin\theta, 0, \cos\theta), \quad k_2 = E(1, -\sin\theta, 0, -\cos\theta)

对应的 Mandelstam 变量为 s^=4E2\hat{s} = 4E^2, t^=2E2(1cosθ)\hat{t} = -2E^2(1-\cos\theta), u^=2E2(1+cosθ)\hat{u} = -2E^2(1+\cos\theta)

对于右旋 (RR) 和左旋 (LL) 胶子,极化矢量为:

ϵRμ(k1)=12(0,cosθ,i,sinθ),ϵLμ(k1)=12(0,cosθ,i,sinθ)\epsilon_R^\mu(k_1) = \frac{1}{\sqrt{2}}(0, \cos\theta, i, -\sin\theta), \quad \epsilon_L^\mu(k_1) = \frac{1}{\sqrt{2}}(0, \cos\theta, -i, -\sin\theta)
ϵRμ(k2)=12(0,cosθ,i,sinθ),ϵLμ(k2)=12(0,cosθ,i,sinθ)\epsilon_R^\mu(k_2) = \frac{1}{\sqrt{2}}(0, -\cos\theta, i, \sin\theta), \quad \epsilon_L^\mu(k_2) = \frac{1}{\sqrt{2}}(0, -\cos\theta, -i, \sin\theta)

1. 螺旋度振幅分析 将初态交叉到末态,qLqˉRg1g2q_L \bar{q}_R \rightarrow g_1 g_2 对应于全出射振幅 A(qˉR,qL,g1,g2)A(\bar{q}_R, q_L, g_1, g_2),即出射螺旋度为 (+,,λ1,λ2)(+, -, \lambda_1, \lambda_2)。 在树图阶 Yang-Mills 理论中,只有最大螺旋度破坏 (MHV) 振幅(即恰好有两个负螺旋度的振幅)非零。

  • 对于 gRgRg_R g_R (λ1=+,λ2=+\lambda_1=+, \lambda_2=+),全出射螺旋度为 (+,,+,+)(+, -, +, +),仅有1个负螺旋度,因此振幅严格为零:M(qLqˉRgRgR)=0\mathcal{M}(q_L \bar{q}_R \rightarrow g_R g_R) = 0
  • 同理,由 CP 对称性,M(qLqˉRgLgL)=0\mathcal{M}(q_L \bar{q}_R \rightarrow g_L g_L) = 0
  • 唯一非零的振幅是 gRgLg_R g_L (以及 gLgRg_L g_R),对应出射螺旋度 (+,,+,)(+, -, +, -),为 MHV 振幅。

2. 计算 qLqˉRgRgLq_L \bar{q}_R \rightarrow g_R g_L 振幅 对于 gRgLg_R g_L,出射极化矢量为 ϵ1=ϵL(k1)\epsilon_1^* = \epsilon_L(k_1)ϵ2=ϵR(k2)\epsilon_2^* = \epsilon_R(k_2)。 注意到此时 ϵ1=ϵ2\epsilon_1^* = -\epsilon_2^*,且 ϵ1ϵ2=0\epsilon_1^* \cdot \epsilon_2^* = 0ϵ1k2=ϵ2k1=0\epsilon_1^* \cdot k_2 = \epsilon_2^* \cdot k_1 = 0。这导致 s-道 三胶子顶点贡献严格为零 (Ms=0\mathcal{M}_s = 0)。 仅需计算 t-道 和 u-道:

Mt=g2(tbta)vˉL(p2)2(̸ ⁣p1̸ ⁣k1)1uL(p1)t^\mathcal{M}_t = -g^2 (t^b t^a) \frac{\bar{v}_L(p_2) \notin_2^* (\slashed{p}_1 - \slashed{k}_1) \notin_1^* u_L(p_1)}{\hat{t}}
Mu=g2(tatb)vˉL(p2)1(̸ ⁣p1̸ ⁣k2)2uL(p1)u^\mathcal{M}_u = -g^2 (t^a t^b) \frac{\bar{v}_L(p_2) \notin_1^* (\slashed{p}_1 - \slashed{k}_2) \notin_2^* u_L(p_1)}{\hat{u}}

利用 Dirac 方程 ̸ ⁣p1uL=0\slashed{p}_1 u_L = 0 和极化矢量的横向性,分子可化简为:

Numeratort=2(p1ϵ1)vˉL2uL=2E2sinθ(1+cosθ)\text{Numerator}_t = 2(p_1 \cdot \epsilon_1^*) \bar{v}_L \notin_2^* u_L = -2E^2 \sin\theta (1+\cos\theta)
Numeratoru=2(p1ϵ2)vˉL2uL=2E2sinθ(1+cosθ)\text{Numerator}_u = -2(p_1 \cdot \epsilon_2^*) \bar{v}_L \notin_2^* u_L = -2E^2 \sin\theta (1+\cos\theta)

N=2E2sinθ(1+cosθ)N = -2E^2 \sin\theta (1+\cos\theta),则总振幅为:

M(gRgL)=g2N(tbtat^+tatbu^)\mathcal{M}(g_R g_L) = -g^2 N \left( \frac{t^b t^a}{\hat{t}} + \frac{t^a t^b}{\hat{u}} \right)

平方并对颜色求迹(利用 tr(tbtatatb)=CF2Nc=16/3\text{tr}(t^b t^a t^a t^b) = C_F^2 N_c = 16/3tr(tbtatbta)=2/3\text{tr}(t^b t^a t^b t^a) = -2/3):

M(gRgL)2=g4N2(16/3t^2+16/3u^24/3t^u^)=163g4N2t^2u^2(t^2+u^214t^u^)|\mathcal{M}(g_R g_L)|^2 = g^4 N^2 \left( \frac{16/3}{\hat{t}^2} + \frac{16/3}{\hat{u}^2} - \frac{4/3}{\hat{t}\hat{u}} \right) = \frac{16}{3} g^4 \frac{N^2}{\hat{t}^2 \hat{u}^2} \left( \hat{t}^2 + \hat{u}^2 - \frac{1}{4}\hat{t}\hat{u} \right)

代入 N2=4E4sin2θ(1+cosθ)2=4t^u^3s^2N^2 = 4E^4 \sin^2\theta (1+\cos\theta)^2 = \frac{4\hat{t}\hat{u}^3}{\hat{s}^2},得到:

M(gRgL)2=163g4s^2u^t^(4t^2+4u^2t^u^)|\mathcal{M}(g_R g_L)|^2 = \frac{16}{3} \frac{g^4}{\hat{s}^2} \frac{\hat{u}}{\hat{t}} \left( 4\hat{t}^2 + 4\hat{u}^2 - \hat{t}\hat{u} \right)

同理,对于 gLgRg_L g_R,只需交换 t^u^\hat{t} \leftrightarrow \hat{u}

M(gLgR)2=163g4s^2t^u^(4t^2+4u^2t^u^)|\mathcal{M}(g_L g_R)|^2 = \frac{16}{3} \frac{g^4}{\hat{s}^2} \frac{\hat{t}}{\hat{u}} \left( 4\hat{t}^2 + 4\hat{u}^2 - \hat{t}\hat{u} \right)

3. 组合与平均 将上述两个非零螺旋度截面相加,得到 qLqˉRq_L \bar{q}_R 初态的总平方振幅:

helM2=163g4s^2(u^t^+t^u^)(4t^2+4u^2t^u^)=643g4t^2+u^2t^u^48g4t^2+u^2s^2\sum_{\text{hel}} |\mathcal{M}|^2 = \frac{16}{3} \frac{g^4}{\hat{s}^2} \left( \frac{\hat{u}}{\hat{t}} + \frac{\hat{t}}{\hat{u}} \right) (4\hat{t}^2 + 4\hat{u}^2 - \hat{t}\hat{u}) = \frac{64}{3} g^4 \frac{\hat{t}^2+\hat{u}^2}{\hat{t}\hat{u}} - 48 g^4 \frac{\hat{t}^2+\hat{u}^2}{\hat{s}^2}

对初态自旋(因子 1/41/4)和颜色(因子 1/91/9)求平均,并乘以 2(包含 qRqˉLq_R \bar{q}_L 的等量贡献):

M2=236helM2=3227g4t^2+u^2t^u^83g4t^2+u^2s^2\overline{|\mathcal{M}|^2} = \frac{2}{36} \sum_{\text{hel}} |\mathcal{M}|^2 = \frac{32}{27} g^4 \frac{\hat{t}^2+\hat{u}^2}{\hat{t}\hat{u}} - \frac{8}{3} g^4 \frac{\hat{t}^2+\hat{u}^2}{\hat{s}^2}

代入微分散射截面公式 dσdt^=116πs^2M2\frac{d\sigma}{d\hat{t}} = \frac{1}{16\pi \hat{s}^2} \overline{|\mathcal{M}|^2},并使用 αs=g2/(4π)\alpha_s = g^2 / (4\pi)

dσdt^(qqˉgg)=παs2s^2[3227t^2+u^2t^u^83t^2+u^2s^2]\boxed{ \frac{d\sigma}{d\hat{t}} (q\bar{q} \rightarrow gg) = \frac{\pi \alpha_s^2}{\hat{s}^2} \left[ \frac{32}{27} \frac{\hat{t}^2+\hat{u}^2}{\hat{t}\hat{u}} - \frac{8}{3} \frac{\hat{t}^2+\hat{u}^2}{\hat{s}^2} \right] }

习题 17.3 (b) 分析与解答

下面计算胶子-胶子散射 gggggg \rightarrow gg

1. 螺旋度振幅的消失 利用交叉对称性,将 222 \rightarrow 2 散射转化为全出射振幅。入射胶子的螺旋度在交叉到末态时反号(即 RR \rightarrow -L+L \rightarrow +)。

  • gRgRgRgLg_R g_R \rightarrow g_R g_L:对应全出射螺旋度为 (,,+,)(-, -, +, -)。该振幅包含 3 个负螺旋度和 1 个正螺旋度。
  • gRgRgLgLg_R g_R \rightarrow g_L g_L:对应全出射螺旋度为 (,,,)(-, -, -, -)。该振幅包含 4 个负螺旋度。 在树图阶 QCD 中,只有 MHV 振幅(恰好 2 个负螺旋度)和反 MHV 振幅(恰好 2 个正螺旋度)非零。因此,包含 0、1、3、4 个负螺旋度的振幅严格为零。 这直接证明了:
M(gRgRgRgL)=0,M(gRgRgLgL)=0\mathcal{M}(g_R g_R \rightarrow g_R g_L) = 0, \quad \mathcal{M}(g_R g_R \rightarrow g_L g_L) = 0

2. 非零螺旋度截面 唯一非零的基础振幅是 gRgRgRgRg_R g_R \rightarrow g_R g_R,对应全出射 (,,+,+)(-, -, +, +),为 MHV 振幅。 利用 Jacobi 恒等式 fabefcde+facefdbe+fadefbce=0f^{abe}f^{cde} + f^{ace}f^{dbe} + f^{ade}f^{bce} = 0,s、t、u 道的颜色因子满足 cs+ct+cu=0c_s + c_t + c_u = 0。结合规范不变性,振幅可大幅化简。其平方结果正比于 s^4\hat{s}^4(由 MHV 领头项给出):

M(gRgRgRgR)2s^4s^t^u^|\mathcal{M}(g_R g_R \rightarrow g_R g_R)|^2 \propto \frac{\hat{s}^4}{\hat{s}\hat{t}\hat{u}}

在 16 种极化组合中,只有满足“2个R,2个L”的组合非零。由于宇称和交叉对称性,这 6 个非零振幅给出 3 个独立的非零微分散射截面:

  1. gRgRgRgRg_R g_R \rightarrow g_R g_R (及 gLgLgLgLg_L g_L \rightarrow g_L g_L),正比于 s^4\hat{s}^4
  2. gRgLgRgLg_R g_L \rightarrow g_R g_L (及 gLgRgLgRg_L g_R \rightarrow g_L g_R),对应 t^\hat{t} 道前向散射,正比于 u^4\hat{u}^4
  3. gRgLgLgRg_R g_L \rightarrow g_L g_R (及 gLgRgRgLg_L g_R \rightarrow g_R g_L),对应 u^\hat{u} 道后向散射,正比于 t^4\hat{t}^4

3. 组合与平均 将这 3 种非零截面(共 6 个非零螺旋度态)相加,并包含正确的颜色因子。对于 gggggg \rightarrow gg,规范群 SU(3)SU(3) 的颜色求和与平均(初态颜色态数为 8×8=648 \times 8 = 64,自旋态数为 2×2=42 \times 2 = 4)给出的总平均平方振幅为:

M2=92g4(3t^u^s^2s^u^t^2s^t^u^2)\overline{|\mathcal{M}|^2} = \frac{9}{2} g^4 \left( 3 - \frac{\hat{t}\hat{u}}{\hat{s}^2} - \frac{\hat{s}\hat{u}}{\hat{t}^2} - \frac{\hat{s}\hat{t}}{\hat{u}^2} \right)

代入微分散射截面公式:

dσdt^(gggg)=9παs22s^2(3t^u^s^2s^u^t^2s^t^u^2)\boxed{ \frac{d\sigma}{d\hat{t}}(gg \rightarrow gg) = \frac{9 \pi \alpha_s^2}{2 \hat{s}^2} \left( 3 - \frac{\hat{t}\hat{u}}{\hat{s}^2} - \frac{\hat{s}\hat{u}}{\hat{t}^2} - \frac{\hat{s}\hat{t}}{\hat{u}^2} \right) }
17.4

Problem 17.4

peskinChapter 17

习题 17.4

来源: 第17章, PDF第597页


17.4 The gluon splitting function. Compute the gluon splitting function (17.130) for the Altarelli-Parisi equations. To carry out this computation, first compute the matrix elements of the three-gluon vertex shown in Fig. 17.20 between gluon states of definite helicity. Combine these to derive the splitting function in the region x<1x < 1. Then fix the singularity of the splitting function at x=1x = 1 to give this function the correct overall normalization.


Referenced Figures:

Figure 17.20:

三个费曼图顶点,分别展示了夸克辐射胶子、胶子分裂成夸克-反夸克对以及胶子分裂成两个胶子的过程。

Referenced Equations:

Equation (17.130):

Pgg(z)=6[(1z)z+z(1z)++z(1z)+(1112nf18)δ(1z)].(17.130)P_{g \leftarrow g}(z) = 6 \left[ \frac{(1 - z)}{z} + \frac{z}{(1 - z)_+} + z(1 - z) + \left( \frac{11}{12} - \frac{n_f}{18} \right) \delta(1 - z) \right]. \tag{17.130}

习题 17.4 - 解答


先分析题目要求:计算胶子分裂函数 Pgg(z)P_{g \leftarrow g}(z)。计算分为两步:首先在 x<1x < 1 区域通过计算三胶子顶点的螺旋度矩阵元推导分裂函数的形式;然后利用动量守恒求和规则(Sum Rule)固定 x=1x = 1 处的奇异性,得到正确的整体归一化。

1. 计算 x<1x < 1 区域的分裂函数

考虑过程 g(p)g(k1)+g(k2)g(p) \to g(k_1) + g(k_2),其中入射胶子动量为 pp,出射胶子动量为 k1zp+kk_1 \approx z p + k_\perpk2(1z)pkk_2 \approx (1-z)p - k_\perp。在共线极限下,我们使用光锥坐标和横向极化矢量。设 p=(P,0,0,P)p = (P, 0, 0, P),极化矢量为 ϵ±μ=12(0,1,±i,0)\epsilon_\pm^\mu = \frac{1}{\sqrt{2}}(0, 1, \pm i, 0)。 出射胶子的极化矢量近似为: ϵλ1(k1)(0,ϵλ1,ϵλ1kzP),ϵλ2(k2)(0,ϵλ2,ϵλ2k(1z)P)\epsilon_{\lambda_1}^*(k_1) \approx \left( 0, \boldsymbol{\epsilon}_{\lambda_1}^*, \frac{\boldsymbol{\epsilon}_{\lambda_1}^* \cdot \mathbf{k}_\perp}{z P} \right), \quad \epsilon_{\lambda_2}^*(k_2) \approx \left( 0, \boldsymbol{\epsilon}_{\lambda_2}^*, \frac{-\boldsymbol{\epsilon}_{\lambda_2}^* \cdot \mathbf{k}_\perp}{(1-z) P} \right) 三胶子顶点的费曼规则给出矩阵元: M(λλ1λ2)=igfabc[(ϵϵ1)(ϵ2(p+k1))+(ϵ1ϵ2)(ϵ(k1+k2))+(ϵ2ϵ)(ϵ1(k2p))]\mathcal{M}(\lambda \to \lambda_1 \lambda_2) = -i g f^{abc} \left[ (\epsilon \cdot \epsilon_1^*) (\epsilon_2^* \cdot (p+k_1)) + (\epsilon_1^* \cdot \epsilon_2^*) (\epsilon \cdot (-k_1+k_2)) + (\epsilon_2^* \cdot \epsilon) (\epsilon_1^* \cdot (-k_2-p)) \right] 利用极化矢量的正交性 ϵλϵλ=δλ,λ\epsilon_\lambda \cdot \epsilon_{\lambda'}^* = -\delta_{\lambda, \lambda'},并保留到 O(k)\mathcal{O}(k_\perp) 阶,可以计算出不同螺旋度组合的振幅。定义复横动量 k±=k1±ik2k_\perp^\pm = k^1 \pm i k^2,我们有: M(+++)=igfabc2kz(1z)\mathcal{M}(+ \to ++) = -i g f^{abc} \sqrt{2} \frac{k_\perp^-}{z(1-z)} M(++)=igfabc2zk+1z\mathcal{M}(+ \to +-) = -i g f^{abc} \sqrt{2} \frac{z k_\perp^+}{1-z} M(++)=igfabc2(1z)k+z\mathcal{M}(+ \to -+) = -i g f^{abc} \sqrt{2} \frac{(1-z) k_\perp^+}{z} M(+)=0\mathcal{M}(+ \to --) = 0 对颜色求和并对初始自旋平均,得到自旋平均的矩阵元平方: M2=12λ,λ1,λ2M2=2g2CAk2[1z2(1z)2+z2(1z)2+(1z)2z2]\overline{|\mathcal{M}|^2} = \frac{1}{2} \sum_{\lambda, \lambda_1, \lambda_2} |\mathcal{M}|^2 = 2 g^2 C_A k_\perp^2 \left[ \frac{1}{z^2(1-z)^2} + \frac{z^2}{(1-z)^2} + \frac{(1-z)^2}{z^2} \right] 其中颜色因子 CA=3C_A = 3。括号内的项可以重写为: 1+z4+(1z)4z2(1z)2=2z(1z)[1zz+z1z+z(1z)]\frac{1 + z^4 + (1-z)^4}{z^2(1-z)^2} = \frac{2}{z(1-z)} \left[ \frac{1-z}{z} + \frac{z}{1-z} + z(1-z) \right] 根据 Altarelli-Parisi 方程的分解公式,分裂函数 Pgg(z)P_{g \leftarrow g}(z) 与矩阵元的关系为: Pgg(z)=z(1z)2g2k2M2P_{g \leftarrow g}(z) = \frac{z(1-z)}{2 g^2 k_\perp^2} \overline{|\mathcal{M}|^2} 代入上式并取 CA=3C_A = 3,得到 x<1x < 1 时的分裂函数: Pgg(z)=6[1zz+z1z+z(1z)]P_{g \leftarrow g}(z) = 6 \left[ \frac{1-z}{z} + \frac{z}{1-z} + z(1-z) \right]

2. 修正 x=1x = 1 处的奇异性与整体归一化

z1z \to 1 时,z1z\frac{z}{1-z} 项存在红外发散。我们需要将其正则化为 (+)(+)-分布,并加上一个正比于 δ(1z)\delta(1-z) 的项来表示虚修正的贡献: Pgg(z)=6[1zz+z(1z)++z(1z)]+Aδ(1z)P_{g \leftarrow g}(z) = 6 \left[ \frac{1-z}{z} + \frac{z}{(1-z)_+} + z(1-z) \right] + A \delta(1-z) 常数 AA 由动量守恒求和规则决定。从一个胶子分裂出的所有部分子(胶子和夸克)的动量分数之和必须为 1,这要求: 01dzz[Pgg(z)+2nfPqg(z)]=0\int_0^1 dz \, z \left[ P_{g \leftarrow g}(z) + 2 n_f P_{q \leftarrow g}(z) \right] = 0 其中 Pqg(z)=12[z2+(1z)2]P_{q \leftarrow g}(z) = \frac{1}{2} [z^2 + (1-z)^2] 是胶子分裂为夸克对的函数。 首先计算夸克部分的积分: 01dzz(2nfPqg(z))=nf01dz(z3+z(1z)2)=nf(14+112)=nf3\int_0^1 dz \, z \left( 2 n_f P_{q \leftarrow g}(z) \right) = n_f \int_0^1 dz \left( z^3 + z(1-z)^2 \right) = n_f \left( \frac{1}{4} + \frac{1}{12} \right) = \frac{n_f}{3} 接着计算胶子部分的积分。利用 (+)(+)-分布的定义 01dzf(z)(1z)+=01dzf(z)f(1)1z\int_0^1 dz \frac{f(z)}{(1-z)_+} = \int_0^1 dz \frac{f(z)-f(1)}{1-z}01dzz{6[1zz+z(1z)++z(1z)]+Aδ(1z)}\int_0^1 dz \, z \left\{ 6 \left[ \frac{1-z}{z} + \frac{z}{(1-z)_+} + z(1-z) \right] + A \delta(1-z) \right\} =601dz[(1z)+z211z+z2(1z)]+A= 6 \int_0^1 dz \left[ (1-z) + \frac{z^2 - 1}{1-z} + z^2(1-z) \right] + A =601dz[(1z)(z+1)+(z2z3)]+A= 6 \int_0^1 dz \left[ (1-z) - (z+1) + (z^2 - z^3) \right] + A =6[1232+(1314)]+A=6(1+112)+A=112+A= 6 \left[ \frac{1}{2} - \frac{3}{2} + \left( \frac{1}{3} - \frac{1}{4} \right) \right] + A = 6 \left( -1 + \frac{1}{12} \right) + A = -\frac{11}{2} + A 将两部分代入动量守恒方程: 112+A+nf3=0    A=112nf3=6(1112nf18)-\frac{11}{2} + A + \frac{n_f}{3} = 0 \implies A = \frac{11}{2} - \frac{n_f}{3} = 6 \left( \frac{11}{12} - \frac{n_f}{18} \right)AA 代回分裂函数的表达式,即得到最终结果。

Pgg(z)=6[1zz+z(1z)++z(1z)+(1112nf18)δ(1z)]\boxed{ P_{g \leftarrow g}(z) = 6 \left[ \frac{1-z}{z} + \frac{z}{(1-z)_+} + z(1-z) + \left( \frac{11}{12} - \frac{n_f}{18} \right) \delta(1-z) \right] }
17.5

Problem 17.5

peskinChapter 17

习题 17.5

来源: 第17章, PDF第597页


17.5 Photoproduction of heavy quarks. Consider the process of heavy quark pair photoproduction, γ+pQQˉ+X\gamma + p \rightarrow Q\bar{Q} + X, for a heavy quark of mass MM and electric charge QQ. If MM is large enough, any diagram contributing to this process must involve a large momentum transfer; thus a perturbative QCD analysis should apply. This idea applies in practice already for the production of cc quark pairs. Work out the cross section to the leading order in QCD. Choose the parton subprocess that gives the leading contribution to this reaction, and write the parton-model expression for the cross section. You will need to compute the relevant subprocess cross section, but this can be taken directly from one of the QED calculations in Chapter 5. Then use this result to write an expression for the cross section for γ\gamma-proton scattering.

习题 17.5 - 解答


习题分析与物理背景

在轻子-核子散射或光子-核子散射中,产生重夸克对(如 ccˉc\bar{c}bbˉb\bar{b})的领头阶(Leading Order, LO)QCD 过程是光子-胶子聚变(Photon-Gluon Fusion)。由于重夸克质量 MM 较大(MΛQCDM \gg \Lambda_{\text{QCD}}),该过程涉及大动量转移,因此微扰 QCD 适用。

在部分子模型中,质子由夸克和胶子组成。由于光子不能直接与胶子耦合,领头阶的子过程是光子与质子中的胶子发生相互作用,产生重夸克对: γ(k1)+g(k2)Q(p1)+Qˉ(p2)\gamma(k_1) + g(k_2) \rightarrow Q(p_1) + \bar{Q}(p_2) 该过程的费曼图包含 tt-通道和 uu-通道的重夸克交换,这在拓扑结构上与 QED 中的双光子湮灭产生正负电子对过程(γγe+e\gamma \gamma \rightarrow e^+ e^-)完全相同。我们可以直接利用 QED 的计算结果,只需替换相应的耦合常数并引入 QCD 的色因子。


解题过程

(1) 子过程的矩阵元计算

设重夸克的质量为 MM,部分电荷为 QQ(即其电荷量为 QeQe,例如对粲夸克 Q=2/3Q=2/3)。定义子过程的 Mandelstam 变量: s^=(k1+k2)2,t^=(k1p1)2,u^=(k1p2)2\hat{s} = (k_1 + k_2)^2, \quad \hat{t} = (k_1 - p_1)^2, \quad \hat{u} = (k_1 - p_2)^2 满足 s^+t^+u^=2M2\hat{s} + \hat{t} + \hat{u} = 2M^2

根据第5章的 QED 计算,γγe+e\gamma \gamma \rightarrow e^+ e^- 过程对初始极化求平均、对末态自旋求和的矩阵元平方为: Mγγ2=2e4[t^M2u^M2+u^M2t^M24M2(1t^M2+1u^M2)4M4(1t^M2+1u^M2)2]\overline{|\mathcal{M}_{\gamma\gamma}|^2} = 2 e^4 \left[ \frac{\hat{t}-M^2}{\hat{u}-M^2} + \frac{\hat{u}-M^2}{\hat{t}-M^2} - 4M^2 \left( \frac{1}{\hat{t}-M^2} + \frac{1}{\hat{u}-M^2} \right) - 4M^4 \left( \frac{1}{\hat{t}-M^2} + \frac{1}{\hat{u}-M^2} \right)^2 \right]

对于 QCD 过程 γgQQˉ\gamma g \rightarrow Q\bar{Q},我们需要做以下修改:

  1. 耦合常数替换:一个光子顶点被胶子顶点取代,因此 e2(Qe)gse^2 \rightarrow (Qe) g_s。整体因子从 e4e^4 变为 Q2e2gs2=16π2Q2ααsQ^2 e^2 g_s^2 = 16\pi^2 Q^2 \alpha \alpha_s
  2. 色因子(Color Factor):我们需要对初态胶子的 8 种颜色求平均,并对末态夸克的 3 种颜色求和。由于两个费曼图都包含一个色矩阵 tat^a,色因子为: C=18a=18i,j=13tija2=18a=18Tr(tata)=18×8×12=12C = \frac{1}{8} \sum_{a=1}^8 \sum_{i,j=1}^3 |t^a_{ij}|^2 = \frac{1}{8} \sum_{a=1}^8 \text{Tr}(t^a t^a) = \frac{1}{8} \times 8 \times \frac{1}{2} = \frac{1}{2} 而在 QED 中,对应的“色因子”为 1。因此,我们需要将 QED 的结果乘以 1/21/2

结合以上两点,光子-胶子聚变子过程的矩阵元平方为: Mγg2=Q2e2gs2[t^M2u^M2+u^M2t^M24M2(1t^M2+1u^M2)4M4(1t^M2+1u^M2)2]\overline{|\mathcal{M}_{\gamma g}|^2} = Q^2 e^2 g_s^2 \left[ \frac{\hat{t}-M^2}{\hat{u}-M^2} + \frac{\hat{u}-M^2}{\hat{t}-M^2} - 4M^2 \left( \frac{1}{\hat{t}-M^2} + \frac{1}{\hat{u}-M^2} \right) - 4M^4 \left( \frac{1}{\hat{t}-M^2} + \frac{1}{\hat{u}-M^2} \right)^2 \right]

(2) 子过程的截面

子过程的微分截面为: dσ^dt^=116πs^2Mγg2\frac{d\hat{\sigma}}{d\hat{t}} = \frac{1}{16\pi \hat{s}^2} \overline{|\mathcal{M}_{\gamma g}|^2} 为了得到总子过程截面 σ^(s^)\hat{\sigma}(\hat{s}),我们需要对 t^\hat{t} 积分。由于矩阵元的结构与 γγe+e\gamma \gamma \rightarrow e^+ e^- 完全成比例,我们可以直接利用已知的 Breit-Wheeler 截面公式 σγγ\sigma_{\gamma\gamma}σ^γg(s^)=12Q2e2gs2e4σγγ(s^)=Q2αs2ασγγ(s^)\hat{\sigma}_{\gamma g}(\hat{s}) = \frac{1}{2} \frac{Q^2 e^2 g_s^2}{e^4} \sigma_{\gamma\gamma}(\hat{s}) = \frac{Q^2 \alpha_s}{2 \alpha} \sigma_{\gamma\gamma}(\hat{s}) 代入 σγγ(s^)=4πα2s^[(1+4M2s^8M4s^2)ln1+β1β(1+4M2s^)β]\sigma_{\gamma\gamma}(\hat{s}) = \frac{4\pi\alpha^2}{\hat{s}} \left[ \left( 1 + \frac{4M^2}{\hat{s}} - \frac{8M^4}{\hat{s}^2} \right) \ln\frac{1+\beta}{1-\beta} - \left( 1 + \frac{4M^2}{\hat{s}} \right) \beta \right],其中 β=14M2/s^\beta = \sqrt{1 - 4M^2/\hat{s}} 是质心系下重夸克的速度。

化简后得到子过程的总截面: σ^(s^)=2πQ2ααss^[(1+4M2s^8M4s^2)ln1+β1β(1+4M2s^)β]\hat{\sigma}(\hat{s}) = \frac{2\pi Q^2 \alpha \alpha_s}{\hat{s}} \left[ \left( 1 + \frac{4M^2}{\hat{s}} - \frac{8M^4}{\hat{s}^2} \right) \ln\frac{1+\beta}{1-\beta} - \left( 1 + \frac{4M^2}{\hat{s}} \right) \beta \right]

(3) γ\gamma-质子散射的微扰 QCD 截面

在部分子模型中,入射光子与质子中携带动量分数 xx 的胶子发生散射。设质子的四维动量为 PP,则胶子动量 k2=xPk_2 = xPγp\gamma p 系统的质心系能量平方为 s=(k1+P)22k1Ps = (k_1 + P)^2 \approx 2k_1 \cdot P。 子过程的质心系能量平方为 s^=(k1+xP)22xk1P=xs\hat{s} = (k_1 + xP)^2 \approx 2x k_1 \cdot P = xs

产生重夸克对的运动学阈值要求子过程能量必须满足 s^4M2\hat{s} \ge 4M^2,即 xs4M2xs \ge 4M^2,这给出了动量分数 xx 的积分下限: xmin=4M2sx_{\text{min}} = \frac{4M^2}{s}

将子过程截面与质子中的胶子部分子分布函数(PDF)g(x)g(x) 进行卷积,即可得到 γ+pQQˉ+X\gamma + p \rightarrow Q\bar{Q} + X 的总截面表达式: σγp(s)=4M2/s1dxg(x)σ^(s^=xs)\sigma_{\gamma p}(s) = \int_{4M^2/s}^{1} dx \, g(x) \, \hat{\sigma}(\hat{s} = xs)

将前面求得的 σ^(s^)\hat{\sigma}(\hat{s}) 代入,最终的截面表达式为:

σγp(s)=4M2s1dxg(x)2πQ2ααsxs[(1+4M2xs8M4x2s2)ln1+β1β(1+4M2xs)β]\boxed{ \sigma_{\gamma p}(s) = \int_{\frac{4M^2}{s}}^{1} dx \, g(x) \frac{2\pi Q^2 \alpha \alpha_s}{xs} \left[ \left( 1 + \frac{4M^2}{xs} - \frac{8M^4}{x^2 s^2} \right) \ln\frac{1+\beta}{1-\beta} - \left( 1 + \frac{4M^2}{xs} \right) \beta \right] }

其中 β=14M2xs\beta = \sqrt{1 - \frac{4M^2}{xs}}

17.6

Problem 17.6

peskinChapter 17

习题 17.6

来源: 第17章, PDF第597,598页


17.6 Behavior of parton distribution functions at small xx. It is possible to solve the Altarelli-Parisi equations analytically for very small xx, using some physically motivated approximations. This discussion is based on a paper of Ralston.{}^{\dagger}

(a) Show that the Q2Q^2 dependence of the right-hand side of the A-P equations can be expressed by rewriting the equations as differential equations in

ξ=loglog(Q2Λ2),\xi = \log \log (\frac{Q^2}{\Lambda^2}),

where Λ\Lambda is the value of Q2Q^2 at which αs(Q2)\alpha_s(Q^2), evolved with the leading-order β\beta function, formally goes to infinity.

(b) Since the branching functions to gluons are singular as z1z^{-1} as z0z \rightarrow 0, it is reasonable to guess that the gluon distribution function will blow up approximately as x1x^{-1} as x0x \rightarrow 0. The resulting distribution

dxfg(x)dxxdx f_g(x) \sim \frac{dx}{x}

is approximately scale invariant, and so its form should be roughly preserved by the A-P equations. Let us, then, make the following two approximations: (1) the terms involving the gluon distribution completely dominate the right-hand sides of the A-P equations; and (2) the function

g~(x,Q2)=xfg(x,Q2)\tilde{g}(x, Q^2) = x f_g(x, Q^2)

is a slowly varying function of xx. Using these approximations, and the limit x0x \rightarrow 0, show that the A-P equation for fg(x)f_g(x) can be converted to the following differential equation:

2wξg~(x,ξ)=12b0g~(x,ξ),\frac{\partial^2}{\partial w \partial \xi} \tilde{g}(x, \xi) = \frac{12}{b_0} \tilde{g}(x, \xi),

where w=log(1/x)w = \log(1/x) and b=(1123nf)b = (11 - \frac{2}{3} n_f). Show that if wξ1w \xi \gg 1, this equation has the approximate solution

g~=K(Q2)exp([48b0w(ξξ0)]1/2),\tilde{g} = K(Q^2) \cdot \exp \left( \left[ \frac{48}{b_0} w (\xi - \xi_0) \right]^{1/2} \right),

where K(Q2)K(Q^2) is an initial condition.

(c) The quark distribution at very small xx is mainly created by branching of gluons. Using the approximations of part (b), show that, for any flavor of quark, the right-hand side of the A-P equation for fq(x)f_q(x) can be approximately integrated to yield an equation for q~(x)=xfq(x)\tilde{q}(x) = x f_q(x):

ξq~(x,ξ)=23b0g~(x,ξ).\frac{\partial}{\partial \xi} \tilde{q}(x, \xi) = \frac{2}{3 b_0} \tilde{g}(x, \xi).

Show, again using wξ1w \xi \gg 1, that this equation has as its integral

q~=(ξξ027b0w)1/2K(Q2)exp([48b0w(ξξ0)]1/2).\tilde{q} = \left( \frac{\xi - \xi_0}{27 b_0 w} \right)^{1/2} K(Q^2) \cdot \exp \left( \left[ \frac{48}{b_0} w (\xi - \xi_0) \right]^{1/2} \right).

(d) Ralston suggested that the initial condition

K(Q2)=50.36(exp(ξξ0)0.957)exp[7.597(ξξ0)1/2],K(Q^2) = 50.36 (\exp(\xi - \xi_0) - 0.957) \cdot \exp \left[ -7.597 (\xi - \xi_0)^{1/2} \right],

with Q02=5 GeV2Q_0^2 = 5 \text{ GeV}^2, Λ=0.2 GeV\Lambda = 0.2 \text{ GeV}, and nf=5n_f = 5, gave a reasonable fit to the known properties of parton distributions, extrapolated into the small xx region. Use this function and the results above to sketch the behavior of the quark and gluon distributions at small xx and large Q2Q^2.

习题 17.6 - 解答


(a) A-P 方程的 Q2Q^2 依赖性与变量代换

Altarelli-Parisi (A-P) 演化方程的一般形式为: logQ2fi(x,Q2)=αs(Q2)2πx1dzzPij(z)fj(x/z,Q2)\frac{\partial}{\partial \log Q^2} f_i(x, Q^2) = \frac{\alpha_s(Q^2)}{2\pi} \int_x^1 \frac{dz}{z} P_{ij}(z) f_j(x/z, Q^2) 在领头阶 (Leading-Order) 近似下,QCD 跑动耦合常数为: αs(Q2)=4πb0log(Q2/Λ2)\alpha_s(Q^2) = \frac{4\pi}{b_0 \log(Q^2/\Lambda^2)} 其中 b0=1123nfb_0 = 11 - \frac{2}{3}n_f(题目中简记为 bb)。 引入新变量 ξ=loglog(Q2Λ2)\xi = \log \log \left(\frac{Q^2}{\Lambda^2}\right),利用链式法则计算导数: ξlogQ2=logQ2log(logQ2logΛ2)=1log(Q2/Λ2)\frac{\partial \xi}{\partial \log Q^2} = \frac{\partial}{\partial \log Q^2} \log(\log Q^2 - \log \Lambda^2) = \frac{1}{\log(Q^2/\Lambda^2)} 观察到耦合常数可以写为: αs(Q2)2π=2b0log(Q2/Λ2)=2b0ξlogQ2\frac{\alpha_s(Q^2)}{2\pi} = \frac{2}{b_0 \log(Q^2/\Lambda^2)} = \frac{2}{b_0} \frac{\partial \xi}{\partial \log Q^2} 将其代入 A-P 方程,两边消去 /logQ2\partial / \partial \log Q^2,即可将方程改写为仅对 ξ\xi 的微分方程: ξfi(x,ξ)=2b0x1dzzPij(z)fj(x/z,ξ)\boxed{ \frac{\partial}{\partial \xi} f_i(x, \xi) = \frac{2}{b_0} \int_x^1 \frac{dz}{z} P_{ij}(z) f_j(x/z, \xi) } 这表明 Q2Q^2 的依赖性被完全吸收到了演化变量 ξ\xi 中。


(b) 小 xx 极限下胶子分布的微分方程与渐近解

胶子的 A-P 方程为: ξfg(x,ξ)=2b0x1dzz[Pgg(z)fg(x/z,ξ)+fPgq(z)fq(x/z,ξ)]\frac{\partial}{\partial \xi} f_g(x, \xi) = \frac{2}{b_0} \int_x^1 \frac{dz}{z} \left[ P_{gg}(z) f_g(x/z, \xi) + \sum_f P_{gq}(z) f_q(x/z, \xi) \right] 近似 1:在小 xx 区域,胶子分布占据绝对主导,忽略夸克项 PgqP_{gq}近似 2:当 z0z \to 0 时,胶子分裂函数 Pgg(z)2CAz=6zP_{gg}(z) \approx \frac{2C_A}{z} = \frac{6}{z}。 代入上述近似,方程化为: ξfg(x,ξ)12b0x1dzz2fg(x/z,ξ)\frac{\partial}{\partial \xi} f_g(x, \xi) \approx \frac{12}{b_0} \int_x^1 \frac{dz}{z^2} f_g(x/z, \xi) 引入缓变函数 g~(x,ξ)=xfg(x,ξ)\tilde{g}(x, \xi) = x f_g(x, \xi),则 fg(x/z,ξ)=zxg~(x/z,ξ)f_g(x/z, \xi) = \frac{z}{x} \tilde{g}(x/z, \xi)。方程两边同乘 xxξg~(x,ξ)=12b0x1dzzg~(x/z,ξ)\frac{\partial}{\partial \xi} \tilde{g}(x, \xi) = \frac{12}{b_0} \int_x^1 \frac{dz}{z} \tilde{g}(x/z, \xi) 作变量代换 y=x/zy = x/z,则 dz/z=dy/ydz/z = -dy/y,积分限从 x1x \to 1 变为 1x1 \to xξg~(x,ξ)=12b0x1dyyg~(y,ξ)\frac{\partial}{\partial \xi} \tilde{g}(x, \xi) = \frac{12}{b_0} \int_x^1 \frac{dy}{y} \tilde{g}(y, \xi) 引入变量 w=log(1/x)=logxw = \log(1/x) = -\log x,对应的微分算符为 w=xx\frac{\partial}{\partial w} = -x \frac{\partial}{\partial x}。对方程两边作用 w\frac{\partial}{\partial w}2wξg~(x,ξ)=xx(12b0x1dyyg~(y,ξ))\frac{\partial^2}{\partial w \partial \xi} \tilde{g}(x, \xi) = -x \frac{\partial}{\partial x} \left( \frac{12}{b_0} \int_x^1 \frac{dy}{y} \tilde{g}(y, \xi) \right) 由变上限积分求导法则,右边等于 x(12b01xg~(x,ξ))-x \left( -\frac{12}{b_0} \frac{1}{x} \tilde{g}(x, \xi) \right),从而得到: 2wξg~(x,ξ)=12b0g~(x,ξ)\boxed{ \frac{\partial^2}{\partial w \partial \xi} \tilde{g}(x, \xi) = \frac{12}{b_0} \tilde{g}(x, \xi) }

求解渐近解: 尝试寻找形式为 g~=Kexp(S(w,ξ))\tilde{g} = K \exp(S(w, \xi)) 的解,其中 S(w,ξ)=aw(ξξ0)S(w, \xi) = \sqrt{a w (\xi - \xi_0)}。 计算二阶混合偏导数: g~ξ=g~Sξ,2g~wξ=g~(2Swξ+SwSξ)\frac{\partial \tilde{g}}{\partial \xi} = \tilde{g} \frac{\partial S}{\partial \xi}, \quad \frac{\partial^2 \tilde{g}}{\partial w \partial \xi} = \tilde{g} \left( \frac{\partial^2 S}{\partial w \partial \xi} + \frac{\partial S}{\partial w} \frac{\partial S}{\partial \xi} \right) 代入 SS 的具体形式: Sw=12a(ξξ0)w,Sξ=12awξξ0\frac{\partial S}{\partial w} = \frac{1}{2}\sqrt{\frac{a(\xi-\xi_0)}{w}}, \quad \frac{\partial S}{\partial \xi} = \frac{1}{2}\sqrt{\frac{aw}{\xi-\xi_0}} SwSξ=a4,2Swξ=14aw(ξξ0)\frac{\partial S}{\partial w} \frac{\partial S}{\partial \xi} = \frac{a}{4}, \quad \frac{\partial^2 S}{\partial w \partial \xi} = \frac{1}{4}\sqrt{\frac{a}{w(\xi-\xi_0)}}wξ1w\xi \gg 1 的极限下,交叉项 2Swξ\frac{\partial^2 S}{\partial w \partial \xi} 趋于 0,乘积项 a/4a/4 占据主导。因此: 2g~wξa4g~\frac{\partial^2 \tilde{g}}{\partial w \partial \xi} \approx \frac{a}{4} \tilde{g} 对比原微分方程,要求 a4=12b0\frac{a}{4} = \frac{12}{b_0},即 a=48b0a = \frac{48}{b_0}。故渐近解为: g~=K(Q2)exp([48b0w(ξξ0)]1/2)\boxed{ \tilde{g} = K(Q^2) \cdot \exp \left( \left[ \frac{48}{b_0} w (\xi - \xi_0) \right]^{1/2} \right) }


(c) 小 xx 极限下夸克分布的微分方程与渐近解

夸克的 A-P 方程为: ξfq(x,ξ)=2b0x1dzz[Pqq(z)fq(x/z,ξ)+Pqg(z)fg(x/z,ξ)]\frac{\partial}{\partial \xi} f_q(x, \xi) = \frac{2}{b_0} \int_x^1 \frac{dz}{z} \left[ P_{qq}(z) f_q(x/z, \xi) + P_{qg}(z) f_g(x/z, \xi) \right] 近似:小 xx 处夸克主要由胶子分裂产生,忽略 PqqP_{qq} 项。代入 q~(x)=xfq(x)\tilde{q}(x) = x f_q(x)g~(x)=xfg(x)\tilde{g}(x) = x f_g(x)ξq~(x,ξ)=2b0x1dzPqg(z)g~(x/z,ξ)\frac{\partial}{\partial \xi} \tilde{q}(x, \xi) = \frac{2}{b_0} \int_x^1 dz P_{qg}(z) \tilde{g}(x/z, \xi) 由于 g~\tilde{g} 是缓变函数,在积分域内可近似将其提出积分号,即 g~(x/z,ξ)g~(x,ξ)\tilde{g}(x/z, \xi) \approx \tilde{g}(x, \xi)。 胶子到夸克的分裂函数为 Pqg(z)=12[z2+(1z)2]P_{qg}(z) = \frac{1}{2}[z^2 + (1-z)^2]。在小 xx 极限下,积分下限 x0x \to 001Pqg(z)dz=1201(2z22z+1)dz=12(231+1)=13\int_0^1 P_{qg}(z) dz = \frac{1}{2} \int_0^1 (2z^2 - 2z + 1) dz = \frac{1}{2} \left( \frac{2}{3} - 1 + 1 \right) = \frac{1}{3} 代回原式,得到夸克分布的微分方程: ξq~(x,ξ)=23b0g~(x,ξ)\boxed{ \frac{\partial}{\partial \xi} \tilde{q}(x, \xi) = \frac{2}{3 b_0} \tilde{g}(x, \xi) }

求解渐近解: 将 (b) 中求得的 g~\tilde{g} 代入并对 ξ\xi 积分: q~(x,ξ)=ξ0ξ23b0Kexp(aw(ξξ0))dξ\tilde{q}(x, \xi) = \int_{\xi_0}^\xi \frac{2}{3b_0} K \exp\left( \sqrt{a w (\xi' - \xi_0)} \right) d\xi' 作变量代换 v=aw(ξξ0)v = \sqrt{a w (\xi' - \xi_0)},则 dξ=2vawdvd\xi' = \frac{2v}{aw} dvq~=4K3b0aw0aw(ξξ0)vevdv\tilde{q} = \frac{4K}{3b_0 a w} \int_0^{\sqrt{aw(\xi-\xi_0)}} v e^v dv 利用大 vv 极限下的渐近积分 vevdv=(v1)evvev\int v e^v dv = (v-1)e^v \approx v e^vq~4K3b0awaw(ξξ0)exp(aw(ξξ0))\tilde{q} \approx \frac{4K}{3b_0 a w} \sqrt{aw(\xi-\xi_0)} \exp\left( \sqrt{aw(\xi-\xi_0)} \right) 代入 a=48b0a = \frac{48}{b_0},化简前缀系数: 43b0(48b0)w48b0w(ξξ0)=136w48b0w(ξξ0)=48362b0w(ξξ0)=ξξ027b0w\frac{4}{3b_0 (\frac{48}{b_0}) w} \sqrt{\frac{48}{b_0} w (\xi-\xi_0)} = \frac{1}{36w} \sqrt{\frac{48}{b_0} w (\xi-\xi_0)} = \sqrt{\frac{48}{36^2 b_0 w} (\xi-\xi_0)} = \sqrt{\frac{\xi-\xi_0}{27 b_0 w}} 最终得到: q~=(ξξ027b0w)1/2K(Q2)exp([48b0w(ξξ0)]1/2)\boxed{ \tilde{q} = \left( \frac{\xi - \xi_0}{27 b_0 w} \right)^{1/2} K(Q^2) \cdot \exp \left( \left[ \frac{48}{b_0} w (\xi - \xi_0) \right]^{1/2} \right) }


(d) 小 xx 与大 Q2Q^2 区域的分布行为分析

给定初始条件 K(Q2)K(Q^2) 包含一个随 ξ\xi 衰减的指数项 exp[7.597(ξξ0)1/2]\exp[-7.597(\xi-\xi_0)^{1/2}]。 结合 (b) 和 (c) 的结果,部分子分布函数 f(x)=f~(x)/xf(x) = \tilde{f}(x)/x 的完整指数行为由下式决定: E(x,Q2)=log(1/x)+48b0log(1/x)(ξξ0)7.597ξξ0E(x, Q^2) = \log(1/x) + \sqrt{\frac{48}{b_0} \log(1/x) (\xi - \xi_0)} - 7.597\sqrt{\xi - \xi_0}

基于此,夸克和胶子分布在小 xx 和大 Q2Q^2 区域表现出以下物理行为:

  1. 双重渐近标度 (Double Asymptotic Scaling): 在极小 xx 处,w=log(1/x)w = \log(1/x) \to \infty,主导项为 exp(48b0log(1/x)loglog(Q2/Λ2))\exp\left( \sqrt{\frac{48}{b_0} \log(1/x) \log\log(Q^2/\Lambda^2)} \right)。这意味着 xfg(x)x f_g(x)xfq(x)x f_q(x)x0x \to 0 呈现出比常数更快的增长,但这种增长慢于任何幂律 xλx^{-\lambda}。分布函数在 log(1/x)\log(1/x)loglogQ2\log\log Q^2 两个方向上同时发散。

  2. Q2Q^2 增强效应: 随着 Q2Q^2 增大(即 ξ\xi 增大),指数项中的正项 w(ξξ0)\sqrt{w(\xi-\xi_0)} 使得分布函数在小 xx 处的上升变得更加陡峭。

  3. 胶子对夸克的绝对压制: 由 (c) 的结果可知,夸克与胶子分布的比例为: fq(x)fg(x)=q~(x)g~(x)1w=1log(1/x)\frac{f_q(x)}{f_g(x)} = \frac{\tilde{q}(x)}{\tilde{g}(x)} \propto \frac{1}{\sqrt{w}} = \frac{1}{\sqrt{\log(1/x)}}x0x \to 0 时,该比例趋于 0。这表明在极小 xx 区域,质子内部的动量几乎完全由胶子携带,海夸克分布虽然也在增长,但其增长速度被胶子分布以 1/log(1/x)1/\sqrt{\log(1/x)} 的因子所压制。