Theory requirements for SM Higgs and EW precision physics at the FCC-ee

High precision experimental measurements of the properties of the Higgs boson at $\sim$ 125 GeV as well as electroweak precision observables such as the W -boson mass or the effective weak leptonic mixing angle are expected at future $e^+e^-$ colliders such as the FCC-ee. This high anticipated precision has to be matched with theory predictions for the measured quantities at the same level of accuracy. We briefly summarize the status of these predictions within the Standard Model (SM) and of the tools that are used for their determination. We outline how the theory predictions will have to be improved in order to reach the required accuracy, and also comment on the simulation frameworks for the Higgs and EW precision program.


Introduction
With the discovery of the Higgs boson, all possible elements of the Standard Model (SM) have been experimentally confirmed and tested in great depth. On the other hand, observational evidence for dark matter and the matterantimatter asymmetry require physics beyond the SM. One promising way to probe such new physics is through precision measurements of the properties of the Higgs boson. A complementary way is to measure electroweak precision (pseudo-)observables (EWPO) with highest precision. These are the avenues pursued by several proposals for a future e + e − collider. In particular, the FCC-ee concept is designed to run at √ s = 250 GeV as a Higgs factory, and at √ s ∼ M Z , 2M W for high-precision EWPO measurements. In this way it can improve (indirect) probes for new physics by several orders of magnitude compared to existing bounds [1][2][3].
The anticipated experimental accuracy of an observable has to be matched with a theory prediction of at least the same level of accuracy to make maximum use of the experimental data. For the accurate study of the properties of the Higgs boson, precise predictions for the various partial decay widths, the branching ratios (BRs) and the Higgs-boson production cross sections along with their theoretical uncertainties are indispensable. For the EWPO on the one hand the SM prediction has to reach the level of the experimental accuracy. Similarly, the extraction of these quantities from experimental data must be equally well under control. All theses types of uncertainties (various theoretical ones and experimental) must be taken into account when deriving constraints on new physics from the data.
Several sources of theory uncertainties have to be distinguished. The intrinsic uncertainties are due to missing higher orders in the perturbative expansion of the SM (or BSM) prediction for an observable. The parametric uncertainties are due to the imperfect experimental knowledge of the SM input parameters (as well as theory uncertainties induced in their extraction from data). The extraction of a quantity from a cross section or an asymmetry requires the theory prediction of this cross section or asymmetry to at least the same order of precision 1 .
In this paper the current status and future implications of theory uncertainties on (SM) Higgs-boson physics and EWPO's will be summarized. While we focus on the FCC-ee, they are valid for all future high precision e + e − colliders running above the HZ threshold (such as ILC, CLIC, or CEPC). We will use anticipated FCC-ee precisions to illustrate the impact of theory uncertainties. We build on [4], where also a discussion of the determination of the SM input parameters can be found. In addition, we review Monte Carlo (MC) event generator frameworks which are indispensable for simulations in the perturbative and non-perturbative regime, building upon the experience gained from the LEP era [5][6][7] and the extensive simulations for CEPC, ILC and CLIC. Work on dedicated FCC-ee MC samples has started.

Pseudo-observables versus realistic observables
The quantities that can be directly measured in experiments are cross sections, line shape observables, forwardbackward asymmetries etc., called "realistic observables" [8,9]. The obtained results depend on the specific set of experimental cuts, however detector inefficiencies are removed and kinematic cuts are brought to simplified form using MC event generators. In order to determine idealized quantities like masses, partial widths or couplings, or the effective electroweak mixing angle, from the primarily measured realistic observables, a "QED unfolding" procedure is applied. This procedure involves manipulations like subtracting photon-exchange or interference terms, subtracting box-diagram contributions, unfolding higher-order QED or QCD corrections, etc. These secondary quantities are therefore called "pseudo-observables". The Higgs-boson observables and EWPOs on which we focus here are pseudoobservables. It should be stressed that the role of EWPOs is to represent and encapsulate experimental data in the model-independent way 2 . Fitting EWPOs to the SM and its extensions comes as an independent next step. In order to reach the required theoretical accuracy special care has to be taken in the application of these unfolding proceduresin particular the procedure of the extraction of EWPOs from data causes certain bias, partial loss of information on physics in the realistic data. The unfolding procedure for the LEP experiments is described in Refs. [8,9]. A discussion of the additional problems in this respect arising in beyond the SM (BSM) models can be found in [10]. In Refs. [8,9] it was proven that the bias induced in the extraction of EWPOs (QED deconvolution) was smaller than the LEP experimental errors. This will not be true in the FCC-ee era where experimental errors will be smaller by up to two orders. The very survival of the technique of the model independent EWPOs as the bridge between the experiment and the theory at FCC-ee will require significant improvement of the precision of the trivial but large QED effects. Ref. [11] covers this issue in a great detail, listing challenges and outlining solutions. In addition, the definition of EWPOs and the algorithm of their extraction will have to be modified. This is elaborated in Sect. 3 in Ref. [12]. Two major changes are anticipated. One is that we may be forced to remove from EWPOs not only QED but also first order pure electroweak corrections, in order to achieve sharper resolution in the search of the BSM effects in the new EWPOs. Another important change will be the increased role of Monte Carlo generators in extracting EWPOs from data -this is because semianalytic programs like ZFITTER and TOPAZ0 of LEP era will not be able to handle precisely enough detector geometry (cut-offs).

Electroweak precision observables
The most important electroweak precision observables (EWPOs) are related to properties of the Z and W bosons. The Z-boson properties are determined from measurements of e + e − → ff on the Z-pole. To isolate the physics of the Z-boson, the typical set of pseudo-observables is defined in terms of the de-convoluted cross-section σ f (s), where the effect of initial-and final-state photon radiation and from s-channel photon and double-boson (box) exchange has been removed. The customary set of pseudo-observables are total and partial Z-boson cross sections with their subsequent decay to quarks or leptons, forward-backward and left-right asymmetries, from which the effective weak leptonic mixing angle, sin 2 θ f eff , is extracted, as well as the W -boson mass, M W , see Ref. [13]. The expected (pure) experimental precision at FCC-ee of M W and sin 2 θ ℓ eff is summarized in Tab. 1 based on Refs. [2,3]. The quantities listed in Tab. 1 can be predicted within the SM in terms of the input parameters G F , α(M Z ), α s (M Z ), M Z , M H and m t . The radiative corrections in these predictions are currently known including complete two-loop corrections In addition, approximate three-and four-loop s ) are available (with α t = y 2 t /(4π) and y t being the top Yukawa coupling), see [14] for a review. The theory uncertainties from missing higher-order corrections are given in the third column of Tab. 1 [15][16][17][18]. Also indicated are the main sources for the respective uncertainties. In order to match the FCC-ee precision the theory predictions will have to be improved substantially. It was estimated in [4,12,19] that the intrinsic uncertainties will be reduced by corrections, and leading fourloop corrections enhanced by the top Yukawa coupling, resulting in the projected intrinsic uncertainties shown in column four of Tab. 1.
A crucial point here is that the determination of the pseudo-observables in Tab. 1 from experimental data also requires theory input for the removal of initial-state and final-state photon radiation and s-channel photon exchange and box contributions. The theory uncertainty from missing higher QED orders is estimated to amount to a few times 0.01% [5,8,9,20] for the Z-peak cross-section and total width measurements. In order not to be dominated by this uncertainty, it will need to be reduced by about a factor 10 for the FCC-ee. This will require the calculation of non-leading log contributions to O(α 3 ) corrections, O(α 3 L 2 ) and O(α 4 L 4 ) contributions (with L ≡ log s/m 2 e ) as well as an improved treatment of fermion pair production from off-shell photons.
The W -boson mass will be determined from a threshold scan near the W -pair threshold, √ s ≈ 161 GeV. It is foreseen that the experimental uncertainty at FCC-ee for this measurement is about 0.5 MeV [1][2][3]. At the point of highest sensitivity, an uncertainty of the cross-section measurement of 0.1% translates to an uncertainty of ∼ 1.5 MeV on M W [21]. Therefore a theoretical prediction for the process e + e − → 4f with an accuracy of ∆σ ∼ 0.01% is desirable, including effects of off-shell W bosons, which become important near threshold.
The currently best calculations are based on complete one-loop results for e + e − → 4f [22] and partial higher-order effects for the total cross section from an effective field theory framework [23,24]. The resulting theory uncertainty on M W is estimated to be about 3 MeV [24]. This result must be improved by complete 2-loop corrections to e + e − → W + W − and to W → ff ′ (based on the effective field theory framework). In addition, a more accurate description of initial-state radiation will be important, which includes universal contributions from soft and collinear photon radiation (see [5,8,20] for a review), as well as hard photon radiation. For the latter, a proper matching and merging procedure needs to be employed to avoid double counting [23,25], see [4] for more details. There it was estimated that a theory induced systematic uncertainty of ∆M W < ∼ 0.60 MeV can be feasible. At the FCC-ee with √ s = 240 GeV (or other e + e − machines near this center-of-mass energy), the Higgs-boson production cross-section is strongly dominated by e + e − → ZH, and e + e − → ννH contributes less than 20% [1,26]. For these two processes full one-loop corrections in the SM are available [27,28]. For the dominating ZH production mode they are found at the level of ∼ 5−10%. Leading two-loop corrections to this cross section were evaluated at O(αα s ) [29,30], which were found to amount to ∼ 1 − 2%. This number has to be compared to the anticipated experimental accuracy of 0.4% [1][2][3]. It becomes clear that with a full two-loop calculation of e + e − → ZH the intrinsic uncertainty will be sufficiently small. Calculational techniques for 2 → 2 processes at the two-loop level exist, and it is reasonable to assume that, if required, this calculation within the SM can be incorporated for the FCC-ee Higgs precision studies (see also [31] for recent progress in this direction.) For WBF production, the calculation of the full two-loop corrections will be significantly more difficult, since this is a 2 → 3 process. However, one may assume that a partial result based on diagrams with closed light-fermion loops and top-quark loops (in a large-m t approximation) can be achieved, which should reduce the intrinsic theory uncertainty to below the 1% level. Given the fact that the WBF process is less crucial than the HZ channel for the Higgs physics program FCC-ee with √ s = 240 GeV, this will probably be adequate for most practical purposes. It would be desirable to have a complete 2 → 4 calculation at two loops, including both processes with their interference, but very likely factorized processes with on-shell projections will be sufficient.
Also the parametric uncertainties can play a non-negligible role for the evaluation of the partial widths. The most important parameters are the bottom quark mass and the strong coupling constant. In Ref. [33] the current uncertainties of α s and m b have been assumed to be δα s = 0.0015 and δm b = 0.03 GeV. Additionally, δm c = Partial width intr. QCD intr. electroweak total para. mq para. αs para.  Table 3. Projected intrinsic and parametric uncertainties for the partial and total Higgs-boson decay width predictions (see text). The last column shows the target of FCC-ee precisions on the respective coupling squared. 0.025 GeV, δm t = 0.85 GeV and δM H = 0.24 GeV [36] have been taken into account. The effect on the various partial widths has been evaluated as in Ref. [32] and is shown in the three right columns of Tab. 2.
When comparing the combined intrinsic and parametric uncertainties with the target precision of FCC-ee [1][2][3], see the last column in Tab. 3, it is clear that improvements are necessary. Concerning the intrinsic theory uncertainty, the available predictions for the ff and γγ channels are already sufficiently precise to match the expected FCC-ee experimental uncertainty. With available calculational techniques, the evaluation of complete two-loop corrections to H → ff can be achieved. This would reduce the uncertainty of the electroweak contributions to less than 0.1%. Similarly, the complete NLO corrections to H → Zγ can be carried out with existing methods, resulting in an estimated precision of about 1%. More theoretical work is needed for H → W W, ZZ, gg, which are currently limited by QCD uncertainties. For H → gg the calculation of massless four-loop QCD diagrams, will be required, which may be within reach [37], reducing the intrinsic uncertainties to the level of about 1%. For H → W W, ZZ, the required QCD corrections are essentially identical to those for e + e − → W W , and it is straightforward to improve them to a practically negligible level. Further significant progress would require the calculation of two-loop electroweak corrections, which for a 1 → 4 process is beyond reach for the forseeable future. Here it should be noted, however, that the HZZ coupling will be mostly constrained by the measurement of the e + e − → HZ production process at FCC-ee with √ s = 240 GeV, rather than the decay H → ZZ * , see the discussion in section 4.1, leading to a remaining intrinsic uncertainty of less than 0.3%.
Also shown in Tab. 3 are the projected parametric uncertainties, assuming FCC-ee precisions. For inputs, we use δα s = 0.0002 and δm t = 50 MeV 3 , δM H ∼ 10 MeV, and δm b ∼ 13 MeV and δm c ∼ 7 MeV, see [4] for details. 4 It becomes clear that the intrinsic and parametric uncertainties will be able to match the high anticipated experimental accuracy at FCC-ee.

Monte-Carlo and exclusive predictions
To determine systematic uncertainties, to properly extract parameters from the Higgs measurements and to find deviations from SM predictions, it is indispensable to have Monte Carlo (MC) simulations. The challenge compared to LEP/SLC times is the strongly increased experimental precision and much higher luminosities which demands much better simulations. The precise description of QED effects, both regarding the correct normalization of cross sections as well as exclusive multi-photon radiation will be one of the highest priorities. In the recent years, there has been steady progress in dealing with the calculation of higher logarithmic and finite orders in the QED radiation, enhanced by the very small electron mass, see Refs. [38][39][40][41][42]. These very advanced calculations are only partially applicable for the exclusive description of photon radiation which is necessary if cuts are to be applied on those photons. At the level of LO matrix elements, there is a matching procedure between explicit matrix element photons and the electron PDFs in collinear factorization in the context of WHIZARD (cf. below) which has been recently shown to work unexpectedly well [43,44].
From the LEP era there are many dedicated MC programs which feature soft photon resummation in the exclusive form. Most of them YFS2/KORALZ [45], BHLUMI [46], BHWIDE [47], YFSWW [48] are based directly on the Yennie-Frautschi-Suura (YSF) resummation scheme [49]. A more powerful coherent exclusive exponentiation (CEEX) at the amplitude level was introduced in [20] and is so far implemented only in KKMC [50] program for two-fermion processes e + e − → ff , f = µ, τ, ν, q. A generic semi-automatic implementation of CEEX for a wider class of processes would be a very desirable development, and might pave the path towards a consistent inclusion of the non-soft higher order contributions over the complete photon phase space, merging together genuine electroweak corrections with QED correction calculated to much higher orders, inclusion of spin correlations, initial-final state interferences and more. The ultimate goal is to match non-soft QED corrections up to forth order and genuine SM electroweak corrections up to second order.
Multi-purpose event generators allow to simulate the full spectrum of SM (and BSM) processes at FCC-ee. They provide (for the moment) a lower level of precision for the two-fermion processes, but allow to simulate them in the same framework as four-, six-and event eight-fermion processes. Due to the much higher precision and better detectors than at SLC/LEP these processes have to be included already at the 250 GeV stage, and not only if triple electroweak resonances or top pairs are kinematically accessible. In this context a lot of experience has been gained from the large full SM event samples generated for ILC, CLIC and CEPC using WHIZARD [51,52]. This supports leading-logarithmic QED collinear factorization for QED photons, automated NLO QCD corrections for arbitrary processes, while NLO electroweak corrections are under way. For experimental studies, crucial observables are the correct numbers of neutral and charged hadrons. These can only be correctly simulated if the color-flow assignments are correctly made, which includes that the invariant masses of the shower systems from the most prominent processes, e + e − → W + W − , ZZ are preserved; in a full electroweak matrix element they are, however, quantum mechanically entangled. The MC hence needs to determine the correct probability of the resonant subprocesses and hand them over to parton shower and hadronization [53,54]. It is important to note that the solution of the long-standing problem of the soft photon resummation to narrow width semi-stable charged resonances like W ± was recently outlined in Ref. [55]. It is also very important to take photon-induced backgrounds from Weizsäcker-Williams QED splittings into account, particularly γγ → hadrons. WHIZARD has a special treatment for the top threshold run [56] that is used by the FCC-ee collaboration. It also supports all different kinds of beam spectra, from simple parameterized files, over Gaussian beam spreads to a sophisticated beamstrahlung setup for linear and circular lepton machines [57].
To search for deviations in the Higgs sector from the SM, WHIZARD supports a large setup of BSM models like supersymmetry, extended Higgs models like 2HDM, Higgs singlet extensions, SMEFT, composite and Little Higgs models, which can be generically extended to its interface to automated Lagrangian model tools, cf. e.g. [58,59] 6 Outlook and deliverables The most important steps to be taken in the future to match the high anticipated accuracy of FCC-ee in the context of Higgs physics and EWPO with SM theory predictions are the following.
-Improved unfolding techniques to go from observables to pseudo-observables. -Improvements by roughly one order of magnitude in the determination of the most important SM parameters (m t , m b , α s , ∆α had , M H , . . . ) -Complete simulations of relevant SM processes at different energy stages in full detector simulation Finally, the following should be kept in mind. The SM constitutes the model in which highest theoretical precision for the predictions of EWPO and Higgs-boson observables can be obtained. This concerns their prediction as well as their extraction. However, as soon as BSM physics will be discovered, an evaluation of the predictions of the EWPO and the Higgs-boson observables in any preferred BSM model will be necessary. The corresponding theory uncertainties, both intrinsic and parametric, can then be larger (see, e.g., [10,60] for the Minimal Supersymmetric SM). A dedicated theory effort (beyond the SM) would be needed in this case.