Impact of $b \to s \ell \ell$ anomalies on rare charm decays in non-universal $Z'$ models

In this work, we study the impact of $b \to s \ell \ell$, $B_s - \bar{B_s}$ mixing and neutrino trident measurements on observables in decays induced by $c \to u $ transition in the context of a non-universal $Z'$ model which generates $C^{\rm NP}_{9}<0$ and $C^{\rm NP}_9 = - \,C^{\rm NP}_{10} $ new physics scenarios at the tree level. We inspect the effects on $D^0 \to \pi^0 \nu \bar{\nu}$, $D^+ \to \pi^+ \nu \bar{\nu}$ and $B_c \to B^+ \nu \bar{\nu} $ decays which are induced by the quark level transition $c \to u \nu \bar{\nu}$. The fact that the branching ratios of these decays are negligible in the standard model (SM) and the long distance effects are relatively smaller in comparison to their charged dileptons counterparts, they are considered to provide genuine null-tests of SM. Therefore the observation of these modes at the level of current as well as planned experimental sensitivities would imply unambiguous signature of new physics. Using the constraints on $Z'$ couplings coming from a combined fit to $b \to s \ell \ell$, $\Delta M_s$ and neutrino trident data, we find that any meaningful enhancement over the SM value is ruled out in the considered framework. The same is true for $D - \bar{D}$ mixing observable $\Delta M_D$ along with $D^0 \to \mu^+ \mu^-$ and $D^+ \to \pi^+ \mu^+ \mu^-$ decay modes which are induced through $c \to u \mu^+ \mu^-$ transition.

= − C NP 10 new physics scenarios at the tree level. We inspect the effects on D 0 → π 0 νν, D + → π + νν and B c → B + νν decays which are induced by the quark level transition c → uνν. The fact that the branching ratios of these decays are negligible in the standard model (SM) and the long distance effects are relatively smaller in comparison to their charged dileptons counterparts, they are considered to provide genuine null-tests of SM. Therefore the observation of these modes at the level of current as well as planned experimental sensitivities would imply unambiguous signature of new physics. Using the constraints on Z couplings coming from a combined fit to b → s , ∆M s and neutrino trident data, we find that any meaningful enhancement over the SM value is ruled out in the considered framework. The same is true for D −D mixing observable ∆M D along with D 0 → µ + µ − and D + → π + µ + µ − decay modes which are induced through c → uµ + µ − transition.

I. INTRODUCTION
The currently running experiments at the LHC have provided several engrossing indirect hints of physics beyond the standard model (SM) of electroweak interactions. Many such signatures of new physics are in observables related to the decays induced by the quark level transition b → s + − (l = e, µ) [1]. The most interesting discrepancies are related to lepton flavour universality (LFU) violation in decays B → (K, K * ) l + l − as within the SM, the coupling of leptons of gauge bosons are flavour independent. The measurements of LFU testing ratios R K ≡ Γ(B + → K + µ + µ − )/Γ(B + → K + e + e − ) [2] and R K * ≡ Γ(B 0 → K * 0 µ + µ − )/Γ(B 0 → K * 0 e + e − ) [3] disagree with the SM prediction which is very close to unity [4][5][6], with small corrections due to the muonelectron mass difference, at the level of 3.1 and 2.5σ, respectively. Very recently such ratios were measured in B 0 → K 0 S l + l − and B + → K * + + − decay modes by the LHCb collaboration using data set corresponding to an integrated luminosity of 9 fb −1 [7]. These measurements have relatively large uncertainties and they are consistent with the SM value at 1.5 and 1.4σ, respectively.
These discrepancies can be attributed to new physics in muon or/and electron sector.
There is another class of discrepancies that can be imputed to new physics only in the muon sector. These are in angular observables constructed from kinematical distribution of differential decay widths of B → K * µ + µ − and B s → φ µ + µ − decays. A series of measurements from different experiments at the LHC along with the Belle experiment have reported deviations from their SM predictions of the angular observable P 5 in B → K * µ + µ − decay at the level of 3σ [8][9][10]. In addition, the measured value of the branching ratio of B s → φ µ + µ − decay shows tension with the SM prediction at the level of 3.5σ [11] .
The piling up of these disparateness with the SM can be considered as signals of new physics. several groups, see for e.g. [12][13][14][15][16][17][18][19][20]. These analyses differ mainly in the treatment of hadronic uncertainties and the statistical approach. Irrespective of the adopted methodology in these analyses, new physics scenarios with non-zero muonic WCs corresponding to vector/axial-vector operators are statistically favoured.
The global analyses of b → s + − data suggest various new physics solutions. Some of these so-lutions are preferred over the SM with a very high significance. For e.g., considering one new physics operator or two related operators at a time and assuming new physics only in the muon sector, the O 9 ≡ (sγ µ P L b)(µγ µ µ) operator as well as a combination of O 9 and O 10 ≡ (sγ µ P L b)(µγ µ γ 5 µ) with C N P 9 = − C N P 10 can account for all b → s + − data. These model independent solutions can be realized in several new physics models. One such simplified model is Z with non-universal couplings which generates b → s µ + µ − decay mode at the tree level, see for e.g. [21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36]. Various versions of these models are widely studied in the context of B-anomalies. However, the Z boson can also contribute to other sectors. The muon coupling to Z is constrained by b → s and neutrino trident data and hence it is possible to look for implications of such couplings in other sector as well. Such a correlation in b → d sector was studied in [37].
In this work we study implications of current b → s and neutrino trident data on several observables in the charm sector in the context of Z models. In particular we study decays induced by the quark level transition c → uνν and c → uµ + µ − . These rare decays provide unique opportunity for probing new physics in the up quark sector. The SM symmetries are utilized to define null-test observables with small theoretical uncertainties in these semi-leptonic decays which are usually dominated by vector resonances [38]. Therefore it would be interesting to investigate the impact of possible new physics in b → s sector on charm decays which can play complementary role in understanding the nature and origin of new physics apparent in the down sector. The theoretical potential embedded in the charm sector can be complemented and utilized by the several ongoing and planned experimental facilities such as LHCb [39], Belle II [40], BES III [41], the FCC-ee [42] and Super Charm-Tau factory [43].
The paper is organized as follows. In section II, we introduce the Z model considered in this work along with its contributions to the effective Hamiltonian of b → s and c → u processes. In the next section we discuss observables which provide useful constraints on Z couplings relevant for the b → s sector. We also provide the fit results. Using the constraints obtained on new physics couplings, in section IV, we provide predictions for several observables in the charm sector. We also show correlations between potential observables in charm and b → sµ + µ − sectors. Finally, the conclusions of this work are illustrated in section V.

II. THE NON-UNIVERSAL Z MODEL
We consider a non-universal Z model where the Z boson is associated with an additional U (1) symmetry. In this model, the b → s µ + µ − transition is generated at the tree level. The Z boson couples to both left and right-handed muons but not to leptons of other generations. Further, the couplings to both left and right-handed quarks are allowed. However, in order to avoid contribution of new chirality flipped operators to flavour changing neutral current (FCNC) decays, the couplings to right-handed quarks are assumed to be flavour-diagonal [44,45]. As we are interested in b → s processes, the change in the Lagrangian density due to the addition of this heavy Z boson can be written as where where g µµ L(R) are the left-handed (right-handed) couplings of the Z boson to muons, and g bs L to quarks. The right handed quark coupling, g bs R , cannot contribute to the b → s processes as these couplings are assumed to be flavour diagonal. Further, P L(R) = (1 ∓ γ 5 )/2, Q i is the i th generation of quark doublet, and L = (ν µ , µ) T is the second generation doublet.
After integrating out the heavy Z , we get the effective four-fermion Hamiltonian which apart from contributing to b → s µ + µ − transition, also induces b → s νν and c → u transitions. The relevant terms in the effective Hamiltonian is given by where h cu L = g bs L V us V * cb . The fact that the down-type quarks in the quark-doublets Q i are taken to be in the mass-flavour diagonal basis, u i → u j transitions are induced by the up-type quarks in the quark doublets.
The first term in eq. (3) induces b → sµ + µ − transition whereas the second term gives rise to B s -B s mixing. The third term contributes to the neutrino trident production ν µ N → ν µ N µ + µ − (N = nucleus). The fourth term generates b → s νν decay whereas the remaining terms induces FCNC transitions in the up quark sector. The contribution to c → uµ + µ − and D 0 -D 0 in this model can be calculated through fifth and sixth terms, respectively. The seventh term corresponds to charm dineutrino decay. The product g bs L g µµ L,R and the individual magnitude |g bs L | are constrained by the b → sµ + µ − data and B s -B s mixing, respectively. The constraints on individual muon couplings g µµ L,R are obtained through the measurement of the neutrino trident production. Observables in c → u sectors depend upon specific combinations of these couplings which can be determined from eq. (3). Thus it would be interesting to see the implications of b → sµ + µ − , B s -B s and neutrino trident data on possible new physics effects in charm decays. We now discuss contributions of Z boson to these processes.

A. b → s transitions
In this subsection, we calculate the contributions of Z boson to b → s transitions. Here we consider b → sµ + µ − and b → sνν decays along with B s −B s mixing.
The effective Hamiltonian for b → sµ + µ − transition in the SM is given by Here G F is the Fermi constant and V ij are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The short distance contributions are encoded in the Wilson coefficients (WC) C i of the

The contributions of operators
Here q 2 is the invariant mass-squared of the final state muon pair. The Z boson contributes to b → sµ + µ − decay mode through the first term in eq. (3). This contribution modifies the SM WCs C SM 9,10 as C 9,10 → C SM 9,10 + C NP 9,10 . The new contributions C NP 9,10 are given as It is well known that the new physics scenario C NP 9 < 0 as well as C NP 9 = − C NP 10 provide a good fit to all b → sµ + µ − data. It is evident from eq. (5) that these one-dimensional (1D) new physics solutions can be generated by substituting g µµ L = g µµ R and g µµ R = 0, respectively. In this work, we consider both scenarios. We denote g µµ L = g µµ R and g µµ R = 0 scenarios as Z 1 and Z 2 , respectively.
is the Inami-Lim function [46] and η B ≈ 0.84 encodes perturbative QCD corrections at two loop [47]. The Z boson generates B s −B s mixing at the tree-level through the second term in eq. (3) which is the same operator as in the SM. The combined contribution is written as [48] ∆M SM+Z for M Z = 1 TeV.

b → sνν decay
Like b → sl + l − transition, b → sνν occurs via one-loop electroweak penguin or box diagrams in the SM. However b → sνν decay rates are predicted with relatively smaller theoretical uncertainties in comparison to the corresponding b → sl + l − decay modes. This is mainly due to the absence of the long-distance electromagnetic interactions. The effective Hamiltonian for b → sνν transition is given by [49] Here B. Neutrino trident production Due to SU (2) L invariance, the Z also couples to the left-handed neutrinos. This leads to the neutrino trident production, ν µ N → ν µ N µ + µ − . The the third term in eq. (3) modifies the cross section σ for neutrino trident production as [50] where v = 246 GeV and s W = sin θ W .
The Z contributes to c → u transition through fifth, sixth and seventh terms in eq. (3). These terms induce c → uµ + µ − decay, D-D mixing and c → uνν transition, respectively.
The c → uµ + µ − process in the SM can be described by the following effective Hamiltonian at where the three contributions correspond to diagram with intermediate quarks d, s, b. It is customary to include the contributions of the states heavier than the charm quark in H peng as The following operators appearing in the above effective Hamiltonian are in particular sensitive to new physics effects The only non-vanishing WCs in the SM are C cu 7,9 . Their values are determined mainly through the effects of QCD renormalization. The value of C cu 7 (m c ) determined by two loop mixing with current-current operators is [51,59]. After inclusion of renormalization group running effects, the value of C cu 9 (m c ) turns out to be small [53]. On the other hand, as the renormalization group running does not affect O 10 , the WC C cu 10 (m c ) is negligibly small [60]. As C cu 10 ≈ 0, the effects of V − A structure of the SM are switched off at the charm scale. This feature makes FCNC charm decays distinct from that of the corresponding B and K decays.

D-D mixing
In the SM, D 0 -D 0 mixing is induced by the box diagrams with d, s and b quarks in the loop.
Due to a strong GIM cancellation, the short-distance contribution is extremely small, ∆M D = O(10 −4 ) ps −1 . Therefore the short-distance contribution cannot explain the measured value of [61]. In particular, the contribution due to b-quark is highly suppressed, O(λ 8 ). The fact that D 0 -D 0 mixing is dominated by the d-and s-quarks, there can be large long-distance contributions, for which there are no reliable estimates at present [62,63]. As long-distance contributions are unknown, in our analysis we focus on the short-distance contributions to the D 0 -D 0 mixing parameter ∆M D .
In Z model, D 0 -D 0 mixing is induced at the tree level and hence expected to provide a much larger contribution in comparison to the short-distance SM contribution. ∆M D generated by the sixth term in eq. (3) is given by

c → uνν decay
Within the SM, the c → uνν transition is induced by box and Z-penguin diagrams. The corresponding short-distance effective Hamiltonian is given by [63] H where Here λ q = V uq V * cq and the structure functions are defined by X (x q ) =D(x q , y )/2, whereD(x q , y ) is given as [46]D (x q , y ) = 1 8 with x q = m 2 q /M 2 W and y l = m 2 l /M 2 W . The Z contributes to c → uνν transition through the seventh term in eq. (3). This contribution modifies the WC as C → C SM + C νν µ , where

III. CONSTRAINTS ON NP COUPLINGS
In the literature there are several works where the rare charm decays have been studied in the context of Z models. However, in most of these works, the constraints coming from b → s µ + µ − sector were not considered. To emphasize the fact that b → s µ + µ − sector can provide valuable constraints on the allowed new physics parameter space and hence restrict large new physics contributions to several observables in the up sector, we perform four types of fit: • F1: only b → sνν data, • F2: b → sνν data and B s −B s mixing, • F3: b → sνν data, B s −B s mixing and neutrino trident,  In the following, we discuss all observables used in the fit. These observables can be classified into following four sectors: These upper bounds are a factor of two to five above the SM predictions. With 50 ab −1 of data, Belle-II experiment is expected to observe these decay modes to an accuracy of 10% on the branching ratio [40,72]. In order to include these observables in the fit, we convert the above upper limits to a branching ratio as (0.5 UB ± 0.625 UB) where UB is the experimental upper bound at 90% C.L. The χ 2 is then written as where the theoretical predictions for B → K ( * ) νν branching ratios are obtained using Flavio [73]. We do not use constraints from the branching ratio of B s → φνν as it is much weaker than B → K ( * ) νν decay modes.
B s −B s sector: Here we consider constraints coming from the mass difference ∆M s . Using eq. (7), the constraints on g bs L can be translated in the following form (g bs L ) 2 = (7.69 ± 12.94) × 10 −6 .
Here we used the value of ∆M SM s /∆M Exp s given in Table I Neutrino trident: The contribution to the total χ 2 coming from neutrino trident production is The theoretical expression of R ν is given in eq. (10) whereas the experimental value is taken from Table I.
b → s µ + µ − sector: A fit to all b → s µ + µ − data leads to constraints on C NP 9 and C NP 9 = − C NP 10 scenarios as given in Table I. The fit value for C NP 9 (C NP 9 = − C NP 10 ) provides constraints on Z 1 (Z 2 ) model. Using the fit values and other inputs from Table I along with eq. (5), we get We use above constraints on (g bs L g µµ L ) in the fit. Hence the χ 2 function can be written as Therefore, the χ 2 for fits F1-F4 can be written as The χ 2 fit is performed using the CERN minimization code MINUIT [74]. The mass of Z is assumed to be 1 TeV in the fits.
FIG. 1: The left and right panels depict 1σ-favoured parameter space of (g bs L , g µµ L ) couplings for Z 1 and Z 2 models, respectively. The black, blue and red regions show allowed parameter space corresponding to F2, F3 and F4 fits, respectively. The allowed regions correspond to M Z = 1 TeV.
For F1 fit, χ 2 is a function of the product (g bs L g µµ L ). The fit result is (g bs L g µµ L ) = −0.0047 ± 0.0051.
The results for F2, F3 and F4 fits for Z 1 and Z 2 models are presented in Fig. 1. It is obvious from the figure that if we only include constraints coming from b → sνν and ∆M s , a relatively large parameter space is allowed. Therefore while correlating c → u decays, in particular c → uνν, with b → s sector if only these constraints are considered, a possibility of large new physics effects may survive.
However, the parameter space shrinks after including neutrino trident in the fit. Finally, after including b → sµ + µ − constraints, the allowed range of new physics couplings reduces considerably.
Therefore it will be interesting to see whether any useful new physics effects in charm sector, particularly decays induced by c → uνν and c → uµ + µ − transitions, is allowed by the current b → s and neutrino trident data.

IV. SEMI-LEPTONIC CHARM DECAYS AND MIXING
Rare or forbidden decays of charm mesons provide a unique channel to probe new physics in the up-type quark sector [38]. This is complementary to new physics searches in the first and third generations of quarks as well as in the down-quark sectors. The fact that the charm quark is not as heavy as bottom, the application of heavy quark effective theory is not as useful as it is for the B-meson sector. Also the mass of charm quark is not small enough to be considered as a light quark. Therefore predictions in the charm sector are usually dominated by the long-distance effects. Within the SM, the short-distance contributions to the purely leptonic and semi-leptonic charm decays along with D −D mixing are highly suppressed, in fact almost negligible, owing to GIM cancellation and CKM suppression. Probing new physics then becomes a challenging task as an unambiguous signal would require enhancements above the level of long-distance contributions.
In the following subsections we discuss the impact of measurements in B meson sector on rare FCNC processes induced by the c → u transition in the context of non-universal Z models. We consider decays induced by c → uνν and c → uµ + µ − transitions along with D −D mixing. The primary motive is to identify observable(s) for which the current b → s data allows short-distance contributions above the long-distance contributions.

A. Branching ratio of dineutrino charm decays
The new physics searches can be exquisitely performed using the null-tests of the SM. In the up-sector, c → u transitions provide such null tests due to strong CKM suppression. In this regard, the decays induced by the quark level transition c → uνν are particularly interesting as the longdistance effects are relatively smaller in comparison to their charged dileptons counterparts. Hence in c → u sector, dineutrino modes can be considered to provide genuine null-tests of SM [38,[75][76][77][78].
The c → uνν transition induces several exclusive decay modes such as D + → π + νν, D 0 → π 0 νν and B + c → B + νν. Any observation of these decays at the level of current experimental sensitivity can be considered as unambiguous signature of beyond SM physics.
On the experimental front, currently there are no experimental upper limits on any of the semileptonic decays included by c → uνν transition. We only have upper limits obtained using the SU (2) L invariance and bounds on the charged lepton decay modes which allow branching ratios to be as high as a few times 10 −5 [75]. These limits can go down if charged leptons bounds are improved. Owing to the clean hadronic environements, the e − e − colliders such as Belle II [40], BES III [41] as well as the future colliders, for instance FCC-ee [42] are well suited for the studies related to the dineurino decay modes. For efficiencies of a permille or better, D decay modes can be observed at the Belle II and FCC-ee experiments provided new physics enhances their branching ratios up to a level of O(10 −6 ) -O(10 −8 ) [76].
In this section we analyse D 0 → π 0 νν, D + → π + νν and B + c → B + νν decays in the context of non-universal Z models. The differential branching ratio of M 1 → M 2 νν mode, where M 1 = D 0 , D + , B c and M 2 = π 0 , π + and B + , respectively is given by where q 2 denotes the invariant mass-squared of the dineutrinos, λ ≡ λ(m 2 M 1 , m 2 M 2 , q 2 ) is the Kallen function defined as λ(x, y, z) = (x + y + z) 2 − 4(xy + yz + zx) and τ M 1 is the lifetime of the M 1 meson. The M 1 → M 2 form-factors are defined as The D → π and B c → B form-factors are calculated within the framework of lattice QCD in refs.
The integrated branching ratio is then given by For D 0 → π 0 νν mode, q 2 min = 0 whereas q 2 max = (m D 0 − m π 0 ) 2 . In D + → π + νν decay, resonant backgrounds through τ -leptons lead to the same final state. This is removed by phase space cuts due to which q 2 min for D + → π + νν mode is (m 2 Using the form-factor calculations in [79,80], we get The short-distance SM predictions for the branching ratios of D 0 → π 0 νν and D + → π + νν are ∼ 2 × 10 −16 and ∼ 8 × 10 −16 , respectively. The long-distance contributions in these decay modes are predicted to be ∼ 10 −16 [52]. For B c → B + νν mode, the short-distance SM contribution is ∼ 1.5 × 10 −16 whereas the long-distance contribution is estimated to be ∼ 10 −16 [77]. Therefore for unambiguous signature of new physics in these decay modes, the short-distance branching ratio should be enhanced well above 10 −16 . The decay mode D + → π + νν in Z models was studied in [76]. It was shown that the branching ratio of D + → π + νν can be enhanced up to ∼ 10 −6 which falls in the same ballpark as the modelindependent upper bound. In [77], it was shown that the branching ratio of B + c → B + νν in the 331 model can be enhanced up to ∼ 10 −11 . It would be interesting to see whether such large enhancements is allowed by the current b → s and neutrino trident data in the non universal Z model considered in this work.   up to an order of magnitude above the SM value. This is also evident from Fig. 2 which illustrates correlations between B(D + → π + νν) and b → sµ + µ − observables R K , P 5 and B(B s → φµ + µ − ) in Z 1 model. The conclusions are almost the same for the Z 2 model. This enhancement is not useful in the sense that it is orders of magnitude smaller than the detection level for these modes in any of the currently running or planned collider facilities.

B. Branching ratio of dimuon charm decays
In this section we consider D 0 → µ + µ − and D + → π + µ + µ − decays which are induced by the quark level transition c → uµ + µ − . In the SM, the short-distance contributions to these decay modes are highly suppressed and the dominant contribution comes from the long-distance effects.
It would be interesting to see whether the new physics contributions coming from the addition of Z can enhance the decay rates above the long-distance contributions.
The short-distance contribution to the branching ratio of D 0 → µ + µ − decay is negligible in the SM. This is due to two fold effect: • Within SM, only C cu 10 WC contributes and its value is negligibly small.
[61], the SM branching ratio is estimated to be [38,57] From the experimental side, we only have an upper bound. The most stringent upper bound is reported by the LHCb experiment in 2013 which is 6.2 × 10 −9 at 90% C.L. using the data-set corresponding to an integrated luminosity of 0.9 fb −1 [83]. Therefore the SM value of the branching ratio including the long-distance effects is well below the current experimental upper limit. Hence one can safely neglect the SM contribution while calculating the branching ration in Z model.
Under this assumption, the branching ratio is given by where f D = 212.0 ± 0.7 MeV [64] and τ D 0 = 4.101 × 10 −13 s [61]. It should be noted that Using constraints obtained on new physics couplings for Z 2 model, we find that B(D 0 → µ + µ − ) Z can be as large as 5 × 10 −13 for the F1 fit. For fits F2, F3 and F4, the upper bounds on B(D 0 → µ + µ − ) Z are 5 × 10 −13 , 2 × 10 −14 and 5 × 10 −16 , respectively. Thus it becomes apparent that the current b → s and neutrino trident data rules out any useful enhancement in The quark level transition c → uµ + µ − generates semi-leptonic decay D + → π + µ + µ − . Within SM, the short-distance contribution to the total branching ratio is highly suppressed, O(10 −12 ) [57]. The dominant contribution comes from vector resonances ρ, ω and φ, decaying to dimuon pair. The effects of these resonances are incorporated assuming naive factorization by addition of q 2 -dependent terms to C cu 9 [57]. The contribution to the total branching ratio coming from vector resonances is O(10 −10 ) [84]. On the experimental front, at present we only have upper bounds. In 2013, LHCb reported following upper limits at 90% C.L. in the low-and high-q 2 bins using data corresponding to an integrated luminosity of 1.0 fb −1 [85] Here low-and high-q 2 bins correspond to (0.0625, 0.276) GeV 2 and (1.56, 4.00) GeV 2 , respectively.

C. Charm Mixing
The observed value of ∆M D is much higher than the short-distance SM value but it is in qualitative accord with the long distance SM predictions. However as no reliable estimate of the later is obtained so far, the possibility of large new physics contributions to D −D mixing is not ruled out. Here we consider such an effect in the context of Z model. In the current framework, the contribution to ∆M D coming from a Z boson is governed by Z bs coupling along with the CKM factor (V us V * cb ). Using the allowed range of g bs L coupling from the current b → s data, we find that (∆M D ) Z < 5 × 10 −5 ps −1 . A correlation between ∆M D and ∆M s is shown in Fig. 3.
Thus we see that a possibility of large new physics contributions to ∆M D coming from a Z boson is disfavoured by the current measurements in the b → s sector.

V. CONCLUSIONS
In this paper we consider a non-universal Z model which generates b → s µ + µ − transition at the tree level. Considering the down-type quarks in the quark-doublets to be in the mass-flavour diagonal basis, u i → u j transitions are induced by the up-type quarks mediated by Z bs and Z µµ couplings along with suitable combinations of the elements of quark mixing matrix. Thus the current b → s data is expected to have potential impact on new physics contributions to several observables in the charm sector. Here we study such an impact on decays induced by c → uνν and c → uµ + µ − transitions along with ∆M D . Performing a combined fit to b → s µ + µ − , ∆M s and neutrino trident data, we obtain following results: • The SM predictions of D 0 → π 0 νν, D + → π + νν and B + c → B + νν decays, including the long-distance contributions, are ∼ 10 −16 . We find that in Z models, the branching ratios of these decay can be enhanced by an order of magnitude above their SM values. However this enhancement is not enough to enable the observation of these decay modes in any of the currently running and planned experimental facilities.
• Within SM, the mixing observable ∆M D is dominated by long-distance contributions. We find that the Z boson contribution is restricted to ∼ 10 −5 ps −1 which is well below the experimental value.
Thus we conclude that the current data in the B sector puts severe constraints on new physics contributions in semi-leptonic charm decays and mixing. Any useful enhancement in c → uνν, c → uµ + µ − decays and ∆M D is ruled out in the context of Z model considered in this work.