New physics effects in charm meson decays involving c→ul+l-(li∓lj±)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c \rightarrow u l^+ l^- (l_i^{\mp } l_j^{\,\pm \,})$$\end{document} transitions

We study the effect of the scalar leptoquark and Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} boson on the rare decays of the D mesons involving flavour changing transitions c→ul+l-(li∓lj±)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c \rightarrow u l^+ l^- (l^{\mp }_i l^{\,\pm \,}_j)$$\end{document}. We constrain the new physics parameter space using the branching ratio of the rare decay mode 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 \rightarrow \mu ^+ \mu ^-$$\end{document} and the D0-D¯0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^0 - {\bar{D}}^0$$\end{document} oscillation data. We compute the branching ratios, forward–backward asymmetry parameters and flat terms in D+(0)→π+(0)μ+μ-\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 \pi ^{+(0)} \mu ^+ \mu ^-$$\end{document} processes using the constrained parameters. The branching ratios of the lepton flavour violating D meson decays, such as D0→μe,τe\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 \mu e, ~\tau e$$\end{document} and D+(0)→π+(0)μ-e+\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 \pi ^{+(0)} \mu ^- e^+$$\end{document} are also investigated.


Introduction
The rare B and D meson decay processes driven by a flavour changing neutral current (FCNC) transitions constitute a subject of great interest in the area of electroweak interactions and provide an excellent testing ground to look for new physics beyond the standard model (SM). The FCNC decays are highly suppressed in the SM and occur only at one-loop level. Of particular interest among the FCNC decays are the rare semileptonic B meson decays involving the transitions b → sl + l − , where several anomalies at the level of few sigma have been observed recently in the LHCb experiment [1][2][3][4]. To complement these results, efforts should also be made towards the search for new physics signal in the up quark sector, mainly in the rare charm meson decays involving c → ul + l − quark level transitions. Recently the LHCb experiment has searched for the branching ratio of the lepton flavour violating (LFV) D 0 → μ ∓ e ± decays and put the limit as BR(D 0 → μ ∓ e ± ) < 1.3 × 10 −8 [5] at 90% confidence level (CL). On the other hand, both the Belle and the BaBar experiments have reported significant deviations on the measured branching fractions ofB → D ( * ) τ ν τ processes from the corresponding SM predictions. The ratio of a e-mail: rmsp@uohyd.ernet.in these branching fractions, the so-called R(D ( * ) ), defined as R(D ( * ) ) = BR(B → D ( * ) τ ν τ )/BR(B → D ( * ) lν l ), where l = e, μ, exceed the SM prediction by 3.5σ [6], thus opening an excellent window to search for new physics (NP) in the up quark sector.
Mixing between a neutral meson and its anti-meson with a specific flavour provides an useful tool to deal with problems in flavour sector. For example, in the past the K 0 −K 0 and B 0 −B 0 oscillations, involving mesons made of up-and down-type quarks, have provided information as regards the charm and top quark mass scale, much before the discovery of these particles in the collider. On the other hand, the D 0 − D 0 system involves mesons with up-type quarks and in the SM the mixing rate is sufficiently small, so that the new physics component might play an important role in this case. The mixing parameters required to describe the D 0 −D 0 mixing are defined by x = M/ and y = /2 , where M ( ) is the mass (width) difference between the mass eigenstates.
In this paper, we focus on the analysis of rare charm meson decays induced by c → uμ + μ − and c → uμ ∓ e ± FCNC transitions. We calculate the branching ratios, forwardbackward asymmetry parameters and the flat terms in D +(0) → π +(0) μ + μ − processes both in the scalar leptoquark (LQ) and the generic Z model. These processes suffer from resonance background through c → u M → ul + l − , where M denotes η ( ) (pseudoscalar), ρ, φ, ω (vector) mesons. However, to reduce the background coming from these resonances, we work in the low and high q 2 regimes, i.e., q 2 ∈ [0.0625, 0.275] GeV 2 and q 2 ∈ [1. 56,4.00] GeV 2 , which lie outside the mass square range of the resonant mesons. We also compute the branching ratios of lepton flavour violating D 0 → μe, τ e and D +(0) → π +(0) μ − e + processes. These LFV processes have negligible contributions from the SM, as they proceed through the box diagrams with tiny neutrino masses in the loop. However, they can occur at tree level in the LQ and Z models and are expected to have significantly large branching ratios. Leptoquarks are hypothetical colour triplet bosonic particles, which couple to quarks and leptons simultaneously and contain both baryon and lepton quantum numbers. It is interesting to study flavour physics with leptoquarks as they allow quark-lepton transitions at tree level, thus explaining several observed anomalies, e.g., the lepton non-universality (LNU) parameter R K | q 2 ∈[1,6] GeV 2 = BR(B → K μ + μ − )/BR(B → K e + e − ) in rare B decays. The existence of the scalar leptoquark is predicted in the extended SM theories, such as grand unified theory [7][8][9][10][11], the Pati-Salam model, the extended technicolour model [12,13] and the composite model [14]. In this work, we consider the model which conserves baryon and lepton numbers and does not allow proton decay. Here we would like to see how this model affects the leptonic and semieptonic decays of the D 0 meson induced by c → ul + l − transitions. The phenomenology of scalar leptoquarks and their implications to the B and D sector has been extensively studied in the literature .
The Z boson is a colour singlet vector gauge boson and it is electrically neutral in nature. By adding an additional U (1) gauge symmetry, the new Z gauge boson could be naturally derived from the extension of electroweak symmetry of the SM, such as superstring theories, grand unified theories and theories with large extra dimensions. The processes mediated via c → u FCNC transitions could be induced by the generic Z model at tree level. The theoretical framework of the heavy new Z gauge boson has been studied in the literature [52,53,56,57]. In this paper, we investigate the Z contribution to the rare D 0 meson decay processes within the parameter space constrained by D 0 −D 0 mixing and D 0 → μ + μ − processes.
The paper is organized as follows. In Sect. 2, we discuss the effective Hamiltonian describing C = 1 transitions i.e., c → ul + l − , and C = 2 transition, which is responsible for D 0 −D 0 mixing. The new physics contribution to c → u transitions and the constraint on leptoquark couplings from the D 0 −D 0 oscillation and the process D 0 → μ + μ − are discussed in Sect. 3. We calculate the constraint on Z couplings from D 0 −D 0 mixing and leptonic D 0 → μ + μ − decays in Sect. 4. In Sect. 5, we compute the branching ratios, forward-backward asymmetry parameters and the flat terms of the process D +(0) → π +(0) μ + μ − in both these models. The lepton flavour violating D +(0) → π +(0) μ − e + and D 0 → μe, τ e processes are discussed in Sects. 6 and 7. Finally we summarize our findings in Sect. 8.

Effective Hamiltonian for C = 1 and C = transitions
Though the rare charm decays are affected by large nonperturbative effects, the short distance structure of FCNC transitions can be investigated well theoretically. The change in charm quantum number for rare FCNC charm meson decays is either of two units or one unit, and hence, they involve either C = 2 or C = 1 transitions. The D 0 −D 0 mixing takes place via a C = 2 transition and the decay processes with C = 1 transitions are c → ul + l − and c → uγ . If we integrate out the heavy degrees of freedom associated with the new interactions at a scale M, an effective Hamiltonian in the form of a series of operators of increasing dimensions can be obtained. However, the operators of dimension d = 6 have important contributions to charm meson decays or mixing. In general, one can write the complete basis of these effective operators in terms of chiral quark fields for both D 0 −D 0 mixing and the process D 0 → l + l − as [52,53] where G has inverse-mass squared dimensions, C i are the Wilson coefficients. 1 The effective operators for D 0 −D 0 mixing at the heavy mass scale M are given by [52,53] where q L(R) = L(R)q are the chiral quark fields with L(R) = (1∓γ 5 )/2 as the projection operators.
In the standard model, the effective weak Hamiltonian for the c → u transitions at the scale μ = m c can be written as the sum of three contributions [49,54,55], where λ q = V uq V * cq is the product of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The explicit form of H peng , which basically is responsible for the c → ul + l − transition is given by where G F is the Fermi constant, the C i 's are the Wilson coefficients evaluated at the charm quark mass scale (μ = m c ) at Next-Next-to-Leading-Order (NNLO) [58]. We use the two loop result of Ref. [59] for the C eff 7 (m c ) Wilson coefficients, V * cb V ub C eff 7 = V * cs V us (0.007 + 0.020i)(1 ± 0.2) and the corresponding effective operators for the c → ul + l − transitions are given by [49] O ( ) The contributions from the primed operators as well as the scalar, pseudoscalar and tensor operators are absent in the SM and arise only in scenarios beyond the standard model. The renormalization group running does not affect the O 10 operator, i.e., C 10 (m c ) = C 10 (M W ) ∝ (m 2 d,s /m 2 W ) and, hence, the Wilson coefficient C 10 is negligible in the SM. In this work we consider the baryon number conserving X = (3, 2, 7/6) scalar leptoquark which induces the interaction between the up-type quarks and charged leptons and, thus, contributes to the semileptonic decay amplitudes. The interaction Lagrangian of the X = (3, 2, 7/6) scalar leptoquark with the SM bilinears is given by [49,50] where˜ = iτ 2 * represents the conjugate state. The transition of weak basis to mass basis divides the Yukawa couplings to two part of couplings pertinent for the upper and lower doublet components. The left-handed quark and lepton doublets are represented by Q and L and u R (l R ) is the right-handed quark (charged-lepton) singlet. We use the basis where CKM and PMNS rotations are assigned to down-type quarks and neutrinos, i.e., Here Y L and Y R are the leptoquark couplings in the mass basis of the up-type quarks and charged leptons. Now writing the leptoquark doublets in terms of its components as = ( (5/3) , (2/3) ) T , where the superscripts denote the electric charge of the LQ components and expanding the terms in Eq. (6), one can obtain the interaction Lagrangian for different components of LQs given as [50] Thus, one can see from (7) where α em is the fine structure constant. After having an idea about the new Wilson coefficients, we now proceed to constrain the combination of LQ couplings using the experimental data on D 0 −D 0 mixing and the process D 0 → l + l − , where l = μ, e.
3.1 Constraint on leptoquark couplings from D 0 −D 0 mixing In the standard model, D 0 −D 0 mixing proceeds through the box diagrams with an internal down-type quarks and W -boson exchange and the weak interaction boxes are suppressed due to the GIM mechanism because of the smallness of the down quark mass in comparison to the weak scale. In the LQ model, there will be a contribution to the D 0 −D 0 mass difference from the box diagrams with the leptoquark and leptons flowing in the loop. Since the SM contribution to the mass difference is very small, we consider its value to be saturated by new physics contributions. Furthermore, the couplings to the left-handed quarks are considered to be zero in order to avoid strict constraints in the down-type quark sector. Thus, considering only right-handed couplings, one can write the effective Hamiltonian due to the leptoquark X (3, 2, 7/6) and charged lepton/neutrinos in the loop as [25,26] where the first term is due to the charged lepton and second term is due to neutrinos in the loop (ignoring the effect of neutrino mixing). The loop function I (x) is given as which is very close to 1, i.e., I (0) = 1. Using the relation we obtain the contribution due to leptoquark exchange as Since M D = 2|M 12 |, we get where l denotes the charged-lepton flavours. In our analysis, the mass of the D 0 meson is taken from [60] [63]. To obtain the bound on the leptoquark coupling, we assume that the individual leptoquark contribution to the mass difference does not exceed the 1σ range of the experimental value. Since we are interested in obtaining the bounds on the Y R μc Y R * μu couplings, here we assume that the leptoquark has dominant coupling to muons and its coupling to electron or tau is negligible. The SM contribution to the mass difference is very small and hence can be neglected. The corresponding experimental value is given by [60] Now comparing the mass difference with the 1σ range of experimental data, the bound on leptoquark coupling for a TeV scale LQ is given by which can be translated with Eq. (8) to give the constraint on the new Wilson coefficients

Constraint from
The rare leptonic D 0 → μ + μ − (e + e − ) processes, mediated by the FCNC transitions c → ul + l − at the quark level, are highly suppressed in the SM due to a negligible C 10 Wilson coefficient and also suffer from CKM suppression. These processes occur only at one-loop level and are considered as some of the most powerful channels to constrain the new physics parameter space in the charm sector. Analogous to the leptonic B meson decay processes, the only non-perturbative quantity involved is the decay constant of the D meson, which can be reliably calculated using non-perturbative methods such as QCD sum rules, lattice gauge theory and so on. The branching ratio of the process D 0 → l + l − is given by [48,49] The process D 0 → μ + μ − has a dominant intermediate γ * γ * state in the SM, which is electromagnetically converted to a μ + μ − pair. After including the contribution of the γ * γ * intermediate state, the predicted branching ratio of this process is BR(D 0 → μ + μ − ) 2.7 × 10 −5 × BR(D 0 → γ γ ) [64]. Using the upper bound BR(D 0 → γ γ ) < 2.2 × 10 −6 at 90% CL reported in [65], the estimated limit on the branching ratio is BR D 0 → μ + μ − SM 10 −10 [49]. The present experimental limits on the branching ratios of the dileptonic decays of the D meson are [60] BR Using the above experimental bounds, the constraint on the leptoquark coupling can be obtained by imposing the condition that the individual leptoquark contribution to the branching ratio does not exceed the experimental limit. In this analysis, we neglect the new physics contribution to the Wilson coefficient C 10 , as the scalar and pseudoscalar Wilson coefficients will be dominating due to the additional multiplication factor M D /m l as noted from Eq. (17). Now, redefining the Wilson coefficients bỹ we show in Fig. 1, the allowed region in theC whereas the bounds obtained from the process D 0 → e + e − are rather weak, i.e., It is obvious that the bounds obtained in Eqs. (20) and (21) could not give us proper information about the bounds on individualC coefficients. Therefore, we consider only one Wilson coefficient at a time to extract the upper bound on individual coefficients. In Table 1, we report the constraint on theC L Q S,P Wilson coefficients obtained from Table 1 The allowed values of the Wilson coefficients obtained from the upper bound of the process D 0 → μ + μ − (e + e − ). The constraint on theC L Q i coefficients can also be applicable toC Wilson coefficients are given in Table 1.
In order to evade the strict bounds in the down-type quark sector, we consider the leptoquark couplings to the lefthanded quarks (Y L ) as zero. Therefore, the only contribution to the rare charm decays comes from theC 9 = −C 10 Wilson coefficients, which are related to the right-handed quark couplings. Now, to include the (pseudo)scalar and (pseudo)tensor Wilson coefficients and to extract the respective upper bound complying with the constraints from B and K physics, we consider a numerically tuned example as discussed in [49]. We assume that the Y R coupling is perturbative, i.e., |Y R | < √ 4π . In particular, we consider a large value for Y R cμ coupling, e.g., Y R cμ = 3.5. We compute the bound on Y R uμ coupling by using the constraint on theC L Q

10
Wilson coefficients from D 0 → μ + μ − process, which is found to be comparatively small, Y R uμ < 8.76 × 10 −3 . Now we instigate a non-zero coupling to the left-handed quark Y L uμ , which along with the large Y R cμ coupling provides nonzero values for the C S,P and C T,T 5 coefficients. However, the process D 0 → μ + μ − imposes a strong bound on the coefficient C S , which, together with the large Y R cμ coupling, limits the left-handed coupling to Y L μu < 1.14 × 10 −3 . Thus, from the above discussion we observe that Our predicted bounds on the leptoquark coupling are in agreement with the constraints obtained in Refs. [48,49] and also with the constraints obtained from B, K physics [66].

New physics contribution in Z model
The new heavy Z gauge boson can exist in many extended SM scenarios and can mediate the FCNC transitions among the fermions in the up quark sector at tree level. The most general Hamiltonian for c → u transition in the Z model is given as [53] H FCNC Analogously, one can write the Hamiltonian for the leptonic sector H l Z as Here g Z i and g Z i are the couplings of the Z boson with the quarks and leptons, respectively, where i = 1 or 2 for the Z μ vector boson coupled to left-handed or right-handed currents.
After having obtained the possible Z couplings with quarks and leptons, we proceed to constrain the new parameter space using the results from charm sector, e.g., the experimental data on D 0 −D 0 mixing and the branching ratios of D 0 → l + l − processes. The constraint on the coupling of Z with the leptonic part is obtained from the upper limit on the branching ratio of the lepton flavour violating τ (μ) − → e − e + e − processes.

Constraint from D 0 −D 0 mixing
In this subsection, we calculate the constraint on the Z couplings from the mass difference of the charm meson mass eigenstates, which characterizes the D 0 −D 0 mixing phenomenon. The D 0 −D 0 oscillation arises from the | C = 2| transition that generates off-diagonal terms in the mass matrix for D 0 andD 0 mesons. The mass difference of D 0 −D 0 mixing at the scale μ = m c is given by [53] M (Z ) At the charm mass scale, the Wilson coefficients in terms of the Z couplings are expressed as where r (m c , M Z ) is the RG factor at the heavy mass scale and r (m c , M Z ) = 0.72 (0.71) for Z mass, M Z = 1(2) TeV [52]. Now we consider two possible cases to constrain the couplings g Z 1 and g Z 2 . One comes with only left-handed coupling present, i.e., (g Z 2 = 0) and the second with both left-handed and right-handed couplings present with equal strengths (g Z 1 = g Z 2 = g Z ). Here, we make the simple assumption that the NP part dominates over the SM contribution in D 0 −D 0 mixing. Thus, for the first case, on substitution of g Z 2 = 0 in Eqs. (25) and (26), the mass difference becomes Now varying the mass difference M D within its 1σ allowed range [60], we obtain and for a representative Z mass M Z = 1 TeV, the value of the coupling is found to be Analogously for the second case, i.e., g Z 1 = g Z 2 = g Z , the constraint is obtained.

Constraint from
The effective Hamiltonian for the process D 0 → l + l − in the Z model is given by [53] H Z c→u where the operators Q 1,2 are and Q 6,7 can be obtained from Q 1,2 by the substitutions of q L → q R and q R → q L . Comparing Eq. (31) with the SM effective Hamiltonian (4) yields the additional contributions to the Wilson coefficients C ( )Z 9,10 as The branching ratio for the process D 0 → μ + μ − in the Z model is given as [53] BR For simplicity, we consider here g Z 2 = 0. Now considering the couplings of the Z boson to the final leptons as the same form as the SM-like diagonal couplings of Z boson to leptons as discussed in [53], i.e., where g is the gauge coupling of Z boson and θ W is the Weinberg mixing angle. Now using the experimental upper limit on the branching ratio BR(D 0 → μ + μ − ) < 6.2 × 10 −9 [60], we obtain For M Z = 1 TeV, the constraint on the g Z 1 coupling is which is rather weak compared with the constraint obtained from D 0 −D 0 mixing. It should be noted that the constraint on Z couplings from the D 0 → μ + μ − decay process and the D 0 −D 0 mixing data have been computed in [53]. Similarly, the constraint on the couplings from D 0 −D 0 oscillation are obtained in Ref. [52]. We found that our constraints are consistent with the above predictions, if we use the updated values of various input parameters.

Constraints on
Considering only the left-handed coupling to the Z boson, the branching ratio of the process μ − → e − e + e − in the Z model is given by [67,68] BR μ − → e − e + e − = τ μ m 5 μ 768π 3 where we have explicitly shown the indices in the couplings.
The experimental upper limit on the branching ratio of this mode is BR(μ − → e − e + e − ) < 10 −12 [60]. For the analysis, we use the mass and lifetime of the muon from [60] and consider the coupling g L ee = g Z 1 as SM-like with its value as presented in Eq. (35). Thus, using the experimental upper limit, we get the bound on the g L μe coupling as Analogously using the branching ratio of τ − → e − e + e − process, BR(τ − → e − e + e − ) < 2.7 × 10 −8 [60], the constraint on the lepton flavour violating g L τ e coupling is found to be |g L τ e | < 0.02

D +(0) → π +(0) μ + μ − process
In this section, we study the rare semileptonic decay process D + → π + μ + μ − , which is mediated by the quark level transition c → uμ + μ − and constitutes a suitable tool to search for new physics. The dominant resonance contributions come from the φ, ρ and ω vector mesons and the effects of η ( ) mesons are comparatively negligible. These decay modes have recently been studied in Refs. [48,49,61] in various new physics scenarios and it is found that the model with scalar/vector leptoquarks and the minimal supersymmetric model with R-parity violation can give significant contributions. The matrix elements of various hadronic currents between the initial D meson and the final π meson can be parametrized in terms of three form factors f 0 , f + and f T [49]: π(k)|ūσ μν (1 ± γ 5 )c|D( p D ) where p D and k are the four momenta of the D and π mesons, respectively, and q = p D − k is the momentum transfer. The form factors for D 0 → π 0 are scaled as f i → f i / √ 2 by isospin symmetry. For the q 2 -dependence of the form factors, we use the parameterization from Refs. [69,70], as where x = q 2 /m 2 pole with m pole = 1.90(8) GeV, a = 0.28 (14) and b = 1.27(17) are the shape parameters [49] measured from D → πlν decay process and f + (0) = 0.67(3) [71]. The parameters in the f T form factor are (4) and a T = 0.18(16) [70]. Thus, one can write the transition amplitude for the process D + → π + μ + μ − as [49,72] where θ is the angle between the D meson and the negatively charged lepton in the rest frame of the dilepton. The functions V , A, S and P are defined in terms of the Wilson coefficients as ) , Here C ( )N P 9,10 , C ( )N P S,P and C N P T,T 5 are the new Wilson coefficients arising from either the scalar leptoquark model or the generic Z model. Using Eq. (44), the double differential decay distribution with respect to q 2 and θ , for the lepton flavour l is given by [49,72] d 2 l dq 2 dcos θ = a l (q 2 ) + b l (q 2 ) cos θ + c l (q 2 ) cos 2 θ, (46) where and Thus, the branching ratio is given by The forward-backward asymmetry ( A F B ) is another useful observable to look for new physics; it is defined by [49] Since the coefficient b l depends only on the scalar and pseudoscalar Wilson coefficients the forward-backward asymmetry is zero in SM. However, the additional new physics contribution can give a non-zero contribution to the forwardbackward asymmetry parameter. Another interesting observable is the flat term, defined as [72] F l H = q 2 max q 2 min dq 2 (a l + c l ) where the uncertainties get reduced due to the cancelation between the numerator and denominator. For numerical evaluation, we take the particle masses and the lifetime of the D meson from [60]. For the CKM matrix elements, we use the Wolfenstein parametrization with values A = 0.814 +0.023 −0.024 , λ = 0.22537 ± 0.00061, ρ = 0.117 ± 0.021 andη = 0.353 ± 0.013 [60]. With these input parameters, we compute the resonant/non-resonant branching ratios of the process D + → π + μ + μ − by integrating the decay distribution with respect to q 2 . We parametrize the contributions from the resonances with the Breit-Wigner shapes for C 9 → C res 9 , for ρ, ω, φ (vector) and C P → C res P for η ( ) (pseudoscalar) mesons [48,49]: Here m M ( M ) denotes the mass (total decay width) of the resonant state M, where M corresponds to η ( ) , ρ, ω, φ mesons. With the approximation of BR(D + → π + M(→ μ + μ − )) BR(D + → π + M)BR(M → μ + μ − ) and considering the experimental upper bound from [60], the magnitudes of the Breit-Wigner parameters are given by [48] a φ = 0.24 +0.05 −0.06 GeV 2 , a ρ = 0.17 ± 0.02 GeV 2 , a ω = a ρ /3, a η = 0.00060 +0.00004 −0.00005 GeV 2 , a η ∼ 0.0007 GeV 2 .
The detailed procedure of SM resonant contributions to the D + → π + μ + μ − process can be found in [48,49,70]. In Fig. 2, we show the q 2 variation of the branching ratio of D + → π + μ + μ − process including the resonant contribution in the SM. The band in the figure is due to the uncertainties associated with the a M parameters as given in (54) and the random variation of relative phases within −π and π . For simplicity we have assumed the same phase for all the resonances. From the figure, one can observe that in the low and high q 2 regions the long distance resonant contributions are approximately one order of magnitude below the current experimental sensitivity, and hence these regions are suitable to look for new physics beyond the SM. Thus, both in the SM and in the leptoquark and Z models, we study the process D + → π + μ + μ − only at the very low and high q 2 regimes. However, it should be emphasized that the uncorrelated variation of the unknown resonant phases affects the branching ratio in the low q 2 region significantly, which makes it quite difficult to infer the possible role of new physics. Fig. 2 The resonant contributions to the branching ratio of D + → π + μ + μ − in the SM. The band arises due to the uncertainties in Breit-Winger parameters and the variation of relative phases. The horizontal black line represents the experimental upper bound from [60] With all the input parameters from [60] along with the SM Wilson coefficients [58,59], we present in Table 2 the predicted values of the branching ratios for the D +(0) → π +(0) μ + μ − processes by integrating the decay distribution in low and high q 2 bins. Here we have used the q 2 regimes as q 2 ∈ [0.0625, 0.275] GeV 2 and q 2 ≥ 1.56 GeV 2 to reduce the background coming from the dominant resonances. The theoretical uncertainties in the SM are associated with the lifetime of the D meson, the CKM matrix elements and the hadronic form factors. In Fig. 3, the variation of the SM branching ratios of the process D + → π + μ + μ − in the very low and high q 2 regimes are shown in red dashed lines and the green bands represent the SM theoretical uncertainties. Now using the constraint on the leptoquark parameter space obtained in Sect. 3, we show in Fig. 3, the q 2 variation of the branching ratio of the process D + → π + μ + μ − in low q 2 (left panel) and high q 2 (right panel) regimes both in the scalar leptoquark and Z models. Here the orange (blue) band represents the contributions from the scalar leptoquark (Z ) model. The 90% CL experimental upper bounds on the branching ratios from [73], are shown in thick black lines. In Table 2, we present the integrated branching ratios of the process D +(0) → π +(0) μ + μ − es in both the low and the high q 2 regions in the leptoquark and Z models. We find that the predicted branching ratios in the leptoquark model have significant deviations from the corresponding SM values due to the effect of the scalar leptoquark and are well below the experimen- Table 2 The predicted branching ratios for D +(0) → π +(0) μ + μ − processes in both the low q 2 and the high q 2 region in the scalar X (3, 2, 7/6) LQ and Z model. This also contains the resonant and non-resonant SM branching ratios Decay process  Fig. 3 The variation of branching ratio of D + → π + μ + μ − with respect to low q 2 (left panel) and high q 2 (right panel). The orange bands represent the contributions from the scalar leptoquark, the blue bands are due to the Z contributions, the red dashed lines are for the non-resonant SM and the cyan bands are for the resonant SM. The green bands stand for the theoretical uncertainties from the input parameters in the SM. The solid black line denotes the 90% CL experimental upper limit [73] tal upper limits. However, the effect of the Z boson to the branching ratios of D +(0) → π +(0) μ + μ − processes is very marginal. In the leptoquark model, the variation of the forwardbackward asymmetry for the process D + → π + μ + μ − in low q 2 (left panel) and high q 2 (right panel) is presented in Fig. 4. The forward-backward asymmetry depends on the combinations of the C ( ) S and C T,T 5 Wilson coefficients, thus they have zero value in the SM. However, the additional contributions of C S,P Wilson coefficients due to scalar leptoquark exchange give non-zero contribution to the forward-backward asymmetry, though it is not so significant. The integrated forward-backward asymmetry for the process D + → π + μ + μ − is given as The Z model provides additional contributions only to the C 9,10 Wilson coefficients, and there are no new contributions to scalar or tensor terms. Thus, the forward-backward asymmetry vanishes in the Z model. In both the LQ and the Z model, the plot for flat term of D + → π + μ + μ − process with respect to low q 2 (left panel) and high q 2 (right panel) is given in Fig. 5 Fig. 4 The variation of the forward-backward asymmetry of D + → π + μ + μ − with respect to low q 2 (left panel) and high q 2 (right panel) in the scalar leptoquark model In addition to the leptoquark and Z models, the rare charm meson decays mediated by the c → u transitions have also been investigated in various new physics models such as the minimal supersymmetric standard model [61,64,74,75], the two Higgs doublet model [74], a warped extra dimensions model [75] and the up vector like quark singlet model [76,77]. In Refs. [48,49], the process D → πμ + μ − is studied in the context of both scalar and vector leptoquark models. Our results on predictions are found to be consistent with the literature.

D +(0) → π +(0) μ − e +
Since the individual lepton flavour number is conserved in the standard model, the observation of lepton flavour violation in the near future will provide unambiguous signal of new physics beyond the SM. The observation of neutrino oscillation implies the violation of lepton flavour in neutral sector and it is expected that there could be FCNC transitions in the charged lepton sector as well, such as l i → l j γ , l i → l j l klk , B → l ∓ i l ± j and B → K ( * ) l ∓ i l ± j etc. The LFV decay modes proceed through box diagrams with tiny neutrino masses in the loop, thus become very rare in the SM. However, these modes can occur at tree level in the leptoquark and Z models, thus can provide observable signature in the high luminosity experiments. In this section, we would like to study the lepton flavour violating semileptonic decay process D + → π + μ − e + . Due to the absence of intermediate states, these LFV processes have no long distance QCD contributions and dominant φ, ω resonance backgrounds. The general expression for the transition amplitude of D + → π + μ − e + process in a generalized new physics model, is given by where the functions F i , i = V, A, S, P are defined as , F A = K N P 10 + K N P 10 , Here the Wilson coefficients (K N P i ) involve the combination in Eq. (8). Now using Eq. (59), the differential decay distribution for the D + → π + μ − e + process with respect to q 2 and cos θ (θ is the angle between the D and μ − in the μ − e rest frame) is given as where and For numerical estimation in the leptoquark model, we use the constrained leptoquark couplings obtained from the process D 0 → μ + μ − and assume that the coupling between different generation of quarks and leptons follow the simple scaling laws, i.e. Y As discussed in [47,78,79], such pattern of ansatz can explain the decay widths of radiative LFV decay μ → eγ . Now using such ansatz, the variation of branching ratio with respect to q 2 for D +(0) → π +(0) μ − e + process in the leptoquark model is shown in left panel of Fig. 6 and the corresponding integrated value is given in Table 3. In this mode, the forward-backward asymmetry depends on K ( )N P 9,10 Wilson coefficients which give non-zero contribution. The left panel of Fig. 7 shows the q 2 variation of the forward-backward asymmetry and the corresponding integrated value is found to be (0.039 → 0.047). The variation of the flat term with respect to q 2 is presented in the left panel of Fig. 8 and the integrated value is (0.137 → 0.33).
For the Z model, we consider the constraint on the coupling of Z boson to the quarks, obtained from the D 0 −D 0 mixing and the process D 0 → μ + μ − as given in Eqs. (29) and (37). For the lepton flavour violating coupling, the constraint is taken from the process μ − → e − e + e − , as discussed in Sect. 4. Thus, using Eqs. (29), (37) and (39), the predicted branching ratio of D +(0) → π +(0) μ − e + process in the Z model is given in Table 3 and the q 2 variation of the process D + → π + μ − e + is shown in Fig. 6 (right panel). The forward-backward asymmetry variation is shown in right panel of Fig. 7 and the predicted value is −1.15 × 10 −3 , which is very small. In Fig. 8 (right panel), we show the plot for q 2 variation of the flat term and the integrated value is 0.158.

D 0 → μ − e + (τ − e + ) LFV decay process
Recently LHCb put the upper limit on branching ratio of the D 0 → μ ∓ e ± lepton flavour violating decay mode [5]: Neglecting the mass of electron, the branching ratio of the process D 0 → μ − e + is given by [48] BR(D 0 → μ − e + ) Fig. 7 The variation of forward-backward asymmetry of the LFV process D + → π + μ − e + in the leptoquark model (left panel) and generic Z model (right panel) with respect to q 2 We use the scaling ansatz as discussed in the previous section to compute the required leptoquark coupling for D 0 → μ − e + process and the predicted branching ratio is found to be BR(D 0 → μ − e + ) = (3.18 − 4.8) × 10 −11 .
The predicted branching ratio is although small, but can be searched at LHCb experiment. The exploration/observation of this decay mode would definitely shed some light in the leptoquark scenarios. Similarly using the new Wilson coefficient generated via leptoquark exchange, the branching ratio for the process D 0 → τ − e + is found to be BR(D 0 → τ − e + ) = (2.84 − 9.75) × 10 −14 .
However, there is no experimental observation of the lepton flavour violating process D 0 → τ − e + . The constraint on Z coupling to tau and electron is obtained from the process τ − → e − e + e − . Using Eq. (40), the branching ratio for the process D 0 → τ − e + in Z model is given as BR(D 0 → τ − e + ) = (0.73 − 1.94) × 10 −15 .

Conclusion
In this paper we have studied the rare decays of the D meson in both scalar leptoquark and generic Z models. We have considered the simple renormalizable scalar leptoquark model with the requirement that proton decay would not be induced in perturbation theory. The leptoquark parameter space is constrained using the present upper limit on branching ratio of D 0 → μ + μ − process and the D 0 −D 0 oscillation data. For the Z model, the constraints on the Z couplings are obtained from the mass difference of D 0 −D 0 mixing, the process D 0 → μ + μ − and the lepton flavour violating τ − (μ − ) → e − e + e − processes. Using the constrained parameter space, we estimated the branching ratios, forward-backward asymmetry parameters and the flat terms in the D +(0) → π +(0) μ + μ − processes. The branching ratios in the LQ model are found to be ∼ O(10 −10 ), which are larger than the corresponding SM predictions in the very low and very high q 2 regimes. If these branching ratios will be observed in near future they would provide indirect hints of leptoquark signal. Furthermore, we estimated the branching ratios of lepton flavour violating D +(0) → π +(0) μ − e + and D 0 → μ(τ ) − e + processes, which are found to be rather small. We also estimated the forward-backward asymmetry parameter and the flat term for the LFV decays.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .