Mass spectroscopy and strong decays of excited open charm D J mesons using relativistic Dirac formalism

.


Introduction
During the past decade, numerous excited charmed states have been observed by the experimental groups L HCb and B AB AR [1][2][3]. These achievements provoked great interest in experimental as well as theoretical studies of charmed mesons. In recent past L HCb has employed the Dalitz Plot (DP) technique to analyse the contributing amplitudes in decay channel B − → D + π − π −, where charmed states were reconstructed through D + → K − π + π + decay process. This analysis was based on data collected by the L HCb detector during 2011 and 2012 when the pp collision center of mass energy was 7 TeV and 8 TeV, respectively. Their a e-mail: bhumispandya@gmail.com (corresponding author) b e-mail: mnshah09@gmail.com c e-mail: p.c.vinodkumar@gmail.com study summarises the resonant contribution coming from D * 2 (2460) 0 , D * 1 (2680) 0 , D * 3 (2760) 0 and D * 2 (3000) 0 states [1]. Their measured masses and widths are listed in Table 1.
In 2013, L HCb has studied the D + π − , D 0 π + and D * + π − channel invariant mass spectra and enriched the spectrum of charmed mesons. They have reported the resonances D 0 1 (2420) in the D * + π − final state and the   sistent for the quantitative study of quarkonium but not appropriate for mesonic systems consisting a heavy-light flavour quarks. For the hadrons consisting a heavy quark, QC D exhibits additional symmetries within the limit that heavy quark mass m Q becomes infinite compared with the typical QC D scale [4]. This symmetry can be realised systematically within relativistic Dirac equation employing the equally mixed scalar plus vector potential. This framework allows to make satisfactory predictions on excited charmed mesons. Hence, for the present study, we consider relativistic Dirac equation with equally mixed four-vector plus scalar powerlaw confinement potential for the single particle bound state energy of the quark and anti-quark. Using these single particle energies along with their j j coupling, the massess of the charmed mesons are computed. We fix the potential parameters for the known ground state and predict the masses of the radial and orbital excited states. Several models are available in literature dealing with decay properties of heavy-light mesons apart from their mass spectra [5][6][7][8][9][10]. The well known 3 P 0 model or quark pair creation model is very effective in strong decay study of mesons [11][12][13]. This model was initially proposed by L. Micu [14] and later developed by Y. Le and collaborators [15][16][17]. The relativized quark model spectroscopy incorpo-rating 3 P 0 model for strong decays [5] supports the D 0 J (2550) and D 0 * J (2600) as the radial excitation states. In their study they also identified D 0 * 1 (2760) and D 0 * J (2760) as the 1 3 D 1 and 1 3 D 3 states. But quoted that identification of D 0 J (2750) state requires further measurements. They concluded that the favourable assignment for states D 0 * J (3000) and D 0 J (3000) as D * 4 (1 3 F 4 ) and D(3 1 S 0 ), respectively. The 3 P 0 is a simplified version of the complicated theory. So, predictability of this model is not accurate. Also,the oscillator parameter affects the shape of wave-functions significantly in 3 P 0 model. More discussion regarding the analysis of uncertainties in 3 P 0 predictions can be available in Refs. [18][19][20]. Another approach to predict the strong decays of heavy-light mesons is based on heavy quark effective theory where P. Colangelo et al. has proposed the framework to classify the hadrons in doublets involving the channels with a light final pseudoscalar mesons [21][22][23][24]. The work exploring the heavy quark effective theory for charmed mesons is covered in Refs. [25,26] by Wang. Very, recently Colangelo et al. [27] incorporated other channels with a light final vector mesons into their previous study. We believe that in the case of relativistic Dirac model mass predictions as an input, Colangelo et al.'s approach is more consistent as compared to 3 P 0 model. In spite of the fact that Heavy quark effective theory contains many unknown phe-nomenological parameters, still HQET combined with chiral perturbation theory is a promising approach in the predictions of the strong decays of heavy-light hadrons [28]. This heavy quark symmetry impose constraints on the range of coupling constant between 0 and 1 for ground state charm mesons [29]. Thus, the present analysis for strong decays is carried out employing the Heavy quark effective theory at the leading order approximation where the mass and the spin degeneracy of the heavy hadrons are treated as approximate internal symmetry of the Lagrangian. The present study is organised as follows. In Sect. 2 we present our framework of relativistic Dirac model for heavy-light system and highlight the strong decays of heavy-light mesons within the framework proposed by Colangelo. Section 2 devoted to the results and analysis of the mass spectra and the strong decay widths. Finally, summary of the present study is included in Sect. 2.

The realtivistic Dirac model for heavy-light mesonic systems
The heavy-light mesons are composed of one heavy quark Q and other light quarkq, so the properties of heavy-light mesons are significantly influenced by relativistic effects. such systems can be systematically emphasised within relativistic Dirac framework. The form of the average flavour independent central potential for quarks and anti-quarks inside a meson is given as [30][31][32] Here, r represents radial distance from the meson centre of mass. Both λ (strength of the potential) and V 0 (depth of the potential) are the phenomenological parameters which are fixed for the ground state. We assume that Dirac equation is obeyed which can be derived from the Lagrangian density [30] The independent quark wave function ψ( r ) satisfies the Dirac equation [30,32] where E is the individual quark binding energy. The solution to the independent-quark wave function (normalised) can be written as [31,32,34,35] ψ nl j (r ) = ψ where and N nl j is the overall normalization constant [31,32,34,36]. The normalized spin angular part is defined by [34] Y l jm (r ) = m l ,m s l, m l , The two spinors χ 1 2 m s are the eigenfunctions of the spin operators and explicitly written as [34] The reduced radial parts g(r ) and f (r ) can be found to satisfy the following equations [6,31,32,35] and with the definition of quantum number κ as [34] Taking the form of V (r ) as given in Eq. (1) and introducing a dimensionless variable ρ = r r 0 with an arbitrary scale factor [31,32,36,37] Eqs. (9) and (10) reduces to the Schrödinger form [31,32,[36][37][38] and and ε represents dimensionless energy eigenvalue given by [31,32,[36][37][38] = Here, E D and m q represent the Dirac confinement energy of the quark having its rest mass m q . The numerical values of extracted Dirac confinement energy (E D ) are presented in Table 15 of Appendix A. For proper choices of κ Eqs. (13) and (14) can be be solved numerically.
The solutions g(ρ) and f (ρ) are normalized to get The wavefunctions for heavy-light charmed meson now can be constructed as the symmetric combinations of the positive energy solution ψ Q and the negative energy solution ψq of Eq. (5) and (6) and the corresponding mass of the Qq system can be written as where E Q/q D are obtained from (15) which include the centrifugal repulsion of the centre of mass also.
The mass of the state M2S+1 L J is obtained by adding the contributions from spin-spin, spin-orbit and tensor interactions of the confined one gluon exchange potential (COGEP) between quark and anti-quark [39,40] to M Qq . Thus, we write The spin-spin coupling part is defined from COGEP as [39,40] where σ represents the j j coupling constant. The angular brackets appearing in the term j 1 j 2 J M|ĵ 1 .ĵ 2 | j 1 j 2 J M contains square of Clebsch-Gordan coefficients.
The tensor and spin-orbit parts of confined one-gluon exchange potential (COGEP) [39,40] have the form (20) where N D , Nq are the normalization constants of quarkantiquark wavefunctions, S Qq = 3(σ Q .r )(σq .r ) − σ Q .σq andr =r Q −rq is the unit vector in the direction of r and where α s is the running strong coupling constant expressed as n f is number of flavours and for charmed mesons we take it as 3. The term λ Q .λq = − 4 3 is the color factor and independent of the flavour contents of the quarks [41].
The confined gluon propagators are given by [39,40,42] as and with α 1 = 1.035, α 2 = 0.3977 GeV, c 0 = 0.3418 GeV, γ = 0.8639 and c 1 = 0.4123 GeV GeV as in the previous study [32,42,43]. Other model parameters employed in the present study are listed in Table 4. Results obtained from the present study are tabulated in Tables 5, 6 and 7. The present results are compared with the available average experimental values reported by PDG [33].

Strong decays of heavy-light mesons
The underlying dynamics of the hadrons involving single heavy quark can be understood systematically by considering the heavy quark mass limit tends to infinity (m Q → ∞) formalised in a heavy quark effective theory (HQET). This H Q limit allows to classify heavy-light (Qq) mesons based upon the decoupling of the heavy quark spin s Q from the total angular momentum s l of the light degrees of freedom. Here, Q corresponds to c quark while q refers to light quarks u and d. The heavy quark spin s Q and total angular momentum of light d.o.f s l are separately conserved in strong interactions. This leads to classification of heavy mesons in doublets as per the different value of s l . Each doublet here has two states, spin partners having total spin J = s l ± 1 2 and parity P = (−1) l+1 where l is the orbital angular momentum of light d.o.f. and s j = l + s q , s q being the spin of light antiquark. and associated with doublets J P = (2 + , 3 + ) (P * 2 , P 3 ) and J P = (3 + , 4 + ) (P 3 , P * 4 ), respectively. This representation is for radial quantum number n = 1. These doublets can be characterized by effective fields (4 × 4 matrices) H a , S a , T a , X a , Y a , Z a and R a [44].
where H a (P, P * ) describes S-wave mesons; S a (P * 0 , P 1 ) and T a (P 1 , P * 2 ) associated with P-wave mesons; X a (P * 1 , P 2 ) Table 6 P-wave D meson (cq, q L HCb L HCb L HCb L HCb and Y a (P 2 , P * 3 ) represents D-wave mesons; Z a (P * 2 , P 3 ) and R a (P 3 , P * 4 ) related with F-wave mesons. Here, a = u, d or s is SU (3) light flavour index and ν corresponds to meson four-velocity. For radial quantum numbers n = 1, 2, 3, 4, . . . , heavy mesons with the same heavy flavour holds identical parity, time-reversal, charge conjugation properties and only differ in the mass. So, it is possible to combine them into effective fields H a , H a , H a . . . S a , S a , S a . . . where , , indicates n = 2, 3, 4 . . . states. Following this analogy, for example, for radial excitation n = 2 one can represent the sates by tilde (P,P * ).
The definition ξ = e iM fπ introduce octet of light pseudoscalar mesons in the theory through the matrix containing light pseudoscalar π , K and η fields. We take f π = 130 MeV for calculations. At the leading order approximation in the heavy quark mass and light quark momentum, the interaction Lagrangian terms L H , L S , L T , L X , L Y ,L Z and L R can be described as [44] where the definitions The vector and axial-vector currents are defined as The chiral symmetry breaking scale are the strong coupling constants. These constants are involved in the strong interactions of higher excited charmed mesons to the ground state +ve and −ve parity charmed mesons emitting light pseudoscalar mesons (π, η, K ) and can be fitted to experimental full widths.
The strong decay width to D ( * ) π , D ( * ) η and D ( * ) K using the chiral Lagrangians where A signifies scattering amplitude; M i and M f are the masses of initial and final charmed mesons,respectively; M p and p p denotes the mass and momentum of light pseudoscalar meson, respectively. Here, J is the total angular momentum of the initial heavy meson while σ shows the summation of all the polarization vectors of the total angular momentum j = 1, 2, 3 or 4. Below we collect the explicit expression of the two body strong decay widths of heavy light mesons for different channels [21,24]: 2. Decaying P wave meson: 3. Decaying P wave meson: 4. Decaying D wave meson: 5. Decaying D wave meson: 6. Decaying F wave meson: 7. Decaying F wave meson: The values of the coefficients C p involved are different for different pseudoscalar mesons: C π ± = C K ± = 1, C π 0 = 1 2 , C η = 2 3 [24]. The strong coupling constants and their notations depends upon the radial quantum number n. For transitions within n i = n f = 1 notations are g H,S,T,X,Y,Z ,R while for n i = 2 and n f = 1 they areg 2 H,S,T,X,Y,Z ,R . Higher order loop corrections are eliminated to bypass the introduction of new coupling constants. For the present study we adopt the approximation A μ = i ∂ μ M f π . If the momenta of the emitted light pseudoscalar mesons are not very small, the additional terms can be added to introduce unknown coupling constants. Moreover, the spin and the flavour violation correction having the order of O( 1 m Q ) to the heavy quark limit could be sizable, and in that case introduction of new coupling constants to the theory may not cancel out in the ratios of the decay widths. But we expect that their contribution would be much less than the leading order contributions. In decay rates the leading order unknown coupling constants can be theoretically predicted or can be evaluated from the available experimental data of decay widths. Successful predictions on coupling constants based on the QCD sum rules [21] and Lattice QCD [45]

Results and discussion
The mass spectra of the charmed mesons (cq); q = u or d are computed by employing the relativistic Dirac framework and listed in Tables 5, 6 [33]. However, if we incorporate the uncertainty of m c m c ≈ 1.5% as per m c = 1.27 ± 0.02 GeV [33] results into just 1% variation in the predicted D meson mass for 1 3 S 1 state and that for 1 1 S 0 state. Here, as up and down quark mass and their uncertainty are negligible (few MeV only) we don't consider it in the uncertainty estimations. The other parameters λ, V 0 and σ are the optimized potential parameters which are fixed to yield the experimental ground state masses. The variation of the parameter λ (potential strength) of 5% and 10% keeping other two fixed leads to the uncertainty of less than 1% and 2% in mass of 1S state. Also, The variation of the parameter σ (coupling strength) of 5% and10% keeping other two fixed gives uncertainty of less than 0.07% and 0.09% in mass of 1S state. The third parameter V 0 variation results into negligible variations in the mass of 1S state. So, the variations are not much compared to the optimized value of these parameters. However, the value of λ is slightly more sensitive than that of V 0 and σ . The detailed results on the sensitivity of charm quark mass and these model parameters are presented Table  16 and 17 of Appendix B.
The states D 2 (2460), D 1 (2420), D 0 (2400) and D 1 (2430) are also well defined 1P multiplets. Our predictions for these states are in excellent agreement with PDG listed values [33]. With the successful description of 1S and 1P multiplets, we are able to assign proper J P values of the newly observed states by L HCb Collaboration (1 3 P 2 , 1 3 P 1 ) to ground state charm meson by emitting a pseudoscalar meson are described by the coupling constant g T . The partial decay widths in terms of the g T are listed in Table 8 for D * 2 (2460). In Table  9 the computed ratio, is compared with those reported by various experimental groups. The average experimental value found to be 2.35 ± 0.6 which is in good agreement with the value of 2.26 predicted by the present study. This indicates the usefulness of the present formalism in the prediction of strong decay widths. From the present study, we found g T = 0.40 ± 0.003 which is in excellent agreement with 0.43 ± 0.05 reported by in [46] and 0.43 ± 0.01 reported by [26]. The details of computing the uncertainty in the estimation of coupling constants are shown in Appendix C.
The L HCb [2] and B AB AR [3] have observed states around the mass region 2550 MeV with unnatural parity and 2600 MeV with natural parity with different labels. But the mass and the widths of these states are so close making them to be considered having the same J P values. Also, their properties and mass range make them suitable candidates for radial excited (2S) states of the charm mesons. In the case of D(2550), the measured total width by L HCb [2] and B AB AR [3] contain very large error bar. The mass difference of L HCb [2] and B AB AR [3] is roughly 40 MeV (see Table 5). From our predictions the mass of 2 1 S 0 is 2543.07 MeV which is in good agreement with the PDG average value of 2564 ± 20 MeV. Similarly for the D(2600) our predicted mass is 2655 MeV. The mass and total width difference between LHCb [2] and BABAR [3] is around again 40 MeV and 47 MeV with large error bars, respectively. The latest effort of L HCb has reported D 1 (2680) with 2681.1 ± 5.6 MeV mass and 186.7 ± 8.5 MeV width. Hence, we identify the D(2580) and D * J (2650) as 2 1 S 0 and 2 3 S 1 ; respectively. The partial decay widths in terms of the effective coupling constants are listed in Table 8 for D J (2580) and D * J (2650) using our computed masses of 2 3 S 1 and 2 1 S 0 . The total decay width (2580) = 1101.39g 2 H MeV and (2650) = 2627.78g 2 H MeV and can be compared with the experimentally measured width to extract the effective coupling constant. Considering D J (2580), D * J (2650) as radial excited states the average effective coupling constantg H is deduced as 0.31 ± 0.017 which is close to 0.28 ± 0.01 reported by [26] but double the value of 0.14 ± 0.03 reported by Ref [46]. This, radial excited state can also decay to 1P states emitting the light pseudoscalar mesons (π, η, K ). Their incorporation requires additional introduction of the unknown coupling constants into the theory. So, to avoid complexity we do not consider them here.
B AB AR and L HCb also reported many states falling within the mass region between 2740 and 2800 MeV [2,3,47,48]. These observed states can be arranged into natural and unnatural parity states. The natural parity states are grouped with designation D * J (2760) (2763.5 ± 3.4 MeV) and unnatural parity states are labeled as D J (2740) (2737 ± 3.5 MeV). This mass range is predicted to be close to 1D and 2P multiplets [5]. The quantum numbers of these states as 2P multiplets are found to be inconsistent with L HCb [47,48] Table 8 for D * (2760) and D J (2740) for both possible assignments. The D + π − mode found to dominant in both 1 3 D 1 and 1 3 D 3 . The ratio, (D + π − ) (D * + π − for 1 3 D 1 is found to be 4.05 whereas for 1 3 D 3 it is 1.93. These suggests that D + π − mode is more dominant in 1 3 D 1 . The B AB AR Collaboration has observed D * (2760) signal in D + π − mode very close in masss to the D(2750) signal observed in D * + π − [3]. The state D * J (2760) was observed in both D * + π − and D + π − decay modes [2]. The LHCb has reported two D(2760) states: One assigned as 1 − [49] and other as 3 − [48] in different analysis. For the case of 1 − assignment, the reported mass and total width are 2781 ± 18 ± 11 ± 6 MeV and 177 ± 32 ± 20 ± 7 MeV respectively. The three quoted errors are statistical, experimental systematic and model uncertainties, respectively. While in the case of 3 − assignment (m = 2798 ± 7 ± 1 ± 7 MeV ; = 105 ± 18 ± 6 ± 23 (Isobar) and m = 2802 ± 11 ± 10 ± 3 Table 8 The strong decay widths (in MeV) for charmed resonance with possible spin-parity assignments. The ratio is calculated from (D * J →D * + π − ) . Fraction (in%) represents the percentage of the partial decay width with respect to total decay width    MeV ; = 154 ± 27 ± 13 ± 9 (K-matrix)), the error is less as compared to previous one [48]. Later, the state D * J (2760) is confirmed by L HCb [1] with J assignment as 3 in their latest attempt. So, we tentatively assign D J (2740) and D * J (2760) as 1 3 D 2 and 1 3 D 3 ; respectively. The average value of effective coupling constant g Y is obtained as 0.49 ± 0.039. This value is consistent with the 0.53 ± 0.13 [46] and 0.42 ± 0.02 [26]. At the same time, we found the ratio D(2760) 0 →D + π − D 0 (2740)→D * + π − equal to 1.40 which is larger than 0.42 ± 0.05 ± 0.11 reported by B AB AR [3]. Although the three states among the four 1D family are observed by experiments but one is still missing. It is very difficult to make precise measurements on the properties of the four overlapping states relying only upon their mass and widths. Also, D * π signal have large contribution from these overlapping states. So, further experimental efforts are required to resolve the above discrepancy and to understand these states more reliably.
The experimental value of the mass and the total decay width of D * J (3000) is 3008.1 ± 4.0 MeV and 110.5 ± 11.5 MeV [2]. The measured mass of the state D J (3000) is 2971.8 ± 8.7 MeV and width is 188.1 ± 44.8 MeV [2]. From the present analysis we may assign these states as the 3S, 2P or 1F multiplets. Considering both of them to be spin partners then the following five assignments of the J P for natural parity state D * J (3000) and unnatural parity state D J (3000) are possible.
Looking into the good agreement for 1S, 2S and 1P masses, we believe that for higher excited states, our predictions are reliable. The predicted masses from the present study for 3 3 S 1 (3143.28 MeV), 1 3 F 2 (3172.58 MeV), 1 3 F 4 (3252.55 MeV) and 3 1 S 0 (3044.31 MeV), 1 1 F 3 (3164.49 MeV), 1 3 F 3 (3208.21 MeV) are high as compared to 3008.1 ± 4 MeV and 2971.8 ± 8 MeV for D * J (3000) and D J (3000). However, the mass of 2 3 P 0 (2947.60 MeV) and 2 1 P 1 (2940.13 MeV) are comparable whereas mass of 2 3 P 2 (3000.90 MeV) and 2 3 P 1 (2972.50 MeV) are close to reported mass of D * J (3000) and D J (3000). The partial decay widths in terms of the effective coupling constants are listed in Tables 10 and 11 for D * (3000) and D J (3000) respectively where we have used our computed masses. Among the five possibilities later two can be completely ruled out as 1 3 F 2 , 1 1 F 3 resulting very small decay width and 1 3 F 4 , 1 3 F 3 have very large decay width which is far from the experimentally predicted widths of the D * J (3000) and D J (3000) [2]. The state D * J (3000) decays to D + π − whereas D J (3000) decays to D * + π − final state [2]. If we identify D * J (3000) as 3 3 S 1 then the dominant mode is D * + π − which is not favoured by experiment. The other modes such as D 0 π 0 , D 0 η are very small. However, the two remaining possibilities leads to the J P of D J (3000) as 1 + . If we consider the D * (3000) 2 3 P 0 then the D * π mode is completely forbidden and Dπ mode is dominant. This fact is consistent with the experimental observation. Also, the D + s K − and D 0 η modes are considerably large which is supported by experimental analysis. Similarly if we assign D J (3000) as 2 1 P 1 , the D + * s K − and D 0 * η modes are sufficiently large as compared to 2 3 P 1 . With all such considerations, we identify D * J (3000) as 2 3 P 1 and D J (3000) as 2 1 P 1 . The previous study Ref. [50] supports this assignments for D * J (3000) and D J (3000) states while the other studies [18] argues that D * J (3000) could be 1 3 F 4 or 1 3 F 2 and Table 10 The strong decay widths (in MeV) for resonance D * J (3000) with possible spin-parity assignments  Table 11 The strong decay widths (in MeV) for resonance D J (3000) with possible spin-parity assignments D J (3000) could be 1F 3 or 2P 1 . Ref. [7] suggests D * J (3000) is a 1 3 F 4 and D J (3000) is a 2P 1 state. All these possibilities also can not be ignored and more precise measurements are needed from the experimental side to clarify the puzzles in D(3000) states. Considering D * J (3000) and D J (3000) as 2P(0 + , 1 + ) 1 2 multiplets, the predicted value of the effective coupling constantg S found to be 0.10 ± 0.015 which would be useful for the future investigations on heavy-light systems. The numerical values of the decay widths and B Rs emitting the light vector meson channels such as Dρ, D s K * and Dω can also be incorporated to make precise predictions based on the heavy quark effective theory framework [27]. Recently, the L HCb group has observed a state labelled as D * 2 (3000) with mass 3214 ± 29 ± 33 MeV and total width 186 ± 38 ± 34 MeV [1]. The labelling of this newly observed state in B − → D + π − π − decays and previously reported state D * J (3000) are resembled to each other. However, the L HCb collaboration assigned J P = 2 + for D * 2 (3000) while D * J (3000) has natural parity. The mass energy difference between D * 2 (3000) and D * J (3000) is 206 MeV which suggests that they both are different states [51]. The mass range 3214 ± 29 ± 33 MeV likely to fall in the mass spectra of 2P, 3P and 1F multiplets. From the present study mass of 1 3 F 4 is close to 3214 ± 29 ± 33 MeV but the quantum numbers are not consistent with this assignment and for 3 3 P 2 the mass difference is more than 240 MeV. So, the two most probable assignments for D * 2 (3000) are 1 3 F 2 or 2 3 P 2 . From the present study the mass of 1 3 F 2 is 3172.58 MeV while the mass of 2 3 P 2 is 3000.90 MeV. Allowed decay channels for D * 2 (3000) are given in terms of coupling constant g Z andg T in Table 12. The two different possibilities bring out different decay widths. The ratio R ( * )0 p of various partial widths and branching fractions can be useful in the identification of the new signals with predicted states in experiments. Using the partial widths such ratios are defined as [27] The ratios for both the possibilities are presented in Table  13. It can be noted that the ratio R 0 K has less sensitivity to the 2P and 1F identification while R 0 π is highly sensitive. In both the cases D + π − mode is found to be dominant. So, to get more insights we also considered the ratio D * 2 (3000)→D * + π − D * 2 (3000)→D + π − for 1 3 F 2 and 2 3 P 2 states. This ratio is found to be 0.95 for 2 3 P 2 while it is 0.38 for 1 3 F 2 . This indicates that D * + π − mode is more dominant in 2 3 P 2 as compared to 1 3 F 2 . On the Table 12 The strong decay widths (in MeV) for resonance D * 2 (3000) with possible spin-parity assignments. The ratio is calculated from (D * J →D * + π − ) . Fraction (in%) represents the percentage of the partial decay width with respect to total decay width Decay mode 2 3 P 2 ( 3 other hand, experimentally D * + π − mode is suppressed in D * 2 (3000) state [1]. This fact is consistent with our predicted decay modes. So, we identify D * 2 (3000) as 1 3 F 2 . As shown in Table 9, D * 2 (3000) being 1 3 F 2 candidate, its total decay width as 12008.54 g 2 Z MeV results into, g Z = 0.12 ± 0.012. After identifying the D * 2 (3000) as 1 3 F 2 ( 5 2 2 + ) one can have further considerations on its spin partner 1 1 F 3 ( 5 2 3 + ). The information about them in literature is very less so it will be helpful for theorists and experimentalists for future investigations on these states. Our predicted mass of 1 1 F 3 is 3164.49 MeV. The Okubo-Zweig-Iizuka (OZI) allowed decay channels for 1 1 F 3 are listed in Table 14. The states 1 1 F 3 can also decays to p-wave mesons through emission of light pseudoscalar mesons but due to small phase space those modes are neglected in the present work. If we consider the value g Z =0.12 then the total decay width of 1 1 F 3 found to be 67.15 MeV. Thus, the total width of 1 1 F 3 is narrower than its spin partner 1 3 F 2 ( = 110.5 ± 11.5 MeV). The branching fraction 44.18% for the decay mode D * + π − suggests as an appropriate channel for the experimental search of 1 1 F 3 . Finally we conclude that the predicted masses of several states of charmed mesons are in good agreement with the experimental values. And the mass spectroscopy is found to be consistent in predictions of strong decays of various Dmesonic states.

Summary
In this article, we obtained the mass spectra of the heavylight charmed mesons within the framework of the relativistic formalism. Our predicted masses are in excellent agreement for the well established low-lying states of charmed mesons. Further, incorporating the heavy quark effective theory at the leading order approximation, we obtained the strong decays of the experimentally observed states and identified D J (2560) as 2 1 S 0 , D * J (2680) as 2 3 S 1 , D J (2740) as 1 3 D 2 , D * J (2760) as 1 3 D 3 , D * J (3000) as 2 3 P 0 , D J (3000) as 2 1 P 1 and D * 2 (3000) as 1 3 F 2 open charm excited states. The effective coupling constants g T ,g H , g Y ,g S and g Z are also computed. We found D * + π − channel is the most suitable mode for experimental search for 1 1 F 3 . The computed decay widths and masses for higher excited states would be useful for the identification of future experimental observation of higher open charm states at B AB AR, L HCb, B E SI I I .

Appendix A: Numerical single particle Dirac confinement energy of the quarks and antiquarks
The Dirac confinement energy (E D ) is extracted from Eq. (15) for quark(Q) and anti-quark(q) in accordance with Eqs. (5), (6) and (13), (14). These values for Q andq are tabulated Table 15.  By considering the charm quark mass, m c = 1.27 + 0.02 MeV as the upper bound and m c = 1.27 − 0.02 MeV as the lower bound, the ground state masses of D− mesons for (1 3 S 1 , 1 1 S 0 ) states are computed. The results are given in Table 16. As one can see from Table 16  The other model parameters λ, V 0 and σ are the optimized potential parameters which are fixed for known ground state. Below we present the effect on 1S mass due to an assumed 5% and 10% variations in the parameters keeping other two fixed. The results are shown in Table 17. One can see from Table 17 that the variations are not much compared to the optimized value of these parameters used in the present study. However, the value of λ is slightly more sensitive than that of V 0 and σ .

Appendix C:Illustration of uncertainty estimation in the coupling constants
In general the total decay width from the present study (say X × g 2 i ; i = T,H , Y,S, Z ) can be compared with experimentally measured width (say E x p ) to extract respective coupling constants.
The uncertainty in the estimations of the strong coupling constants as appear in Table 8 Column 4 due to uncertainty in the experimental width as shown in Table 8 Column 7 are computed as below.
Taking the log and the derivative on both the sides we get, Here, the δ E x P is the uncertainty in E x P and δ X that in X . In the present case δ X = 0 results into As a sample case, consider the total strong decay width (2460) = 286.39g 2 T MeV from Table 8 and is equated to the experimentally measured width of (2460) = 47.0 ±0.8 MeV, we get  The same process is followed for g z . For computing the gH coupling constant we have two relations as shown in Table 8 . It follows that, the total decay width Thus we obtain gH = 0.31 ± 0.017 The same method is followed for g Y and gS. The estimated uncertainty in g i for i = T,H , Y,S, Z are given in Table 18.