-
因与氢原子非常相似,单重介子系统的研究一直受人们关注。二十多年来,重介子家族不断扩大,这得益于B工厂实验的BaBar、Belle、CLEO和最近的LHCb合作组[1]发现了大量激发态粲介子和一系列激发态底介子。尽管如此,一些粲介子和底介子激发态的结构仍不清楚,需要人们去确定。例如,
$ B_{sJ}^ * \left( {5850} \right) $ 可以理解为$ s\bar b $ 态,但需进一步确认。$ {D_J}{\left( {3000} \right)^0} $ [2-3]的情况与之类似,其结构仍未确定。D0合作组[4-5]已观测到$ X{\left( {5568} \right)^ \pm } $ 态,但目前该粒子仍未被LHCb合作组[4-5]所发现。一些奇异粲介子,例如最近发现的$ D_{s0}^ * \left( {2317} \right) $ ,$ {D_{s1}}\left( {2460} \right) $ 和$ {X_{0,1}}\left( {2900} \right) $ [6],被认为含有奇异的成分。为此,包括格点模拟[7-8]在内的很多方法[9-21]被用来研究重轻介子谱,详情见文献[22-23]等。在夸克模型中[24],激发态粲介子和底介子的窄质量劈裂一直有争议。在本文中,我们结合线性Regge轨迹和相对论夸克模型重新考察了自旋宇称$ {J^ + } $ ($ J = 0 $ ,$ 1 $ ,$ 2 $ )的单重(SH)介子(D,$ {D_s} $ ,B和$ {B_s} $ )的P波质量,并解释它们的窄质量劈裂主要是由介子内轻夸克的相对论效应引起。对P波单重介子的质量计算表明,底介子B和$ {B_s} $ 的两个实验缺失态分别是质量5 659 MeV的$ {0^ + } $ 态和质量为5 788 MeV的$ {0^ + } $ 态。与实验数据相比有质量平移的介子$ {D_{s0}}\left( {2317} \right) $ 和$ {D_{s1}}\left( {2460} \right) $ 可能存在耦合道效应。在表1和表2中,我们列出了粲介子和底介子的S波和P波质量的实验值[1]以及典型的夸克模型预言[9]。
State $ {J^P} $ 介子 质量/MeV GI[9]/MeV 介子 质量/MeV GI[9]/MeV $ 1{}^1{S_0} $ $ {0^ - } $ $ {D^ \pm } $ 1 869.6 1 880 $ {D_s} $ 1 968.4 1 980 $ 1{}^3{S_1} $ $ {1^ - } $ $ {D^ * }{\left( {2\;010} \right)^ \pm } $ 2 010.3 2 040 $ D_s^ * {\text{[}}{{\text{?}}^{\text{?}}}{\text{]}} $ 2 112.2 2 130 $ 1{}^3{P_0} $ $ {0^ + } $ $ D_0^*\left( {2\;300} \right) $ 2 343(10) 2 400 $ D_{s0}^ * {\left( {2\;317} \right)^ \pm } $ 2 317.8 2 480 $ 1{P_1} $ $ {1^ + } $ $ {D_1}{\left( {2\;430} \right)^0} $ 2 412(9) 2 440 $ {D_{s1}}{\left( {2\;460} \right)^ \pm } $ 2 459.5 2 530 $ 1{P_1} $ $ {1^ + } $ $ {D_1}{\left( {2\;420} \right)^ \pm } $ 2 422.1 2 490 $ {D_{s1}}{\left( {2\;536} \right)^ \pm } $ 2 535.1 2 570 $ 1{}^3{P_2} $ $ {2^ + } $ $ D_2^ * \left( {2\;460} \right) $ 2 461.1 2 500 $ D_{s2}^ * \left( {2\;573} \right) $ 2 569.1 2 590 ?:表示安排于此的量子数为夸克模型预言但未被实验确认;这里小于5 MeV的误差未列出。 State $ {J^P} $ 介子 质量/MeV GI[9]/MeV 介子 质量/MeV GI[9]/MeV $ 1{}^1{S_0} $ $ {0^ - } $ $ {B^0} $ 5 279.7 1 880 $ {B_s} $ 5 366.9 5 390 $ 1{}^3{S_1} $ $ {1^ - } $ $ {B^ * } $ 5 324.7 2 040 $ B_s^ * $ 5 415.4 5 450 $ 1{}^3{P_0} $ $ {0^ + } $ $ 1{P_1} $ $ {1^ + } $ $B_J^ * \left( {5\;732} \right){\text{[} }{ {\text{?} }^{\text{?} } }{\text{]} }$ 5 698(8) $B_{sJ}^ * \left( {5\; 850} \right){\text{[} }{ {\text{?} }^{\text{?} } }{\text{]} }$ 5 853(15) $ 1{P_1} $ $ {1^ + } $ ${B_1}{\left( {5\; 721} \right)^0}$ 5 725.9 ${B_{s1} }{\left( {5\;830} \right)^0}$ 5 828.7 $ 1{}^3{P_2} $ $ {2^ + } $ $B_2^ * {\left( {5\; 747} \right)^0}$ 5 737.2 5 800 $B_{s2}^*{\left( {5 \,840} \right)^0}$ 5 839.9 5 880 ?:表示安排于此的量子数为夸克模型预言但未被实验确认;这里小于5 MeV的误差未列出。 -
我们将P波单重介子
$ \bar qQ $ ($Q = c,\,b,\,{\text{ }}q = {u \mathord{\left/ {\vphantom {u d}} \right. } d},\,s$ )的质量分为两个部分,即$ M = \bar M + \Delta M $ ,其中$ \bar M $ 为自旋无关部分的质量,$ \Delta M $ 为自旋相关部分的质量。为计算$ \bar M $ ,我们采用文献[25-26]中的线性Regge关系,它由QCD弦模型所推导出并成功应用于单重强子中[27-28]。设$ {\bar M_L} $ 是轨道激发为$ L $ 的单重介子的质量,则Regge轨迹关系式为[25-26]其中:a是重夸克和轻反夸克之间的QCD弦张力系数;重夸克质量为
$ {M_Q} $ ;轻反夸克质量为$ {m_q} $ 。 当$ Q = c\left( b \right) $ 时,$ {m_{{\text{bare}}Q}} = 1.275\left( {4.18} \right)\;{\text{GeV}} $ ,这里$ {m_{{\text{bare}}Q}} $ 表示重夸克Q的裸质量。对于单重介子,用轻反夸克$ \bar q $ 的有效质量$ {m_q} $ 来代替文献[25]中轻的双夸克质量$ {m_d} $ , 即可得到方程(1)。 在方程(1)中,截距$ {a_I} = \left( {{m_q} + } \right.{M_Q} - $ ${\left. {{{m_{{\text{bare}}Q}^2} \mathord{\big/ {\vphantom {{m_{{\text{bare}}Q}^2} {{M_Q}}}} } {{M_Q}}}} \right)^2}$ 和反夸克质量$ {m_q} $ 与重夸克非相对论动能${{P_Q^2} \mathord{/ {\vphantom {{P_Q^2} {{M_Q}}}} } {{M_Q}}}$ 有关,因为截距也可以写为$ {a_I} = \left( {{m_q} + } \right.P_Q^2 $ $ {\left. {{{} \mathord{\left/ {\vphantom {{} {{M_Q}}}} \right. } {{M_Q}}}} \right)^2} $ ,由此得到$ {P_Q} \equiv {M_Q}{v_Q} $ ,其中${v_Q} = {\left( {1 - {{m_{{\text{bare}}Q}^2} \mathord{\big/ {\vphantom {{m_{{\text{bare}}Q}^2} {M_Q^2}}} } {M_Q^2}}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}$ 为重夸克速度,它在重夸克极限下守恒。当L=0时,方程(1)可简化为重夸克对称性[29]预言的基态质量方程:注意,重夸克极限
$ {m_{{\text{bare}}Q}} \to \infty $ 意味着,$ {M_Q} \cong {m_{{\text{bare}}Q}} $ $ \to \infty \left( {{v_Q} \to 0} \right) $ ,而$ {P_Q} \equiv {M_Q}{v_Q} $ 是守恒的。结合单重重子和单重介子的实验质量可得到相应的夸克质量和弦张力系数$ a $ ,结果在文献[30]中给出,如表3所列,其中$ {a_{c\bar n}} $ ,$ {a_{c\bar s}} $ ,$ {a_{b\bar n}} $ 和$ {a_{b\bar s}} $ 分别代表D介子、$ {D_s} $ 介子、B介子和$ {B_s} $ 介子的弦张力系数。Parameters $ {M_c}/{\rm GeV} $ $ {M_b}/{\rm GeV} $ $ {m_n}/{\rm GeV} $ $ {m_s}/{\rm GeV} $ $ a\left( {c\bar n} \right) $ $ a\left( {c\bar s} \right) $ $ a\left( {b\bar n} \right) $ $ a\left( {b\bar s} \right) $ Ref. [30] 1.44 4.48 0.23 0.328 0.223 0.249 0.275 0.313 EFG[31] 1.55 4.88 0.33 0.5 0.64/0.58 0.68/0.64 1.25/1.21 1.25/1.21 $ a $的单位为$ {\text{Ge}}{{\text{V}}^2} $,RMS误差为$ {\chi _{{\text{RMS}}}} = {\text{0}}{\text{.001 GeV}} $。 由表3中列出的参数,代入方程(1)求得单重介子D(
$ {D_s} $ )和B($ {B_s} $ )的P波自旋平均质量。例如,粲介子D和$ {D_s} $ 的自旋平均质量$ \bar M\left( {{D \mathord{\left/ {\vphantom {D {{D_S}}}} \right. } {{D_S}}}} \right) $ 的计算结果为其中
$v_c^2 = 1 - {{m_{{\text{bare}}c}^2} \mathord{\big/ {\vphantom {{m_{{\text{bare}}c}^2} {M_c^2}}} } {M_c^2}} = 1 - \frac{{{{1.275}^2}}}{{{{1.44}^2}}} = 0.216$ 。类似地,我们可以计算出底介子B和$ {B_s} $ 的自旋平均质量($ v_b^2 = 1 - $ ${\left( {{{4.18} \mathord{\left/ {\vphantom {{4.18} {4.48}}} \right. } {4.48}}} \right)^2} = 0.129~4$ )。所有单重介子自旋平均质量的计算结果已在表4中列出。
GeV Means $ {\left( D \right)_{\rm 1S}} $ $ {\left( D \right)_{\rm 1P}} $ $ {\left( {{D_s}} \right)_{\rm 1S}} $ $ {\left( {{D_s}} \right)_{\rm 1P}} $ $ {\left( B \right)_{\rm 1S}} $ $ {\left( B \right)_{\rm 1P}} $ $ {\left( {{B_s}} \right)_{\rm 1S}} $ $ {\left( {{B_s}} \right)_{\rm 1P}} $ Exp. 1.975 2.429 2.076 2.512 5.313 5.733 5.403 5.836 Regge 1.981 2.437 2.079 2.532 5.290 5.713 5.388 5.825 由表1和表2中的实验值所计算的S波自旋平均质量(GeV)亦列在表中。 -
由于存在色库仑项
$ V\left( r \right) $ ,等式(7)的谐振子近似仅适用于$ V\left( r \right) $ 较小的情况,即r很大($r \gg {1 \mathord{\big/ {\vphantom {1 {\sqrt b }}} } {\sqrt b }}$ )时。本文将$ V\left( r \right) $ 看成一阶微扰来计算自旋平均质量$ \bar M $ 。对短程夸克间极限下($ r \to 0 $ )的径向波函数我们单独求解,这样就可以用它来计算由短程势$ V\left( r \right) $ 所引起的自旋相关质量$ \Delta M $ 的贡献,具体结果将在第五部分给出。在长程范围内,由谐振子波函数
${\psi ^N} = R_{{n_r}L}^N\left( {\alpha ,\,r} \right)$ ${Y_{Lm}}\left( {\theta ,\,\phi } \right)$ 可给出单重介子的波函数$ {\Psi _M} $ 。其中,量子数$ N \equiv 2{n_r} + L $ ,$\alpha = \sqrt {\mu \omega } = {\left( {{{{b^2}\mu } \mathord{/ {\vphantom {{{b^2}\mu } \nu }} } \nu }} \right)^{{1 \mathord{/ {\vphantom {1 4}} } 4}}}$ 是逆谐振子长度,$ {Y_{Lm}} $ 为球谐函数,$ R_{{n_r}L}^N\left( {\alpha ,r} \right){\sim}{\left( {\alpha r} \right)^L}L_{{n_r}}^{L + {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}\left[ {{{\left( {\alpha r} \right)}^2}} \right] $ $ {{\rm e}^{{{ - {{\left( {\alpha r} \right)}^2}} \mathord{\left/ {\vphantom {{ - {{\left( {\alpha r} \right)}^2}} 2}} \right. } 2}}} $ 为谐振子的径向波函数,$ L_{{n_r}}^{L + {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}} $ 为对应的拉盖尔多项式。对于
$ {n_r} = 0 $ 的纯轨道态,可以用近似等式${\mu _q} \simeq \sqrt {m_q^2 + {{\left\langle {{p^2}} \right\rangle }_N}} $ 和$\nu \simeq b\sqrt {{{\left\langle {{r^2}} \right\rangle }_N}} $ 来将等式(6)写为其中用到了谐振子态
$ {\psi ^N} $ 下的$ {p^2} $ 和$ {r^2} $ 平均值。由此可以写出辅助场变量$\left\{ {\mu ,\,\nu } \right\}$ 所满足的两个非线性方程:借助数值迭代方法(附录),我们可以在长程和短程区域内分别求解方程(10)。选取输入值
$ {m_q} $ ,$ b $ 和$ {\alpha _s} $ 如表5中所列,其中的夸克质量源自文献[30],则可解出辅助场变量$\left\{ {\mu ,\,\nu } \right\}$ 和${\mu _q} = {\mu \mathord{\left/ {\vphantom {\mu {\left( {1 - {\mu \mathord{\left/ {\vphantom {\mu {{m^{}_Q}}}} \right. } {{m^{}_Q}}}} \right)}}} \right. } {\left( {1 - {\mu \mathord{\left/ {\vphantom {\mu {{m^{}_Q}}}} \right. } {{m^{}_Q}}}} \right)}}$ 的解如表6所列,其中还列出了计算自旋耦合参数所需的其他参数$ {\varepsilon _{V,S}} $ 。在图1中,我们展示了单重介子系统(P波)径向波函数,即哈密顿量(7)的本征函数数值解,其中输入值如表5所列。State $ {m_q} $ b/GeV2 $ {\alpha _s} $ $ {C_0} $ ${\varepsilon ^{}_V}$ ${\varepsilon ^{}_S}$ $ D $ 0.230 0.219 0.761 −0.590 −0.70 0.62 $ {D_s} $ 0.328 0.246 0.680 −0.660 −2.30 2.10 $ B $ 0.230 0.280 0.680 −0.428 −0.65 0.56 $ {B_s} $ 0.328 0.317 0.730 −0.450 −0.75 0.82 重夸克质量$ {m_c} = 1440\;{\text{MeV}} $,$ {m_b} = 4480\;{\rm{MeV}} $,轻夸克质量$ {m_q} $取自文献[30]。禁闭势参数$ b $在相应介子的Regge轨迹张力$ a $附近取值。 State $ {\mu _d} $ $ \mu $ $ \nu $ $ {\mu _{dH}} $ $ {\mu _H} $ ${a_B}{\text{/} }{\rm GeV}^{-1}$ $ {a_1} $ $ {a_2} $ $ b $ $ c $ $ D $ 684 464 809 234 201 11.84 34.46 32.08 26.78 10.48 Ds 760 497 854 332 270 9.95 85.23 35.82 27.36 11.52 D 835 704 830 234 222 12.00 28.52 15.49 11.96 6.99 Ds 920 763 877 334 311 8.09 14.14 13.53 10.21 6.65 由式(20)计算的自旋耦合参数$ {a_1} $,$ {a_2} $,$ b $和$ c $,由式(10)和P波其他相关参数求得的辅助场($ \mu $,$ \nu $)。除了$ {a_B} $之外单位都为MeV。 在短程范围,
$ V $ 相对于谐振子势中占主导地位,我们可以通过忽略后者来求解哈密顿量(7)的本征方程。小r范围内径向波函数[记为$ R_{{n_r}L}^H\left( r \right) $ ]的结果是知道的,即类氢原子解,其中有效玻尔半径${a_B} = {1 \mathord{\left/ {\vphantom {1 {\left( {{\mu ^{}_H}{k_s}} \right)}}} \right. } {\left( {{\mu ^{}_H}{k_s}} \right)}}$ ,约化质量${\mu ^{}_H} = {{{\mu ^{}_{qH}}} \mathord{\left/ {\vphantom {{{\mu ^{}_{qH}}} {\left( {1 + {{{\mu ^{}_{qH}}} \mathord{\left/ {\vphantom {{{\mu ^{}_{qH}}} {{m_Q}}}} \right. } {{m^{}_Q}}}} \right)}}} \right. } {\left( {1 + {{{\mu ^{}_{qH}}} \mathord{\left/ {\vphantom {{{\mu ^{}_{qH}}} {{m^{}_Q}}}} \right. } {{m^{}_Q}}}} \right)}}$ ,$ {\mu _{qH}} = {\left\langle {\sqrt {m_q^2 + {p^2}} } \right\rangle _H} $ 用到了类氢波函数${\psi ^H} = R_{{n_r}L}^H\left( r \right){Y_{Lm}}\left( {\theta ,\,\phi } \right)$ 的平均值。$ {n_r} = 0 $ 时,在态$ {\psi ^H} $ 中的动量平均值为$ {\left\langle {{p^2}} \right\rangle _H} = {{{\mu _H}{k_s}} \mathord{\left/ {\vphantom {{{\mu _H}{k_s}} {\left[ {{a_B}{{\left( {L + 1} \right)}^2}} \right]}}} \right. } {\left[ {{a_B}{{\left( {L + 1} \right)}^2}} \right]}} = $ $ {\left[ {{{{\mu _H}{k_s}} \mathord{\left/ {\vphantom {{{\mu _H}{k_s}} {\left( {L + 1} \right)}}} \right. } {\left( {L + 1} \right)}}} \right]^2} $ ,而$ {\mu _{qH}} $ 近似等于$ \sqrt {m_q^2 + {{\left\langle {{p^2}} \right\rangle }_H}} $ ,结果发现$ {\mu _{qH}} $ 满足自洽方程,对这个方程,我们可以用迭代法求解,
$ {\mu _{qH}} $ 和$ {\mu _H} $ 的数值结果列在表6中。于是,单重介子的自旋无关质量可写为
其中
在长程(N标记)和短程(H标记)范围,单粒子动能
$ E_i^{N,H} = {\left[ {m_i^2 + {{\left\langle {{p^2}} \right\rangle }_{N,H}}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} $ ,可以分别写成参数
$\left\{ {\mu ,\,\nu } \right\}$ 以及其他参数都在表5和表6中列出。我们用式(12)和式(13)来计算P波单重介子D,$ {D_s} $ ,$ B $ 和$ {B_s} $ 的自旋无关质量,其结果如表4所列。表4中还列出了介子自旋平均质量的实验值和由Regge轨迹所预测的结果,以供比较。 -
与类氢原子类似,重介子中的夸克-反夸克之间也存在色自旋轨道势和张量势,这是针对重夸克系统通过洛伦兹不变的树图级夸克-反夸克散射振幅推导的[34]。当应用于单重介子系统
$ \bar qQ $ 时,它的形式是[9, 24, 34]其中
$ V $ 和$ S $ 分别是重夸克和反夸克之间的矢量势和标量势,$ V' $ ,$ S' $ 和$ V'' $ 是它们的导数,由于轻的反夸克(
${m_q} \ll {\varLambda _{\rm QCD}}$ )的运动是相对性的,因此它的瞬时势(15)必须包含对应的相对论效应。注意,在手征极限$ {m_q} \to 0 $ 下,磁矩项${\sim}{{{{\boldsymbol{S}}_q}} \mathord{\left/ {\vphantom {{{{\boldsymbol{S}}_q}} {{m_q}}}} \right. } {{m_q}}}$ 的存在导致自旋相互作用(15)是发散的。一种简单的发散处理方法是进行以下替换[9]来调节瞬时势
$ {V^{{\text{quasi - static}}}} $ ,其中$ {\varepsilon _{V,S}} $ 是待定的小参数,动量为${\boldsymbol P}\left( {{{\boldsymbol P}_q} = - {{\boldsymbol P}_Q}} \right)$ 的夸克$i\left( {i = \bar q,\,Q} \right)$ 的单粒子动能为${E_i} = {\left[ {m_i^2 + {{\left| {\boldsymbol{P}} \right|}^2}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}$ 。注意,比值${{{m_i}} \mathord{\left/ {\vphantom {{{m_i}} {{E_i}}}} \right. } {{E_i}}} = {{{m_i}} \mathord{\left/ {\vphantom {{{m_i}} {\sqrt {m_i^2 + {{\left| {\boldsymbol{P}} \right|}^2}} }}} \right. } {\sqrt {m_i^2 + {{\left| {\boldsymbol{P}} \right|}^2}} }} = $ $ \sqrt {1 - v_i^2}$ 是一个相对论因子。当夸克$ i $ 作相对论性运动时,相对论因子$ {{{m_i}} \mathord{\left/ {\vphantom {{{m_i}} {{E_i}}}} \right. } {{E_i}}} $ 随着相对论长度收缩$r \to r\sqrt {1 - {v^2}} = $ $ r\left( {{m \mathord{\left/ {\vphantom {m E}} \right. } E}} \right) $ 而变大。显然,相对论因子在夸克的低动量极限($ {p \mathord{\left/ {\vphantom {p {{m_i}}}} \right. } {{m_i}}} \to 0 $ )下是回归到1。式(17)的相对论处理通过用$ {1 \mathord{\left/ {\vphantom {1 {{E_i}}}} \right. } {{E_i}}} = {\left( {m_i^2 + {p^2}} \right)^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} $ 替换$ {1 \mathord{\left/ {\vphantom {1 {{m_i}}}} \right. } {{m_i}}} $ 来消除式(15)的发散($ {{{\sim}1} \mathord{\left/ {\vphantom {{{\sim}1} {{m_i}}}} \right. } {{m_i}}} $ ),这确保了与轻夸克$ q $ 相关的自旋相互作用不会像$ {1 \mathord{\left/ {\vphantom {1 {{m_i}}}} \right. } {{m_i}}} $ 一样在手征极限下发散。相反,它们应该具有固定的有限大小极限$ \left\langle {{p^{ - 1}}} \right\rangle $ (和口袋模型一样),这个极限又由紧闭尺度$ \varLambda _{\rm QCD}^{ - 1} $ 所控制。通过将式(17)正则化,自旋相互作用(15)变为
其中四个自旋耦合参数为
这里的量子平均值是在
$ \bar qQ $ 系统的径向波函数$ {\Psi _M} $ 上平均得到的,而类似的关系也适用于
$ \tilde V'' $ 和$ {\nabla ^2}\tilde V $ 。当$ \Delta M = \left\langle {\Delta {H^{{\text{quasi - static}}}}} \right\rangle $ 已知时,即可计算单重介子的质量$ M = \bar M + $ $ \Delta M $ 。我们在表7中列出了粲介子和底介子P波的质量的计算结果。
Mass Splitting of Excited Charmed and Bottom Mesons in 1P Wave
doi: 10.11804/NuclPhysRev.39.2022001
- Received Date: 2022-01-02
- Rev Recd Date: 2022-03-03
- Publish Date: 2022-03-01
-
Key words:
- regge trajectory /
- phenomenological model /
- heavy-light mesons
Abstract: We apply linear Regge trajectory and relatived quark model to revisit the low-lying excited singly-heavy mesons and explain the narrow mass splitting of these excited states to arise mainly from relativistic effect of the light quarks. The P-wave masses of the singly heavy mesons (D,
Citation: | Qing MEI, Yongxin SONG, Duojie JIA. Mass Splitting of Excited Charmed and Bottom Mesons in 1P Wave[J]. Nuclear Physics Review, 2022, 39(1): 30-36. doi: 10.11804/NuclPhysRev.39.2022001 |