Probing top-quark operators with precision electroweak measurements

  • In the standard model effective field theory, operators involving the top quark are generally difficult to probe and can generate sizable loop contributions to electroweak precision observables measured by past and future lepton colliders. Could the high precision of electroweak measurements compensate for loop suppression and provide competitive reaches on these operators? Would the inclusion of these contributions introduce too many additional parameters for a meaningful global electroweak analysis to be performed? In this paper, we perform a detailed phenomenological study to address these two important questions. Focusing on eight dimension-6 operators that generate anomalous couplings between electroweak gauge bosons and third-generation quarks, we calculate their one loop contributions to $ e^+e^- \to f\bar{f} $ processes, both on and off the Z-pole, and the $ e^-e^+ \to WW $ process. A global analysis is performed with these eight operators and those that contribute to the above processes at tree level using measurements at the LEP, SLC, and several low energy experiments. We find that although current electroweak precision measurements are sensitive to the one-loop effects of top-quark operators, it is difficult to separate them from the operators that contribute at tree level, making a global analysis rather challenging. Under further assumptions (for instance, new physics contributes to only third generation quark operators and the S and T parameters), competitive reaches may be obtained in a global fit. Another important finding of our study is that the two operators that generate the dipole interactions of the bottom quark have a significant impact on the Z-pole measurements and should not be omitted. We also discuss the implications of the recently reported W-boson mass measurement at the CDF for our results. Finally, we estimate the reaches of future lepton colliders in probing top-quark operators with precision electroweak measurements.
  • 加载中
  • [1] A. V. Manohar, Lect. Notes Phys. 479, 311-362 (1997), arXiv:hep-ph/9606222
    [2] I. Z. Rothstein, TASI lectures on effective field theories, 8 2003. hep-ph/0308266
    [3] D. B. Kaplan, Five lectures on effective field theory, 10 2005. nucl-th/0510023
    [4] C. P. Burgess, Ann. Rev. Nucl. Part. Sci. 57, 329-362 (2007) doi: 10.1146/annurev.nucl.56.080805.140508
    [5] S. Weinberg, PoS CD09, 001 (2009), arXiv:0908.1964
    [6] W. Buchmuller and D. Wyler, Nucl. Phys. B 268, 621-653 (1986) doi: 10.1016/0550-3213(86)90262-2
    [7] B. Grzadkowski, M. Iskrzyński, M. Misiak et al., Jou. of High. Ener. Phys. 10, 085 (2010), arXiv:1008.4884 doi: 10.1007/JHEP10(2010)085
    [8] I. Brivio and M. Trott, Phys. Rept. 793, 1-98 (2019), arXiv:1706.08945 doi: 10.1016/j.physrep.2018.11.002
    [9] A. Falkowski and F. Riva, JHEP 02, 039 (2015), arXiv:1411.0669 doi: 10.1007/JHEP02(2015)039
    [10] A. Efrati, A. Falkowski, and Y. Soreq, JHEP 07, 018 (2015), arXiv:1503.07872 doi: 10.1007/JHEP07(2015)018
    [11] A. Falkowski, M. Gonzalez-Alonso, A. Greljo et al., Phys. Rev. Lett. 116(1), 011801 (2016), arXiv:1508.00581 doi: 10.1103/PhysRevLett.116.011801
    [12] G. Durieux, C. Grojean, J. Gu et al., JHEP 09, 014 (2017), arXiv:1704.02333 doi: 10.1007/JHEP09(2017)014
    [13] T. Barklow, K. Fujii, S. Jung et al., Phys. Rev. D 97(5), 053003 (2018), arXiv:1708.08912 doi: 10.1103/PhysRevD.97.053003
    [14] J. Ellis, C. W. Murphy, V. Sanz et al., JHEP 06, 146 (2018), arXiv:1803.03252 doi: 10.1007/JHEP06(2018)146
    [15] G. Durieux, M. Perelló, M. Vos et al., JHEP 10, 168 (2018), arXiv:1807.02121 doi: 10.1007/JHEP10(2018)168
    [16] G. Durieux, J. Gu, E. Vryonidou et al., Chin. Phys. C 42(12), 123107 (2018), arXiv:1809.03520 doi: 10.1088/1674-1137/42/12/123107
    [17] A. Falkowski and D. Straub, JHEP 04, 066 (2020), arXiv:1911.07866 doi: 10.1007/JHEP04(2020)066
    [18] J. de Blas et al., JHEP 01, 139 (2020), arXiv:1905.03764 doi: 10.1007/JHEP01(2020)139
    [19] J. De Blas, G. Durieux, C. Grojean et al., JHEP 12, 117 (2019), arXiv:1907.04311 doi: 10.1007/JHEP12(2019)117
    [20] G. Durieux, A. Irles, V. Miralles et al., The electro-weak couplings of the top and bottom quarks — Global fit and future prospects, JHEP 12 (2019) 98, arXiv: 1907.10619, [Erratum: JHEP 01, 195 (2021)]
    [21] N. P. Hartland, F. Maltoni, E. R. Nocera et al., JHEP 04, 100 (2019), arXiv:1901.05965 doi: 10.1007/JHEP04(2019)100
    [22] I. Brivio, S. Bruggisser, F. Maltoni et al., JHEP 02, 131 (2020), arXiv:1910.03606 doi: 10.1007/JHEP02(2020)131
    [23] S. Dawson, S. Homiller, and S. D. Lane, Phys. Rev. D 102(5), 055012 (2020), arXiv:2007.01296 doi: 10.1103/PhysRevD.102.055012
    [24] J. Ellis, M. Madigan, K. Mimasu et al., JHEP 04, 279 (2021), arXiv:2012.02779
    [25] S. Jung, J. Lee, M. Perelló et al., Phys. Rev. D 105(1), 016003 (2022), arXiv:2006.14631 doi: 10.1103/PhysRevD.105.016003
    [26] J. J. Ethier, G. Magni, F. Maltoni et al. (SMEFiT Collaboration), JHEP 11, 089 (2021), arXiv:2105.00006
    [27] E. d. S. Almeida, A. Alves, O. J. P. Éboli (SMEFiT Collaboration), Phys. Rev. D 105(1), 013006 (2022), arXiv:2108.04828 doi: 10.1103/PhysRevD.105.013006
    [28] V. Miralles, M. M. López, M. M. Llácer (SMEFiT Collaboration), JHEP 02, 032 (2022), arXiv:2107.13917
    [29] M. Cepeda et al., CERN Yellow Rep. Monogr. 7, 221-584 (2019), arXiv:1902.00134
    [30] M. McCullough, An Indirect Model-Dependent Probe of the Higgs Self-Coupling, Phys. Rev. D 90(1), 015001 (2014), [Erratum: Phys. Rev. D 92, 039903 (2015)] arXiv: 1312.3322
    [31] G. Degrassi, P. P. Giardino, F. Maltoni et al., JHEP 12, 080 (2016), arXiv:1607.04251 doi: 10.1007/JHEP12(2016)080
    [32] S. Di Vita, C. Grojean, G. Panico et al., JHEP 09, 069 (2017), arXiv:1704.01953 doi: 10.1007/JHEP09(2017)069
    [33] S. Di Vita, G. Durieux, C. Grojean et al., JHEP 02, 178 (2018), arXiv:1711.03978
    [34] E. Vryonidou and C. Zhang, JHEP 08, 036 (2018), arXiv:1804.09766
    [35] R. Contino, T. Kramer, M. Son et al., JHEP 05, 074 (2007), arXiv:hep-ph/0612180
    [36] R. Barbieri, D. Buttazzo, F. Sala et al., JHEP 07, 181 (2012), arXiv:1203.4218
    [37] R. Barbieri, D. Buttazzo, F. Sala et al., JHEP 05, 069 (2013), arXiv:1211.5085
    [38] A. Falkowski, M. González-Alonso, and K. Mimouni, JHEP 08, 123 (2017), arXiv:1706.03783
    [39] L. Alasfar, J. de Blas, and R. Gröber, Higgs probes of top quark contact interactions and their interplay with the Higgs self-coupling, arXiv: 2202.02333
    [40] T. Aaltonen et al. (CDF Collaboration), Science 376(6589), 170-176 (2022) doi: 10.1126/science.abk1781
    [41] S. Schael et al. (ALEPH, DELPHI, L3, OPAL, LEP Electroweak Collaboration), Phys. Rept. 532, 119-244 (2013), arXiv:1302.3415 doi: 10.1016/j.physrep.2013.07.004
    [42] P. A. Zyla et al. (Particle Data Group Collaboration), PTEP 2020(8), 083C01 (2020)
    [43] K. Hagiwara, S. Ishihara, R. Szalapski et al., Phys. Rev. D 48, 2182-2203 (1993)
    [44] Z. Han and W. Skiba, Phys. Rev. D 71, 075009 (2005), arXiv:hep-ph/0412166 doi: 10.1103/PhysRevD.71.075009
    [45] L. Berthier and M. Trott, JHEP 05, 024 (2015), arXiv:1502.02570
    [46] L. Berthier and M. Trott, JHEP 02, 069 (2016), arXiv:1508.05060
    [47] M. E. Peskin and T. Takeuchi, Phys. Rev. D 46, 381-409 (1992)
    [48] C. Zhang, N. Greiner, and S. Willenbrock, Phys. Rev. D 86, 014024 (2012), arXiv:1201.6670 doi: 10.1103/PhysRevD.86.014024
    [49] B. Yan and C. P. Yuan, Phys. Rev. Lett. 127(5), 051801 (2021), arXiv:2101.06261 doi: 10.1103/PhysRevLett.127.051801
    [50] B. Yan, Z. Yu, and C. P. Yuan, Phys. Lett. B 822, 136697 (2021), arXiv:2107.02134 doi: 10.1016/j.physletb.2021.136697
    [51] H. T. Li, B. Yan, and C. P. Yuan, Jet Charge: A new tool to probe the anomalous \begin{document}$Zb\bar{b}$\end{document} couplings at the EIC, arXiv: 2112.07747
    [52] H. Dong, P. Sun, B. Yan et al., Phys. Lett. B 829, 137076 (2022), arXiv:2201.11635 doi: 10.1016/j.physletb.2022.137076
    [53] J. Alwall, R. Frederix, S. Frixione et al., JHEP 07, 079 (2014), arXiv:1405.0301
    [54] A. van Hameren, Computer Physics Communications 182(11), 2427-2438 (2011) doi: 10.1016/j.cpc.2011.06.011
    [55] S. Dawson and P. P. Giardino, Phys. Rev. D 101(1), 013001 (2020), arXiv:1909.02000 doi: 10.1103/PhysRevD.101.013001
    [56] S. Dawson and P. P. Giardino, Flavorful Electroweak Precision Observables in the Standard Model Effective Field Theory, arXiv: 2201.09887
    [57] J. V. Allaby, U. Amaldi, G. Barbiellini et al., Zeitschrift für Physik C Particles and Fields 36(4), 611-628 (1987)
    [58] A. Blondel, P. Böckmann, H. Burkhardt et al., Zeitschrift für Physik C Particles and Fields 45(3), 361-379 (1990)
    [59] K. S. McFarland et al., The European Physical Journal C-Particles and Fields 1(3), 509-513 (1998)
    [60] G. Zeller, K. McFarland, T. Adams et al., Phys. Revi. Lett. 88(9), 091802 (2002) doi: 10.1103/PhysRevLett.88.091802
    [61] C. Wood, S. Bennett, D. Cho et al., Science 275(5307), 1759-1763 (1997) doi: 10.1126/science.275.5307.1759
    [62] M. T. Gericke and Q. Collaboration, The q weak p experiment a test for physics beyond the standard model via a precision measurement of the proton weak charge, in “AIP Conference Proceedings”, vol. 1149, pp. 237-240, American Institute of Physics. 2009
    [63] D. Wang, K. Pan, R. Subedi et al., Nature 506(7486), 67 (2014) doi: 10.1038/nature12964
    [64] D. Geiregat, G. Wilquet, U. Binder et al., Phys. Lett. B 245(2), 271-275 (1990) doi: 10.1016/0370-2693(90)90146-W
    [65] A. Collaboration, L. Collaboration, L. E. W. G. Collaboration, O. Collaboration, D. Collaboration et al., Electroweak Measurements in Electron-Positron Collisions at W-Boson-Pair Energies at LEP, Physics Reports 532, no. DESY-2014-00955, 119-244 (2013)
    [66] The LEP Collaborations: ALEPH Collaboration, DELPHI Collaboration, L3 Collaboration, OPAL Collaboration, and the LEP Electroweak Working Group, A Combination of Preliminary Electroweak Measurements and Constraints on the Standard Model, arXiv e-prints, Dec 2006, arXiv: hep-ex/0612034.
    [67] Y. Inoue, A. Miyamoto, E. Nakano et al., The European Physical Journal C-Particles and Fields 18(2), 273-282 (2000)
    [68] M. Baak, J. Cúth, J. Haller et al. (Gfitter Group Collaboration), Eur. Phys. J. C 74, 3046 (2014), arXiv:1407.3792
    [69] M. Kobel et al. (Two Fermion Working Group Collaboration), Two-Fermion Production in Electron-Positron Collisions, in “Proceedings, Monte Carlo Workshop: Report of the working groups on precision calculation for LEP-2 physics: CERN, Geneva, Switzerland, March 12-13, June 25-26, October 12-13 Oct 1999”. 2000. arXiv: hep-ph/0007180
    [70] M. Caffo, H. Czyz, and E. Remiddi, Nuovo Cim. A 110, 515-536 (1997), arXiv:hep-ph/9704443
    [71] TOPAZ Collaboration, Y. Inoue et al., Eur. Phys. J. C 18, 273-282 (2000), arXiv:hep-ex/0012033
    [72] L. Berthier and M. Trott, “Consistent constraints on the standard model effective field theory”, Journal of High Energy Physics 2016 Feb (2016).
    [73] P. Artoisenet, R. Frederix, O. Mattelaer et al., Automatic spin-entangled decays of heavy resonances in monte carlo simulations, Journal of High Energy Physics 2013 Mar (2013)
    [74] M. W. Grunewald et al., JHEP, (2000), arXiv:hep-ph/0005309
    [75] S. Jadach, W. Placzek, M. Skrzypek et al., Comput. Phys. Commun. 140, 475-512 (2001), arXiv:hep-ph/0104049 doi: 10.1016/S0010-4655(01)00296-X
    [76] S. Gori, J. Gu, and L.-T. Wang, JHEP 04, 062 (2016), arXiv:1508.07010
    [77] R. Barbieri, A. Pomarol, R. Rattazzi et al., Nucl. Phys. B 703, 127-146 (2004), arXiv:hep-ph/0405040 doi: 10.1016/j.nuclphysb.2004.10.014
    [78] J. D. Wells and Z. Zhang, JHEP 01, 123 (2016), arXiv:1510.08462
    [79] A. Falkowski and D. Straub, Flavourful smeft likelihood for higgs and electroweak data, Journal of High Energy Physics 2020 Apr (2020)
    [80] M. Farina, G. Panico, D. Pappadopulo et al., Phys. Lett. B 772, 210-215 (2017), arXiv:1609.08157 doi: 10.1016/j.physletb.2017.06.043
    [81] E. Bagnaschi, J. Ellis, M. Madigan et al., SMEFT Analysis of mW , arXiv: 2204.05260
    [82] J. Gu, Z. Liu, T. Ma et al., Speculations on the W-Mass Measurement at CDF, arXiv: 2204.05296
    [83] J. Fan, L. Li, T. Liu et al., W-Boson Mass, Electroweak Precision Tests and SMEFT, arXiv: 2204.04805
    [84] C.-T. Lu, L. Wu, Y. Wu et al., Electroweak Precision Fit and New Physics in light of W Boson Mass, arXiv: 2204.03796
    [85] A. Strumia, Interpreting electroweak precision data including the W-mass CDF anomaly, arXiv: 2204.04191
    [86] R. S. Gupta, Running away from the T-parameter solution to the W mass anomaly, arXiv: 2204.13690
    [87] E. d. S. Almeida, A. Alves, O. J. P. Eboli et al., Impact of CDF-Ⅱ measurement of MW on the electroweak legacy of the LHC Run Ⅱ, arXiv: 2204.10130
    [88] R. Balkin, E. Madge, T. Menzo et al., On the implications of positive W mass shift, arXiv: 2204.05992
    [89] M. Endo and S. Mishima, New physics interpretation of W-boson mass anomaly, arXiv: 2204.05965
    [90] P. Asadi, C. Cesarotti, K. Fraser et al., Oblique Lessons from the W Mass Measurement at CDF Ⅱ, arXiv: 2204.05283
    [91] M. Dong et al. (CEPC Study Group Collaboration), CEPC Conceptual Design Report: Volume 2 - Physics & Detector”, arXiv: 1811.10545
    [92] A. Abada et al. (FCC Collaboration), Eur. Phys. J. C 79(6), 474 (2019) doi: 10.1140/epjc/s10052-019-6904-3
    [93] V. Bresó-Pla, A. Falkowski, and M. González-Alonso, JHEP 08, 021 (2021), arXiv:2103.12074
    [94] M. Diehl and O. Nachtmann, Z. Phys. C 62, 397-412 (1994) doi: 10.1007/BF01555899
    [95] J. Erler and M. J. Ramsey-Musolf, Prog. Part. Nucl. Phys. 54, 351-442 (2005), arXiv:hep-ph/0404291 doi: 10.1016/j.ppnp.2004.08.001
    [96] V. Shtabovenko, J. Phys. Conf. Ser. 762(1), 012064 (2016), arXiv:1604.06709
    [97] R. Mertig, M. Bohm, and A. Denner, Comput. Phys. Commun. 64, 345-359 (1991) doi: 10.1016/0010-4655(91)90130-D
    [98] J. Kublbeck, M. Bohm, and A. Denner, Comput. Phys. Commun. 60, 165-180 (1990) doi: 10.1016/0010-4655(90)90001-H
    [99] T. Han, Collider phenomenology: Basic knowledge and techniques, in “Physics in D >= 4. Proceedings, Theoretical Advanced Study Institute in elementary particle physics, TASI 2004, Boulder, USA, June 6-July 2, 2004”, pp. 407-454. 2005. arXiv: hep-ph/0508097
  • 加载中

Figures(14) / Tables(7)

Get Citation
Yiming Liu, Yuhao Wang, Cen Zhang, Lei Zhang and Jiayin Gu. Probing Top-quark Operators with Precision Electroweak Measurements[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac82e1
Yiming Liu, Yuhao Wang, Cen Zhang, Lei Zhang and Jiayin Gu. Probing Top-quark Operators with Precision Electroweak Measurements[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac82e1 shu
Milestone
Received: 2022-05-15
Article Metric

Article Views(896)
PDF Downloads(37)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Probing top-quark operators with precision electroweak measurements

  • 1. Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
  • 2. School of Physics, Nanjing University, Nanjing 210093, China
  • 3. Department of Physics and Center for Field Theory and Particle Physics, Fudan University, Shanghai 200438, China; Key Laboratory of Nuclear Physics and Ion-beam Application (MOE), Fudan University, Shanghai 200433, China
  • 4. School of Physical Sciences, University of Chinese Academy of Science, Beijing 100049, China
  • 5. Center for High Energy Physics, Peking University, Beijing 100871, China

Abstract: In the standard model effective field theory, operators involving the top quark are generally difficult to probe and can generate sizable loop contributions to electroweak precision observables measured by past and future lepton colliders. Could the high precision of electroweak measurements compensate for loop suppression and provide competitive reaches on these operators? Would the inclusion of these contributions introduce too many additional parameters for a meaningful global electroweak analysis to be performed? In this paper, we perform a detailed phenomenological study to address these two important questions. Focusing on eight dimension-6 operators that generate anomalous couplings between electroweak gauge bosons and third-generation quarks, we calculate their one loop contributions to $ e^+e^- \to f\bar{f} $ processes, both on and off the Z-pole, and the $ e^-e^+ \to WW $ process. A global analysis is performed with these eight operators and those that contribute to the above processes at tree level using measurements at the LEP, SLC, and several low energy experiments. We find that although current electroweak precision measurements are sensitive to the one-loop effects of top-quark operators, it is difficult to separate them from the operators that contribute at tree level, making a global analysis rather challenging. Under further assumptions (for instance, new physics contributes to only third generation quark operators and the S and T parameters), competitive reaches may be obtained in a global fit. Another important finding of our study is that the two operators that generate the dipole interactions of the bottom quark have a significant impact on the Z-pole measurements and should not be omitted. We also discuss the implications of the recently reported W-boson mass measurement at the CDF for our results. Finally, we estimate the reaches of future lepton colliders in probing top-quark operators with precision electroweak measurements.

    HTML

    I.   INTRODUCTION
    • The standard model (SM), despite its enormous success, is generally considered an effective theory with a cutoff that can be as low as a few TeVs. Considerable effort is being devoted to constructing and studying extensions of the SM that predict new particles with masses of approximately the TeV scale. To date, direct searches for such new particles at colliders have been unsuccessful. A complementary approach to direct searches is indirect searches, where precise measurements of SM processes are compared with SM predictions, and an observed deviation from the latter is a strong indication that the process may receive virtual contributions of heavy new particles. A powerful model-independent framework to identify, constrain, and parametrize potential deviations with respect to SM predictions is the standard model effective field theory (SMEFT) [15]. Assuming electroweak symmetry breaking is linearly realized as in the SM, and the new physics scale (usually denoted as Λ) is significantly larger than the electroweak scale, the SM Lagrangian is augmented by a series of higher dimensional operators suppressed by powers of $ 1/\Lambda $. Here, we focus on the effects of the leading operators (denoted as $ Q_i $) that preserve baryon and lepton numbers, which are of dimension six.

      $ \begin{equation} \mathcal{L}_\mathrm{SMEFT}=\mathcal{L}_\mathrm{SM}+ \sum\limits_{i}\frac{c_{i}Q_{i}}{\Lambda^{2}}, \end{equation} $

      (1)

      where Λ is the energy scale of new physics, and $ c_{i} $ are the dimensionless Wilson coefficients. $ Q_{i} $ spans the entire space of dimension-6 operators [68]. If experiments detect significant deviations from SM predictions, the SMEFT could help characterize their possible origin and guide direct searches for new physics. In the absence of any significant deviation, the SMEFT can be used to systematically constrain the scales of different BSM physics scenarios.

      Many dimension-6 operators have been extensively studied in various experiments. Among them, operators that modify electroweak processes are stringently constrained by the precision measurements of the Z and W bosons at lepton colliders. Global analyses have been performed for these operators at both current and future colliders [911]. Although similar analyses have also been conducted for the Higgs and top sectors (and their combinations with the EW sector) [1228], several corresponding operators are not well constrained for a number of reasons. A well-known example is the operator $ (H^\dagger H)^3 $, which modifies the trilinear Higgs coupling. Even at the high luminosity LHC, its coefficient is only probed at the order-one level with the measurement of the double-Higgs process [29]. Interestingly, the measurement of single-Higgs processes, where the $ (H^\dagger H)^3 $ operator contributes at the one-loop order, offers competitive reaches on it due to their better measurement precision (in particular, the Higgsstrahlung process at future lepton colliders) [3033]. Similarly, operators that generate anomalous gauge couplings of the top quark are generally less-well constrained at the tree level and contribute to many Higgs and electroweak processes at the one-loop order [16, 25, 34]. Some of these contributions are also enhanced by the large top mass. From the UV point of view, sizable contributions to top-quark operators are generated in many plausible new physics models. A typical example is the composite Higgs model with partial compositeness [35]. A large mixing between the top quark and strong sector is required to generate a large top mass, which also tends to generate large contributions to 3rd-generation quark operators. Strong constraints from flavor observables can be avoided by imposing approximate $ U(2)^3 $ symmetry on the first two generation quarks [36], which can also be naturally embedded in composite Higgs models [37]. Therefore, with the high precision of electroweak measurements, it is entirely possible that one-loop effects from top-quark operators are non-negligible, despite loop suppression, and should be included in the electroweak analysis. This calls for two important questions: First, are the one-loop contributions of these top-quark operators sufficiently sizable for the electroweak measurements to have a sensitivity that is comparable to, or even better than, that from the LHC, which probes them directly? Second, by introducing additional degrees of freedom to the electroweak analysis, are we still able to obtain meaningful bounds in a global framework? In other words, is it possible to separate the effects of top-quark operators from those of tree-level electroweak ones?

      In this paper, we perform a comprehensive global analysis with current precision electroweak data to answer these two questions. We focus on the effects of eight dimension-6 operators that generate anomalous couplings between electroweak gauge bosons and third-generation quarks. These include four operators that modify the SM gauge couplings of top and bottom quarks, and four operators that introduce dipole interactions. We include Z-pole measurements from the LEP and SLD, the measurements of the $ e^+e^-\to WW $ process at LEP2, and measurements of several low energy scattering processes that are also sensitive to 4-fermion operators [38]. We perform global analyses with the eight operators above and other tree-level operators that contribute to these processes. For the latter, we impose the flavor universality condition to reduce the size of the parameter space. Interestingly, we also find that the tree level effects of the bottom dipole operators, although suppressed by the bottom mass, may be even larger than the one-loop contributions of the other operators, and these tree level effects must be included for consistency. As a first step toward a more complete global analysis, we do not consider the effects of any top-quark-related four-fermion operators in our study. These operators introduce additional degrees of freedom and, in many cases, are notoriously difficult to separate from other top operators [15, 20, 28]. Recent studies also showed that they have a significant impact on Higgs processes [39].

      The remainder of this paper is organized as follows: In Sec. II, we lay out the theoretical framework of our study, including the operators we consider, the corresponding tree-level and one-loop contributions to the electroweak processes, and details of the Monte Carlo simulation we use to obtain some of the results. In Sec. III, we provide a detailed description of the experimental inputs used in our analysis. Our results are presented in Sec. IV. We consider both the general framework and a more restrictive "semi-universal" scenario. Results with the new CDF-II W mass measurement [40] are presented, and its implication is also discussed. Our projections for future lepton colliders are provided in Sec. V. Finally, we draw conclusion in Sec. VI. More details of our analysis and additional results are presented in Appendix A.

    II.   THEORY FRAMEWORK

      A.   SMEFT framework

    • We work in a global SMEFT framework and consider all dimension-6 operators that contribute to current electroweak measurements (listed in Sec. III). Our standpoint is that third generation quarks are special. Operators involving third generation quarks are thus separated from the other dimension-6 operators, and their one-loop contributions are considered in addition to possible tree-level contributions. For the other operators, we consider only the tree-level contributions and impose $U(2)_{u} \bigotimes U(2)_{d} \bigotimes U(2)_{q} \bigotimes$$ U(3)_ {l} \bigotimes U(3)_e $ flavor symmetry. This setup allows us to investigate the impact of third-generation-quark operators while maintaining a relatively small parameter space. As mentioned earlier, four-fermion operators involving the top quark are not included in our study. We leave a more general analysis with four-fermion top-quark operators and fewer flavor assumptions to future studies.

      Under the above assumptions, the dimension-6 operators involved in our study are summarized in Table 1, where the Warsaw basis is used [7]. These operators are divided into two classes. The first class involves third generation quarks, which are

      $ \psi ^ { 2 } \varphi ^ { 3 } $$ X^{3} $$ \varphi ^ { 4 } D ^ { 2 } $
      $ Q_{u\varphi}^{ i j}=(\varphi^\dagger \varphi)(\bar q_i u_j \widetilde\varphi) $$ Q_W=\epsilon^{IJK} W_\mu^{I\nu} W_\nu^{J\rho} W_\rho^{K\mu} $$ Q_{\varphi D}=\left(\varphi^\dagger D^\mu\varphi\right)^\star \left(\varphi^\dagger D_\mu\varphi\right) $
      $ Q_{d\varphi}^{ i j}=(\varphi^\dagger \varphi)(\bar q_i d_j \varphi) $
      $ \psi ^ { 2 } \varphi ^ { 2 } D $$ \psi ^ { 2 } X \varphi $$ X ^ { 2 } \varphi ^ { 2 } $
      $ Q_{\varphi l}^{ij(1)}=\left(\varphi^{\dagger}i\stackrel{\leftrightarrow}{D}_{\mu}\varphi\right)\left(\overline{l}_{i}\gamma^{\mu}l_{j} \right) $$ Q_{uW}^{i j}=(\bar q_{i}\sigma^{\mu\nu}u_{j})\tau^{I}\widetilde{\varphi}W_{\mu\nu}^{I} $$ Q_{\varphi WB}=\varphi^{\dagger}\tau^{I} \varphi W^{I}_{\mu\nu} B^{\mu\nu} $
      $ Q_{\varphi l}^{ij(3)}=\left(\varphi^{\dagger}i\stackrel{\leftrightarrow}{D}_{\mu}^{I}\varphi \right)\left(\overline {l}_{i}\tau^{I}\gamma^{\mu}l_{j} \right) $$ Q_{uB}^{ij}=(\bar q_{i}\sigma^{\mu\nu}u_{j}) \widetilde\varphi B_{\mu\nu} $
      $ Q_{\varphi e}^{ i j}=\left(\varphi ^ { \dagger } i \stackrel { \leftrightarrow } { D } _ { \mu } \varphi \right)(\bar e_i \gamma^\mu e_j) $$ Q_{dW}^{ i j}=(\bar q_i \sigma^{\mu\nu} d_j) \tau^I \varphi\, W_{\mu\nu}^I $
      $ Q_{\varphi q}^{ i j(1)}=\left(\varphi ^ { \dagger } i \stackrel { \leftrightarrow } { D } _ { \mu } \varphi \right)(\bar q_i \gamma^\mu q_j) $$ Q_{dB }^{i j}=(\bar q_i \sigma^{\mu\nu} d_j) \varphi B_{\mu\nu} $
      $ Q_{\varphi q}^{ i j(3)}=\left(\varphi ^ { \dagger } i \stackrel { \leftrightarrow } { D } _ { \mu } ^ { I } \varphi \right)(\bar q_i \tau^I \gamma^\mu q_j) $
      $ Q_{\varphi u}^{ i j}=\left(\varphi ^ { \dagger } i \stackrel { \leftrightarrow } { D } _ { \mu } \varphi \right)(\bar u_i \gamma^\mu u_j) $
      $ Q_{\varphi d}^{ i j}=\left(\varphi ^ { \dagger } i \stackrel { \leftrightarrow } { D } _ { \mu } \varphi \right)(\bar d_i \gamma^\mu d_j) $
      $ Q_{\varphi u d}^{ i j}=i(\widetilde\varphi^\dagger D_\mu \varphi)(\bar u_i \gamma^\mu d_j) $
      $ (\overline { L } L) (\overline { L } L) $$ (\overline { R } R) (\overline { R } R) $$ (\overline { L } L) (\overline { R } R) $
      $ Q_{ll}^{ p r st }=(\bar l_p \gamma_\mu l_r)(\bar l_s \gamma^\mu l_t) $$ Q_{ee}^{ p r s t}=(\bar e_p \gamma_\mu e_r)(\bar e_s \gamma^\mu e_t) $$ Q_{le}^{ p r s t}=(\bar l_p \gamma_\mu l_r)(\bar e_s \gamma^\mu e_t) $
      $ Q_{lq }^{p r s t(1)}=(\bar l_p \gamma_\mu l_r)(\bar q_s \gamma^\mu q_t) $$ Q_{eu}^{ p r s t}=(\bar e_p \gamma_\mu e_r)(\bar u_s \gamma^\mu u_t) $$ Q_{lu}^{ p r s t}=(\bar l_p \gamma_\mu l_r)(\bar u_s \gamma^\mu u_t) $
      $ Q_{lq }^{p r s t(3)}=(\bar l_p \gamma_\mu \tau^I l_r)(\bar q_s \gamma^\mu \tau^I q_t) $$ Q_{ed}^{ p r s t}=(\bar e_p \gamma_\mu e_r)(\bar d_s\gamma^\mu d_t) $$ Q_{ld }^{p r s t}=(\bar l_p \gamma_\mu l_r)(\bar d_s \gamma^\mu d_t) $
      $ Q_{qe }^{p r s t}=(\bar q_p \gamma_\mu q_r)(\bar e_s \gamma^\mu e_t) $

      Table 1.  Operators in the Warsaw basis [8] that are used in our study. The indices i, j, and p, r, s, t label the fermion generation.

      $ \begin{equation} Q^{(1)}_{\varphi Q} \,,\; \; \; \; Q^{(3)}_{\varphi Q} \,,\; \; \; \; Q_{\varphi t} \,,\; \; \; \; Q_{\varphi b} \,,\; \; \; \; Q_{tW} \,,\; \; \; \; Q_{tB} \,,\; \; \; \; Q_{bW} \,,\; \; \; \; Q_{bB} \,, \end{equation} $

      (2)

      where, instead of presenting the flavor indices $ ij=33 $, we make the replacements $ q\to Q $, $ u\to t $, and $ d\to b $. Here, $ Q,\; t,\; b $ denote the $ S U(2) $ doublet $ (t,\,b)_L $ and the singlets $ t_R $and $ b_R $, respectively. The corresponding Wilson coefficients $ c_i/\Lambda^2 $ follow the same notation. The four operators $ Q^{(1)}_{\varphi Q} $, $ Q^{(3)}_{\varphi Q} $, $ Q_{\varphi t} $, and $ Q_{\varphi b} $ modify the SM gauge couplings between 3rd generation quarks and the W and Z bosons. There is another operator, $ Q_{\varphi tb} = {\rm i}(\widetilde\varphi^\dagger D_\mu \varphi)(\bar t_i\gamma^\mu b_j) $, which generates a right-handed $ Wtb $ coupling. It contributes only to the W-boson self energy, and this contribution is strongly suppressed by the bottom mass. For this reason, we do not include it in our analysis. The remaining four operators $ Q_{tW} $, $ Q_{tB} $, $ Q_{bW} $, and $ Q_{bB} $ generate dipole interactions between 3rd generation quarks and the gauge bosons. Owing to their different helicities, a fermion mass insertion is required to generate an interference term with the SM amplitude. For the top quark, this contribution may be sizable and should be included in the analysis. The impact of the bottom-quark dipole operators is more subtle — while their contributions to the $ Zb\bar{b} $ vertex are suppressed by the small (but non-negligible) bottom mass, these are tree level contributions, which turn out to be comparable to, or even larger than, the one-loop contributions of top operators. As such, they are also included in our analysis. For all eight operators, their leading contributions are included into all the observables in a consistent manner. That is, if an operator already contributes to a process at the tree level (such as the bottom-quark operators to $ e^+e^-\to b\bar{b} $), we shall only consider its tree-level contribution; otherwise, its one-loop contribution is included.

      The second class contains all other operators that contribute to the $ e^+e^- \to f \bar{f} $ and $ e^+e^- \to WW $ processes at the tree level, including those that modify the propagators of Z and W bosons, their couplings to fermions (excluding 3rd generation quarks), and triple gauge couplings,

      $ \begin{aligned}[b] &Q^{(1)}_{\varphi q} \,,\; \; \; \; Q^{(3)}_{\varphi q} \,,\; \; \; \; Q_{\varphi u} \,,\; \; \; \; Q_{\varphi d} \,,\; \; \; \; Q^{(1)}_{\varphi l} \,,\; \; \; \; Q^{(3)}_{\varphi l} \,,\; \; \; \; Q_{\varphi e} \,, \\ & Q'_{ll} \,,\; \; \; \; Q_{\varphi D} \,,\; \; \; \; Q_{\varphi WB} \,,\; \; \; \; Q_W\,, \end{aligned} $

      (3)

      as well as 4-fermion operators that directly contribute to $ e^+e^-\to f\bar{f} $ processes and several low energy scattering processes,

      $ \begin{aligned}[b] &Q_{qe} \,,\; \; \; \; Q_{eu} \,,\; \; \; \; Q_{ed} \,,\; \; \; \; Q^{(1)}_{lq} \,,\; \; \; \; Q^{(3)}_{lq} \,,\; \; \; \; Q_{lu} \,,\; \; \; \; Q_{ld} \,, \\ & Q_{ll} \,,\; \; \; \; Q_{ee} \,,\; \; \; \; Q_{le} \,, \end{aligned} $

      (4)

      where the flavor indices are omitted due to the flavor symmetries we impose. Note that we distinguish $ Q'_{ll} \equiv Q^{1221}_{ll} $ (which contributes to the μ decay) from $ Q_{ll} \equiv Q^{11ii}_{ll} $ (which contributes to $ e^+e^- \to l^+ l^- $), even though they are not independent in the universal flavor case. This is because many new physics models that contribute to the latter (such as a flavor-diagonal $ Z' $ boson) do not contribute to the former. In addition, we include the 4-fermion operators in the $ e^+e^- \to b\bar{b} $ process to the second class and only consider their tree-level effects. They are

      $ \begin{equation} Q^{(1)}_{lQ}, \,\; \; \; Q^{(3)}_{lQ}, \,\; \; \; Q_{lb} \,,\; \; \; \; Q_{eQ} \,,\; \; \; \; Q_{eb} \,, \end{equation} $

      (5)

      where we again make the replacements $ q\to Q $, $ u\to t $, and $ d\to b $ instead of presenting the 3rd generation flavor indices. For convenience, we also use the following combinations of operators instead of the original ones for 3rd generation quarks:

      $ \begin{aligned}[b] Q^{(+)}_{\varphi Q} \equiv Q^{(1)}_{\varphi Q} + Q^{(3)}_{\varphi Q} \,, & \quad Q^{(+)}_{l Q} \equiv Q^{(1)}_{l Q} + Q^{(3)}_{l Q} \,, \\ Q^{(-)}_{\varphi Q} \equiv Q^{(1)}_{\varphi Q} - Q^{(3)}_{\varphi Q} \,, & \quad Q^{(-)}_{l Q} \equiv Q^{(1)}_{l Q} - Q^{(3)}_{l Q} \,, \end{aligned} $

      (6)

      and their Wilson coefficients follow the same labeling system. $ Q^{(-)}_{l Q} $ contributes only to a contact $ e^+e^- t\bar{t} $ interaction and is not considered in our analysis. In total, 33 Wilson coefficients are included in our analysis, and we consider their leading contributions for each observable, which are at the one-loop level if there is no tree level contribution.

      It is well known that in the Warsaw basis, the operator coefficients that contribute to Z-pole observables at the tree level exhibit flat directions such that a global SMEFT fit with only the Z-pole observables (and the W mass measurement) cannot be closed. These flat directions are lifted by the measurement of the $ e^+e^-\to WW $ process (see, for example, [9]). However, because the measurement precision of $ e^+e^-\to WW $ at the LEP is significantly worse than that of the Z-pole, large correlations remain among many operator coefficients. This may obscure the impact of the 3rd-generation-quark operators investigated in this study. To resolve this issue, in most parts of our analysis, we use a slightly different basis obtained by replacing the operators $ Q_{\varphi D} $ and $ Q_{\varphi WB} $ with

      $ \begin{aligned}[b] Q_{D \varphi W} =& {\rm i} D_{\mu} \phi^{\dagger} \sigma_{a} D_{v} \phi W^{a \mu v}\,, \\ Q_{D \varphi B} =& {\rm i} D_{\mu} \phi^{\dagger} D_{v} \phi B^{\mu v} \,, \end{aligned} $

      (7)

      which do not contribute to the Z-pole observables. The translation to the Warsaw basis is given by

      $ \begin{aligned}[b] Q_{D \varphi B} \to& -\frac{g^{\prime}}{4} Q_{\varphi B}+\frac{g^{\prime}}{2} \sum_{\psi} Y_{\psi} Q_{\varphi \psi}^{(1)}+\frac{g^{\prime}}{4} Q_{\varphi \square}\\&+g^{\prime} Q_{\varphi D}-\frac{g}{4} Q_{\varphi W B} \,, \\ Q_{D \varphi W} \to & \; \frac{g}{4} \sum_{F} Q_{\varphi F}^{(3)}+\frac{g}{4}\left(3 Q_{\varphi \square}+8 \lambda_{\phi} Q_{\varphi}-4 \mu_{\phi}^{2}\left(\phi^{\dagger} \phi\right)^{2}\right) \\ &+\frac{g}{2}\left(y_{i j}^{e}\left(Q_{e \varphi}\right)_{i j}+y_{i j}^{d}\left(Q_{d \varphi}\right)_{i j}+y_{i j}^{u}\left(Q_{u \varphi}\right)_{i j}+\rm { h.c. }\right) \\ &-\frac{g^{\prime}}{4} Q_{\varphi W B}-\frac{g}{4} Q_{\varphi W} \,. \end{aligned} $

      (8)

      We use the following input parameters in our analysis [42]:

      $ \begin{aligned}[b]& \alpha=\frac{1}{127.9},\,\; \; \; \; m_{Z}=91.1876\,~{\rm{Ge}}{\rm{V}}, \,\; \; \; \; \\& G_{F}=1.166379\times 10^{-5}\,~{\rm{Ge}}{\rm{V}}, \\ & m_{b}=4.7\,~{\rm{Ge}}{\rm{V}},\,\; \; \; \; m_{t}=172.5\,~{\rm{Ge}}{\rm{V}} \,. \end{aligned} $

      (9)
    • B.   Loop contributions

    • The tree level contributions of the SMEFT dimension-6 operators to the $ e^+e^- \to f\bar{f} $ and $ e^+e^- \to WW $ processes have been well studied in the past [911, 4346], and we follow these references in our analysis. Here we focus on the one-loop contributions of the 3rd-generation-quark operators in Eq. (2). First, they generate universal contributions to the gauge boson self-energies, as shown in Fig. 1. They have both direct and indirect effects. The indirect effects are the contributions to the measurements of the three input parameters α, $ m^2_Z $, and $ G_F $. As with the tree level dimension-6 contributions, they change the "inferred SM values" of these parameters. Their effects can be parameterized as

      Figure 1.  Gauge boson self-energy correction caused by the operators in Eq. (2). The black dot is the dimension-6 vertex.

      $ \begin{aligned}[b] \alpha&=\alpha _{O}\left(1+\Pi^{'}_{\gamma\gamma}(0)\right) \,, \\ m^{2}_{Z}&=m^{2} _{ZO}+\Pi_{ZZ}(m^{2}_{Z}) \,, \\ G_{F}&=G _{FO}\left(1-\frac{\Pi_{WW}(0)}{m_{W}^{2}}\right)\,, \end{aligned} $

      (10)

      where α, $ m^{2}_{Z} $, and $ G_{F} $ take the values in Eq. (9), $ \alpha _{O} $, $ m^{2} _{ZO} $, and $ G _{FO} $ are the renormalized SM parameters, which are entered into the calculations of EW observables, and $ \Pi^{(i)}_{XY}(q^2) $ parameterizes the one-loop corrections to the gauge boson self-energies [47], for which we only consider the contributions from the operators in Eq. (2). Their definitions are shown in Sec. A.3. A renormalized $ s^2_W \equiv \sin^2\theta_W $ can be defined as

      $ \begin{aligned}[b] s^{2} _{WO}&=\frac{1}{2}\left(1-\sqrt{1-\frac{4\pi\alpha _{O}}{\sqrt{2}G _{FO}m _{ZO}^{2}}}\right)\\ &=s_{W}^{2}\left(1-\frac{c^{2}_{W}}{c^{2}_{W}-s^{2}_{W}}\left(\Pi^{'}_{\gamma\gamma}(0)+\frac{\Pi_{WW}(0)}{m^{2}_{W}}-\frac{\Pi_{ZZ}(m^{2}_{Z})}{m^{2}_{Z}}\right)\right), \end{aligned} $

      (11)

      where $ s_{W}^{2}=\dfrac{1}{2}\left(1-\sqrt{1-\dfrac{4\pi\alpha}{\sqrt{2}G_{F}m_{Z}^{2}}}\right) $. The W and Z self-energies also directly enter the observables. The contributions to $ e^+e^- \to f\bar{f} $ as well as the W and Z decay rates can be characterized by the tree-level neutral and charged current interactions. At tree level, they are given by (assuming no tree-level contributions from dimension-6 operators)

      $ \begin{aligned}[b] M_{NC}&=e^{2}_{O}\frac{QQ^{'}}{q^{2}}+\frac{e^{2}_{O}}{s_{WO}^{2}c_{WO}^{2}}(I_{3}-s_{WO}^{2}Q)\frac{1}{q^{2}-m_{ZO}^{2}}(I^{'}_{3}-s_{WO}^{2}Q^{'}) \,, \\ M_{CC}&=\frac{e^{2}_{O}}{2s_{WO}^{2}}I_+\frac{1}{q^{2}-m_{WO}^{2}}I_- \,, \end{aligned} $

      (12)

      where $ I^{(')}_{3} $and $ Q^{(')} $ are the $ S U(2) $ and electric charges of the external fermions, respectively, $ I_{+}, I_{-} $ are the isospin-raising and isospin-lowering matrices, respectively, and $ \alpha _{O} $, $ m^{2} _{ZO} $, $ G _{FO} $, and $ s^{2} _{WO} $ are the renormalized SM parameters in Eqs. (10) and (11). With the direct contributions from the W and Z self-energies, Eq. (12) is modified to

      $ \begin{aligned}[b] M^{1}_{NC}&=e_{*}^{2}\frac{QQ^{'}}{q^{2}}+\frac{e_{*}^{2}}{s_{W^{*}}^{2}c_{W^{*}}^{2}}(I_{3}-s_{W^{*}}^{2}Q)\frac{Z_{Z^{*}}}{q^{2}-m_{Z^{*}}^{2}}(I^{'}_{3}-s_{W^{*}}^{2}Q^{'}) \,, \\ M^{1}_{CC}&=\frac{e_{*}^{2}}{2s_{W^{*}}^{2}}I_+\frac{Z_{W^{*}}}{q^{2}-m_{W^{*}}^{2}}I_- \,, \end{aligned} $

      (13)

      where $ m^2_{W*} (q^2) $, $ m^2_{Z*} (q^2) $, $ Z_{W*} (q^2) $, $ Z_{Z*} (q^2) $, $ s^2_{W*} (q^2) $, and $ e^2_* (q^2) $ generally depend on $ q^2 $ of the propagators, and their expressions are also listed in Sec. A.3 [48]. In addition to the W and Z self-energies, the 3rd-generation-quark operators also directly modify the $ Zb\bar{b} $ vertex, as shown in Fig. 2, thus generating non-universal effects. Finally, they also contribute to the processes $ e^+e^-\to WW $ and $ e^{+}e^{-}\to e^{+}e^{-} $. The calculation of these effects are rather complicated, and we rely on Monte Carlo integration with MadGraph5_aMC@NLO [53]. The method of reweighting [54] is used to generate weighted events. The SM input parameters used in Monte Carlo are the same as those in Eq. (9). We also check that the statistical uncertainties due to simulation sample size are sufficiently small and can be neglected. The numerical results of our calculation can be found in the GitHub repository.

      Figure 2.  Diagram of top loop modification to $ bb $ production.

      The one-loop contributions to EW observables from dimension-6 operators were also recently studied in Refs. [55, 56]. In Ref. [55], the one-loop order QCD and electroweak corrections to the Z and W pole observables are computed under the flavor universality assumption. In contrast, our study focuses on 3rd generation quarks. By also computing the one-loop contributions of the first two generations of quarks and then imposing the flavor universality assumption, we manage to reproduce the results in Ref. [55] on the contributions of $ O_{\varphi t} $, $ O_{tW} $, and $ O_{tB} $ to the $ Z\rightarrow bb $ decay width. For $ O_{\varphi Q}^{(-)} $, our result turns out to be approximately four times larger than that in Ref. [55], which is likely due to the different choice of cutoff energy. In Ref. [56], the one-loop corrections of flavor non-universal 4-fermion interactions were investigated, which complements our study.

    III.   EXPERIMENT
    • The data used in this study are mainly from the LEP experiment. Data from low energy experiments, such as CHARM [57], CDHS [58], CCFR [59], NuTeV [60], APV [61], QWEAK [62], and PVDIS [63] are also included. These can be divided into the following categories:

      ● Precision electroweak measurement: We use the Z-pole measurement from the LEP experiment, W mass measurement taken from the combined results of the PDG Group and branching ratio (BR) information from the LEP.

      $ e^+ e^- \to f \bar{f} $: Measurements of electron collision to quark and lepton pairs from the LEP experiment are collected. Several observables measured by TRISTAN at 58 GeV are also included.

      $ e^+ e^- \to W^+ W^- \to 4f $: The total cross-section and differential cross-section measurements of $ e^+ e^- \to W^+ W^- \to 4f $ from the LEP experiment are included in this study.

      ● Low energy measurement: The neutral-current parameters for ν-hadron and $ \nu-e $ processes measured using deep inelastic scattering (DIS) experiments, including CDHS and CHARM in CERN, and CCFR and NuTeV in FermiLab. The parameters for electron scattering are measured using the Cs atomic parity violation (APV) experiment, electron-deuteron DIS (eDIS) experiment PVDIS, and Qweak, which measures the weak charge of protons. Additionally, the cross-section measurement of trident production ($ \nu_\mu \gamma^* \rightarrow \nu_\mu \mu^+ \mu^- $) over its SM prediction by CHARM and CDHS is included.

      For the data mentioned above, we take the center values, uncertainties, and correlations. The correlations between different categories and experiments are taken as zero. For those with unidentifiable correlations, we make the assumption that they are not correlated with each other. Table 2 shows the varieties of observables included in this study and their references.

      ExperimentObservablesReference
      Low energyCHARM/CDHS/CCFR/NuTeV/
      APV/QWEAK/PVDIS
      Effective couplings[42, 64]
      Z-poleLEP/SLCTotal decay width $\Gamma_Z$
      Hadronic cross-section $\sigma_{\rm had}$
      Ratio of decay width $R_f$
      Forward-Backward asymmetry $A_{\rm FB}^f$
      Polarized asymmetry $A_f$
      [10]
      W-poleLHC/Tevatron/LEP/SLCTotal decay width $ \Gamma_W $
      W branching ratios ${\rm Br}(W\rightarrow l v_l)$
      Mass of W boson $ M_W $
      [10, 42]
      $ ee\rightarrow qq $LEP/TRISTANHadronic cross-section $\sigma_{\rm had}$
      Ratio of cross-section $ R_f $
      Forward-Backward asymmetry for b/$c A_{\rm FB}^f$
      [65, 66, 67]
      $ ee \rightarrow ll $LEPCross-section $ \sigma_f $
      Forward-Backward asymmetry $A_{\rm FB}^f$
      Differential cross-section ${ {\rm d}\sigma_f}/{ {\rm d}{\rm cos}\theta}$
      [65]
      $ ee\rightarrow WW $LEPCross-section $ \sigma_{WW} $
      Differential cross-section ${ {\rm d}\sigma_{WW} }/{ {\rm d}{\rm cos}\theta}$
      [65]

      Table 2.  Observables used in this study.

      Table A1 in the appendix lists the type-I operators that have an impact on the above observables. Besides these, five operators ($ Q_{\varphi WB} $, $ Q^{ijji}_{ll} $, $ Q_{\varphi D} $, $ Q^{11(3)}_{\varphi l} $, and $ Q^{22(3)}_{\varphi l} $) have an impact on all the obersvables because they affect the electroweak parameters g, $ g^{'} $, and v. In the following part, details about several of the categories involved in the experiments will be discussed. These details also affect the generation of Monte Carlo events. Additionally, the uncertainties on SM calculation and experimental measurement will be discussed.

      ProcessesObservablesWilson Coefficient
      Neutrino DIS and APV$ g _ { L V } ^ { \nu e },g _ { L A } ^ { \nu e } $$ c^{(1)}_{\varphi l},c_{\varphi e},c_{ll},c_{le} $
      $g _ { A V } ^ { e u } + 2g _ { A V } ^ { e d } $
      $ 2g _ { A V } ^ { e u } - g _ { A V } ^ { e d }$
      $ 2g _ { V A } ^ { e u } - g _ { V A } ^ { e d }$
      $c^{(1)}_{\varphi q}$,$c^{(3)}_{\varphi q},c_{\varphi u}$,$c_{\varphi d}$,$c^{(1)}_{\varphi l}$
      $c_{\varphi e}$,$c^{(1)}_{lq}$,$c^{(3)}_{lq}$,$c_{lu},c_{ld}$, $c_{qe}$,$c_{eu}$,$c_{ed}$
      $ g _ { V A } ^ { e e } $$c _ { \varphi l } ^ { (1)},c_{\varphi e}, c _ { e e },c _ { l l}$
      $2 \frac{g_{LV}^{\nu_{\mu}\mu,\mathrm{SM}}\delta g_{LV}^{\nu_{\mu}\mu}+g_{LA}^{\nu_{\mu}\mu,\mathrm{SM}}\delta g_{LA}^{\nu_{\mu}\mu}}{\left(g_{LV}^{\nu_{\mu}\mu,\mathrm{SM}} \right)^{2}+\left(g_{LA}^{\nu_{\mu}\mu,\mathrm{SM}} \right)^{2}} $$c_{\varphi l}^{(1)}$,$c_{\varphi e}$,$c_{le}$, $c_{ll}$
      Z-pole$ \Gamma_{Z} $$ c_ { \varphi l } ^ { (1) },c _ {\varphi e },c_ { \varphi q } ^ { (3) },c_ { \varphi q } ^ { (1) },c_{\varphi u},c _ { \varphi d },c_ { \varphi Q} ^ { (+) },c_ { \varphi b }, c_{bB}, c_{bW}$
      $ \sigma_{had} $$ c_ { \varphi l } ^ { (1) },c _ { \varphi e },c_ { \varphi q } ^ { (3) },c_ { \varphi q } ^ { (1) },c_{\varphi u},c _ {\varphi d },c_ { \varphi Q} ^ { (+) },c_ { \varphi b }, c_{bB}, c_{bW}$
      $ R_{e} $$ c_ { \varphi l } ^ { (1) },c _ { \varphi e },c_ {\varphi q } ^ { (3) },c_ { \varphi q } ^ { (1) },c_ { \varphi u },c_{\varphi d},c_ { \varphi Q} ^ { (+) },c _ {\varphi b }, c_{bB}, c_{bW}$
      $ R_{\mu} $$ c_ { \varphi l } ^ { (1) },c _ { \varphi e },c_ { \varphi q } ^ { (3) },c_ { \varphi q } ^ { (1) },c_ { \varphi u },c_{\varphi d},c_ { \varphi Q} ^ { (+) },c _ {\varphi b }, c_{bB}, c_{bW}$
      $ R_{\tau} $$ c_ { \varphi l } ^ { (1) },c _ { \varphi e },c_ { \varphi q } ^ { (3) },c_ { \varphi q } ^ { (1) },c_ { \varphi u },c_{\varphi d},c_ { \varphi Q} ^ { (+) },c _ {\varphi b }, c_{bB}, c_{bW}$
      $ A^{o,e}_{FB} $$c _ { \varphi l } ^ { (1) },c_ { \varphi e }$
      $ A^{o,\mu}_{FB} $$c _ { \varphi l } ^ { (1) },c_ { \varphi e }$
      $ A^{o,\tau}_{FB} $$c _ { \varphi l } ^ { (1) },c_ { \varphi e }$
      Continued on next page

      Table A1.  List of the corresponding coefficients of type I operators that contribute to different observables. Besides these coefficients, there are also four Wilson coefficients that contribute to all observables: $ c^{(3)}_{\varphi l}, c^{'}_{ll}, c_{\varphi D}, c_{\varphi WB} $.

      Table A1-continued from previous page
      ProcessesObservablesWilson coefficient
      $ R_{b} $$ c_ { \varphi q } ^ { (3) },c_ { \varphi q } ^ { (1) },c _ { \varphi u },c _ { \varphi d},c_ { \varphi Q} ^ { (+) },c_{\varphi b}, c_{bB}, c_{bW} $
      $ R_{c} $$ c_ { \varphi q } ^ { (3) },c_ { \varphi q } ^ { (1) },c _ { \varphi u },c _ { \varphi d},c_ { \varphi Q} ^ { (+) },c_{\varphi b}, c_{bB}, c_{bW}$
      $ A^{b}_{FB} $$c_ { \varphi l } ^ { (1) },c_ { \varphi e },c_ {\varphi Q } ^ { (+) },c _ { \varphi b }, c_{bB}, c_{bW}$
      $ A^{c}_{FB} $$ c _ { \varphi l } ^ { (1) },c_ { \varphi e },c_ { \varphi q } ^ { (1) },c _ { \varphi q } ^ { (3) },c _ { \varphi u },$
      $ A_{e} $$c _ { \varphi l } ^ { (1) },c_ { \varphi e }$
      $ A_{\mu} $$c _ { \varphi l } ^ { (1) },c_ { \varphi e }$
      $ A_{\tau} $$c _ { \varphi l } ^ { (1) },c_ { \varphi e }$
      $ A_{b} $$ c_ { \varphi Q } ^ { (+) },c _ { \varphi b }, c_{bB}, c_{bW}$
      $ A_{c} $$ c_ { \varphi q } ^ { (1) },c_ { \varphi q } ^ { (3) },c _ { \varphi u }$
      $ A_{s} $$ c_ { \varphi q } ^ { (1) },c_ { \varphi q } ^ { (3) },c _ { \varphi d }$
      W-pole$ M_{W} $
      $Br(W\to e\nu_{e})$ $Br(W\to\mu\nu_{\mu})$
      $Br(W \to\tau\nu_{\tau})$$\Gamma_{W}$
      $c _ { \varphi q } ^ { (3) }$
      $ee\to q\bar{q}$$ \sigma_{eeqq} $$ c _ { \varphi l } ^ { (1) }$,$c _ { \varphi e }$,$ c _ { \varphi q } ^ { (3) }$,$c _ {\varphi q } ^ { (1) }$,$c _ { \varphi u }$,$c _ { \varphi d }$,$c _ { \varphi Q} ^ { (+) }$,$c_ { \varphi b }$,$c _ { l q } ^ { (1) }$,$c _ { l q } ^ { (3) }$,
      $c_{ l u }$,$c _ { l d }$,$c_ { l Q } ^ { (+) }$,$c _ { l b}$,$c _{ q e }$,$c_{eu}$,$c_ { e d }$,$c _ { eQ}$,$c _ { e b}, c_{bB}, c_{bW}$
      $ee\to b\bar{b}$$ \sigma_{eebb} $$ c _ { \varphi l } ^ { (1) },c _ { \varphi e },c _ { \varphi Q } ^ { (+) },c _ { \varphi b },c_ { l Q}^ { (+) },c _ { l b},c _ { eQ},c_ { e b}, c_{bB}, c_{bW}$
      $ ee\to c\bar{c} $$ \sigma_{eecc} $$c _ { \varphi l } ^ { (1)}$,$c _ { \varphi e }$,$c _ { \varphi q } ^ { (1) }$,$c _ { \varphi q } ^ { (3) }$,$c_ { \varphi u }$,$c_{ l q } ^ { (1) }$,$c_ { l q } ^ { (3) }$,$c _ { l u }$,$c_ { q e }$,$c_ { e u }$
      $ ee\to \mu^{+}\mu^{-} $$ \sigma_{ee\mu\mu} $$c_ { \varphi l } ^ {(1) }$,$c_ { \varphi e }$,$c_ { ll}$,$c _ { l e }$,$c _ { e e } $
      $ ee\to\tau^{+}\tau^{-} $$ \sigma_{ee\tau\tau} $$c_ { \varphi l } ^ { (1) }$,$c_ { \varphi e }$,$c_ { ll}$,$c _ { l e }$,$c _ { e e } $
      $ ee\to e^{+}e^{-} $$ \sigma_{eeee} $$c_ { \varphi l } ^ { (1)}$,$c_ { \varphi e }$,$c_ { ll}$,$c _ { l e }$,$c _ { e e } $
      $ ee\to W^{+}W^{-} $$ \sigma_{eeww} $$c _ { \varphi l } ^ { (1)}$,$c_ { \varphi e }$,$c _ { W }$
    • A.   Z, W pole observables and low energy couplings

    • The theoretical predictions of Z pole and W pole observables are taken from Table 10.5 of [42] and Table 2 of [11]. The theoretical errors on the Z and W pole observables and low energy couplings, except for $ m_{W} $, are ignored (the theoretical prediction error om $ m_{W} $ is 0.004 GeV [42, 68]).

    • B.   Fermion-pair production

    • The main Feynman diagram of $ e^+ e^- \rightarrow f \bar{f} $ contains an s-channel photon and Z annihilation (Fig. 3(a)). However, there is an additional t-channel photon and Z exchange in the $ e^+ e^- \rightarrow e^+ e^- $ process (Fig. 3(b)), which have a small contribution at the Z-pole. We make the assumption of zero electron and muon mass in our model. As a result, there is only a negligible difference between the contribution of new operators to the observables of the $ e^+ e^- \rightarrow e^+ e^- $ and $ e^+ e^- \rightarrow \mu^+ \mu^- $ processes, and we take the result of $ e^+ e^- \rightarrow \mu^+ \mu^- $ as that of $ e^+ e^- \rightarrow e^+ e^- $ at the Z-pole. At other energies, both t-channel and s-channel production of $ e^+ e^- \rightarrow e^+ e^- $ are considered.

      Figure 3.  Feynman diagrams for the $ e^+ e^- \to $$ f \bar{f} $ process at the Born level.

      For quark pair production, the decay of quarks is not considered. Additional observables in the low energy region for b and c quark production are included, which helps us obtain a better constraint on the 3rd-generation-quark operators.

      For fermion-pair production, especially at energies above Z resonance, the QED radiative corrections are very large. This is caused by the initial-state radiation of photons, which lowers the center-of-mass energy, $ \sqrt{s} $, to a value of $ \sqrt{s'} $. In the LEP experiment, the measurement is performed using events with a small amount of initial state radiation, that is, large $ \sqrt{s/s'} $. Therefore, in Monte Carlo integration, initial-state radiation is not considered.

      ● In the $ e^{+}e^{-}\to f\bar{f} $($ f=\mu, \tau, q $) processes, the experimental values of cross-sections (in pb) and forward-backward asymmetries are reported in Table 3.4 of [65]. The theoretical prediction uncertainties on $ \sigma(q\bar{q}), \sigma(\mu^{+}\mu^{-}) $, and $ \sigma(\tau^{+}\tau^{-}) $ and leptonic forward-backward asymmetries are $ 0.26 $%, $ 0.4 $%, $ 0.4 $%, and $ 0.4 $%, respectively [65, 69]. Their theoretical prediction uncertainties can be neglected because their experimental uncertainties are at least four times larger.

      ● The experimental and theoretical values of the $ e^{+}e^{-}\to\mu^{+}\mu^{-} $and $ e^{+}e^{-}\to\tau^{+}\tau^{-} $ differential cross-sections are reported in Tables 3.8 and 3.9 of [65]. In the processes $ e^{+}e^{-}\to\mu^{+}\mu^{-} $ and $ e^{+}e^{-}\to\tau^{+}\tau^{-} $, we assume that the theoretical uncertainty on the differential distribution is 0.4%, which can be neglected compared with experimental values.

      ● The experimental and theoretical values of the $ e^{+}e^{-}\to e^{+}e^{-} $ differential cross-section are reported in Tables 3.11 and 3.12 of [65]. The theoretical uncertainty on large-angle Bhabha scattering for $ \sqrt{s} $ ranging from 189 to 207 GeV is approximately 0.5% [69, 70]. The theoretical uncertainty on the differential distribution is also taken as 0.5%.

      ● For the process $ e^{+}e^{-}\to q\bar{q} $ at $ \sqrt{s}=58 $ GeV, Table 8 in [71] show the experimental values for $ R_{b} $ and $ R_{c} $, which are the bottom and charm quark pair production cross-section ratios to the total hadronic cross-section. In [72], their theoretical uncertainties are estimated as 1% and can be ignored.

      ● For the process $ e^{+}e^{-}\to q\bar{q} $ at $ \sqrt{s} $ ranging from 189 to 207 GeV, Tables 8.9 and 8.10 in [66] show the experimental values for $ R_{b} $ and $ R_{c} $ without theoretical uncertainties. We assume the theoretical uncertainties are negligible.

    • C.   Boson-pair production

    • Boson-pair production is also very important in LEP experiments. In our study, the $ e^+ e^- \rightarrow W^+ W^- $ process is included. There are three main Feynman diagrams for this process, shown in Fig. 4, which are known as CC03 (Charged Current).

      Figure 4.  Main Feynman diagrams for $ e^+ e^- \rightarrow W^+ W^- $ at the Born level.

      Figure 13.  (color online) Comparison of the current precision of EW observables and the projections for CEPC and FCC-ee.

      The W boson will decay into a quark-antiquark pair or lepton-neutrino pair. The BR of hadronic decay is approximately 67%, and that of leptonic decay is approximately 33%. $ WW $ events are classified into fully hadronic, semi-leptonic, and fully leptonic events. In the Monte Carlo simulation of the total cross-section of $ WW $ production, we ignore the decay of W bosons and only produce events that take the W pair as the final state for convenience. The angular distribution of $ WW $ events ($ {\rm d}\sigma_{WW}/{\rm dcos}\theta_{W^-} $) is also included in this study, where $ \theta_{W^-} $ is the polar angle between $ W^- $ and the $ e^- $ beam direction. In the differential cross-section measurement, the decay of $ WW $ is selected as $ q\bar{q}e\nu $ and $ q\bar{q}\mu\nu $ because these two types of leptonic decays provide the charge tag and a clean background compared with the τ channel. There is a constraint on the experiment that the charged lepton should be $ 20^\circ $ from the beam, $ |\theta_l|>20^\circ $. This angular cut brings about a 93% cut efficiency in the experiment, [65] which matches our Monte Carlo simulation with the SM model. This angular requirement corresponds to the experimental acceptance of the four LEP experiments and also greatly reduces the difference between the full 4f cross-section and the CC03 cross-section by reducing the contribution of t-channel diagrams in the $ q\bar{q}e\nu $ final state. (As stated in the LEP result, the difference of the 4f and CC03 cross-sections decreased from 24.0% to approximately 3.5% with the angular cut.) However, its effect is no longer the same for the contribution from new physics. Based on our simulation, the angular cut efficiency of some operators is shown in Fig. 5. The cut efficiency of operators that contribute to the differential cross-section at different energies can be found in the appendix. As a result, the decay of W cannot be ignored in the simulation of W angular distribution. Taking the four fermions as the final state, the process $ e^+ e^- \rightarrow W^+ W^- \rightarrow q \bar{q} l \nu $ would be a subset of the process $ e^+ e^- \rightarrow q \bar{q} l \nu $. If we simulate the process $ e^+ e^- \rightarrow W^+ W^- \rightarrow q \bar{q} l \nu $ and allow the W boson to be off-shell, the gauge invariance would be broken. However, simulating the process $ e^+ e^- \rightarrow q \bar{q} l \nu $ would bring about an additional Feynman diagram. Finally, we decided to simulate the process $e^+ e^- \rightarrow W^+ W^- \rightarrow q \bar{q} l \nu$ using MadSpin [73] to place W on shell decay to guarantee gauge invariance.

      Figure 5.  (color online) Impact of the 20 degree angular cut on the operators $c_{\varphi D} $ and $ c_{\varphi e}$ at 183 GeV. The red line shows the differential cross-section of the operators before the angular cut, and the blue line shows the differential cross-section after the cut. The overall cut efficiency is given by the total cross-section after the cut divided by that before the cut.

      In [74, 75], the theoretical uncertainties on the $ W^{-} $ angular distribution are approximately $ 0.5 $% when $ \sqrt{s} $ ranges from 180 to 210 GeV. Because the experimental precision of $ W^{-} $ differential angular cross-section are at least seven times larger than the theoretical ones, we can neglect the theoretical prediction uncertainties on $ W^{-} $ differential angular cross-section. Although the theoretical uncertainties on the total cross-section of W pair production are also approximately $ 0.5 $% for $ \sqrt{s} $ ranging from 180 to 210 GeV, they are non-negligible because the experimental uncertainties are comparable with the theoretical values.

    IV.   FIT RESULT

      A.   Methodology

    • The SMEFT Wilson coefficients are estimated using the least squares method under the assumption of Gaussian errors. The $ \chi^2 $ function is constructed as

      $ \begin{equation} \chi^2=(\vec{y}-A\vec{c})^TV^{-1}(\vec{y}-A\vec{c}) \,, \end{equation} $

      (14)

      where $ \vec{y} $ is the vector of the difference between the experimental result and SM prediction of the observables, A is the contribution matrix ($ A_{ij} $ is the i-th operator contribution to j-th observables), $ \vec{c} $ is the vector of the Wilson coefficient, and V is the covariance matrix of the observables.

      By letting $ \nabla\chi^2=0 $, we can obtain the least square estimator $ \hat{\vec{c}}\; $ that minimizes $ \chi^2 $ and its covariance matrix U,

      $ \begin{equation} \hat{\vec{c}}=(A^TV^{-1}A)^{-1}A^TV^{-1}\vec{y}=B\vec{y}, \end{equation} $

      (15)

      $ \begin{equation} U=BVB^T=(A^TV^{-1}A)^{-1}. \end{equation} $

      (16)

      The element of the covariance matrix, $ U_{ij} $, represents the covariance of the estimators $ \hat{c_i} $ and $ \hat{c_j} $. The one-sigma bound of $ c_i $ can be obtained from the diagonal element of this matrix.

    • B.   Global analysis

    • We apply the aforementioned fit strategy to obtain the constraints of the EFT operator coefficients. Two types of bounds are shown. The first is the marginalized bound, derived allowing all operator coefficients to float. When conducting this type of fit, the correlations between the coefficients also contain useful information. The second is the individual bound, which is obtained by considering only one operator coefficient at a time while fixing all others to zero. In the following results, all the bounds are given in the 68% confidence level (CL). The new physics scale Λ of the Wilson coefficients $ \dfrac{c_i}{\Lambda^2} $ is fixed to be 1 TeV throughout this section.

      We consider a total of 33 operators listed in Section II, which are constrained by the measurements in Sec. III. As mentioned in Sec. II, we trade $ \dfrac{c_{\varphi D}}{\Lambda^{2}}Q_{\varphi D} $ and $ \dfrac{c_{\varphi WB}}{\Lambda^{2}}Q_{\varphi WB} $ in the Warsaw basis for $ \dfrac{c_{D\varphi B}}{\Lambda^{2}}{\rm i}D^{\mu}\varphi^{\dagger}D^{\nu}\varphi B_{\mu\nu} $ and $ \dfrac{c_{D\varphi W}}{\Lambda^{2}}{\rm i}D^{\mu}\varphi^{\dagger}\sigma_{a}\varphi W^{a}_{\mu\nu} $ to disentangle several large correlations, whereas the other operators in the Warsaw basis remain unchanged. Our main results are presented in Fig. 6 in terms of the one-sigma (corresponding to a 68% confidence level (CL)) bounds of the Wilson coefficients. To understand the impact of the 3rd generation quark operators, we also consider two additional scenarios for comparison: The first considers only tree-level contributions from dimension-6 operators, whereas the second also includes the loop contributions of the bottom-quark operators (that is, excluding $ Q^{(-)}_{\varphi Q} $, $ Q_{\varphi t} $, $ Q_{tW} $, and $ Q_{tB} $). The correlations between the coefficients are shown in Fig. 7 for the tree-level-only scenario and Fig. 8 for the full scenario. The numerical results can be found in Table A2 of the appendix, along with the Fisher information matrix for the full scenario.

      Figure 6.  (color online) 68% CL bounds of operator coefficients in the modified Warsaw basis for three different scenarios. (blue) Tree-level contribution only. (green) Tree-level contribution and loop contributions from the bottom-quark operators. (red) Tree-level contribution and the loop contributions of all third-generation-quark operators listed in Sec. II. The dark red bar is the individual bound in this scenario. See Table A2 in the appendix for the best-fitted central values. The value of Λ is set to 1 TeV for all coefficients.

      Figure 7.  (color online) Correlation matrix between operator Wilson coefficients in the modified Warsaw basis. Only the tree-level contribution of these operators is considered.

      Figure 8.  (color online) Correlation matrix between operator Wilson coefficients in the modified Warsaw basis. The top and bottom loop contributions are considered in addition to the tree-level contribution.

      OperatorsMarginalized,
      tree-only
      Marginalized,
      b loop considered
      Marginalized,
      t/b loop considered
      Individual,
      t/b loop considered
      $ c_{tB} $$ 30 \pm 35 $$ -0.13 \pm 0.18 $
      $ c_{le} $$ 0.0019 \pm 0.018 $$ 0.0014 \pm 0.018 $$ 0.036 \pm 0.036 $$ -0.0015 \pm 0.016 $
      $ c_{ee} $$ -0.019 \pm 0.014 $$ -0.019 \pm 0.014 $$ -0.0082 \pm 0.018 $$ -0.012 \pm 0.012 $
      $ c_{\varphi l}^{(3)} $$ -0.023 \pm 0.023 $$ -0.027 \pm 0.018 $$ -0.42 \pm 2 $$ -0.0036 \pm 0.0033 $
      $ c_{\varphi q}^{(3)} $$ -0.079 \pm 0.07 $$ -0.11 \pm 0.074 $$ -0.47 \pm 2 $$ 0.0025 \pm 0.0076 $
      $ c_{W} $$ -7.9 \pm 5.1 $$ -1.8 \pm 3.7 $$ -4.5 \pm 4.5 $$ 0.41 \pm 0.37 $
      $ c_{D\varphi B} $$ -13 \pm 10 $$ -7.9 \pm 11 $$ -12 \pm 17 $$ -0.2 \pm 1.9 $
      $ c_{D\varphi W} $$ 35 \pm 21 $$ 12 \pm 17 $$ 23 \pm 20 $$ 0.82 \pm 1.1 $
      $ c_{\varphi l} $$ 0.024 \pm 0.013 $$ 0.021 \pm 0.013 $$ -0.26 \pm 0.67 $$ 0.0014 \pm 0.0042 $
      $ c_{\varphi e} $$ 0.0051 \pm 0.014 $$ -2.5e-05 \pm 0.013 $$ -0.57 \pm 1.3 $$ -0.0018 \pm 0.0054 $
      $ c_{\varphi q} $$ -0.042 \pm 0.13 $$ -0.053 \pm 0.14 $$ 0.085 \pm 0.28 $$ -0.015 \pm 0.014 $
      $ c_{\varphi u} $$ 0.086 \pm 0.17 $$ 0.11 \pm 0.17 $$ 0.49 \pm 0.91 $$ -0.011 \pm 0.021 $
      $ c'_{ll} $$ -0.043 \pm 0.03 $$ -0.055 \pm 0.029 $$ -0.069 \pm 0.094 $$ 0.0035 \pm 0.0051 $
      $ c_{ll} $$ 0.054 \pm 0.033 $$ 0.066 \pm 0.031 $$ 0.1 \pm 0.093 $$ 0.0037 \pm 0.011 $
      $ c_{tW} $$ 4.2 \pm 27 $$ -0.3 \pm 0.22 $
      $ c_{\varphi b} $$ -0.34 \pm 0.62 $$ -0.27 \pm 0.63 $$ -12 \pm 14 $$ -0.17 \pm 0.11 $
      $ c_{bB} $$ 10 \pm 15 $$ 9.2 \pm 14 $$ 12 \pm 18 $$ 0.25 \pm 0.23 $
      $ c_{bW} $$ 16 \pm 27 $$ 9.7 \pm 25 $$ 17 \pm 32 $$ 0.1 \pm 0.13 $
      $ c_{\varphi Q}^{(+)} $$ -3.5 \pm 5.4 $$ -2.3 \pm 5 $$ -6.1 \pm 5.5 $$ 0.017 \pm 0.02 $
      $ c_{lQ}^{(+)} $$ -1.5 \pm 2.1 $$ -2 \pm 2.4 $$ -4.7 \pm 5.5 $$ -0.048 \pm 0.042 $
      $ c_{lb} $$ 1.2 \pm 2.3 $$ 1.4 \pm 2.5 $$ -1.2 \pm 5.9 $$ 0.0074 \pm 0.23 $
      $ c_{eQ} $$ 0.64 \pm 1.8 $$ 0.56 \pm 1.8 $$ -1.5 \pm 4.8 $$ 0.035 \pm 0.11 $
      $ c_{eb} $$ 3.6 \pm 5 $$ 5.1 \pm 5.6 $$ 11 \pm 12 $$ -0.1 \pm 0.1 $
      $ c_{\varphi Q}^{(-)} $$ 31 \pm 2.7e+02 $$ -0.043 \pm 0.14 $
      $ c_{\varphi t} $$ 3.8 \pm 1.9e+02 $$ 0.011 \pm 0.12 $
      $ c_{\varphi d} $$ -0.7 \pm 0.8 $$ -0.92 \pm 0.86 $$ -0.69 \pm 0.99 $$ -0.038 \pm 0.028 $
      $ c_{lq}^{(1)} $$ 2.3 \pm 1.2 $$ 2.7 \pm 1.2 $$ 2.6 \pm 1.2 $$ 0.014 \pm 0.014 $
      $ c_{lq}^{(3)} $$ 0.74 \pm 0.43 $$ 0.88 \pm 0.43 $$ 0.84 \pm 0.44 $$ 0.029 \pm 0.018 $
      $ c_{lu} $$ -1.1 \pm 0.83 $$ -1.3 \pm 0.82 $$ -1.3 \pm 0.83 $$ 0.017 \pm 0.026 $
      $ c_{ld} $$ -4.8 \pm 2.4 $$ -5.8 \pm 2.4 $$ -5.2 \pm 2.4 $$ 0.042 \pm 0.03 $
      $ c_{eq} $$ 2.6 \pm 1.6 $$ 3 \pm 1.6 $$ 2.8 \pm 1.6 $$ -0.019 \pm 0.014 $
      $ c_{eu} $$ -2.2 \pm 1.1 $$ -2.6 \pm 1.1 $$ -2.4 \pm 1.1 $$ -0.039 \pm 0.023 $
      $ c_{ed} $$ -3.6 \pm 2.1 $$ -4.2 \pm 2 $$ -4.1 \pm 2.1 $$ -0.022 \pm 0.029 $

      Table A2.  Numerical fit results of the operator coefficients in the modified Warsaw basis for three different scenarios. The error bound is given in the 68% CL.

      Overall, our results demonstrate the relevance of 3rd generation quarks in electroweak measurements. The coefficients of the four top operators $ Q^{(-)}_{\varphi Q} $, $ Q_{\varphi t} $, $ Q_{tW} $, and $ Q_{tB} $ are constrained in the range of $ 0.1 $ to $ 1 $ for the individual fit, despite the fact that their contributions only enter at the one-loop order. Assuming an order-one coupling, this corresponds to a new physics scale of a few TeVs. After comparing these results with those in Ref. [25], which uses LHC data to probe these four top operators, it is found that the electroweak data have greater constraining power on $ Q^{(-)}_{\varphi Q} $ and $ Q_{\varphi t} $, whose advantage is estimated to be of the order of 10. Conversely, for $ Q_{tW} $ and $ Q_{tB} $, the LHC data are more powerful. Their 95% CL individual bounds are listed in Table 3 for comparison.

      $ c_{\varphi t} $$ c_{\varphi Q}^{(-)} $$ c_{tW} $$ c_{tB} $
      Electroweak0.2860.3360.8220.592
      LHC data2.2751.220.060.145

      Table 3.  95% CL individual bounds of top quark operator coefficients using different sets of observables.The electroweak results are from this study, whereas the LHC data are the results in [26].

      In a global fit, however, the marginalized bounds of the top-quark operator coefficients become significantly looser, and some of them possibly exceed the range of EFT validity. This is not surprising because the introduction of additional degrees of freedom tends to bring additional flat directions, making the fit difficult to converge. This can also be verified by comparing the correlations in Fig. 7 and Fig. 8; the increase in correlation is visible in the latter. The inclusion of top operators also significantly degrades the reaches of several other operator coefficients. Among them, the most notable are the leptonic operators $ c_{l e} $, $ c_{e e} $, $ c^{(3)}_{\varphi l} $, $ c_{\varphi l} $, and $ c_{\varphi e} $ (as well as $ c^{(3)}_{\varphi q} $), which were previously well constrained and are thus most sensitive to additional degrees of freedom. Note that $ c_{l e} $and $ c_{e e} $ contribute directly to 4-lepton processes, and their contributions can be distinguished from those with a gauge boson propagator by measuring the processes at several different energies. On the other hand, when only the loop contributions of the bottom-quark operators are included (green bar), the overall reach of the global fit is not significantly degraded compared with that of the tree-level fit, and the reach on some of the operator coefficients are even improved. This is because these operators already contribute at the tree level; hence, the number of operators is not increased, whereas their dependence on the observables changes with the inclusion of loop effects.

      It should be noted that, as pointed out in Section II, the effects of the two bottom-quark dipole operators $ Q_{bW} $ and $ Q_{bB} $ are non-negligible and should be included in the fit. Even at the tree level, this introduces large flat directions with the operators $ Q^{(+)}_{\varphi Q} $ and $ Q_{\varphi b} $because the four operators are mainly constrained by only two observables, $ R_b $ and $ A_b \,(A^b_{\rm FB}) $. Among them, we observe particular strong correlations among $ Q_{bW} $, $ Q_{bB} $, and $ Q^{(+)}_{\varphi Q} $ ($ >0.98 $), whereas their correlations with $ Q_{\varphi b} $ are relatively smaller ($ <0.9 $). This is because the former three operators mainly contribute to $ R_b $, while $ Q_{\varphi b} $ modifies the $ Zb_R\bar{b}_R $ coupling and is more sensitive to asymmetry (see, for example, Ref. [76]).

      Additional results are also provided in Table A2 in the appendix, which shows both the marginalized bound with central values and individual bounds. Figure A1 shows the impact of different sets of measurements on the Wilson coefficients, which is given by Fisher information.

      Figure A1.  (color online) Impact of different sets of measurements on the Wilson coefficients in the Warsaw basis, measured using Fisher information. The larger the number is, the bigger impact of the set have on the constraints of Wilson coefficients.

    • C.   ${\boldsymbol {S\,T \, t\,b} }$ scenario

    • In addition to the general global-fitting framework, we also consider a special case denoted here as the "$ S\,T \, t\,b $" scenario, where we assume that all new physics effects, apart from those in Eq. (2), can be parameterized by the two oblique parameters S and T [47]. This scenario is motivated by a large class of models with top/bottom partners, which generally mixes with 3rd generation quarks and can also contribute to oblique parameters at the one-loop order. For convenience, we work with the modified parameters $ \hat{S} $ and $ \hat{T} $ [77], which are related to the original parameters via

      $ \begin{equation} \hat{S} = \frac{\alpha}{4 s^2_w} S \,, \quad \hat{T} = \alpha T \,. \end{equation} $

      (17)

      $ \hat{S} $and $ \hat{T} $ are defined based on [78]. The measurements considered are still the same as in Sec. IV.B.

      For comparison, we first consider the two-parameter fit of $ \hat{S} $ and $ \hat{T} $. Their 68% CL preferred region and correlation are shown in Fig. 9. The numerical results of the individual and marginalized bound are

      Figure 9.  (color online) 68% CL preferred region of $ \hat{S} $ and $ \hat{T} $.

      $ \begin{equation} \begin{matrix} \;\;\;\;\;\;\hat{S}=&(-2.8\pm2.5)\times10^{-4} \\ \;\;\;\;\;\;\hat{T}=&(2.7\pm1.9)\times10^{-4} \end{matrix} \quad \rm{(individual bound)}\,, \end{equation} $

      (18)

      and

      $ \begin{equation} \begin{matrix} \;\;\;\;\;\; \hat{S}=&(4.0\pm6.9)\times10^{-4} \\ \;\;\;\;\;\; \hat{T}=&(5.5\pm5.1)\times10^{-4} \end{matrix} \quad \rm{(marginalized bound)}\,, \end{equation} $

      (19)

      with the correlation

      $ \begin{equation} corr(\hat{S},\hat{T})=0.93. \end{equation} $

      (20)

      These results are comparable with those in Ref. [79], which uses the same set of measurements as we do.

      With the inclusion of the 3rd-generation quark operators, we perform a 10-parameter fit, and the results are presented in Fig. 10. It should be noted that part of the loop contributions of the 3rd-generation quark operators is universal and can be absorbed in the original definition of the S and T parameters. Here, by $ \hat{S}, \hat{T} $, we simply denote the contributions of the corresponding tree-level operators. From Fig. 10, we observe a significantly better reach on the top quark operators compared with the general case in Sec. IV.B. However, the constraining power to the $ \hat{S} $ and $ \hat{T} $ parameters with the electroweak data is significantly reduced from the 2-parameter case, as expected.

      Figure 10.  (color online) Fit results ((left) 68% CL marginalized bound, (right) correlation matrix) of the $ \hat{S} $ and $ \hat{T} $ parameters and the eight operator coefficients in Eq. (2).

      One could also consider a more general set of universal corrections [78] (in addtion to the 3rd generation quark operators), which also requires additional measurements. For instance, the W and Y parameters are strongly constrained by the high energy Drell-Yan measurement at the LHC [80]. A detailed analysis in this direction is left for future studies.

    • D.   Impact of the new W-mass measurement at the CDF

    • The CDF experiment recently announced their new W mass measurement [40], which deviated from the SM prediction by seven standard deviations. This measurement will greatly affect the fit results. Therefore, in this section, the fit results with this new measured W mass are shown in contrast with those using the previously measured value. Fig. 11 shows the individual bound changes with the new measurement.

      Figure 11.  (color online) Individual bounds of operator coefficients in the modified Warsaw basis with the previous (blue) and new (red) W mass measurement. Here, both the tree-level contribution and the loop contributions of all third-generation-quark operators are considered.

      With the new CDF W-mass measurement, the results of the $ \hat{S} $-$ \hat{T} $ two-parameter fit are

      $ \begin{matrix} \;\;\;\;\;\;\hat{S}=&(-13.0\pm2.3)\times10^{-4} \\ \;\;\;\;\;\; \hat{T}=&(10.9\pm1.6)\times10^{-4} \end{matrix} \quad \rm{(individual bound)}\,, $

      (21)

      and

      $ \begin{equation} \begin{matrix} \;\;\;\;\;\; \hat{S}=&(13.2\pm6.8)\times10^{-4} \\ \;\;\;\;\;\; \hat{T}=&(19.7\pm4.8)\times10^{-4} \end{matrix} \quad \rm{(marginalized bound)}\,. \end{equation} $

      (22)

      For the "$ S\,T \, t\,b $" scenario, Fig. 12 shows the marginalized bound with the new W mass measurement. The impact of the shifted W mass value is greatly absorbed by the 3rd-generation-quark operators.

      Figure 12.  (color online) 68% CL marginalized bound of operator coefficients in the "$ S\,T \, t\,b $" scenario with the new W mass measurement (red) and previous value (blue).

      There are several papers [8190] that also study the impact of the new $ m_W $ measurement on new physics scenarios, especially the SMEFT fit and oblique parameters. These results are generally in good agreement with our results.

    V.   FUTURE COLLIDER
    • Recently, there have been several proposals on the construction of future colliders that can reach higher luminosities and energies and will thus significantly increase the precision of electroweak measurements. It is interesting to see how the improvement of measurement affects the constraints on these EFT operators. As a result, the projections of electroweak observables from FCC-ee and CEPC are collected to conduct a chi-square fit. For the other two future collider proposals, the ILC and CLIC, because their Giga-Z proposals are still under discussion and currently there is little information on their expected performance for electroweak observables, they are not included in this study.

      In Table 4, we show a comparison of the current precision of electroweak observables and the projections for the CEPC and FCC. Because both experiments are designed to be tera-Z factories, the projections of observables at the Z-pole are taken as a common value. For the CEPC, the projections of luminosity at the WW threshold and 240 GeV are taken from Ref. [91]. For the FCC-ee, the predictions are taken from Ref. [92]. The projections of the constraints on the EFT operators involved in this study are shown in Fig. 13. In general, the future collider will greatly improve the constraining power on these EFT operators; most of the coefficients will obtain at least 10 times better constraints, both on individual and marginalized bounds. Though 3rd-generation quark operators can only be probed at loop level in current projections, they still have significantly better constraints. For several operators, such as $ Q_W $ and $ Q_{eq} $, their constraints experience relatively small improvement because their constraints mainly originate from the $ e^+e^-\to W^+W^- $ process or low energy coupling measurements, which lack the projection of their precision for future colliders. However, it is believed that they can also benefit from the high precision measurements of future colliders. Finally, a high-energy lepton collider running at and above the top threshold will be able to measure the $ e^+e^-\to t\bar{t} $ process to an unprecedented precision, providing the ultimate probe for 3rd-generation quark operators [15, 20, 25]. The interplay between the top, Higgs, and EW measurements in probing 3rd-generation quark operators is an important topic for future collider studies; however, this is beyond the scope of our current study.

      ObservablesSMCurrent precisionCEPCFCC-ee
      $ \Gamma_Z $2.4942 G${\rm e}{\rm V}$2.3 M${\rm e}{\rm V}$0.1 M${\rm e}{\rm V}$0.1 M${\rm e}{\rm V}$
      $\sigma_{\rm had}$41.481 nb0.037 nb0.004 nb0.004 nb
      $ R_e $20.7370.050.0010.001
      $ R_\mu $20.7370.0330.0010.001
      $ R_\tau $20.7820.0450.0010.001
      $ R_b $0.215820.000660.000060.00006
      $ R_c $0.172210.0030.000260.00026
      $A_{\rm FB}^b$0.1030.00030.00030.0003
      $ A_\tau $0.14720.002160.00020.0002
      $ M_W $80.379 G${\rm e}{\rm V}$12 M${\rm e}{\rm V}$1 M${\rm e}{\rm V}$0.5M${\rm e}{\rm V}$
      $ \Gamma_W $2.085 G${\rm e}{\rm V}$42 M${\rm e}{\rm V}$2.8 M${\rm e}{\rm V}$1.2 M${\rm e}{\rm V}$
      $\sigma_{ee\to ff}^{160\;{\rm GeV} }$$ \sqrt{\dfrac{\sigma_{\rm SM}}{2.6ab^{-1}}} $$ \sqrt{\dfrac{\sigma_{\rm SM}}{12ab^{-1}}} $
      $\sigma_{ee\to ff}^{240\;{\rm GeV} }$$ \sqrt{\dfrac{\sigma_{\rm SM}}{5.6ab^{-1}}} $$ \sqrt{\dfrac{\sigma_{\rm SM}}{5ab^{-1}}} $
      $\sigma_{ee\to \mu\mu/\tau\tau, \rm differential}^{240\;{\rm GeV} }$$ \sqrt{\dfrac{\sigma_{\rm SM}}{5.6ab^{-1}}} $$ \sqrt{\dfrac{\sigma_{\rm SM}}{5ab^{-1}}} $

      Table 4.  Comparison of the projections of electroweak observables for the CEPC [91] and FCC-ee [92].

    VI.   CONCLUSION
    • Precision measurements of electroweak processes offer important probes for physics beyond the SM. Many SMEFT analyses of electroweak measurements focus on the tree-level contributions of dimension-6 operators. However, given the outstanding precision of these measurements (especially at future lepton colliders), they could be sensitive to many important loop contributions of new physics, which are not captured by simple tree-level treatment. In this paper, we attempt to extend the tree-level framework by including the one-loop contributions of operators involving third generation quarks and study their impact in a global analysis of electroweak measurements. This is motivated by many new physics scenarios in which 3rd generation quarks play a special role. We include the measurements of the $ e^+e^-\to f\bar{f} $ processes around the Z-pole and at several other energies, the measurements of $ e^+e^-\to WW $ at LEP2, and a collection of low energy scattering processes. We find that the 3rd-generation quark operators, and especially top quark ones, offer significant contributions to electroweak processes. In individual fits, where only one operator coefficient is considered at a time, we obtain competitive reaches on these operators, which are all constrained to be at least approximately 1 TeV, with the order-one coupling assumption. On the other hand, in a global framework, where all tree-level operator contributions are also included, the reaches on the operator coefficients become significantly worse because it is difficult to separate the loop effects of the 3rd-generation quark operators from other tree-level effects. However, one should not conclude that the considerations of these loop contributions are meaningless in the SMEFT framework. One important goal of the SMEFT is to provide a bridge between the experimental constraints and the parameters in the UV model, and the likelihood from the SMEFT global analysis could be directly translated to the bounds on the UV model, even if the fit in the SMEFT is not close. In a particular UV model, we usually expect a significantly smaller parameter space, and a global fit with loop effects becomes much more feasible. As a demonstration, we perform a fit in a more constrained scenario, where the only tree level contributions are parameterized by the two oblique parameters S and T. In this case, better constraints are obtained. We also apply our analysis to the electroweak measurements at future lepton colliders. Another important finding of our study is that the tree-level contributions of the bottom dipole operators to the electroweak processes are non-negligible, and their effects are generally difficult to separate from the modifications of $ Zb\bar{b} $ couplings.

      Our study is one of the many first steps toward a more complete loop-level SMEFT global analysis, for which many improvements are still required. Throughout our study, the loop effects of 4-fermion operators involving 3rd generation quarks are not considered. These contributions could be comparable to those considered in our study and should, in principle, also be included. However, they would introduce more degrees of freedom, and additional measurements would be required to discriminate their effects. Similarly, it is also desirable to remove the flavor assumptions imposed on our study, which also significantly increases the size of the parameter space and requires additional measurements [93]. It is also important to study the complementarity between the direct probes of the top quark operators, either at hadron colliders or a future lepton collider with higher center-of-mass energies, and the indirect probes studied here. Third generation quark operators also enter Higgs processes, and a combined Higgs and electroweak analysis is particularly relevant for future lepton colliders in this framework. Previously, an optimal-observable analysis of $ e^+e^- \to WW $ with information on W decay angles was shown to be useful in probing the corresponding tree-level operators [19, 94], which could be extended to include loop effects. However, this requires additional effort in calculating the one-loop contributions to the full differential cross section. We leave these many possible extensions of our current analysis to future studies.

    ACKNOWLEDGMENTS
    • Cen Zhang, who was our collaborator, friend, and Yiming Liu's advisor, unexpectedly passed away while leading the work of this paper. Jiayin Gu joined at a later stage to help complete this paper.

    APPENDIX A

      A.1.   List of operator contributions to observables

      A.2.   Additional fit results

    • The Fisher information of set i on coefficient $ c_j $ is calculated using

      $ f_i=\frac{\dfrac{\partial^2 \chi^2_i}{\partial c_j^2}}{\dfrac{\partial^2 \chi^2_{all}}{\partial c_j^2}}, \tag{A1} $

      where $ \chi^2_i $ indicates the chi-squared calculated with only the data of set i, and $ \chi^2_{\rm all} $ indicates the chi-squared calculated with all sets of measurements.

    • A.3.   Definition of renormalized parameters

    • $ m^{2}_{W^{*}}(q^{2})=(1-Z_{W})q^{2}+Z_{W}\bigg(m^{2} _{WO}+\Pi_{WW}(q^{2})\bigg), \tag{A2}$

      $ m^{2}_{Z^{*}}(q^{2})=(1-Z_{Z})q^{2}+Z_{Z}\bigg(m^{2} _{ZO}+\Pi_{ZZ}(q^{2})\bigg), \tag{A3}$

      $ Z_{W}=1+\frac{\mathrm{d}\Pi_{WW}(q^{2})}{\mathrm{d}q^{2}}\Big|_{q^{2}=m^{2}_{W}}, \tag{A4} $

      $ Z_{Z}=1+\frac{\mathrm{d}\Pi_{ZZ}(q^{2})}{\mathrm{d}q^{2}}|_{q^{2}=m^{2}_{Z}}, \tag{A5} $

      $ Z_{W^{*}}(q^{2})=1+\frac{\mathrm{d}\Pi_{WW}(q^{2})}{\mathrm{d}q^{2}}|_{q^{2}=m^{2}_{W}}-\Pi^{'}_{\gamma\gamma}(q^{2})-\frac{c_{W}}{s_{W}}\Pi^{'}_{\gamma Z}(q^{2}), \tag{A6} $

      $ \begin{aligned}[b] Z_{Z^{*}}(q^{2})=&1+\frac{\mathrm{d}\Pi_{ZZ}(q^{2})}{\mathrm{d}q^{2}}|_{q^{2}=m^{2}_{Z}}-\Pi^{'}_{\gamma\gamma}(q^{2})\\&-\frac{c^{2}_{W}-s^{2}_{W}}{s_{W}c_{W}}\Pi^{'}_{\gamma Z}(q^{2}),\end{aligned} \tag{A7}$

      $ s^{2}_{W^{*}}(q^{2})=s^{2} _{WO}-s_{W}c_{W}\Pi^{'}_{\gamma Z}(q^{2}), \tag{A8} $

      $ e^{2}_{*}(q^{2})=e^{2} _{O}+e^{2}\Pi^{'}_{\gamma\gamma}(q^{2}), \tag{A9} $

      $ \Pi^{'}_{XY}(q^{2})=\frac{\big(\Pi_{XY}(q^{2})-\Pi_{XY}(0)\big)}{q^{2}}, \tag{A10} $

      $ \Pi_{XY}(q^{2})=\sum_{i} c_{i}\Pi^{(i)}_{XY}(q^{2}), \tag{A11} $

      $ Q_{\varphi q}^{33(3)} $

      $ \begin{aligned}[b] \Pi^{(1)}_{WW}=-N_c\frac{g^2}{4\pi^2}\frac{v^2}{\Lambda^2} \Bigg[ \left(\frac{1}{6}q^2-\frac{1}{4}(m_t^2+m_b^2)\right)E -q^2b_2(m_t^2,m_b^2,q^2)+\frac{1}{2}\Big(m_b^2b_1(m_t^2,m_b^2,q^2)+m_t^2b_1(m_b^2,m_t^2,q^2)\Big) \Bigg], \end{aligned}\tag{A12} $

      $ \begin{aligned}[b] \Pi^{(1)}_{ZZ}=&-N_c\frac{g^2}{\cos^2\theta_W}\frac{1}{4\pi^2}\frac{v^2}{\Lambda^2} \Bigg[ \Bigg(\frac{1}{6}(1-\sin^2\theta_W)q^2-\frac{1}{4}(m_t^2+m_b^2)\Bigg)E -q^2\Bigg(\left(\frac{1}{2}-\frac{2}{3}\sin^2\theta_W\right)b_2(m_t^2,m_t^2,q^2)\\&+\left(\frac{1}{2}-\frac{1}{3}\sin^2\theta_W\right)b_2(m_b^2,m_b^2,q^2)\Bigg) \left. +\frac{1}{4}\left(m_t^2b_0(m_t^2,m_t^2,q^2)+m_b^2b_0(m_b^2,m_b^2,q^2)\right) \right], \end{aligned}\tag{A13} $

      $ \Pi^{(1)}_{\gamma\gamma}=0, \tag{A14} $

      $ \begin{aligned}[b] \Pi^{(1)}_{\gamma Z}=-N_cg^2\frac{\sin\theta_W}{\cos\theta_W}\frac{1}{8\pi^2}\frac{v^2}{\Lambda^2} \Bigg[\frac{1}{6}E-\frac{2}{3}b_2(m_t^2,m_t^2,q^2)-\frac{1}{3}b_2(m_b^2,m_b^2,q^2)\Bigg]q^2, \end{aligned}\tag{A15} $

      $ Q_{\varphi q}^{33(1)} $

      $ \Pi^{(2)}_{WW}=0, \tag{A16} $

      $ \begin{aligned}[b] \Pi^{(2)}_{ZZ}=&N_c\frac{g^2}{\cos^2\theta_W}\frac{1}{4\pi^2}\frac{v^2}{\Lambda^2} \left[-\left(\frac{1}{4}m_t^2-\frac{1}{4}m_b^2+\frac{1}{18}q^2\sin^2\theta_W\right)E \right.\left. -q^2\left(\left(\frac{1}{2}-\frac{2}{3}\sin^2\theta_W\right)b_2(m_t^2,m_t^2,q^2) \right.\right.\\&-\left.\left(\frac{1}{2}-\frac{1}{3}\sin^2\theta_W\right)b_2(m_b^2,m_b^2,q^2)\right) \left. +\frac{1}{4}\left(m_t^2b_0(m_t^2,m_t^2,q^2)-m_b^2b_0(m_b^2,m_b^2,q^2)\right) \right], \end{aligned}\tag{A17} $

      $ \Pi^{(2)}_{\gamma\gamma}=0, \tag{A18} $

      $ \Pi^{(2)}_{\gamma Z}=N_cg^2\frac{\sin\theta_W}{\cos\theta_W}\frac{1}{8\pi^2}\frac{v^2}{\Lambda^2} \left[\frac{1}{18}E-\frac{2}{3}b_2(m_t^2,m_t^2,q^2)+\frac{1}{3}b_2(m_b^2,m_b^2,q^2)\right]q^2, \tag{A19} $

      $ Q^{33}_{\varphi u} $

      $ \Pi^{(3)}_{WW}=0, \tag{A20} $

      $ \Pi^{(3)}_{ZZ}=N_c\frac{g^2}{\cos^2\theta_W}\frac{1}{4\pi^2}\frac{v^2}{\Lambda^2} \left[\left(\frac{1}{4}m_t^2-\frac{1}{9}q^2\sin^2\theta_W\right)E \right.\left. -\left(\frac{1}{4}m_t^2b_0(m_t^2,m_t^2,q^2)-\frac{2}{3}q^2\sin^2\theta_Wb_2(m_t^2,m_t^2,q^2)\right) \right], \tag{A21}$

      $ \Pi^{(3)}_{\gamma\gamma}=0, \tag{A22} $

      $ \Pi^{(3)}_{\gamma Z}=N_cg^2\frac{\sin\theta_W}{\cos\theta_W}\frac{1}{12\pi^2}\frac{v^2}{\Lambda^2} \left(\frac{1}{6}E-b_2(m_t^2,m_t^2,q^2)\right)q^2, \tag{A23}$

      $ Q^{33}_{\varphi d} $

      $ \Pi^{(4)}_{WW}=0, \tag{A24}$

      $ \Pi^{(4)}_{ZZ}=N_c\frac{g^2}{\cos^2\theta_W}\frac{1}{4\pi^2}\frac{v^2}{\Lambda^2} \left[-\left(\frac{1}{4}m_b^2-\frac{1}{18}q^2\sin^2\theta_W\right)E \right.\left. +\left(\frac{1}{4}m_b^2b_0(m_b^2,m_b^2,q^2)-\frac{1}{3}q^2\sin^2\theta_Wb_2(m_b^2,m_b^2,q^2)\right) \right], \tag{A25} $

      $ \Pi^{(4)}_{\gamma\gamma}=0, \tag{A26} $

      $ \Pi^{(4)}_{\gamma Z}=-N_cg^2\frac{\sin\theta_W}{\cos\theta_W}\frac{1}{24\pi^2}\frac{v^2}{\Lambda^2} \left(\frac{1}{6}E-b_2(m_b^2,m_b^2,q^2)\right)q^2, \tag{A27} $

      $ Q^{33}_{\varphi ud} $

      $ \Pi^{(5)}_{WW}=-N_cg^2\frac{1}{16\pi^2}\frac{v^2}{\Lambda^2}m_tm_b\left(E-b_0(m_t^2,m_b^2,q^2)\right), \tag{A28} $

      $ \Pi^{(5)}_{ZZ}=0, \tag{A29} $

      $ \Pi^{(5)}_{\gamma\gamma}=0, \tag{A30} $

      $ \Pi^{(5)}_{\gamma Z}=0, \tag{A31} $

      $ Q^{33}_{uW} $

      $ \Pi^{(6)}_{WW}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_t}{\Lambda^2}\left(\frac{1}{2}E-b_1(m_b^2,m_t^2,q^2)\right)q^2, \tag{A32} $

      $ \Pi^{(6)}_{ZZ}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_t}{\Lambda^2}\left(\frac{1}{2}-\frac{4}{3}\sin^2\theta_W\right) \left(E-b_0(m_t^2,m_t^2,q^2)\right)q^2,\tag{A33} $

      $ \Pi^{(6)}_{\gamma\gamma}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_t}{\Lambda^2}\frac{4}{3}\sin^2\theta_W\left(E-b_0(m_t^2,m_t^2,q^2)\right)q^2, \tag{A34} $

      $ \Pi^{(6)}_{\gamma Z}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_t}{\Lambda^2}\frac{\sin\theta_W}{\cos\theta_W} \left(\frac{11}{12}-\frac{4}{3}\sin^2\theta_W\right)\left(E-b_0(m_t^2,m_t^2,q^2)\right)q^2, \tag{A35}$

      $ Q^{33}_{dW} $

      $ \Pi^{(7)}_{WW}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_b}{\Lambda^2}\left(\frac{1}{2}E-b_1(m_t^2,m_b^2,q^2)\right)q^2, \tag{A36} $

      $ \Pi^{(7)}_{ZZ}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_b}{\Lambda^2}\left(\frac{1}{2}-\frac{2}{3}\sin^2\theta_W\right) \left(E-b_0(m_b^2,m_b^2,q^2)\right)q^2, \tag{A37} $

      $ \Pi^{(7)}_{\gamma\gamma}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_b}{\Lambda^2}\frac{2}{3}\sin^2\theta_W\left(E-b_0(m_b^2,m_b^2,q^2)\right)q^2, \tag{A38} $

      $ \Pi^{(7)}_{\gamma Z}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_b}{\Lambda^2}\frac{\sin\theta_W}{\cos\theta_W} \left(\frac{7}{12}-\frac{2}{3}\sin^2\theta_W\right)\left(E-b_0(m_b^2,m_b^2,q^2)\right)q^2, \tag{A39}$

      $ Q^{33}_{uB} $

      $ \Pi^{(8)}_{WW}=0, \tag{A40} $

      $ \Pi^{(8)}_{ZZ}=N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_t}{\Lambda^2}\frac{\sin\theta_W}{\cos\theta_W}\left(\frac{1}{2}-\frac{4}{3}\sin^2\theta_W\right) \left(E-b_0(m_t^2,m_t^2,q^2)\right)q^2, \tag{A41} $

      $ \Pi^{(8)}_{\gamma\gamma}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_t}{\Lambda^2}\frac{4}{3}\sin\theta_W \cos\theta_W\left(E-b_0(m_t^2,m_t^2,q^2)\right)q^2, \tag{A42} $

      $ \Pi^{(8)}_{\gamma Z}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_t}{\Lambda^2} \left(\frac{1}{4}-\frac{4}{3}\sin^2\theta_W\right)\left(E-b_0(m_t^2,m_t^2,q^2)\right)q^2, \tag{A43}$

      $ Q^{33}_{dB} $

      $ \Pi^{(9)}_{WW}=0, \tag{A44} $

      $ \Pi^{(9)}_{ZZ}=-N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_b}{\Lambda^2}\frac{\sin\theta_W}{\cos\theta_W}\left(\frac{1}{2}-\frac{2}{3}\sin^2\theta_W\right) \left(E-b_0(m_b^2,m_b^2,q^2)\right)q^2,\tag{A45}$

      $ \Pi^{(9)}_{\gamma\gamma}=N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_b}{\Lambda^2}\frac{2}{3}\sin\theta_W \cos\theta_W\left(E-b_0(m_b^2,m_b^2,q^2)\right)q^2, \tag{A46}$

      $ \Pi^{(9)}_{\gamma Z}=N_cg\frac{\sqrt{2}}{4\pi^2}\frac{vm_b}{\Lambda^2} \left(\frac{1}{4}-\frac{2}{3}\sin^2\theta_W\right)\left(E-b_0(m_b^2,m_b^2,q^2)\right)q^2, \tag{A47} $

      where $ \theta_W $ is the weak angle, $ N_c=3 $ is the number of colors, $ E=\dfrac{2}{4-d}-\gamma+\ln 4\pi $, and the functions $ b_i $ are given by

      $ b_0(m_1^2,m_2^2,q^2)=\int_0^1\ln\frac{(1-x)m_1^2+xm_2^2-x(1-x)q^2}{\mu^2}{\rm d}x, \tag{A48} $

      $ b_1(m_1^2,m_2^2,q^2)=\int_0^1x\ln\frac{(1-x)m_1^2+xm_2^2-x(1-x)q^2}{\mu^2}{\rm d}x, \tag{A49} $

      $ b_2(m_1^2,m_2^2,q^2)=\int_0^1x(1-x)\ln\frac{(1-x)m_1^2+xm_2^2-x(1-x)q^2}{\mu^2}{\rm d}x, \tag{A50} $

      where μ is the 't Hooft mass. They have the following analytical expressions:

      $ b_0(m_1^2,m_2^2,q^2)=-2+\log\frac{m_1m_2}{\mu^2}+\frac{m_1^2-m_2^2}{q^2}\log\left(\frac{m_1}{m_2}\right) +\frac{1}{q^2}\sqrt{|(m_1+m_2)^2-q^2||(m_1-m_2)^2-q^2|}f(m_1^2,m_2^2,q^2), \tag{A51} $

      where

      $ \begin{aligned} f(m_1^2,m_2^2,q^2)=\left\{ \begin{array}{ll} \log\dfrac{\sqrt{(m_1+m_2)^2-q^2}-\sqrt{(m_1-m_2)^2-q^2}}{\sqrt{(m_1+m_2)^2-q^2}+\sqrt{(m_1-m_2)^2-q^2}} &\quad q^2\leq(m_1-m_2)^2\\ 2\arctan\sqrt{\dfrac{q^2-(m_1-m_2)^2}{(m_1+m_2)^2-q^2}} &\quad (m_1-m_2)^2<q^2<(m_1+m_2)^2\\ \log\dfrac{\sqrt{q^2-(m_1-m_2)^2}+\sqrt{q^2-(m_1+m_2)^2}}{\sqrt{q^2-(m_1-m_2)^2}-\sqrt{q^2-(m_1+m_2)^2}} &\quad q^2\geq(m_1+m_2)^2\\ \end{array}\right., \end{aligned} \tag{A52}$

      and

      $ \begin{aligned}[b] b_1(m_1^2,m_2^2,q^2)=&-\frac{1}{2}\left[\frac{m_1^2}{q^2}\left(\log\frac{m_1^2}{\mu^2}-1\right)-\frac{m_2^2}{q^2}\left(\log\frac{m_2^2}{\mu^2}-1\right)\right]+\frac{1}{2}\frac{m_1^2-m_2^2+q^2}{q^2}b_0(m_1,m_2,q),\\ b_2(m_1^2,m_2^2,q^2)=&\frac{1}{18}+\frac{1}{6}\left[\frac{m_1^2(2m_1^2-2m_2^2-q^2)}{(q^2)^2}\log\frac{m_1^2}{\mu^2} +\frac{m_2^2(2m_2^2-2m_1^2-q^2)}{(q^2)^2}\log\frac{m_2^2}{\mu^2}\right]\\ &-\frac{1}{3}\left(\frac{m_1^2-m_2^2}{q^2}\right)^2-\frac{1}{6}\left[2\left(\frac{m_1^2-m_2^2}{q^2}\right)^2-\left(\frac{m_1^2+m_2^2+q^2}{q^2}\right)\right]b_0(m_1,m_2,q). \end{aligned}\tag{A53} $

    • A.4.   Examples of tree and loop level contributions

    • According to [48], we can calculate the loop contribution of the following observables:

      1. The decay width of $ Z\to e^{+}e^{-} $, $ \Gamma_{Z\to ee} $. The numerical expressions are

      $ \delta \Gamma^\mathrm{tree}_{Z\to ee}= (10.7 c^{(1)}_{\varphi l}-1.15 c^{(3)}_{\varphi l}-2.96 c_{\varphi D}-9.40 c_{\varphi e}-2.08 c_{\varphi WB}+5.93 c^{'}_{ll})\times 10^{-3}. \tag{A54} $

      $ \delta\Gamma^\mathrm{loop}_{Z\to ee}= (-124c^{(-)}_{\varphi Q}+3.24 c^{(+)}_{\varphi Q}+140c_{\varphi t} -4.56c_{\varphi b}-2.11c_{\varphi tb} +12.3 c_{tW}+1.33c_{bW} +17.6 c_{tB}+4.45c_{bB})\times10^{-6}. \tag{A55} $

      2. The decay width of $ W\to l\nu_{l} $, $ \Gamma_{W\to l\nu_{l}} $.

      $ \delta \Gamma^\mathrm{tree}_{W\to l\nu_{l}}=(-1.77c^{(3)}_{\varphi l}-1.45c_{\varphi D}-3.21c_{\varphi WB}+2.23c^{'}_{ll})\times 10^{-2}, \tag{A56} $

      $ \delta\Gamma^\mathrm{loop}_{W\to l\nu_{l}}= (-537c^{(-)}_{\varphi Q}+126c^{(+)}_{\varphi Q}+629c_{\varphi t}-34.3c_{\varphi b} -12.8c_{\varphi tb} -64.7c_{tW}+8.53c_{bW}+245c_{tB} +41.5c_{bB})\times10^{-6}. \tag{A57} $

      3. The coupling between the axial-vector current and Z boson and the coupling between the vector current and Z boson in the $ \nu_{\mu}-e $ scattering process at low energy, $ g^{\nu_{\mu}e}_{LV} $, $ g^{\nu_{\mu}e}_{LA} $ [95].

      $ \delta g^{\nu_{\mu}e}_\mathrm{LVtree}= (-2.83c^{(1)}_{\varphi l}+5.32c^{(3)}_{\varphi l}+2.13c_{\varphi D}-3.03c_{\varphi e} +9.63 c_{\varphi WB}-3.03c_{le}-6.06c_{ll}-4.27c^{'}_{ll})\times10^{-2},\tag{A58} $

      $ \delta g^{\nu_{\mu}e}_\mathrm{LAtree}=(1.51c_{\varphi D} + 3.03 c_{\varphi e} + 3.03c_{le} - 6.06c_{ll} - 3.03c^{'}_{ll})\times10^{-2}. \tag{A59} $

      $ \delta g^{\nu_{\mu} e}_\mathrm{LVloop}= (814c^{(-)}_{\varphi Q}+6.73 c^{(+)}_{\varphi Q}-896c_{\varphi t}+213c_{\varphi b} +15.2c_{\varphi tb}-715c_{t W}+111 c_{b W}-747c_{t B}-170c_{bB})\times10^{-6}, \tag{A60}$

      $ \delta g^{\nu_{\mu} e}_\mathrm{LAloop}=(647c^{(-)}_{\varphi Q}+87.1c^{(+)}_{\varphi Q}-729c_{\varphi t}-5.5072c_{\varphi b}+10.8c_{\varphi tb})\times10^{-6}. \tag{A61} $

    • A.5.   $\boldsymbol {e^{+}e^{-}\to W^{+}W^{-}}$ and $\boldsymbol {e^{+}e^{-}\to W^{+}W^{-}\to l\nu_{l}ud}$

    • The W pair production cross-section and $ W^{-} $ angular distribution are measured at LEP-II. However, MadGraph cannot calculate the change in space phase induced by $ \delta m^{2}_{W}(m^{2}_{W}) $, $ \delta m_{W} $, and $ \delta\Gamma_{W} $. $ \delta\Gamma_{W} $ is the third-generation quark loop and tree level contribution to the decay width of W and $ \Gamma_{W} $, respectively. Therefore, the calculation of the effect induced by $ \delta m^{2}_{W}(m^{2}_{W}) $, $ \delta m_{W} $, and $ \delta\Gamma_{W} $ is performed with Feyncalc [96, 97] and Feynarts [98]. The remaining contributions from dimension-6 operators to the W pair production cross-section and $ W^{-} $ angular distribution can be calculated with MadGraph. Next, we demonstrate how to calculate these contributions with Feynarts and Feyncalc in detail.

      $ \delta\sigma_{\rm{on-shell}} $: $ \delta\sigma_{\rm{on-shell}} $ is the contribution of $ \delta m^{2}_{W}(m^{2}_{W}) $ and $ \delta m_{W} $ to $ \sigma_{\rm{on-shell}} $. $ \sigma_{\rm{on-shell}} $ is the W pair production cross-section in the process $ e^{+}e^{-}\to W^{+}W^{-} $ at the SM tree level.

      $ \delta\sigma_{\rm{on-shell}}=\partial_{m_{W}}\sigma_{\rm{on-shell}}\Bigg(\frac{\delta m^{2}_{W}(m^{2}_{W})}{2m_{W}}+\delta m_{W}\Bigg).\tag{A62} $

      $\dfrac{{\rm d}^{i} \delta\sigma_{\rm{off-shell}}}{{\rm d} \cos\theta}$: $ \dfrac{{\rm d}^{i} \delta\sigma_{\rm{off-shell}}}{{\rm d} \cos\theta} $ is the contribution of $ \delta m^{2}_{W}(m^{2}_{W}) $, $ \delta m_{W} $, and $ \delta\Gamma_{W} $ to $ \dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta} $. $ \dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta} $ is the $ W^{-} $ angular distribution (θ is the polar angle between the $ W^{-} $ and $ e^{-} $ beams) in the ith bin at the SM tree level in the process $ e^{+}e^{-}\to W^{+}W^{-}\to l\nu_{l}ud $.

      We can use the phase space recursion relation and narrow width approximation [99] to calculate $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d} \cos\theta}$(this approximation will have only a 1% deviation from the total cross-section).

      $ \frac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta}\approx \frac{{\rm d}^{i}\sigma_{\rm{on-shell}}}{{\rm d}\cos\theta} \frac{\Gamma_{W\to l\nu_{l}}}{\Gamma_{W}}\frac{\Gamma_{W\to ud}}{\Gamma_{W}}, \tag{A63}$

      where $ \Gamma_{W\to ud} $ is the decay width of $ W\to ud $ at the SM tree level, and $\dfrac{{\rm d}^{i}\sigma_{\rm{on-shell}}}{{\rm d}\cos\theta}$ is the $ W^{-} $ angular distribution in the ith bin at the SM tree level in the process $ e^{+}e^{-}\to W^{+}W^{-} $. Then, $\dfrac{{\rm d}^{i} \delta\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta}$ is

      $ \begin{aligned}[b] \frac{{\rm d}^{i} \delta\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta}=&\frac{{\rm d}^{i} \delta\sigma_{\rm{on-shell}}}{{\rm d} \cos\theta}\frac{\Gamma_{W\to l\nu_{l}}}{\Gamma_{W}}\frac{\Gamma_{W\to ud}}{\Gamma_{W}}\\&+\frac{{\rm d}\sigma_{\rm{on-shell}}}{{\rm d}\cos\theta} \frac{\delta\Gamma_{W\to l\nu_{l}}}{\Gamma_{W}}\frac{\Gamma_{W\to ud}}{\Gamma_{W}}\\&+\frac{{\rm d}\sigma_{\rm{on-shell}}}{{\rm d}\cos\theta} \frac{\Gamma_{W\to l\nu_{l}}}{\Gamma_{W}}\frac{\delta\Gamma_{W\to ud}}{\Gamma_{W}}\\&-2\frac{{\rm d}\sigma_{\rm{on-shell}}}{{\rm d}\cos\theta} \frac{\Gamma_{W\to l\nu_{l}}}{\Gamma_{W}}\frac{\Gamma_{W\to ud}}{\Gamma_{W}}\frac{\delta\Gamma_{W}}{\Gamma_{W}}, \end{aligned}\tag{A64} $

      where $ \delta\Gamma_{W\to l\nu_{l}} $ and $ \delta\Gamma_{W\to ud} $ are the $ \delta m^{2}_{W}(m^{2}_{W}) $ and $ \delta m_{W} $ contributions to $ \Gamma_{W\to l\nu_{l}} $and $ \Gamma_{W\to ud} $, respectively. $ \dfrac{{\rm d}^{i} \delta\sigma_{\rm{on-shell}}}{{\rm d} \cos\theta} $ is the $ \delta m^{2}_{W}(m^{2}_{W}) $ and $ \delta m_{W} $ contributions to $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d} \cos\theta}$. $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d} \cos\theta}$ is the $ W^{-} $ angular distribution in the ith bin in the process $ e^{+}e^{-}\to W^{+}W^{-} $ at the SM tree level.

      However, the complete $ \mathcal{O}(\alpha) $ electroweak corrections cannot be ignored. We must calculate the contribution of dimension-6 operators including complete $ \mathcal{O}(\alpha) $ electroweak corrections.

    • A.5.1.   Contribution of $ \delta m^{2}_{W}(m^{2}_{W}), \delta m_{W} $ to $ \sigma_{\rm{on-shell}}^\mathrm{NLO} $
    • In the process $ e^{+}e^{-}\to W^{+}W^{-} $, our observables are $ \sigma_{\rm{on-shell}}^\mathrm{NLO} $ at different $ \sqrt{s} $. $ \sigma_{\rm{on-shell}}^\mathrm{NLO} $ is the total cross-section of W pair production including complete $ \mathcal{O}(\alpha) $ electroweak correction in the process $ e^{+}e^{-}\to W^{+}W^{-} $. The assumption is made that the operators' contribution to $ \sigma_{\rm{on-shell}}^\mathrm{NLO} $ can be calculated in the following way:

      $ \delta\sigma_{\rm{on-shell}}^\mathrm{NLO}\approx\frac{\sigma_{\rm{on-shell}}^\mathrm{NLO}}{\sigma_{\rm{on-shell}}}\delta\sigma_{\rm{on-shell}} \tag{A65} $

      where $ \delta\sigma_{\rm{on-shell}}^\mathrm{NLO} $ is the contribution of $ \delta m^{2}_{W}(m^{2}_{W}) $ and $ \delta m_{W} $ to $ \sigma_{\rm{on-shell}}^\mathrm{NLO} $.

    • A.5.2.   Contribution of $ \delta m^{2}_{W}(m^{2}_{W}) $, $ \delta m_{W} $, and $ \delta\Gamma_{W} $ to $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d} \cos\theta}$
    • In the process $ e^{+}e^{-}\to W^{+}W^{-}\to l\nu_{l}ud $, the observables are $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d}\cos\theta}$. $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d}\cos\theta}$ is the $ W^{-} $ angular distribution in the ith bin including complete $ \mathcal{O}(\alpha) $ electroweak correction. Because $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d}\cos\theta}$ changes very slowly in every certain bin ($ \cos\theta_{i}=(i-1)\times0.2 $), it can be regarded as $\dfrac{{\rm d}\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d}\cos\theta}\Big|_{\cos\theta=\textstyle\frac{\cos\theta_{i}+\cos\theta_{i+1}}{2}}$. The reweighting method is applied to calculate $\dfrac{{\rm d}^{i}\delta\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d} \cos\theta}$ from $ \dfrac{{\rm d}^{i}\delta\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta} $. $ \dfrac{{\rm d}^{i}\delta\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta} $ and $\dfrac{{\rm d}^{i}\delta\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d} \cos\theta}$ are the contributions of $ \delta m^{2}_{W}(m^{2}_{W}) $, $ \delta m_{W} $, and $ \delta\Gamma_{W} $ to $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d} \cos\theta}$ and $\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}^{\mathrm{NLO}}}{{\rm d} \cos\theta}$, respectively.

      $ \frac{{\rm d}^{i}\delta\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d} \cos\theta}\approx\frac{\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}^\mathrm{NLO}}{{\rm d}\cos\theta}}{\dfrac{{\rm d}^{i}\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta}}\frac{{\rm d}^{i}\delta\sigma_{\rm{off-shell}}}{{\rm d}\cos\theta}. \tag{A66} $

      The numerical result of ${\rm d}^{1}\delta\sigma_{\rm{off-shell}}^\mathrm{NLO}$ when $ \sqrt{s} $ is 182.66 GeV.

      $ \begin{aligned}[b] {\rm d}^{1}\delta\sigma_{\rm{off-shell}}^\mathrm{NLO}=&(-114.082c^{(-)}_{\varphi Q}+27.4928 c^{(+)}_{\varphi Q} + 132.446c_{\varphi t} \\&- 3615.79c^{(3)}_{\varphi l} -361.279c^{(3)}_{\varphi q}-7.22494c_{\varphi b} \\&- 2.66921c_{\varphi tb} - 5.95838c_{t W}+2.23408c_{b W}\\&+ 51.7414c_{tB} + 8.73857c_{bB} - 3260.26c_{\varphi D} \\ &-7201.74c_{\varphi WB}+ 1446.62c^{'}_{ll})\times 10^{-5}. \end{aligned}\tag{A67} $

Reference (99)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return