Prospects of discovering new physics in rare charm decays

The LHCb bounds on the branching ratio of the rare decay 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 constraints on the branching ratio of D+→π+μ+μ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \mu ^+ \mu ^-$$\end{document} in the nonresonant regions enable us to improve constraints on new physics contributions. Using the effective Lagrangian approach we determine the sizes of the Wilson coefficients allowed by the existing LHCb bounds on rare charm decays. Then we discuss contributions to rare charm meson decay observables in several models of new physics: a model with an additional spin-1 weak triplet, leptoquark models, Two Higgs doublets model of type III, and a Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z'$$\end{document} model. Here we complement the discussion by 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} oscillations data. Among the considered models, only leptoquarks can significantly modify the Wilson coefficients. Assuming that the differential decay width for D+→π+μ+μ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \mu ^+ \mu ^-$$\end{document} receives a NP contribution, while the differential decay width for D+→π+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^+ \rightarrow \pi ^+ e^+ e^-$$\end{document} is Standard Model-like, we find that lepton flavor universality can be violated and might be observed at high dilepton invariant mass.


Introduction
Processes with charmed mesons and top quarks offer an excellent opportunity to search for new physics (NP) in the up-type quark sector. In contrast to B meson physics, which is convenient to search for NP due to good exposure of the short-distance effects, charm quark systems are dominated by large long-distance quantum chromodynamics contributions. Such effects then screen the short-distance contributions of interest. Within the Standard Model (SM) the shortdistance physics in rare charm decays is strongly affected by the Glashow-Iliopoulos-Maiani (GIM) mechanism [1]. Namely, box or penguin diagram amplitudes get contributions from down-type quarks which are approximately massless from the weak scale perspective, and this warrants a very a e-mail: svjetlana.fajfer@ijs.si b e-mail: nejc.kosnik@ijs.si effective GIM cancellation. Flavor changing neutral current (FCNC) processes with charm mesons might change charm quantum number for two or one unit (| C| = 2 or | C| = 1 transitions). The | C| = 2 transition occurs in D 0 -D 0 oscillations and leads to strong constraints on NP from the measured observables as pointed out in [2,3]. There are two possibilities for NP in the | C| = 2 transition: the transition might occur at tree level, in which case a new neutral scalar or a vector boson possesses FCNC couplings to u and c quarks, or at loop level via NP degrees of freedom affecting the box diagrams. The processes with | C| = 1 on the quark level are c → uγ and c → u + − [4][5][6][7][8][9]. Both transitions can be approached in the familiar effective Lagrangian formalism [3]. Additional constraints on NP arise from the down-type quark sector whenever new bosons couple to lefthanded quark doublets [10,11]. Since NP is very constrained by the current experimental results coming from B and K physics [12] the only chance to observe NP in rare charm decays seems to be when new bosons are coupled to weak singlets. This then allows one to avoid the strong flavor constraints in the down-type quark sector.
On the experimental side the LHCb experiment succeeded to improve the bound on the rates of | C| = 1 decays by almost two orders of magnitude with respect to previous bounds. For the dileptonic decay the best bound to date is [13] BR(D 0 → μ + μ − ) < 7.6 × 10 −9 . (1) The above limit as well as other quoted limits in the following, unless stated otherwise, correspond to the 95 % CL upper bounds. In the decay D + → π + μ + μ − the LHCb experiment focused on kinematic regions of dilepton mass, q 2 = (k − + k + ) 2 , that are below or above the dominant resonant contributions due to vector resonances in the range m 2 ρ q 2 m 2 φ . The measured total branching ratio, obtained by extrapolating spectra over the resonant region, is [14] BR(D + → π + μ + μ − ) < 8.3 × 10 −8 , while separate branching fractions in the low-and high-q 2 bins were bounded as [14] 1 BR(π + μ + μ − ) I ≡ BR(D + → π + μ + μ − ) q 2 ∈[0.0625,0.276] GeV 2 < 2.5 × 10 −8 Motivated by these improved bounds we consider several NP models and either derive constraints on their flavor parameters and masses, or for the models that are severely bounded from alternative flavor observables (e.g. D 0 -D 0 mixing, K , or B physics), we comment on the prospects of observing their signals in rare charm decays. To this end, we use the effective Lagrangian encoding the short-distance NP contributions in a most general way. Namely, the experimental results (1) and (3) give us a possibility to constrain NP in c → u + − also in a model independent way. In the case of b → s + − transitions, LHCb has recently observed a large departure of the experimentally determined lepton flavor universality (LFU) ratio [15]. This value was found to be R LHCb K = 0.745 +0.090 −0.074 ± 0.036, lower than the SM prediction R SM K = 1.0003 ± 0.0001 [16]. This surprising result of LHCb indicates possible violation of LFU in the μ-e sector. Due to the importance of this result, we investigate whether analogous tests in the μ-e LFU can be carried out in c → u + − processes.
The outline of this article is as follows. In Sect. 2 we describe effective Lagrangian of | C| = 1 transition and determine bounds on the Wilson coefficients coming from the experimental limits on BR(D + → π + μ + μ − ) and BR(D 0 → μ + μ − ). Section 3 contains an analysis in the context of specific theoretical models of new physics, contributing to the c → u + − and related processes. Section 4 discusses lepton flavor universality violation. Finally, we summarize the results and present conclusions in Sect. 5.

Effective Hamiltonian for c → u + −
The relevant effective Hamiltonian at scale μ c ∼ m c is split into three contributions corresponding to diagrams with intermediate quarks q = d, s, b [9,17] where each of them is weighted by an appropriate combination λ q = V uq V * cq of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Virtual contributions of states heavier than charm quark are by convention contained within The operators appearing in the above Hamiltonian have thus enhanced sensitivity to new physics contributions: The chiral projectors are defined as P L ,R = (1 ∓ γ 5 )/2, F μν is the electromagnetic field strength tensor. For each of the operators O 7,9,10,S,P we introduce the corresponding counterpart O 7,9,10,S,P with opposite chiralities of quarks. Within the SM the Wilson coefficients C i result from the perturbative dynamics of the electroweak interactions and QCD renormalization. The latter effect determines the value of C 7 (m c ) by two-loop mixing with current-current operators and was found to be V * cb V ub C SM 7 = V * cs V us (0.007 + 0.020i)(1 ± 0.2) [4,8]. On the other hand the value of the C 9 Wilson coefficient was found to be small after including renormalization group running effects, as shown in [7] and confirmed in [6], while C 10 is negligible in the SM [18].
In order to analyze NP effects in D + → π + μ + μ − one needs to evaluate the hadronic transition matrix elements of currentsūγ μ P L ,R c andūσ μν P L ,R c. The standard parametrization expresses these matrix elements in terms of three form factors: π(k)|ūσ μν (1 ± γ 5 )c|D( p) where q = p − k is the dilepton four-momentum. For the f +,0 (q 2 ) form factors we use the Bečirević-Kaidalov (BK) parametrization [19]: with the shape parameters m pole and a determined by measurements of D → π ν decay spectra. We make an average of four experimental fits to the shape parameters, by taking as input the CLEO-c tagged [20] and untagged analysis [21], BES III [22], and Babar [23] results, all compiled by the HFAG [24]. The fitted shape parameters are m pole = 1.90(8) GeV and a = 0.28 (14). For the normalization of the form factor we rely on the lattice result f + (0) = 0.67(3) calculated by the HPQCD Collaboration [25]. The shape parameter b = 1.27 (17) has also been extracted in lattice simulations [26]. For the tensor current form factor we rely on the fit of lattice data to BK shape as in [26]: (4), and a T = 0.18 (16). Based on the effective Hamiltonian (4), the most general expression for the short-distance amplitude can be written as [27] Here θ is defined as the angle between the three-momenta of B and − in the rest frame of lepton pair whereas V, A, S, P, T , and T 5 are q 2 -dependent functions expressed in terms of hadronic form factors and Wilson coefficients, We have employed a shorthand notation The decay spectrum can be expressed in terms of q 2 -dependent angular coefficients as whereas the angular coefficients are The coefficients a and c enter then the q 2 distribution of branching ratio whereas b is proportional to the forwardbackward asymmetry: Contributions of the vector resonances ρ, ω, and φ, decaying to μ + μ − , is due to the first two terms in the effective Hamiltonian (4) and electromagnetic interaction. Effects of vector resonances to the spectrum can be treated assuming naïve factorization by adding a q 2 -dependent piece to C 9 that contains the vector current of leptons. Analogously, the scalar contribution of η feeds into C S . The procedure is described in detail in Ref. [28] for the contribution of D + → π + ρ 0 (ω) and updated for D + → π + φ → π + μ + μ − in Ref. [26]. The current experimental upper bound outside the resonance region indicates that the long-distance contribution is very suppressed. One might expect that at high invariant dilepton mass bin some excited states of vector mesons might give an additional long-distance contribution. However, it was shown in [29,30] that contributions of these states are negligible in comparison with the leading long-distance contributions. We parametrize the resonances with the Breit-Wigner shapes, The magnitude of unknown parameters a X (X = ρ, ω, φ, η), can be fitted to the measured resonant branching ratios, given in Table 1 [31]. The corresponding values of |a X | are given in the second row in Table 1. We treat the relative phases as free parameters. Alternatively, for the relative phases and magnitudes of a X one can use flavor structure arguments [18]. In the left-hand panel in Fig. 1 we present the long-distance contributions to the differential branching ratio for D + → π + μ + μ − as a function of dilepton invariant mass for a representative set of parameters |a X | from the 1σ region (Table 1) and random phases of a X . On the righthand panel in Fig. 1 we also indicate the interpretation of experimental upper bounds (3) in the case where the total amplitude would be constant, namely in the case where all angular coefficient functions a , b , c would be independent of q 2 . We also estimate the saturation of these bounds by the  (2) total resonant decay branching ratio and find for the lowand high-q 2 bin contributions to be smaller than 7.3 × 10 −9 and 5.3 × 10 −9 , respectively. On the other hand, the shortdistance contribution to the total branching ratio of the SM due to the quoted value of C 7 is of the order 10 −12 and thus negligible. The branching ratio for D 0 → μ + μ − can be written in its most general form as where the decay constant of a D meson, f D = 209 (3) MeV, has been averaged over N f = 2 + 1 lattice simulations [32][33][34]. In the SM this decay is dominated by the intermediate , and, together with the upper bound BR(D 0 → γ γ ) < 2.2 × 10 −6 at 90 % CL [35], this leads to the limit BR(D 0 → μ + μ − ) SM 10 −10 .

Constraints on the Wilson coefficients
In this section we interpret the experimental bounds in kinematical regions I and II given in Eq.
(3) as constraints on benchmark scenarios with NP contributions affecting individual Wilson coefficients. In the nonresonant regions of D + → π + μ + μ − the long-distance resonant contributions are one order of magnitude below the current experimental sensitivity. This allows us to saturate experimental results for the differential decay width distribution at the low/high dilepton invariant mass bins by the contributions of the effective Wilson coefficients. In Fig. 2   the time, where we have neglected the SM resonant contributions. A strong kinematical dependence of spectra in the cases of tensor and especially pseudoscalar coefficients suggests they will be better constrained in the high-than in the low-q 2 bin. For EM dipole and (axial-)vector interactions the enhancement at low-q 2 bin is hindered by relatively smaller phase space devoted to that bin. We allow only one Wilson coefficient at a time to have a real nonzero value and extract its upper bound. This is repeated for each choice of random phases and moduli of the Breit-Wigner parameters a η,ρ,ω,φ , where the latter are sampled in their 1σ regions (90 % CL bound for |a ω |); cf. Table 1. The most relaxed bound obtained in this way is then reported in Table 2, where we use the notationC i = V ub V * cb C i . At the same time the branching ratio of D 0 → μ + μ − can give bounds on the Wilson coefficients C 10 , C S , and C P . It turns out that the upper bound on BR(D 0 → μ + μ − ) is more restrictive for the C S,P,10 Wilson coefficients than any of the invariant dilepton mass bins of D + → π + μ + μ − .
The high invariant dilepton mass bin is more restrictive than the low dilepton invariant mass bin. Due to the parity conservation in D → π transition the bounds forC j , j = 7, 9, 10, S, P are the same as forC j .
In specific cases the angular distribution with respect to cos θ can be a good discriminant between the resonant and genuine short-distance contributions. It was shown that the forward-backward asymmetry (FBA) can be enhanced toward the larger end of the q 2 -spectrum in models with tensor and scalar Wilson coefficients (or pseudoscalar and pseudotensor) simultaneously present [36]. In principle such a scenario can be realized by a nonchiral leptoquark to be discussed in the following section. As a numerical example we chooseC S = 0.049, allowed by the D 0 → μ − μ + , and in additionC T = 0.2, which results in BR(D + → π + μ + μ − ) II < 10 −8 and is therefore hard to distinguish from the resonant background. On the other hand, the FBA Table 2 Maximal allowed values of the Wilson coefficient moduli, |C i | = |V ub V * cb C i |, calculated in the nonresonant regions of D + → π + μ + μ − in the low lepton invariant mass region (q 2 ∈ [0.0625, 0.276] GeV 2 ), denoted by I, in the high invariant mass region (q 2 ∈ [1.56, 4.00] GeV 2 ), denoted by II, and from the upper bound BR(D 0 → μ + μ − ) < 7.6 × 10 −9 [13]. The last row gives the maximal value for the case whereC 9 = ±C 10 . All the quoted bounds have been derived for real C i . The bounds forC i apply also to the chirally flipped coefficientsC j  in this case is strongly enhanced in the high-q 2 region, as shown in Fig. 3. We turn to the discussion of specific models in the next section.

Spin-1 weak triplet
Introducing an additional vector particle that transforms as a triplet under SU (2) L affects a plethora of flavor observables. It has been shown recently that a model of this type explains the current B sector anomalies (R D ( * ) , R K ) even in the scenario with U (2) q × U (2) flavor symmetry [11]. The relevant effective Lagrangian that follows from integrating out the vector triplet at tree level reads representing a contact interaction between vector currents of left-handed quark and lepton doublets: The indices of λ q and λ denote the mass-eigenstates of down-type quarks and charged leptons. The Hermitian matrices λ q, i j are conventionally normalized to λ q 33 = λ 33 = 1, with subleading entries in λ 3i , whereas entries involving only the first two generations of quarks are severely suppressed. This hierarchy is a direct consequence of the imposed flavor symmetry. The dominant contributions to the processes involving only the first two generations are induced by λ q bb and accompanying CKM rotations. Following [11] for the rare charm decays these are where the quark sector couplings originate dominantly from the CKM mixing, namely (V λ q V † ) uc ≈ V ub V cb * . A purely left-handed current effective interaction generates the following pair of Wilson coefficients: where R 0 = (g q g m 2 W )/(g 2 m 2 V ) is directly related to the LFU τ/ ratio R D ( * ) in semileptonic B meson decays whose experimental value requires R 0 = 0.14±0.04. The constraint from τ → 3μ implies λ μμ = (0.013 ± 0.011)(0.15/R 0 ) g q g , while the constraint on the | C| = 2 operators from CP violation in D 0 -D 0 mixing results in the inequality g /g q > 1.26R 0 . From these ingredients one can estimate the maximum value of C 9 , which is unfortunately too small to have a detectable effect in D → πμ + μ − or in D 0 → μ + μ − .

Leptoquarks
There exist several scalar and vector leptoquark (LQ) states which may leave imprint on c → u + − transitions [37]. Among all scalar LQs we will consider only the baryon number conserving state (3, 2, 7/6), which comes with a rich set of couplings that are in general severely constrained by B and K physics [38]. Then, among the two baryon number conserving vector LQs we will focus on the state in the representation (3, 1, 5/3) whose phenomenology is limited to the up-type quarks and charged leptons. (3, 2, 7/6) The renormalizable LQ couplings for the state (3, 2, 7/6) are [38]

Scalar leptoquark
The LQ Yukawa matrices Y L and Y R are written in the mass basis of up-type quarks and charged leptons with the CKM and PMNS rotations present in the down-type quarks and neutrinos. Thus, the couplings of LQ component with charge 5/3 are The tree-level amplitude induced by a nonchiral LQ state (5/3) involves both chiralities of fermions and is matched onto the set of (axial)vector, (pseudo)scalar, and (pseudo)tensor operators: In the minimal numerical scenario, strict bounds in the downtype quark sector can be evaded completely by putting to zero the couplings to the left-handed quarks. In this case we are allowed to have significant contributions to rare charm decays via the C 9 = −C 10 contributions for which the bound from the last line of Table 2 applies. The contribution to the D 0 -D 0 mixing amplitude is matched onto the effective Hamiltonian We have assumed that the leptoquark does not couple to electrons or tau leptons. Hadronic matrix element of the above operator in mixing is customarily expressed as where the bag parameter in the MS scheme B D (3 GeV) = 0.757 (27) (4) has been computed on the lattice by the ETM Collaboration with 2 + 1 + 1 dynamical fermions [39]. The SM part of the mixing amplitude is poorly known due to its nonperturbative nature and the only robust bound on the LQ couplings is obtained by requirement that the mixing frequency (in the absence of CP violation) has to be smaller than the world average x = 2|M 12 |/ = (0.49 +0.14 −0.15 )% as quoted by the HFAG [24], where r = 0.76 is a renormalization factor due to running of C 6 from scale m = 1 TeV down to 3 GeV [3]. Finally, we find a bound on C 9 slightly stronger than but comparable to the one obtained from D 0 → μ + μ − : |C 6 (m )| < 2.5 × 10 −13 GeV −2 ⇒ |C 9 ,C 10 | < 0.34.
One can imagine an extension of this scenario which would include also scalar and tensor operators. Namely, we consider a numerically tuned example with m = 1 TeV and large Y R cμ = 3. The bound on C 10 from D 0 → μ + μ − would then impose the smallness of the coupling Y R uμ , Y R uμ < 0.007. Bounds of similar strengths are expected from D 0 -D 0 mixing. Now one can introduce a nonzero coupling to the lefthanded quark doublet Y L μu that would, together with large Y R cμ contribute to the Wilson coefficients C S,P and C T,T 5 . However, a very strong bound on C S now emerges from D 0 → μ + μ − and limits the left-handed coupling, Y L μu < 1.2 × 10 −3 . Thus we can realize −C 10 =C 9 = 0.63, together with small enough Y L μu = 1.2×10 −3 to comply with the constraints from B, K physics and four fermion operator constraints [40]. leptoquark (3, 1, 5/3) The interactions of the vector LQ state V (5/3) (3, 1, 5/3) with the SM fermions are contained in a single term at the renormalizable level:

Vector
Generation indices are denoted i, j. Integrating out V (5/3) results in the right-handed current operators: On the other hand, the same combination of couplings enters the D 0 -D 0 mixing. We employ the same type of Hamiltonian as in the preceding section, this time with the Wilson coefficient: A consequence of the bound (27) is that the rare decay Wilson coefficients are limited: The above knowledge of C 9,10 implies that the branching ratio of D → πμ + μ − in the high-q 2 bin is at most 1.4 × 10 −8 , where the long-distance uncertainties have been taken into account. The effect is twice smaller than the existing experimental bound.

Two Higgs doublet model type III
In the Two Higgs doublet model of type III (THDM III) the neutral Higgses have flavor changing couplings to the fermions. The spectrum includes two neutral scalars, h and H , one pseudoscalar, A, and two charged scalars, H ± . In the scenario with MSSM-like scalar potential their masses and mixing angles are related [41], where β, tan β = v u /v d , is the angle that diagonalizes the mass matrix of the charged states, α is the mixing angle of neutral scalars. The vacuum expectation values are normalized to the electroweak vacuum expectation value, v/ GeV. The part of the interaction Lagrangian responsible for FCNCs in the up-type quarks and charged leptons is [41] and the neutral Yukawa couplings for the charged leptons and up-type quarks are respectively. The flavor off-diagonal terms f i , u f i are free parameters of the model. The coefficients x k q for H k = (H, h, A) are determined by the mixing angles of the neutral scalars and the VEVs [41] x k u = (− sin α, − cos α, i cos β), For the transitions c → u + − the driving flavor changing parameter is u 12 , which induces scalar and pseudoscalar Wilson coefficients, while we assume that 22 is negligible [41]: The best upper bounds on C P , C S , or C P , C S pairs are obtained from BR(D 0 → μ + μ − ) and read |C S −C S | ≤ 0.05 and |C P −C P | ≤ 0.05, which makes them very difficult to probe in D → πμ + μ − decay, unless the cancellation between C S (C P ) and C S (C P ) in D 0 → μ + μ − is arranged by fine-tuning.

Flavor specific Z extension
An additional neutral gauge boson appears in many extensions of the SM. Current searches for Z at the LHC are well motivated by many extensions of the SM; see e.g. [42,43]. Even more, a Z boson can explain the B → K * μ + μ − angular asymmetries puzzle, as presented in e.g. [44,45]. Assuming as in [43] flavor nonuniversal couplings of Z to fermions, we allow Z to couple only to the paircu and cū. Such a model in the most general way has been considered by the authors of [3]. In order to avoid constraints coming from the down-type quark sector which will affect left-handed quark couplings, we allow only right-handed couplings of L q Z = C u (ūγ μ P R c)Z μ . This assumption leads to the same effective operator H eff = C 6 (ūγ μ P R c)(ūγ μ P R c) as already discussed in the case of leptoquarks. The effective Wilson coefficient describing the D 0 -D 0 transition is now The bound on C 6 (27) leads to |C u | < 7.1 × 10 −4 (m Z / 1 TeV). Allowing Z to couple to muons as in the SM with g L = (g/ cos θ W )(−1/2 + sin 2 θ W ) and g R = g sin 2 θ W / cos θ W , we obtain and For m Z ∼ 1 TeV this amounts to |C 9 | 8 and |C 10 | 100 (|C 9 | < 10 −3 , and |C 10 | < 0.014); this induces negligible effects in D → πμ + μ − and D → μ + μ − decays.

Lepton flavor universality violation
Lepton flavor universality was checked in the case of B → K + − with = e, μ by the LHCb experiment [15] in the low dilepton invariant mass region, q 2 ∈ [1, 6] GeV 2 . The disagreement between the measurement and the value predicted within the SM is 2.6 σ [46]. This disagreement might be a result of NP, as first pointed out in Ref. [46]. Many subsequent studies found a number of models which can account for the observed discrepancy. In the following we assume that the amplitude for D + → π + e + e − receives SM contributions only, while in the case of the π + μ + μ − mode, there can be NP contributions, similarly to what was assumed for R K in Ref. [47]. We define LFU ratios in the low-and high-q 2 regions by and In the SM the departure of the above ratios from 1 comes entirely from the lepton mass differences. We find R I,SM π = 0.87 ± 0.09 in the low-q 2 and R II,SM π = 0.999 ± 0.001 in the high-q 2 region, where in the latter region both leptons are effectively massless. In Table 3 we quote ranges for the ratio R II π for the maximal allowed values of Wilson coefficients by rare charm decays considered in the previous sections. Generally we find that with currently allowed Wilson coefficients and assuming no NP contribution in the electronic modes these ratios could become much larger. The spread in these predictions is large because of the unknown relative phases in the resonant part of the spectrum, i.e., BR(D + → Note that large enhancements are allowed in the scenarios which are currently constrained by D + → π + μ + μ − . In the low-q 2 region the interference terms in R I π are even more pronounced since the effect of nearby ρ resonance is interfering either in positive or in negative direction, and thus we cannot conclude the sign of the deviation from the SM value of R I π .

Summary and outlook
Motivated by the great improvement of the bounds on rare charm decays by the LHCb experiment we determine the bounds on the effective Wilson coefficients. Existing data implies upper bounds on the effective Wilson coefficients as presented in Table 2. The strongest constraints on C 10 , C P , C S , and C 10 , C P , C S are obtained from the bound on the branching fraction of D 0 → μ + μ − decay. The nonresonant differential decay width distribution gives bounds on C i , i = 7, 9, 10, S, P, T, T 5 as well as on the coefficients of the operators of opposite chirality. The constraints are stricter in the high dilepton invariant mass bin than in the low dilepton invariant mass bin, and this statement applies in particular to the contributions of the scalar and pseudoscalar operators. The forward-backward asymmetry is sensitive to the combination of scalar and tensor coefficients at high-q 2 . Next, we have investigated new physics models in which the effective operators may be generated. We have found that the presence of a leptoquark, which is either a scalar and weak doublet, (3, 2, 7/6), or has spin-1 and is a weak singlet, (3, 1, 5/3), can lead to sizable contributions to the Wilson coefficients C 9 and C 10 . The sensitivity to the LQ scenarios is similar in the high-q 2 bin of D + → π + μ + μ − and D 0 → μ + μ − , while D 0 -D 0 mixing results in a somewhat stronger constraint. For the Two Higgs doublet model of type III the presence of scalar and pseudoscalar operators enhances the sensitivity in D 0 → μ + μ − and therefore results in small effects in D + → π + μ + μ − . We have also discussed a SM extension by a Z gauge boson where the tree-level amplitude in D 0 -D 0 mixing is a dominant constraint and leaves no possibility of signals in rare charm decays.
Our study indicates the possibility to check whether lepton flavor universality between muonic and electronic channels is valid by means of studying ratios of widths of D + → π + + − at low or high dilepton invariant mass bins, R I,II π . In the SM the two ratios are close to 1, especially in the highq 2 bin. Assuming the electronic decay is purely SM-like, we find that in the high-q 2 bin the ratio R II π is in most cases significantly increased with respect to the SM prediction, while there is no clear preference between higher and lower values at low-q 2 bin ratio R I π . In the leptoquark models studied in this paper the ratio may be greatly increased, but a slight decrease cannot be excluded, presently due to the unknown interplay of weak phases with the phases of the resonant spectrum. Chances to observe new physics in rare charm decays are present in models where the connection to the stringent constraints stemming from B and K flavor physics are hindered. New physics models which fulfill this condition are main candidates to be exposed experimentally by future progress in bounding the rare charm decays D → πμ + μ − and D 0 → μ + μ − , as well as by more precise studies of D 0 -D 0 mixing with the potential NP contributions. Alternatively, experimental tests of lepton flavor universality in rare charm decays might point toward the presence of new physics in the charm sector, which can easily be hidden in the case of existing experimental observables.
Note While we were finishing this paper another work [48] appeared in which the authors studied rare charm decays.