-
因与氢原子非常相似,单重介子系统的研究一直受人们关注。二十多年来,重介子家族不断扩大,这得益于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]。
表 1 粲介子和奇异粲介子的实验质量[1]
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的误差未列出。 表 2 底介子和奇异底介子的实验质量[1]
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]$$ {\left( {{{\bar M}_L} - {M_Q}} \right)^2} = \pi aL + {\left( {{m_q} + {M_Q} - \frac{{m_{{\text{bare}}Q}^2}}{{{M_Q}}}} \right)^2} \text{,} $$ (1) 其中: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]预言的基态质量方程:$$ \bar M\left( {1S} \right) = {M_Q} + {m_q} + \frac{{P_Q^2}}{{{M_Q}}} 。 $$ (2) 注意,重夸克极限
$ {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} $ 介子的弦张力系数。表 3 通过Regge轨迹方程(1)和相对论夸克模型确定的夸克有效质量
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) $ 的计算结果为$$ \begin{aligned}{\bar M_{L = 1}}\left( D \right) =& 1.44 + {\left[ {0.223\pi + {{\left( {0.230 + 1.44v_c^2} \right)}^2}} \right]^{{1 {\left/ {\vphantom {1 2}} \right. } 2}}}\\ =& 2.436~6\;{\text{GeV}} \text{,}\end{aligned} $$ $$\begin{split} {\bar M_{L = 1}}\left( {{D_s}} \right) =& 1.44 + {\left[ {0.249\pi + {{\left( {0.328 + 1.44v_c^2} \right)}^2}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} \\=& 2.531~2\;{\text{GeV}} \text{,} \end{split}$$ 其中
$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$ )。$$ \begin{split}{\bar M_{L = 1}}\left( B \right) =& 4.48 + {\left[ {0.275\pi + {{\left( {0.230 + 4.48v_b^2} \right)}^2}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} \\=& 5.712~7\;{\text{GeV}} \text{,} \end{split}$$ $$\begin{split} {\bar M_{L = 1}}\left( {{B_s}} \right) =& 4.48 + {\left[ {0.313\pi + {{\left( {0.328+ 4.48v_b^2} \right)}^2}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} \\=& 5.824~3\;{\text{GeV}} 。 \end{split}$$ 所有单重介子自旋平均质量的计算结果已在表4中列出。
表 4 由Regge关系计算所得的P波单重介子自旋平均质量
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)亦列在表中。 -
为了计算P波单重介子的质量并求出相应的波函数
$ {\Psi _M} $ ,我们考虑文献[9]中提出的相对论夸克模型,其中重夸克$ Q $ 和轻反夸克$ \bar q $ 之间的相互作用$ V + S $ 为色库仑势$ \left( V \right) $ 加上线性势$ \left( S \right) $ ,$$ {H^{SI}} = {\left[ {m_Q^2 + {{\left| {\boldsymbol{P}} \right|}^2}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} + {\left[ {m_q^2 + {{\left| {\boldsymbol{P}} \right|}^2}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} + V\left( r \right) + S\left( r \right) \text{,} $$ (3) $$ V\left( r \right) = - \frac{{4{\alpha _s}}}{3}\frac{{F\left( r \right)}}{r} \text{,} \quad S\left( r \right) = br + {C_0} \text{,} $$ (4) 其中
${\boldsymbol{P}}\left( { - {\boldsymbol{P}}} \right)$ 为夸克$ \bar q\left( Q \right) $ 的三动量,$ {\alpha _s} $ 为强耦合参数,且$$ F\left( r \right) \equiv 0.25{\rm Erf}\left( r \right) + 0.15{\rm Erf}\left( {{r \mathord{\left/ {\vphantom {r {10}}} \right. } {10}}} \right) + 0.2{\rm Erf}\left( {{r \mathord{\left/ {\vphantom {r {10}}} \right. } {100~0}}} \right) \text{,} $$ 其中函数
${\rm Erf}\left( x \right) = \left( {{2 \mathord{\big/ {\vphantom {2 {\sqrt \pi }}} } {\sqrt \pi }}} \right)\int_0^x {{\rm d}t{{\rm e}^{ - {t^2}}}}$ 是描述短程范围内轻反夸克$ \bar q $ 非点状性质的形状因子[9]。这里的$ b $ 是与Regge轨迹斜率$ a $ 相关的禁闭参数,$ {C_0} $ 是常数。如表1和表2所示,文献[9]预言低能质量与实验值相当一致,但对激发态质量的预言,比实验数据偏大。我们借助辅助场(AF)方法[32, 33]来求解本征方程
$ {H^{SI}}{\Psi _M} = \bar M{\Psi _M} $ ,其中$ \bar M $ 是我们所考虑的介子质量。AF方法的基本思想是通过引入一个辅助场$ \lambda $ 将非线性算符$ \sqrt O $ 重新表述为$ \sqrt O = {\text{mi}}{{\text{n}}_{\lambda > 0}}\left[ {{O \mathord{\left/ {\vphantom {O {\left( {2\lambda } \right) + {\lambda \mathord{\left/ {\vphantom {\lambda 2}} \right. } 2}}}} \right. } {\left( {2\lambda } \right) + {\lambda \mathord{\left/ {\vphantom {\lambda 2}} \right. } 2}}}} \right] $ ,其中最小化仅当$ \lambda = \sqrt O $ 时实现。假设$ \lambda $ 是一个慢变量,它满足一个半经典的、类Hartree方程$ \lambda = \left\langle {\sqrt O } \right\rangle $ 。据此,我们为第二个动量项引入一个辅助场$ {\mu _q} $ ,为S中的线性势$ br = \sqrt {{{\left( {br} \right)}^2}} $ 引入另一个辅助场变量$ \nu $ ,从而将哈密顿量(3)重写为$$ \begin{split} {H^{SI}} =& {m_Q^{}} + \frac{{{p^2}}}{{2{m_Q}}} + \frac{{m_q^2 + {p^2}}}{{2{\mu _q^{}}}} + \frac{{{{\left( {br} \right)}^2}}}{{2\nu }} +\\& \frac{{{\mu _q} + \nu }}{2} - \frac{{4{\alpha _s}}}{3}\frac{{F\left( r \right)}}{r} + {C_0} , \end{split} $$ (5) 其中两个辅助场满足方程
$$ {\mu _q} = \left\langle {\sqrt {m_q^2 + {p^2}} } \right\rangle \text{,} \quad \nu = b\left\langle r \right\rangle \text{,} $$ (6) 这里
$ {p^2} = {\left| {\boldsymbol{P}} \right|^2} $ ,且对第一个动量项做了重夸克($ {m_Q} \gg \left| {\boldsymbol{P}} \right| $ ) 展开。 在某种意义上,AF方法用类似于多体系统的平均场方法,将哈密顿量(3)线性化。显然,引入介子系统的有效约化质量$\mu \equiv {{{\mu _q}} \mathord{\left/ {\vphantom {{{\mu _q}} {\left( {1 + {{{\mu _q}} \mathord{\left/ {\vphantom {{{\mu _q}} {{m_Q}}}} \right. } {{m_Q}}}} \right)}}} \right. } {\left( {1 + {{{\mu _q}} \mathord{/ {\vphantom {{{\mu _q}} {{m_Q}}}} } {{m_Q}}}} \right)}}$ 后,方程(5)就约化为$$ \begin{split} {H^{{SI}}} =& \frac{{{p^2}}}{{2\mu }} + \frac{1}{2}\mu {\omega ^2}{r^2}{\text{ }} - \frac{{4{\alpha _s}}}{3}\frac{{F\left( r \right)}}{r} +\\& {m_Q} + \frac{{m_q^2}}{{2{\mu _q}}} + \frac{{{\mu _q} + \nu }}{2} + {C_0} 。 \end{split} $$ (7) 若忽略色库仑相互作用,它就是一个三维谐振子(HO)的非相对论模型,谐振子频率为
$\omega = {b \mathord{\big/ {\vphantom {b {\sqrt {\mu \nu } }}} } {\sqrt {\mu \nu } }}$ 。注意,虽然哈密顿方程(7)看起来像HO的哈密顿量,
$\bar M - {m_Q} \sim \omega \left( {L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}} \right)$ ,但其能级不是均匀分布的,这是因为谐振子频率$\omega \sim {1 \mathord{\big/ {\vphantom {1 {\sqrt {\mu \nu } }}} } {\sqrt {\mu \nu } }}$ 与$ \sqrt {b\left( {L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}} \right)} $ 成反比。这样,就与方程(1)所示的Regge关系$\bar M - {m_Q} \sim$ $\sqrt {L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}} $ 相一致。 -
由于存在色库仑项
$ 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)写为$$ {\mu _q} = {\left[ {m_q^2 + {\alpha ^2}\left( {L + \frac{3}{2}} \right)} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} = {\left[ {m_q^2 + \sqrt {\frac{\mu }{\nu }b} \left( {L + \frac{3}{2}} \right)} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} \text{,} $$ (8) $$ \nu = b{\left[ {\frac{{L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}}}{{\mu \omega }}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} = {\left[ {b\left( {L + \frac{3}{2}} \right)\sqrt {\frac{\nu }{\mu }} } \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} 。 $$ (9) 其中用到了谐振子态
$ {\psi ^N} $ 下的$ {p^2} $ 和$ {r^2} $ 平均值。由此可以写出辅助场变量$\left\{ {\mu ,\,\nu } \right\}$ 所满足的两个非线性方程:$$\begin{split}& \mu \equiv \frac{{{{\left[ {m_q^2 + b\sqrt {{\mu \mathord{\left/ {\vphantom {\mu \nu }} \right. } \nu }} \left( {L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}} \right)} \right]}^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}}}{{{{1 + {{\left[ {m_q^2 + b\sqrt {{\mu \mathord{\left/ {\vphantom {\mu \nu }} \right. } \nu }} \left( {L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}} \right)} \right]}^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}} \mathord{\left/ {\vphantom {{1 + {{\left[ {m_q^2 + b\sqrt {{\mu \mathord{\left/ {\vphantom {\mu \nu }} \right. } \nu }} \left( {L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}} \right)} \right]}^{{1 {/ {\vphantom {1 2}} } 2}}}} {{m^{}_Q}}}} \right. } {{m^{}_Q}}}}} \text{,} \\[1mm]& \nu = \sqrt {b\left( {L + {3 \mathord{\left/ {\vphantom {3 2}} \right. } 2}} \right)} {\left( {{\nu \mathord{\left/ {\vphantom {\nu \mu }} \right. } \mu }} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. } 4}}} 。\end{split} $$ (10) 借助数值迭代方法(附录),我们可以在长程和短程区域内分别求解方程(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所列。表 5 P波单重介子哈密顿量(7)的模型参数的取值
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 $附近取值。 表 6 自旋耦合参数
$ {a_1} $ ,$ {a_2} $ ,$ b $ 和$ c $ 及辅助场($ \mu $ ,$ \nu $ )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 ^{}_{qH}}{k_s}} \mathord{\left/ {\vphantom {{{\mu _{qH}}{k_s}} {\left[ {{a^{}_B}} \right]}}} \right. } {\left[ {{a^{}_B}} \right]}} = {\left[ {m_q^2 + \frac{1}{{{{\left( {L + 1} \right)}^2}}}\frac{{{{\left( {{\mu ^{}_{qH}}{k_s}} \right)}^2}}}{{{{\left( {1 + {{{\mu ^{}_{qH}}} \mathord{\left/ {\vphantom {{{\mu ^{}_{qH}}} {{m^{}_Q}}}} \right. } {{m^{}_Q}}}} \right)}^2}}}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} 。 $$ (11) 对这个方程,我们可以用迭代法求解,
$ {\mu _{qH}} $ 和$ {\mu _H} $ 的数值结果列在表6中。于是,单重介子的自旋无关质量可写为
$$\begin{split} \bar M\left( {qQ} \right) =& \omega \left( {L + \frac{3}{2}} \right) - \frac{{4{\alpha _s}}}{3}{\left\langle {\frac{{F\left( r \right)}}{r}} \right\rangle _N} + {m_Q} + \frac{{m_q^2}}{{2{\mu _q}}} + \frac{{{\mu _q} + \nu }}{2} + {C_0} \\=& {m_Q} + \frac{b}{{\sqrt {\mu \nu } }}\left( {L + \frac{3}{2}} \right) - \frac{4}{3}{\alpha _s}{F_L} + E\left( {{\mu _q},\,\nu } \right) \text{,}\\[-15pt] \end{split}$$ (12) 其中
$$\begin{split} {F_L} \equiv & \frac{{4{\alpha _s}}}{3}{\left\langle {\frac{{F\left( r \right)}}{r}} \right\rangle _N} = \int {{\rm d}r{r^2}} {\left| {R_{0L}^H\left( r \right)} \right|^2}\left( {\frac{{F\left( r \right)}}{r}} \right) \text{,} \\& E\left( {{\mu _q},\nu } \right) \equiv \frac{1}{2}\left( {{\mu _q} + \nu } \right) + \frac{{m_q^2}}{{2{\mu _q}}} + {C_0} 。 \end{split}$$ (13) 在长程(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}}} $ ,可以分别写成$$ E_q^H = {\mu _{qH}} \text{,} \quad E_Q^H = {\left[ {m_Q^2 + \frac{{{{\left( {{\mu _H}{k_s}} \right)}^2}}}{{{{\left( {L + 1} \right)}^2}}}} \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} \text{,} $$ $$ E_{q,\,Q}^N = {\left[ {m_{q,\,Q}^2 + a\left( {L + \frac{3}{2}} \right)\sqrt {{\mu \mathord{\left/ {\vphantom {\mu \nu }} \right. } \nu }} } \right]^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}} 。 $$ (14) 参数
$\left\{ {\mu ,\,\nu } \right\}$ 以及其他参数都在表5和表6中列出。我们用式(12)和式(13)来计算P波单重介子D,$ {D_s} $ ,$ B $ 和$ {B_s} $ 的自旋无关质量,其结果如表4所列。表4中还列出了介子自旋平均质量的实验值和由Regge轨迹所预测的结果,以供比较。 -
与类氢原子类似,重介子中的夸克-反夸克之间也存在色自旋轨道势和张量势,这是针对重夸克系统通过洛伦兹不变的树图级夸克-反夸克散射振幅推导的[34]。当应用于单重介子系统
$ \bar qQ $ 时,它的形式是[9, 24, 34]$$ \begin{aligned} \Delta {H^{{\text{quasi - static}}}} =& \left( {\frac{{V' - S'}}{r}} \right){\boldsymbol L} \boldsymbol\cdot \left( {\frac{{{{\boldsymbol S}_q}}}{{2m_q^2}} + \frac{{{{\boldsymbol S}_Q}}}{{2m_Q^2}}} \right) + \\&\left( {\frac{{V'}}{r}} \right){\boldsymbol L} \boldsymbol\cdot \left( {\frac{{{{\boldsymbol S}_q} + {{\boldsymbol S}_Q}}}{{{m_q} + {m_Q}}}} \right) + \frac{1}{{3{m_q}{m_Q}}}\left( {\frac{{V'}}{r} - V''} \right){S_{12}} +\\& \frac{2}{{3{m_q}{m_Q}}}\left( {{\nabla ^2}V} \right){{\boldsymbol S}_q} \boldsymbol\cdot {{\boldsymbol S}_Q} 。\\[-15pt] \end{aligned} $$ (15) 其中
$ V $ 和$ S $ 分别是重夸克和反夸克之间的矢量势和标量势,$ V' $ ,$ S' $ 和$ V'' $ 是它们的导数,$$\begin{split}& V' = \frac{{4{\alpha _s}}}{3}\left( {{F \mathord{\left/ {\vphantom {F {{r^2}}}} \right. } {{r^2}}} - {{F'} \mathord{\left/ {\vphantom {{F'} r}} \right. } r}} \right) \text{,} \\& V'' = \frac{{4{\alpha _s}}}{3}\left( {2{{F'} \mathord{\left/ {\vphantom {{F'} {{r^2}}}} \right. } {{r^2}}} - 2{F \mathord{\left/ {\vphantom {F {{r^3}}}} \right. } {{r^3}}} - {{F''} \mathord{\left/ {\vphantom {{F''} r}} \right. } r}} \right) \\& S' = a \text{,} \quad {\nabla ^2}V = \frac{{4{\alpha _s}}}{3}\left[ { - 4\pi {\delta ^3}\left( {\boldsymbol{r}} \right) + {{F''} \mathord{\left/ {\vphantom {{F''} r}} \right. } r}} \right] 。\end{split} $$ (16) 由于轻的反夸克(
${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]$$\begin{split}& V\left( r \right) \to \tilde V\left( r \right) = {\left( {\frac{{{m_q}{m_Q}}}{{{E_q}{E_Q}}}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2} + {{{\varepsilon _V}} \mathord{\left/ {\vphantom {{{\varepsilon _V}} 2}} \right. } 2}}}V\left( r \right){\left( {\frac{{{m_q}{m_Q}}}{{{E_q}{E_Q}}}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2} + {{{\varepsilon _V}} \mathord{\left/ {\vphantom {{{\varepsilon _V}} 2}} \right. } 2}}} , \\& S\left( r \right) \to \tilde S\left( r \right) = {\left( {\frac{{{m_q}{m_Q}}}{{{E_q}{E_Q}}}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2} + {{{\varepsilon _S}} \mathord{\left/ {\vphantom {{{\varepsilon _S}} 2}} \right. } 2}}}S\left( r \right){\left( {\frac{{{m_q}{m_Q}}}{{{E_q}{E_Q}}}} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2} + {{{\varepsilon _S}} \mathord{\left/ {\vphantom {{{\varepsilon _S}} 2}} \right. } 2}}} , \end{split}$$ (17) 来调节瞬时势
$ {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)变为
$$ \Delta {H^{\rm{quasi - static}}} = {a_1}{\boldsymbol L} \boldsymbol\cdot {{\boldsymbol S}_q} + {a_2}{\boldsymbol L} \boldsymbol\cdot {{\boldsymbol S}_Q} + b{{S}_{12}} + c{{\boldsymbol S}_q} \boldsymbol\cdot {{\boldsymbol S}_Q} \text{,} $$ (18) $$ {S_{12}} = 3{{{\boldsymbol S}}_q} \boldsymbol\cdot \hat r{{\boldsymbol S}_Q} \boldsymbol\cdot \hat r - {{\boldsymbol S}_q} \boldsymbol\cdot {{\boldsymbol S}_Q} \text{,} $$ (19) 其中四个自旋耦合参数为
$$ \begin{split}&{a_1} = \frac{1}{{2{m_q}}}\left\langle {\frac{{\tilde V' - \tilde S'}}{{{m_q}r}} + \frac{{2\tilde V'}}{{{m_Q}r}}} \right\rangle \text{,} \\& {a_2} = \frac{1}{{2{m_Q}}}\left\langle {\frac{{\tilde V' - \tilde S'}}{{{m_Q}r}} + \frac{{2\tilde V'}}{{{m_q}r}}} \right\rangle \text{,} \\& b = \frac{1}{{3{m_q}{m_Q}}}\left\langle {{{\tilde V'} \mathord{\left/ {\vphantom {{\tilde V'} r}} \right. } r} - \tilde V''} \right\rangle \text{,} \quad c = \left\langle {\frac{{2{\nabla ^2}\tilde V}}{{3{m_q}{m_Q}}}} \right\rangle 。\end{split} $$ (20) 这里的量子平均值是在
$ \bar qQ $ 系统的径向波函数$ {\Psi _M} $ 上平均得到的,而$$ \tilde V'\left( r \right) = V'\left( r \right){\left( {\frac{{{m_q}{m_Q}}}{{E_q^HE_Q^H}}} \right)^{1 + {\varepsilon ^{}_V}}} \text{,} \tilde S'\left( r \right) = S'\left( r \right){\left( {\frac{{{m_q}{m^{}_Q}}}{{E_q^NE_Q^N}}} \right)^{1 + {\varepsilon ^{}_S}}} \text{,} $$ 类似的关系也适用于
$ \tilde V'' $ 和$ {\nabla ^2}\tilde V $ 。当$ \Delta M = \left\langle {\Delta {H^{{\text{quasi - static}}}}} \right\rangle $ 已知时,即可计算单重介子的质量$ M = \bar M + $ $ \Delta M $ 。我们在表7中列出了粲介子和底介子P波的质量的计算结果。表 7 P波粲介子和底介子的计算结果,平均值
$ \bar M $ 如表4所列
-
摘要: 应用线性Regge轨迹和相对论夸克模型,研究了低激发态单重介子谱,并解释了这些激发态的窄质量劈裂源自轻夸克的相对论效应。计算了自旋宇称为
$ {J^ + } $ ($ J = 0 $ ,$ 1 $ ,$ 2 $ )的单重介子(D, Ds, B, Bs)P波的质量,并建议存在$ {0^ + } $ 态的,未被实验发现的$ B $ 介子和$ {B_s} $ 介子,其质量分别为5 659 MeV和5 788 MeV。计算表明,$ {D_{s{\text{0}}}}\left( {{\text{2317}}} \right) $ 包含70 MeV左右的质量下移而$ {D_{s{\text{1}}}}\left( {{\text{2460}}} \right) $ 则没有显著的质量移动。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,$ {D_s} $ , B,$ {B_s} $ ) with spin-parity$ {J^ + } $ ($ J = 0 $ ,$ 1 $ ,$ 2 $ ) have been computed, suggesting the unseen states of the mesons B and$ {B_s} $ to be the$ {0^ + } $ state with mass 5 659 and 5 788 MeV, respectively. Our calculations show that$ {D_{s{\text{0}}}}\left( {{\text{2\;317}}} \right) $ contains a mass downshift of around 70 MeV while$ {D_{s{\text{1}}}}\left( {{\text{2\;460}}} \right) $ has no significant mass shift.-
Key words:
- regge trajectory /
- phenomenological model /
- heavy-light mesons
-
表 1 粲介子和奇异粲介子的实验质量[1]
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的误差未列出。 表 2 底介子和奇异底介子的实验质量[1]
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的误差未列出。 表 3 通过Regge轨迹方程(1)和相对论夸克模型确定的夸克有效质量
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}} $。 表 4 由Regge关系计算所得的P波单重介子自旋平均质量
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)亦列在表中。 表 5 P波单重介子哈密顿量(7)的模型参数的取值
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 $附近取值。 表 6 自旋耦合参数
$ {a_1} $ ,$ {a_2} $ ,$ b $ 和$ c $ 及辅助场($ \mu $ ,$ \nu $ )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。 表 7 P波粲介子和底介子的计算结果,平均值
$ \bar M $ 如表4所列 -
[1] ZYLA P A, BARNETT R M, BERINGER J, et al. Prog Theor Exp Phys, 2020, 2020(8): 083C01. doi: 10.1093/ptep/ptaa104 [2] AAIJ R, ADEVA B, ADINOLFI M, et al. J High Energ Phys, 2013, 2013(9): 145. doi: 10.1007/JHEP09(2013)145 [3] AAIJ R, ADEVA B, ADINOLFI M, et al. Phys Rev D, 2016, 94(7): 072001. doi: 10.1103/PhysRevD.94.072001 [4] AAIJ R, ADEVA B, ADINOLFI M, et al. Phys Rev Lett, 2016, 117(15): 152003. doi: 10.1103/PhysRevLett.117.152003 [5] AABOUD M, AAD G, ABBOTT B, et al. Phys Rev Lett, 2018, 120(20): 202007. doi: 10.1103/PhysRevLett.120.202007 [6] AAIJ R, ABELLÁN BETETA C, ACKERNLEY T, et al. Phys Rev D, 2020, 102(11): 112003. doi: 10.1103/PhysRevD.102.112003 [7] MOIR G, PEARDON M, RYAN S M, et al. J High Energ Phys, 2013, 2013(5): 21. doi: 10.1007/JHEP05(2013)021 [8] LEWIS R, WOLOSHYN R M. Phys Rev D, 2000, 62(11): 114507. doi: 10.1103/PhysRevD.62.114507 [9] GODFREY S, ISGUR N. Phys Rev D, 1985, 32(1): 189. doi: 10.1103/PhysRevD.32.189 [10] ROSNER J L. Comments Nucl Part Phys, 1986, 16: 109. [11] GODFREY S, KOKOSKI R. Phys Rev D, 1991, 43: 1679. doi: 10.1103/PhysRevD.43.1679 [12] GODFREY S, MOATS K. Phys Rev D, 2016, 93(3): 034035. doi: 10.1103/PhysRevD.93.034035 [13] KE H W, ZHOU J H, LI X Q. Eur Phys J C, 2015, 75(1): 28. doi: 10.1140/epjc/s10052-015-3264-5 [14] SONG Q T, CHEN D Y, LIU X, et al. Eur Phys J C, 2015, 75(1): 30. doi: 10.1140/epjc/s10052-015-3265-4 [15] SONG Q T, CHEN D Y, LIU X, et al. Phys Rev D, 2015, 92(7): 074011. doi: 10.1103/PhysRevD.92.074011 [16] SUN Y, LIU X, MATSUKI T. Phys Rev D, 2013, 88(9): 094020. doi: 10.1103/PhysRevD.88.094020 [17] WANG Z G. Phys Rev D, 2011, 83(1): 014009. doi: 10.1103/PhysRevD.83.014009 [18] VAN BEVEREN E, RUPP G. Phys Rev D, 2011, 81(11): 118101. doi: 10.1103/PhysRevD.81.118101 [19] BADALIAN A M, BAKKER B L G. Phys Rev D, 2011, 84(3): 034006. doi: 10.1103/PhysRevD.84.034006 [20] CLOSE F E, THOMAS C E, LAKHINA O, et al. Phys Lett B, 2007, 647(2): 159. doi: 10.1016/j.physletb.2007.01.052 [21] GODFREY S, JARDINE I T. Phys Rev D, 2014, 89(7): 074023. doi: 10.1103/PhysRevD.89.074023 [22] OLSEN S L, SKWARNICKI T, ZIEMINSKA D. Rev Mod Phys, 2018, 90(1): 015003. doi: 10.1103/RevModPhys.90.015003 [23] CHEN H X, CHEN W, LIU X, et al. Rept Prog Phys, 2017, 80(7): 076201. doi: 10.1088/1361-6633/aa6420 [24] CAHN R N, JACKSON J D. Phys Rev D, 2003, 68(3): 037502. doi: 10.1103/PhysRevD.68.037502 [25] JIA D J, DONG W C. Eur Phys J Plus, 2019, 134(3): 123. doi: 10.1140/epjp/i2019-12474-8 [26] CHEN B, WEI K W, ZHANG A L. Eur Phys J A, 2015, 51: 82. doi: 10.1140/epja/i2015-15082-3 [27] JIA D J, DONG W C, HOSAKA A. JPS Conf Proc, 2019, 26: 022021. doi: 10.7566/JPSCP.26.022021 [28] CHEN K, DONG Y, LIU X, et al. Eur Phys J C, 2018, 78(1): 20. doi: 10.1140/epjc/s10052-017-5512-3 [29] MANOHAR A V, WISE M B. Heavy Quark Physics[M]. Cambridge: Cambridge Univ Press, 2000: 191. [30] JIA D J, LIU W N, HOSAKA A. Phys Rev D, 2020, 101(3): 034016. doi: 10.1103/PhysRevD.101.034016 [31] EBERT D, FAUSTOV R N, GALKIN V O. Eur Phys J C, 2010, 66: 197. doi: 10.1140/epjc/s10052-010-1233-6 [32] SEMAY C, SILVESTRE-BRAC B, NARODETSKII I M. Phys Rev D, 2004, 69(1): 014003. doi: 10.1103/PhysRevD.69.014003 [33] SILVESTRE-BRAC B, SEMAY C, BUISSERET F. J Phys A, 2008, 41(27): 275301. doi: 10.1088/1751-8113/41/27/275301 [34] YNDURÁIN F J. The Theory of Quark and Gluon Interactions[M]. Berlin:Springer-Verlag, 1983: 241. [35] KALASHNIKOVA Y S, NEFEDIEV A V, SIMONOV Y A. Phys Rev D, 2001, 64(1): 014037. doi: 10.1103/PhysRevD.64.014037 [36] ALI KHAN A, BHATTACHARYA T, COLLINS S, et al. Phys Rev D, 2000, 62(5): 054505. doi: 10.1103/PhysRevD.62.054505 [37] KALASHNIKOVA Y S, NEFEDIEV A V. Phys Lett B, 2000, 492(1): 91. doi: 10.1016/S0370-2693(00)01063-7 [38] LACOURSE D, OLSSON M G. Phys Rev D, 1989, 39(9): 2751. doi: 10.1103/PhysRevD.39.2751