NNLO EW Fermionic Corrections

• NNLO electroweak fermionic corrections are order α² contributions with closed fermion loops that refine theoretical predictions for processes like ZH production and Drell–Yan.
• Semi-numerical two-loop techniques combine Feynman parameters, dispersion relations, and Passarino–Veltman reduction to achieve stable numerical integration and precision.
• Robust UV and IR subtraction along with strategic renormalization schemes significantly reduce theoretical uncertainties, aligning predictions with experimental targets.

Next-to-next-to-leading order (NNLO) electroweak fermionic corrections refer to the contributions at order $\mathcal{O}(\alpha^2)$ in the Standard Model (SM) perturbative expansion to scattering amplitudes and cross sections, involving Feynman diagrams with at least one closed fermion loop in the electroweak sector. These corrections are a crucial component in achieving the precision required for theoretical predictions at future lepton colliders and for benchmark processes such as Drell–Yan production and Higgsstrahlung ($e^+e^- \to ZH$). The evaluation and numerical stability of these corrections require advanced semi-numerical two-loop techniques, control of ultraviolet (UV) and infrared (IR) divergences, and detailed renormalization-scheme studies.

1. Perturbative Expansion and Classification of Corrections

The perturbative expansion of the unpolarized cross section for processes such as $e^+e^- \to ZH$ or Drell–Yan ($e^+e^- \to \ell^+\ell^- / q\bar{q}$) is organized as

$$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$

where

• $\sigma^{(0)}$ is the tree-level (Born) cross section,
• $\sigma^{(1)}$ is the full one-loop (next-to-leading order, NLO) electroweak correction (composed of fermionic and bosonic contributions),
• $\sigma^{(2)}$ includes the NNLO contributions.

At the NNLO, “fermionic corrections” are defined as those with at least one closed fermion loop in the two-loop amplitude. These can be further subdivided according to the number of closed loops ($N_f = 1, 2$). For instance, in $e^+e^- \to ZH$, the decomposition is

$e^+e^- \to ZH$0

Typical topologies include reducible diagrams (self-energy insertions), two-loop vertices, box diagrams, and triangle sub-loops, all characterized by at least one (possibly two) closed fermion (lepton or quark) loop(s) (Freitas et al., 2022, Freitas et al., 2023, Freitas et al., 17 Dec 2025).

2. Semi-Numerical Two-Loop Evaluation Techniques

The computation of genuine two-loop electroweak diagrams with fermion loops is performed using semi-numerical strategies that combine analytical and numerical integration. The key steps are:

• Introduction of Feynman parameters for one sub-loop (e.g., $e^+e^- \to ZH$1-loop).
• Use of dispersion relations for the same sub-loop leading to an integral over a spectral parameter $e^+e^- \to ZH$2,
• Application of Passarino–Veltman reduction for the remaining sub-loop (e.g., $e^+e^- \to ZH$3), and
• Final numerical integration over the relevant Feynman and dispersion parameters.

For a generic two-loop vertex diagram, one encounters integrals of the form: $e^+e^- \to ZH$4 which are reduced, after a Feynman-parameter shift and dispersion for the $e^+e^- \to ZH$5-loop, to a well-controlled integral over a finite region. The method systematically isolates UV subdivergences and global UV divergences through analytic subtraction terms. Tensor reduction and vacuum-integral subtraction further yield fully finite quantities that are then computed numerically with 3–4 digit precision per diagram; in cases of strong cancellations, quadruple precision is employed (Freitas et al., 2022, Freitas et al., 17 Dec 2025).

3. UV and IR Subtraction and Renormalization Schemes

Two-loop fermionic corrections feature nested UV divergences, as well as potential soft-photon IR singularities, necessitating a robust subtraction strategy. Three levels of subtraction are implemented:

• Sub-loop UV subtraction: Achieved by dropping external momenta in the propagators of the divergent sub-loop and treating the divergence analytically, e.g., using $e^+e^- \to ZH$6 or $e^+e^- \to ZH$7 master integrals.
• Global UV subtraction: Remaining divergences are canceled by vacuum diagrams (with all external momenta set to zero).
• Infrared subtraction: IR divergences due to soft internal photons are removed by process-specific eikonal counterterms following the Yennie–Frautschi–Suura scheme, ensuring the two-loop amplitude is both UV- and IR-finite (Freitas et al., 17 Dec 2025).

Renormalization can be performed in the on-shell (OS) mass scheme for all fields, with two main input definitions for $e^+e^- \to ZH$8:

• The $e^+e^- \to ZH$9-scheme, using the Thomson limit value of $e^+e^- \to ZH$0, and requiring inclusion of $e^+e^- \to ZH$1,
• The $e^+e^- \to ZH$2 scheme, using Fermi constant $e^+e^- \to ZH$3 from muon decay, absorbing leading universal corrections ($e^+e^- \to ZH$4) in the input parameters (Freitas et al., 2023).

4. Numerical Results and Polarization Dependence

At $e^+e^- \to ZH$5 GeV for $e^+e^- \to ZH$6, the numerical evaluation in the $e^+e^- \to ZH$7 scheme finds: $e^+e^- \to ZH$8 resulting in a net shift at NNLO of $e^+e^- \to ZH$9 fb or approximately $e^+e^- \to \ell^+\ell^- / q\bar{q}$0 relative to the NLO result (Freitas et al., 2022).

For polarized beams:

• $e^+e^- \to \ell^+\ell^- / q\bar{q}$1: NNLO correction increases cross-section by $e^+e^- \to \ell^+\ell^- / q\bar{q}$2,
• $e^+e^- \to \ell^+\ell^- / q\bar{q}$3: NNLO correction decreases cross-section by $e^+e^- \to \ell^+\ell^- / q\bar{q}$4 (Freitas et al., 2023).

In neutral-current Drell–Yan-like processes ($e^+e^- \to \ell^+\ell^- / q\bar{q}$5 at $e^+e^- \to \ell^+\ell^- / q\bar{q}$6GeV): $e^+e^- \to \ell^+\ell^- / q\bar{q}$7 thus the NNLO fermionic correction is $e^+e^- \to \ell^+\ell^- / q\bar{q}$8 at this energy. For quark final states, effects are $e^+e^- \to \ell^+\ell^- / q\bar{q}$9 (Freitas et al., 17 Dec 2025). Modifications of differential distributions (e.g., $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$0) are of similar relative size, and the shape shows mild forward–backward asymmetries.

5. Renormalization-Scheme Dependence and Theory Uncertainties

The inclusion of NNLO fermionic corrections leads to a notable reduction of renormalization-scheme dependence. For $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$1, the scheme difference between $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$2 and $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$3:

• At LO: $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$4,
• At NLO: $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$5,
• At NNLO (fermionic): $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$6,
• With mixed QCD–EW and NNLO fermionic: $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$7 (Freitas et al., 2023).

Estimates of missing bosonic NNLO two-loop corrections indicate that the associated theory error is below $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$8, dominated by the square of the NLO bosonic amplitude. Thus, the full fermionic NNLO contribution is essential to reach sub-percent-level precision.

6. Phenomenological Impact and Future Directions

Given experimental precision targets of $$\sigma = \sigma^{(0)} + \alpha \sigma^{(1)} + \alpha^2 \sigma^{(2)} + \ldots$$9 (FCC-ee), $\sigma^{(0)}$0 (CEPC), and $\sigma^{(0)}$1 (ILC) for the $\sigma^{(0)}$2 cross section, the inclusion of NNLO electroweak fermionic corrections is indispensable to match the anticipated measurement uncertainty (Freitas et al., 2022, Freitas et al., 2023). For Drell–Yan and $\sigma^{(0)}$3 processes, these corrections are similarly essential for FCC-ee, ILC, CEPC physics programs and high-luminosity LHC analyses requiring per-mille accuracy (Freitas et al., 17 Dec 2025).

A plausible implication is that the methods and computational tools developed (dispersive semi-numerical approaches, robust subtraction frameworks) will form the cornerstone for future calculations of complete NNLO electroweak corrections, including the outstanding bosonic two-loop pieces and real emission. This will further reduce theory uncertainties and enable sensitivity to new-physics deviations at future colliders. Moreover, polarization-dependent effects at the per-mille to percent level may provide additional handles to disentangle possible new-physics contributions in electroweak observables (Freitas et al., 2023).