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 $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.

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±1) MeV +δ∆ S , M B ′ s1 = (5763 ± 1) MeV + δ∆ S with δ∆ 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 ∼ 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. 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 .

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: ρ, ω, K * , φ. Since the constituent strange quark is heavier than 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 low-lying 3 P 0 qq states simultaneously. If the light scalar nonet is identified with the Pwave 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.
• 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 I). 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 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], while the neutral one is from BaBar [3], Belle [4] and FOCUS [5]. 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 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 [8], has the quantum numbers J P C = 1 ++ [9]. X(3872) cannot be a pure charmonium as it cannot be identified as χ 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 [10]. 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 [11,12].
The above discussions suggest that X(3872) is most likely an admixture of the S-wave DD * molecule and the P -wave charmonium as first advocated in [11] |X(3872) = c 1 |cc P−wave + c 2 |D 0 D * 0 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. [18,19]). 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). 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 [6,7]. It was first stressed and proposed in [20] 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 [21]). This conjecture was realized in both QCD sum rule [22,23] and lattice [24][25][26] calculations. For example, when the contribution from the DK continuum is other. As a result, the world averaged masses for the neutral and charged D * 0 (2400) shown in Table I 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 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). I: 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 [6,7]. "Large" means a broad width of order 100 MeV, while "small" implies a narrow width of order 10 MeV. included in QCD sum rules, it has been shown that this effect will significantly lower the mass of the D * s0 state [22]. Recent lattice calculations using cs, DK and D * K interpolating fields show the existence of D * s0 (2317) below the DK threshold [24] and D ′ s1 (2460) below the D * K threshold [25]. 2 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 [28]. 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) [29][30][31][32]. 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 [33][34][35]. Both f 0 (980) and a 0 (980) have the strong couple channel KK. They are often viewed as KK molecules, which accounts for their near degeneracy with 2m K . Schematically, the self-energy KK 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π, KK, and KK, respectively [36]. Alternatively, it is well known that the tetraquark picture originally advocated by Jaffe [37] provides a simple solution 2 A recent lattice calculation with N f = 2 + 1 + 1 optimal domain-wall fermions [27] yields a mass of 2317 ± 15 ± 5 MeV for D * s0 and 2463 ± 13 ± 9 MeV for D ′ s1 (2460).
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 4-quark 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 [38]. The fall-apart decays f 0 (500) → ππ, K * 0 (800) → πK, and f 0 (980), a 0 (980) → KK 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 KK threshold.
In [39] 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 [40] who followed the framework of [41] 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 [40], 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 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 ′+ heavy mesons have been considered in the literature [28,[39][40][41][42][43][44][45]. 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 [46][47][48]. The relevant Lagrangian is [49] 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): 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 Σ = exp(i √ 2φ/f π ) with f π = 93 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 A = i 2 (ξ † ∂ µ ξ − ξ∂ µ ξ † ). In Eq. (2.4), 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 FIG. 1: Self-energy contributions to scalar and axial-vector heavy mesons.
so that 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 [49,50] 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 [41]. The masses of heavy mesons can be expressed as where ∆ a and∆ a denote the residual mass contributions to odd-and even-parity mesons, respectively. and 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 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 channels D 0 K + , D + K 0 , D + s η and 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 mass shifts from above equations. Assuming that the bare mass M is the one obtained in the quark model, then from Eq. (2.10) we have For the self-energy of the axial-vector meson, we consider D ′ s1 as an illustration which receives contributions from D * 0 K + , D * + K 0 , D * + s η and D * 0 0 K + , D * + 0 K 0 , D * + s0 η intermediate states (see Fig.  1 where and [40,41] Note that the function F (−m/ω) can be recast to 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 the conventional practice, it is often to choose Λ ∼ Λ χ , the chiral symmetry breaking scale of order 1 GeV, to get numerical estimates of chiral loop effects. However, as pointed out in [39], 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 higher-order heavy quark expansion needs to be taken into account to justify the use of Λ ∼ Λ χ .
In our previous study [39] 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 , we find numerically that Π ′ (ω ′ ) contributes destructively to the mass shifts. Moreover, the selfenergy 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 . Take the quark model of [6] 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 ′ = 0. 25 Table I): 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] 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) [53].
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 [6] as an illustration. If we use the quark model predictions from [7], 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 D * s0 and D * 0 0 masses . The measured mass of D * ± 0 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 ′ | 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 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 Sec. IV below.

III. COMPARISON WITH OTHER WORKS IN HMCHPT
Self-energy hadronic loop corrections to 0 + and 1 ′+ heavy mesons have been discussed in the literature within the framework of HMChPT [28,[39][40][41][42][43][44][45]. Guo, Krewald and Meißner (GKM) [28] 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 [28] 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 [39]. Contributions from axial-vector heavy mesons to the self-energy of scalar mesons were also not considered in this work. Mehen and Springer [41] 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  Table II except that the propagator i/(2v · k) is used for all heavy meson states and v · k is set to be zero inside the loop integral calculation.

Meson
Bare and The results of calculations are summarized in Table III. It is clear that while near degeneracy is achieved for scalar charmed mesons, it is not the case for scalar B mesons.
It is clear from Table IV 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 [39], 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 spin-averaged 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 I,  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 will obtain (in MeV) M B * s0 = 5716.2 ± 0.6, M B ′ s1 = 5762.7 ± 1.8 . Our results for B * 0 0 and B ′0 1 masses are slightly different from those shown in the last row of Table  IV. 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 [62] to assume 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 [62]. Evidently, mass degeneracy in the scalar D sector will repeat itself in the B system through heavy quark symmetry. A scrutiny of Table IV shows that this degeneracy is respected only in few of the predictions such as [39], [40], [62], [72], and [81].
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 λ S(bq) 2 = λ S(cq) 2 is given by [82] 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 [39] for details) While the parameter λ H 1 is calculable, the other parameter λ S 1 is unknown. In [39] the 1/m Q correction is estimated to be (4.14) to be compared with the estimate of −50 ± 25 MeV in [41].
As pointed out in [39], so long as δ∆ S is not large in magnitude, the 1/m Q correction to the relation ∆ S [see Eq. (4.13)], to a very good approximation, amounts to lowering the masses of B * (s)0 and B ′ (s)1 by an equal amount of |δ∆ S |: 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 the relation (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 Table IV). (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 nonstrange sectors satisfy the relation (see [39]) is not respected by some of the models listed in Table  IV. (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 IV). 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 BK 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.

V. 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.