Heavy-quark box-loop corrections to $q\bar q \to Zγ$ at two loops in QCD

This paper presents a numerical computation of two-loop QCD corrections to $Z\gamma$ production at the LHC, incorporating both light- and heavy-quark box-loop contributions via a Monte Carlo integration method that simultaneously handles loop and phase space integrations while validating results against known benchmarks.

The Gist

Imagine you are trying to predict exactly how a specific type of traffic jam will form on a massive, multi-lane highway (the Large Hadron Collider, or LHC). You know the basic rules of the road, but you want to know the exact probability of a crash happening when two cars (quarks) collide and produce a specific pair of vehicles: a Z-boson (a heavy, mysterious truck) and a photon (a flash of light).

This paper is about calculating the "traffic rules" for these collisions with extreme precision, specifically looking at the tiny, invisible detours that heavy particles take during the crash.

Here is the breakdown of what the authors did, using simple analogies:

1. The Problem: The "Ghost" Detours

In the quantum world, when particles collide, they don't just bounce off each other. They briefly create "ghost" particles that pop in and out of existence before vanishing. These are called loops.

• The Light Quarks: Imagine these ghosts are like tiny, fast ants. They are everywhere, but they are so light they barely change the traffic flow.
• The Heavy Quarks: Now imagine these ghosts are like giant, heavy elephants. They are rare, but when they appear in the loop, they weigh down the whole system and change the outcome significantly.

For a long time, physicists could calculate the "ant" loops perfectly. But calculating the "elephant" loops (heavy quarks like the Top and Bottom quarks) in a two-step collision (two loops) was like trying to solve a Rubik's cube while blindfolded and juggling. It was too complex for standard math tools.

2. The Solution: A New "Traffic Simulation"

The authors developed a new numerical method (a computer simulation) to handle this complexity. Instead of trying to solve the whole equation at once (which is impossible), they broke it down:

• Local Subtraction (The "Traffic Cop"): In their simulation, certain parts of the calculation blow up to infinity (like a traffic jam that never ends). The authors invented a way to put up "Traffic Cop" signs at the exact spot where the jam starts. These signs cancel out the infinite parts locally, right where they happen, leaving behind a manageable, finite number.
• The "Loop-Tree" Transformation: They used a clever trick called "Loop-Tree Duality." Imagine taking a tangled ball of yarn (the complex loop) and cutting it open to lay it flat like a tree branch. This makes it much easier to count and measure without getting lost in the knots.
• Monte Carlo Integration (The "Dice Roll"): Since the math is still too messy to solve with a pen and paper, they used a method called Monte Carlo. Imagine throwing millions of darts at a target to estimate its shape. Their computer threw "darts" (random numbers) at the collision scenario billions of times to find the average result.

3. What They Found

They ran their simulation for three scenarios:

1. Massless Quarks (The Ants): They checked their work against known answers. Their simulation matched perfectly, proving their "Traffic Cop" and "Yarn" methods work.
2. Bottom Quarks (The Medium Elephants): They calculated the effect of these heavier particles. The results were slightly different from the "ants," showing that even medium-sized ghosts matter.
3. Top Quarks (The Giant Elephants): This was the big discovery. Because the Top quark is so heavy, it changes the physics entirely. The authors found that for these heavy particles, certain "traffic jams" (singularities) simply disappear because the elephant is too heavy to fit through the narrow door. This resulted in a completely different numerical outcome.

4. Why It Matters

The authors didn't just stop at the math; they applied it to the real world. They combined their results with data about how protons are made (Parton Distribution Functions) to predict what will happen in actual collisions at the LHC.

The Bottom Line:
Think of this paper as upgrading the GPS for the world's most powerful particle collider.

• Before: The GPS could tell you the route for light cars (massless quarks) but got confused by heavy trucks (heavy quarks).
• Now: The authors have built a new GPS that can handle both light cars and heavy elephants, even when they are driving in complex, multi-lane loops.

This allows scientists to predict the "Z-boson + Photon" signal with much higher accuracy. If the real data from the LHC matches their new, heavy-quark-inclusive prediction, it confirms our understanding of the universe. If it doesn't match, it might mean there is a new, undiscovered particle hiding in the traffic!

In short: They built a super-accurate calculator that can finally weigh the "elephants" in the quantum world, ensuring our predictions for future experiments are rock solid.

Technical Summary

Here is a detailed technical summary of the paper "Heavy-quark box-loop corrections to $\bar{q}q \to Z\gamma$ at two loops in QCD" by Dario Kermanschah and Matilde Vicini.

1. Problem Statement

The computation of two-loop multi-leg scattering amplitudes involving multiple mass scales remains a significant challenge in Quantum Chromodynamics (QCD). Specifically, the authors address the calculation of two-loop QCD corrections to the process $q\bar{q} \to Z\gamma$ (quark-antiquark annihilation into a Z boson and a photon) at the Large Hadron Collider (LHC).

The primary difficulties addressed are:

• Mass Scales: The process involves massless incoming quarks and photons, but a massive $Z$ boson and potentially massive quarks (bottom and top) circulating in the loop.
• Singularities: The calculation must handle infrared (IR) and ultraviolet (UV) divergences, as well as threshold singularities arising from on-shell intermediate states in the loop integrals.
• Lack of Benchmarks: While massless quark loop contributions have analytic solutions, there were no available analytic benchmarks for the heavy-quark box-loop contributions (specifically top and bottom quarks) in this specific process.

2. Methodology

The authors employ a fully numerical approach based on the pipeline developed in their previous works [1, 2]. The method avoids analytic integration where possible, instead performing numerical integration over spatial loop momenta.

Key Technical Components:

• Loop-Tree Duality & Local Subtraction:
  • The calculation constructs finite two-loop amplitudes directly in loop momentum space.
  • IR and UV Subtraction: Local counterterms are introduced to render individual Feynman diagrams finite in four dimensions before integration. This involves subtracting collinear IR limits and applying UV renormalization via the $R$-operation.
  • Threshold Subtraction: Threshold singularities (where internal particles go on-shell) are identified using Cutkosky cuts. These are subtracted locally using a counterterm method based on the $i\epsilon$ causal prescription, avoiding the need for complex contour deformation into the complex plane.
• Loop-Energy Integration:
  • The energy components of the loop momenta are integrated out analytically using a custom FORM script.
  • This utilizes partial fraction decomposition to isolate simple poles, effectively converting the loop integrals into representations equivalent to partially time-ordered perturbation theory or causal loop-tree duality.
• Numerical Integration:
  • The remaining integrals over spatial loop momenta and phase space are performed using a multi-channel Monte Carlo method with adaptive importance sampling (VEGAS algorithm).
  • The pipeline handles the simultaneous integration of loop momenta and parton distribution functions (PDFs) for proton-proton collisions.
• Massive vs. Massless Handling:
  • The framework treats the $Z$ boson as an off-shell photon ($\gamma^*$) for the kinematic structure, adjusting only for vector couplings.
  • Heavy quark masses ($m_b, m_t$) are explicitly included in the propagators. The code dynamically activates threshold counterterms based on the kinematic region (e.g., whether $\hat{s} < 4m_q^2$).

3. Key Contributions

1. First Numerical Calculation of Heavy-Quark Box Loops: The paper presents the first numerical results for the two-loop squared matrix elements of $q\bar{q} \to Z\gamma$ mediated by bottom and top quark box loops.
2. Validation: The massless quark loop results were validated against known analytic results from Ref. [3], confirming the accuracy of the numerical pipeline.
3. Extension to Massive Scales: The work demonstrates the flexibility of the numerical approach in handling processes with multiple mass scales (massless external legs, massive $Z$, and massive internal quarks) without requiring new analytic derivations.
4. Phenomenological Application: The authors performed the full convolution with Parton Distribution Functions (PDFs) to compute the double-virtual corrections to the $pp \to Z\gamma$ cross-section at $\sqrt{s} = 13$ TeV.

4. Results

The authors evaluated the helicity-summed squared matrix elements and the integrated cross-section corrections for three distinct cases: massless quarks, bottom quarks ($m_b = 4.18$ GeV), and top quarks ($m_t = 172$ GeV).

• Massless Validation: The numerical results for massless quarks matched the analytic benchmarks within statistical errors (relative errors $\sim 0.1\% - 0.3\%$).
• Heavy-Quark Effects:
  • Bottom Quarks: Due to the small mass, the bottom-quark contribution is numerically close to the massless result, with differences primarily visible in non-planar diagrams (percent-level effects).
  • Top Quarks: The large top-quark mass significantly alters the analytic structure. Specifically, final-state threshold singularities are absent (since $\hat{s} < 4m_t^2$ in the relevant kinematic regions), leaving only $s$-channel thresholds. This leads to distinct numerical values compared to the massless case.
• Cross-Section Corrections (Table 3):
  • The double-virtual corrections were computed with high precision (targeting $<1\%$ total error, though individual components required higher precision due to cancellations).
  • Massless: $\sigma^{(2)} \approx 1.70 \times 10^{-5}$ (arbitrary units).
  • Bottom: $\sigma^{(2)} \approx 3.95 \times 10^{-5}$.
  • Top: $\sigma^{(2)} \approx 8.88 \times 10^{-6}$.
  • The results show that heavy-quark loops provide non-negligible corrections to the $Z\gamma$ production cross-section.
• Performance: The evaluation times per sample (Table 4) are in the millisecond range, demonstrating the efficiency of the optimized integrand expressions (using Spenso and Symbolica).

5. Significance

• Phenomenological Precision: These calculations are crucial for achieving Next-to-Next-to-Leading Order (NNLO) precision in $Z\gamma$ production at the LHC. Combining these double-virtual corrections with real-emission contributions (available in other literature) will allow for highly accurate theoretical predictions to match experimental data.
• Methodological Robustness: The paper validates a purely numerical framework capable of handling complex, multi-scale two-loop amplitudes where analytic solutions are either unknown or intractable. This paves the way for computing similar corrections for other processes involving heavy quarks (e.g., Higgs production, top-quark pair production).
• Heavy-Quark Physics: By providing the first benchmarks for heavy-quark box contributions in this channel, the work fills a gap in the Standard Model precision tests, particularly regarding the role of the top quark in electroweak processes.

In conclusion, Kermanschah and Vicini have successfully extended a state-of-the-art numerical integration pipeline to compute previously unavailable two-loop heavy-quark corrections for $Z\gamma$ production, demonstrating both the feasibility of the method and its necessity for high-precision collider phenomenology.