Masses of Scalar and Axial-Vector B Mesons Revisited

TL;DR

This paper investigates the near degeneracy of scalar charmed mesons and predicts the masses of related bottom mesons using heavy meson chiral perturbation theory, highlighting the importance of coupled channel effects and heavy quark symmetry.

Contribution

It provides a qualitative and quantitative analysis of meson mass degeneracy through self-energy effects and predicts bottom meson masses with corrections from heavy quark symmetry and QCD effects.

Findings

Near degeneracy explained by coupled channel effects.

Predicted B meson masses with specific values and corrections.

Mass differences larger than quark model predictions.

Abstract

The SU(3) quark model encounters a great challenge in describing even-parity mesons. Specifically, the $q\bar q$ quark model has difficulties in understanding the light scalar mesons below 1 GeV, scalar and axial-vector charmed mesons and $1^+$ charmonium-like state $X(3872)$. A common wisdom for the resolution of these difficulties lies on the coupled channel effects which will distort the quark model calculations. In this work, we focus on the near mass degeneracy of scalar charmed mesons, $D_{s0}^*$ and $D_0^{*0}$, and its implications. Within the framework of heavy meson chiral perturbation theory, we show that near degeneracy can be qualitatively understood as a consequence of self-energy effects due to strong coupled channels. Quantitatively, the closeness of $D_{s0}^*$ and $D_0^{*0}$ masses can be implemented by adjusting two relevant strong couplings and the renormalization scale appearing in the loop diagram. Then this in turn implies the mass similarity of $B_{s0}^*$ and $B_0^{*0}$ mesons. The $P_0^* P'_1$ interaction with the Goldstone boson is crucial for understanding the phenomenon of near degeneracy. Based on heavy quark symmetry in conjunction with corrections from QCD and $1/m_Q$ effects, we obtain the masses of $B^*_{(s)0}$ and $B'_{(s)1}$ mesons, for example, $M_{B_{s0}^*}= (5715\pm1)\,{\rm MeV}+\delta\Delta_S$, $M_{B'_{s1}}=(5763\pm1)\,{\rm MeV}+\delta\Delta_S$ with $\delta\Delta_S$ being $1/m_Q$ corrections. We find that the predicted mass difference of 48 MeV between $B'_{s1}$ and $B_{s0}^*$ is larger than that of $20\sim 30$ MeV inferred from the relativistic quark models, whereas the difference of 15 MeV between the central values of $M_{B'_{s1}}$ and $M_{B'_1}$ is much smaller than the quark model expectation of $60-100$ MeV.

Tables

Table 1. Table 1: Measured masses and widths of even-parity charmed mesons. The four p 𝑝 p -wave charmed meson states are denoted by D 0 ∗ , D 1 ′ , D 1 superscript subscript 𝐷 0 subscript superscript 𝐷 ′ 1 subscript 𝐷 1 D_{0}^{*},D^{\prime}_{1},D_{1} and D 2 ∗ superscript subscript 𝐷 2 D_{2}^{*} , respectively. In the heavy quark limit, D 1 ′ subscript superscript 𝐷 ′ 1 D^{\prime}_{1} has j = 1 / 2 𝑗 1 2 j=1/2 and D 1 subscript 𝐷 1 D_{1} has j = 3 / 2 𝑗 3 2 j=3/2 with j 𝑗 j being the total angular momentum of the light degrees of freedom. The data are taken from the Particle Data Group PDG . The last two columns are the predictions from the quark model Godfrey ; Di Pierro . “Large” means a broad width of order 100 MeV, while “small” implies a narrow width of order 10 MeV.

$J^P$	Meson	Mass (MeV)	$\Gamma$ (MeV)	Mass (MeV)	$\Gamma$
$0^{+}$	$D_0^{\ast }{\left(2400\right)}^0$	$2318\pm 29$	$267\pm 40$	2340$-$2410	large
	$D_0^{\ast }{\left(2400\right)}^{\pm }$	$2351\pm 7$	$230\pm 17$	2340$-$2410	large
$1^{\prime +}$	$D_1^{\prime }{\left(2430\right)}^0$	$2427\pm 26\pm 25$	${384}_{-75}^{+107}\pm 74$	2470$-$2530	large
$1^{+}$	$D_1{\left(2420\right)}^0$	$2420.8\pm 0.5$	$31.7\pm 2.5$	2417$-$2434	small
	$D_1{\left(2420\right)}^{\pm }$	$2432.2\pm 2.4$	$25\pm 6$	2417$-$2434	small
$2^{+}$	$D_2^{\ast }{\left(2460\right)}^0$	$2460.57\pm 0.15$	$47.7\pm 1.3$	2460$-$2467	small
	$D_2^{\ast }{\left(2460\right)}^{\pm }$	$2465.4\pm 1.3$	$46.7\pm 1.2$	2460$-$2467	small
$0^{+}$	$D_{s0}^{\ast }{\left(2317\right)}^{\pm }$	$2317.7\pm 0.6$		2400$-$2510	large
$1^{\prime +}$	$D_{s1}^{\prime }{\left(2460\right)}^{\pm }$	$2459.5\pm 0.6$		2528$-$2536	large
$1^{+}$	$D_{s1}{\left(2536\right)}^{\pm }$	$2535.10\pm 0.06$	$0.92\pm 0.05$	2543$-$2605	small
$2^{+}$	$D_{s2}^{\ast }{\left(2573\right)}^{\pm }$	$2569.1\pm 0.8$	$16.9\pm 0.8$	2569$-$2581	small

Table 2. Table 2: Mass shifts ( δ ​ M ≡ M − ℳ 𝛿 𝑀 𝑀 ℳ \delta M\equiv M-{\cal M} ) of heavy scalar mesons calculated in HMChPT. The renormalization scale is taken to be Λ = 1.3 Λ 1.3 \Lambda=1.3 (1.2) GeV for scalar D 𝐷 D ( B 𝐵 B ) mesons. Bare masses are taken from Godfrey . All masses and widths are given in MeV and only the central values are listed here.

Meson	Bare mass $ℳ$	$\delta M$	$M$
$D_{s0}^{\ast }$	2480	$-162$	2318
$D_0^{\ast }$	2400	$-79$	2321
$B_{s0}^{\ast }$	5831	$-136$	5694
$B_0^{\ast }$	5756	$-55$	5701
$D_{s1}^{\prime }$	2550	$-98$	2452
$D_1^{\prime }$	2460	$-41$	2419
$B_{s1}^{\prime }$	5857	$-85$	5772
$B_1^{\prime }$	5777	$-34$	5743

Table 3. Table 3: Same as Table 2 except that the propagator i / ( 2 ​ v ⋅ k ) 𝑖 ⋅ 2 𝑣 𝑘 i/(2v\cdot k) is used for all heavy meson states and v ⋅ k ⋅ 𝑣 𝑘 v\cdot k is set to be zero inside the loop integral calculation.

Meson	Bare mass $ℳ$	$\delta M$	$M$
$D_{s0}^{\ast }$	2480	$-117$	2363
$D_0^{\ast }$	2400	$-44$	2356
$B_{s0}^{\ast }$	5831	$-63$	5768
$B_0^{\ast }$	5756	$-27$	5729
$D_{s1}^{\prime }$	2550	$-48$	2502
$D_1^{\prime }$	2460	$-22$	2438
$B_{s1}^{\prime }$	5857	$-40$	5817
$B_1^{\prime }$	5777	$-9$	5768

Table 4. Table 4: Predicted masses (in MeV) of B ( s ) ​ 0 ∗ subscript superscript 𝐵 𝑠 0 B^{*}_{(s)0} and B ( s ) ​ 1 ′ subscript superscript 𝐵 ′ 𝑠 1 B^{\prime}_{(s)1} mesons in the literature.

		$B_0^{\ast }$	$B_1^{\prime }$	$B_{s0}^{\ast }$	$B_{s1}^{\prime }$
Relativistic quark model
Di Pierro et al.	Di Pierro	5706	5742	5804	5842
Godfrey et al. (2016)	Godfrey:2016	5720	5738	5805	5822
Lakhina et al.	Swanson	5730	5752	5776	5803
Liu et al.	Liu	5749	5782	5815	5843
Ebert et al.	Ebert	5749	5774	5833	5865
Sun et al.	Zhu	5756	5779	5830	5858
Godfrey et al. (1991)	Godfrey	5756	5777	5831	5857
Nonrelativistic quark model
Lu et al.	Lu	5683	5729	5756	5801
Ortega et al.	Ortega			5741	5858
HMChPT
Bardeen et al.	Bardeen	$5627\pm 35$	$5674\pm 35$	$5718\pm 35$	$5765\pm 35$
Colangelo et al.	Colangelo	$5708.2\pm 22.5$	$5753.3\pm 31.1$	$5706.6\pm 1.2$	$5765.6\pm 1.2$
Cheng et al.	Cheng:2014	$5715\pm 22+\delta {\Delta }_S$	$5752\pm 31+\delta {\Delta }_S$	$5715\pm 1+\delta {\Delta }_S$	$5763\pm 1+\delta {\Delta }_S$
Unitarized ChPT
Torres-Rincon et al.	Rincon	5530	5579	5748	5799
Guo et al.	Guo:dynamic ; Guo:Ds1	Two poles		$5725\pm 39$	$5778\pm 7$
Cleven et al.	Cleven:2014			$5625\pm 45$	$5671\pm 45$
Altenbuchinger et al.	Altenbuchinger			$5726\pm 28$	$5778\pm 26$
Cleven et al.	Cleven:2010aw			$5696\pm 40$	$5742\pm 40$
Lattice QCD
Gregory et al. (HPQCD)	Gregory			$5752\pm 30$	$5806\pm 30$
Koponen et al. (UKQCD)	Koponen			$5760\pm 9$	$5807\pm 9$
Lang et al.	Lang:Bs			$5713\pm 22$	$5750\pm 26$
Chiral loop corrections
I. W. Lee et al.	Lee	5637	5673	5634	5672
Albaladejo et al.	Albaladejo			$5709\pm 8$	$5755\pm 8$
QCD sum rules
Z. G. Wang	ZGWang	$5720\pm 50$	$5740\pm 50$	$5700\pm 60$	$5760\pm 60$
Others
Kolomeitsev et al.	Lutz:2003fm	5526	5590	5643	5690
Matsuki et al.	Matsuki	5592	5649	5617	5682
Orsland et al.	Orsland	5592	5671	5667	5737
Vijande et al.	Vijande	5615		5679	5713
Badalian et al.	Badalian	$5675\pm 20$	$5725\pm 20$	$5710\pm 15$	$5730\pm 15$
Lähde et al.	Riska	5678	5686	5781	5795
Dmitrašinović	Dmitrasinovic:2012zz	$5718\pm 25$	$5732\pm 25$	$5719\pm 25$	$5765\pm 25$

Equations

$$|X(3872)⟩$$

$$\frac{i}{2v\cdot k-\Pi (v\cdot k)},$$

$$2v\cdot \tilde{k}-\mathrm{Re}\Pi (v\cdot \tilde{k})=0.$$

$$m=m_0+v\cdot \tilde{k}=m_0+\frac{1}{2}\mathrm{Re}\Pi (v\cdot \tilde{k}).$$

$$\mathcal{L}$$

$$H_a=\frac{1+v/}{2}[P_{a\mu }^{\ast }{\gamma }^{\mu }-P_a{\gamma }_5],S_a=\frac{1}{2}(1+v/)[{P^{\prime }}_{1a}^{\mu }{\gamma }_{\mu }{\gamma }_5-P_{0a}^{\ast }],$$

$$⟨M_H⟩\equiv \frac{3M_{P^{\ast }}+M_P}{4},⟨M_S⟩\equiv \frac{3M_{P_1^{\prime }}+M_{P_0^{\ast }}}{4},$$

$${\Delta }_S=⟨M_S⟩-⟨M_H⟩.$$

$${\mathcal{L}}_{1/m_Q}=\frac{1}{2m_Q}\left\{{\lambda }_2^H\mathrm{Tr}[{\bar{H}}_a{\sigma }^{\mu \nu }{H}_a{\sigma }_{\mu \nu }]-{\lambda }_2^S\mathrm{Tr}[{\bar{S}}_a{\sigma }^{\mu \nu }{S}_a{\sigma }_{\mu \nu }]\right\},$$

$${\lambda }_2^H=\frac{1}{4}(M_{P^{\ast }}^2-M_P^2),{\lambda }_2^S=\frac{1}{4}(M_{P_1^{\prime }}^2-M_{P_0^{\ast }}^2),$$

$$M_{P_a}$$

$$M_{P_{0a}^{\ast }}$$

$$\frac{i}{2(v\cdot k+\frac{3}{4}\Delta M_P-{\Delta }_a)+iϵ},\frac{-i(g_{\mu \nu }-v_{\mu }{v}_{\nu })}{2(v\cdot k-\frac{1}{4}\Delta M_P-{\Delta }_a)+iϵ},$$

$$\frac{i}{2(v\cdot k-{\Delta }_S+\frac{3}{4}\Delta M_S-{\tilde{\Delta }}_a)+iϵ},\frac{-i(g_{\mu \nu }-v_{\mu }{v}_{\nu })}{2(v\cdot k-{\Delta }_S-\frac{1}{4}\Delta M_S-{\tilde{\Delta }}_a)+iϵ},$$

$${\Pi }_{DK}({\omega }_D)$$

$$\frac{i}{2(v\cdot k-{\Delta }_S+\frac{3}{4}\Delta M_S-{\tilde{\Delta }}_s)-[2{\Pi }_{DK}({\omega }_D)+\frac{2}{3}{\Pi }_{D_s\eta }({\omega }_{D_s})+2{\Pi }_{D_1^{\prime }K}^{\prime }({\omega }_{D_1^{\prime }}^{\prime })+\frac{2}{3}{\Pi }_{D_{s1}^{\prime }\eta }^{\prime }({\omega }_{D_{s1}^{\prime }}^{\prime })]}$$

$${\Pi }_{D_1^{\prime }K}^{\prime }({\omega }_{D_1^{\prime }}^{\prime })$$

$$\displaystyle i\Bigg{[}2(v\cdot k-\Delta_{S}+{3\over 4}\Delta M_{S}-\tilde{\Delta}_{u})-\Big{[}{3\over 2}\,\Pi_{D\pi}(\omega_{D})+{1\over 6}\,\Pi_{D\eta}(\omega_{D})+\Pi_{D_{s}K}(\omega_{D_{s}})$$

$$\displaystyle\qquad+{3\over 2}\,\Pi^{\prime}_{D^{\prime}_{1}\pi}(\omega^{\prime}_{D^{\prime}_{1}})+{1\over 6}\,\Pi^{\prime}_{D^{\prime}_{1}\eta}(\omega^{\prime}_{D^{\prime}_{1}})+\Pi^{\prime}_{D^{\prime}_{s1}K}(\omega^{\prime}_{D^{\prime}_{s1}})\Big{]}\Bigg{]}^{-1}.$$

$${\Delta }_S-\frac{3}{4}\Delta M_S+{\tilde{\Delta }}_s={\mathcal{M}}_{D_{s0}^{\ast }}-M_0={\mathcal{M}}_{D_{s0}^{\ast }}-M_D-\frac{3}{4}\Delta M_D+{\Delta }_u.$$

$$F(v\cdot k{)}_{D_{s0}^{\ast }}$$

$$F(v\cdot k{)}_{D_0^{\ast }}$$

$$M_{D_{(s)0}^{\ast }}=M_0+v\cdot \tilde{k}=M_D+\frac{3}{4}\Delta M_D-{\Delta }_u+v\cdot \tilde{k}.$$

$$\frac{-i(g_{\mu \nu }-v_{\mu }{v}_{\nu })}{2(v\cdot k-{\Delta }_S-\frac{1}{4}\Delta M_S-{\tilde{\Delta }}_s)-[2{\Pi }_{D^{\ast }K}({\omega }_{D^{\ast }})+\frac{2}{3}{\Pi }_{D_s^{\ast }\eta }({\omega }_{D_s^{\ast }})+2{\Pi }_{D_0^{\ast }K}^{\prime }({\omega }_{D_0^{\ast }}^{\prime })+\frac{2}{3}{\Pi }_{D_{s0}^{\ast }\eta }^{\prime }({\omega }_{D_{s0}^{\ast }}^{\prime })]},$$

$$\Pi (\omega )$$

$${\Pi }^{\prime }(\omega )$$

$$\displaystyle F\left({1\over x}\right)=\cases{{\sqrt{x^{2}-1}\over x}\,{\rm ln}(x+\sqrt{x^{2}-1}),&$|x|\geq 1$\cr-{\sqrt{1-x^{2}}\over x}\left[{\pi\over 2}-{\rm tan}^{-1}\left({x\over\sqrt{1-x^{2}}}\right)\right],&$|x|\leq 1$\,.}$$

$$F\left(-\frac{m}{\omega }\right)=\frac{1}{\omega }G(\omega ,m)$$

$$\displaystyle G(\omega,m)=\cases{\sqrt{\omega^{2}-m^{2}}\,[{\rm cosh}^{-1}({\omega\over m})-i\pi],&$\omega>m$\cr\sqrt{m^{2}-\omega^{2}}\,{\rm cos}^{-1}(-{\omega\over m}),&$\omega^{2}<m^{2}$\cr-\sqrt{\omega^{2}-m^{2}}\,{\rm cosh}^{-1}(-{\omega\over m}),&$\omega<-m$\,.}$$

$$\Pi (\omega )$$

$$M_{D_{s0}^{\ast }}={\mathcal{M}}_{D_{s0}^{\ast }}+(v\cdot k{)}_{D_{s0}^{\ast }},M_{D_0^{\ast }}={\mathcal{M}}_{D_0^{\ast }}+(v\cdot k{)}_{D_0^{\ast }},$$

$$(v\cdot k{)}_{D_{s0}^{\ast }}$$

Full text

September, 2017

Masses of Scalar and Axial-Vector $B$ Mesons Revisited

Hai-Yang Cheng1, Fu-Sheng Yu2

1 Institute of Physics, Academia Sinica

Taipei, Taiwan 115, Republic of China

2 School of Nuclear Science and Technology, Lanzhou University

Lanzhou 730000, People’s Republic of China

Abstract

The SU(3) quark model encounters a great challenge in describing even-parity mesons. Specifically, the $q\bar{q}$ quark model has difficulties in understanding the light scalar mesons below 1 GeV, scalar and axial-vector charmed mesons and $1^{+}$ charmonium-like state $X(3872)$. A common wisdom for the resolution of these difficulties lies on the coupled channel effects which will distort the quark model calculations. In this work, we focus on the near mass degeneracy of scalar charmed mesons, $D_{s0}^{*}$ and $D_{0}^{*0}$, and its implications. Within the framework of heavy meson chiral perturbation theory, we show that near degeneracy can be qualitatively understood as a consequence of self-energy effects due to strong coupled channels. Quantitatively, the closeness of $D_{s0}^{*}$ and $D_{0}^{*0}$ masses can be implemented by adjusting two relevant strong couplings and the renormalization scale appearing in the loop diagram. Then this in turn implies the mass similarity of $B_{s0}^{*}$ and $B_{0}^{*0}$ mesons. The $P_{0}^{*}P^{\prime}_{1}$ interaction with the Goldstone boson is crucial for understanding the phenomenon of near degeneracy. Based on heavy quark symmetry in conjunction with corrections from QCD and $1/m_{Q}$ effects, we obtain the masses of $B^{*}_{(s)0}$ and $B^{\prime}_{(s)1}$ mesons, for example, $M_{B_{s0}^{*}}=(5715\pm 1)\,{\rm MeV}+\delta\Delta_{S}$, $M_{B^{\prime}_{s1}}=(5763\pm 1)\,{\rm MeV}+\delta\Delta_{S}$ with $\delta\Delta_{S}$ being $1/m_{Q}$ corrections. We find that the predicted mass difference of 48 MeV between $B^{\prime}_{s1}$ and $B_{s0}^{*}$ is larger than that of $20\sim 30$ MeV inferred from the relativistic quark models, whereas the difference of 15 MeV between the central values of $M_{B^{\prime}_{s1}}$ and $M_{B^{\prime}_{1}}$ is much smaller than the quark model expectation of $60-100$ MeV. Experimentally, it is important to have a precise mass measurement of $D_{0}^{*}$ mesons, especially the neutral one, to see if the non-strange scalar charmed meson is heavier than the strange partner as suggested by the recent LHCb measurement of the $D_{0}^{*\pm}$.

I Introduction

Although the SU(3) quark model has been applied successfully to describe the properties of hadrons such as pseudoscalar and vector mesons, octet and decuplet baryons, it often encounters a great challenge in understanding even-parity mesons, especially scalar ones. Take vector mesons as an example and consider the octet vector ones: $\rho,\omega,K^{*},\phi$. Since the constituent strange quark is heavier than up or down quark by 150 MeV, one will expect the mass hierarchy pattern $m_{\phi}>m_{K^{*}}>m_{\rho}\sim m_{\omega}$ which is borne out by experiment. However, this quark model picture faces great challenges in describing the even-parity meson sector:

• •

Many scalar mesons with masses lower than 2 GeV have been observed and they can be classified into two nonets: one nonet with mass below or close to 1 GeV, such as $f_{0}(500)$ (or $\sigma$), $K_{0}^{*}(800)$ (or $\kappa$), $f_{0}(980)$ and $a_{0}(980)$ and the other nonet with mass above 1 GeV such as $K_{0}^{*}(1430)$, $a_{0}(1450)$ and two isosinglet scalar mesons. Of course, the two nonets cannot be both low-lying ${}^{3}P_{0}$ $q\bar{q}$ states simultaneously. If the light scalar nonet is identified with the P-wave $q\bar{q}$ states, one will encounter two major difficulties: First, why are $a_{0}(980)$ and $f_{0}(980)$ degenerate in their masses? In the $q\bar{q}$ model, the latter is dominated by the $s\bar{s}$ component, whereas the former cannot have the $s\bar{s}$ content since it is an $I=1$ state. One will expect the mass hierarchy pattern : $m_{f_{0}(980)}>m_{K_{0}^{*}(800)}>m_{a_{0}(980)}\sim m_{f_{0}(500)}$. However, this pattern is not seen by experiment. In contrast, it is $m_{a_{0}(980)}\approx m_{f_{0}(980)}>m_{K_{0}^{*}(800)}>m_{f_{0}(500)}$ experimentally. Second, why are $f_{0}(500)$ and $K_{0}^{*}(800)$ so broad compared to the narrow widths of $a_{0}(980)$ and $f_{0}(980)$ even though they are all in the same nonet?

• •

In the scalar meson sector above 1 GeV, $K_{0}^{*}(1430)$ with mass $1425\pm 50$ MeV PDG is almost degenerate in masses with $a_{0}(1450)$ which has a mass of $1474\pm 19$ MeV PDG despite having one strange quark for the former.

• •

In the even-parity charmed meson sector, we compare the experimentally measured masses and widths with what are expected from the quark model (see Table 1). There are some prominent features from this comparison: (i) The measured masses of $D_{0}^{*},D^{\prime}_{1},D_{s0}^{*}$ and $D^{\prime}_{s1}$ are substantially smaller than the quark model predictions. (ii) The physical $D_{s0}^{*}$ mass is below the $DK$ threshold, while $D^{\prime}_{s1}$ is below $DK^{*}$. This means that both of them are quite narrow, in sharp contrast to the quark model expectation of large widths for them. (iii) $D_{0}^{*}(2400)^{0}$ and $D_{s0}^{*}(2317)$ are almost equal in their masses, while $D_{0}^{*}(2400)^{\pm}$ is heavier than $D_{s0}^{*}$ even though the latter contains a strange quark. 111Strictly speaking, the masses of the neutral and charged states of $D_{0}^{*}(2400)$ are not consistently determined due mainly to its broadness. The mass of the charged one is primarily from the LHCb measurements LHCb:D0pm , while the neutral one is from BaBar BaBar:D0 , Belle Belle:D0 and FOCUS FOCUS:D0 . It is worthwhile to notice that only FOCUS has measured both the neutral and charged $D_{0}^{*}(2400)$, and their masses are quite similar with a small difference of a few MeV. All the other three groups have not reported the masses for both the neutral and charged $D_{0}^{*}(2400)$. In addition, the masses reported by different groups are very different from each other. As a result, the world averaged masses for the neutral and charged $D_{0}^{*}(2400)$ shown in Table 1 are very different as well. The difference is well beyond the expectation from isospin splitting. Thus, the two averaged masses given by the Particle Data Group (PDG) PDG are not consistent with each other. In this work, we consider both PDG masses for $D_{0}^{*}(2400)$. In Sec. II we focus on the closeness of $D_{s0}^{*}$ and $D_{0}^{*}$ masses and discuss its implication to the scalar $B$ sector. In Sec. III we derive two different sets of scalar $B$ meson masses, corresponding to two different PDG masses for $D_{0}^{*}(2400)$. (iv) The masses of $D_{1},D_{2}^{*},D_{s1}$ and $D_{s2}^{*}$ predicted by the quark model are consistent with experiment. These four observations lead to the conclusion that $0^{+}$ and $1^{\prime+}$ charmed mesons have very unusual behavior not anticipated from the quark model.

• •

The first $X\!Y\!Z$ particle, namely $X(3872)$, observed by Belle in 2003 in $B^{\pm}\to K^{\pm}+(J/\psi\pi^{+}\pi^{-})$ decays Belle:X , has the quantum numbers $J^{PC}=1^{++}$ LHCb:X . $X(3872)$ cannot be a pure charmonium as it cannot be identified as $\chi_{c1}(1^{3}P_{1})$ with a mass 3511 MeV PDG or $\chi_{c1}(2\,^{3}P_{1})$ with the predicted mass of order 3950 MeV Godfrey:1985xj . Moreover, a pure charmonium for $X(3872)$ cannot explain the large isospin violation observed in $X(3872)\to J/\psi\omega,J/\psi\rho$ decays. The extreme proximity of $X(3872)$ to the threshold suggests a loosely bound molecule state $D^{0}\bar{D}^{*0}$ for $X(3872)$. On the other hand, $X(3872)$ cannot be a pure $D\!\bar{D}^{*}$ molecular state either for the following reasons: (i) It cannot explain the prompt production of $X(3872)$ in high energy collisions Suzuki:2005ha ; Braaten:2009wk . (ii) The ratio $R_{1}\equiv\Gamma(B^{0}\to K^{0}X(3872))/\Gamma(B^{-}\to K^{-}X(3872))$ is predicted to be much less than unity in the molecular scenario, while it was measured to be $0.50\pm 0.30\pm 0.05$ and $1.26\pm 0.65\pm 0.06$ by Belle Belle:Xcs . (iii) For the ratio $R_{2}\equiv\Gamma(X(3872)\to\psi(2S)\gamma)/\Gamma(X(3872)\to J/\psi(1S)\gamma)$, the molecular model leads to a very small value of order $3\times 10^{-3}$ Swanson:2006st ; Dong:2009uf ; Ferretti:2014xqa while the charmonium model predicts $R_{2}$ to be order of unity. The LHCb measurement yields $R_{2}=2.46\pm 0.64\pm 0.29$ LHCb:Xgamma . Hence, $X(3872)$ cannot be a pure $D\!\bar{D}^{*}$ molecular state.

The above discussions suggest that $X(3872)$ is most likely an admixture of the $S$-wave $D\!\bar{D}^{*}$ molecule and the $P$-wave charmonium as first advocated in Suzuki:2005ha

[TABLE]

More specifically, the charmonium is identified with $\chi_{c1}(2\,^{3}P_{1})$. Some calculations favor a larger $c\bar{c}$ component over the $D^{0}\overline{D}^{*0}$ component (see e.g. Matheus ; Kang ). Then the question is how to explain the mass of $X(3872)$ through the charmonium picture.

In short, the $q\bar{q}$ quark model has difficulties in describing light scalar mesons below 1 GeV, $0^{+}$ and $1^{\prime+}$ charmed mesons and $1^{+}$ charmonium-like state $X(3872)$. A common wisdom for the resolution of aforementioned difficulties lies on the coupled channel effects which will distort the quark model calculations.

In the quark potential model, the predicted masses for $D_{s0}^{*}$ and $D_{0}^{*0}$ are higher than the measured ones by order 160 MeV and 70 MeV, respectively Godfrey ; Di Pierro . It was first stressed and proposed in vanBeveren that the low mass of $D_{s0}^{*}(2317)$ ($D_{0}^{*}(2400)^{0}$) arises from the mixing between the $0^{+}$ $c\bar{s}$ ($c\bar{q}$) state and the $DK$ ($D\pi$) threshold (see also Browder ). This conjecture was realized in both QCD sum rule Dai ; Dai:D0 and lattice Mohler:Ds0 ; Lang ; Mohler:D0 calculations. For example, when the contribution from the $DK$ continuum is included in QCD sum rules, it has been shown that this effect will significantly lower the mass of the $D_{s0}^{*}$ state Dai . Recent lattice calculations using $c\bar{s}$, $DK$ and $D^{*}K$ interpolating fields show the existence of $D_{s0}^{*}(2317)$ below the $DK$ threshold Mohler:Ds0 and $D^{\prime}_{s1}(2460)$ below the $D^{*}K$ threshold Lang . 222A recent lattice calculation with $N_{f}=2+1+1$ optimal domain-wall fermions Chen:2017kxr yields a mass of $2317\pm 15\pm 5$ MeV for $D_{s0}^{*}$ and $2463\pm 13\pm 9$ MeV for $D^{\prime}_{s1}(2460)$.

All these results indicate that the strong coupling of scalars with hadronic channels will play an essential role of lowering their masses.

By the same token, mass shifts of charmed and bottom scalar mesons due to self-energy hadronic loops have been calculated in Guo . The results imply that the bare masses of scalar mesons calculated in the quark model can be reduced significantly. Mass shifts due to hadronic loops or strong coupled channels have also been studied in different frameworks to explain the small mass of $D_{s0}^{*}(2317)$ Hwang:2004cd ; Zhou ; TLee ; Ortega:2016mms . In the same spirit, even if $X(3872)$ is dominated by the $c\bar{c}$ component, the mass of $\chi_{c1}^{\prime}$ can be shifted down due to its strong coupling with $D\!\bar{D}^{*}$ channels Li:2009ad ; Danilkin:2010cc ; Achasov:2015oia .

Both $f_{0}(980)$ and $a_{0}(980)$ have the strong couple channel $K\overline{K}$. They are often viewed as $K\overline{K}$ molecules, which accounts for their near degeneracy with $2m_{K}$. Schematically, the self-energy $K\overline{K}$ loop diagram of $a_{0}(980)$ will shift its mass to the physical one. In the unitarized chiral perturbation theory, light scalar mesons $f_{0}(500)$, $K_{0}^{*}(800)$, $f_{0}(980)$ and $a_{0}(980)$ can be dynamically generated through their strong couplings with $\pi\pi$, $K\pi$, $K\overline{K}$, and $K\overline{K}$, respectively lightS . Alternatively, it is well known that the tetraquark picture originally advocated by Jaffe Jaffe provides a simple solution to the mass and width hierarchy problems in the light scalar meson sector. The tetraquark structure of light scalars accounts for the mass hierarchy pattern $m_{a_{0}(980)}\approx m_{f_{0}(980)}>m_{K_{0}^{*}(800)}>m_{f_{0}(500)}$, Moreover, the $S$-wave 4-quark nonet can be lighter than the $P$-wave $q\bar{q}$ nonet above 1 GeV due to the absence of the orbital angular momentum barrier and the presence of strong attraction between the diquarks $(qq)_{\bf 3^{*}}$ and $(\bar{q}\bar{q})_{\bf 3}$ Jaffe2 . The fall-apart decays $f_{0}(500)\to\pi\pi$, $K_{0}^{*}(800)\to\pi K$, and $f_{0}(980),a_{0}(980)\to K\overline{K}$ are all OZI-superallowed. This explains the very broad widths of $f_{0}(500)$ and $K_{0}^{*}(800)$, and the narrowness of $f_{0}(980)$ and $a_{0}(980)$ owing to the very limited phase space available as they are near the $K\overline{K}$ threshold.

In Cheng:2014 we have studied near mass degeneracy of scalar charmed and bottom mesons. Qualitatively, the approximate mass degeneracy can be understood as a consequence of self-energy effects due to strong coupled channels which will push down the mass of the heavy scalar meson in the strange sector more than that in the non-strange partner. However, we showed that it works in the conventional model without heavy quark expansion, but not in the approach of heavy meson chiral perturbation theory (HMChPT) as mass degeneracy and the physical masses of $D_{s0}^{*}$ and $D_{0}^{*}$ cannot be accounted for simultaneously. Mass shifts in the strange charm sector are found to be largely overestimated. It turns out that the conventional model works better toward the understanding of near mass degeneracy.

Our previous work was criticized by Alhakami Alhakami who followed the framework of Mehen to write down the general expression of HMChPT and fit the unknown low-energy constants in the effective Lagrangian to the experimentally measured odd- and even-parity charmed mesons. Using the results from the charm sector, Alhakami predicted the spectrum of odd- and even-parity bottom mesons. He concluded that the near degeneracy of nonstrange and strange scalar $B$ mesons is confirmed in the predictions using HMChPT. He then proceeded to criticize that we should use physical masses instead of bare masses to evaluate the hardronic loop effects and that we have missed the contributions from axial-vector heavy mesons to the self-energy of scalar mesons.

Motivated by the above-mentioned criticisms Alhakami , in this work we shall re-examine our calculations within the framework of HMChPT. We show that the closeness of $D_{s0}^{*}$ and $D_{0}^{*0}$ masses can be achieved by taking into account the additional contribution, which was missing in our previous work, from axial-vector heavy mesons to the self-energy diagrams of scalar mesons by adjusting two relevant strong couplings and the renormalization scale $\mu$ appearing in the loop diagram. Then we proceed to confirm that near degeneracy observed in the charm sector will imply the similarity of $B_{s0}^{*}$ and $B_{0}^{*}$ masses in the $B$ system.

This work is organized as follows. In Sec. II we consider the self-energy corrections to scalar and axial-vector heavy mesons in HMChPT. In the literature, the self-energy loop diagrams were sometimes evaluated in HMChPT by neglecting the corrections from mass splittings and residual masses to the heavy meson’s propagator. We shall demonstrate in Sec. III that the calculation in this manner does not lead to the desired degeneracy in both charm and $B$ sectors simultaneously. The masses of $B_{s0}^{*}$ and $B_{0}^{*}$ are discussed in Sec. IV with focus on the predictions based on heavy quark symmetry and possible $1/m_{Q}$ and QCD corrections. Sec. V comes to our conclusions.

II Mass shift of scalar and axial-vector heavy mesons due to hadronic loops

Self-energy hadronic loop corrections to $0^{+}$ and $1^{\prime+}$ heavy mesons have been considered in the literature Becirevic ; Fajfer ; Brdnik ; Mehen ; Alhakami ; Ananth ; Guo ; Cheng:2014 . Since chiral loop corrections to heavy scalar mesons have both finite and divergent parts, it is natural to consider the framework of HMChPT where the divergences and the renormalization scale dependence arising from the chiral loops induced by the lowest-order tree Lagrangian can be absorbed into the counterterms which have the same structure as the next-order tree Lagrangian.

The heavy meson’s propagator in HMChPT has the expression

[TABLE]

where $v$ and $k$, respectively, are the velocity and the residual momentum of the meson defined by $p=vm_{0}+k$, and $\Pi(v\cdot k)$ is the 1PI self-energy contribution. In general, $\Pi(v\cdot k)$ is complex as its imaginary part is related to the resonance’s width. The particle’s on-shell condition is then given by

[TABLE]

The physical mass reads

[TABLE]

Consider the self-energy diagrams depicted in Fig. 1 for scalar and axial-vector heavy mesons. We will evaluate the loop diagrams in the framework of HMChPT in which the low energy dynamics of hadrons is described by the formalism in which heavy quark symmetry and chiral symmetry are synthesized Yan ; Wise ; Donoghue . The relevant Lagrangian is Casalbuoni

[TABLE]

where $H$ denotes the odd-parity spin doublet $(P,P^{*})$ and $S$ the even-parity spin doublet $(P_{0}^{*},P^{\prime}_{1})$ with $j=1/2$ ($j$ being the total angular momentum of the light degrees of freedom):

[TABLE]

with $a=u,d,s$, for example, $P_{a}=(D^{0},D^{+},D_{s}^{+})$. The nonlinear chiral symmetry is realized by making use of the unitary matrix $\Sigma={\rm exp}(i\sqrt{2}\phi/f_{\pi})$ with $f_{\pi}=93$ MeV and $\phi$ being a $3\times 3$ matrix for the octet of Goldstone bosons. In terms of the new matrix $\xi=\Sigma^{1/2}$, the axial vector field ${\cal A}$ reads ${\cal A}={i\over 2}(\xi^{\dagger}\partial_{\mu}\xi-\xi\partial_{\mu}\xi^{\dagger})$. In Eq. (5), the parameter $\Delta_{S}$ is the residual mass of the $S$ field; it measures the mass splitting between even- and odd-parity doublets and can be expressed in terms of the spin-averaged masses

[TABLE]

so that

[TABLE]

There exist two corrections to the chiral Lagrangian (5): one from $1/m_{Q}$ corrections and the other from chiral symmetry breaking. The $1/m_{Q}$ corrections are given by Casalbuoni ; Boyd

[TABLE]

with

[TABLE]

where $\lambda_{H}$ ($\lambda_{S}$) is the mass splitting between spin partners, namely, $P^{*}$ and $P$ ($P^{\prime}_{1}$ and $P^{*}_{0}$) of the pseudoscalar (scalar) doublet. We will not write down the explicit expressions for chiral symmetry breaking terms and the interested reader is referred to Mehen . The masses of heavy mesons can be expressed as

[TABLE]

where $\Delta_{a}$ and $\tilde{\Delta}_{a}$ denote the residual mass contributions to odd- and even-parity mesons, respectively. Note that $\lambda_{2}^{H}/m_{Q}\approx{1\over 2}(M_{P^{*}}-M_{P})\equiv{1\over 2}\Delta M_{P}$ and $\lambda_{2}^{S}/m_{Q}\approx{1\over 2}(M_{P^{\prime}_{1}}-M_{P_{0}^{*}})\equiv{1\over 2}\Delta M_{S}$ in the heavy quark limit. The propagators for $P_{a}(v)$, $P_{a}^{*}(v)$, $P_{0a}^{*}(v)$ and $P^{\prime}_{1a}(v)$ read

[TABLE]

and

[TABLE]

respectively.

Consider the hadronic loop contribution to $D_{s0}^{*}$ in Fig. 1(a) with the intermediate states $D^{0}$ and $K^{+}$. The self-energy loop integral is

[TABLE]

where $m$ is the mass of the Goldstone boson. The residual momentum $k^{\prime}$ of the heavy meson in the loop is given by $k^{\prime}=p+q-vM_{D}=q+k+v({\cal M}_{D_{s0}^{*}}-M_{D})$, and $\omega_{D}=v\cdot k+{\cal M}_{D_{s0}^{*}}-M_{D}+{3\over 4}\Delta M_{D}-\Delta_{u}$. The calligraphic symbol has been used to denote the bare mass. The full $D_{s0}^{*}$ propagator becomes

[TABLE]

after taking into account the contributions from the channels $D^{0}K^{+},D^{+}K^{0},D_{s}^{+}\eta$ and $D_{1}^{\prime 0}K^{+},D_{1}^{\prime+}K^{0},D^{\prime+}_{s1}\eta$. In Eq. (II),

[TABLE]

where $\omega^{\prime}_{D^{\prime}_{1}}=v\cdot k+{\cal M}_{D_{s0}^{*}}-M_{D^{\prime}_{1}}-(M_{D^{\prime}_{1}}-M_{D})+{3\over 4}\Delta M_{D}-\Delta_{u}$, and use of Eq. (II) has been made. Likewise, the full $D_{0}^{*}$ propagator reads

[TABLE]

Since many parameters such as $\Delta_{S}$, $\Delta M_{S}$, $\tilde{\Delta}_{s}$ and $\tilde{\Delta}_{u}$ in Eqs. (II) and (II) are unknown, we are not able to determine mass shifts from above equations. Assuming that the bare mass ${\cal M}$ is the one obtained in the quark model, then from Eq. (II) we have

[TABLE]

With

[TABLE]

the on-shell conditions read $F(v\cdot\tilde{k})_{D_{s0}^{*}}=0$ for $D_{s0}^{*}$ and $F(v\cdot\tilde{k})_{D_{0}^{*}}=0$ for $D_{0}^{*}$. The physical masses are then given by

[TABLE]

Since $\Delta_{u}$ is of order 1 MeV Stewart , it can be neglected in practical calculations. Note that in the above equation, one should not replace $M_{0}$ by the bare mass ${\cal M}_{D_{(s)0}^{*}}$. Indeed, in the absence of chiral loop corrections, $M_{D_{s0}^{*}}=M_{0}+v\cdot\tilde{k}=M_{0}+\Delta_{S}-{3\over 4}\Delta M_{S}+\tilde{\Delta}_{s}={\cal M}_{D_{s0}^{*}}$, as it should be.

For the self-energy of the axial-vector meson, we consider $D^{\prime}_{s1}$ as an illustration which receives contributions from $D^{*0}K^{+},D^{*+}K^{0},D_{s}^{*+}\eta$ and $D_{0}^{*0}K^{+},D_{0}^{*+}K^{0},D^{*+}_{s0}\eta$ intermediate states (see Fig. 1(b)). The full $D^{\prime}_{s1}$ propagator reads

[TABLE]

where $\omega_{D^{*}}=v\cdot k+M_{D^{\prime}_{s1}}-M_{D^{*}}-{1\over 4}\Delta M_{D}$ and $\omega^{\prime}_{D_{0}^{*}}=v\cdot k+M_{D^{\prime}_{s1}}-M_{D^{*}_{0}}-(M_{D_{0}^{*}}-M_{D})+{3\over 4}\Delta M_{D}$.

The loop integrals in Eqs. (14) and (16) have the expressions Falk ; Boyd ; Stewart

[TABLE]

and Mehen ; Alhakami

[TABLE]

respectively, with

[TABLE]

Note that the function $F(-m/\omega)$ can be recast to the form

[TABLE]

with

[TABLE]

The parameter $\Lambda$ appearing in Eqs. (22) and (23) is an arbitrary renormalization scale. In the dimensional regularization approach, the common factor ${2\over\epsilon}-\gamma_{E}+{\rm ln}4\pi+1$ with $\epsilon=4-n$ can be lumped into the logarithmic term ${\rm ln}(\Lambda^{2}/m^{2})$. In the conventional practice, it is often to choose $\Lambda\sim\Lambda_{\chi}$, the chiral symmetry breaking scale of order 1 GeV, to get numerical estimates of chiral loop effects. However, as pointed out in Cheng:2014 , contrary to the common wisdom, the renormalization scale has to be larger than the chiral symmetry breaking scale of order 1 GeV in order to satisfy the on-shell conditions. In general, there exist two solutions for $v\cdot\tilde{k}$ due to two intercepts of the curve with the $v\cdot k$ axis. We shall consider the smaller solution for $v\cdot\tilde{k}$ as the other solution will yield too large masses. It could be that higher-order heavy quark expansion needs to be taken into account to justify the use of $\Lambda\sim\Lambda_{\chi}$.

In our previous study Cheng:2014 we argued that near mass degeneracy and the physical masses of $D_{s0}^{*}$ and $D_{0}^{*}$ cannot be accounted for simultaneously in the approach of HMChPT. In this work we show that near mass degeneracy can be implemented by taking into account the additional contributions from axial-vector heavy mesons to the self-energy diagram Fig. 1(a) of scalar mesons. Since $\omega_{P}\sim{\cal M}_{P_{0}^{*}}-M_{P}+\cdots$, $\omega^{\prime}_{P^{\prime}_{1}}\sim{\cal M}_{P_{0}^{*}}-M_{P^{\prime}_{1}}-(M_{P^{\prime}_{1}}-M_{P})+\cdots$ and ${\cal M}_{P_{0}^{*}}>M_{P}$, ${\cal M}_{P_{0}^{*}}<M_{P^{\prime}_{1}}$, we find numerically that $\Pi^{\prime}(\omega^{\prime})$ contributes destructively to the mass shifts. Moreover, the self-energy of $D_{0}^{*}$ ($D_{s0}^{*}$) is sensitive (insensitive) to the coupling $g^{\prime}$. Therefore, we can adjust the couplings $h$, $g^{\prime}$ and the renormalization scale $\mu$ to get a mass degeneracy of $D_{s0}^{*}$ and $D_{0}^{*}$. Take the quark model of Godfrey as an illustration which predicts the bare masses of 2480 and 2400 MeV for $D_{s0}^{*}$ and $D_{0}^{*}$, respectively. We find that $h=0.51$, $g^{\prime}=0.25$ and $\Lambda=1.27$ GeV will lead to $M_{D_{s0}^{*}}=2222$ MeV and $M_{D_{0}^{*}}=2225$ MeV. By fixing $h$ and $\Lambda$ and varying the coupling $g^{\prime}$, we obtain $(M_{D_{s0}^{*}},M_{D_{0}^{*}})=(2338,2241)$, $(2277,2235)$ and $(2222,2225)$ MeV for $g^{\prime}=0$, 0.15 and 0.25, respectively. This explains why we cannot obtain near degeneracy of $D_{s0}^{*}$ and $D_{0}^{*}$ in our previous work with $g^{\prime}=0$. Hence, although the coupling $g^{\prime}$ which characterizes the interaction of $P^{*}_{0}P^{\prime}_{1}$ with the Goldstone boson is small, it does play a crucial role for achieving near degeneracy. The absolute mass is not an issue as we can scale up the unknown mass $M_{0}$ by an amount of 96 MeV for scalar and axial-vector heavy mesons, i.e. $P_{0a}^{*}$ and $P^{\prime}_{1a}$ in Eq. (II), so that we have $M_{D_{s0}^{*}}=2318$ MeV and $M_{D_{0}^{*}}=2321$ MeV.

Using the same set of values for the parameters $h$ and $g^{\prime}$ but $\Lambda=1.18$ GeV for even-parity $B$ mesons, we see from Table 2 that the closeness of $B_{s0}^{*}$ and $B_{0}^{*}$ masses also holds in the scalar $B$ sector. The predicted masses for axial-vector charmed mesons $M_{D^{\prime}_{s1}}=2452$ MeV and $M_{D^{\prime}_{1}}=2419$ MeV are in agreement with the respective measured masses (see Table 1): $2459.5\pm 0.6$ MeV and $2427\pm 36$ MeV.

The coupling constant $h$ at tree level can be extracted from the measured widths of $D_{0}^{*0},D_{0}^{*\pm}$ and ${D^{\prime}}_{1}^{0}$. We find the values $h=0.60\pm 0.07$ from $D_{0}^{*}(2400)^{0}$, $0.514\pm 0.017$ from $D_{0}^{*}(2400)^{\pm}$ and $0.79\pm 0.17$ from $D^{\prime}_{1}(2430)^{0}$. It is obvious that $D_{0}^{*}(2400)^{\pm}$ gives the best determination of $h$ as its mass and width measured recently by LHCb LHCb:D0pm are more accurate than $D_{0}^{*0}$ and ${D^{\prime}}_{1}^{0}$. Our result of $h=0.508$ agrees well with the LHCb experiment. As for the coupling $g^{\prime}$, its magnitude and even the sign relative to the coupling $g$ in Eq. (5) are unknown. A recent lattice calculation yields $h=0.84(3)(2)$ and $g^{\prime}=-0.122(8)(6)$ Blossier .

Several remarks are in order. (i) The chiral loop calculations presented here are meant to demonstrate that the mass difference between strange heavy scalar mesons and their non-strange partners can be substantially reduced by self-energy contributions. Near mass degeneracy in the scalar charm sector will imply the same phenomenon in the scalar $B$ system. (ii) For the bare masses of even-parity heavy mesons, we have used the relativistic quark model of Godfrey as an illustration. If we use the quark model predictions from Di Pierro , near degeneracy can be implemented by employing $h=0.52$ and $g^{\prime}=0.28$ and similar $\Lambda$ given before. (iii) So far we have focused on the closeness of $D_{s0}^{*}$ and $D_{0}^{*0}$ masses . The measured mass of $D_{0}^{*\pm}$ by LHCb is greater than $D_{0}^{*0}$ and $D_{s0}^{*}$ by order 30 MeV. More precise measurement of the $D_{0}^{*0}$ mass is certainly needed. If it turns out that the mass of $D_{0}^{*0}$ is larger than that of $D_{s0}^{*}$, then we need a larger $|g^{\prime}|$ to render $D_{0}^{*0}$ heavier than $D_{s0}^{*}$. (iv) There exist other contributions which are not considered in Fig. 1, for example, the intermediate states $D^{0}\overline{K}^{*0}$, $D_{s}\eta^{\prime},\cdots$ for the self-energy diagram of $D_{s0}^{*}$. (v) The calculated masses are sensitive to the choice of the renormalization scale $\Lambda$. Of course, physics should be independent of the renormalization scale. However, this issue cannot be properly addressed in the phenomenological model discussed here. In view of this, we should not rely on the chiral loop calculations to make quantitative predictions on the scalar meson masses. For this purpose, we will turn to heavy quark symmetry in Sec. IV below.

III Comparison with other works in HMChPT

Self-energy hadronic loop corrections to $0^{+}$ and $1^{\prime+}$ heavy mesons have been discussed in the literature within the framework of HMChPT Becirevic ; Fajfer ; Brdnik ; Mehen ; Alhakami ; Ananth ; Guo ; Cheng:2014 . Guo, Krewald and Meißner (GKM) Guo considered three different models for calculations. Models I and III correspond to non-derivative and derivative couplings of the scalar meson with two pseudoscalar mesons, while Model II is based on HMChPT. However, our HMChPT results are different from GKM. This can be traced back to the loop integral $\Pi(v\cdot k)$ which has the following expression in Guo

[TABLE]

where use of Eq. (25) has been made. Comparing with Eq. (22), we see that the $m^{2}$ coefficient in Eq. (27) should be replaced by $\omega^{2}$ and a term $(m^{2}/\omega){\rm ln}(\Lambda^{2}/m^{2})$ in $[\cdots]$ is missing. Contrary to the claim made by GKM, the self-energy contribution should not vanish in the chiral limit. Some other detailed comparisons are referred to Cheng:2014 . Contributions from axial-vector heavy mesons to the self-energy of scalar mesons were also not considered in this work.

Mehen and Springer Mehen have systematically studied chiral loop corrections to the masses of scalar and axial-vector heavy mesons. For the propagators of heavy mesons, they only keep the $i/(2v\cdot k)$ terms and neglect all mass splittings such as $\Delta_{S}$, $\Delta M_{P,S}$ and residual masses such as $\Delta_{a}$ and $\tilde{\Delta}_{a}$ in Eqs. (12) and (13). Then they set $v\cdot k=0$ inside the loop integral. As a consequence, $\omega=M_{\rm ext}-M_{\rm int}$, the mass difference between external and internal heavy meson states. If we follow this prescription and repeat the calculations presented in the last section, we will obtain different results. The physical masses of scalar charmed mesons read

[TABLE]

with (see Eq. (II))

[TABLE]

and $\omega_{D}={\cal M}_{D_{(s)0}^{*}}-M_{D}$, $\omega_{D^{\prime}_{1}}={\cal M}_{D_{(s)0}^{*}}-M_{D^{\prime}_{1}}$, $\cdots$ etc. The results of calculations are summarized in Table 3. It is clear that while near degeneracy is achieved for scalar charmed mesons, it is not the case for scalar $B$ mesons.

Note that in the work of Mehen and Alhakami , the bare mass ${\cal M}$ in the mass relation $M={\cal M}+\Pi$ is expressed in terms of the unknown low-energy constants in the effective chiral Lagrangian which are fitted to the measured masses of odd- and even-parity charmed mesons. In our work, we have $M={\cal M}+{1\over 2}{\rm Re}\Pi$ (see Eq. (4)) with ${\cal M}$ being inferred from the quark model. We use full propagators to evaluate the mass shifts. Based on the HMChPT framework of Mehen and Springer Mehen , Alhakami Alhakami has shown that the nonstrange and strange scalar $B$ mesons are nearly degenerate and the mass difference between $B_{0}^{*}$ and $B_{s0}^{*}$ is $\sim$ 8 MeV. However, the renormalization scale is chosen to be $\Lambda=317$ MeV in Alhakami which is too small according to the spirit of ChPT.

IV Masses of $B^{*}_{(s)0}$ and $B^{\prime}_{(s)1}$ mesons

The states $B^{*}_{(s)0}$ and $B^{\prime}_{(s)1}$ are yet to be observed. The predictions of their masses in the literature are collected in Table 4. They can be classified into the following categories: (i) relativistic quark potential model Di Pierro ; Godfrey ; Swanson ; Liu ; Ebert ; Zhu ; Godfrey:2016 , (ii) nonrelativistic quark model Lu ; Ortega , (iii) heavy meson chiral perturbation theory Bardeen ; Colangelo ; Cheng:2014 , (iv) unitarized chiral perturbation theory Rincon ; Cleven:2010aw ; Guo:dynamic ; Guo:Ds1 ; Altenbuchinger ; Cleven:2014 , (v) lattice QCD Gregory ; Koponen ; Lang:Bs , (vi) potential model with one loop corrections Lee ; Swanson ; Albaladejo , (vii) QCD sum rules ZGWang , and (viii) others, such as the nonlinear chiral SU(3) model Lutz:2003fm , the semi-relativistic quark potential model Matsuki , the MIT bag model Orsland , the mixture of conventional $P$-wave quark-antiquark states with four-quark components Vijande , the chiral quark-pion Lagrangian with strong coupled channels Badalian , heavy quark symmetry and the assumption of flavor independence of mass differences between $0^{+}$ and $0^{-}$ states in Dmitrasinovic:2012zz .

It is clear from Table 4 that the mass difference between strange and non-strange scalar $B$ mesons predicted by the relativistic quark models is of order $60-110$ MeV, and that $B_{s0}^{*}$ is above the $BK$ threshold. As shown in Sect. II, the closeness of $B_{0}^{*}$ and $B_{s0}^{*}$ masses is expected in view of the near degeneracy observed in the charm sector.

Since this section based on heavy quark symmetry was already discussed in details in our previous work Cheng:2014 , we shall recapitulate the main points and update the numerical results as the masses of charged and neutral $D_{0}^{*}$ under the current measurement are somewhat different.

In the heavy quark limit, the two parameters $\Delta_{S}$ and $\lambda_{2}^{S}$ defined in Eqs. (8) and (10), respectively, are independent of the heavy quark flavor, where $\Delta_{S}$ measures the spin-averaged mass splitting between the scalar doublet $(P^{\prime}_{1},P^{*}_{0})$ and the pseudoscalar doublet $(P^{*},P)$ and $\lambda_{S}$ is the mass splitting between spin partners, namely $P^{\prime}_{1}$ and $P^{*}_{0}$, of the scalar doublet. From the data listed in Table 1, we are led to

[TABLE]

for the scalar $c\bar{u}$ meson,

[TABLE]

for the scalar $c\bar{d}$ meson and

[TABLE]

for the $c\bar{s}$ meson, where we have assumed that the charged $D^{\prime\pm}_{1}$ has the same mass as $D^{\prime 0}_{1}$. The large errors associated with $\Delta_{S}$ and $\lambda_{2}^{S}$ in the non-strange scalar meson sector reflect the experimental difficulty in identifying the broad states. The light quark flavor dependence of $\Delta_{S}$ and $\lambda_{2}^{S}$ shown in Eqs. (30)$-$(32) indicates SU(3) breaking effects.

Note that the parameter $\Delta_{S}$ can be decomposed in terms of the mass differences between $0^{+}$ and $0^{-}$ states and between $1^{\prime+}$ and $1^{-}$ states defined by

[TABLE]

We have

[TABLE]

If we follow Dmitrasinovic:2012zz to assume the heavy flavor independence of $\Delta_{0}$ and $\Delta_{1}$, namely,

[TABLE]

we will obtain (in MeV)

[TABLE]

Our results for $B^{*0}_{0}$ and $B^{\prime 0}_{1}$ masses are slightly different from those shown in the last row of Table 4. We see that the central values of $B^{*0}_{0}$ and $B^{*}_{s0}$ are not very close. It appears that the assumption of heavy quark flavor independence of $\Delta_{0}$ and $\Delta_{1}$ is probably too strong. Hence, we shall follow Colangelo to assume the heavy quark flavor independence of $\Delta_{S}$ and $\lambda_{2}^{S}$,

[TABLE]

In this case, we find the results (in MeV)

[TABLE]

Predictions along this approach were first obtained in Colangelo . Evidently, mass degeneracy in the scalar $D$ sector will repeat itself in the $B$ system through heavy quark symmetry. A scrutiny of Table 4 shows that this degeneracy is respected only in few of the predictions such as Cheng:2014 , Alhakami , Colangelo , Lee , and Dmitrasinovic:2012zz .

Though the relations in Eq. (37) are valid in the heavy quark limit, they do receive $1/m_{Q}$ and QCD corrections. The leading QCD correction to the relation $\lambda_{2}^{S(b\bar{q})}=\lambda_{2}^{S(c\bar{q})}$ is given by Amoros:1997rx

[TABLE]

and a similar relation with the replacement of $\bar{q}$ by $\bar{s}$. We then obtain

[TABLE]

and the masses (in MeV)

[TABLE]

Comparing with Eq. (IV) we see that QCD corrections will mainly enhance the masses of $B_{0}^{*}$ and $B_{s0}^{*}$ by an amount of a few MeV.

We next proceed to the corrections to $\Delta_{S}^{(b)}=\Delta_{S}^{(c)}$. In heavy quark effective theory, it follows that (see Cheng:2014 for details)

[TABLE]

While the parameter $\lambda_{1}^{H}$ is calculable, the other parameter $\lambda_{1}^{S}$ is unknown. In Cheng:2014 the $1/m_{Q}$ correction is estimated to be

[TABLE]

to be compared with the estimate of $-50\pm 25$ MeV in Mehen .

As pointed out in Cheng:2014 , so long as $\delta\Delta_{S}$ is not large in magnitude, the $1/m_{Q}$ correction to the relation $\Delta_{S}^{(b)}=\Delta_{S}^{(c)}$ [see Eq. (42)], to a very good approximation, amounts to lowering the masses of $B^{*}_{(s)0}$ and $B^{\prime}_{(s)1}$ by an equal amount of $|\delta\Delta_{S}|$:

[TABLE]

Evidently, the mass degeneracy between $B_{s0}^{*}$ and $B_{0}^{*}$ shown in Eq. (IV) implied by heavy quark symmetry is not spoiled by $1/m_{Q}$ and QCD corrections.

Several remarks are in order. (i) There exists an empirical mass relation

[TABLE]

for charmed mesons. Therefore, one will expect the mass difference relation

[TABLE]

also valid in the $B$ sector. The mass difference between $B^{\prime}_{s1}$ and $B_{s0}^{*}$ in Eq. (IV) is independent of $\delta\Delta_{S}$; its numerical value $48\pm 1$ MeV indeed respects the relation (46). It is interesting that this mass relation is not obeyed by any of the existing relativistic quark models which predict a smaller mass difference $M_{B^{\prime}_{s1}}-M_{B^{*}_{s0}}\sim 20-30$ MeV (see Table 4). (ii) In the charm sector we have the hierarchy $M_{D^{\prime}_{s1}}-M_{D^{*}_{s0}}>M_{D^{\prime}_{1}}-M_{D^{*}_{0}}$. It is natural to generalize this to the $B$ sector: $M_{B^{\prime}_{s1}}-M_{B^{*}_{s0}}>M_{B^{\prime}_{1}}-M_{B^{*}_{0}}$. Indeed, the mass splittings of the $j={1\over 2}$ doublet in strange and nonstrange sectors satisfy the relation (see Cheng:2014 )

[TABLE]

The relation $M_{B^{\prime}_{s1}}-M_{B^{*}_{s0}}>M_{B^{\prime}_{1}}-M_{B^{*}_{0}}$ is not respected by some of the models listed in Table 4. (iii) The difference 15 MeV between the central values of $M_{B^{\prime}_{s1}}$ and $M_{B^{\prime}_{1}}$ inferred from Eq. (IV) is much smaller than the quark model expectation of $60-100$ MeV (cf Table 4). 333The splitting between strange and nonstrange axial-vector bottom mesons is of order 19 MeV in the HHChPT approach of Alhakami . Just as the scalar $B$ mesons, the mass splitting is reduced by the self-energy effects due to strong coupled channels. (iv) The current world-averaged masses of neutral and charged $D_{0}^{*}(2400)$ mesons are not the same; they differ by an amount of $33\pm 29$ MeV and the charged one is better measured. On the experimental side, it is thus important to have a more precise mass measurement of the scalar charmed meson. If it turns out the $D_{0}^{*0}$ mass is also of order 2350 MeV, we find from Eq. (IV) that the $B_{0}^{*}$ mass is larger than that of $B_{s0}^{*}$ by an order of 13 MeV.

Since the masses of $B_{s0}^{*}$ and $B^{\prime}_{s1}$ are below the $BK$ and $B^{*}K$ thresholds, respectively, their widths are expected to be very narrow with the isospin-violating strong decays into $B_{s}\pi^{0}$ and $B_{s}^{*}\pi^{0}$. Experimentally, it will be even more difficult to identify the non-strange $B_{0}^{*}$ and $B^{\prime}_{1}$ mesons owing to their broad widths.

V Conclusions

The $q\bar{q}$ quark model encounters a great challenge in describing even-parity mesons. Specifically, it has difficulties in understanding the light scalar mesons below 1 GeV, scalar and axial-vector charmed mesons and $1^{+}$ charmonium-like state $X(3872)$. Many studies indicate that quark model results are distorted by the interaction of the even-parity meson with strong coupled channels. In this work, we focus on the near mass degeneracy of scalar charmed mesons, $D_{s0}^{*}$ and $D_{0}^{*0}$, and its implications to the $B$ sector.

• •

We work in the framework of HMChPT to evaluate the self-energy diagrams for $P_{0}^{*}$ and $P^{\prime}_{1}$ heavy mesons. The approximate mass degeneracy can be qualitatively understood as a consequence of self-energy hadronic loop corrections which will push down the mass of the heavy scalar meson in the strange sector more than that in the non-strange one. Quantitatively, the closeness of $D_{s0}^{*}$ and $D_{0}^{*0}$ can be implemented by adjusting the strong couplings $h$ and $g^{\prime}$ and the renormalization scale $\Lambda$. This in turn implies the mass similarity of $B_{s0}^{*}$ and $B_{0}^{*0}$. The $P_{0}^{*}P^{\prime}_{1}$ interaction with the Goldstone boson characterized by the strong coupling $g^{\prime}$ is crucial for understanding the phenomenon of near degeneracy.

• •

The self-energy loop diagram has been evaluated in the literature by neglecting mass splittings and residual masses in the heavy meson’s propagator. If we follow this prescription, we find near degeneracy in the scalar charm sector but not in the corresponding $B$ system.

• •

The masses of $B^{*}_{(s)0}$ and $B^{\prime}_{(s)1}$ mesons (i.e. Eq. (IV)) can be deduced from the charm spectroscopy and heavy quark symmetry with corrections from QCD and $1/m_{Q}$ effects, for example, $M_{B_{s0}^{*}}=(5715\pm 1)\,{\rm MeV}+\delta\Delta_{S}$, $M_{B^{\prime}_{s1}}=(5763\pm 1)\,{\rm MeV}+\delta\Delta_{S}$ with $\delta\Delta_{S}$ being $1/m_{Q}$ corrections. We found that the mass difference of 48 MeV between $B^{\prime}_{s1}$ and $B_{s0}^{*}$ is larger than the result of $20\sim 30$ MeV predicted by the quark models, whereas the difference of 15 MeV between the central values of $M_{B^{\prime}_{s1}}$ and $M_{B^{\prime}_{1}}$ is much smaller than the quark model expectation of $60-100$ MeV.

• •

The current world-averaged masses of neutral and charged $D_{0}^{*}(2400)$ mesons are not the same; they differ by an amount of $33\pm 29$ MeV and the charged one is more well measured. Experimentally, it is thus important to have a more precise mass measurement of the scalar charmed meson, especially the neutral one. If it turns out the $D_{0}^{*0}$ mass is of order 2350 MeV, this means $D^{*}_{0}$ is heavier than $D_{s0}^{*}$ even though the latter contains a strange quark. In this case, we find that the $B_{0}^{*}$ mass is larger than that of $B_{s0}^{*}$ by an order of 13 MeV.

Acknowledgements.

We wish to thank Ting-Wai Chiu and Keh-Fei Liu for insightful discussions. This research was supported in part by the Ministry of Science and Technology of R.O.C. under the Grant No. 104-2112-M-001-022 and by the National Science Foundation of China under the Grant No. 11347027.