-
在BLFQ中,我们通过求解定态薛定谔方程[69]得到相对论性束缚态的可观测量,
$$ P^{\mu}P_{\mu} \vert \varPsi \rangle = M^2\vert \varPsi \rangle, $$ (1) 其中:
$ P^{\mu} $ 是光前坐标系下能量动量算符,表示为$ P^{\mu} = (P^+,{{P}}^{\perp},P^-) $ [69],这里$ P^{+} = P^0 + P^3 $ 、$ P^{-} = $ $ \frac{({{P}}^{\perp})^2+M^2}{P^+} $ 和$ {{P}}^{\perp} $ 分别表示纵向动量、光前哈密顿量和横向动量。由于QCD的相互作用,物理介子可以看作是正反夸克与胶子云的集合体,介子的波函数可以在Fock空间表示为以下展开形式,
$$ |{\rm{meson}} \rangle = a\vert {\rm{q}}{\bar{\rm{q}}} \rangle+b |{\rm{q}}{\bar{\rm{q}}}{\rm{g}} \rangle+c \vert {{\rm{q}}{\bar{\rm{q}}}{\rm{gg}} } \rangle+..., $$ (2) 其中,不同Fock空间的粒子数不同,但动量守恒,并且总角动量在纵向的投影守恒。由于计算资源有限,在本文中我们将Fock空间截断至前两阶
$ \vert {\rm{q}}{\bar{\rm{q}}} \rangle+ \vert {{\rm{q}}{\bar{\rm{q}}}{\rm{g}} } \rangle $ ,在后续工作中,我们可以将其扩展至更高阶。 -
我们将波函数用给定的基矢展开,其纵向我们采用了平面波
$ {\rm{e}}^{- {\rm{i}}\frac{1}{2}p^+_ix^-} $ [52],并采用箱归一化:$$ p_i^ + = \frac{{2\pi {k_i}}}{L},\quad \left\{ {\begin{array}{*{20}{l}} {{k_i} = 1,2,3...({\text{玻色子}}),}\\ {{k_i} = 1/2,3/2...({\text{费米子}}),} \end{array}} \right.$$ 其中,纵向总动量为
$ P^+ = 2\pi K/L $ ,$ L $ 为箱归一化中箱子的长度,由纵向动量守恒有$ K = \sum_i k_i $ ,下标$ i $ 为粒子标号,这里,玻色子采用周期性边界条件对应$ k $ 的值为整数(忽略零模$ k_i = 0 $ ),而费米子采用半周期性边界条件对应$ k $ 的值为半整数,规定粒子的纵向动量分数为$ x_i\equiv k_i/K $ ($ x_1 $ 、$ x_2 $ 和$ x_3 $ 分别表示夸克、反夸克和胶子的纵向动量分数)。显然,当纵向截断参数K趋于无穷时,纵向动量将趋于连续。 -
我们采用了单粒子坐标系[52]下二维谐振子的波函数作为横向基矢,通过求解如下二维谐振子动量空间的薛定谔方程,可以得到二维谐振子波函数[52],
$$ \begin{split} & P_{\varOmega}^{-} = \frac{\varOmega}{2} \left(\frac{{{p}}^2_{\perp}}{ p_i^+\varOmega} + p_i^+ \varOmega {{r}}^2_{\perp} \right),\\ &E = (2n +|m|+1)\varOmega , \end{split} $$ $$ \begin{split} &\phi_{nm}(b' = \sqrt{x_i}b,{{p}} _{\perp}) = \frac{1}{b'} \sqrt{\frac{4\pi n!}{(n+|m|)!}}\times\\ &\quad (\frac{p_{\perp}}{b'})^{|m|} {\rm{e}}^{-\frac{p^{2}_{\perp}}{2b'^2}} L^{|m|}_{n}(\frac{p^2_{\perp}}{b'^2}){\rm{e}}^{ {\rm{i}} m \theta}, \end{split} $$ (3) 其中:第一行为二维谐振子的哈密顿量;第二行和第三行分别对应二维谐振子的能量本征值和波函数。这里
$ b $ 是特征动量标度,取为$ b = \sqrt{P^+\varOmega} $ ,$ \varOmega $ 为二维谐振子的振动频率,二维谐振子波函数记作$ \phi_{nm}(\sqrt{x_i}b,{{p}}_{\perp}) $ ,$ n $ 和$ m $ 分别是径向量子数和角向量子数,而$ p_{\perp} $ 和$ \theta $ 分别对应$ p_{\perp} = |{{p}}_{\perp}| $ 和$ \theta = {\rm{arg}}({{p}}_{\perp}) $ ,$ L^{|m|}_{n}(z) $ 是拉盖尔多项式。以$ n = 5 $ 为例,不同$ m $ 取值下的二维谐振子波函数$ \phi_{nm}(b' = 1,\;{{p}}_{\perp}) $ 如图1所示。在本文中,我们选取
$ x_i $ 依赖的特征动量标度,对应的横向基矢记作$ \phi_{nm}(\sqrt{x_i}b,\;{{p}}_{\perp}) $ 。同样地,由于计算资源有限,我们对横向基矢进行能量截断$ \sum_{i} (2n_i+|m_i|+1) \leqslant N_{\rm{max}} $ ($ i $ 为粒子标号)。理论上,随着横向截断参数$ N_{\rm{max}} $ 趋于无穷,横向基矢趋于完备。 -
在BLFQ中,由于基矢空间被截断,总角动量算符
$ J $ 与哈密顿量算符$ P^- $ 不对易,空间旋转对称性只是近似满足。因此,总角动量量子数$ J $ 不是我们基矢中的好量子数,但总角动量在纵向的投影$ M_J $ 是基矢的好量子数,满足:$$ M_J = \sum\limits_i m_i +s_i, $$ (4) 其中:
$ m_i $ 为横向基矢的角向量子数;$ s_i $ 为粒子的自旋量子数[52]。在颜色空间我们选取色单态的基矢。因为在
$ | {\rm{q}}{\bar{\rm{q}}}\rangle $ 与$ |{{\rm{q}}{\bar{\rm{q}}}{\rm{g}} }\rangle $ Fock空间中各只包含有一个色单态[53],我们省略基矢的颜色量子数。综上,我们采用四个量子数$ k, n, m, s $ 对基矢粒子进行表述,其中,$ k $ 为纵向量子数,$ n $ 和$ m $ 为横向量子数,$ s $ 为自旋量子数。这些量子数满足以下对称性与截断条件:$$ \left\{ {\begin{array}{*{20}{l}} {K = \sum\nolimits_i {{k_i}} ,}\\ {{N_{{\rm{max}}}} \geqslant\sum\nolimits_i {(2{n_i} + |{m_i}| + 1)} ,}\\ {{M_J} = \sum\nolimits_i {({m_i} + {s_i})} ,} \end{array}} \right.$$ 其中:对于领头Fock空间;
$ i $ 取1和2,对应于夸克和反夸克;而对于次领头Fock空间,$ i $ 分别对应于夸克与反夸克,还包括胶子。 -
我们的光前哈密顿量
$ P^- = P^-_{\rm{QCD}}+P^-_{\rm{C}} $ ,由QCD哈密顿量$ P^-_{\rm{QCD}} $ 与唯象禁闭势$ P^-_{\rm{C}} $ 两部分构成,其中,$ P^-_{\rm{QCD}} $ 可由QCD拉氏量通过勒让德变换并考虑Fock空间的截断得到[69]:$$ \begin{split} P^-_{\rm{QCD}} =& \int {} {\rm{d}}^2 x^{\perp} {\rm{d}} x^- \frac{1}{2} \bar{\varPsi} \gamma^+ \frac{m_{\rm{q0}}^2+( {\rm{i}}\partial^{\perp})^2}{ {\rm{i}}\partial^+} \varPsi -\\& \frac{1}{2}A^i_a ( {\rm{i}}\partial^{\perp})^2 A^i_a +g \bar{\varPsi} \gamma_{\mu} T^a A^{\mu}_{a} \varPsi +\\& \frac{1}{2}g^2 \bar{\varPsi} \gamma^+ T^a \varPsi \frac{1}{( {\rm{i}}\partial^+)^2} \bar{\varPsi} \gamma^+ T^a \varPsi, \end{split} $$ (5) 这里,
$ P^-_{\rm{QCD}} $ 中前两项对应着哈密顿量的动能部分,分别是夸克和胶子动能项,后两项是相互作用部分记作$ P^-_{\rm{int}} $ ,包括矢量顶点相互作用和瞬时胶子交换相互作用。其中,$ m_{\rm{q0}} $ 为夸克的裸质量,$ x^{\perp} $ 和$ x^- $ 为横向平面和纵向坐标,$ \gamma^+ = \gamma^0+\gamma^3 $ [52],$ g $ 为强相互作用常数,$ T^a $ 为$ {SU(3)} $ 规范群的八个伴随矩阵,$ \varPsi $ 和$ A $ 分别为夸克场和胶子场。 -
考虑到夸克禁闭,我们在纵向和横向唯象地引入禁闭势
$ P^-_{\rm{C}} $ [58]为$$ P^-_{\rm{C}} P^{+} = \kappa_T^4 {{\xi}}^2_{\perp}-\frac{\kappa^4_L}{(m_{\rm{q}}+m_{\bar{\rm{q}}})^2}\partial_{x_1}(x_1(1-x_1)\partial_{x_1}), $$ (6) 其中,第一项为横向禁闭势,用
$ \kappa_T $ 描述横向禁闭势的强弱,而$ {{\xi}}_{\perp}\equiv \sqrt{x_1(1-x_1)}{{r}}_{\perp} $ 是Brodsky和de Téramond提出的全息变量[69](其中$ {{r}}_{\perp} = {{r}}_{1\perp}- {{r}}_{2\perp} $ 为相对坐标,$ {{r}}_1 $ 和$ {{r}}_2 $ 为单粒子坐标系下正反夸克的坐标);第二项为纵向禁闭势,用$ \kappa_L $ 描述纵向禁闭势的强度,这里我们取$ \kappa_T = \kappa_L = \kappa $ ,而$ m_{\rm{q}} $ 和$ m_{\bar{\rm{q}}} $ 分别为夸克和反夸克的物理质量。为简单起见,这里的禁闭势
$ P^-_{\rm{C}} $ 只作用在领头Fock空间。而对于次领头Fock空间,其哈密顿量我们目前只考虑了夸克-胶子相互作用,我们将在以后的工作中考虑次领头Fock空间中的禁闭相互作用。 -
在单粒子坐标系中,通过引入质心哈密顿量,我们将波函数中的质心运动与相对运动因子化:
$$ P'^{-} P^+ = P^-P^+ +\lambda(P^-_{\rm{cm}}P^+-2b^2I), $$ (7) 其中,质心哈密顿量
$ P^-_{\rm{cm}} $ 为[52],$$ P^-_{\rm{cm}}P^+ = {{P}}^2_{\perp}+b^4 {{R}}^2_{\perp} = (\sum\limits_i {{p}}_{\perp i})^2+b^4 (\sum\limits_i x_i {{r}}_{\perp i})^2 \;{\text{ 。}} $$ (8) 这里,我们选择一个足够大的拉格朗日乘子
$ \lambda $ 将质心激发的能谱移到高激发态,从而在低能级部分得到不含质心激发的能谱。 -
已知光前QCD哈密顿量为
$ P^{-}_{\rm{QCD}} $ ,将其作用在给定的BLFQ基矢上,可得对应的哈密顿矩阵。这里,QCD哈密顿量中夸克-胶子相互作用$ P^{-}_{\rm{int}} $ (对应式(5)右边的后两项)的形式如图2所示。为应对求解束缚态过程中遇到的发散,我们考虑了如下的重整化与正规化手续:
(1) 夸克(反夸克)质量重整化。在BLFQ的基矢中,动力学胶子不仅会贡献夸克-反夸克间的相互作用,它也会对单个夸克(反夸克)产生自能修正,如图3所示。因此,我们参照Karmanov等[71]提出的“Fock空间依赖的重整化方案”,在领头Fock空间中密顿量的动能部分,引入了质量抵消项
$ \delta m_{{\rm{q}}({\bar{\rm q}})} = m_{{\rm{q}} ({\bar{\rm q}})0}- $ $m_{{\rm{q}} ({\bar{\rm q}})} $ 来重整化夸克(反夸克)的自能发散。其具体操作与BLFQ在物理电子[72-73]和电子偶素[52, 74]的重整化类似。(2) 其它发散的正规化。为使瞬时胶子交换相互作用中瞬时胶子的截断与基矢空间中动力学胶子的截断一致,我们在瞬时胶子交换相互作用中,引入了一个宽度可调的高斯函数(
$ {\rm{e}}^{-{{p}}^2_{\perp}/b^2_ {\rm{inst}}} $ )对瞬时胶子相互作用的横向动量进行软截断[74]。此外,参照单胶子交换相互作用的有效模型[52, 58],我们在次领头Fock空间的胶子动能项中引入了一个胶子质量$ m_{\rm{g}} $ 来消除交换小动量胶子引起的红外发散。 -
在Fock空间中,
$ {{\pi}} $ 介子的光前波函数可表示为$$ \begin{split} |\varPsi_{\pi}(P) \rangle =& \sum\limits_{s_1s_2} \int {} [ {\rm{d}}^3p_1] [ {\rm{d}}^3p_2] \varPsi_{2}^{s_1s_2}(p_1,p_2)\times \\& b^{\dagger}_{s_1}(p_1)d^{\dagger}_{s_2 } (p_2)|0\rangle+\\ & \sum\limits_{s_1s_2s_3} \int {} [{\rm{d}}^3p_1] [{\rm{d}}^3p_2] [{\rm{d}}^3p_3] \varPsi_{3}^{s_1s_2s_3}(p_1,p_2,p_3)\times\\& b^{\dagger}_{s_1}(p_1)d^{\dagger}_{s_2 } (p_2)a_{s_3 }^{\dagger}(p_3)|0\rangle, \end{split} $$ (9) 通过对角化哈密顿量矩阵,可以得到
$ {{\pi}} $ 介子领头阶与次领头阶Fock空间的波函数,记作$ \varPsi_2^{{s_1}{s_2}}({p_1},{p_2}) $ 和$ \varPsi_{3}^{s_1s_2s_3}(p_1,p_2,p_3) $ 。约定$ [ {\rm{d}}^3p_i]\!\!\equiv\!\!\frac{ {\rm d} x_{i}{\rm d}^2{{p}}_{i\perp}}{(2\pi)^3} $ ,$ p_i\!\!\equiv\!\! (x_i,{{p}}_{i\perp}) $ ,$ b^{\dagger} $ 、$ d^{\dagger} $ 和$ a^{\dagger} $ 分别为夸克、反夸克和胶子的产生算符,$ |0\rangle $ 为真空态。 -
从强子到真空的矩阵元与其对应的衰变常数为
$$ \langle 0|\bar{\varPsi}(0) \gamma^+ \gamma_5 \varPsi(0)|\varPsi_{\pi}(P)\rangle = {\rm{i}} P^{+} f_{P}, $$ (10) $$ \frac{f_{P}}{2\sqrt{2N_{C}}} = \int {} [ {\rm{d}}^3p_1][ {\rm{d}}^3 p_2]\varPsi_{2}^{\uparrow\downarrow-\downarrow\uparrow}(p_1,p_2), $$ (11) 其中,式(10)左边
$ \varPsi(0) $ 表示在$ r = 0 $ 处的夸克场算符;而式(11)左边$ f_{P} $ 为衰变常数[75],它来自于$ \pi^+ \rightarrow \mu^+\; \nu $ ($ \pi^- \rightarrow \mu^-\; \bar{\nu} $ )的弱衰变,等式右边是对领头Fock空间中波函数的动量积分。约定$\varPsi_{2}^{\uparrow\downarrow-\downarrow\uparrow} (p_1,p_2)\equiv $ $ \frac{1}{\sqrt{2}}[\varPsi_{2}^{\frac{1}{2}\frac{-1}{2}}(p_1,p_2) -\varPsi_{2}^{\frac{-1}{2}\frac{1}{2}}(p_1,p_2)] $ ,色因子$ N_C = 3 $ 。 -
$ {{\pi}} $ 介子的形状因子(FFs)定义为[76]$$ I_{0,0}(Q^2)\triangleq \frac{1}{2P^+}\langle \varPsi^{J*}_{\pi}(P')|j^+(0)|\varPsi^{J}_{\pi}(P)\rangle, $$ (12) 这里,
$ J $ 为对应的总角动量对于$ {{\pi}} $ 介子$ J = 0 $ ,$ P $ 和$ P' $ 分别为对应$ {{\pi}} $ 介子的初态和末态动量,而$ \Delta = P'-P $ 为动量转移,$ Q^2 = -\Delta^2 $ ,其$j^+(0) = $ $ \bar{\varPsi}(0) \gamma^+ \varPsi(0) $ 来自于矢量流$ j^{\mu}(0) $ 。在本工作中基于领头Fock空间的波函数
$ \varPsi_2^{{s_1}{s_2}}({p_1},{p_2}) $ ,我们计算了电磁形状因子(FF)与电磁半径,它们的表达式分别为[76]$$ F(Q^2) = \frac{1}{N_ {{\rm{q}}{\bar{\rm{q}}}}}I_{0,0}(Q^2), $$ (13) $$ \langle r^2_c\rangle = -6\lim\limits_{{Q^2\rightarrow 0}}\frac{\rm d}{{\rm d} Q^2} F(Q^2), $$ (14) 其中,领头Fock空间所占的概率为
$ N_ {{\rm{q}}{\bar{\rm{q}}}} = \sum\nolimits_{{s_1}{s_2}} $ $ {\int {[{{\rm{d}}^3}{p_1}][{{\rm{d}}^3}{p_2}]} } \varPsi _2^{{s_1}{s_2}*}({p_1},{p_2})\varPsi _2^{{s_1}{s_2}}({p_1},{p_2}) $ ,电磁半径与电磁形状因子在$ Q^2 = 0 $ 处的导数相关。我们通过调禁闭势强度、强相互作用耦合常数、夸克质量、胶子质量、特征动量标度和正规化参数等参数,使得
$ {{\pi}} $ 介子的质量、衰变常数和电磁半径等物理量尽量与粒子物理数据手册上的值[77]匹配。这里,$ {{\pi}} $ 介子的物理量(质量、衰变常数和电磁半径)是输入参数,所调参数(禁闭势强度、强相互作用耦合常数、夸克质量、胶子质量、特征动量标度和正规化参数)是在经验值[62]周围遍历的。在横向截断为
$ N_{\rm{max}} = 8 $ 、纵向截断$ K = 9 $ 和$ M_J = 0 $ 时,所调的参数取值:禁闭势强度均为$ \kappa = 0.643\; {\rm{GeV}} $ ,强相互作用耦合常数$ g = 1.82 $ ,$ u $ ($ \bar{d} $ )夸克质量取$ m_{{\rm{q}}({\bar{\rm q}})} = 0.2\; {\rm{GeV}} $ ,胶子质量$ m_{\rm{g}} = 0.01\; {\rm{GeV}} $ ,特征动量标度$ b = 0.4\; {\rm{GeV}} $ ,正规化参数$ b_ {\rm{inst}} = 9.8\; {\rm{GeV}} $ ;得到的输入参数如表1所列,包括$ {{\pi}} $ 介子的质量、衰变常数和电磁半径。表 1
$ {{\pi}}$ 介子的质量、衰变常数以及电磁半径来源 $m_{\pi^+}$/MeV $f_{\pi^+}$/MeV $\sqrt{\langle r^2_c\rangle}|\pi^+ $/fm 本工作 $139.57$ $206.1$ $0.643$ PDG[77] $775.26$ $130.2\pm1.7$ $0.672\pm0.008$ $ {{\pi}} $ 介子的质量、衰变常数和(基于领头Fock空间的)电磁半径的数值计算结果(见表1),与粒子数据手册(PDG)的结果[77]相比,本结果在一定程度上描述了$ {{\pi}} $ 介子的基本性质。本工作通过匹配输入量($ {{\pi}} $ 介子的质量、衰变常数和电磁半径),得到了$ {{\pi}} $ 介子的波函数,该波函数可以计算很多可观测量,通过该波函数计算的可观测量是我们的输出结果,本工作给出了两个输出结果,一个是部分子分布函数(见下一小节),另一个是电磁形状因子。电磁形状因子的数值计算结果,如图4所示,在小
$ Q^2 $ 区域与实验数据[78-83]接近,在大$ Q^2 $ 区域我们的结果相比于实验数据有一定程度的偏离。在我们的计算中,我们考虑了$ {{\pi}} $ 介子的领头阶$ | {{\rm{q}}{\bar{\rm q}} }\rangle $ 和次领头阶$ | {{\rm{q}} {\bar{\rm q}} {\rm{g}} } \rangle $ ,其中领头阶Fock空间占的概率为$ N_ {{\rm{q}}{\bar{\rm{q}}}} = 0.887 $ ,由于领头阶占主导地位,在本工作中我们根据领头阶所占概率对波函数重新进行归一,并根据归一后的波函数计算电磁形状因子和部分子分布函数,对于次领头阶的贡献我们将在今后的工作中研究。 -
部分子分布函数(PDF)是强子中发现一个携带纵向动量分数为
$ x $ 的部分子 (夸克或胶子)的概率密度分布。由于领头阶的贡献达到$ 88.7\% $ ,为方便与之前的结果[49-50]进行对比,在本工作中,我们只计算领头Fock空间贡献的部分子分布,并根据领头阶所占概率对波函数重新进行了归一化处理,我们将在以后的工作中考虑次领头Fock空间的贡献。参考文献[58],领头Fock空间的部分子分布计算为:$$ \begin{split} f(x) = & \frac{1}{N_ {{\rm{q}}{\bar{\rm{q}}}}}\sum\limits_{s_1s_2}\int {} [ {\rm{d}}^2 {{p}}_{1 \perp}][ {\rm{d}}^2{{p}}_{2 \perp}]\times \\ &\varPsi_{2}^{s_1s_2 *}(p_1,\,p_2) \varPsi_{2}^{s_1s_2}(p_1,\,p_2), \end{split} $$ (15) 其中:
$ N_ {{\rm{q}}{\bar{\rm{q}}}} $ 为领头Fock空间所占的概率。我们对波函数的模平方横向动量进行积分自旋部分求和并归一,即可得到基于领头Fock空间的初始部分子分布。约定$ [ {\rm{d}}^2 {{p}}_{i \perp}] = \frac{ {\rm{d}}^2{{p}}_{i \perp}}{(2\pi)^2} $ 。如图5所示我们发现,在一定程度上目前的部分子分布比之前的结果宽。为进一步深入地对比,我们在下面引入了QCD演化方程,对初始部分子分布进行动力学演化。
-
在QCD中,夸克可以辐射出胶子,而胶子既可以辐射出胶子,也可以湮灭产生正反夸克对。因此,知道初始部分子分布后,可以通过QCD演化得到实验能标下的部分子分布情况。该QCD演化由Dokshitzer[84]、Gribov等[85]、Altarelli等[86]提出,简称DGLAP演化,
$$ \begin{split} \frac{\partial }{{\partial {\rm{ln}}{\mu ^{\rm{2}}}}}\!\left(\! {\begin{array}{*{20}{c}} {q(x,{\mu ^2})}\\ {g(x,{\mu ^2})} \end{array}}\! \right) \!= &\frac{{{\alpha _s}({\mu ^2})}}{{2\pi }} \int_x^1 \frac{{{\rm{d}}y}}{y} \!{\left( \!\begin{array}{*{20}{c}} {{P_{{\rm{qq}}}}(\frac{x}{y})}&{{P_{{\rm{qg}}}}(\frac{x}{y})}\\ {{P_{{\rm{gq}}}}(\frac{x}{y})}&{{P_{{\rm{gg}}}}(\frac{x}{y})} \end{array} \!\right)}\!\times \\& \left( \!{\begin{array}{*{20}{c}} {q(y,{\mu ^2})}\\ {g(y,{\mu ^2})} \end{array}}\! \right), \end{split} $$ (16) 其中:
$ P_ {\rm{qq}}(z) $ ,$ P_ {\rm{qg}}(z) $ ,$ P_ {\rm{gq}}(z) $ 和$ P_ {\rm{gg}}(z) $ 为演化核,也称劈裂函数,$ q(x,\mu^2) $ 和$ g(x,\mu^2) $ 分别是夸克和胶子的部分子分布,而$ \alpha_s(\mu^2) $ 是强相互作用中的跑动耦合函数,对应的能标由$ \mu^2 $ 表示。为与之前的工作[49-50]一致,这里采用高阶微扰部分子演化包(HOPPET)[87]来数值求解上式DGLAP演化方程,其阶数也考虑到次次领头阶(NNLO)。将领头Fock空间里的(归一化后的)初始部分子分布作为输入(
$ q(x,\mu^2_0) = f(x) $ ),由于其初始价夸克携带全部纵向动量,选择起始能标为$ \mu_0^2 = 0.24\; {{\rm{GeV}}^2} $ [49-50],演化至实验能标$ \mu^2 = 16\; {{\rm{GeV}}^2} $ ,如图6所示,紫色实线是目前的结果,其它部分可参考之前的文章[49-50]。我们发现,演化后胶子和海夸克的部分子分布与原先结果接近,而价夸克分布相比原先结果对E-0615的实验结果[6]符合得更好。在实验能标下,价夸克的部分子分布在大$ x $ 区的下降趋势为$ (1-x)^{1.3} $ ,相比原先结果$ (1-x)^{1.44} $ 更倾向于支持线性下降[1-5]。基于演化后的部分子分布,我们可以计算Drell-Yan过程
$ \pi^- N \rightarrow \mu^+\mu^- X $ 的微分散射截面。 -
在强子与强子的碰撞中,湮灭一对分别来自两个强子的夸克与反夸克对并产生一对轻子的过程,称为Drell-Yan过程,如图7所示,该过程是量子色动力学中最基本的过程之一。
在该过程中,
$ p_{1,2} $ 和$ l_{1,2} $ 分别是入射强子和出射轻子的四动量,进而该过程可由如下物理量描述[22, 88-93]:$$ \begin{split} & s = {({p_1} + {p_2})^2},\quad \quad \quad \quad q = {l_1} + {l_2},\\ &{m^2} = {q^2},\quad \quad \quad \quad \quad Y = \frac{1}{2} \ln \frac{{{q_0} + {q_3}}}{{{q_0} - {q_3}}} ,\;\;\\ &{x_{\rm{F}}} = {x_1} - {x_2},\quad \quad \quad \quad \quad \tau = \frac{{{m^2}}}{s},\\ & z = \frac{{{m^2}}}{{\hat s}} = \frac{\tau }{{{x_1}{x_2}}},\quad y = \frac{{\frac{{{x_1}}}{{{x_2}}}{{\rm e}^{ - 2Y}} - z}}{{(1 - z)(1 + \frac{{{x_1}}}{{{x_2}}}{{\rm e}^{ - 2Y}})}} , \end{split} $$ 其中:
$ \sqrt{s} $ 为质心系能量;$ q $ 为虚光子的四动量;$ m $ 为轻子对的不变质量;$ Y $ 为虚光子的快度;$ x_{\rm{F}} $ 为费曼变量;$ \sqrt{\tau} $ 为轻子对不变质量与质心系能量之比;$ \sqrt{\hat{s}} $ 是湮灭的夸克对的质心系能量;而$ z $ 则是轻子对不变质量与湮灭的夸克对质心系能量之比;$ y $ 是描述该过程引入的一个辅助变量。那么,湮灭的强子对的纵向动量可表示为$$ \begin{split} x_{1} =& \sqrt{\frac{\tau}{z}\frac{1-(1-y)(1-z)}{1-y(1-z)}} {\rm{e}}^{Y}, \nonumber\\ x_{2} =& \sqrt{\frac{\tau}{z}\frac{1-y(1-z)}{1-(1-y)(1-z)}} {\rm{e}}^{-Y} \;{\text{ 。}} \end{split} $$ 因此,我们可以定义固定阶数的Drell-Yan微分散射截面计算公式为[22, 90-92]:
$$ \begin{split} \frac{m^3 {\rm{d}}^2{\bf \sigma}}{ {\rm{d}} m\, {\rm{d}} Y} =& \frac{8\pi\alpha^2}{9 }\frac{m^2}{s}\sum\limits_{ij} \int {} {\rm{d}} x_1 {\rm{d}} x_2 \widetilde{C}_{ij}(x_1,\,x_2,\,m^2,\,\mu^2)\times\\ & f_{i/\pi}(x_1,\mu^2) f_{j/N}(x_2,\,\mu^2), \end{split} $$ (17) 其中:
$ \alpha $ 为电磁相互作用中的耦合常数;$ f_{i/\pi} $ 和$ f_{j/N} $ 分别为$ {{\pi}} $ 介子的部分子分布和对应原子核的部分子分布;$ \widetilde{C}_{ij} $ 则是相应能标下的微扰系数[90-92]。这里我们计算到次领头阶,与原先计算同阶[50],求和$ ij $ 包括$ {\rm{q}}{\bar{\rm{q}}} $ 、$ {{\bar{\rm q}}{\rm{q}}} $ 、$ {{\rm{qg}}} $ 、gq、$ {{\bar{\rm q}}{\rm{g}}} $ 和$ {{\rm{g}}{\bar{\rm q}}} $ ,其详细表达式可参考上篇文章的附录[50]。同样地,我们采用了2015年的nCTEQ原子核部分子分布[94]作为对应靶核部分子分布的输入。在图8~10中,紫色实线为目前的结果,黑色虚线为BLFQ-NJL模型的计算结果①[50],实验数据包括:费米国家实验室(FNAL)的E-0615实验数据[6]、E-0326实验数据[95]和E-0444实验数据[96],欧洲核子中心(CERN)的NA-003实验数据[7]、NA-010实验数据[97]、WA-011实验数据[98]和WA-039实验数据[99],实验数据对应的
$ \pi^- $ 束流能量、靶核和质心系能量$ \sqrt{s} $ 见表2。表 2
$ \pi^- N \rightarrow \mu^+\mu^- X$ 的实验数据信息表进行初步的对比,我们发现,目前的结果与BLFQ-NJL模型的结果[50]相近,都能很好地描述实验数据。进行更详细的对比,我们发现,在较大轻子对不变质量 (
$ \sqrt{\tau} $ 、$ \tau $ 和$ m $ )区,目前的结果比原结果有一定程度的增大,并且与E-0615、NA-003、WA-011和WA-039实验数据更接近,但与NA-010、E-0326和E-0444实验数据偏离稍大。
-
摘要:
为研究
介子的性质,通过基矢光前量子化方法获得
介子的光前波函数。基失光前量子化是一种在哈密顿量体系下基于量子场论的非微扰方法。在哈密顿量中我们考虑了动能项、基于全息色动力学的横向禁闭势、与横向禁闭势互补的纵向禁闭势和基于QCD的夸克-胶子相互作用。我们的基矢空间包括最低阶的两个Fock空间,即领头阶
与次领头阶
。根据所得的光前波函数我们计算了
介子衰变常数以及(基于领头Fock空间的)电磁半径,这些结果与粒子数据手册(PDG)上的结果相近。此外,我们计算了(基于领头Fock空间的)
介子部分子分布,QCD演化后,与原先的结果相近(蓝江山等,Phys Rev Lett, 2019, 122: 172001.),能够很好地描述费米国家实验室(FNAL)与欧洲核子中心(CERN)的实验数据。
-
关键词:
- ${{\pi}}$介子 /
- 基矢光前量子化 /
- 部分子分布 /
- QCD演化
Abstract:We study the properties of the pion through the light-front wave function (LFWF) obtained from the basis light-front quantization (BLFQ) approach. BLFQ is a nonperturbative approach to quantum field theory based on the Hamiltonian formalism. Our Hamiltonian contains the kinetic energy terms, a transverse confining potential motivated by holographic quantum chromodynamics (QCD), a complementary longitudinal confining potential, and the quark-gluon interactions based on QCD. Our basis space includes the lowest two Fock sectors, namely
and
. The obtained pion decay constant and electromagnetic radius (at the leading order Fock sector) are comparable to those from the particle data group (PDG). Next, we calculate the parton distribution function (PDF) of the pion based on the leading Fock sector LFWF. After QCD evolution, the resulting PDF is similar to the previous result (Lan et al., Phys Rev Lett, 2019, 122: 172001.), and agrees with the experimental data from the Fermi National Accelerator Laboratory (FNAL), as well as those from the European Organization for Nuclear Research (CERN).
-
Key words:
- pion /
- basis light-front quantization /
- parton distribution function /
- QCD evolution
-
图 2 相互作用部分哈密顿量示意图
左上费曼图为领头Fock空间的瞬时胶子交换相互作用,对应哈密顿矩阵中的对角部分;右上费曼图为领头Fock空间与次领头Fock空间交叉作用部分,具体为夸克到夸克-胶子的矢量顶点相互作用,对应哈密顿矩阵中的上角部分;同样地,左下费曼图为相应夸克-胶子到夸克的矢量顶点相互作用,对应哈密顿矩阵中的下角部分;右下费曼图为0,是由于根据“规范原理”[70]其相互作用涉及到$ | {{\rm{q}}{\bar{\rm q}}{\rm{gg}}}\rangle$Fock空间,而我们的基矢空间只包括领头阶与次领头阶Fock展开。
表 1
$ {{\pi}}$ 介子的质量、衰变常数以及电磁半径来源 $m_{\pi^+}$ /MeV$f_{\pi^+}$ /MeV$\sqrt{\langle r^2_c\rangle}|\pi^+ $ /fm本工作 $139.57$ $206.1$ $0.643$ PDG[77] $775.26$ $130.2\pm1.7$ $0.672\pm0.008$ -
[1] BORDALO P, BUSSON P, KLUBERG L, et al. Phys Lett B, 1987, 193: 368. doi: 10.1016/0370-2693(87)91253-6 [2] FREUDENREICH K. Int J Mod Phys A, 1990, 5: 3643. doi: 10.1142/S0217751X90001586 [3] SUTTON P J, MARTIN A D, ROBERTS R G. Phys Rev D, 1992, 45: 2349. doi: 10.1103/PhysRevD.45.2349 [4] GLUCK M, REYA E, SCHIENBEIN I. Eur Phys J C, 1999, 10: 313. doi: 10.1007/s100529900124 [5] WIJESOORIYA K, REIMER P E, HOLT R J. Phys Rev C, 2005, 72: 065203. doi: 10.1103/PhysRevC.72.065203 [6] CONWAY J S, ADOLPHSEN C E, ALEXANDER J P, et al. Phys Rev D, 1989, 39: 92. doi: 10.1103/PhysRevD.39.92 [7] BADIER J, BOUCROT J, BOUROTTE J, et al. Z Phys C, 1983, 18: 281. doi: 10.1007/BF01573728 [8] AICHER M, SCHAFER A, VOGELSANG W. Phys Rev Lett, 2010, 105: 252003. doi: 10.1103/PhysRevLett.105.252003 [9] WATANABE A, SANWADA T, KAO C W. Phys Rev D, 2018, 97: 074015. doi: 10.1103/PhysRevD.97.074015 [10] HECHT M B, ROBERTS C D, SCHMIDT S M. Phys Rev C, 2001, 63: 025213. doi: 10.1103/PhysRevC.63.025213 [11] NAM S I. Phys Rev D, 2012, 86: 074005. doi: 10.1103/PhysRevD.86.074005 [12] DETMOLD W, MELNITCHOUK W, THOMAS A W. Phys Rev D, 2003, 68: 034025. doi: 10.1103/PhysRevD.68.034025 [13] HOLT R J, ROBERTS C D. Rev Mod Phys, 2010, 82: 2991. doi: 10.1103/RevModPhys.82.2991 [14] PUMPLIN J, STUMP D R, HUSTON J, et al. JHEP, 2002, 0207: 012. doi: 10.1088/1126-6708/2002/07/012 [15] BALL R D, BERTONE V, CARRAZZA S, et al. Eur Phys J C, 2017, 77: 663. doi: 10.1140/epjc/s10052-017-5199-5 [16] ALEKHIN S, BLUMLEIN J, MOCH S, et al. Phys Rev D, 2017, 96: 014011. doi: 10.1103/PhysRevD.96.014011 [17] DULAT S, HOU T J, GAO J, et al. Phys Rev D, 2016, 93: 033006. doi: 10.1103/PhysRevD.93.033006 [18] HARLAND-LANG L A, MARTIN A D, MOTYLINSKI P, et al. Eur Phys J C, 2015, 75: 204. doi: 10.1140/epjc/s10052-015-3397-6 [19] BEDNAR K D, CLOET I C, TANDY P C. arXiv: 1811.12310[nucl-th]. [20] AGUILAR A C, AHMED Z, AIDALA C, et al. Eur Phys J A, 2019, 55: 190. doi: 10.1140/epja/i2019-12885-0 [21] DING M, RAYA K, BINOSI D, et al. arXiv: 1905.05208[nuclth]. [22] BARRY P C, SATO N, MELNITCHOUK W, et al. Phys Rev Lett, 2018, 121: 152001. doi: 10.1103/PhysRevLett.121.152001 [23] DAVIDSON R M, RUIZ ARRIOLA E. Acta Phys Polon B, 2002, 33: 1791. [24] NGUYEN T, BASHIR A, ROBERTS C D, et al. Phys Rev C, 2011, 83: 062201. doi: 10.1103/PhysRevC.83.062201 [25] RANDYUSHKIN A V. Phys Rev D, 2017, 95: 056020. doi: 10.1103/PhysRevD.95.056020 [26] BRONIOWSKI W, RUIZ ARRIOLA E. Phys Lett B, 2017, 773: 385. doi: 10.1016/j.physletb.2017.08.055 [27] XU S S, CHANG L, ROBERTS C D, et al. Phys Rev D, 2018, 97: 094014. doi: 10.1103/PhysRevD.97.094014 [28] HUTAURUK P T P, CLOET I C, THOMAS A W. Phys Rev C, 2016, 94: 035201. doi: 10.1103/PhysRevC.94.035201 [29] CHEN C, CHANG L, ROBERTS C D, et al. Phys Rev D, 2016, 93: 074021. doi: 10.1103/PhysRevD.93.074021 [30] SHI C, MEZRAG C, ZONG H S. Phys Rev D, 2018, 98(5): 054029. doi: 10.1103/PhysRevD.98.054029 [31] LIN H W, NOCERA E R, OLNESS F, et al. Prog Part Nucl Phys, 2018, 100: 107. doi: 10.1016/j.ppnp.2018.01.007 [32] BRONIOWSKI W, RUIZ A E, GOLEC-BIERNAT K. Phys Rev D, 2008, 77: 034023. doi: 10.1103/PhysRevD.77.034023 [33] GUTSCHE T, LYUBOVITSKIJ V E, SCHMIDT I, et al. J Phys G, 2015, 42: 095005. doi: 10.1088/0954-3899/42/9/095005 [34] GUTSCHE T, LYUBOVITSKIJ V E, SCHMIDT I, et al. Phys Rev D, 2014, 89: 054033. Erratum: [Phys Rev D, 2015, 92: 019902]. [35] DE TERAMOND F G, LIU T, SUFIAN R S, et al. [HLFHS Collaboration]. Phys Rev Lett, 2018, 120: 182001. doi: 10.1103/PhysRevLett.120.182001 [36] BROMMEL D, DIEHL M, GOCKELER M, et al. PoS LATTICE, 2007, 2007: 140. doi: 10.22323/1.042.0140 [37] MARTINELLI G, SACHRAJDA C T. Nucl Phys B, 1988, 306: 865. doi: 10.1016/0550-3213(88)90445-2 [38] ABDEL-REHIM A, ALEXANDROU C, CONSTANTINOU M, et al. Phys Rev D, 2015, 92: 114513. Erratum: [Phys Rev D, 2016, 93: 039904]. [39] SUFIAN R S, KARPIE J, EHERER C, et al. Phys Rev D, 2019, 99: 074507. doi: 10.1103/PhysRevD.99.074507 [40] AHMADY M, MONDAL C, SANDAPEN R. Phys Rev D, 2018, 98: 034010. doi: 10.1103/PhysRevD.98.034010 [41] OEHM M, ALEXANDROU C, CONSTANTINOU M, et al. Phys Rev D, 2019, 99: 014508. doi: 10.1103/PhysRevD.99.014508 [42] FREDERICO T, MILLER A G. Phys Rev D, 1994, 50: 210. doi: 10.1103/PhysRevD.50.210 [43] SHIGETANI T, SUZUKI K, TOKI H. Phys Lett B, 1993, 308: 383. doi: 10.1016/0370-2693(93)91302-4 [44] YUAN F. Phys Rev D, 2004, 69: 051501. doi: 10.1103/PhysRevD.69.051501 [45] BRODSKY S J, YUAN F. Phys Rev D, 2006, 74: 094018. doi: 10.1103/PhysRevD.74.094018 [46] BERGER E L, BRODSKY S J. Phys Rev Lett, 1979, 42: 940. doi: 10.1103/PhysRevLett.42.940 [47] FARRAR G R, JACKSON D R. Phys Rev Lett, 1979, 43: 246. doi: 10.1103/PhysRevLett.43.246 [48] MELNITCHOUK W. Eur Phys J A, 2003, 17: 223. doi: 10.1140/epja/i2003-10006-6 [49] LAN J, MONDAL C, JIA S, et al. Phys Rev Lett, 2019, 122: 172001. doi: 10.1103/PhysRevLett.122.172001 [50] LAN J, MONDAL C, JIA S, et al. Phys Rev D, 2020, 101: 034024. doi: 10.1103/PhysRevD.101.034024 [51] VARY J P, HONKANEN H, LI J, et al. Phys Rev C, 2010, 81: 035205. doi: 10.1103/PhysRevC.81.035205 [52] WIECKI P, LI Y, ZHAO X, et al. Phys Rev D, 2015, 91: 105009. doi: 10.1103/PhysRevD.91.105009 [53] LI J. Light front Hamiltonian and Its Application in QCD[D]. Iowa: Iowa State University, 2009. https://lib.dr.iastate.edu/etd/11067 [54] LI Y, MARIS P, ZHAO X, et al. Phys Lett B, 2016, 758: 118. doi: 10.1016/j.physletb.2016.04.065 [55] BRODSKY S J, PAULI H C, PINSKY S S. Phys Rept, 1998, 301: 299. doi: 10.1016/S0370-1573(97)00089-6 [56] ZHAO X, ILDERTON A, MARIS P, et al. Phys Lett B, 2013, 726: 856. doi: 10.1016/j.physletb.2013.09.030 [57] ZHAO X, HONKANEN H, MARIS P, et al. Phys Lett B, 2014, 737: 65. doi: 10.1016/j.physletb.2014.08.020 [58] LI Y, MARRIS P, VARY J P. Phys Rev D, 2017, 96: 016022. doi: 10.1103/PhysRevD.96.016022 [59] LEITAO S, LI Y, MARIS P, et al. Eur Phys J C, 2017, 77: 696. doi: 10.1140/epjc/s10052-017-5248-0 [60] LI M, LI Y, MARIS P, et al. Phys Rev D, 2018, 98: 034024. doi: 10.1103/PhysRevD.98.034024 [61] TANG S, LI Y, MARIS P, et al. Phys Rev D, 2018, 98: 114038. doi: 10.1103/PhysRevD.98.114038 [62] JIA S, VARY J P. Phys Rev C, 2019, 99: 035206. doi: 10.1103/PhysRevC.99.035206 [63] MONDAL C, XU S, LAN J, et al. PoS DIS, 2019, 2019: 190. doi: 10.22323/1.352.0190 [64] MONDAL C, XU S, LAN J, et al. arXiv: 1911.10913[hep-ph]. [65] DU W, LI Y, ZHAO X, et al. arXiv: 1911.10762[nucl-th]. [66] LAN J, MONDAL C, LI M, et al. arXiv: 1911.11676[nucl-th]. [67] TANG S, LI Y, MARRIS P, et al. arXiv: 1912.02088[nucl-th]. [68] HU B, ILDERTON A, ZHAO X. arXiv: 1911.12307[nucl-th]. [69] BRODSKY S J, DE TERAMOND G F, ERLICH J, et al. Phys Rept, 2015, 584: 1. doi: 10.1016/j.physrep.2015.05.001 [70] TANG A C, BRODSKY S J, PAULI H C. Phys Rev D, 1991, 44: 1842. doi: 10.1103/PhysRevD.44.1842 [71] LI Y, KARMANOV V A, MARRIS P, et al. Phys Lett B, 2015, 748: 278. doi: 10.1016/j.physletb.2015.07.014 [72] ZHAO X, HONKANEN H, MARRIS P, et al. Few Body Syst, 2012, 52: 339. doi: 10.1007/s00601-011-0273-2 [73] WIECKI P, LI Y, ZHAO X, et al. Few Body Syst, 2015, 56: 489. doi: 10.1007/s00601-015-0962-3 [74] FU K, ZHAO H, ZHAO X, et al. Manuscript in Preparation, 2020. [75] LEPAGE G P, BRODSKY S J. Phys Rev D, 1980, 22: 2157. doi: 10.1103/PhysRevD.22.2157 [76] ADHIKARI L, LI Y, ZHAO X, et al. Phys Rev C, 2016, 93: 055202. doi: 10.1103/PhysRevC.93.055202 [77] TANABASHI M, HAGIWARA K, HIKASA K, et al. Phys Rev D, 2018, 98: 030001. doi: 10.1103/PhysRevD.98.030001 [78] AMENDOLIA S R, ARIK M, BADELEK B, et al. Nucl Phys B, 1986, 277: 168. doi: 10.1016/0550-3213(86)90437-2 [79] BEBEK C J, BROWN C N, HERZLINGER M, et al. Phys Rev D, 1974, 9: 1229. doi: 10.1103/PhysRevD.9.1229 [80] BEBEK C J, BROWN C N, HERZLINGER M, et al. Phys Rev D, 1976, 13: 25. doi: 10.1103/PhysRevD.13.25 [81] BEBEK C J, BROWN C N, HOLMES S D, et al. Phys Rev D, 1978, 17: 1693. doi: 10.1103/PhysRevD.17.1693 [82] VOLMER J, ABBOTT D, ANKLIN H, et al. Phys Rev Lett, 2001, 86: 1713. doi: 10.1103/PhysRevLett.86.1713 [83] HORN T, ANIOL K, ARRINGTON J, et al. Phys Rev Lett, 2006, 97: 192001. doi: 10.1103/PhysRevLett.97.192001 [84] DOKSHITZER Y L. Sov Phys JETP, 1997, 46: 641. Erratum: [Zh Eksp Teor Fiz 1977, 73: 1216]. http://inspirehep.net/record/126153. [85] GRIBOV V N, LIPATOV L N. Sov J Nucl Phys, 1972, 15: 438. http://inspirehep.net/record/73449. [86] ALTARELLI G. PARISI G Nucl Phys B, 1977, 126: 298. doi: 10.1016/0550-3213(77)90384-4 [87] SALAM G P, ROJO J. Comput Phys Commun, 2009, 180: 120. doi: 10.1016/j.cpc.2008.08.010 [88] DRELL S D, YAN T M. Phys Rev Lett, 1970, 25: 316. Erratum: [Phys Rev Lett 1970, 25: 902]. doi: 10.1103/PhysRevLett.25.316, 10.1103/PhysRevLett.25.902.2. [89] MCGAUGHEY P L, MOSS J M, ALDE D M, et al. Phys Rev D, 1994, 50: 3038. Erratum: [Phys Rev D, 1999, 60: 119903]. doi: 10.1103/PhysRevD.50.3038, 10.1103/PhysRevD.60.119903. [90] ANASTASIOU C, DIXON L J, MELNIKOV K, et al. Phys Rev Lett, 2003, 91: 182002. doi: 10.1103/PhysRevLett.91.182002 [91] ANASTASIOU C, DIXON L J, MELNIKOV K, et al. Phys Rev D, 2004, 69: 094008. doi: 10.1103/PhysRevD.69.094008 [92] BECHER T, NEUBERT M, XU G. JHEP, 2008, 0807: 030. doi: 10.1088/1126-6708/2008/07/030 [93] PASQUINI B, SCHWEITZER P. Phys Rev D, 2014, 90: 014050. doi: 10.1103/PhysRevD.90.014050 [94] KOVARIK K, KUSINA A, JEZO T, et al. Phys Rev D, 2016, 93: 085037. doi: 10.1103/PhysRevD.93.085037 [95] GREENLEE H B, FRISCH H J, JOHNSON K F, et al. Phys Rev Lett, 1985, 55: 1555. doi: 10.1103/PhysRevLett.55.1555 [96] ANDERSON K J, COLEMAN R N, HOGAN G E, et al. Phys Rev Lett, 1979, 42: 944. doi: 10.1103/PhysRevLett.42.944 [97] BETEV B, BLAISING J J, BORDALO P, et al. [NA10 Collaboration]. Z Phys C, 1985, 28: 9. doi: 10.1007/BF01550243 [98] BARATE R, BAREYRE P, BONAMY P, et al. Phys Rev Lett, 1979, 43: 1541. doi: 10.1103/PhysRevLett.43.1541 [99] CORDEN M, DOWELL J D, GARVEY J, et al. Phys Lett, 1980, 96B: 417. doi: 10.1016/0370-2693(80)90800-X