Electric transitions of the charmed-strange mesons in a relativistic quark model

In the present work, we adopt a relativistic constituent quark model to depict the charmed strange meson spectroscopy, in which $D_{s0}(2317)$ and $D_{s1}(2460)$ are considered as the $1^3P_0$ and $1P_1^\prime$ charmed strange mesons, respectively. By using the wave function obtained from the relativistic quark model, we further investigate the electric transitions between charmed strange mesons. We find the long wave length approximation is reasonable for the charmed strange meson radiative decay by comparing the results with different approximations. The estimated partial widths are all safely under the upper limits of the experimental data. Moreover, we find the branching ratio of $D_{s1}(2536) \to D_s^\ast \gamma/D_s \gamma$ are large enough to be detected, which could be searched by further experiments in Belle II and LHCb.


I. INTRODUCTION
Charmed-strange meson is one of important members of meson family. The ground S −wave states, D s and D * s , were first observed more than 40 years ago in the e + e − annihilation process by DASP Collaboration [1] 1 . Later, in theνN collisions, a new state, D s1 (2536), was observed in the D * s γ invariant mass spectrum [2], which could be a ground P−wave state. The second P− wave states D * s2 (2537) was observed in the DK and D * K modes in the B meson decay processes by CLEO Collabortaion [3].
Nearly ten years later after the observation of D s2 (2537), the rest two ground P−wave state candidates, D * s0 (2317) and D s1 (2460), were discovered [4,5]. The former one was first observed in the D s π 0 invariant mass spectrum of B decay process by BaBar Collaboration [4] and the later one was reported in a similar process but in the D * s π 0 invariant mass spectrum by Belle Collaboration [5]. These two states are particular interesting since their observed masses are much lower than the quark model expectation [6] and several tens MeV below the threshold of DK and D * K, respectively. Thus, these two states were ever considered as DK and D * K molecular states due to their particular properties [7][8][9][10][11][12][13][14]. However, considering the coupled channel effects and the fact that there are no additional states around quark model predicted masses, the authors in Refs. [15][16][17][18][19][20][21][22][23] assigned these two states as P− wave charmed-strange mesons. In this case, the ground P−wave charmed strange mesons are established.
In 2006, the BaBar Collaboration reported two new charmed-strange meson in the DK invariant mass spectrum of B meson decay [24]. The narrow one is D sJ (2860) and the broader one is D * s1 (2700). The theoretical estimation indicate that the D * s1 (2700) could be a good candidate of 2 3 S 1 state [23,[25][26][27]. In 2014, the LHCb Collaboration analyzed theDK invariant mass spectrum of B 0 s →DK − π + process D s (1968) FIG. 1: The history of charmed-strange meson discovery [1-5, 24, 28]. Here, the masses of the charmed-mesons are taken from the Review of Particle Physics [35]. and find the structure around 2860 MeV announced by BaBar Collaboration consist of two particle with spin-1 and spin-3 [28], which were named as D * s1 (2860) and D * s3 (2860). As indicated in Refs. [29][30][31][32], these two states could be good candidates of D−wave charmed-strange mesons 1 3 D 1 and 1 3 D 3 , respectively. To data, the observed heaviest charmed strange meson is D sJ (3040), which was discovered in D * K invariant mass spectrum of B decay process by BaBar Collaboration [33], which can be assigned as 2P 1 state as indicated in Refs. [23,34].
In Fig. 1, we present the history of the observation of charmed-strange mesons, where we find most excited charmed-strange mesons were observed during the year of 2003-2014. Moreover, from the discovery history one can find most of the charmed-strange mesons are firstly observed in the bottom or bottom-strange meson decays. With the running of Belle II and LHCb, more excited charmed strange meson are expected to be discovered in the bottom or bottom-strange meson decays, which will make the charm-strange family abundant.
Besides the observed resonance parameters, i.e., the mass and width, the decay behaviors of the observed states are also crucial to understand their inner structures. In particular, the electromagnetic transitions can be well described by Quantum electrodynamics in the quark level, which is unlike the non- Possibly Seen perturbative strong interactions in the hadron energy. Thus, the electromagnetic transitions could reflect the inner structure in a more comprehensive manner. On the experimental side, there are some experimental measurements for the radiative transitions between charmed mesons, the corresponding experimental information are collected in Table I. Thus the investigation of the radiative decays of charmed strange mesons are interesting and necessary.
In the present work, the charmed meson spectroscopy is depicted by a relativistic quark model [20], where the masses of the charmed strange mesons are well reproduced. With the wave functions estimated in the quark model, we estimate the electric transitions between the charmed-strange mesons, which could not only further test the relativistic quark model by comparing the theoretical estimations with the experimental measurements but also provide some useful predictions.
This work is organized as follows, in Section II, we present a short review of the relativistic quark model, by which the mass spectroscopy of charmed strange mesons are estimated. In Section III, we present the formula of the electric transitions between charmed-strange mesons, and the numerical results and discussions are presented in Section IV. A short summary is given in Section V.

II. REVIEW OF MASS SPECTROSCOPY OF CHARMED STRANGE MESONS
Relativistic quark model is usually adopted to depict the mass spectroscopy of hadrons since the non-perturbative properties of QCD in hadron energy. Such kind of quark model was proposed to investigate the meson spectroscopies systematically in 1985 by Godfrey and Isgur [6].
In this model, the mass spectroscopy and wave functions of the mesons can be determined by solving the relativistic Schrödinger equation, where the spin independent Hamiltonian can be, where r, p are the coordinates and the momentum of quark in the center-of-mass frame, respectively. m 1 and m 2 are the masses of the quark and the antiquark, respectively. V(r) is the effective spin-independent potential between the quark and the antiquark, including a Coulomb term and a linear confin-ing term [6], which is As for the spin-dependent part H , it includes the spin-spin interaction and spin-orbital interactions, which is, and the concrete form of spin-spin and spin-orbital interactions are where the functions f (r), g(r), h 1 (r), h 2 (r) can be found in Ref. [20]. With these spin dependent terms, the S -D mixings and spin-singlet and spin-triplet mixings have been included. In this model, the mass spectroscopy of the charmed strange mesons can be well reproduced, thus, in the present work, we adopt the same model parameters as those in Ref. [20] to investigate the electric radiative decays of the charmed strange mesons. Before the estimations of the radiative decays, we present the mass spectroscopy of the charmedstrange mesons in Table II, where the theoretical estimations from Ref. [34,36] and experimental data [35] are also listed for comparison. The theoretical estimated mass of 1 3 P 0 and 1P 1 states are 2317 and 2425 MeV, respectively, which are more consistent with the experimental measurements comparing to other works [34,36]. In the quark level, the quark-photon electromagnetic interaction can be written as where ψ j and e j represent the j−th quark fields and its charges in the charmed-strange meson, respectively. The k is three momentum of the emitted photon. After performing some algebra estimation as shown in Appendix A, the amplitude of the electromagnetic transition can be expressed, where |i and | f are the wave functions of initial and final states, respectively. ω is the energy of the emitted photon.
In the present work, we mainly focus on the electric transition processes, and the helicity amplitude is where the initial and final hadron wave functions can be estimated by the relativistic quark model. In the estimation, we can choose the photon momentum direction along the z axis, i.e., k = kẑ, and the photon polarization vector is in right-hand form, which is In this case, e −ik·r j can be expanded as, then the helicity amplitude for the angular momentum l can be [37], and then the decay width of the electric transition between Qq can be estimated as 3 where |n 2S +1 L J and |n 2S +1 L J represent the final and initial states, respectively. C f i is a coefficient related to involved states, which is Considering the lowest order of the electric transition, the terms related to j 2 (kr) can be ignored, then the electric transition width can be In the literatures, the zeroth order spherical Bessel function j 0 (kr) is usually expanded as j 0 (kr) = 1 + O(x 2 ), keeping the lowest order, one can get the partial width as where e M = (eQmq − e q m Q )/(m q + m Q ). The approximation in above formula corresponds to the long wave length approximation, where e ik·r ∼ 1.

IV. NUMERICAL RESULTS AND DISCUSSIONS
As indicated in the last section, the partial widths of electric transition can be estimated with different approximations, hereafter, we use Mode I, Mode II and Mode III to refer the estimations with Eqs.(10), (12) and (13), respectively and further check the reliability of different approximations. With the wave functions estimated from the relativistic quark model and the formula in above section, we can get the partial widths of the electric transitions, which are listed in Tables III-VI. In Table III, we present the electric transitions for P → S γ processes, where P and S indicate the P− and S − wave charmed-strange mesons, respectively. In addition, we also listed the theoretical results from other groups [38][39][40][41] for comparison. From the table, one can find the estimation from different approximations are almost the same, which indicates that the approximation from Mode I to Mode III are still reliable and long wave length approximation in the considered electric transitions of charmed-strange mesons is reasonable. Our estimation indicates that most of our results are of the same order as those in Refs. [38][39][40][41]. In particular, we find the partial widths of D s (2 1 S 0 ) → D s1 (1P 1 )γ and D s (2 1 S 0 ) → D s1 (1P 1 )γ from different literatures are very different. Our estimation shows that the partial width of D s (2 1 S 0 ) → D s1 (1P 1 )γ is 0.07 keV, which is of same order as those in Ref. [38,39], but the estimation in Ref. [40,41] are about two order larger than the present estimation. As for D s (2 1 S 0 ) → D s1 (1P 1 )γ, our estimation is of the same order as the one in Ref. [40], but much smaller than those in Ref. [41]. The estimations in the present work and in Refs. [38][39][40] are all based on relativistic quark model. But it should be notice that the estimated mass spectroscopy in Refs. [38,39] are similar to the present one, where the masses of D s0 (2317) and D s1 (2460) were well reproduced. Thus the meson wave functions should be similar and so do the electric transition widths. III: Electric transition widths for S → Pγ processes, where P and S are P− and S −wave charmed-strange mesons, respectively. For comparison, we also present the theoretical estimations from Refs. [38][39][40][41]. The results of Ref. [40] are estimated with the mixing angle θ 1P = −38 • [42,43] Initial Final Decay Width (keV) Mode I Mode II Mode III Ref. [38] Ref. [39] Ref. [40] Ref. [ [38] Ref. [39] Ref. [44] Ref. [45] Ref. [40] Ref. [ As for Ref. [41], the estimated mass spectroscopy are much different with the present one and the estimated masses of D s0 (2317) and D s1 (2460) are far above the measured values, then the meson wave functions and electric transition widths are much different. As for Ref. [40], the estimations are based on heavy quark limit, which should be more reliable for bottom mesons.
As for the radiative decay of 3S states, we find the that Γ(D * s (3 3 S 1 ) → D * s0 (2 3 P 0 )γ) = (13 ∼ 14) keV. As for D * s (3 3 S 1 ), it is far above the threshold of DK, and it dominantly decay into a charmed meson and a strange meson, and its total width are estimated to be around 100 MeV [23], and with such a large width, the branching ratio of D * s (3 3 S 1 ) → D * s0 (2 3 P 0 )γ is of order 10 −4 .
In Table IV, we present our estimated widths for P → S γ. Our estimation indicates that the partial widths for 1P → 1S γ vary from several keV to 10 keV, which is consistent with those in the previous literatures [38-41, 44, 45]. The partial width of D s0 (1 3 P 0 ) → D * s γ is estimated to be around 2 keV. The measured upper limits of Γ D s0 (1 3 P 0 ) and B(D s0 (1 3 P 0 ) → D * s γ) are 3.5 MeV and 6%, respectively. Thus, the upper limit of the partial width of D s0 (1 3 P 0 ) → D * s γ is 210 keV, which indicates our estimation is safely under the upper limit of the experimental values.
As for D s1 (2460), the widths of D s γ and D * s γ modes are 3.61 and 4.79 keV, respectively, which are both safely under the upper limits of the experimental values. Moreover, from our estimation, we find that the partial width of D * s γ mode is a bit larger than the one of D s γ, which is similar to those in Refs. [38][39][40][41], but different with the experimental measurements, which are B(D s1 (2460) → D s γ) = (18 ± 4)% and B(D s1 (2460) → D * s γ) < 8%. It should be notice that the D s1 (2460) state has the components with both S = 0 (corresponding to 1 P 1 state) and S = 1 (corresponding to 3 P 1 state), while in the electric transitions, the spin of the initial and final states should be the same, thus, the electric transitions involves D s1 (nP 1 ) and D s1 (nP 1 ) states are sensitive to the spin singlet and triplet mixing.
Our estimation for P → Dγ and D → Pγ process are listed in Table V and VI. As for P → Dγ processes, the largest one is D * s2 (2 3 P 2 ) → D s3 (1 3 D 3 ), which are 3.22 keV. As for D → Pγ processes, the partial widths of D * s (1 3 D 1 ) → D * s0 (1 3 P 0 )γ, D * s (1 3 D 3 ) → D * s0 (1 3 P 2 )γ and D s2 (2D 2 ) → D s1 (2P 1 )γ processes are greater than 10 keV. These highly excited states are far above the threshold of DK and D * K, and they dominantly decay into a charmed meson and a strange meson, their width should be of order 100 MeV. Thus the branching ratios of these radiative decays should be of order of 10 −4 .

V. SUMMARY
The radiative decay is one of important decay modes of charmed strange mesons, especially for the low lying charmed strange mesons. In the present work, we adopt a relativistic constituent quark model to depict the mass spectroscopy of the charmed meson, in which D s0 (2317) and D s1 (2460) are considered as 1 3 P 0 and 1P 1 charmed strange mesons, respectively, while D sJ (3040) is assigned as D s1 (2P 1 ) states.
With the wave function estimated by the relativistic quark model, we evaluate the electric transitions between the charmed strange mesons. By comparing the transition widths obtained with different approximations, we find that the long wave length approximation is reasonable for most cases of the electric transitions between charmed-strange mesons. Our estimation indicates that the partial widths of D s0 (1 3 P 0 ) → D * s γ, D s1 (1P 1 ) → D * s γ and D s1 (1P 1 ) → D s γ are all safely under the upper limits of the experimental data. As for D s1 (1P 1 ) → D s γ and D s1 (1P 1 ) → D * s γ, our estimation find that the branching ratios of these processes are large enough to be detected, which could be searched in further experiments in Belle II and LHCb. As for P → Dγ and D → Pγ processes, the width of some channels can reach up to 10 keV, which may be tested by further experimental measurements.
By replacing the quark fieldψ with ψ † , one can use matrix α instead of the γ matrix in Eq. (5). Then the electromagnetic transition matrix elements for a radiative decay process becomes, Considering the fact that the involve mesons are composite systems and the relativistic Hamiltonian is, we have the following identity, Then, the electromagnetic transition matrix can be expressed as, with h e = j e j r j · (1 − α j ·k)e −ik·r j , g e = j e j r j · e −ik·r j .
E i , E f and ω γ in Eq. (A4) are the energies of the initial meson, the final meson and the emitted photon, respectively. Thus, E i − E f − ω γ ≡ 0 due to the conservation of energy. Thus, one has, M = −iω γ f |h e | i (A6) Following the procedures used in Refs. [46,47], one can get the non-relativistic expansion of h e , which is, h e j e j r j · − e j 2m j σ j · ×k e −ik·r j , where the first and the second terms are corresponding to electric and magnetic transitions, respectively.