Masses of scalar and axial-vector B mesons revisited

The SU(3) quark model encounters a great challenge in describing even-parity mesons. Specifically, the qq¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q\bar{q}$$\end{document} quark model has difficulties in understanding the light scalar mesons below 1 GeV, scalar and axial-vector charmed mesons and 1+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1^+$$\end{document} 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, Ds0∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_{s0}^*$$\end{document} and D0∗0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_0^{*0}$$\end{document}, 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 Ds0∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_{s0}^*$$\end{document} and D0∗0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_0^{*0}$$\end{document} 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 Bs0∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{s0}^*$$\end{document} and B0∗0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_0^{*0}$$\end{document} mesons. The P0∗P1′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_0^* P'_1$$\end{document} 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/mQ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/m_Q$$\end{document} effects, we obtain the masses of B(s)0∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^*_{(s)0}$$\end{document} and B(s)1′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B'_{(s)1}$$\end{document} mesons, for example, MBs0∗=(5715±1)MeV+δΔS\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{B_{s0}^*}= (5715\pm 1)\,\mathrm{MeV}+\delta \Delta _S$$\end{document}, MBs1′=(5763±1)MeV+δΔS\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{B'_{s1}}=(5763\pm 1)\,\mathrm{MeV}+\delta \Delta _S$$\end{document} with δΔS\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta \Delta _S$$\end{document} being 1/mQ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/m_Q$$\end{document} corrections. We find that the predicted mass difference of 48 MeV between Bs1′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B'_{s1}$$\end{document} and Bs0∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{s0}^*$$\end{document} is larger than that of 20–30 MeV inferred from the relativistic quark models, whereas the difference of 15 MeV between the central values of MBs1′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{B'_{s1}}$$\end{document} and MB1′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{B'_1}$$\end{document} is much smaller than the quark model expectation of 60–100 MeV. Experimentally, it is important to have a precise mass measurement of D0∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_0^*$$\end{document} 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 D0∗±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_0^{*\pm }$$\end{document}.


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 evenparity mesons, especially scalar ones. Take vector mesons as an example and consider the octet vector ones: ρ, ω, K * , φ. Since the constituent strange quark is heavier than the up or down quark by 150 MeV, one will expect the mass hierarchy pattern m φ > m K * > m ρ ∼ m ω , 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 σ ), K * 0 (800) (or κ), 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 lowlying 3 P 0 qq states simultaneously. If the light scalar nonet is identified with the P-wave qq states, one will encounter two major difficulties: first, why are a 0 (980) and f 0 (980) degenerate in their masses? In the qq model, the latter is dominated by the ss component, whereas the former cannot have the ss 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) ∼ m f 0 (500) . However, this pattern is not seen by experiment. In contrast, it is m a 0 (980) ≈ 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 ± 50 MeV [1] is almost degenerate in masses with a 0 (1450), which has a mass of 1474 ± 19 MeV [1] despite having one strange quark for the former. Table 1 Measured masses and widths of even-parity charmed mesons. The four p-wave charmed meson states are denoted by D * 0 , D 1 , D 1 and D * 2 , respectively. In the heavy quark limit, D 1 has j = 1/2 and D 1 has j = 3/2 with j being the total angular momentum of the light degrees of freedom. The data are taken from the Particle Data Group [1]. The last two columns are the predictions from the quark model [7,8]. "Large" means a broad width of order 100 MeV, while "small" implies a narrow width of order 10 MeV • 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 1 , D * s0 and D s1 are substantially smaller than the quark model predictions. (ii) The physical D * s0 mass is below the DK threshold, while D 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) ± is heavier than D * s0 even though the latter contains a strange quark. 1 (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 + charmed mesons have very unusual behavior not anticipated from the quark model. • The first XYZ particle, namely X (3872), observed by Belle in 2003 in B ± → K ± + (J/ψπ + π − ) decays [9], has the quantum numbers J PC = 1 ++ [10]. X (3872) cannot be a pure charmonium as it cannot be identified as 1 Strictly 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 [2,3], while the neutral one is from BaBar [4], Belle [5] and FOCUS [6]. It is worthwhile to notice that only FOCUS has measured both the neutral and the 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 the 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) [1] are not consistent with each other. In this work, we consider both PDG masses for D * 0 (2400). In Sect. 2 we focus on the closeness of D * s0 and D * 0 masses and discuss its implication to the scalar B sector. In Sect. 3 we derive two different sets of scalar B meson masses, corresponding to two different PDG masses for D * 0 (2400). χ c1 (1 3 P 1 ) with a mass 3511 MeV [1] or χ c1 (2 3 P 1 ) with the predicted mass of order 3950 MeV [11]. Moreover, a pure charmonium for X (3872) cannot explain the large isospin violation observed in X (3872) → J/ψω, J/ψρ decays. The extreme proximity of X (3872) to the threshold suggests a loosely bound molecule state D 0D * 0 for X (3872). On the other hand, X (3872) cannot be a pure DD * molecular state either for the following reasons: (i) It cannot explain the prompt production of X (3872) in high-energy collisions [12,13]. (ii) The ratio is predicted to be much less than unity in the molecular scenario, while it was measured to be 0.50 ± 0.30 ± 0.05 and 1.26 ± 0.65 ± 0.06 by Belle [14,15]. (iii) For the ratio R 2 ≡ (X (3872) → ψ(2S)γ )/ (X (3872) → J/ψ(1S)γ ), the molecular model leads to a very small value of order 3 × 10 −3 [16][17][18], while the charmonium model predicts R 2 to be of order of unity. The LHCb measurement yields R 2 = 2.46 ± 0.64 ± 0.29 [19]. Hence, X (3872) cannot be a pure DD * molecular state. The above discussions suggest that X (3872) is most likely an admixture of the S-wave DD * molecule and the Pwave charmonium as first advocated in [12], More specifically, the charmonium is identified with χ c1 (2 3 P 1 ). Some calculations favor a larger cc component over the D 0 D * 0 component (see e.g. [20,21]).
Then the question is how to explain the mass of X (3872) through the charmonium picture.
In short, the qq quark model has difficulties in describing light scalar mesons below 1 GeV, 0 + and 1 + charmed mesons and 1 + charmonium-like state X (3872). Common wisdom for the resolution of aforementioned difficulties lies in 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 and 70 MeV, respectively [7,8]. It was first stressed and proposed in [22] that the low mass of D * s0 (2317) (D * 0 (2400) 0 ) arises from the mixing between the 0 + cs (cq) state and the DK (Dπ ) threshold (see also [23]). This conjecture was realized in both QCD sum rule [24,25] and lattice [26][27][28][29] 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 [24]. Recent lattice calculations using cs, DK and D * K interpolating fields show the existence of D * s0 (2317) below the DK threshold [26,27] and D s1 (2460) below the D * K threshold [28]. 2 All these results indicate that the strong coupling of scalars with hadronic channels will play the 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 [31]. 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) [32][33][34][35]. In the same spirit, even if X (3872) is dominated by the cc component, the mass of χ c1 can be shifted down due to its strong coupling with DD * channels [36][37][38][39][40].
Both f 0 (980) and a 0 (980) have the strong couple channel K K . They are often viewed as K K molecules, which accounts for their near degeneracy with 2m K . Schematically, the self-energy K 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 ππ, K π , K K , and K K , respectively [41][42][43][44][45]. Alternatively, it is well known that the tetraquark picture originally advocated by Jaffe [46] 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) ≈ m f 0 (980) > m K * 0 (800) > m f 0 (500) , Moreover, the S-wave 4quark nonet can be lighter than the P-wave qq nonet above 1 GeV due to the absence of the orbital angular momentum barrier and the presence of strong attraction between the diquarks (qq) 3 * and (qq) 3 [47]. The fall-apart decays f 0 (500) → ππ, K * 0 (800) → π K , and f 0 (980), a 0 (980) → K 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 K threshold.
In [48] we have studied near mass degeneracy of scalar charmed and bottom mesons. Qualitatively, the approximate mass degeneracy can be understood as a consequence of selfenergy 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 [49] who followed the framework of [50] 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 non-strange 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 hadronic loop effects and that we have missed the contributions from axialvector heavy mesons to the self-energy of scalar mesons.
Motivated by the above-mentioned criticisms [49], 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 μ 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 Sect. 2 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 Sect. 3 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 Sect. 4 with focuses on the predictions based on heavy quark symmetry and possible 1/m Q and QCD corrections. Section 5 comes to our conclusions.
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 where v and k, respectively, are the velocity and the residual momentum of the meson defined by p = vm 0 +k, and (v·k) is the 1PI self-energy contribution. In general, (v · k) is complex as its imaginary part is related to the resonance's width. The particle's on-shell condition is then given by The physical mass reads 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 [55][56][57]. The relevant Lagrangian is [58] where H denotes the odd-parity spin doublet (P, P * ) and S the even-parity spin doublet (P * 0 , P 1 ) with j = 1/2 ( j being the total angular momentum of the light degrees of freedom): The nonlinear chiral symmetry is realized by making use of the MeV and φ being a 3 × 3 matrix for the octet of Goldstone bosons. In terms of the new matrix ξ = 1/2 , the axial-vector field A reads , the parameter 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, There exist two corrections to the chiral Lagrangian (2.4): one from 1/m Q corrections and the other from chiral symmetry breaking. The 1/m Q corrections are given by [58,59] where λ H (λ S ) is the mass splitting between spin partners, namely, P * and P (P 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 [50]. The masses of the heavy mesons can be expressed as where a and˜ a denote the residual mass contributions to odd-and even-parity mesons, respectively. Note that 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 , (2.13) where m is the mass of the Goldstone boson. The residual momentum k of the heavy meson in the loop is given by The calligraphic symbol has been used to denote the bare mass. The full D * s0 propagator becomes (2.14) after taking into account the contributions from the chan- In Eq. (2.14), where and use of Eq. (2.10) has been made. Likewise, the full D * 0 propagator reads (2.16) Since many parameters such as S , M S ,˜ s and˜ u in Eqs. (2.14) and (2.16) are unknown, we are not able to determine the mass shifts from the above equations. Assuming that the bare mass M is the one obtained in the quark model, then from Eq. (2.10) we have The physical masses are then given by , as it should. For the self-energy of the axial-vector meson, we consider D s1 as an illustration which receives contributions from Fig. 1b). The full D s1 propagator reads where and [49,50] ( respectively, with (2.23) Note that the function F(−m/ω) can be recast into the form (2.25) The parameter appearing in Eqs. (2.21) and (2.22) is an arbitrary renormalization scale. In the dimensional regularization approach, the common factor 2 − γ E + ln4π + 1 with = 4 − n can be lumped into the logarithmic term ln( 2 /m 2 ). In conventional practice, one often chooses ∼ χ , the chiral symmetry breaking scale of order 1 GeV, to get numerical estimates of chiral loop effects. However, as pointed out in [48], 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 ·k due to two intercepts of the curve with the v · k axis. We shall consider the smaller solution for v ·k as the other solution will yield too large masses. It could be that a higherorder heavy quark expansion needs to be taken into account to justify the use of ∼ χ .
In our previous study [48] 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. 1a of scalar mesons. Since ω P ∼ M P * 0 − M P + · · · , ω P 1 ∼ M P * 0 − M P 1 − (M P 1 − M P ) + · · · and M P * 0 > M P , M P * 0 < M P 1 , we find numerically that (ω ) contributes destructively to the mass shifts. Moreover, the self-energy of D * 0 (D * s0 ) is sensitive (insensitive) to the coupling g . Therefore, we can adjust the couplings h, g and the renormalization scale μ to get a mass degeneracy of D * s0 and D * 0 . We take the quark model of Godfrey et al. [7] as an illustration; it predicts the bare masses of 2480 and 2400 MeV for D *  near degeneracy of D * s0 and D * 0 in our previous work with g = 0. Hence, although the coupling g , which characterizes the interaction of P * 0 P 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 1a in Eq. (2.10), 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 but = 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 s1 = 2452 MeV and M D 1 = 2419 MeV are in agreement with the respective measured masses (see Table 1): 2459.5 ± 0.6 MeV and 2427 ± 36 MeV.
The coupling constant h at tree level can be extracted from the measured widths of D * 0 0 , D * ± 0 and D 0 1 . We find the values h = 0.60 ± 0.07 from D * 0 (2400) 0 , 0.514 ± 0.017 from D * 0 (2400) ± and 0.79 ± 0.17 from D 1 (2430) 0 . It is obvious that D * 0 (2400) ± gives the best determination of h as its mass and width measured recently by LHCb [2,3] are more accurate than D * 0 0 and D 0 1 . Our result of h = 0.508 agrees well with the LHCb experiment. As for the coupling g , its magnitude and even the sign relative to the coupling g in Eq. (2.4) are unknown. A recent lattice calculation yields h = 0.84(3)(2) and g = −0.122(8)(6) [62].
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 selfenergy 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 [7] as an illustration. If we use the quark model predictions from [8], near degeneracy can be implemented by employing h = 0.52 and g = 0.28 and similar given before. (iii) So far we have focused on the closeness of the masses of D * s0 and D * 0 0 . The measured mass of D * ± 0 by LHCb is greater than D * 0 0 and D * s0 by an amount of the order of 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 | 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 K * 0 , D s η , · · · for the self-energy diagram of D * s0 . (v) The calculated masses are sensitive to the choice of the renormalization scale . 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 Sect. 4.

Comparison with other work in HMChPT
Self-energy hadronic loop corrections to 0 + and 1 + heavy mesons have been discussed in the literature within the framework of HMChPT [31,[48][49][50][51][52][53][54]. Guo, Krewald and Meißner (GKM) [31] 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 (v · k), which has the following expression in [31]: where use of Eq. (2.24) has been made. Comparing with Eq. (2.21), we see that the m 2 coefficient in Eq. (3.1) should be replaced by ω 2 and a term (m 2 /ω)ln( 2 /m 2 ) in [· · · ] 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 [48]. Contributions from axial-vector heavy mesons to the self-energy of scalar mesons were also not considered in this work. Mehen and Springer [50] 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 · k) terms and neglect all mass splittings such as S , M P,S and residual masses such as a and˜ a in Eqs. (2.11) and (2.12). Then they set v · k = 0 inside the loop integral. As a consequence, ω = M ext − M 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 the scalar charmed mesons read and 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 Springer [50] and Alhakami [49], the bare mass M in the mass relation M = M+ 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 = M + 1 2 Re (see Eq. (2.3)) with 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 [50], Alhakami [49] has shown that the non-strange and strange scalar B mesons are nearly degenerate and the mass difference between B * 0 and B * s0 is ∼8 MeV. However, the renormalization scale is chosen to be = 317 MeV in [49], which is too small according to the spirit of ChPT.
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 B K threshold. As shown in Sect. 2, the closeness of the 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 detail in our previous work [48], 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 S and λ S 2 defined in Eqs. (2.7) and (2.9), respectively, are independent of the heavy quark flavor, where S measures the spinaveraged mass splitting between the scalar doublet (P 1 , P * 0 ) and the pseudoscalar doublet (P * , P) and λ S is the mass splitting between spin partners, namely P 1 and P * 0 , of the scalar doublet. From the data listed in Table 1, we are led to for the scalar cd meson and  Note that the parameter S can be decomposed in terms of the mass differences between 0 + and 0 − states and between 1 + and 1 − states defined by We have If we follow [93] assuming the heavy flavor independence of 0 and 1 , namely, we will obtain (in MeV) Our results for B * 0 0 and B 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 0 and 1 is probably too strong. Hence, we shall follow [71,72] assuming the heavy quark flavor independence of S and λ S 2 , (4.8) In this case, we find the results (in MeV) Predictions along this approach were first obtained in [71,72]. Evidently, mass degeneracy in the scalar D sector will repeat itself in the B system through heavy quark symmetry. Scrutiny of Table 4 shows that this degeneracy is respected only in few of the predictions such as [48,49,71,72,83,93]. Though the relations in Eq. (4.8) are valid in the heavy quark limit, they do receive 1/m Q and QCD corrections. The leading QCD correction to the relation λ is given by [94] Comparing with Eq. (4.9) 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 S . In heavy quark effective theory, it follows that (see [48] for details) (4.13) While the parameter λ H 1 is calculable, the other parameter, λ S 1 , is unknown. In [48] the 1/m Q correction is estimated to be (4.14) to be compared with the estimate of −50 ± 25 MeV in [50].
As pointed out in [48], so long as δ S is not large in magnitude, the 1/m Q correction to the relation Evidently, the mass degeneracy between B * s0 and B * 0 shown in Eq. (4.9) 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 for charmed mesons. Therefore, one will expect the mass difference relation also valid in the B sector. The mass difference between B s1 and B * s0 in Eq. (4.15) is independent of δ S ; its numerical value 48 ± 1 MeV indeed respects Eq. (4.17). 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 s1 − M B * s0 ∼20-30 MeV (see Table 4). (ii) In the charm sector we have the hierarchy It is natural to generalize this to the B sector: Indeed, the mass splittings of the j = 1 2 doublet in strange and non-strange sectors satisfy the relation (see [48]) is not respected by some of the models listed in Table 4. (iii) The difference 15 MeV between the central values of M B s1 and M B 1 inferred from Eq. (4.15) is much smaller than the quark model expectation of 60-100 MeV (cf. Table 4). 3 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 ± 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. (4.15) 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 s1 are below the B K and B * K thresholds, respectively, their widths are expected to be very narrow with the isospin-violating strong decays into B s π 0 and B * s π 0 . Experimentally, it will be even more difficult to identify the non-strange B * 0 and B 1 mesons owing to their broad widths.

Conclusions
The qq 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 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 and the renormalization scale . This in turn implies the mass similarity of B * s0 and B * 0 0 . The P * 0 P 1 interaction with the Goldstone boson characterized by the strong coupling g 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 (s)1 mesons (i.e. Eq. (4.15)) 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 ± 1) MeV + δ S , M B s1 = (5763 ± 1) MeV + δ S with δ S being 1/m Q corrections. We found that the mass difference of 48 MeV between B s1 and B * s0 is larger than the result of 20 ∼ 30 MeV predicted by the 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.
• The current world-averaged masses of neutral and charged D * 0 (2400) mesons are not the same; they differ by an amount of 33±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.