The assignments of the B s mesons within the screened potential model and 3 P 0 model

We investigate the mass spectrum and the decay properties of the B s mesons within the screened nonrelativistic quark model and the 3 P 0 model. Our results suggest that the B sJ ( 6064 ) and B sJ ( 6114 ) states, as the ﬁrst solution of the recently LHCb measurements, could be explained as the B s ( 1 3 D 3 ) and B s ( 1 3 D 1 ) , respectively. In addition, the B sJ ( 6109 ) and B sJ ( 6158 ) states, as the second solution of the LHCb measurements, could be explained as the B (cid:2) s 2 ( 1 D ) and B s 1 ( 2 P ) , respectively. Meanwhile, the B s 1 ( 5830 ) could be interpreted as the candidate of the B s 1 ( 1 P ) . We also calculated the decay properties of the other excited B s mesons with the predicted masses, which should be helpful for the experimental searching in future.


I. INTRODUCTION
Recently, the LHCb Collaboration observed an excess structure 300 MeV above the B ± K ∓ threshold in the B ± K ∓ mass spectrum in the proton-proton collisions, which could be described by a two-peak hypothesis [1].
However, if a decay through B * ± K ∓ with a missing photon from the B * ± → B ± γ decay is assumed, the masses * Electronic address: haowei2020@itp.ac.cn † Electronic address: ylu@ucas.ac.cn ‡ Electronic address: wangen@zzu.edu.cn and widths will shift to be, B sJ (6158) : M = 6158 ± 4 ± 5 MeV, Γ = 72 ± 18 ± 25 MeV.(8) These observations of the LHCb have enriched the bottom-strange spectrum.We have tabulated the experimental information of all the B s mesons in Table I.
Although these two states were suggested to be the Dwave orbital excited bottom-strange mesons, their masses are significantly lower than the quark model predictions [2][3][4][5].There have been some theoretical works to study these states [6][7][8][9] (A recent review about these two states can be found in Ref. [10]).Based on the nonrelativistic linear potential model, Ref.For the B sJ (6064), there are also two possible interpre-tations in Ref. [7].According to the first interpretation, the narrow structure around 6064 MeV is mainly caused by the B sJ (6109) resonance, which is regarded as the In addition, authors of Ref. [9] have studied the B and B s mesons using the heavy quark effective theory, which suggests that the B sJ (6064) could be the candidate of the B s (2 Besides the conventional B s mesons explanation, the two states are also regarded as bsq q tetraquark states in Ref. [11], or as the BK * molecular state with quantum numbers (I)J P = (0)1 + in Ref. [12].Thus, one can find that the natures of these two states are still in debate, and more efforts are needed to shed light on their internal structures.
As we known, although the quenched quark models have obtained lots of success in the last decades, the effects of the sea quarks and gluons interactions are not taken into account.The unquenched models, considering all kinds of additional effects, have been developed, and widely used to describe the hadron spectra, such as the coupled channel model [13][14][15][16][17][18][19] and the screened potential model [20][21][22][23][24][25].In the quenched potential model, the potentials mainly contain a coulomb term at short distances and the linear confining interaction at large distances.However this is not appreciate in the large mass range, since the linear potential, which is expected to be dominant in large mass region, will be screened or softened by the vacuum polarization effects of dynamical fermions [26,27], i.e., the unquenched effects reflecting the sea quarks or gluons contributions to some extent.
Clearly, the unquenched effects can lead to important influence for higher radial and orbital excited hadrons, which means that the predicted masses of the higher ex-cited states will be smaller than the ones of the general liner potential models.Comparing with the coupled channel model, the screened potential model is simpler, and have been successfully used to describe the spectra of the charmed-strange meson [21,28], charmed meson [20], excited ρ mesons [29,30], bottom mesons [25], charmonium [23,24,27], and bottomonium [22,31].B sJ (6158) 6158 ± 4 ± 5 72 ± 18 ± 25 0(? ? ) [1] In this paper, we use the screened nonrelativistic quark This article is organized as follows.In Sec.II, we give a brief introduction about the screened nonrelativestic quark model and the 3 P 0 model.In Sec.III, the numerical results and the discussions are presented.Finally, the summary is given in Sec.IV.

II. THEORETICAL MODELS A. Screened nonrelativistic quark model
The nonrelativistic quark model mainly includes the confinement term , the spin-dependent term, and the oneloop correction for the spin-dependent terms [33][34][35], and the Hamiltonian for a q q meson system is defined as [2,36], where H 0 is the zeroth-order Hamiltonian, H sd is the spin-dependent Hamiltonian, and C q q is a constant, which will be fixed to experimental data.The H 0 can be compressed as, where the confinement interaction includes the standard Coulomb potential −4α s /3r and the linear scalar potential br.The last term is the hyperfine interaction that could be treated nonperturbatively.p is quark momentum in the system of q q meson, r = | r | is the q q separation, M r = 2m q m q /(m q + m q ), m q (m q ) and S q (S q ) are the reduced mass of the q q system, the mass and spin of the constituent quark q (antiquark q), respectively.
The spin-dependent term H sd is, with where S ± = S q ± S q , L is the relative orbital angular momentum of the q q system.We take Euler constant The screening effects are introduced by the following replacement, where V scr (r) behaves like br at short distances and constant b/β at large distance [20,21], β is the parameter which controls the power of the screening effects.One can find that the screened potential approximates the liner potential for a small distance r, and will be softened for large distance r.Since the distance between the quarks in the excited bottom-strange mesons is larger than the one of the ground bottom-strange meson, it is expected that the spectrum of the screened potential model is more sensitive for the excited bottom-strange mesons.
The spin-orbit term H sd can be decomposed into symmetric part H sym and antisymmetric part H anti , which can be expressed as [2] The antisymmetric part H anti gives rise to the the spinorbit mixing of the heavy-light mesons with different total spins but with the same total angular momentum, such as B s (n 3 L L ) and B s (n 1 L L ).Hence, the mixing of the two physical states B sL (nL) and B ′ sL (nL) can be expressed as, where the θ nL is the mixing angles.
With above formalism, one can solve the Schrödinger equation with Hamiltonian H of Eq. ( 9) to obtain the mass spectrum and the meson wave functions, where the wave functions will be used to calculate the strong decays of excited bottom-strange mesons in the 3 P 0 model.

B. The 3 P0 model
The 3 P 0 model was proposed by Micu [37] and further developed by Le Yaouanc [38][39][40][41], and it has been widely used to calculate the OZI allowed decay processes [2, 20-25, 36, 42-52].In this model, the meson decay occurs through the regroupment between the q q of the initial meson and the another q q pair created from vacuum with the quantum numbers J P C = 0 ++ .The transition operator T of the decay A → BC in the 3 P 0 model is given by where ) is solid harmonic polynomial in the momentum space of the created quarkantiquark pair.χ 34  1,−m , φ 34 0 , and ω 34 0 are the spin, flavor, and color wave functions, respectively.The paramtere γ is the quark pair creation strength parameter for uū and d d pairs, and for ss we take γ ss = γ mu ms [41].The parameter γ can be determined by fitting to the experimental data.The partial wave amplitude M LS (P ) of the decay A → BC is be given by, where M MJ A MJ B MJ C (P ) is the helicity amplitude, Here, |A , |B , and |C denote the mock meson states.
Then, the decay width Γ(A → BC) can be expressed as where

III. RESULTS AND DISCUSSIONS
In the calculation, the screened parameter β = 0.025 GeV and the constant C q q = 0.1035 GeV were obtained by fitting the well-known states B s (1 1 S 0 ), B * s (1 3 S 1 ), and B * s2 (5840)(1 3 P 2 ).The other parameters are taken from Ref. [35].The 3 P 0 model parameter γ = 0.354 is obtained by fitting the total decay width of the B * s2 (5840), which is regarded as the B s (1 3 P 2 ).With these parameters, the predicted ratio of the B * s2 (5840) decay modes, which is consistent with LHCb experimental data of 0.093 ± 0.013 ± 0.012 [53].
The predicted masses of the B s mesons are listed in  As we discussed in the introduction, if a decay through B * ± K ∓ with a missing photon from the B * ± → B ± γ decay is assumed, the masses and widths of the two states observed by LHCb will shift, and the two states are named as B sJ (6109) and B sJ (6158) [1].In this case, the mass and width of the B sJ (6109) are close to the predicted mass (6132 MeV) and width (43 MeV) of B ′ s2 (1D), which implies that B sJ (6109) could be regarded as the B ′ s2 (1D) state.On the other hand, the mass and width of the B sJ (6158) are close to the predicted mass (6194 MeV) and width (75 MeV) of B s1 (2P ), respectively, which supports the assignment of the B sJ (6158) as the B s1 (2P ) state.It should be stressed that the two solutions could be not distinguished according to the present LHCb measurements.Thus, the future precise measurements of their masses, widths, and the quantum numbers of the spin-parity would be helpful to shed light on this problem, and deepen our understanding the spectra of the bottom-strange mesons.According to the predicted mass spectrum of Table II, the predicted mass 5840 MeV of the B ′ s1 (1P ) is very close to the one of the B s1 (5830).With the mixing angle of −55.8 • , the decay widths of the B ′ s1 (1P ) are calculated, as shown in Table III, and the predicted total decay width is very small, in good agreement with the experimental value 0.5 ± 0.3 ± 0.3 MeV of the B s1 (5830).We show the total decay widths of the B ′ s1 (1P ) versus the mixing angle in Fig. 1, where one can find the total decay width is still consistent with the experimental data with the mixing angle in the range of −60 • ∼ −50 • .
In addition, we also predict the decay widths of the other excited B s mesons with the predicted masses of Table II.For the B s (2 3 P 0 ), its total decal width is predicted to be 66.8  In addition, the B sJ (6109) and B sJ (6158) states, as the second solution of the LHCb measurements, could be explained as the B ′ s2 (1D) and B s1 (2P ), respectively.Since those two solutions could be not distinguished according to the present LHCb measurements, thus the future precise measurements of their masses, widths, and the quantum numbers of the spin-parity would be helpful to shed light on this problem, and deepen our understanding the spectra of the bottom-strange mesons.
In addition, we suggest that the state B s1 (5830) could be explained as the candidate of the B ′ s1 (1P ) state.With the predicted masses of the other excited B s states, we also predict their decay widths as well as the dominant   It should be stressed that there are already many theoretical studies about the family of the bottom-strange mesons.Comparing with those works, we have adopted the nonrelativistic quark model by taking into account the screening effects, which play an important role for the higher radial and orbital excited mesons, and our results could give better descriptions for all the existed bottom-strange mesons.
[7] explained the B sJ (6114) as the B s (2 3 S 1 ) − B s (1 3 D 1 ) mixing state, and obtained the mass M = 6114 MeV and width Γ = 95 ± 15 MeV with the mixing angle θ = −(45 ± 16) • ,which is supported by the results based on the screened potential model[8].In addition, Ref.[7] also supports the interpretation of the B sJ (6114) as a pure B s (1 3 D 1 ) state if there is small mixing between B s (2 3 S 1 ) and B s (1 3 D 1 ).
mixing state.For the second interpretation, the B sJ (6064) could be explained as a pure B s (1 3 D 3 ) with the predicted mass M = 6067 MeV and width Γ = 13 MeV.

model and the 3 P
0 model to study the spectrum and the strong decay properties of the B s mesons, and also to explore the possible assignments of the two resonances recently observed by the LHCb Collaboration.
and M C are the masses of the mesons A, B, and C, respectively.The spatial wave functions of the mesons in the 3 P 0 model are obtained by solving the Schrödinger equation in Eq. (9).
) and B s (1 3 D 3 ) are predicted to be 6061 MeV and 6117 MeV, in good agreement with the experimental results of the B sJ (6064) (6063.5 ± 1.2 ± 0.8 MeV) and B sJ (6114) (6114 ± 3 ± 5 MeV), respectively, which indicates that the two states could be the possible candidates of the B s (1 3 D 1 ) and B s (1 3 D 3 ) .Of course, only the mass information is not enough to establish these assignments.We also calculate the strong decay widths of the B s mesons, as shown in Table III.One can find that the predicted width of 23 MeV for B s (1 3 D 3 ) is in good agreement with the measured width 26 ± 4 ± 4 MeV of the B sJ (6064), which supports the B s (1 3 D 3 ) assignment of the B sJ (6064).In addition, the predicted width for B s (1 3 D 1 ) is 127 MeV, reasonably consistent with the one of B sJ (6114) if taking into account the large experimental uncertainties.Thus, the B sJ (6114) could be explained as the B s (1 3 D 1 ) state, and more precise measurements will be helpful to pin down this assignment.

FIG. 1 :
FIG. 1: Total decay widths of the Bs1(5830) as the 1P ′ versus the mixing angle.The blue band denotes the uncertainties of the experimental data.
MeV, and the dominant decay mode is BK with the branching fraction 92%.The total decay width of B ′ s1 (2P ) is 113.9MeV, while the dominant decay modes are BK * /B * K * .The predicted width and the dominant decay mode of the B s (2 3 P 2 ) are 171.6MeV and B * K * , respectively.In addition, the decay widths of the B s2 (1D) are 127.3MeV, and the dominant decay mode is B * K. Our results should be helpful to search for them in experiments, such as LHC.IV.SUMMARY Recently, the LHCb Collaboration has observed two resonances B sJ (6064) and B sJ (6114) assuming they decay directly to the BK final states.However, their masses and widths will shift if a decay through B * K with a missing photon from the B * decay, and those two resonances are named as B sJ (6109) and B sJ (6158).Motivated by the recently LHCb measurements, in this paper we calculate the spectrum of the B s mesons within the screened nonrelativistic quark model, and also investigate the strong decay properties of these mesons with the 3 P 0 model.By comparing with the experimental data, it is found that the B sJ (6064) and B sJ (6114) states, as the first solution of the recently LHCb measurements, could be explained as the B s (1 3 D 3 ) and B s (1 3 D 1 ), respectively.
3 S 1 ) with the predicted mass 6033.0 ± 2.4 MeV and width 170 ± 1.5 MeV, and the B sJ (6114) could be the candidate of the B s (1 3 D 3 ) with predicted mass 6247.0±2.4MeV and width 82.0±1.0MeV.However, the predicted width for B sJ (6064) and mass for B sJ (6114) are larger than experimental values.

Table II
, where we also show the predictions of other theoretical works for comparison.The B s [B s (1 1 S 0 )], B * s [B s (1 3 S 1 )], and B * s2 (5840)[B s (1 3 P 2 )] can be well described in the spectrum.The masses of the B s (1 3 D 1

TABLE III :
The decay widths of the Bs mesons in units of MeV.The label '−' means that the channel is forbidden or there is no experimental information.