Top-quark flavor-changing $tqZ$ couplings and rare $\Delta F=1$ processes

TL;DR

This paper investigates how anomalous top-quark $tqZ$ couplings influence rare flavor-changing processes, revealing significant effects on certain Kaon and B-meson decays and providing constraints based on current experimental bounds.

Contribution

It offers a model-independent analysis of $tqZ$ couplings' impacts on low-energy flavor physics, highlighting their effects on rare decays and CP violation, with updated bounds from ATLAS data.

Findings

$tuZ$-coupling can significantly affect rare $K$ and $B$ decays.

Enhanced $BR(K_S o \mu^+ \mu^-)$ by 11% over SM.

Constraints on $BR(t o uZ)$ impact rare decay predictions.

Abstract

We model-independently study the impacts of anomalous $tqZ$ couplings ($q=u,c$), which lead to the $t\to q Z$ decays, on low energy flavor physics. It is found that the $tuZ$-coupling effect can significantly affect the rare $K$ and $B$ decays, whereas the $tcZ$-coupling effect is small. Using the ATLAS's branching ratio (BR) upper bound of $BR(t\to uZ) < 1.7\times 10^{-4}$, the influence of the anomalous $tuZ$-coupling on the rare decays can be found as follows: (a) The contribution to the Kaon direct CP violation can be up to $Re(\epsilon'/\epsilon) \lesssim 0.8 \times 10^{-3}$; (b) $BR(K^+\to \pi^+ \nu \bar \nu) \lesssim 12 \times 10^{-11}$ and $BR(K_L \to \pi^0 \nu\bar \nu)\lesssim 7.9 \times 10^{-11}$; (c) the BR for $K_S \to \mu^+ \mu^-$ including the long-distance effect can be enhanced by $11\%$ with respect to the standard model result, and (d) $BR(B_d \to \mu^+ \mu^-) \lesssim 1.97 \times 10^{-10}$. In addition, although $Re(\epsilon'/\epsilon)$ cannot be synchronously enhanced with $BR(K_L\to \pi^0 \nu \bar\nu)$ and $BR(K_S\to \mu^+ \mu^-)$ in the same region of the CP-violating phase, the values of $Re(\epsilon'/\epsilon)$, $BR(K^+ \to \pi^+ \nu \bar\nu)$, and $BR(B_d \to \mu^+ \mu^-)$ can be simultaneously increased.

Tables

Table 1. Table 1: Inputs for the numerical estimates.

$m_s=0.109$ GeV	$m_d=5.10$ MeV	$m_c=1.3$ GeV	$m_t\left(m_t\right)=165$ GeV
$m_t^{\mathrm{pole}}=172$ GeV	$m_W=80.38$ GeV	${\Gamma }_t=1.43$ GeV	$m_K=0.498$ GeV
$m_{B_d}=5.28$ GeV	$V_{ud,tb,cs}\approx 1$	$V_{td}=0.0088e^{-i{23}^{\circ }}$	$V_{ts}=-0.041$
$V_{us}=0.225$	$V_{cd}=-0.225$	${\sin }^2{\theta }_W=0.23$	$f_{\pi }=0.13$ GeV
$f_K=0.16$ GeV	$f_B=0.191$ GeV	$\left|{ϵ}_K\right|=2.228\times {10}^{-3}$	${\tau }_{K_S\left(B\right)}=89.5\left(1.52\right)\times {10}^{-12}$ s

Equations

$$BR(t\to uZ)$$

$$BR(t\to cZ)$$

$$Re({ϵ}^{\prime }/ϵ)=1.38(5.15)(4.59)\times 10^{-4},$$

$$Re({ϵ}^{\prime }/ϵ{)}_{\mathrm{SM}}=(1.9\pm 4.5)\times 10^{-4}.$$

$$BR(K^{+}\to {\pi }^{+}\nu \bar{\nu })$$

$$BR(K_L\to {\pi }^0\nu \bar{\nu })$$

$$-{\mathcal{L}}_{tqZ}$$

$$\sqrt{{\left|{\zeta }_q^L\right|}^2+{\left|{\zeta }_q^R\right|}^2}$$

$$C_{tqZ}$$

$$\sqrt{|{\zeta }_u^L{|}^2+|{\zeta }_u^R{|}^2}<0.019,$$

$$\sqrt{|{\zeta }_c^L{|}^2+|{\zeta }_c^R{|}^2}<0.022.$$

$${\mathcal{H}}_{tqZ}$$

$$Q_3$$

$$Q_7$$

$$Q_9$$

$$y_3^Z$$

$$y_9^Z$$

$$I_Z(x_t)=-\frac{1}{4}+\frac{x_t\ln x_t}{2(x_t-1)}\approx 0.693.$$

$$Re\left(\frac{{ϵ}^{\prime }}{ϵ}\right)\approx -\frac{\omega }{\sqrt{2}|{ϵ}_K|}\left[\frac{Im{A}_0}{Re{A}_0}-\frac{Im{A}_2}{Re{A}_2}\right],$$

$$y_6^Z(m_c)$$

$$y_8^Z(m_c)$$

$$Re{\left(\frac{{ϵ}^{\prime }}{ϵ}\right)}_P^Z$$

$$a_8^{(3/2)}$$

$$⟨Q_8{⟩}_2=\sqrt{2}{\left(\frac{m_K^2}{m_s(\mu )+m_d(\mu )}\right)}^2{f}_{\pi }{B}_8^{3/2}.$$

$${\mathcal{H}}_{d_i\to d_j\ell \bar{\ell }}$$

$$-\frac{G_F{V}_{t{d}_j}^{\ast }{V}_{t{d}_i}}{\sqrt{2}}\frac{\alpha }{\pi }{\bar{d}}_j{\gamma }_{\mu }{P}_L{d}_i\left[C_9^Z\bar{\ell }{\gamma }^{\mu }\ell +C_{10}^Z\bar{\ell }{\gamma }^{\mu }{\gamma }_5\ell \right],$$

$$C_L^Z=C_{10}^Z$$

$$BR(K_L\to {\pi }^0\nu \bar{\nu })$$

$$BR(K^{+}\to {\pi }^{+}\nu \bar{\nu })$$

$$X_{\mathrm{eff}}={\lambda }_t\left(X_L^{\mathrm{SM}}-s_W^2{C}_L^{Z\ast }\right),$$

$$BR(K_S\to {\mu }^{+}{\mu }^{-}{)}_{\mathrm{SD}}$$

$$BR(B_d\to {\mu }^{+}{\mu }^{-})$$

$$\times {\left|V_{td}^{\ast }{V}_{tb}\left(C_{10}^{\mathrm{SM}}+C_{10}^Z\right)\right|}^2.$$

$$\mathcal{L}={\mathcal{L}}_{\mathrm{SM}}^{(4)}+\frac{1}{\Lambda }\sum _k{C}_k^{(5)}{Q}_k^{(5)}+\frac{1}{{\Lambda }^2}\sum _k{C}_k^{(6)}{Q}_k^{(6)}+...,$$

$$\mathcal{L}$$

$$\displaystyle\left.+\left(\mbox{${\varphi^{\dagger}i\,\raisebox{5.69054pt}{\boldmath${}^{\leftrightarrow}$}D_{\mu}\,\varphi}$}\right)\left[U_{R}\gamma^{\mu}C_{\phi u}U_{R}+D_{R}\gamma^{\mu}C_{\phi d}D_{R}\right]+\left(\widetilde{\varphi}^{\dagger}iD_{\mu}\varphi\right)\left(U_{R}\gamma^{\mu}C_{\phi ud}D_{R}\right)\right\}\,,$$

$$\displaystyle\mbox{${\varphi^{\dagger}i\,\raisebox{5.69054pt}{\boldmath${}^{\leftrightarrow}$}D_{\mu}\,\varphi}$}=i\varphi^{\dagger}\left(D_{\mu}-\mbox{${\raisebox{5.69054pt}{\boldmath${}^{\leftarrow}$}D}$}_{\mu}\right)\varphi\,,~{}\mbox{${\varphi^{\dagger}i\,\raisebox{5.69054pt}{\boldmath${}^{\leftrightarrow}$}D_{\mu}^{\,I}\,\varphi}$}=i\varphi^{\dagger}\left(\tau^{I}D_{\mu}-\mbox{${\raisebox{5.69054pt}{\boldmath${}^{\leftarrow}$}D}$}_{\mu}\tau^{I}\right)\varphi\,.$$

Full text

I Introduction

Top-quark flavor changing neutral currents (FCNCs) are extremely suppressed in the standard model (SM) due to the Glashow-Iliopoulos-Maiani (GIM) mechanism Glashow:1970gm . The branching ratios (BRs) for the $t\to q(g,\gamma,Z,h)$ decays with $q=u,c$ in the SM are of order of $10^{-12}-10^{-17}$ AguilarSaavedra:2004wm ; Abbas:2015cua , and these results are far below the detection limits of LHC, where the expected sensitivity in the high luminosity (HL) LHC for an integrated luminosity of 3000 fb*-1* at $\sqrt{s}=14$ TeV is in the range $10^{-5}-10^{-4}$ ATLAS:2013hta ; ATLAS:2016 . Thus, the top-quark flavor-changing processes can serve as good candidates for investigating the new physics effects. Extensions of the SM, which can reach the HL-LHC sensitivity, can be found in Abraham:2000kx ; Eilam:2001dh ; AguilarSaavedra:2002kr ; Agashe:2006wa ; Fox:2007in ; Li:2011fza ; Li:2011af ; Gong:2013sh ; Yuan:2010vk ; tqZ ; Gaitan:2017tka ; Mandrik:2018gud .

Using the data collected with an integrated luminosity of 36.1 fb*-1* at $\sqrt{s}=13$ TeV, ATLAS reported the current strictest upper limits on the BRs for $t\to qZ$ as Aaboud:2018nyl :

[TABLE]

Based on the current upper bounds, we model-independently study the implications of anomalous $tqZ$ couplings in the low energy flavor physics. It is found that the $tqZ$ couplings through the $Z$-penguin diagram can significantly affect the rare decays in $K$ and $B$ systems, such as $\epsilon^{\prime}/\epsilon$, $K\to\pi\nu\bar{\nu}$, $K_{S}\to\mu^{+}\mu^{-}$, and $B_{d}\to\mu^{+}\mu^{-}$. Since the gluon and photon in the top-FCNC decays are on-shell, the contributions from the dipole-operator transition currents are small. In this study we thus focus on the $t\to qZ$ decays, especially the $t\to uZ$ decay.

From a phenomenological perspective, the importance of investigating the influence of these rare decays are stated as follows: The inconsistency in $\epsilon^{\prime}/\epsilon$ between theoretical calculations and experimental data was recently found based on two analyses: (i) The RBC-UKQCD collaboration obtained the lattice QCD result with Blum:2015ywa ; Bai:2015nea :

[TABLE]

where the numbers in brackets denote the errors. (ii) Using a large $N_{c}$ dual QCD Buras:1985yx ; Bardeen:1986vp ; Bardeen:1986uz ; Bardeen:1986vz ; Bardeen:1987vg , the authors in Buras:2015xba ; Buras:2015yba obtained:

[TABLE]

Note that the authors in Gisbert:2017vvj could obtain $Re(\epsilon^{\prime}/\epsilon)=(15\pm 7)\times 10^{-4}$ when the short-distance (SD) and long-distance (LD) effects are considered. Both RBC-UKQCD and DQCD results show that the theoretical calculations exhibit an over $2\sigma$ deviation from the experimental data of $Re(\epsilon^{\prime}/\epsilon)_{\rm exp}=(16.6\pm 2.3)\times 10^{-4}$, measured by NA48 Batley:2002gn and KTeV AlaviHarati:2002ye ; Abouzaid:2010ny . Based on the results, various extensions of the SM proposed to resolve the anomaly can be found in Buras:2015qea ; Buras:2015yca ; Buras:2015kwd ; Buras:2015jaq ; Tanimoto:2016yfy ; Buras:2016dxz ; Kitahara:2016otd ; Endo:2016aws ; Bobeth:2016llm ; Cirigliano:2016yhc ; Endo:2016tnu ; Bobeth:2017xry ; Crivellin:2017gks ; Bobeth:2017ecx ; Haba:2018byj ; Buras:2018lgu ; Chen:2018ytc ; Chen:2018vog ; Matsuzaki:2018jui ; Haba:2018rzf ; Aebischer:2018rrz ; Aebischer:2018quc ; Aebischer:2018csl ; Chen:2018dfc ; Chen:2018stt . We find that the direct Kaon CP violation arisen from the $tuZ$-coupling can be $\epsilon^{\prime}/\epsilon\lesssim 0.8\times 10^{-3}$ when the bound of $BR(t\to uZ)<1.7\times 10^{-4}$ is satisfied.

Unlike $\epsilon^{\prime}/\epsilon$, which strongly depends on the hadronic matrix elements, the calculations of $K^{+}\to\pi^{+}\nu\bar{\nu}$ and $K_{L}\to\pi^{0}\nu\bar{\nu}$ are theoretically clean and the SM results can be found as Bobeth:2016llm :

[TABLE]

where the QCD corrections at the next-to-leading-order (NLO) Buchalla:1993bv ; Misiak:1999yg ; Buchalla:1998ba and NNLO Gorbahn:2004my ; Buras:2005gr ; Buras:2006gb and the electroweak corrections at the NLO Buchalla:1997kz ; Brod:2008ss ; Brod:2010hi have been calculated. In addition to their sensitivity to new physics, $K_{L}\to\pi^{0}\nu\bar{\nu}$ is a CP-violating process and its BR indicates the CP-violation effect. The current experimental situations are $BR(K^{+}\to\pi^{+}\nu\bar{\nu})^{\rm exp}=(17.3^{+11.5}_{-10.5})\times 10^{-11}$ Artamonov:2008qb and $BR(K_{L}\to\pi^{0}\nu\bar{\nu})^{\rm exp}<2.6\times 10^{-8}$ Ahn:2009gb . The NA62 experiment at CERN is intended to measure the BR for $K^{+}\to\pi^{+}\nu\bar{\nu}$ with a $10\%$ precision Rinella:2014wfa ; Moulson:2016dul , and the KOTO experiment at J-PARC will observe the $K_{L}\to\pi^{0}\nu\bar{\nu}$ decay Komatsubara:2012pn ; Beckford:2017gsf . In addition, the KLEVER experiment at CERN starting in Run-4 could observe the BR of $K_{L}\to\pi^{0}\nu\bar{\nu}$ to $20\%$ precision Moulson:2018mlx . Recently, NA62 reported its first result using the 2016 taken data and found that one candidate event of $K^{+}\to\pi^{+}\nu\bar{\nu}$ could be observed, where the corresponding BR upper bound is given by $BR(K^{+}\to\pi^{+}\nu\bar{\nu})<14\times 10^{-10}$ at a $95\%$ confidence level (CL) Velghe:2018erh . We will show that the anomalous $tuZ$-coupling can lead to $BR(K^{+}\to\pi^{+}\nu\bar{\nu})\lesssim 12\times 10^{-11}$ and $BR(K_{L}\to\pi^{0}\nu\bar{\nu})\lesssim 7.9\times 10^{-11}$. It can be seen that NA62, KOTO, and KLEVER experiments can further constrain the $tuZ$-coupling using the designed sensitivities.

Another important CP violating process is $K_{S}\to\mu^{+}\mu^{-}$, where the SM prediction including the long-distance (LD) and short-distance (SD) effects is given as $BR(K_{S}\to\mu^{+}\mu^{-})=(5.2\pm 1.5)\times 10^{-12}$ Ecker:1991ru ; Isidori:2003ts ; DAmbrosio:2017klp . The current upper limit from LHCb is $BR(K_{S}\to\mu^{+}\mu^{-})<0.8(1.0)\times 10^{-9}$ at a 90%(95%) CL. It is expected that using the LHC Run-2 data, the LHCb sensitivity can be improved to $[4,\,200]\times 10^{-12}$ with 23 fb*-1* and to $[1,\,100]\times 10^{-12}$ with 100 fb*-1* Dettori:UKF2017 . Although the $tuZ$-coupling can significantly enhance the SD contribution of $K_{S}\to\mu^{+}\mu^{-}$, due to LD dominance, the increase of $BR(K_{S}\to\mu^{+}\mu^{-})_{\rm LD+SD}$ can be up to $11\%$.

It has been found that the $tuZ$-coupling-induced $Z$-penguin can significantly enhance the $B_{d}\to\mu^{+}\mu^{-}$ decay, where the SM prediction is given by $BR(B_{d}\to\mu^{+}\mu^{-})=(1.06\pm 0.09)\times 10^{-10}$ Bobeth:2013uxa . From the data, which combine the full Run I data with the results of 26.3 fb*-1* at $\sqrt{s}=13$ TeV, ATLAS reported the upper limit as $BR(B_{d}\to\mu^{+}\mu^{-})<2.1\times 10^{-10}$ Aaboud:2018mst . In addition, the result combined CMS and LHCb was reported as $BR(B_{d}\to\mu^{+}\mu^{-})=(3.9^{+1.6}_{-1.4})\times 10^{-10}$ CMS:2014xfa , and LHCb recently obtained the upper limit of $BR(B_{d}\to\mu^{+}\mu^{-})<3.4\times 10^{-10}$ Aaij:2017vad . It can be seen that the measured sensitivity is close to the SM result. We find that using the current upper limit of $BR(t\to uZ)$, the $BR(B_{d}\to\mu^{+}\mu^{-})$ can be enhanced up to $1.97\times 10^{-10}$, which is close to the ATLAS upper bound.

The paper is organized as follows: In Sec. II, we introduce the effective interactions for $t\to qZ$ and derive the relationship between the $tqZ$-coupling and $BR(t\to qZ)$. The $Z$-penguin FCNC processes induced via the anomalous $tqZ$ couplings are given in Sec. III. The influence on $\epsilon^{\prime}/\epsilon$ is shown in the same section. The $tqZ$-coupling contribution to the other rare $K$ and $B$ decays is shown in Sec. IV. A summary is given in Sec. V.

II Anomalous $tqZ$ couplings and their constraints

Based on the prescription in AguilarSaavedra:2004wm , we write the anomalous $tqZ$ interactions as:

[TABLE]

where $g$ is the $SU(2)_{L}$ gauge coupling; $c_{W}=\cos\theta_{W}$ and $\theta_{W}$ is the Weinberg angle; $P_{L(R)}=(1\mp\gamma_{5})/2$, and $\zeta^{L(R)}_{q}$ denote the dimensionless effective couplings and represent the new physics effects. In this study, we mainly concentrate the impacts of the $tqZ$ couplings on the low energy flavor physics, in which the rare $K$ and $B$ decays are induced through the penguin diagram. The rare $D$-meson processes, such as $D-\bar{D}$ mixing and $D\to\ell\bar{\ell}$, can be induced through the box diagrams; however, the processes in $D$ system can always be suppressed by taking one of the involved anomalous couplings, e.g. $tcZ$, to be small. Thus, in the following analysis, we focus on the study in the rare $K$ and $B$ decays. In order to study the influence on the Kaon CP violation, we take $\zeta^{L,R}_{q}$ as complex parameters, and the new CP violating phases are defined as $\zeta^{\chi}_{q}=|\zeta^{\chi}_{q}|e^{-i\theta^{\chi}_{q}}$ with $\chi=L,R$. The top anomalous couplings in Eq. (6) can basically arise from the dimension-six operators in the SM effective field theory (EFT), where the theory with new physics effects obeys the $SU(2)_{L}\times U(1)_{Y}$ gauge symmetry. For clarity, we show the detailed analysis for the left-handed quark couplings in Appendix. It can be found that the couplings in Eq. (6), which are generated from the SM-EFT, are not completely excluded by the low-energy flavor physics when the most general couplings are applied. The case with the strict constraints can be found in Fox:2007in . In addition to the SM-EFT Buchmuller:1985jz ; AguilarSaavedra:2009mx ; Grzadkowski:2010es , the top anomalous $tqZ$ couplings can be induced from the lower dimensional operators in the extension of the SM, such as $SU(2)$ singlet vector-like up-type quark model AguilarSaavedra:2002kr , extra dimensions Agashe:2006wa , and generic two-Higgs-doublet model Gaitan:2017tka . Hence, in this study, we take $\zeta^{\chi}_{q}$ are the free parameters and investigate the implications of the sizable $\zeta^{\chi}_{q}$ effects without exploring their producing mechanism.

Using the interactions in Eq. (6), we can calculate the BR for $t\to qZ$ decay. Since our purpose is to examine whether the anomalous $tqZ$-coupling can give sizable contributions to the rare $K$ and $B$ decays when the current upper bound of $BR(t\to qZ)$ is satisfied, we express the parameters $\zeta^{L,R}_{q}$ as a function of $BR(t\to qZ)$ to be:

[TABLE]

For the numerical analysis, the relevant input values are shown in Table 1. Using the numerical inputs, we obtain $C_{tqZ}\approx 0.40$. When $BR(t\to u(c)Z)<1.7(2.3)\times 10^{-4}$ measured by ATLAS are applied, the upper limits on $\sqrt{|\zeta^{L}_{u(c)}|^{2}+|\zeta^{R}_{u(c)}|^{2}}$ can be respectively obtained as:

[TABLE]

Since the current measured results of the $t\to(u,c)Z$ decays are close each other, the bounds on $\zeta^{\chi}_{u}$ and $\zeta^{\chi}_{c}$ are very similar. We note that BR cannot determine the CP phase; therefore, $\theta^{\chi}_{u}$ and $\theta^{\chi}_{c}$ are free parameters.

III Anomalous $tqZ$ effects on $\epsilon^{\prime}/\epsilon$

In this section, we discuss the $tqZ$-coupling contribution to the Kaon direct CP violation. The associated Feynman diagram is shown in Fig. 1, where $q=u,c$; $q^{\prime}$ and $q^{\prime\prime}$ are down type quarks, and $f$ denotes any possible fermions. That is, the involved rare $K$ and $B$ decay processes in this study are such as $K\to\pi\pi$, $K\to\pi\nu\bar{\nu}$, and $K_{S}(B_{d})\to\ell^{+}\ell^{-}$. It is found that the contributions to $K_{L}\to\pi\ell^{+}\ell^{-}$ and $B\to\pi\ell^{+}\ell^{-}$ are not significant; therefore, we do not discuss the decays in this work.

Based on the $tqZ$ couplings shown in Eq. (6), the effective Hamiltonian induced by the $Z$-penguin diagram for the $K\to\pi\pi$ decays at $\mu=m_{W}$ can be derived as:

[TABLE]

where $\lambda_{t}=V^{*}_{ts}V_{td}$; the operators $Q_{3,7,9}$ are the same as the SM operators and are defined as:

[TABLE]

with $e_{q^{\prime}}$ being the electric charge of $q^{\prime}$-quark, and the effective Wilson coefficients are expressed as:

[TABLE]

with $\alpha=e^{2}/4\pi$, $x_{t}=m^{2}_{t}/m^{2}_{W}$, and $s_{W}=\sin\theta_{W}$. The penguin-loop integral function is given as:

[TABLE]

Since $W$-boson can only couple to the left-handed quarks, the right-handed couplings $\zeta^{R}_{u,c}$ in the diagram have to appear with $m_{u(c)}$ and $m_{t}$, in which the mass factors are from the mass insertion in the quark propagators inside the loop. When we drop the small factors $m_{c,u}/m_{W}$, the effective Hamiltonian for $K\to\pi\pi$ only depends on $\zeta^{L}_{u,c}$. Since $|V_{ud}/V_{td}|$ is larger than $|V_{cs}/V_{ts}|$ by a factor of $4.67$, the dominant contribution to the $\Delta S=1$ processes is from the first term of $\eta_{Z}$ defined in Eq. (11). In addition, $V_{ud}$ is larger than $|V_{cd}|$ by a factor of $1/\lambda\sim 4.44$; therefore, the main contribution in the first term of $\eta_{Z}$ comes from the $V_{ud}\zeta^{L*}_{u}/V_{td}$ effect. That is, the anomalous $tuZ$-coupling is the main effect in our study.

Using the isospin amplitudes, the Kaon direct CP violating parameter from new physics can be estimated using Buras:2015yba :

[TABLE]

where $\omega=ReA_{2}/ReA_{0}\approx 1/22.35$ denotes the $\Delta I=1/2$ rule, and $|\epsilon_{K}|\approx 2.228\times 10^{-3}$ is the Kaon indirect CP violating parameter. It can be seen that in addition to the hadronic matrix element ratios, $\epsilon^{\prime}/\epsilon$ also strongly depends on the Wilson coefficients at the $\mu=m_{c}$ scale. It is known that the main new physics contributions to $\epsilon^{\prime}/\epsilon$ are from the $Q^{(\prime)}_{6}$ and $Q^{(\prime)}_{8}$ operators Buras:2014sba ; Buras:2015yca . Although these operators are not generated through the $tqZ$ couplings at $\mu=m_{W}$ in our case, they can be induced via the QCD radiative corrections. The Wilson coefficients at the $\mu=m_{c}$ scale can be obtained using the renormalization group (RG) evolution Buchalla:1995vs . Thus, the induced effective Wilson coefficients for $Q_{6,8}$ operators at $\mu=m_{c}$ can be obtained as:

[TABLE]

It can be seen that $y^{Z}_{6}(m_{c})$ is much smaller than $y^{Z}_{8}(m_{c})$; that is, we can simply consider the $Q_{8}$ operator contribution.

According to the $K\to\pi\pi$ matrix elements and the formulation of $Re(\epsilon^{\prime}/\epsilon)$ provided in Buras:2015yba , the $O_{8}$ contribution can be written as:

[TABLE]

where $r_{2}=\omega G_{F}/(2|\epsilon_{K}|)\approx 1.17\times 10^{-4}$ GeV*-2*, $B^{(3/2)}_{8}\approx 0.76$; $ReA^{\rm exp}_{2(0)}\approx 1.21(27.04)\times 10^{-8}$ GeV PDG , and the matrix element of $\langle Q_{8}\rangle_{2}$ is defined as:

[TABLE]

Although the $Q_{8}$ operator can contribute to the isospin $I=0$ state of $\pi\pi$, because its effect is a factor of $15$ smaller than the isospin $I=2$ state, we thus neglect its contribution.

Since the $t\to(u,c)Z$ decays have not yet been observed, in order to simplify their correlation to $\epsilon^{\prime}/\epsilon$, we use $BR(t\to qZ)\equiv{\rm Min}(BR(t\to cZ),\,BR(t\to uZ))$ instead of $BR(t\to u(c)Z)$ as the upper limit. The contours for $Re(\epsilon^{\prime}/\epsilon)^{Z}_{P}$ ( in units of $10^{-3}$) as a function of $BR(t\to qZ)$ and $\theta^{L}_{u}$ are shown in Fig. 2, where the solid and dashed lines denote the results with $\theta^{L}_{c}=-\theta^{L}_{u}$ and $\zeta^{L}_{c}=0$, respectively, and the horizontal dashed line is the current upper limit of $BR(t\to qZ)$. It can be seen that the Kaon direct CP violation arisen from the anomalous $tuZ$-coupling can reach $0.8\times 10^{-3}$, and the contribution from $tcZ$-coupling is only a minor effect. When the limit of $t\to qZ$ approaches $BR(t\to qZ)\sim 0.5\times 10^{-4}$, the induced $\epsilon^{\prime}/\epsilon$ can be as large as $Re(\epsilon^{\prime}/\epsilon)^{Z}_{P}\sim 0.4\times 10^{-3}$.

IV Z-penguin induced (semi)-leptonic $K$ and $B$ decays and numerical analysis

The same Feynman diagram as that in Fig. 1 can be also applied to the rare leptonic and semi-leptonic $K(B)$ decays when $f$ is a neutrino or a charged lepton. Because $|V_{us}/V_{ts}|\ll|V_{cs}/V_{ts}|\sim|V_{us}/V_{td}|\ll|V_{ud}/V_{td}|$, it can be found that the anomalous $tu(c)Z$-coupling contributions to the $b\to s\ell\bar{\ell}$ ($\ell=\nu,\ell^{-})$ processes can deviate from the SM result being less than $7\%$ in terms of amplitude. However, the influence of the $tuZ$ coupling on $d\to s\ell\bar{\ell}$ and $b\to d\ell\bar{\ell}$ can be over $20\%$ at the amplitude level. Accordingly, in the following analysis, we concentrate the study on the rare decays, such as $K\to\pi\nu\bar{\nu}$, $K_{S}\to\mu^{+}\mu^{-}$, and $B_{d}\to\mu^{+}\mu^{-}$, in which the channels are sensitive to the new physics effects and are theoretically clean.

According to the formulations in Bobeth:2017ecx , we write the effective Hamiltonian for $d_{i}\to d_{j}\ell\bar{\ell}$ induced by the $tuZ$ coupling as:

[TABLE]

where we have ignored the small contributions from the $tcZ$-coupling; $d_{i}\to d_{j}$ could be the $s\to d$ or $b\to d$ transition, and the effective Wilson coefficients are given as:

[TABLE]

Because $-1+4s^{2}_{W}\approx-0.08$, the $C^{Z}_{9}$ effect can indeed be neglected.

Based on the interactions in Eq. (17), the BRs for the $K_{L}\to\pi^{0}\nu\bar{\nu}$ and $K^{+}\to\pi^{+}\nu\bar{\nu}$ decays can be formulated as Buras:2015yca :

[TABLE]

where $\lambda_{c}=V^{*}_{cs}V_{cd}$, $\Delta_{EM}=-0.003$; $P_{c}(X)=0.404\pm 0.024$ denotes the charm-quark contribution Isidori:2005xm ; Mescia:2007kn ; the values of $\kappa_{L}$ and $\kappa_{+}$ are respectively given as $\kappa_{L}=(2.231\pm 0.013)\times 10^{-10}$ and $\kappa_{+}=(5.173\pm 0.025)\times 10^{-11}$ , and $X_{\rm eff}$ is defined as:

[TABLE]

with $X^{\rm SM}_{L}=1.481\pm 0.009$ Buras:2015yca . Since $K_{L}\to\pi^{0}\nu\bar{\nu}$ is a CP violating process, its BR only depends on the imaginary part of $X_{\rm eff}$. Another important CP violating process in $K$ decay is $K_{S}\to\mu^{+}\mu^{-}$, where its BR from the SD contribution can be expressed as Bobeth:2017ecx :

[TABLE]

with $C^{\rm SM}_{10}\approx-4.21$. Including the LD effect Ecker:1991ru ; Isidori:2003ts , the BR for $K_{S}\to\mu^{+}\mu^{-}$ can be estimated using $BR(K_{S}\to\mu^{+}\mu^{-})_{\rm LD+SD}\approx 4.99_{\rm LD}\times 10^{-12}+BR(K_{S}\to\mu^{+}\mu^{-})_{\rm SD}$ DAmbrosio:2017klp . Moreover, it is found that the effective interactions in Eq. (17) can significantly affect the $B_{d}\to\mu^{+}\mu^{-}$ decay, where its BR can be derived as:

[TABLE]

Because $B_{d}\to\mu^{+}\mu^{-}$ is not a pure CP violating process, the BR involves both the real and imaginary part of $V^{*}_{td}V_{tb}\left(C^{\rm SM}_{10}+C^{Z}_{10}\right)$. Note that the associated Wilson coefficient in $B_{d}\to\mu^{+}\mu^{-}$ is $C^{Z}_{10}$, whereas it is $C^{Z*}_{10}$ in the $K$ decays.

After formulating the BRs for the investigated processes, we now numerically analyze the $tuZ$-coupling effect on these decays. Since the involved parameter is the complex $\zeta^{L}_{u}=|\zeta^{L}_{u}|e^{-i\theta^{L}_{u}}$, we take $BR(t\to uZ)$ instead of $|\zeta^{L}_{u}|$. Thus, we show $BR(K_{L}\to\pi^{0}\nu\bar{\nu})$ (in units of $10^{-11})$ as a function of $BR(t\to uZ)$ and $\theta^{L}_{u}$ in Fig. 3(a), where the CP phase is taken in the range of $\theta^{L}_{u}=[-\pi,\pi]$; the SM result is shown in the plot, and the horizontal line denotes the current upper limit of $BR(t\to uZ)$. It can be clearly seen that $BR(K_{L}\to\pi^{0}\nu\bar{\nu})$ can be enhanced to $7\times 10^{-11}$ in $\theta^{L}_{u}>0$ when $BR(t\to uZ)<1.7\times 10^{-4}$ is satisfied. Moreover, the result of $BR(K_{L}\to\pi^{0}\nu\bar{\nu})\approx 5.3\times 10^{-11}$ can be achieved when $BR(t\to uZ)=0.5\times 10^{-4}$ and $\theta^{u}_{L}\ =2.1$ are used. Similarly, the influence of $\zeta^{L}_{u}$ on $BR(K^{+}\to\pi^{+}\nu\bar{\nu})$ is shown in Fig. 3(b). Since $BR(K^{+}\to\pi^{+}\nu\bar{\nu})$ involves the real and imaginary parts of $X_{\rm eff}$, unlike the $K_{L}\to\pi^{0}\nu\bar{\nu}$ decay, its BR cannot be enhanced manyfold due to the dominance of the real part. Nevertheless, the BR of $K^{+}\to\pi^{+}\nu\bar{\nu}$ can be maximally enhanced by $38\%$; even, with $BR(t\to uZ)=0.5\times 10^{-4}$ and $\theta^{u}_{L}=2.1$, the $BR(K^{+}\to\pi^{+}\nu\bar{\nu})$ can still exhibit an increase of $15\%$. It can be also found that in addition to $|\zeta^{L}_{u}|$, the BRs of $K\to\pi\nu\bar{\nu}$ are also sensitive to the $\theta^{L}_{u}$ CP-phase. Although the observed $BR(K\to\pi\nu\bar{\nu})$ cannot constrain $BR(t\to uZ)$, the allowed range of $\theta^{L}_{u}$ can be further limited.

For the $K_{S}\to\mu^{+}\mu^{-}$ decay, in addition to the SD effect, the LD effect, which arises from the absorptive part of $K_{S}\to\gamma\gamma\to\mu^{+}\mu^{-}$, predominantly contributes to the $BR(K_{S}\to\mu^{+}\mu^{-})$. Thus, if the new physics contribution is much smaller than the LD effect, the influence on $BR(K_{S}\to\mu^{+}\mu^{-})_{\rm LD+SD}=BR(K_{S}\to\mu^{+}\mu^{-})_{\rm LD}+BR(K_{S}\to\mu^{+}\mu^{-})_{\rm SD}$ from new physics may not be so significant. In order to show the $tuZ$-coupling effect, we plot the contours for $BR(K_{S}\to\mu^{+}\mu^{-})_{\rm LD+SD}$ ( in units of $10^{-12}$) in Fig. 3(c). From the result, it can be clearly seen that $BR(K_{S}\to\mu^{+}\mu^{-})_{\rm LD+SD}$ can be at most enhanced by $11\%$ with respect to the SM result, whereas the BR can be enhanced only $\sim 4.3\%$ when $BR(t\to uZ)=0.5\times 10^{-4}$ and $\theta^{L}_{u}=2.1$ are used. We note that the same new physics effect also contributes to $K_{L}\to\mu^{+}\mu^{-}$. Since the SD contribution to $K_{L}\to\mu^{+}\mu^{-}$ is smaller than the SM SD effect by one order of magnitude, we skip to show the case for the $K_{L}\to\mu^{+}\mu^{-}$ decay.

As discussed earlier that the $tcZ$-coupling contribution to the $B_{s}\to\mu^{+}\mu^{-}$ process is small; however, similar to the case in $K^{+}\to\pi^{+}\nu\bar{\nu}$ decay, the BR of $B_{d}\to\mu^{+}\mu^{-}$ can be significantly enhanced through the anomalous $tuZ$-coupling. We show the contours of $BR(B_{d}\to\mu^{+}\mu^{-})$ ( in units of $10^{-10})$ as a function of $BR(t\to uZ)$ and $\theta^{L}_{u}$ in Fig. 3(d). It can be seen that the maximum of the allowed $BR(B_{d}\to\mu^{+}\mu^{-})$ can reach $1.97\times 10^{-10}$, which is a factor of 1.8 larger than the SM result of $BR(B_{d}\to\mu^{+}\mu^{-})^{\rm SM}\approx 1.06\times 10^{-10}$. Using $BR(t\to uZ)=0.5\times 10^{-4}$ and $\theta^{L}_{u}=2.1$, the enhancement factor to $BR(B_{d}\to\mu^{+}\mu^{-})^{\rm SM}$ becomes $1.38$. Since the maximum of $BR(B_{d}\to\mu^{+}\mu^{-})$ has been close to the ATLAS upper bound of $2.1\times 10^{-10}$, the constraint from the rare $B$ decay measured in the LHC could further constrain the allowed range of $\theta^{L}_{u}$

V Summary

We studied the impacts of the anomalous $tqZ$ couplings in the low energy physics, especially the $tuZ$ coupling. It was found that the anomalous coupling can have significant contributions to $\epsilon^{\prime}/\epsilon$, $BR(K\to\pi\nu\bar{\nu})$, $K_{S}\to\mu^{+}\mu^{-}$, and $B_{d}\to\mu^{+}\mu^{-}$. Although these decays have not yet been observed in experiments, with the exception of $\epsilon^{\prime}/\epsilon$, their designed experiment sensitivities are good enough to test the SM. It was found that using the sensitivity of $BR(t\to uZ)\sim 5\times 10^{-5}$ designed in HL-LHC, the resulted $BR(K\to\pi\nu\bar{\nu})$ and $BR(B_{d}\to\mu^{+}\mu^{-})$ can be examined by the NA62, KOTO, KELVER, and LHC experiments.

According to our study, it was found that we cannot simultaneously enhance $Re(\epsilon^{\prime}/\epsilon)$, $BR(K_{L}\to\pi^{0}\nu\bar{\nu})$, and $BR(K_{S}\to\mu^{+}\mu^{-})$ in the same region of the CP violating phase, where the positive $Re(\epsilon^{\prime}/\epsilon)$ requires $\theta^{L}_{u}<0$, but the large $BR(K_{L}\to\pi^{0}\nu\bar{\nu})$ and $BR(K_{S}\to\mu^{+}\mu^{-})$ have to rely on $\theta^{L}_{u}>0$. Since $BR(K^{+}\to\pi^{+}\nu\bar{\nu})$ and $BR(B_{d}\to\mu^{+}\mu^{-})$ involve both real and imaginary parts of Wilson coefficients, their BRs are not sensitive to the sign of $\theta^{L}_{u}$. Hence, $Re(\epsilon^{\prime}/\epsilon)$, $BR(K^{+}\to\pi^{+}\nu\bar{\nu})$ and $BR(B_{d}\to\mu^{+}\mu^{-})$ can be enhanced at the same time.

Appendix A Anomalous gauge couplings from the SM-EFT

If we take the SM as an effective theory at the electroweak scale, the new physics effects should appear in terms of higher dimensional operators when the heavy fields above electroweak scale are integrated out. Thus, the effective Lagrangian with respect to the SM gauge symmetry can be generally expressed as Buchmuller:1985jz ; AguilarSaavedra:2009mx ; Grzadkowski:2010es :

[TABLE]

where ${\cal L}^{(4)}_{\rm SM}$ is the original SM; $Q^{(n)}_{k}$ are the dimension-$n$ effective operators, and $C^{(n)}_{k}$ are the associated Wilson coefficients. The top flavor-changing anomalous couplings can be generated from the dimension-6 operators, where based on the notations in Grzadkowski:2010es , the relevant operators in our study can be written as Grzadkowski:2010es :

[TABLE]

where $\varphi$ denotes the SM Higgs doublet, $Q^{T}_{L}=(U_{L},D_{L})$ is left-handed quark doublet, $D_{\mu}\varphi$ is the covariant derivative acting on $\varphi$, $\tau^{I}$ are the Pauli matrices; $\widetilde{\varphi}=i\tau_{2}\varphi^{*}$, \varphi^{\dagger}\mbox{{\raisebox{5.69054pt}{\boldmath${}^{\leftarrow}$}D}}_{\mu}\varphi=(D_{\mu}\varphi)^{\dagger}\varphi, and

[TABLE]

The flavor indices are suppressed; therefore, the Wilson coefficients $\{C_{i}\}$ are $3\times 3$ matrices. Since the top anomalous gauge couplings in this study are mainly related to the left-handed couplings, in the following discussions, we focus on the couplings to the left-handed quarks.

After electroweak symmetry breaking, the relevant $Z$ and $W$ gauge couplings to the quark weak eigenstates in Eq. (24) can be formulated as:

[TABLE]

where $<\varphi>=v/\sqrt{2}$ is the vacuum expectation value (VEV) of $\varphi$. It can be seen that the $Z$ couplings to the down-type quarks can be removed if we assume $C^{(1)}_{\phi q}=-C^{(3)}_{\phi q}\equiv-C_{qL}$. Under such circumstance, the FCNCs at the tree level could only occur in the up-type quarks. In order to use the physical quark states to express Eq. (26), we introduce the unitary matrices $U^{u,d}_{L,R}$ to diagonalize the quark mass matrices. Thus, defining $C^{\prime}_{qL}=V^{u}_{L}C_{qL}V^{u\dagger}_{L}$, Eq. (26) can be written as:

[TABLE]

where $\xi_{qL}=2C^{\prime}_{qL}+2C^{\prime\dagger}_{qL}$, and $V=V^{u}_{L}V^{d\dagger}_{L}$ is the Cabibbo-Kobayashi-Maskawa (CKM) matrix. It can be seen that the anomalous gauge couplings in the neutral current interactions are strongly correlated with those in the charged-current interactions.

It is known that the CKM matrix has a hierarchical structure, such as $V_{11(22,33)}\sim 1$, $|V_{12(21)}|\sim\lambda$, $|V_{23(32)}|\sim\lambda^{2}$, and $|(V_{13(31)}|\sim\lambda^{3}$, where $\lambda\approx 0.22$ is the Wolfenstein parameter Wolfenstein:1983yz . Since each CKM matrix element is measured well, it is necessary to examine if the sizable $t\to qZ$ FCNCs are excluded by the experimental measurements, which are dictated by the charged current interactions. Thus, in the following analysis, we concentrate on the modifications of $V_{ub}$, $V_{ts}$, and $V_{td}$. First, we consider $(\xi_{qL}V)_{ub}$ for the $b\to u$ transition effect and decompose it as:

[TABLE]

where the $V_{ub,cb}$ terms in the second line are dropped due to $V_{ub,cb}\ll V_{tb}$. In order to obtain a small effect in the $b\to u$ transition, we have to require $v^{2}(\xi_{qL})_{ut}/(\Lambda^{2})$ to be much less than $0.02$, which is the current upper limit shown in Eq. (8). Similarly, the $(\xi_{qL}V_{\rm CKM})_{ts(td)}$ factors can be expressed in terms of $\lambda$ as:

[TABLE]

If we take $(\xi_{qL})_{tu}\sim\lambda(\xi_{qL})_{tc}-\lambda^{3}(\xi_{qL})_{tt}$, i.e., the $(\xi_{qL}V)_{td}$ effect is suppressed, $(\xi_{qL}V)_{ts}$ can be rewritten as:

[TABLE]

Because $(\xi_{qL})_{tu,tc,tt}$ are taken as the free parameters, we have the degrees of freedom to obtain $v^{2}|(\xi_{qL}V)_{ts}|/(2\Lambda^{2})<|V_{ts}|\sim\lambda^{2}$ without $(\xi_{qL})_{tu(tc)}\ll 1$. Using the result, we can obtain $|\zeta^{L}_{u}|=v^{2}|(\xi_{qL})_{tu}|/\Lambda^{2}<0.021$, where the upper limit is consistent with that shown in Eq. (8). Hence, although $v^{2}(\xi_{qL})_{ut}/\Lambda^{2}$ in a general SM-EFT is bounded by the measured CKM matrix elements, $\zeta^{L}_{u}=v^{2}(\xi_{qL})_{tu}/\Lambda^{2}$ could be a free parameter and $\zeta^{L}_{u}<0.021$ is still allowed.

Acknowledgments

This work was partially supported by the Ministry of Science and Technology of Taiwan, under grants MOST-106-2112-M-006-010-MY2 (CHC).