Towards the two-loop electroweak corrections to the Drell-Yan process: the complete fermionic contributions

This paper presents the complete set of ultraviolet-renormalized and infrared-subtracted two-loop virtual fermionic contributions to the Drell-Yan process ($u\bar u\to \mu^+\mu^-$) using an automated methodology, providing a crucial building block for next-to-next-to-leading order electroweak simulations.

The Gist

Imagine the universe as a giant, incredibly complex machine where tiny particles like quarks and electrons are constantly colliding and interacting. Physicists try to predict exactly what happens during these collisions using a set of rules called the "Standard Model." However, these rules aren't perfect; they are like a map that works well for a city but gets fuzzy when you try to zoom in on every single crack in the sidewalk. To get a truly accurate map, scientists have to calculate "corrections"—tiny adjustments that account for the chaotic quantum noise happening in the background.

This paper is about the team taking a massive step forward in drawing that ultra-precise map for a specific event: when a quark and an anti-quark smash together to create a pair of muons (heavy cousins of electrons). This event is known as the Drell-Yan process.

Here is the breakdown of what they did, using everyday analogies:

1. The Goal: The "Two-Loop" Challenge

Think of calculating a particle collision like trying to predict the weather.

• Level 1 (Tree Level): You look at the sky and say, "It's sunny." (This is the basic, simple prediction).
• Level 2 (One-Loop): You realize, "Oh, there are some clouds and a breeze." You add those details.
• Level 3 (Two-Loop): This is the level this paper tackles. It's like realizing the breeze is causing a cloud to swirl, which is creating a tiny rain shower that affects the temperature, which in turn changes the wind. It's a second layer of complexity.

The authors calculated the complete set of "fermionic" corrections for this two-loop level. In plain English, they tracked every possible way that a closed loop of matter particles (fermions) could wiggle in and out of existence during the collision and change the outcome. They didn't just guess; they calculated the entire set of these specific loops.

2. The Messy Parts: Cleaning Up the Math

When you try to do these calculations, the math often explodes into infinity. It's like trying to measure a room with a ruler that keeps stretching to infinity. To fix this, the team had to perform two major "clean-up" operations:

• UV Renormalization (The "Infinite" Fix): This is about removing the "ultraviolet" infinities. Imagine you are building a house, but your blueprint has a section where the walls are infinitely tall. You have to rewrite the blueprint to make the walls a sensible height without changing the actual shape of the house. The authors developed a rigorous method to cut out these infinities and replace them with real, measurable numbers (like the mass of the Z boson).
• IR Subtraction (The "Glitch" Fix): This is about removing "infrared" glitches. Imagine your camera lens has a smudge that makes the picture blurry. In particle physics, this blurriness comes from particles that are too soft to be detected but still mess up the math. The team created a "cleaning cloth" (mathematical subtraction) to wipe away these blurs so they could see the clear picture of the collision.

3. The "Chiral" Puzzle (The $\gamma_5$ Problem)

One of the biggest headaches in this field is a mathematical object called $\gamma_5$. Think of it as a special 4-dimensional gear in a machine. When the physicists try to run their calculations in a slightly different dimension (a mathematical trick called "dimensional regularization" used to handle the infinities), this gear doesn't fit right. It's like trying to put a square peg in a round hole.

There are different ways to force the gear to fit, but they all break a different rule of the machine. The authors used a specific strategy (the Kreimer scheme) that keeps the gears turning smoothly by accepting that you have to look at the machine from a specific angle (breaking "cyclicity") to make the math work. They proved that no matter which angle they looked from, the final result was the same.

4. The Automation: The "Robot Factory"

Calculating these diagrams by hand is impossible. There are thousands of them. The authors built a "robot factory" (automated computer code) that:

1. Generates all the possible diagrams (the blueprints).
2. Calculates the complex integrals (the measurements).
3. Applies the renormalization and subtraction rules (the clean-up crew).
4. Checks for errors (the quality control).

They tested this robot extensively to ensure it didn't make mistakes, verifying that the infinities canceled out perfectly and that the results were consistent.

5. The Result: A Sharper Lens

The paper presents the final, finite numbers that remain after all the cleaning and fixing. These numbers represent the "fermionic" contribution to the two-loop corrections for creating muon pairs.

Why does this matter?
The Large Hadron Collider (LHC) and future colliders are becoming incredibly precise. They can measure things with an accuracy of less than one part in a thousand. To match this precision, the theoretical predictions must be just as sharp. This paper provides a crucial "building block" for those predictions. Without these specific calculations, the theoretical map would be too blurry to compare with the high-definition photos taken by the experiments.

In summary: The authors built a highly sophisticated, automated mathematical machine to calculate the second-order "wiggles" of matter particles in a specific collision. They solved the problem of infinite numbers, fixed the tricky 4-dimensional gears, and delivered a clean, precise result that helps physicists understand the universe with unprecedented accuracy.

Technical Summary

Technical Summary: Towards the two-loop electroweak corrections to the Drell-Yan process: the complete fermionic contributions

Problem Statement
The precision of experimental measurements at the LHC (High-Luminosity phase) and future colliders like FCC-ee demands theoretical predictions for the Drell-Yan process ($u\bar{u} \to \mu^+\mu^-$) with sub-per mille accuracy. While Next-to-Leading Order (NLO) electroweak (EW) corrections and mixed QCD-EW corrections at Next-to-Next-to-Leading Order (NNLO) are available, the complete set of second-order (NNLO) pure EW virtual corrections remains a major challenge in Quantum Field Theory. Specifically, the evaluation of the complete subset of fermionic contributions—defined by the presence of closed fermionic loops in Feynman diagrams or counterterms—is required to reduce theoretical uncertainties below the current percent level in hadron-hadron collisions and match expected experimental precision. This calculation involves complex theoretical issues, including UV renormalisation with unstable particles, IR subtraction, and the handling of chiral couplings ($\gamma_5$) in dimensional regularisation.

Methodology
The authors present an automated, semi-analytical framework to compute the complete set of second-order fermionic EW virtual corrections. The methodology is structured as follows:

1. Renormalisation Scheme: The calculation employs the Complex Mass Scheme (CMS) for gauge boson masses to preserve gauge invariance in the presence of unstable particles, combined with the $\overline{\text{MS}}$ scheme for the electromagnetic coupling. The authors utilize the Background Field Gauge (BFG) approach, which ensures QED-like Ward identities for Green's functions involving background photon fields.
2. UV Renormalisation: A comprehensive set of counterterms is derived up to $\mathcal{O}(\alpha^2)$. This includes mass, wavefunction, and coupling counterterms. The authors explicitly address the "sub-renormalisation" diagrams (insertion of one-loop counterterms into one-loop diagrams) and the products of first-order counterterms arising from the expansion of renormalised vertex functions. Tadpole contributions are handled via the Parameter-Renormalised Tadpole Scheme (PRTS), ensuring the vanishing of the Higgs tadpole at the renormalised level.
3. Infrared (IR) Subtraction: The authors define a finite remainder by subtracting IR divergences from the UV-renormalised amplitude. The subtraction terms are derived from the abelianisation of QCD soft and cusp anomalous dimensions, adapted for the specific fermionic contributions. The subtraction operator $I^{(0,2)}_{fer}$ accounts for the emission of photons with fermionic corrections and tree-level fermion pair emissions.
4. $\gamma_5$ Prescription: To handle chiral couplings in dimensional regularisation, the authors adopt the Kreimer scheme (implemented in the Abiss package). This approach preserves the anti-commutativity of $\gamma_5$ with Dirac matrices and the trace condition involving the Levi-Civita tensor, while abandoning trace cyclicity. A specific "reading point" is introduced for each distinct trace to manage the non-cyclic nature.
5. Integral Reduction and Evaluation: The two-loop Feynman integrals are reduced to Master Integrals (MIs) using the KIRA and FireFly frameworks within the Oceann automation suite. Numerical evaluation of these MIs is performed using AMFlow and SeaSyde. To mitigate computational complexity, the muon mass is set to zero in diagrams where it does not generate final-state collinear singularities, retaining it only where necessary to regulate such divergences.

Key Contributions

• Complete Fermionic Subset: The paper provides the first evaluation of the complete set of second-order fermionic EW virtual corrections to the Drell-Yan process ($u\bar{u} \to \mu^+\mu^-$) for arbitrary kinematics. This subset is defined by the presence of at least one closed fermionic loop in the diagrams or counterterms.
• Automated Workflow: The authors demonstrate a fully automated workflow for the calculation of two-loop amplitudes, integrating UV renormalisation, IR subtraction, and the specific treatment of $\gamma_5$.
• Validation of Tools: The work serves as a benchmark for the Oceann, Abiss, AMFlow, and SeaSyde toolchains. The authors explicitly verify the cancellation of UV and IR poles and the independence of the final results from the specific $\gamma_5$ prescription (reading point and fixed-point conventions) used during intermediate steps.
• Ward Identity Verification: The calculation verifies the validity of BFG Ward identities for the renormalised Green's functions, ensuring the gauge invariance of the amplitude.

Results

• Finite Amplitude: The authors present the UV-renormalised and IR-subtracted finite contribution to the scattering amplitude. Numerical results for the coefficients of the Laurent expansion in $\epsilon$ (where $d=4-2\epsilon$) are provided for a specific kinematic point ($s=10000$ GeV$^2$, $t=-2500$ GeV$^2$).
• Pole Cancellation: The sum of the diagrammatic part, the UV counterterms, and the IR subtraction terms results in the cancellation of all $1/\epsilon^n$ poles to a relative precision of $10^{-7}$. The remaining imperfections are attributed to the approximation of the muon mass being set to zero in specific diagrams, which introduces terms proportional to $m_\mu^2/\mu_W^2$.
• Structure of Divergences: The analysis confirms that UV divergences cancel within specific subsets of diagrams (e.g., vertex and wavefunction corrections satisfying QED-like Ward identities) and that the cancellation of IR divergences relies on the correct combination of virtual and real-emission contributions.

Significance and Claims
The authors position this work as a critical building block toward the complete evaluation of second-order EW corrections to the Drell-Yan process. They claim that:

1. The calculation addresses major conceptual and technical challenges, including the renormalisation of unstable particles, the systematic subtraction of IR divergences, and the consistent treatment of chiral couplings.
2. The gauge invariance of the results is enforced through the CMS and the verification of Ward identities.
3. The universality of the derived counterterms and subtraction operators allows these ingredients to be applied to other scattering processes.
4. While this paper establishes the virtual fermionic corrections, the authors note that the combination of these results with real-emission contributions and the treatment of perturbative unitarity violations inherent in the CMS (regarding the complex mass of stable particles) will be addressed in future publications.

The work does not claim to provide the final NNLO EW cross-section prediction for the LHC but rather delivers the necessary virtual component and validates the computational infrastructure required to achieve such predictions.