Prospects of discovering new physics in rare charm decays

The LHCb bounds on the branching ratio of rare decay $D^0 \to \mu^+ \mu^-$ and the constraints on the branching ratio of $D^+ \to \pi^+ \mu^+ \mu^-$ in the nonresonant regions enable us to improve constraints on new physics contributions. Using the effective Lagrangian approach we determine 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'$ model. Here we complement the discussion by $D^0 - \bar D^0$ oscillations data. Among considered models, only leptoquarks can significantly modify Wilson coefficients. Assuming that the differential decay width for $D^+ \to \pi^+ \mu^+ \mu^-$ receives NP contribution, while the differential decay width for $D^+ \to \pi^+ e^+ e^-$ is Standard Model-like, we find that lepton flavor universality can be violated and might be observed at high dilepton invariant mass.


I. 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 short distance 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 very 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 |∆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 then down-type quark sector whenever new bosons couple to left-handed 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 to avoid the strong flavor constraints in the down-type quark sector.
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 large departure of the experimentally determined lepton flavor universality (LFU) ratio R K = BR(B → Kµ + µ − ) q 2 ∈[1,6]GeV 2 /BR(B → Ke + e − ) q 2 ∈[1,6]GeV 2 from the expected SM value [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 Section 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 → µ + µ − ). Sec. 3 contains analysis in the context of specific theoretical models of new physics, contributing to the c → u + − and related processes. Sec. 4 discusses lepton flavor universality violation. Finally, we summarize the results and present conclusions in Sec. 5.

II. OBSERVABLES AND MODEL INDEPENDENT CONSTRAINTS
A. 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 is 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 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: 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] f where x = q 2 /m 2 D * , f T (0) = 0.46(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]: Vū / pv + Aū / pγ 5 v + (S + T cos θ)ūv + (P + T 5 cos θ)ūγ 5 v .
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 introduced a shorthand notation λ = λ(m 2 D , m 2 π , q 2 ), where λ(x, y, z) = (x+y+z) 2 −4(xy+yz+zx), as well as β = β (q 2 ) = 1 − 4m 2 /q 2 . 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 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 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 the 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 additional long distance contribution. However, it was shown in [29] and [30] that contributions of these states is 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 Tab. I [31]. The corresponding values of |a X | are given in the second row in Tab. I. 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 (Tab. I) and random phases of a X . On the right-hand 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 total resonant decay branching ratio and find for the low-and high-q 2 bin contributions to be smaller than 7.3 × 10 −9 and 5.3 × 10 −9 , respectively. On the other hand, the short distance 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.    Figure 1. SM resonant contributions in D + → π + µ + µ − shown in orange. On the right-hand side panel cyan regions correspond to a scenario with constant decay amplitude that would saturate LHCb bounds (3).

III. 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 we show the kinematical effect of setting to 1 individual Wilson coefficients one at the time, where we have neglected the SM resonant contributions. 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 Breit-Wigner parameters a η,ρ,ω,φ , where the latter are sampled in their 1σ regions (90% CL bound for |a ω |), c.f. Tab. I. The most relaxed bound obtained in this way is then reported in Tab. II, where we use notationC i = V ub V * cb C i . At the same time the branching ratio of D 0 → µ + µ − can give bound on Wilson coefficients C 10 , C S , and C P . It turns out that the upper bound on BR(D 0 → µ + µ − ) is more restrictive for 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 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 towards 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 in this case is strongly enhanced in the high-q 2 region, as shown in Fig. 3.   Table II. Maximal allowed values of the Wilson coefficient moduli, |Ci| = |V ub V * cb Ci|, 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 whereC9 = ±C10. All the quoted bounds have been derived for real Ci. The bounds for Ci apply also to the chirally flipped coefficientsC j .

IV. IMPACT ON SPECIFIC MODELS
A. 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, ij 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 * . 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. Constraint from τ → 3µ implies λ µµ = (0.013±0.011)(0.15/R 0 ) gq g , while the constraint on the |∆C| = 2 operators from CP violation in D 0 −D 0 mixing results in an inequality g /g q > 1.26R 0 . From these ingredients one can estimate the maximum value of C 9 , C 9 50 0.013 ± 0.011 R 0 10 , which is unfortunately too small to have a detectable effect in D → πµ + µ − or in D 0 → µ + µ − .

B. 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. The renormalizable LQ couplings for the state ∆(3, 2, 7/6) are [38] 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 down-type 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 Tab. II applies. The contribution to D 0 −D 0 mixing amplitude is matched onto the effective Hamiltonian We have assumed that leptoquark does not couple to electrons or tau leptons. Hadronic matrix element of the above operator in mixing is customarily expressed 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 but comparable to the one obtained from D 0 → µ + µ − : 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 coupling Y R uµ , Y R uµ < 0.007. Bounds of similar strength are expected from D 0 −D 0 mixing. Now one can introduce a nonzero coupling to left-handed 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]. 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: Generation indices are denoted by 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 the Wilson coefficient: 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.

C. 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/ √ 2 = v 2 u + v 2 d = 174 GeV. The part of the interaction Lagrangian responsible for FCNCs in the up-type quarks and charged leptons is [41] L = y 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 β) , x k d = (− cos α, sin α, i sin β) .
For the transitions c → u + − the driving flavor changing parameter is u 12 that 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.

D. 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 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 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 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), and induces negligible effects in D → πµ + µ − and D → µ + µ − decays.

V. 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 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 π + µ + µ − mode, there can be NP contributions, similarly to what was assumed for R K in Ref. [47]. We define LFU ratios in the lowand high-q 2 regions as and 0.74 ∼ 6-60 C9 = ±C10 0.63 ∼ 3-60 C 9 = −C 10 LQ(3,2,7/6) 0.34 ∼ 1-20 Table III. The LFU ratio R II π at high dilepton invariant mass bin and maximal value of each Wilson coefficient (applies also for the primed coefficients,C i ). It is assumed that NP contributes only to the muonic mode. The SM value of R II π is given in the first row.
In the SM the departure of the above ratios from 1 comes entirely from 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 Tab. III 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 electronic modes these ratios could become much larger. The spread in these predictions is large because of unknown relative phases in the resonant part of the spectrum, i.e., BR(D + → π + e + e − ) ≈ BR(D + → π + µ + µ − ) ≈ (0.5-5.3) × 10 −9 . 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 deviation from the SM value of R I π .

VI. SUMMARY AND OUTLOOK
Motivated by the great improvement of bounds on rare charm decays by the LHCb experiment we determine bounds on the effective Wilson coefficients. Existing data implies upper bounds on the effective Wilson coefficients as presented in Tab.II. 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. Forward-backward asymmetry is sensitive to the combination of scalar and tensor coefficients at high-q 2 .
Then, 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 scalar and weak doublet, (3, 2, 7/6), or has spin-1 and is a weak singlet, (3, 1, 5/3), can lead to sizeable contributions to the Wilson coefficients C 9 and C 10 . Sensitivity to the LQ scenarios is similar in high-q 2 bin of D + → π + µ + µ − and D 0 → µ + µ − , while D 0 −D 0 mixing results in somewhat stronger constraint. For the Two Higgs doublet model of type III the presence of scalar and pseudoscalar operators enhances 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 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 a 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 high-q 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 slight decrease cannot be excluded, presently due to unknown interplay of weak phases with the phases of resonant spectrum. Chances to observe new physics in rare charm decays are possible 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 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 towards 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.