NNLO QCD corrections to $\gamma \gamma \rightarrow Q\bar{Q}$ from Local Unitarity combined with Coulomb resummation and NLO EW effects

This paper presents state-of-the-art NNLO QCD predictions for heavy-quark pair production in direct photon fusion by applying the Local Unitarity formalism to handle infrared singularities and combining these results with NLO electroweak corrections and NLP Coulomb resummation for top, bottom, and charm quarks.

The Gist

Imagine you are trying to predict exactly how much "stuff" (heavy particles like top, bottom, or charm quarks) will be created when two beams of light (photons) smash into each other. This isn't just a simple collision; it's a chaotic dance of quantum mechanics where particles pop in and out of existence, and the math to describe it is notoriously difficult.

This paper is like a master chef presenting a brand-new, ultra-precise recipe for calculating this collision, specifically for the most complex version of the recipe yet: Next-to-Next-to-Leading Order (NNLO).

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

1. The Problem: The "Infinite Noise" of Quantum Math

In the world of particle physics, when you try to calculate what happens in a collision, your equations often blow up. They produce "infinities" (singularities) because particles can be infinitely soft or travel in perfectly straight lines.

• The Old Way: Traditionally, physicists had to do a two-step dance. First, they would calculate the "virtual" particles (loops in the math) separately from the "real" particles flying out. Then, they would manually subtract the infinities from one to cancel out the infinities in the other. It was like trying to balance a checkbook by calculating income and expenses on separate pieces of paper and hoping they match up at the end.
• The New Way (Local Unitarity): The authors used a new method called Local Unitarity (LU). Think of this as doing the calculation in a single, unified room where the income and expenses are balanced instantly as you write them down. The "noise" (infinities) cancels out locally, right where it happens, before you even finish the calculation. This allows them to use a computer to simulate the collision directly, without needing to solve impossible algebraic puzzles first.

2. The Experiment: Smashing Light into Light

The paper focuses on Photon-Photon collisions ($\gamma\gamma \to Q\bar{Q}$).

• The Setup: Imagine two flashlights shining at each other. In the real world, these photons come from huge machines like the Large Hadron Collider (LHC) or future electron-positron colliders.
• The Target: They are looking for the creation of heavy quarks (the building blocks of protons and neutrons). They studied three types:
  • Top Quarks: The "giants" of the quark world (very heavy).
  • Bottom Quarks: The "medium" sized ones.
  • Charm Quarks: The "light" ones (but still heavy compared to the energy scale).

3. The Challenge: The "Coulomb Glue"

When heavy particles are created, they move slowly at first. Because they are charged and moving slowly, they stick together like magnets before flying apart. This is called the Coulomb effect.

• The Analogy: Imagine two magnets snapping together. If you try to calculate their energy just by looking at them moving fast, you get it wrong. You have to account for the "glue" pulling them together.
• The Solution: The authors didn't just calculate the collision; they also "resummed" (added up) these glue effects. For the heavy Top quark, this glue is a small detail. But for the lighter Bottom and Charm quarks, this glue is a massive force that changes the result significantly. If you ignore it, your prediction is wrong.

4. The Results: A New Standard of Precision

The authors combined their new "Local Unitarity" math with the "Coulomb glue" correction and some Electroweak (electromagnetic) tweaks.

• Top Quarks: The new calculation shows that previous estimates were slightly off. The new "NNLO" correction adds about 6% to the predicted amount of top quarks. It's a small but crucial adjustment for future experiments.
• Bottom & Charm Quarks: Here, the story is wilder. The "glue" (Coulomb effect) actually cancels out a lot of the extra particles predicted by the new math. It's like adding a new ingredient to a cake that makes it rise, but then adding a second ingredient that makes it sink, resulting in a perfect balance.
• The "Discrepancy" Mystery: For years, experiments at the old LEP collider measured more Bottom quarks than the theory predicted. This paper provides the most precise theory yet. While it doesn't fully solve the mystery (the data is still a bit higher than the theory), it narrows the gap and gives scientists a much sharper tool to investigate why the difference exists.

5. The Toolkit: PHIQUE

The authors didn't just write a paper; they built a public software tool called PHIQUE.

• The Metaphor: Think of this as releasing a high-end, open-source calculator app. Before this, only a few experts with supercomputers could do these complex calculations. Now, any physicist can download PHIQUE, type in "I want to know the collision rate at the LHC," and get a precise answer that includes all the latest math corrections.

Summary

This paper is a triumph of computational physics.

1. New Method: They used a revolutionary way to handle the math (Local Unitarity) that avoids the usual headaches of infinities.
2. High Precision: They calculated the collision of light creating heavy matter with the highest precision ever achieved (NNLO).
3. Real World Application: They applied this to real-world scenarios at the LHC and future colliders, providing the "gold standard" predictions that experimentalists need to know if they are seeing new physics or just standard behavior.

In short, they built a better microscope to see how light turns into matter, revealing that the "glue" holding particles together plays a much bigger role than we thought, especially for the lighter heavy particles.

Technical Summary

1. Problem Statement

The paper addresses the calculation of heavy-quark pair production ($\gamma\gamma \to Q\bar{Q}$, where $Q = c, b, t$) in photon-photon collisions. While leading-order (LO) and next-to-leading-order (NLO) calculations are well-established, the Next-to-Next-to-Leading Order (NNLO) QCD corrections for this process had remained unknown.

Key challenges in this calculation include:

• Complexity of NNLO: Traditional approaches require separating loop (virtual) and phase-space (real emission) integrals, necessitating complex infrared (IR) subtraction schemes and dimensional regularization. For processes involving massive colored particles, this often involves calculating two-loop amplitudes with complicated special functions (e.g., elliptic integrals).
• Threshold Singularities: Near the production threshold ($s \approx 4m_Q^2$), the cross-section is plagued by Coulomb singularities ($\alpha_s/\beta_Q$) that require resummation to obtain reliable predictions.
• Resolved Photon Ambiguity: At NNLO, the distinction between "direct" photon interactions and "resolved" photon interactions (where the photon fluctuates into hadrons) becomes scheme-dependent. The authors focus on the direct channel ($\gamma\gamma \to Q\bar{Q}$), which is unambiguously defined if interactions between photons and massless quarks are ignored.

2. Methodology: Local Unitarity (LU)

The core innovation of this work is the application of the Local Unitarity (LU) formalism to a previously unknown NNLO process.

• Integrand-Level Cancellation: Unlike traditional methods that integrate separately and subtract divergences, LU combines loop and phase-space integrals at the integrand level. It utilizes the Kinoshita–Lee–Nauenberg (KLN) theorem to ensure that final-state IR singularities cancel locally before integration.
• Forward-Scattering (FS) Diagrams: The cross-section is expressed as a sum over Forward-Scattering diagrams (obtained by gluing the amplitude and its conjugate). The method applies Cutkosky cutting rules to these diagrams.
• Loop-Tree Duality (LTD): The method employs LTD to express loop integrals in terms of a common measure over spatial loop momenta.
• Causal Flow: A "causal flow" is constructed to parametrize the phase space, correlating IR singularities and allowing for the local cancellation of thresholds without dimensional regularization.
• Renormalization: UV singularities are subtracted locally using the BPHZ formalism.
• Handling Thresholds: For the specific process $\gamma\gamma \to Q\bar{Q}$, the authors demonstrate that for the non-singlet contributions (Classes A, B, D), the LU integrand is naturally free of threshold singularities. For the singlet contributions (Class C), which involve double-gluon cuts, they employ a subtraction strategy using a separate calculation of $\gamma\gamma \to gg$ to isolate the physical contribution.
• Coulomb Resummation: The fixed-order NNLO results are combined with Non-Relativistic QCD (NRQCD) and potential NRQCD (pNRQCD) to resum Coulomb gluon exchanges up to Next-to-Leading Power (NLP) in the threshold region.
• Electroweak (EW) Effects: NLO EW corrections are computed separately using the MadGraph5_aMC@NLO framework and combined with the QCD results.

3. Key Contributions

• First NNLO Calculation: This is the first computation of NNLO QCD corrections for direct heavy-quark pair production in photon fusion.
• 138 Diagrams: The authors enumerated and computed contributions from 138 distinct forward-scattering diagrams, categorized into four gauge-invariant classes (A, B, C, D).
• Proof of Feasibility: The work serves as a major proof-of-principle for the LU formalism, demonstrating its ability to handle complex, semi-inclusive processes with massive particles at NNLO without analytic amplitude reduction.
• PHIQUE Code: The authors released a public code, PHIQUE (PHoton-Induced heavy QUark pair procEss), which integrates the LU results, performs Coulomb resummation, and convolves with photon fluxes for various collider scenarios.
• Comprehensive Phenomenology: The study provides state-of-the-art predictions for top, bottom, and charm quark production in:
  • Ultraperipheral Collisions (UPCs) at the LHC (p-p, p-Pb, Pb-Pb) and future FCC-hh.
  • $e^+e^-$ colliders (FCC-ee, CEPC, ILC, CLIC).

4. Key Results

The paper presents total cross-sections and K-factors (ratio of higher-order to LO) for $Q=t, b, c$.

• Top Quark ($t\bar{t}$):

  • Convergence: Perturbative convergence is excellent. NLO QCD increases the LO cross-section by ~21% (LHC), while NNLO adds another ~6%.
  • Coulomb Effects: Coulomb resummation effects are marginal (~1% reduction) because the top quark is produced with high velocity ($\beta_t$) in most kinematic regions.
  • EW vs. QCD: NLO EW corrections reduce the cross-section by ~5.5%, partially canceling the NNLO QCD enhancement.
  • Observability: Observation of $\gamma\gamma \to t\bar{t}$ is expected at the High-Luminosity LHC (HL-LHC) and future colliders.

• Bottom Quark ($b\bar{b}$):

  • Large Corrections: Due to the lower mass scale, QCD corrections are larger. NLO increases the cross-section by ~35-40%, and NNLO adds another ~17-23%.
  • Coulomb Resummation: The effect is significant (~10-14% reduction), partially canceling the NNLO QCD enhancement.
  • Tension with Data: The results revisit the historical discrepancy between LEP measurements and NLO predictions. The new NNLO+Resummation predictions are crucial for resolving this tension.

• Charm Quark ($c\bar{c}$):

  • Perturbative Instability: Due to the charm mass being close to $\Lambda_{QCD}$, the perturbative series converges slowly. NNLO corrections increase the cross-section by >40%, and scale uncertainties actually increase at NNLO compared to NLO.
  • Coulomb Cancellation: Coulomb resummation reduces the cross-section by ~40%, nearly compensating for the NNLO increase, thereby stabilizing the prediction.
  • Comparison with OPAL: The combined NNLO QCD + NLP Coulomb + NLO EW prediction agrees with OPAL data at LEP within one standard deviation, whereas LO and NLO alone show larger deviations.

• Scale Uncertainties:

  • For top quarks, scale uncertainties drop from $\pm 2\%$ (NLO) to $\pm 1\%$ (NNLO) and sub-percent with resummation.
  • For bottom and charm, uncertainties remain larger, highlighting the need for resummation and potentially higher-order calculations.

5. Significance

• Theoretical Advancement: The paper validates the Local Unitarity approach as a viable, powerful alternative to traditional semi-analytical NNLO methods, capable of handling massive particles and complex topologies directly in 4D.
• Phenomenological Impact: It provides the most precise theoretical predictions for heavy-flavor production in photon-photon collisions, essential for:
  • Testing the universality of heavy-quark fragmentation.
  • Probing color-singlet exchange mechanisms.
  • Searching for Beyond the Standard Model (BSM) physics (e.g., anomalous couplings, extra dimensions).
  • Reducing uncertainties in background estimates for other processes.
• Future Readiness: The results and the PHIQUE code are prepared for the high-luminosity phase of the LHC and future lepton colliders (FCC-ee, CEPC), where these processes will be measured with high precision.

In summary, this work represents a milestone in high-energy physics calculations, successfully applying a novel numerical integration technique to a difficult NNLO problem and delivering precise, resummed predictions for heavy-quark production across a wide range of collider energies.