New Physics effects in charm meson decays involving $c \to u l^+ l^- (l_i^\mp l_j^\pm)$ transitions

We study the effect of scalar leptoquark and $Z^\prime$ boson on the rare decays of $D$ mesons involving flavour changing transitions $c \to u l^+ l^- (l^\mp_i l^\pm_j)$. We constrain the new physics parameter space using the branching ratio of the rare decay mode $D^0 \to \mu^+ \mu^-$ and the $D^0-\bar D^0$ oscillation data. We compute the branching ratios, forward-backward asymmetry parameters and flat terms in $D^{+(0)} \to \pi^{+(0)} \mu^+ \mu^-$ processes using the constrained parameters. The branching ratios of the lepton flavour violating $D$ meson decays, such as $D^0 \to \mu e, ~\tau e$ and $D^{+(0)} \to \pi^{+(0)} \mu^- e^+$ are also investigated.


I. INTRODUCTION
The rare B and D meson decay processes driven by 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]. To compliment 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 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 [2] at 90% confidence level (CL). On the other hand, both Belle and BaBar experiments have reported significant deviations on the measured branching fractions ofB → D ( * ) τ ν τ processes from the corresponding SM predictions. The ratio of 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σ [3], thus open 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 up of down type quarks, had provided information about the charm and top quark mass scale, much before the discovery of these particles in the collider. On the other hand, 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, forward-backward asymmetry parameters and the flat terms in D +(0) → π +(0) µ + µ − processes both in the scalar leptoquark (LQ) and generic Z model. These processes suffer from resonance back-ground through c → uM → 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 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 color 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 explain several observed anomalies, e.g., lepton non-universality (LNU) parameter R K | q 2 ∈[1,6] GeV 2 = BR(B → Kµ + µ − )/BR(B → Ke + e − ) in rare B decays. The existence of scalar leptoquark is predicted in the extended SM theories, such as grand unified theory [4,5], Pati-Salam model, extended technicolor model [6] and the composite model [7]. 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 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 [8][9][10][11][12][13][14][15][16][17].
The Z boson is a color singlet vector gauge boson and 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 generic Z model at tree level. The theoretical framework of the heavy new Z gauge boson has been studied in the literature [18,19,21]. 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 section II, 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 D 0 −D 0 oscillation and D 0 → µ + µ − process are discussed in section III. We calculate the constraint on Z couplings from D 0 −D 0 mixing and leptonic D 0 → µ + µ − decays in section IV. In section V, we compute the branching ratios, forwardbackward asymmetry parameters and the flat terms of D +(0) → π +(0) µ + µ − process in both these models. The lepton flavour violating D +(0) → π +(0) µ − e + and D 0 → µe, τ e processes are discussed in sections VI and VII. Finally we summarize our findings in section VIII.

II. EFFECTIVE HAMILTONIAN FOR ∆C = 1 AND ∆C = 2 TRANSITIONS
Though the rare charm decays are affected by large non-perturbative 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 or one unit, and hence, involve either ∆C = 2 or ∆C = 1 transitions. The D 0 −D 0 mixing takes place via ∆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 D 0 → l + l − process as [18,19] 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 [18,19] where q L(R) = L(R)q are the chiral quark fields with L(R) = (1 ∓ γ 5 )/2 as the projection operators. 1 We denote the Wilson coefficients for ∆C = 2 operators as c i and those for ∆C = 1 operators as C i through out this work In the standard model, the effective weak Hamiltonian for c → u transitions at the scale µ = m c , can be written as the sum of three contributions as [15,20] where λ q = V uq V * cq is the product of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The explicit form of H peng , which basically responsible for the c → ul + l − transition is given by where G F is the Fermi constant, C i 's are the Wilson coefficients evaluated at the charm quark mass scale (µ = m c ) at Next-Next-to-Leading-Order (NNLO) [22]. We use the two loop result of Ref. [23] 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 c → ul + l − transitions are given as [15] O ( ) The contributions from the primed operators as well as the scalar, pseudoscalar and tensor  (3, 2, 7/6) [15,16], only the leptoquark with multiplet (3, 2, 7/6) conserves both baryon and lepton numbers and thus, avoids rapid proton decay at the electroweak scale. Similarly out of the vector multiplets (3, 3, 2/3), (3, 1, 5/3), (3, 2, 1/6) and (3, 2, −5/6), only first two leptoquark states don't allow baryon number violation and can be considered to study the observed anomalies in flavour sector.
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 X = (3, 2, 7/6) scalar leptoquark with the SM bilinears is given by [15,16] L Thus, one can see from (7), that only ∆ (5/3) component mediates the interaction between up-type quarks and charged lepton. Now applying the Fierz transformation, we obtain additional contributions to the SM Wilson coefficients for c → uµ + µ − transition as [15] C LQ 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 D 0 → l + l − process, where l = µ, e.
A. 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 GIM mechanism because of the smallness of down-quark mass in comparison to the weak scale. In the LQ model, there will be 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 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 [10] H where the first term is due 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 D 0 meson is taken from [24], the value of the decay constant f D = 222.6 ± 16.7 +2.3 −2.4 MeV [26] and B D (3 GeV) = 0.757 (27)(4) [27]. To obtain the bound on the leptoquark coupling, we assume that individual leptoquark contribution to the mass difference does not exceed the 1σ range of the experimental value. Since we are interested to obtain the bounds on Y R µc Y R * µu couplings, here we assume that 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 [24] ∆M D = 0.0095 +0.0041 −0.0044 ps −1 .
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 Eqn. (8) to give the constraint on new Wilson coefficients as The rare leptonic D 0 → µ + µ − (e + e − ) processes, mediated by FCNC transitions c → ul + l − at the quark level, are highly suppressed in SM due to 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 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 D 0 → l + l − process is given by [14,15] The D 0 → µ + µ − process has dominant intermediate γ * γ * state in the SM, which is electromagnetically converted to a µ + µ − pair. After including the contribution of γ * γ * intermediate state, the predicted branching ratio of this process is BR CL reported in [29], the estimated limit on branching ratio is BR [15]. The present experimental limits on the branching ratios of dileptonic decays of D meson are [24] BR Using the above experimental bounds, the constraint on leptoquark coupling can be obtained by imposing the condition that 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 Eqn. (17). Now, redefining the Wilson coefficients as we show in Fig. 1, the allowed region inC LQ from Eqn. (8). From the figure, we found the allowed range for the above combinations of Wilson coefficients from D 0 → µ + µ − process as whereas the bounds obtained from D 0 → e + e − process is rather weak, i.e., It is obvious that the bounds obtained in Eqns. (20) and (21)  In Table I, we report the constraint onC LQ S,P Wilson coefficients obtained from the experimental bound on the branching fraction of D 0 → µ + µ − (e + e − ) process. The bounds onC LQ i Wilson coefficients will be same as those forC LQ i . If we impose chirality on scalar leptoquarks i.e., they couple to either left-handed or righthanded quarks, but not to both then the C  Table I.
In order to evade the strict bounds in the down type quark sector, we consider the leptoquark couplings to the left handed quarks (Y L ) as zero. Therefore, the only contribution to the rare charm decays comes fromC 9 = −C 10 Wilson coefficients, which are related to the right-handed quark couplings. Now to include (pseudo)scalar and (pseudo)tensor Wilson coefficients and to extract respective upper bound complying with the constraints from B and K physics, we consider a numerically tuned example as discussed in [15]. We assume 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 onC LQ 10 Wilson coefficients from D 0 → µ + µ − process, which is found to be small comparatively, Y R uµ < 8.76 × 10 −3 . Now we instigate a nonzero coupling to the left handed quark Y L uµ , which along with the large Y R cµ coupling provides nonzero values for C S,P and C T,T 5 coefficients. However, the D 0 → µ + µ − process imposes strong bound on C S coefficient, which together with large Y R cµ coupling, limits the left handed coupling as Y L µu < 1.14 × 10 −3 . Thus, from the above discussion we observe that Our predicted bound on leptoquark coupling are in agreement with the constraints obtained in Refs. [14,15] and also with the constraints obtained from B, K physics [30].

IV. 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 [19] 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 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 knowing 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 branching ratio of lepton flavour violating τ (µ) − → e − e + e − processes.
At the charm mass scale, the Wilson coefficients in terms of 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 [18].
Now we consider two possible cases to constrain the couplings g Z 1 and g Z 2 . One with only left handed coupling present, i.e., (g Z 2 = 0) and the second where both left handed and right handed couplings are present with equal strength (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, substitution of g Z 2 = 0 in Eqns. (25) and (26), the mass difference becomes Now varying the mass difference ∆M D within its 1σ allowed range [24], 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 obtained as The effective Hamiltonian for D 0 → l + l − process in the Z model is given by [19] H 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 Eqn. (31) with the SM effective Hamiltonian (4), yields the additional contributions to Wilson coefficients C ( )Z 9,10 as The branching ratio for D 0 → µ + µ − process in the Z model is given as [19] BR For simplicity, we consider here g Z 2 = 0. Now considering the couplings of Z boson to the final leptons as the same form as the SM-like diagonal couplings of Z boson to leptons as discussed in [19], i.e., where g is the gauge coupling of Z boson and θ W is the weak mixing angle. Now using the experimental upper limit on branching ratio BR(D 0 → µ + µ − ) < 6.2 × 10 −9 [24], we obtain For M Z = 1 TeV, the constraint on g Z 1 coupling is which is rather weak compared to 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 [19]. Similarly, the constraint on the couplings from D 0 −D 0 oscillation are obtained in Ref. [18]. We found that our constraints are consistent with the above predictions, if we use the updated values of various input parameters.

C. Constraints on
Considering only the left handed coupling to the Z boson, the branching ratio of µ − → e − e + e − process in the Z model is given by [31,32] BR µ − → e − e + e − = τ µ m 5 where we have explicitly shown the indices in the couplings. The experimental upper limit on branching ratio of this mode is BR(µ − → e − e + e − ) < 10 −12 [24]. For the analysis, we use the mass and lifetime of muon from [24] and consider the coupling g L ee = g Z 1 as SM-like with its value as presented in Eqn. (35). Thus, using the experimental upper limit, we get the bound on g L µe coupling as Analogously using the branching ratio of τ − → e − e + e − process, BR(τ − → e − e + e − ) < 2.7 × 10 −8 [24], the constraint on lepton flavour violating g L τ e coupling is found to be |g L τ e | < 0.02 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 are recently studied in Refs. [14,15,25] in various new physics scenarios and it is found that the model with scalar/vector leptoquarks and 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 [15] as 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. [33,34], 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 [15] measured from D → πlν decay process and f + (0) = 0.67(3) [35]. The parameters in the f T form factor are: x T = q 2 /M 2 D * , f T (0) = 0.46(4) and a T = 0.18 (16) [34]. Thus, one can write the transition amplitude for D + → π + µ + µ − process as [15,36] 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 Eqn. (44), the double differential decay distribution with respect to q 2 and θ, for the lepton flavour l is given by [15,36] where with and (49) Thus, the branching ratio is given by The forward-backward asymmetry (A F B ) is another useful observable to look for new physics, which is defined as [15] A . (51) Since the coefficient b l depends only on scalar and pseudoscalar Wilson coefficients the forward-backward asymmetry is zero in SM. However, the additional new physics contribution can give non-zero contribution to the forward-backward asymmetry parameter. Another interesting observable is the flat term, defined as [36] F where the uncertainties get reduced due to the cancelation between the numerator and denominator.
The detailed procedure of SM resonant contributions to D + → π + µ + µ − process can be found in [14,15,34]. In Fig. 2, we show the q 2 variation of 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)  With all the input parameters from [24] along with the SM Wilson coefficients [22,23], we present in Table II, the predicted values of 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 D meson, CKM matrix elements and the hadronic form factors. In Fig.   3, the variation of SM branching ratios of D + → π + µ + µ − process 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 section III, we show in Fig. 3, the q 2 variation of branching ratio of D + → π + µ + µ − process in low q 2 (left panel) and high q 2 (right panel) 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 [37] BR(D + → π + µ + µ − )| low q 2 < 2.0 × 10 −8 , are shown in thick black lines. In Table II, we present the integrated branching ratios of D +(0) → π +(0) µ + µ − processes in both the low and 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 scalar leptoquark and are well below the experimental upper limits. However the effect of Z boson to the branching ratios of D +(0) → π +(0) µ + µ − processes is very marginal.
In the leptoquark model, the variation of forward-backward asymmetry for D + → π + µ + µ − process 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 C ( ) S and C T,T 5 Wilson coefficients, thus 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 forwardbackward asymmetry, though it is not so significant. The integrated forward-backward asymmetry for D + → π + µ + µ − process 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 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. The predicted values in low q 2 range are and in the region of high q 2 F H | SM = 0.03 ± 0.005, In addition to the leptoquark and Z models, the rare charm meson decays mediating by the c → u transitions have also been investigated in various new physics models such as, Minimal Supersymmetric Standard Model [25,28,38,39], two Higgs doublet model [38], warped extra dimensions model [39] and the up vector like quark singlet model [40]. In the Ref. [14,15], the D → πµ + µ − process is studied in the context of both scalar and vector leptoquark models. Our predicted results are found to be consistent with the literature. The horizontal black line represents the experimental upper bound from [24].
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 The predicted branching ratios for D +(0) → π +(0) µ + µ − processes in both the low q 2 and high q 2 region in the scalar X(3, 2, 7/6) LQ and Z model. This also contains the resonant and nonresonant SM branching ratios.
Decay process 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 Here the Wilson coefficients (K N P i ) involve the combination of LQ couplings as Y L µc Y in Eqn. (8). Now using Eqn. (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 D 0 → µ + µ − process and assume that the coupling between different generation of quarks and leptons follow the simple scaling laws, i.e. Y As discussed in [13,41], 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 III. In this mode, the forward backward asymmetry depends on K ( )N P 9,10 Wilson coefficients which give nonzero 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  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 D 0 → µ + µ − process as given in Eqn. (29) and (37).
For the lepton flavour violating coupling, the constraint is taken from µ − → e − e + e − process, as discussed in section IV. Thus, using Eqn (29), (37) and (39), the predicted branching ratio of D +(0) → π +(0) µ − e + process in the Z model is given in Table III and the q 2 variation of D + → π + µ − e + process 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.
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.
However, there is no experimental observation of lepton flavour violating D 0 → τ − e + process.
The constraint on Z coupling to tau and electron is obtained from the τ − → e − e + e − process.

VIII. CONCLUSION
In this paper we have studied the rare decays of 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 Z couplings are obtained from the mass difference of D 0 −D 0 mixing, D 0 → µ + µ − process 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 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.