This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

WW-Boson Mass, Electroweak Precision Tests and SMEFT

JiJi Fan [email protected] Department of Physics, Brown University, Providence, RI, 02912, USA Brown Theoretical Physics Center, Brown University, Providence, RI, 02912, USA    Lingfeng Li [email protected] Department of Physics, Brown University, Providence, RI, 02912, USA    Tao Liu [email protected] Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong S.A.R., PRC    Kun-Feng Lyu [email protected] School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, USA
Abstract

Recently the CDF collaboration at the Tevatron reported a significant discrepancy between the direct measurement of the WW-boson mass and its Standard Model (SM) prediction based on electroweak precision tests (EWPTs). In this paper we explore the potential origin of this discrepancy from physics beyond the SM. Explicitly, we work on a set of six-dimensional operators in the SM effective field theory (SMEFT) which are relevant to the EWPTs. By fitting to the data, we demonstrate that an upward shift in mWm_{W} is driven by the operator 𝒪T=12(HDμH)2\mathcal{O}_{T}=\frac{1}{2}(H^{\dagger}\overset{\text{$\leftrightarrow$}}{D}_{\mu}H)^{2} with a coefficient cT(TeV/Λ)20.01c_{T}({\rm TeV}/\Lambda)^{2}\gtrsim 0.01. This suggests that the new physics scale favored by the CDF data should be multiple TeV for tree-level effects and sub TeV for loop-level effects. One simple example is to introduce a hypercharge-free electroweak triplet scalar which can raise the cTc_{T} value at tree level. We also study the potential to further test the relevant SMEFT by measuring Higgs-coupling, mWm_{W} and other EWPTs at future circular ee+e^{-}e^{+} colliders.

I Introduction

The discovery of the Higgs boson Aad et al. (2012); Chatrchyan et al. (2012) marks a great success of the Standard Model (SM) in particle physics. It opens the door to explore mass origin of the massive elementary particles and the underlying physics for electroweak symmetry breaking (EWSB), and thus has far-reaching impacts for the development of particle physics. The measurements of the WW boson can also play a significant role in this direction. As a direct consequence of EWSB, the WW-boson mass (mWm_{W}) could be determined by the ZZ-boson mass (mZm_{Z}) and the Weinberg angle (θW\theta_{W}). Thanks to the high precision achieved in electroweak precision tests (EWPTs), mWm_{W} has been well-constrained in the SM. We denote this measured value as mWEWm_{W}^{\rm EW}. Then any deviation from mWEWm_{W}^{\rm EW} in the direct measurements of mWm_{W} may serve as a signal of physics beyond the SM (BSM).

Recently, the Collider-Detector-at-Fermilab (CDF) collaboration at Tevatron reported its updated measurement of mWm_{W} Aaltonen et al. (2022a). With 4.24×106\sim 4.24\times 10^{6} semileptonic WW-decay events, this measurement yields mW=80.4335±0.0094m_{W}=80.4335\pm 0.0094 GeV Aaltonen et al. (2022a). Here a combination of the almost quadrupled statistics and the improvements in systematic effects such as parton distribution function and lepton resolution leads to a greatly reduced uncertainty of 9.4\sim 9.4 MeV. This precision is approximately halved compared to that of the previous CDF measurement Aaltonen et al. (2012) and also exceeds the current experimental average (12\sim 12 MeV) Zyla et al. (2020). Such a WW-boson mass results in a tension of 7σ\sim 7\sigma with the latest EW global fit of mWEW=80.3545±0.0059m_{W}^{\rm EW}=80.3545\pm 0.0059 GeV de Blas et al. (2021). Caution should be taken when we attend to make an interpretation of this tension. There is a significant tension between the new CDF result and the direct measurements performed by the D0 collaboration at the Tevatron Abazov et al. (2012) and the ATLAS/LHCb collaborations at the Large Hadron Collider (LHC) Aaboud et al. (2018); Aaij et al. (2022), while the latter ones are in good agreement with mWEWm_{W}^{\rm EW}. This indicates that a better understanding of the uncertainties in these measurements is needed, to ensure a solid interpretation of the data. Notably, dedicated (HL-)LHC runs with a low instantaneous luminosity and upgraded detectors may also measure mWm_{W} with a precision 10\lesssim 10 MeV ATL (2018a); Azzi et al. (2019), falsifying or strengthening the new CDF result.

With these subtleties in mind, we would explore the potential origin of this discrepancy from the BSM physics. For this purpose, we will simply assume that the direct measurements of mWm_{W} will later converge to a value close to the CDF one Aaltonen et al. (2022a). It is already known that the BSM scenarios exist where the mWm_{W} value can deviate from the SM prediction without heavily disturbing other electroweak precision observables (EWPOs) Berthier and Trott (2015); Bjørn and Trott (2016); Corbett et al. (2021). The examples include the singlet López-Val and Robens (2014) or doublet Lopez-Val and Sola (2013) extensions of the SM Higgs sector where the corrections to mWm_{W} enter at one-loop level. A significant correction to mWm_{W} can be also achieved in supersymmetry with the sfermion-loop contributions Heinemeyer et al. (2013); Diessner and Weiglein (2019); Bagnaschi et al. (2022). Moreover, a positive correction to mWm_{W} is possible via a tree-level mixing between the SM ZZ boson and the BSM “dark photon” Curtin et al. (2015); Algueró et al. (2022). In composite Higgs theories, mWm_{W} can be shifted away from its SM prediction due to the UV-scale suppressed interactions Bellazzini et al. (2014).

In this paper, we will try to explore this discrepancy in the framework of the SM effective field theory (SMEFT) and make a prediction on its impacts for the ongoing and upcoming measurements at colliders. In the SMEFT, we parametrize the leading effects of the BSM physics above the EW scale with a set of dimension-six (6D) operators

eff=SM+iciΛ2𝒪i,\mathcal{L}_{\text{eff}}=\mathcal{L}_{\text{SM}}+\sum_{i}\frac{c_{i}}{\Lambda^{2}}\mathcal{O}_{i}\ , (1)

where SM\mathcal{L}_{\text{SM}} defines the SM, and cic_{i} and Λ\Lambda denote dimensionless Wilson coefficients and the cutoff scale of the BSM physics, respectively. Instead of presenting a comprehensive study on the whole set of 6D operators, we will focus on the key ones which are closely related to the WW-boson physics. We would view this as a first step for more refined analyses later.

We organize this paper in the following way. We will introduce the analysis formalism in Section II, which was originally developed in Chiu et al. (2018). The analysis results of the SMEFT and the predictions for the ongoing and future tests will be presented in Section III. We will conclude in Section IV.

II Analysis Formalism

The observables to be considered in this paper are summarized in Tab. 1, which include mWm_{W}, the EWPOs and one Higgs observable. Below is the list of the EWPOs:

Rb=ΓbΓhad,R=ΓhadΓ,Af=2gVfgVf+gAf,sin2θefflep=14(1gVlgAl),ΓZ,σhad0=12πmZ2ΓeΓhadΓZ2,AFBf=34AeAf(f=b,),ΓW,BRWhad.\begin{split}&R_{b}=\dfrac{\Gamma_{b}}{\Gamma_{\text{had}}},\quad R_{\ell}=\dfrac{\Gamma_{\text{had}}}{\Gamma_{\ell}},\quad A_{f}=\dfrac{2g_{V}^{f}}{g_{V}^{f}+g_{A}^{f}},\\ &\sin^{2}\theta_{\text{eff}}^{\text{lep}}=\dfrac{1}{4}\left(1-\dfrac{g^{l}_{V}}{g^{l}_{A}}\right),\quad\Gamma_{Z},\quad\sigma_{\rm had}^{0}=\frac{12\pi}{m_{Z}^{2}}\frac{\Gamma_{e}\Gamma_{\rm had}}{\Gamma_{Z}^{2}},\\ &A_{FB}^{f}=\dfrac{3}{4}A_{e}A_{f}\quad(f=b,\ell),\quad\Gamma_{W},\quad\text{BR}_{W\to{\rm had}}\ .\end{split} (2)

Here Γf=Γ(Zff¯)\Gamma_{f}=\Gamma(Z\rightarrow f\bar{f}), Γhad=Γ(Zhadrons)\Gamma_{\rm had}=\Gamma(Z\to\text{hadrons}) and ΓZ\Gamma_{Z} are the partial and total decay widths of the ZZ boson. gV(A)fg_{V(A)}^{f} is the (axial-)vector couplings of fermion ff. ΓW\Gamma_{W} and BRWhad\text{BR}_{W\to{\rm had}} denote the total decay width and branching fraction of hadronic decays of WW boson. Note that we take the experimental average at the LHC and Tevatron as the value of sin2θefflep\sin^{2}\theta_{\text{eff}}^{\text{lep}}, which is thus uncorrelated with the asymmetry observables (AA’s) measured at LEP.

Experimental Measurement Precisions at Future ee+e^{-}e^{+} Theoretical Prediction
EW Inputs mZm_{Z}(GeV) 91.1875±0.002191.1875\pm 0.0021Schael et al. (2006) ±1.0×104\pm 1.0\times 10^{-4}Abada et al. (2019) /
GF(1010GeV2)G_{F}(10^{-10}\text{GeV}^{-2}) 1166378.7±0.61166378.7\pm 0.6 Patrignani et al. (2016) ±0.6\pm 0.6 /
α1(mZ)\alpha^{-1}(m_{Z}) 127.952±0.009127.952\pm 0.009  Zyla et al. (2020) ±0.003\pm 0.003 Abada et al. (2019) /
Observables mWm_{W}(GeV) 80.4335 ±\pm 0.0094 Aaltonen et al. (2022b) ±3.9×104\pm 3.9\times 10^{-4} Abada et al. (2019) 80.3545 ±\pm 0.0059 (±\pm 0.0014)
AbA_{b} 0.923 ±\pm 0.020 Zyla et al. (2020); Schael et al. (2006) ±2.9×103\pm 2.9\times 10^{-3} Abada et al. (2019) 0.934727 (±0.000014\pm 0.000014)
AA_{\ell} (PτP_{\tau}) 0.1465 ±\pm 0.0033 Zyla et al. (2020); Schael et al. (2006) ±4.3×105\pm 4.3\times 10^{-5} Abada et al. (2019) 0.14692 (±0.00018\pm 0.00018)
AA_{\ell} (SLD) 0.1513 ±\pm 0.0021 Zyla et al. (2020); Schael et al. (2006) / 0.14692 (±0.00018\pm 0.00018)
RbR_{b} 0.21629 ±\pm 0.00066 Zyla et al. (2020); Schael et al. (2006) ±6.0×105\pm 6.0\times 10^{-5} Abada et al. (2019) 0.21588 (±0.00005\pm 0.00005)
RR_{\ell} 20.767 ±\pm 0.025 Zyla et al. (2020); Schael et al. (2006) ±1.0×103\pm 1.0\times 10^{-3} Abada et al. (2019) 20.7464 (±0.0020\pm 0.0020)
AFBbA_{\rm FB}^{b} 0.0996 ±\pm 0.0016 Zyla et al. (2020); Schael et al. (2006) ±3.0×104\pm 3.0\times 10^{-4} Abada et al. (2019) 0.10300 (±\pm 0.00013)
AFBA_{\rm FB}^{\ell} 0.0171 ±\pm 0.0010 Zyla et al. (2020); Schael et al. (2006) ±2.0×104\pm 2.0\times 10^{-4} Abada et al. (2019) 0.01619 (±0.00040\pm 0.00040)
ΓZ\Gamma_{Z}(GeV) 2.4955 ±\pm 0.00230.0023 Zyla et al. (2020) ±2.5×105\pm 2.5\times 10^{-5} Abada et al. (2019) 2.49414 (±0.00019\pm 0.00019)
σhad0\sigma_{\rm had}^{0}(nb) 41.480 ±\pm 0.0330.033 Janot and Jadach (2020) ±4.0×103\pm 4.0\times 10^{-3} Abada et al. (2019) 41.4929 (±0.0032\pm 0.0032)
ΓW\Gamma_{W}(GeV) 2.085 ±\pm 0.0420.042 Zyla et al. (2020) ±1.2×103\pm 1.2\times 10^{-3} Abada et al. (2019) 2.08782 (±0.00011\pm 0.00011)
BRW→had 0.6741 ±\pm 0.00210.0021 Zyla et al. (2020) ±6.4×105\pm 6.4\times 10^{-5} d’Enterria and Srebre (2016) 0.6748 (±0.00010\pm 0.00010)
sin2θefflep(105)\sin^{2}\theta_{\text{eff}}^{\text{lep}}(10^{-5}) (23143±25)(23143\pm 25) Aad et al. (2015); Aaij et al. (2015); Sirunyan et al. (2018); ATL (2018b); Aaltonen et al. (2018) ±0.31\pm 0.31 Abada et al. (2019) 23153.4 (±2.3\pm 2.3)
μgghγγ/μSM\mu_{ggh}^{\gamma\gamma}/\mu_{\rm SM} 1.02 ±\pm 0.11 Rossi (2022) / 1
BR(hγγh\to\gamma\gamma)/BR(hγγh\to\gamma\gamma)SM / ±0.0032\pm 0.0032 Abada et al. (2019) 1
Table 1: Input parameters, EWPOs and Higgs observable for the analyses with current data and the sensitivity projections at a future circular ee+e^{-}e^{+} collider. For the Higgs observable, we use μgghγγ/μSM\mu_{ggh}^{\gamma\gamma}/\mu_{\rm SM} in the former case and BR(hγγh\to\gamma\gamma)/BR(hγγh\to\gamma\gamma)SM in the latter case. These two are equivalent in our framework. For the uncertainties of theoretical predictions, we consider the one of mWm_{W} only in the former case, but the ones for all except the Higgs observable in the latter case. Their values have been listed outside and inside the brackets in the last column, respectively.

Totally six CP-even 6D operators will be considered in this study, which include

𝒪WB\displaystyle\mathcal{O}_{WB} =\displaystyle= ggHσaHWμνaBμν,\displaystyle gg^{\prime}H^{\dagger}{\sigma}^{a}HW^{a}_{\mu\nu}B^{\mu\nu},
𝒪T\displaystyle\mathcal{O}_{T} =\displaystyle= 12(HDμH)2,\displaystyle\frac{1}{2}(H^{\dagger}\overset{\text{$\leftrightarrow$}}{D}_{\mu}H)^{2}\,,
𝒪L(3)l\displaystyle\mathcal{O}^{(3)l}_{L} =\displaystyle= (iHσaDμH)(L¯LγμσaLL),\displaystyle(iH^{\dagger}{\sigma}^{a}\overset{\text{$\leftrightarrow$}}{D}_{\mu}H)(\bar{L}_{L}{\gamma}^{\mu}{\sigma}^{a}L_{L}),
𝒪LL(3)l\displaystyle\mathcal{O}^{(3)l}_{LL} =\displaystyle= (L¯LγμσaLL)(L¯LγμσaLL),\displaystyle(\bar{L}_{L}{\gamma}_{\mu}{\sigma}^{a}L_{L})(\bar{L}_{L}{\gamma}^{\mu}{\sigma}^{a}L_{L}),
𝒪Ll\displaystyle\mathcal{O}^{l}_{L} =\displaystyle= (iHDμH)(L¯LγμLL),\displaystyle(iH^{\dagger}\overset{\text{$\leftrightarrow$}}{D}_{\mu}H)(\bar{L}_{L}{\gamma}^{\mu}L_{L}),
𝒪Re\displaystyle\mathcal{O}^{e}_{R} =\displaystyle= (iHDμH)(l¯RγμlR).\displaystyle(iH^{\dagger}\overset{\text{$\leftrightarrow$}}{D}_{\mu}H)(\bar{l}_{R}{\gamma}^{\mu}l_{R})\,.

Here we work in the Warsaw basis Grzadkowski et al. (2010).111Our 𝒪T\mathcal{O}_{T} is slightly different from 𝒪HD=(HDμH)(HDμH)\mathcal{O}_{HD}=(H^{\dagger}D_{\mu}H)^{*}(H^{\dagger}D_{\mu}H) in the original Warsaw basis. They are related as 𝒪T=2𝒪HD12μ(HH)μ(HH)\mathcal{O}_{T}=-2\mathcal{O}_{HD}-\frac{1}{2}\partial_{\mu}\left(H^{\dagger}H\right)\partial^{\mu}\left(H^{\dagger}H\right). These operators can affect the EW physics by renormalizing wave functions, shifting the definition of EW parameters and the SM couplings, and introducing new vertices. We will follow the linear formalism developed in Chiu et al. (2018) to analyze these effects. Strictly speaking, there are additional CP-even 6D operators that could affect EWPTs such as Higgs coupling to quarks Ellis et al. (2021). Yet our choice already suffices for the leading-order discussions on the potential BSM physics to explain the new CDF result, and explore other possible measurable consequences.

To derive the relative shifts caused by these operators in the observables listed Eq. (II), we need three basic input parameters: the ZZ-boson mass mZm_{Z}, the Fermi constant GFG_{F}, and the fine-structure constant α1(mZ)\alpha^{-1}(m_{Z}). Their experimental values and uncertainties are presented in Tab. 1. For the other input parameters that contribute to the EW predictions in the SM at one or higher-loop levels, namely mhm_{h}, mtm_{t}, αs(mZ)\alpha_{s}(m_{Z}), their values are taken to be the same as those in de Blas et al. (2021). Following the procedures described in Fan et al. (2015); d’Enterria and Srebre (2016); Dubovyk et al. (2018) with the updated parameter values, we calculate the SM values for the EWPOs which agree with de Blas et al. (2021) decently, with a difference 1σ\lesssim 1\sigma. As for theoretical uncertainties, only the one for mWm_{W} is relevant. The overall uncertainty of 5.9\sim 5.9 MeV for the mWm_{W} prediction is evenly contributed by the uncertainty of mtm_{t} and the high-order corrections. By taking this effect into consideration, the predicted uncertainties for these observables also match well with de Blas et al. (2021).

To show how the 6D operators affect these EWPOs, let us consider mWm_{W} as an example. The WW-boson mass receives the contributions from these operators via the EW-parameter shifts only, which yields

ΔmWmW=δgZgZsinθWcosθWδθW12δGFGF.\displaystyle\dfrac{\Delta m_{W}}{m_{W}}=\dfrac{\delta g_{Z}}{g_{Z}}-\dfrac{\sin\theta_{W}}{\cos\theta_{W}}\delta\theta_{W}-\dfrac{1}{2}\dfrac{\delta G_{F}}{G_{F}}\ . (3)

Here θW\theta_{W} and gZg_{Z} are defined as

sin2θW\displaystyle\sin 2\theta_{W} =\displaystyle= (4πα2GFmZ2)1/2,\displaystyle\left(\dfrac{4\pi\alpha}{\sqrt{2}{G}_{F}m^{2}_{Z}}\right)^{1/2},
gZ\displaystyle g_{Z} =\displaystyle= gcosθW=4παsin2θW=2(2GFmZ2)1/2,\displaystyle\dfrac{g}{\cos\theta_{W}}=\dfrac{4\sqrt{\pi\alpha}}{\sin 2\theta_{W}}=2(\sqrt{2}G_{F}m_{Z}^{2})^{1/2}\ , (4)

where gg is the SM SU(2)SU(2) coupling. Then we have

δθW\displaystyle\delta\theta_{W} =\displaystyle= sinθWcosθW2(cos2θWsin2θW)(δααδGFGF2δmZmZ),\displaystyle\dfrac{\sin\theta_{W}\cos\theta_{W}}{2(\cos^{2}\theta_{W}-\sin^{2}\theta_{W})}\left(\dfrac{\delta\alpha}{\alpha}-\dfrac{\delta G_{F}}{G_{F}}-\dfrac{2\delta m_{Z}}{m_{Z}}\right),
δgZgZ\displaystyle\dfrac{\delta g_{Z}}{g_{Z}} =\displaystyle= 12δGFGF+δmZmZ.\displaystyle\dfrac{1}{2}\dfrac{\delta G_{F}}{G_{F}}+\dfrac{\delta m_{Z}}{m_{Z}}\ . (5)

The parameter shifts depend on the Wilson coefficients of the 6D operators as

δmZmZ\displaystyle\dfrac{\delta m_{Z}}{m_{Z}} =\displaystyle= δZZ+cTv22Λ2,δGFGF=2(cLL(3)lcL(3)l)v2Λ2,\displaystyle-\delta Z_{Z}+\dfrac{c_{T}v^{2}}{2\Lambda^{2}},\quad\dfrac{\delta G_{F}}{G_{F}}=\dfrac{2(c^{(3)l}_{LL}-c_{L}^{(3)l})v^{2}}{\Lambda^{2}},
δαα\displaystyle\dfrac{\delta\alpha}{\alpha} =\displaystyle= 2δZA,\displaystyle-2\delta Z_{A}\ , (6)

where v=246v=246GeV is the SM Higgs VEV. δZZ\delta Z_{Z} and δZA\delta Z_{A} are field renormalization factors induced by 𝒪WB\mathcal{O}_{WB}. They are given by

δZZ=v2Λ2cosθWsinθWggcWB,δZA=v2Λ2cosθWsinθWggcWB,\begin{split}\delta Z_{Z}&=\dfrac{v^{2}}{\Lambda^{2}}\cos\theta_{W}\sin\theta_{W}gg^{\prime}c_{WB},\\ \delta Z_{A}&=-\dfrac{v^{2}}{\Lambda^{2}}\cos\theta_{W}\sin\theta_{W}gg^{\prime}c_{WB}\ ,\end{split} (7)

where gg^{\prime} is the SM U(1)YU(1)_{Y} coupling. Combining these derivations, one will find the dependence of mWm_{W} on the Wilson coefficients of these 6D operators at linear level. Generalizing this discussion, we have

Δ𝒪𝒪\displaystyle\frac{\Delta\mathcal{O}}{\mathcal{O}} =\displaystyle= C𝒪WBcWBΛ2+C𝒪TcTΛ2+C𝒪3LcL(3)lΛ2\displaystyle C_{\mathcal{O}}^{WB}\dfrac{c_{WB}}{\Lambda^{2}}+C_{\mathcal{O}}^{T}\dfrac{c_{T}}{\Lambda^{2}}+C_{\mathcal{O}}^{3L}\dfrac{c_{L}^{(3)l}}{\Lambda^{2}} (8)
+\displaystyle+ C𝒪3LLcLL(3)lΛ2+C𝒪LcLlΛ2+C𝒪RcReΛ2\displaystyle C_{\mathcal{O}}^{3LL}\dfrac{c^{(3)l}_{LL}}{\Lambda^{2}}+C_{\mathcal{O}}^{L}\dfrac{c^{l}_{L}}{\Lambda^{2}}+C_{\mathcal{O}}^{R}\dfrac{c^{e}_{R}}{\Lambda^{2}}

for all EWPOs. The coefficients C𝒪C_{\mathcal{O}}’s are summarized in Tab. 2.

𝒪\mathcal{O} C𝒪WBC_{\mathcal{O}}^{WB} C𝒪TC_{\mathcal{O}}^{T} C𝒪3LC_{\mathcal{O}}^{3L} C𝒪3LLC_{\mathcal{O}}^{3LL} C𝒪LC_{\mathcal{O}}^{L} C𝒪RC_{\mathcal{O}}^{R}
mWm_{W} -0.0111 0.0433 -0.0264 0.0264 / /
AbA_{b} -0.00781 0.0142 -0.0285 0.0285 / /
RbR_{b} 0.00189 -0.00345 0.00691 -0.00691 / /
RR_{\ell} -0.00969 0.0177 -0.159 0.0353 -0.124 0.109
AFBA_{FB}^{\ell} -1.01 1.84 -2.41 3.69 1.28 1.46
AA_{\ell} -0.583 1.06 -1.38 2.13 0.739 0.843
AFBbA_{FB}^{b} -0.625 1.14 -1.50 2.28 0.784 0.894
ΓZ\Gamma_{Z} -0.0112 0.079 -0.121 0.158 -0.0113 -0.0113
σhad0\sigma_{\rm had}^{0} 0.00142 -0.00259 0.0572 -0.00519 0.152 -0.0895
ΓW\Gamma_{W} -0.0322 0.126 -0.174 0.193 / /
BRW→had / / -0.0200 / / /
sin2θefflep\sin^{2}\theta^{\text{lep}}_{\text{eff}} 0.0483 -0.0881 0.115 -0.176 -0.0612 -0.0698
μgghγγ/μSM\mu_{ggh}^{\gamma\gamma}/\mu_{\rm SM} 5.8 / / / / /
Table 2: The coefficients for calculating the 6D-operator contributions to all the observables Chiu et al. (2018). We have also neglected the effects of renormalization group running.

Notably, the main EWPOs are sensitive to four linear combinations of the six Wilson coefficients only: ξ0=1.1cWB+2cT4cL(3)l+4cLL(3)l\xi_{0}=-1.1c_{WB}+2c_{T}-4c^{(3)l}_{L}+4c^{(3)l}_{LL}, ξ±=cL(3)l±cLl\xi_{\pm}=c_{L}^{(3)l}\pm c^{l}_{L} and cRec_{R}^{e}. Explicitly, σhad\sigma_{\rm had} is sensitive to all of them; AbA_{b} and RbR_{b} depend on ξ0\xi_{0} only; AFBb,A_{\text{FB}}^{b,\ell}, AA_{\ell}, and sin2θefflep\sin^{2}\theta_{\text{eff}}^{\text{lep}} have the same dependence on ξ0\xi_{0}, ξ+\xi_{+} and cRec_{R}^{e}; and RR_{\ell} depends on ξ0\xi_{0}, ξ+\xi_{+} and cRec_{R}^{e} uniquely. This necessarily leaves two (approximately) runaway directions. These directions could be partly lifted by ΓZ\Gamma_{Z} and the three WW observables, which have different dependences on the variables beyond ξ0,±\xi_{0,\pm} and cRec_{R}^{e}. In addition, we will include the relative Higgs diphoton signal strength in the gluon fusion channel, namely μgghγγ/μSM\mu_{ggh}^{\gamma\gamma}/\mu_{\rm SM}, to enhance the constraint on 𝒪WB\mathcal{O}_{WB}. Only the ATLAS measurements Rossi (2022) will be considered here. The CMS results are similar CMS (2020) and combining them will not qualitatively change the conclusion. With this set of observables, all runaway directions can be lifted, which we will discuss more in next section.

In this paper, we plan to study the potential to further test the relevant SMEFT at future circular ee+e^{-}e^{+} colliders such as FCC-ee Abada et al. (2019) and CEPC Dong et al. (2018) also. These machines, with extraordinary integrated luminosities and state-of-the-art detector technologies, would provide unprecedented precisions for measuring mWm_{W} and other EWPOs. For concreteness, we consider the ones at FCC-ee Abada et al. (2019) and list them in Table 1 also. These projected precisions receive contributions from both the systematic and statistical uncertainties. The projected mWm_{W} precision goes below 0.50.5 MeV, improved by more than one order of magnitude compared to the CDF measurement, thanks to the dedicated beam energy scanning around the W+WW^{+}W^{-} threshold. The other EWPTs would also be improved substantially at future ee+e^{-}e^{+} colliders. Notably, we will use BR(hγγh\to\gamma\gamma)/BR(hγγh\to\gamma\gamma)SM to replace μgghγγ/μSM\mu_{ggh}^{\gamma\gamma}/\mu_{\rm SM} here as they are mutually equivalent given the six SMEFT operators in Eq. (II) while the latter is not directly measurable at ee+e^{-}e^{+} colliders. Such a choice also allows us to combine the HL-LHC and future ee+e^{-}e^{+} collider data for a further precision improvement Fan et al. (2015).

At a future ee+e^{-}e^{+} collider, the theoretical uncertainties for the considered EWPOs become more relevant. Although not playing a significant role at the current stage, these theoretical uncertainties may become comparable to the experimental ones as the latter are expected to be improved more by the time of running a future ee+e^{-}e^{+} collider. For most of the relevant EWPOs, the major factors determining their theoretical uncertainties by that time, namely the expected input parameter precisions and the intrinsic uncertainties in relation to higher-order corrections, have been discussed in d’Enterria and Srebre (2016); Freitas et al. (2019). The exceptions are σhad0\sigma_{\rm had}^{0}, ΓW\Gamma_{W} and BRW→had. For these observables, no future projections for their intrinsic uncertainties from the higher-order corrections are available. So we adopt a conservative assumption that their future values will be reduced by half compared to current ones. The projected theoretical uncertainties are summarized in Table 1.

III Analysis of SMEFT

III.1 General Analysis

Below we will analyze the impacts of the CDF mWm_{W} value and other EWPTs on the set of 6D operators listed in Eq. (II). We will fit the data in three cases, where (1) only 𝒪WB\mathcal{O}_{WB} and 𝒪T\mathcal{O}_{T}, (2) all of the six operators except 𝒪Re\mathcal{O}_{R}^{e}, and (3) all of the six operators are turned on. Predictions of all the observables at the best-fit points are summarized in Tab. 3. The values of the Wilson coefficients favored by the data are shown in Fig. 1. For completeness, we also present the 2D posterior distributions of these Wilson coefficients in Fig. 2 (case (2)) and Fig. 3 (case (3)).

Refer to caption
Refer to caption
Figure 1: Upper: Posterior distributions of the Wilson coefficients cWBc_{WB} and cTc_{T} in five fitting scenarios. The dashed and solid contours are defined at 68% and 95% C.L., respectively. The contours have been marginalized in the five- and six-operator fitting scenarios. For the six-operator case, we have also plotted the future ee+e^{-}e^{+} collider projections as the dark orange contours. Bottom: Marginalized constraints for individual Wilson coefficients at 68% C.L. in the fitting scenarios of the two-, five-, six-operator with current precision, and of six-operator at a future ee+e^{-}e^{+} collider.
Observables Case (1) Case (2) Case (3)
mWm_{W}(GeV) 80.4182 80.4335 80.4335
AbA_{b} 0.934895 0.93481 0.934944
AA_{\ell} (PτP_{\tau}) 0.14889 0.14744 0.14736
AA_{\ell} (SLD) 0.14889 0.14744 0.14736
RbR_{b} 0.21587 0.21588 0.21587
RR_{\ell} 20.7510 20.7592 20.7634
AFBbA_{\rm FB}^{b} 0.10448 0.10340 0.10335
AFBA_{\rm FB}^{\ell} 0.01657 0.01629 0.01627
ΓZ\Gamma_{Z}(GeV) 2.49818 2.49515 2.49537
σhad0\sigma_{\rm had}^{0}(nb) 41.4915 41.4729 41.4771
ΓW\Gamma_{W}(GeV) 2.09262 2.09109 2.09261
BRW→had 0.6748 0.6748 0.6749
sin2θefflep(105)\sin^{2}\theta_{\text{eff}}^{\text{lep}}(10^{-5}) 23127.7 23146.6 23147.7
μgghγγ/μSM\mu_{ggh}^{\gamma\gamma}/\mu_{\rm SM} 1.11 1.03 1.02
χ2/D.O.F\chi^{2}/{\rm D.O.F} 1.38 1.20 1.34
Table 3: Values of the EWPOs and Higgs observable at the best-fit points and corresponding χ2/D.O.F\chi^{2}/{\rm D.O.F}.

Case (1): only 𝒪WB\mathcal{O}_{WB} and 𝒪T\mathcal{O}_{T} are turned on. The results are presented as black (with the CDF mWm_{W} value Aaltonen et al. (2022a)) and gray (with the global average value Zyla et al. (2020)) contours in Fig. 1. This is equivalent to the standard EW fit with new physics encoded in the oblique parameters SS and TT Kennedy and Lynn (1989); Holdom and Terning (1990); Peskin and Takeuchi (1990, 1992); Golden and Randall (1991). Here cWBc_{WB} and cTc_{T} are related to the SS and TT parameters as

S\displaystyle S =\displaystyle= 16πcWBv2Λ23cWB(TeVΛ)2,\displaystyle 16\pi\frac{c_{WB}\,v^{2}}{\Lambda^{2}}\approx 3c_{WB}\left(\frac{{\rm TeV}}{\Lambda}\right)^{2},
T\displaystyle T =\displaystyle= cTαv2Λ27.7cT(TeVΛ)2,\displaystyle\frac{c_{T}}{\alpha}\frac{v^{2}}{\Lambda^{2}}\approx 7.7c_{T}\left(\frac{{\rm TeV}}{\Lambda}\right)^{2}\ , (9)

where v=246v=246 GeV. It is clear that with the global average of mWm_{W}, the contours are approximately centered at {0,0}\{0,0\} as expected, while the contours with the CDF mWm_{W} value take a right-upward offset from the {0,0}\{0,0\} point such that a positive cTc_{T} is strongly favored. The latter yields cT(TeV/Λ)20.015c_{T}({\rm TeV}/\Lambda)^{2}\gtrsim 0.015 at 68% C.L. Combining Tab. 2 and Eq. (8), we have

ΔmWmW=0.0111cWBΛ2+0.0433cTΛ20.0264cL(3)lΛ2+0.0264cLL(3)lΛ2.\begin{split}\dfrac{\Delta m_{W}}{m_{W}}=&-0.0111\dfrac{c_{WB}}{\Lambda^{2}}+0.0433\dfrac{c_{T}}{\Lambda^{2}}\\ &-0.0264\dfrac{c^{(3)l}_{L}}{\Lambda^{2}}+0.0264\dfrac{c^{(3)l}_{LL}}{\Lambda^{2}}\end{split}\ . (10)

A shift in mWm_{W} required for explaining the CDF measurement is thus essentially driven by the operator 𝒪T\mathcal{O}_{T}. This conclusion can be extended to other cases also.

Case (2): all of the six operators except 𝒪Re\mathcal{O}_{R}^{e} are turned on. The results are presented as blue contours in the left panel of Fig. 1 and the light-shaded bars in its right panel. In the left panel, we have marginalized the other three coefficients to obtain the contours in the (cWB,cT)(c_{WB},c_{T}) plane. Compared to the 2D fit in the previous case, the contours expand as expected, corresponding to a larger allowed parameter space. Yet an even larger positive cTc_{T} is needed to explain the CDF mWm_{W} value. This point can also be seen from the right panel. We also provide all the 2D posterior distributions of the five Wilson coefficients in Fig. 2.

Case (3): all of the six operators are turned on. The results are presented as orange contours in the left panel of Fig. 1 and the dark shaded bars in its right panel. The marginalized contours in the 2D space and the 1D bars continue to expand in this case. However, at 68% C.L., cT(TeV/Λ)20.01c_{T}({\rm TeV}/\Lambda)^{2}\gtrsim 0.01 is favored after a full marginalization, as it happens to the other cases. The 2D posterior distributions of the six Wilson coefficients are presented in Fig. 3. We could see correlations between several Wilson coefficients, such as cRec_{R}^{e}, cTc_{T} and cL3lc_{L}^{3l}.

Refer to caption
Figure 2: 2D posterior distributions of Wilson coefficients (Case (2)), obtained from the marginalized χ2\chi^{2} analysis. The contours based on the current precisions and the precisions from a future ee+e^{-}e^{+} collider are drawn in cyan and dark cyan, respectively.
Refer to caption
Figure 3: 2D posterior distributions of Wilson coefficients (Case (3)), obtained from the marginalized χ2\chi^{2} analysis. The contours based on the current precisions and the precisions from a future ee+e^{-}e^{+} collider are drawn in orange and dark orange, respectively.

III.2 Discussions

In general, cT(TeV/Λ)20.01c_{T}({\rm TeV}/\Lambda)^{2}\gtrsim 0.01 is needed to accommodate the large mWm_{W} value and its high precision reported by CDF. These results may not be easily explained by the BSM physics, even if we ignore the tension between the CDF result and other collider measurements. Firstly, the relevant cutoff scale tends to be low. If this operator is generated at one-loop level in a BSM scenario, the new mass scale is 𝒪(100)\sim{\cal O}(100) GeV, assuming order-one couplings. If this operator is generated at tree level, the new mass scale could be raised to a few TeV. In either case, the new physics scale could be well within the reach of current and near-future LHC searches. Secondly, cTc_{T} has to be positive to increase mWm_{W}. This may not always be the case in BSM scenarios that could change mWm_{W}. As an example, depending on the symmetry breaking pattern, some little Higgs models shift mWm_{W} in the opposite direction Han and Skiba (2005); Giudice et al. (2007); Chiu et al. (2018). Combining these two considerations, it is easier to seek a tree-level explanation for the large mWm_{W}. One example is augmenting the SM with a single electroweak triplet scalar field with zero hypercharge. In this case, mWm_{W} receives a positive contribution at tree level Khandker et al. (2012); Ellis et al. (2021), while the contribution to cWBc_{WB} arises at one-loop order. More work needs to be done to exhaust all possible models and check other constraints on them, which we will not explore further in this article. Note that our analysis does not apply to low-scale BSM scenarios where the SMEFT fails.

Separately, the physical origin of this mWm_{W} anomaly can be further tested indirectly with the observables which are sensitive on these operators. Such observables typically exist for the EW processes at colliders, such as di-boson productions and vector-boson fusion scatterings. Moreover, one may consider the partial decay of Higgs boson to diphoton. As shown in the data fits (especially in case (1)), cWBc_{WB} is allowed to vary at 𝒪(0.01)\mathcal{O}(0.01). Because of its high-sensitivity to cWBc_{WB} (since such decays arise from one-loop level in the SM; see Tab. 2 also), such cWBc_{WB} values would lead to a deviation in Γhγγ\Gamma_{h\to\gamma\gamma} at 10%\sim 10\% level. This is already comparable to the current precision at the LHC. At HL-LHC, this precision is expected to reach a level 4%\sim 4\% in the gluon fusion channel Cepeda et al. (2019). Such measurements certainly will probe the BSM physics along the cWBc_{WB} direction (constraining |cWB|(TeV/Λ)20.007|c_{WB}|({\rm TeV}/\Lambda)^{2}\lesssim 0.007) and hence narrow the allowed space for explaining the mWm_{W} anomaly.

Compared to the LHC and HL-LHC, the future ee+e^{-}e^{+} colliders can provide a more decisive test of the relevant SMEFT. To demonstrate this point, we present the marginalized 1D constraints and 2D posterior distributions at FCC-ee for the Wilson coefficients in Case (3) in Fig. 1 (right panel) and Fig. 3 (also see Fig. 1 (left panel)), respectively. We also present the 2D posterior distributions at FCC-ee for the Wilson coefficients in Case (2) in Fig. 2. Here the future experimental central values for the relevant observables have been assumed to have no shift from their current ones (see Tab. 1). The FCC-ee constraints are much stronger than the current ones, thanks to the improvements of both the experimental and theoretical uncertainties. In particular, from Fig. 1 (right panel), one can see that |cWB||c_{WB}| is tightly constrained to be 103\lesssim 10^{-3}, while cTc_{T}, as the potential main cause of the MWM_{W} anomaly, can be tested with more than 10σ10\sigma away from its null limit. At last, we expect comparable sensitivities to be achieved at CEPC.

IV Summary

Our paper explores the possible new physics origin of the recent CDF mWm_{W} measurement, focusing on alleviating the tension between the CDF result and EWPTs. We carry out a model-independent analysis using a subset of 6D operators which EWPOs are most sensitive to, in the SMEFT framework. We implement three different fits and show that to accommodate the CDF result, new physics has to generate 𝒪T=12(HDμH)2\mathcal{O}_{T}=\frac{1}{2}(H^{\dagger}\overset{\text{$\leftrightarrow$}}{D}_{\mu}H)^{2} with a coefficient cT(TeV/Λ)20.01c_{T}({\rm TeV}/\Lambda)^{2}\gtrsim 0.01. This suggests a new energy scale of multiple TeV for tree-level effects and sub TeV for loop-level effects, which could be searched for either through direct searches or other indirect probes.

While the CDF result clearly needs to be cross-checked with the upcoming LHC measurements, it certainly urges the particle physics community to prepare better for understanding the potential physical cause of the anomalies from indirect precision measurements. This could be crucial for the planning of future collider projects. As we have shown, a future ee+e^{-}e^{+} collider covering the ZZ-pole, WWWW threshold, Higgs factory, and tt¯t\bar{t} threshold modes could be highly relevant or even decisive in this regard.

Acknowledgements.
We thank Ayres Freitas, Ashutosh Kotwal, Kirtimaan Mohan, Simone Pagan Griso, and Keping Xie for useful discussions. J. Fan and L. Li are supported by the DOE grant DE-SC-0010010. T. Liu is supported partly by the Area of Excellence (AoE) under the Grant No. AoE/P-404/18-3, and partly by the General Research Fund (GRF) under Grant No. 16305219. Both the AoE and GRF grants are issued by the Research Grants Council of Hong Kong S.A.R. K. Lyu is partially supported by the DOE grant DE-SC0022345.

References

  • Aad et al. (2012) G. Aad et al. (ATLAS), Phys. Lett. B716, 1 (2012), eprint 1207.7214.
  • Chatrchyan et al. (2012) S. Chatrchyan et al. (CMS), Phys. Lett. B716, 30 (2012), eprint 1207.7235.
  • Aaltonen et al. (2022a) T. Aaltonen et al. (CDF), Science 376, 170 (2022a), eprint https://www.science.org/doi/pdf/10.1126/science.abk1781, URL https://www.science.org/doi/abs/10.1126/science.abk1781.
  • Aaltonen et al. (2012) T. Aaltonen et al. (CDF), Phys. Rev. Lett. 108, 151803 (2012), eprint 1203.0275.
  • Zyla et al. (2020) P. A. Zyla et al. (Particle Data Group), PTEP 2020, 083C01 (2020).
  • de Blas et al. (2021) J. de Blas, M. Ciuchini, E. Franco, A. Goncalves, S. Mishima, M. Pierini, L. Reina, and L. Silvestrini (2021), eprint 2112.07274.
  • Abazov et al. (2012) V. M. Abazov et al. (D0), Phys. Rev. Lett. 108, 151804 (2012), eprint 1203.0293.
  • Aaboud et al. (2018) M. Aaboud et al. (ATLAS), Eur. Phys. J. C 78, 110 (2018), [Erratum: Eur.Phys.J.C 78, 898 (2018)], eprint 1701.07240.
  • Aaij et al. (2022) R. Aaij et al. (LHCb), JHEP 01, 036 (2022), eprint 2109.01113.
  • ATL (2018a) (2018a).
  • Azzi et al. (2019) P. Azzi et al., CERN Yellow Rep. Monogr. 7, 1 (2019), eprint 1902.04070.
  • Berthier and Trott (2015) L. Berthier and M. Trott, JHEP 05, 024 (2015), eprint 1502.02570.
  • Bjørn and Trott (2016) M. Bjørn and M. Trott, Phys. Lett. B 762, 426 (2016), eprint 1606.06502.
  • Corbett et al. (2021) T. Corbett, A. Helset, A. Martin, and M. Trott, JHEP 06, 076 (2021), eprint 2102.02819.
  • López-Val and Robens (2014) D. López-Val and T. Robens, Phys. Rev. D 90, 114018 (2014), eprint 1406.1043.
  • Lopez-Val and Sola (2013) D. Lopez-Val and J. Sola, Eur. Phys. J. C 73, 2393 (2013), eprint 1211.0311.
  • Heinemeyer et al. (2013) S. Heinemeyer, W. Hollik, G. Weiglein, and L. Zeune, JHEP 12, 084 (2013), eprint 1311.1663.
  • Diessner and Weiglein (2019) P. Diessner and G. Weiglein, JHEP 07, 011 (2019), eprint 1904.03634.
  • Bagnaschi et al. (2022) E. Bagnaschi, M. Chakraborti, S. Heinemeyer, I. Saha, and G. Weiglein (2022), eprint 2203.15710.
  • Curtin et al. (2015) D. Curtin, R. Essig, S. Gori, and J. Shelton, JHEP 02, 157 (2015), eprint 1412.0018.
  • Algueró et al. (2022) M. Algueró, A. Crivellin, C. A. Manzari, and J. Matias (2022), eprint 2201.08170.
  • Bellazzini et al. (2014) B. Bellazzini, C. Csáki, and J. Serra, Eur. Phys. J. C 74, 2766 (2014), eprint 1401.2457.
  • Chiu et al. (2018) W. H. Chiu, S. C. Leung, T. Liu, K.-F. Lyu, and L.-T. Wang, JHEP 05, 081 (2018), eprint 1711.04046.
  • Schael et al. (2006) S. Schael et al. (SLD Electroweak Group, DELPHI, ALEPH, SLD, SLD Heavy Flavour Group, OPAL, LEP Electroweak Working Group, L3), Phys. Rept. 427, 257 (2006), eprint hep-ex/0509008.
  • Abada et al. (2019) A. Abada et al. (FCC), Eur. Phys. J. C 79, 474 (2019).
  • Patrignani et al. (2016) C. Patrignani et al. (Particle Data Group), Chin. Phys. C40, 100001 (2016).
  • Aaltonen et al. (2022b) T. Aaltonen et al. (CDF), Science 376, 170 (2022b).
  • Janot and Jadach (2020) P. Janot and S. Jadach, Phys. Lett. B 803, 135319 (2020), eprint 1912.02067.
  • d’Enterria and Srebre (2016) D. d’Enterria and M. Srebre, Phys. Lett. B 763, 465 (2016), eprint 1603.06501.
  • Aad et al. (2015) G. Aad et al. (ATLAS), JHEP 09, 049 (2015), eprint 1503.03709.
  • Aaij et al. (2015) R. Aaij et al. (LHCb), JHEP 11, 190 (2015), eprint 1509.07645.
  • Sirunyan et al. (2018) A. M. Sirunyan et al. (CMS), Eur. Phys. J. C 78, 701 (2018), eprint 1806.00863.
  • ATL (2018b) (2018b).
  • Aaltonen et al. (2018) T. A. Aaltonen et al. (CDF, D0), Phys. Rev. D 97, 112007 (2018), eprint 1801.06283.
  • Rossi (2022) E. Rossi (ATLAS), PoS EPS-HEP2021, 580 (2022).
  • Grzadkowski et al. (2010) B. Grzadkowski, M. Iskrzynski, M. Misiak, and J. Rosiek, JHEP 10, 085 (2010), eprint 1008.4884.
  • Ellis et al. (2021) J. Ellis, M. Madigan, K. Mimasu, V. Sanz, and T. You, JHEP 04, 279 (2021), eprint 2012.02779.
  • Fan et al. (2015) J. Fan, M. Reece, and L.-T. Wang, JHEP 09, 196 (2015), eprint 1411.1054.
  • Dubovyk et al. (2018) I. Dubovyk, A. Freitas, J. Gluza, T. Riemann, and J. Usovitsch, Phys. Lett. B 783, 86 (2018), eprint 1804.10236.
  • CMS (2020) (2020).
  • Dong et al. (2018) M. Dong et al. (CEPC Study Group) (2018), eprint 1811.10545.
  • Freitas et al. (2019) A. Freitas et al. (2019), eprint 1906.05379.
  • Kennedy and Lynn (1989) D. C. Kennedy and B. W. Lynn, Nucl. Phys. B 322, 1 (1989).
  • Holdom and Terning (1990) B. Holdom and J. Terning, Phys. Lett. B 247, 88 (1990).
  • Peskin and Takeuchi (1990) M. E. Peskin and T. Takeuchi, Phys. Rev. Lett. 65, 964 (1990).
  • Peskin and Takeuchi (1992) M. E. Peskin and T. Takeuchi, Phys. Rev. D 46, 381 (1992).
  • Golden and Randall (1991) M. Golden and L. Randall, Nucl. Phys. B 361, 3 (1991).
  • Han and Skiba (2005) Z. Han and W. Skiba, Phys. Rev. D 72, 035005 (2005), eprint hep-ph/0506206.
  • Giudice et al. (2007) G. F. Giudice, C. Grojean, A. Pomarol, and R. Rattazzi, JHEP 06, 045 (2007), eprint hep-ph/0703164.
  • Khandker et al. (2012) Z. U. Khandker, D. Li, and W. Skiba, Phys. Rev. D 86, 015006 (2012), eprint 1201.4383.
  • Cepeda et al. (2019) M. Cepeda et al., CERN Yellow Rep. Monogr. 7, 221 (2019), eprint 1902.00134.