New physics signature in D0(D¯0)→f\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^0 (\bar{D}^0)\rightarrow f$$\end{document} effective width asymmetries

Violation of charge conjugation-parity (CP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm CP$$\end{document}) symmetry plays a major rule in the dominance of matter in our universe. A kind of CP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm CP$$\end{document} violation results from the asymmetry of the life time measured in M0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M^0$$\end{document} and M¯0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{M}}^0$$\end{document}, here M is a heavy meson, decays to final states which is referred in the literature as AΓf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_{\Gamma }^f$$\end{document}. In this paper, we give an estimation of the upper bound on |AΓf|\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|A_{\Gamma }^f|$$\end{document} for the Cabibbo Favored D0→K-π+\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 K^- \pi ^+$$\end{document} decay process in different models. We show that in the standard model, |AΓf|≲O(10-10)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|A_{\Gamma }^f| \lesssim {\mathcal {O}} (10^{-10})$$\end{document}. Recently a bound on AΓf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_{\Gamma }^f$$\end{document} has been obtained: (AΓf)Exp.=(1.6±1)×10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(A^f_{\Gamma })^{Exp.}= (1.6 \pm 1)\times 10^{-4}$$\end{document}. This result motivates further studies on AΓf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_{\Gamma }^f$$\end{document} in beyond standard model physics. In the framework of two Higgs doublet model with generic Yukawa structure, we show that |AΓf|≲O(10-7)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|A_\Gamma ^{f}|\lesssim {\mathcal {O}} (10^{-7})$$\end{document} which is several orders of magnitude smaller than the current experimental value. Finally, in the framework of left-right symmetric models in which the mixing between the left and the right gauge bosons is allowed and the left-right symmetry is not manifest at unification scale, we find that AΓf\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_{\Gamma }^f$$\end{document} can be as large as |AΓf|≲O(10-5)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|A_{\Gamma }^f|\lesssim {\mathcal {O}} (10^{-5})$$\end{document} which is one order of magnitude smaller than the experimentally measured value by LHCb collaborators.


Introduction
Symmetries play an important rule in particle physics. In the standard model (SM), masses of quarks, charged leptons and weak gauge bosons can be attributed to the breaking of the electroweak symmetry. On the other hand, the difference between the decay rates of particle and its antiparticle can be an indication of direct violation of charge-parity (CPV) symmetry. Weak decays of hadrons can serve as a probe for a e-mail: delepine@fisica.ugto.mx b e-mail: gaberfaisel@sdu.edu.tr c e-mail: jpjdramirez@yahoo.com CPV. This remark can be understood as in SM CPV originates from the presence of complex couplings in the Cabibbo-Kobayashi-Maskawa (CKM) matrix which appears only in the quark sector in the interactions of quarks and the charged weak gauge bosons W ± [1,2].
Direct CPV has been confirmed in the weak decays of K and B mesons [3][4][5][6]. On the other hand, the remarkable experimental progress in D mesons has leaded to the observation of D 0 − D 0 meson mixing [7][8][9][10] and measurements of direct CP asymmetries in D mesons decays, with precision of O(10 −3 ) [11]. A sensitive probe of CP violation in the weak decays of D 0 meson is given by the direct CP asymmetry difference, A C P , between D 0 → π + π − and D 0 → K + K − which can be expressed as The first observation of A C P was reported in 2011 by the LHCb Collaboration [12] and later confirmed by CDF [13] and Belle [14]. Recently, the LHCb collaboration has presented new measurements at Moriond 2019 and the combined value with previous LHCb results leads to [15] A C P = (−15.6 ± 2.9) × 10 −4 (2) which is 5.3 standard deviations away from zero and thus confirm direct CPV in these particular weak decays of D 0 mesons. This progress motivates further search and further studies of CP violation in D meson decays. Indirect CP violation has been searched also in the decays of D 0 andD 0 to final states K + K − , π + π − and K − π + modes [16]. This kind of CP violation results from the asymmetry of the life time measured in D 0 andD 0 decays to the same final states or equivalently the asymmetry in effective decay widths and usually denoted by A . The latest measurements are given as [16] The first two processes are single Cabibbo suppressed (SCS) and the Standard Model (SM) contributions to these asymmetries are expected to be of order 10 −4 [17][18][19][20]. On the other hand, the last one is Cabibbo favored (CF) with very suppressed direct CP asymmetry in the framework of the SM [21]. So any observation of CPV in these CF channels will be a strong hint for New Physics. At present time, all results are compatible with no direct or indirect CPV [16,[22][23][24][25][26][27].
In the literature, the study of the interference between direct CPV and mixing has been performed through the introduction of a non-universal weak phase defined as [28][29][30]. It is important to notice that this weak phase is irrelevant for the direct CPV as direct CPV is proportional to |A f | 2 − |Ā f | 2 . To get non vanishing direct CPV it is necessary to write the amplitude, A f , as a sum of at least two amplitudes with different relative weak and strong phases.
The mass eigenstates of the neutral D mesones, denoted as |D 1,2 > with masses (total widths) m 1,2 ( 1,2 ), are linear combinations of the flavor eigenstates |D 0 > and |D 0 > and can be defined as follows: with the imposed normalization condition | p| 2 + |q| 2 = 1. We consider the decay modes with f ≡ K − π + andf ≡ K + π − . These modes are examples of D 0 andD 0 decays to final non-CP eigenstate modes. In the following we denote the decay amplitudes as, Moreover, we follow Ref. [31] and express the amplitudes as here A T f e iφ T f and A T f e i( f +φ T f ) represent SM tree-level contribution, the phases φ T f and φ T f represent weak CP violating phases while f represent strong CP conserving one all generated at tree-level. It should be noted that, upon neglecting small terms of order |( , the quantities r f and rf express the relative magnitudes of subleading contributions that can arise from new physics or from SM amplitudes with suppressed CKM factors. Moreover, the phases φ f, f and δ f, f denote the relative weak and strong CP violating phases respectively that account for the difference between the phases generated by the subleading contributions and the tree-levels ones.
Using the general formalism for D 0 − D 0 mixing, it is possible to compute the widths as a function of time. The time-dependent decay rates can be expressed as [31,32] here = ( 1 + 2 )/2 is the mean D 0 width and The experimental values for the mixing and CPV parameters in D neutral mesons are given as [33] x = 0.41 +0.14 −0.15 % ( 8 ) where the fit assuming all floating parameters is used and The expressions of the time-dependent decay rates into a final statef can be obtained via the substitutions f →f in the above expressions [31]. Due to the small values of the mixing parameters x and y and |λ f | 1 and |λ −1 f | 1, these approximations are experimentally confirmed for the decay modes under consideration [34], one can expand the expressions of the time-dependent decay rates of D 0 → f and D 0 →f and keep the terms up to first order in time.
Thus, we get Upon defining [28,31] where The phases δφ λ f and δφ λ f , to first order in r f and r f , are given as [28] Using the above definitions of λ f and λ f , we find that here φ λ f and φ λ f are the arguments of −λ f and −λ f respectively. Now, one can define the following CP observable, the asymmetry in effective decay widths A f , for the D 0 and D 0 decays to final two-body non-CP eigenstate mode f [31]: Thus, using Eq. (16), we obtain Upon substitution of the expressions of φ λ f and φ λ f we get In the models where r f and r f are so small, δφ λ f δφ λ f 0 and hence A f reduces to This expression is in agreement with the result in first line of Eq. (23) in Ref. [31]. At tree-level, the amplitudes of the decay processes under concern have no CP violating weak phases and thus φ T The quantities δφ λ f and δφ λ f are expected to be small and thus we can expand A f in the preceding equation and keep terms up to linear order in δφ λ f and δφ λ f . Thus, we obtain In the SM, the contributions to the amplitudes of the CF D 0 → K − π + and the DCS D 0 → K + π − decays originate from integrating out the W ± boson mediating the treelevel diagrams. These contributions are proportional to the Fermi coupling constant, G F , and CKM matrix elements For the scenarios in which subleading contributions to the tree-level amplitudes arise from new physics with particles heavier than m W or from SM amplitudes with suppressed CKM factors, one upper bounds on r f and r f assuming no suppressions from the weak and the strong phases in a treatment similar to the one adopted in Section IV in Ref. [31] although the treatment there is for direct CP asymmetry.
As discussed in Ref. [31], and in the presence of new contributions to the QCD penguin and dipole operators, one can use the QCD factorization as a framework to obtain order-of-magnitude estimates for the amplitudes. However, for hadronic D decays, the 1/m c expansion is not expected to work very well [31]. This can be understood as the mass of the charm quark is of order 1.5 GeV which is not heavy enough to allow for a sensible heavy quark expansion like the case of 1/m b expansion in B meson decay. Thus, as proposed in Ref. [31], see appendix A for details, one can ignore O(α s ) corrections to the matrix elements, as they are negligible compared to the overall theoretical uncertainties and work primarily at leading order in QC D /m c , using naive factorization (NF) for tree and QCD penguin operators in the effective Hamiltonian governs the decay process. Thus, in the following we adopt NF in our analysis to give an estimation of the upper bounds on r f and r f for the CF decay mode D 0 → K − π + in the framework of the SM and classes of NP models.

The effective time-integrated CP asymmetry in the SM
In the SM, the total amplitudes of D 0 → K − π + and D 0 → K + π − decay processes can be expressed as [21,40] A SM with X P 1 P 2 P 3 is given by here f P is the P meson decay constant and F P 2 P 3 0 is the form factor. In Eq. (23) a 1 = c 1 +c 2 /N C and a 2 = −(c 2 +c 1 /N C ) where N C is the color number account for the tree-level contributions to the amplitudes. These coefficients originate from integrating out the W ± boson mediating the tree-level diagrams. On the other hand, and in the same equation, a 1,2 and a 1,2 express the contributions to the amplitudes resulting from integrating out the W ± boson mediating the box and di-penguin diagrams. These loop contributions are essential for generating the weak phase required for having nonvanishing A f as the tree-level contributions are real. Their expressions are given as where κ = (m 2 D + m 2 K )/2 + 3m 2 π /4 and The quantities B x and P g originate from the box and dipenguin diagrams respectively and their expressions are given as where [41] f (x, y) = 7x y − 4 4(1 − x)(1 − y) x y 4 and the Inami function E 0 (x) is given as Turning now to a 1,2 we find that their expressions are given as where the quantities B x , P g can be obtained from the expressions of B x , P g , given in Eq. (27), via the replacement d ↔ s.

The effective time-integrated CP asymmetry in NP models
In this section we consider two particular extensions of the SM based on their potentials to enhance CP violation due to the presence of new complex couplings. The first model is based on extending the scalar sector of the SM to include new Higgs doublet. The other model is based on extending the gauge symmetry of the SM to include new gauge group.
In both models, the new interactions can provide new sources for the weak CP violating phases essential for CP violation as we showed in our earlier studies in Refs. [21,40]. Based on the studies and due to the strong constraints on the parameter space of the two models, r f are expected to be small compared to r f and thus in the following we give an estimation of the upper bound on r f only.

Models with charged Higgs contributions
Two Higgs doublet models (2HDM) are simple extensions of the SM. In 2HDM, only the scalar sector of the SM is extended to include extra Higgs doublet [46,47]. In the literature, 2HDM have been classified, according to their couplings to quarks and leptons, into: 2HDM type I, II or III (for a review see Ref. [48]). The 2HDM type III (2HDM III) is of a particular interest to our study due to the presence of complex couplings of Higgs to quarks which are relevant for generating the desired CP violating weak phases. The model has five physical mass eigenstates; heavy CP-even Higgs (H 0 ), light CP-even Higgs (h 0 ), CP-odd Higgs ( A 0 ) and finally the charged Higgs (H ± ). In the model also, both Higgs doublets can couple to up-type and down-type quarks implying that the couplings of the neutral Higgs mass eigenstates can lead to flavor violation in neutral currents at tree-level. As a result, flavor changing neutral current processes can be used to strongly restraint these couplings [49,50]. We turn now to the charged Higgs couplings to the quarks. The interaction Lagrangian in this case is given as [49,50] where (31) where v u and v d denote the vacuum expectations values of the neutral component of the Higgs doublets, tan β = v u /v d and V is the CKM matrix. Extensive study of all possible constraints that can be imposed on the parameters u,d j f has been carried in Ref. [50]. We also have studied the constraints on u,d j f relevant to the process D 0 → K − π + in a previous work in Ref. [21]. Based on our study in Ref. [21], the total amplitude, including Higgs contribution, can be written as Keeping only the dominant contributions to a H ± 1,2 we find that where A recent analysis has set the bound m H ± 600 GeV independent of tan β in 2HDM II [51]. This result has been obtained after considering the most recent constraints from flavour physics and direct charged and neutral Higgs boson searches at LEP and the LHC. It should be noted that the obtained bound must be respected also for the charged Higgs mass in 2HDM III [49], Thus, for tan β = 50 and m H ± = 600 GeV we find that where we kept only the imaginary parts required for generating the weak phases and neglected the real parts of a H ± 1 and a H ± 2 as they are much smaller than the SM contributions and they are not relevant for generating weak phases. The most dominant constraints on I m( u 22 ) arise from the electric dipole moment of the neutron [50].

A new charged gauge boson as left right models
Possible extensions of the SM include models based on the gauge group SU (2) L ×SU (2) R ×U (1) B−L [52][53][54][55][56][57][58][59][60][61]. In these class of models, new complex couplings can arise due to the interactions of quarks and leptons with the new charged boson. In turn, this can affects CP violation in meson and lepton sectors. Previous analyses showed that large direct CP violation can be generated in the Charm and muon sectors if the mixing between the left and the right gauge bosons is allowed and the left-right symmetry is not manifest at unification scale [21,62,63]. Motivated by this finding, we study the impact of the new complex couplings in such particular setup of LRS model on A f of the decay process D 0 → K − π + .
The charged current mixing matrix can be parameterized as [59,62,64] here ξ is a mixing angle, W ± 1 and W ± 2 denote the mass eigenstates and ω is a weak CP violating phase. This mixing results in interactions between charged quarks and charged W bosons that reads whereV R = e iω V R . Upon integrating out W 1 in the usual way and neglecting W 2 contributions, given its mass is much higher, we can express the total amplitude of D 0 → K − π + as with The measurement of the muon decay parameter ρ, which governs the shape of the overall momentum spectrum, performed by the TWIST collaboration [65,66] can set constraint on the left right mixing angle ξ . The ρ parameter can be linked to ξ via [65]: Upon defining ζ = g R g L ξ and using the TWIST value, from their latest global fit given in Table VII in Ref. [66], ρ = 0.74960 ± 0.00019 we obtain 3.7 × 10 −3 ζ 2.3 × 10 −2 (41) which represents the allowed 2σ range of the mixing parameter ζ . We move now to discuss the allowed values ofV R ud ,V R cs . The real parts of these quark mixing matrix elements will be always suppressed by a factor ζ and thus can be neglected compare to the SM contributions. This will be the case also for the imaginary parts ofV R ud ,V R cs where they are also suppressed by the same factor ζ . However, they provide new source of the desired weak CP violating phases and thus can not be neglected.
Recently, the authors of Ref. [67] have investigated the possible bounds that can be imposed on the complex couplings of the W ± boson to right-handed quarks using lowenergy precision measurements, flavor physics and collider physics. These bounds can be applied to the couplings in general left-right symmetric model that allows mixing between the charged gauge bosons of the SU (2) R and SU (2) L as the one we consider here. The findings of the study in Ref. [67], imply that the experimental value of ( / ) K and the stringent bounds on the electric dipole moment of the neutron can allow I m(V R ud ) to be as large as 9 × 10 −4 . Moreover, the dominant constraint on ζ I m(V R cs ) arise from the process K L → π 0 e + e − and can allow ζ I m(V R cs ) to have a maximum value 7 × 10 −3 . Consequently, with the range of ζ in Eq. (41), I m(V R cs ) can have a value O(1) without violating the imposed constraints from the process K L → π 0 e + e − . Taking these values into account, we obtain | a L RS 1 | O(10 −2 ) and | a L RS 2 | O(10 −1 ). As a consequence, we find that r f O(10 −2 ) and hence, in this class of NP models |A f | O(10 −5 ). The result is one order of magnitude smaller than the experimentally measured value by LHCb collaborators.

Conclusion
In this work we have studied the CP asymmetry in the timeintegrated effective widths, A f , for the Cabibbo Favored D 0 → K − π + decay process within different models. This asymmetry is very sensitive to both the scale and the weak phases of the amplitudes generated from the radiative corrections to the SM tree-level amplitude or from the New Physics contributions. In the SM, due to the suppression of the radiative corrections to the tree-level amplitude, we have shown that |A f | O(10 −10 ) which is very suppressed compared to the recently measured value (A f ) E x p. = (1.6 ± 1.0) × 10 −4 .
It should be noted that, this experimental result shows only a 1.6 σ deviation from zero and the mode served only as a control channel in Ref. [16], where its consistence with zero has been used as a justification of the main result of that paper. Regarding the prediction in the framework of 2HDM III, we have found that |A f SM+H ± | O(10 −7 ) which is several orders of magnitude smaller than the current experimental value. Finally, in a general left-right symmetric models, allowing the mixing between the left and the right gauge bosons and adopting the scenario that the left-right symmetry is not manifest at unification scale can lead to a value |A f | O(10 −5 ) which is one order of magnitude smaller than the experimentally measured value by LHCb collaborators.