D meson spectroscopy and their decay properties using Martin potential in a relativistic Dirac formalism

For the present study, we have used the Martin-like potential for the quark confinement. Our predicted states in the S-wave, 23S1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2\ ^3S_1$$\end{document} (2605.86 MeV) and 21S0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2\ ^1S_0$$\end{document} (2521.72 MeV), are in very good agreement with experimental results of 2608±2.4±2.5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2608\pm 2.4 \pm 2.5$$\end{document} MeV and 2539.4±4.5±6.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2539.4\pm 4.5\pm 6.8$$\end{document} MeV, respectively, reported by the BABAR Collaboration. The calculated P-wave D meson states, 13P2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1^3P_2$$\end{document} (2462.50 MeV), 13P1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1^3P_1$$\end{document} (2407.56 MeV), 13P0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1^3P_0$$\end{document} (2373.82 MeV) and 11P1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1^1P_1$$\end{document} (2423.28 MeV), are in close agreement with experimental average (Particle Data Group) values of 2462.6±0.7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2462.6 \pm 0.7 $$\end{document} MeV, 2427±26±25\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2427 \pm 26 \pm 25$$\end{document} MeV, 2318±29\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2318 \pm 29 $$\end{document} MeV and 2421.3±0.6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2421.3 \pm 0.6 $$\end{document} MeV, respectively. The pseudoscalar decay constant (fP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_P$$\end{document}= 202.57 MeV) of the D meson is in very good agreement with the experiment as well as with the lattice predictions. The Cabibbo favoured nonleptonic decay branching ratios, BR(D0→K-π+)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(D^0\rightarrow K^- \pi ^+)$$\end{document} of 4.071%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.071\,\%$$\end{document} and BR (D0→K+π-)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(D^0\rightarrow K^+ \pi ^-)$$\end{document} of 1.135×10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.135 \times 10^{-4} $$\end{document}, are also in very good agreement with the respective experimental values of 3.91±0.08%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ 3.91 \pm 0.08\,\%$$\end{document} and (1.48±0.07)×10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(1.48\pm 0.07) \times 10^{-4}$$\end{document} reported by CLEO Collaboration. The mixing parameters of the D0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^0 $$\end{document}–D¯0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \bar{D}^0$$\end{document} oscillation, xq\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_q$$\end{document} (5.14 ×10-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times 10^{-3}$$\end{document}), yq\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_q$$\end{document} (6.02 ×10-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times 10^{-3}$$\end{document}) and RM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_M$$\end{document} (3.13 ×10-5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times 10^{-5}$$\end{document}), are in very good agreement with BaBar and Belle Collaboration results.


Introduction
Very recently, experiments at LHCb [1] have reported a large number of D J resonances in the mass range of 2.0 GeV/c 2 to 4.0 GeV/c 2 , of which many belong to natural excited states of the D meson, while quite a number of them belong to unnatural states [1]. It is important and necessary to exhaust the possible conventional descriptions of qQ excitations [2,3] before resorting to more exotic interpretations [4][5][6]. a e-mail: mnshah09@gmail.com b e-mail: azadpatel2003@gmail.com c e-mail: p.c.vinodkumar@gmail.com Further theoretical efforts are still required in order to explain satisfactorily the recent experimental data concerning these open-charm states.
Apart from the challenges posed by the exotics, there are also many states which are admixtures of their nearby natural states. For example, the discoveries of new resonances of D states such as D(2550) [7], D(2610) [7], D(2640) [8], D(2760) [7] etc. have further generated considerable interest towards the spectroscopy of these open-charm mesons. The study of the D meson is of special interest as it is a hadron with two open flavours (c,ū ord) which restricts its decay via strong interactions. The ground state (D, D * ) mesons provide us with a clean laboratory to study weak decay and are useful to study the electromagnetic transitions. The masses of low-lying 1S and 1P J states of the D mesons are recorded both experimentally [2,3] and theoretically [9][10][11][12][13][14]. Though lattice QCD and QCD sum rules are quite successful, their predictions for the excited states of the open flavour mesons in the heavy sector are very few. However, recent experimental data on excited D-states are partially inconclusive and require a more detailed analysis involving their decay properties. The understanding of the weak transition form factors of heavy mesons is important for a proper extraction of the quark mixing parameters, for the analysis of nonleptonic decays and CP-violating effects. The QCD sum rule (QSR) [15][16][17][18][19] is a non-perturbative approach to evaluate the hadron properties by using the correlator of the quark currents over the physical vacuum and it is implemented with the operator product expansion (OPE). Lattice QCD (LQCD) [20][21][22] is also a non-perturbative approach to use a discrete set of spacetime points (lattice) to reduce the analytically intractable path integrals of the continuum theory to a very difficult numerical computation. QCD sum rules are suitable for describing the low q 2 region of the form factors; lattice QCD gives good predictions for high q 2 . As a result these methods do not provide a full picture of the form factors and, more significant, for the relations between the various decay channels. Potential models provide such relations and give the form factors in the full q 2 -range. However, the potential models are not derived from the basic principles of QCD. In particular, the Martin potential is not really inspired from QCD unlike the case of the Cornell potential. Yet it gives spectacular phenomenological predictions of the hadron spectra. In the case of heavy-light flavour systems (qQ/q Q), its connection with heavy quark effective theory (HQET) remains unclear. Thus it is important to check its consistency for the successful predictions of the properties of other well recorded open flavour hadronic systems. With this perspective we employ the Martin-like potential for the study of open-charm mesonic properties.
Thus any attempts towards the understanding of these newly observed states become very important for our understanding of the light quark/antiquark dynamics within qQ/Qq bound states. So, a successful theoretical model aims to provide important information as regards the quarkantiquark interactions and the behaviour of QCD within the doubly open flavour hadronic system. Though there exist many theoretical models [9][10][11] in the study of the hadron properties based on its quark structure, the predictions for low-lying states are off by 60-90 MeV with respect to the respective experimental values. Moreover, the issues related to the hyperfine and fine structure splitting of the mesonic states, their intricate dependence with the constituent quark masses and the running strong coupling constant are still unresolved. Though the validity of nonrelativistic models is very well established and had significant success in the description of heavy quarkonia, disparities exist in the description of meson containing light flavour quarks or antiquarks.
For any successful attempt to understand these states we should not only be able to satisfactorily predict the mass spectra but also be able to predict their decay properties. For better predictions of the decay widths, many models have incorporated additional contributions such as radiative and higher order QCD corrections [13,[23][24][25][26][27]. Thus, in this paper we make an attempt to study properties like the mass spectrum, decay constants and other decay properties of the D meson based on a relativistic Dirac formalism. We investigate the heavy-light mass spectra of the D meson in this framework with a Martin-like confinement potential as in the case of D s mesons studied recently [28].
Along with the mass spectra, the pseudoscalar decay constants of the heavy-light mesons have also been estimated in the context of many QCD-motivated approximations. The predictions of such methods spread over a wide range of values [29,30]. It is important thus to have reliable estimate of the decay constant as it is an important parameter in many weak processes such as quark mixing, CP violation etc. The leptonic decay of charged meson is another important annihilation channel through the exchange of virtual W boson.
Though this annihilation process is rare, we find clear experimental signatures due to the presence of a highly energetic leptons in the final state. The leptonic decays of mesons entails an appropriate representation of the initial state of the decaying vector mesons in terms of the constituent quark and antiquark with their respective momenta and spin. The bound constituent quark and antiquark inside the meson are in definite energy states having no definite momenta. However, one can find the momentum distribution amplitude for the constituent quark and antiquark inside the meson just before their annihilation to a lepton pair. Thus, it is appropriate to compute the leptonic branching ratio and compare our result with the experimental values as well as with the predictions based on other models.

Theoretical framework
From the first principles of QCD, the non-perturbative multigluon mechanism is unfeasible to estimate theoretically, but this mechanism is produced in the framework of the quark confining interaction of the meson. On the other hand there exists ample experimental support for the quark structure of the hadrons. This is the origin of phenomenological models which are proposed to understand the properties of hadrons and quark dynamics at the hadronic scale. To a first approximation, the confining part of the interaction is believed to provide the zeroth-order quark dynamics inside the meson through the quark Lagrangian density, In the present study, we assume that the constituent quarkantiquark inside a meson is independently confined by an average potential of the form [28,31] where λ is the potential strength. Here we use the index ν = 0.1, the form of the potential is Martin-like. In the stationary case, the spatial part of the quark wave functions ψ(r) satisfies the Dirac equation given by The two component (positive and negative energies in the zeroth order) form represents the solution of the Dirac equation as where (6) and N nl j is the overall normalisation constant. The normalised spin angular part is expressed as Here the spinors χ 1 2 m s are eigenfunctions of the spin operators, The reduced radial part g(r ) of the upper component and f (r ) of the lower component of Dirac spinor ψ nl j (r ) are the solutions of the equations given by and It can be transformed into a convenient dimensionless form given by [32] and where ρ = (r/r 0 ) is a dimensionless variable with the arbitrary scale factor chosen conveniently as and is a corresponding dimensionless energy eigenvalue defined as Here, it is suitable to define a quantum number κ by Equations (11) and (12) now can be solved numerically [33] for each choice of κ.
The solutions g(ρ) and f (ρ) are normalised to get Using Eqs. (5) and (6), the wavefunction for a D(cq) meson now can be constructed and the corresponding mass of the quark-antiquark system can be written as where E D Q/q = E Q/q are obtained using Eqs. (14), (15), and (16). The expression for E Q/q thus contains the centrifugal repulsion term which includes the centre-of-mass correction also. For the spin triplet (vector) and spin singlet (pseudoscalar) state, the choices of ( j 1 , j 2 ) are l 1 + 1 2 , l 2 + 1 2 and l 1,2 + 1 2 , l 2,1 − 1 2 , respectively. The previous work of the independent quark model within the Dirac formalism by [28,31] has been extended here by incorporating the spin-orbit and tensor interactions of the confined one-gluon exchange potential (COGEP) [34,35], in addition to the j-j coupling of the quarkantiquark. Finally, the mass of the specific 2S+1 L J states of the Qq system is expressed as The spin-spin part is defined here as where σ is the j-j coupling constant. The expectation value of j 1 j 2 J M|ĵ 1 .ĵ 2 | j 1 j 2 J M contains the ( j 1 . j 2 ) coupling and the square of the Clebsch-Gordan coefficients. The tensor and spin-orbit parts of the confined one-gluon exchange potential (COGEP) [34,35] are given as where 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 strong coupling constant; it is computed as with n f = 3 and QCD = 0.250 GeV obtained by fixing the strong running coupling strength α s (M 2 z ) at the Z-boson mass (M Z = 91 GeV) scale to be equal to 0.11838. In Eq. (22) the spin-orbit term has been split into symmetric (σ Q + σ q ) and antisymmetric (σ Q − σ q ) terms. It is to be noted that in the relativistic Dirac formalism, the relativistic mass of the quark/antiquark contains contribution from its rest mass (m Q/q ) as well as from the kinetic and interaction part of the Hamiltonian. Hence (E Q/q + m Q/q ) appears here for its inertial mass. Though m represents the current quark mass in the MS scheme, (E Q/q + m Q/q ) corresponds to the constituent quark mass in its nonrelativistic reduction. For the present calculations, we have taken the quark mass parameters as quoted in the particle data group (PDG) (2014).
We have adopted the same parametric form of the confined gluon propagators as are given by [34,35] and with α 1 = 0.011, α 2 = 0.036, c 0 = 0.1017 GeV, c 1 = 0.1522 GeV, γ = 0.0109 as in our earlier study [28]. Other optimised model parameters employed in the present study are listed in Table 1. The current charm quark mass of 1.27 GeV is taken from the particle data group (PDG) [2,3]. In the case of l = 0 orbitally excited states, we find small variations in the choice of λ for the l = 0 states due to the centrifugal repulsion from the centre of mass of the bound system, which is proportional to (n + l). This centrifugal repulsion thus incorporates the centre-of-mass correction.
The computed S-wave masses and other P-wave and Dwave masses of D meson states are listed in Tables 2 and  3, respectively. We have also calculated the mixed state of 3 P 1 and 1 P 1 and compared with the experimental results and other theoretical results in Table 4. A statistical analysis of the sensitivity of the model parameters (i.e. potential strength (λ) and j-j coupling strength σ in the present case) shows about 0.76 % variations in the binding energy with 5 % changes in the parameters λ and σ . Figure 1 shows the energy level diagram of the D meson spectra along with available experimental results.

Magnetic (M1) transitions of open-charm meson
The decay widths of energetically allowed radiative transitions ( A → B + γ ) of vector and pseudoscalar states of D meson are computed from the spectroscopic studies. The magnetic transition correspond to spin flip and hence the vector meson decay to pseudoscalar V → Pγ represents a typical M1 transition. Such transitions are experimentally important for the identification of newly observed states. Assuming that such transitions are single vertex processes governed mainly by photon emission from independently confined Semi-relativistic model [36], quasi-potential approach [37], relativistic quark-antiquark potential (Coulomb plus power) model [14], nonrelativistic constituent quark model [38], relativistic quark model [1], lattice QCD [LQCD] [20], QCD sum rule [QSR] [15,18]  quark and antiquark inside the meson, the S-matrix elements in the rest frame of the initial meson is written in the form The common choice of the photon field A μ (x) is made here in Coulomb-gauge with (k, λ) as the polarisation vector of the emitted photon having an energy momentum (k 0 = |k|, k) in the rest frame of A. The quark field operators find possible expansions in terms of the complete set of positive and negative energy solutions given by Eqs. (5) and (6) in the form where the subscript q stands for the quark flavour and ζ represents the set of Dirac quantum numbers. Here b qζ and b † qζ are the quark annihilation and the antiquark creation operators corresponding to the eigenmodes ζ . After some standard calculations (the details of the calculations can be found in Refs. [40][41][42]), the S-matrix elements can be expressed as Here Jq mm (k, λ) = e q d 3 r exp(−ik · r) One can reduce the above equations to the simple forms of and where K = k × (k, λ). Equation (28) further can be simplified to where μ q (k) is expressed as where j 1 (kr) is the spherical Bessel function and the energy of the outgoing photon in the case of a vector meson undergoing a radiative transition to its pseudoscalar state, for instance, D * → Dγ is given by The relevant transition magnetic moment is expressed as Now, the magnetic (M1) transition width of D * → Dγ can be obtained: The computed transition widths of the low-lying S-wave states are tabulated in Table 5 and are compared with other model predictions.

Decay constant of the D meson
In the study of leptonic or nonleptonic weak decay processes, the decay constant of a meson is an important parameter. The decay constant ( f p ) of the pseudoscalar state is obtained by parameterizing the matrix elements of weak current between the corresponding meson and the vacuum as [43] 0|qγ μ γ 5 c|P μ = i f p P μ .
It is possible to express the quark-antiquark eigenmodes in the ground state of the meson in terms of the corresponding momentum distribution amplitudes. Accordingly, the eigenmodes, ψ (+) A in the state of definite momentum p and spin projection s p , can be expressed as where U q ( p, s p ) is for the usual free Dirac spinors.
In the relativistic quark model, the decay constant can be expressed through the meson wave function G q ( p) in the momentum space [41,44] Here M p is mass of the pseudoscalar meson and I p and J p are defined as respectively. Here, and E p i = k i 2 + m 2 q i . The computed decay constants of the D meson from 1S to 4S states are tabulated in Table 6. The present result for 1S state is compared with experimental as well as other model predictions. There are no model predictions available for a comparison of the decay constants of the 2S to 4S states.

Leptonic decay of the D meson
The leptonic decays of open flavour mesons belong to rare decay [46,47], they have clear experimental signatures due to the presence of a highly energetic lepton in the final state. Such decays are very clean due to the absence of hadrons in  the final state [51]. Charged mesons produced from a quark and antiquark can decay to a charged lepton pair when these objects annihilate via a virtual W ± boson as given in Fig. 2. The leptonic width of the D meson is computed using the relation given by [2,3] (46) in complete analogy to π + → l + ν. These transitions are helicity suppressed; i.e., the amplitude is proportional to m l , the mass of the lepton l.
where τ is the experimental lifetime of the respective D meson state. The computed leptonic widths are tabulated in Table 7 along with other model predictions as well as with

Nonleptonic decays of the D meson
The study of flavour changing decays of heavy flavour quarks is useful for determining the parameters of the Standard Model and for testing phenomenological models which include strong effects. The interpretation of the nonleptonic decays of the c-meson within a hadronic state is complicated by the effects of the strong interaction and by its interplay with the weak interaction. The nonleptonic decays of heavy mesons can be understood in this model and we assume that Cabibbo favoured nonleptonic decays proceed via the basic process (c → q + u +d; q ∈ s, d), and the decay widths are given by [43] for q = s and Here, C f is the colour factor and (|V cs |, |V cd |, |V us |) are the CKM matrices. f π is the decay constant of π meson and its value is taken as 0.130 GeV. Here, f + (q 2 ) is the form factor and the factor λ(M 2 D , M 2 K + , M 2 π ) can be computed as The renormalised colour factor without the interference effect due to QCD is given by (C 2 A + C 2 B ). The coefficients C A and C B are further expressed as [43] where and where M W is the mass of W meson. Consequently, the form factors f ± (q 2 ) correspond to the D final state are related to the Isgur-Wise function as [43] The Isgur-Wise function, ξ(ω) can be evaluated according to the relation given by [65] where E q is the binding energy of decaying meson and ω is given by In a good approximation the form factor f − (q 2 ) does not contribute to the decay rate, so we have neglected it here. The heavy flavour symmetry provides a model-independent normalisation of the weak form factors f ± (q 2 ) either at q = 0 or q = q max , and we have applied q = q max in Eqs. (48) and (49) for nonleptonic decay. From the computed exclusive semileptonic and hadronic decay widths, the branching ratios are obtained: here the lifetime (τ ) of D (τ D + = 1.040 ps −1 and τ D 0 = 0.410 ps −1 ) is taken as the world average value reported by Particle Data Group [2,3]. The decay widths and their branching ratios are listed in Table 8 along with the known experimental and other theoretical predictions for comparison.  [76] 7 Mixing parameters of the D 0 -D 0 oscillation A different D 0 decay channel [66][67][68][69][70] has been reported by three experimental groups as evidence of the D 0 -D 0 oscillation. We discuss here the mass oscillation of the neutral open-charm meson and the integrated oscillation rate using our spectroscopic parameters deduced from the present study. In the standard model, the transitions D 0 -D 0 andD 0 -D 0 occur through the weak interaction. The neutral D meson mixes with the antiparticle leading to oscillations between the mass eigenstates [2,3]. In the following, we adopt the notation introduced in [2,3], and assume CPT conservation in our calculations. If CP symmetry is violated, the oscillation rates for meson produced as D 0 andD 0 can differ, further enriching the phenomenology. The study of CP violation in D 0 oscillation may lead to an improved understanding of possible dynamics beyond the standard model [71][72][73].
The time evolution of the neutral D−meson doublet is described by a Schrödinger equation with an effective 2 × 2 Hamiltonian given by [43,74] i d dt where the M and matrices are Hermitian, and they are defined as CPT invariance imposes The off-diagonal elements of these matrices describe the dispersive and absorptive parts of the D 0 -D 0 mixing [77]. The two eigenstates D 1 and D 2 of the effective Hamiltonian matrix (M − i 2 ) are given by is small, and a power expansion of |q/ p| 2 yields Therefore, considering both Eqs. (71) and (72), the CPviolating parameter given by is expected to be very small: ∼ O(10 −3 ) for the D 0 −D 0 system. In the approximation of negligible C P violation in mixing, the ratio / m is equal to the small quantity | 12 /M 12 | of Eq. (72); it is hence independent of the CKM matrix elements, i.e., the same for the D 0 -D 0 system. Theoretically, the hadron lifetime (τ D 0 ) is related to 11 (τ D 0 = 1/ 11 ), while the observables m and are related to M 12 and 12 as [2,3] and = 2| 12 |.
The gluonic correction can be found from a different model, like the Wilson coefficient and the evolution of the Wilson coefficient from the new physics scale [73]. We have used the values of the gluonic correction (η D 0 = 0.86; η D 0 = 0.21) from [81,82]. The bag parameter B D 0 = 1.34 is taken from the lattice result of [83], while the pseudoscalar mass (M D 0 ) and the pseudoscalar decay constant ( f D ) of the D mesons are the values obtained from our present study using a relativistic independent quark model using a Martinlike potential. The values of m s (0.1 GeV), M W (80.403 GeV) and the CKM matrix elements V cs (1.006) and V us (0.2252) are taken from the Particle Data Group [2,3]. The resulting mass oscillation parameter m are tabulated in Table 9 with the latest experimental results.
The integrated oscillation rate (χ q ) is the probability to observe aD meson in a jet initiated by ac quark. As the mass difference m D is a measure of the frequency of the change from a D 0 into aD 0 or vise versa. This change is reflected in either the time-dependent oscillations or in the time-integrated rates corresponding to the di-lepton events having the same sign. The time evolution of the neutral states from the pure |D 0 phys or |D 0 phys state at t = 0 is given by which means that the flavour states remain unchanged (g + ) or oscillate into each other (g − ) with time-dependent probabilities proportional to Starting at t = 0 with an initially pure D 0 , the probability for finding a D 0 (D 0 ) at time t = 0 is given by |g + (t)| 2 (|g − (t)| 2 ). Taking |q/ p| = 1, one gets Conversely, from an initially pureD 0 at t = 0, the probability for finding aD 0 (D 0 ) at time t = 0 is also given by |g + (t)| 2 (|g − (t)| 2 ). The oscillation of D 0 orD 0 as shown by Eq. (81) gives m directly. Integrating |g ± (t)| 2 from t = 0 to t = ∞, we get  [86] 0 . 0 4 +0.7 −0.6 ×10 −3 [87] 0 . 0 2 ± 0.47 ± 0.14 ×10 −3 where = D = ( 1 + 2 )/2. The ratio where reflects the change of a pure D 0 into aD 0 , or vice versa. The time-integrated mixing rate relative to the timeintegrated right-sign decay rate for semileptonic decays [2,3] is In the standard model, CP violation in charm mixing is small and |q/ p| ≈ 1.
For the present estimation of these mixing parameters, x q , y q and χ q , we employ our predicated m values and the experimental average lifetime of PDG [2,3] of the D-meson.

Results and discussion
We have studied here the mass spectra and decay properties of the D meson in the framework of the relativistic independent quark model. Our computed D meson spectral states are in good agreement with the reported PDG values of the known states. The predicted masses of the S-wave D meson state 2 3 S 1 (2605.86 MeV) and 2 1 S 0 (2521.72 MeV) are in very good agreement with the respective experimental results of 2608 ± 2.4 ± 2.5 MeV [39] and 2539.4 ± 4.5 ± 6.8 MeV [39] by the BABAR Collaboration. The expected results of other S-wave excited states of the D meson are also in good agreement with other reported values [14,[36][37][38]. The predicted P-wave D meson states, 1 3 P 2 (2462.50 MeV), 1 3 P 1 (2407.56 MeV), 1 3 P 0 (2373.82 MeV) and 1 1 P 1 (2423.28 MeV), are in good agreement with experimental [2,3] results of 2462.6 ± 0.7 MeV, 2427 ± 26 ± 25 MeV, 2318 ± 29 MeV and 2421.3 ± 0.6 MeV, respectively. We have also compared lattice QCD and QCD sum rule results with our predicted results where the numerical values in Table 3 for lattice QCD results are extracted from the energy level diagram available in [20]. With reference to the available experimental masses of D-mesonic states, we observe that the LQCD predictions [20] are off by a standard deviation of ±58.52 and those found by the QCD sum rule [15,18] predictions are off by ±59.22, while the predicted calculations show a standard deviation of ±21.88.
In the limit of one heavy quark, the 1 + resonance is expected to show a very simple mixing pattern [88,89]. The angular momentum j q = s q + L of a light quark is a good quantum number, which is a conserved quantum number, while the angular momentum s Q of a heavy quark is not a good quantum number, so the physically observed two 1 + states are mixed states of the 3 P 1 and 1 P1. Accordingly, the mixed states can be expressed as [88,89], Further, we can write the masses of these states in terms of the predicted masses of 3 P 1 and 1 P 1 states as which falls within the error bars of the experimentally observed state D(2427 ± 26 ± 25) and which is very close to D 1 (2420±0.6). The calculated mixing states of 3 P 1 and 1 P 1 are listed in Table 4. We look forward to see more precise measurements of the mass of the D 1 (2430) state as it supports the right contributions from the tensor and spin-orbit interactions.
In the relativistic Dirac formalism, the spin degeneracy is primarily broken; therefore to compare the spin average mass, we employ the relation of The spin average or the centre of weight masses M CW are calculated from the known values of the different meson states and are compared with other model predictions [14,37] in Table 10. The table also contains the different spin dependent contributions for the observed state. Mass splittings of   [45] and other model predictions [14,37] in Table 11. The precise experimental measurements of the masses of the D meson states provided a real test for the choice of the hyperfine and the fine structure interactions adopted in the study of the D meson spectroscopy. A recent study of the D meson mass splitting in lattice QCD [LQCD] [45] using 2 ± 1 flavour configurations generated with the Clover-Wilson fermion action by the PACS-CS Collaboration [45] has been listed for comparison. The present results as seen in Table 8 are in very good agreement with the respective experimental values over the lattice results [45]. In this table, the present results, on average, are in agreement with the available experimental value within 12 % variations, while the lattice QCD predictions [45] show 30 % variations.
The magnetic transitions (M1) can probe the internal charge structure of hadrons, and therefore they will likely play an important role in determining the hadronic structures of the D meson. The present M1 transitions widths of the D meson states as listed in Table 5 are in accordance with the model prediction of [49] while the upper bound provided by PDG [2,3] is very wide. We do not find any theoretical predictions for M1 transition width of excited states for comparison. Thus we only look forward to see future experimental support to our predictions.
The calculated pseudoscalar decay constant ( f P ) of the D meson is listed in Table 6 along with other model predictions as well as experimental results. The value of f D (1S) = 202.57 MeV obtained in our present study is in very good agreement with other theoretical predictions for 1S state. The predicted f D for higher S-wave states are found to increase with energy. However, there are no experimental or theoretical values available for a comparison. The computed vector decay constant ( f V ) and P-wave decay constant ( f i ) of the D meson are listed in Tables 12 and 13 along with other model predictions. Another important property of the D meson studied in the present case is the leptonic decay widths. The present branching ratios for D → τν τ (9.73 × 10 −4 ) and D → μν μ (3.846 × 10 −4 ) are in accordance with the exper-   Table 7. The large experimental uncertainty in the electron channel makes it difficult to reach any reasonable conclusion. The Cabibbo favoured nonleptonic branching ratios BR (D 0 → K − π + ) and BR (D 0 → K + π − ) obtained, respectively, as 4.071 % and 1.135 × 10 −4 , are in very good agreement with experimental values of 3.91 ± 0.08 % and (1.48 ± 0.07) × 10 −4 , [76].
We obtained the CP-violation parameter in mixing |q/ p| (0.9996) in this case, and the D 0 andD 0 decays show no evidence for CP violation and provides the most stringent bounds on the mixing parameters. The mixing parameter x q , y q , and mixing rate (R M ) are in very good agreement with BaBar, Belle and other Collaborations, as shown in Table  9. However, due to a larger uncertainty in the experimental values it is difficult for us to draw a conclusion on this mixing parameter. Thus, the present study of the mixing parameters of the neutral open-charm meson is found to be one of the successful attempts to extract the effective quark-antiquark interaction in the case of heavy-light flavour mesons. Thus the present study is an attempt to indicate the importance of spectroscopic (strong interaction) parameters in the weak decay processes.
Finally we look forward to see future experimental support of many of our predictions on the spectral states and decay properties of the open-charm meson.