Complete NLO BFKL impact factors for quarkonium hadroproduction in NRQCD: the case of ${}^1S_0^{[1]}$, ${}^1S_0^{[8]}$, and ${}^3S_1^{[8]}$ states

This paper presents the first complete next-to-leading-order calculation of BFKL impact factors for the hadroproduction of specific NRQCD quarkonium states by combining virtual corrections with real-emission contributions, demonstrating the cancellation of soft divergences and the compatibility of surviving collinear singularities with factorization to enable future next-to-leading-logarithmic precision studies.

The Gist

Imagine you are trying to predict the weather, but instead of clouds and rain, you are predicting how tiny, heavy particles called quarkonia (think of them as "heavy atoms" made of a charm or bottom quark and its anti-particle) are created when two protons smash into each other at the fastest speeds in the universe (like at the Large Hadron Collider).

This paper is a massive breakthrough in the "mathematical weather forecast" for these particles. Here is the story of what they did, explained simply.

1. The Problem: The "High-Energy Traffic Jam"

When protons collide at incredibly high speeds, they don't just bounce off; they create a chaotic storm of other particles. Physicists use a set of rules called QCD (Quantum Chromodynamics) to predict what happens.

Usually, when things move slowly, the math is easy. But when they move near the speed of light, the math gets messy because of "logarithmic corrections." Imagine trying to count the number of cars in a traffic jam where every car is honking, and the honks are getting louder and louder the faster you go. If you ignore the honks, your prediction is wrong.

To fix this, physicists use a special tool called BFKL. Think of BFKL as a "super-calculator" that sums up all those annoying honks (logarithms) to give an accurate prediction for high-speed collisions.

2. The Missing Piece: The "Impact Factor"

In the BFKL calculator, there are two main parts:

1. The Green's Function: This is the "highway" where the particles travel and interact.
2. The Impact Factor: This is the "entrance ramp." It describes exactly how a particle enters the highway from the proton.

For a long time, physicists had a very rough sketch of the entrance ramp (Leading Order). They knew the general shape, but they didn't know the details. To get a truly accurate prediction, they needed the Next-to-Leading Order (NLO) version. This is like upgrading from a hand-drawn sketch to a 3D blueprint with every pothole and speed bump accounted for.

The Problem: Until now, no one had finished the blueprint for the entrance ramp specifically for quarkonia (the heavy atoms). The math was too hard, and the calculations were incomplete.

3. The Solution: Finishing the Blueprint

This paper by Fucilla, Lansberg, Nefedov, Szymanowski, and Wallon is the first time anyone has completed the NLO blueprint for three specific types of quarkonium states.

They did this by tackling two types of "traffic" that happen during the collision:

• Virtual Corrections (The Ghosts): These are particles that pop in and out of existence for a split second, changing the energy slightly.
• Real Emissions (The New Cars): These are actual new particles (gluons) that get shot out during the crash.

The Magic Trick:
In physics, these two types of traffic often create "infinities" (mathematical errors where numbers go to infinity). It's like a calculator saying "Error: Divide by Zero."

• The authors calculated the "Ghost" traffic and found it created a huge infinity.
• They calculated the "New Car" traffic and found it created an opposite infinity.
• The Result: When they added them together, the infinities canceled each other out perfectly! This proved their math was consistent and the "blueprint" was solid.

4. Why This Matters: The "Color" of the Car

The paper also had to deal with a tricky concept called "Color." In this world, particles have a "charge" called color (Red, Green, Blue).

• Color Singlet: The particle is "color-neutral" (like a white car). It's stable and easy to handle.
• Color Octet: The particle is "color-charged" (like a neon car). It's messy and interacts violently with the environment.

The authors had to figure out how to build the entrance ramp for both the "white cars" and the "neon cars." They discovered that the "neon cars" (Octet states) had an extra, sneaky source of chaos (a soft divergence) that they had to specifically subtract out to make the math work. They successfully did this, creating a complete guide for both types.

5. The Big Picture: What Can We Do Now?

Before this paper, if you wanted to predict how often these heavy atoms are made in a specific direction (forward or backward) at the LHC, you had to guess or use approximations.

Now, thanks to this paper:

• Precision: Physicists can now make "Next-to-Leading Logarithmic" predictions. This is the difference between guessing the weather and having a satellite image with a 99% accuracy forecast.
• New Experiments: This allows scientists to study "associated production" (creating a quarkonium atom and a jet of particles at the same time, far apart from each other). This is a perfect test for the BFKL framework.
• Future Tech: It paves the way for understanding the "gluon saturation" (a state where protons are so packed with particles they act like a single fluid), which is crucial for understanding the early universe.

Summary Analogy

Imagine you are building a bridge (the BFKL framework) to cross a river of high-energy physics.

• Previous work: You had the pillars (the Green's function) and a rough idea of the road surface (Leading Order Impact Factors).
• This paper: The authors poured the concrete, installed the guardrails, and fixed the potholes for the entrance ramps (Impact Factors) for three specific types of vehicles. They even figured out how to handle the "neon cars" that vibrate the bridge.

Now, the bridge is safe to drive on, and scientists can finally cross over to the other side to discover new physics about how the universe works at its most fundamental level.

Technical Summary

Here is a detailed technical summary of the paper "Complete NLO BFKL impact factors for quarkonium hadroproduction in NRQCD: the case of $1S_0^{[1]}$,$1S_0^{[8]}$, and$3S_1^{[8]}$ states."

1. Problem Statement and Motivation

The paper addresses a critical gap in the theoretical description of heavy quarkonium production at high-energy hadron colliders (like the LHC).

• Context: In the high-energy limit ($s \gg Q^2 \gg \Lambda_{QCD}^2$), large logarithmic corrections of the form $\alpha_s^n \ln^n(s/Q^2)$ become significant. The Balitsky-Fadin-Kuraev-Lipatov (BFKL) framework resums these logarithms to all orders.
• Current Limitation: While Next-to-Leading Logarithmic (NLL) accuracy has been achieved for processes like Mueller-Navelet dijets, quarkonium production within the BFKL formalism has been restricted to Leading Logarithmic (LL) accuracy or specific approximations (e.g., fragmentation limits).
• Specific Challenge: The authors aim to compute the Next-to-Leading Order (NLO) corrections to the Impact Factors (IFs) for quarkonium production. Impact factors are the process-dependent coefficients in the BFKL cross-section formula that describe the coupling of the BFKL Pomeron (reggeized gluons) to the produced particles.
• Target States: The study focuses on the S-wave NRQCD states:
  • Color-Singlet (CS): $1S_0^{[1]}$(pseudoscalar, e.g.,$\eta_c, \eta_b$).
  • Color-Octet (CO): $1S_0^{[8]}$and$3S_1^{[8]}$(vector, e.g.,$J/\psi, \Upsilon$).
• Goal: To provide the first complete NLO calculation of these impact factors, including both virtual and real-emission corrections, thereby enabling NLL-precision phenomenology for forward-backward quarkonium production.

2. Methodology

The authors employ a sophisticated combination of theoretical frameworks and computational techniques:

• Factorization Frameworks:

  • BFKL Formalism: Used to resum high-energy logarithms. The cross-section is expressed as a convolution of two impact factors and a BFKL Green's function.
  • NRQCD Factorization: Used to separate the short-distance production of a heavy quark-antiquark pair ($Q\bar{Q}$) from the long-distance hadronization into a quarkonium state. The calculation is performed for specific $Q\bar{Q}$ intermediate states ($m$) with defined quantum numbers ($^{2S+1}L_J^{[1,8]}$).
  • Gauge-Invariant EFT: The authors utilize the Effective Field Theory (EFT) for Multi-Regge processes in QCD (Lipatov's EFT) to handle the off-shell reggeized gluons and ensure gauge invariance.

• Calculation Strategy:

  • Virtual Corrections: The authors build upon their previous work (JHEP 12 (2024) 129) which computed the one-loop virtual corrections.
  • Real Emission: They compute the tree-level squared matrix elements for the process $R + a \to Q\bar{Q}[m] + a$ (where $R$ is a reggeized gluon and $a$ is a parton).
  • Subtraction Method: To handle infrared (IR), collinear, and soft divergences in the real-emission contributions, they use a subtraction approach. They construct a subtraction term $J$ that mimics the singular behavior of the exact matrix element in the soft/collinear limits. This allows the finite part to be integrated numerically while the divergent parts are extracted analytically.
  • Divergence Cancellation: The soft and collinear divergences from the real emission are shown to cancel against the corresponding poles ($1/\epsilon, 1/\epsilon^2$) in the virtual corrections and the BFKL counterterms.
  • Rapidity Factorization: A crucial step involves defining a rapidity-factorization scheme to separate the impact factor from the BFKL Green's function. The authors derive the necessary counterterms to transform the result from the EFT scheme to the standard BFKL scheme.

• Technical Tools:

  • Dimensional Regularization ($D = 4-2\epsilon$) for UV and IR divergences.
  • Integration-by-Parts (IBP) identities to reduce loop integrals to master integrals.
  • Symbolic computation tools (FeynArts, FeynCalc, FORM) for generating and manipulating Feynman diagrams and algebraic expressions.

3. Key Contributions

The paper makes several distinct technical breakthroughs:

1. First Complete NLO Impact Factors: This is the first calculation of complete NLO impact factors for quarkonium production within the BFKL framework. Previous works were limited to LL or specific kinematic approximations.
2. Handling Color-Octet States: The authors explicitly address the additional soft divergences arising in Color-Octet (CO) states ($1S_0^{[8]}$and$3S_1^{[8]}$). Unlike Color-Singlet states, CO states allow for the radiation of soft gluons by the total color charge of the final state, requiring a specific subtraction term ($\Delta J^{(m, CO)}$) to isolate these singularities.
3. Analytic Derivation of Real Emission: They provide the full analytic expressions for the real-emission contributions, including the extraction of distributions (delta functions and $+$-distributions) necessary for numerical implementation.
4. Scheme Transformation: They derive the specific counterterms required to match the EFT calculation to the BFKL rapidity-factorization scheme, ensuring the universality of the Green's function.
5. Validation of Factorization: The work demonstrates the cancellation of soft divergences between real and virtual contributions and confirms that the remaining collinear singularities are compatible with standard factorization up to one loop.

4. Key Results

The paper presents the final expressions for the NLO impact factors $V^{(m,1)}_a$ for the gluon ($a=g$) and quark ($a=q$) channels.

• Structure of the Result: The NLO impact factor is decomposed into:
  $$V^{(m,1)} = V^{(m,1, \text{fin})} + V^{(m,1, \text{analyt})}$$

  • Finite Part ($V^{(m,1, \text{fin})}$): Obtained by integrating the subtracted real-emission matrix element (Eq. 3.1).
  • Analytic Part ($V^{(m,1, \text{analyt})}$): Contains the virtual corrections, the BFKL kernel contributions, and the counterterms. It is expressed in terms of distributions ($\delta(1-z)$, $(1-z)_+$) and logarithms.

• Specific Findings:

  • **For $1S_0^{[1]}$(CS):** The result is given in Eq. (3.12), combining the subtracted real part and the analytic part (Eq. 3.8). The virtual part includes finite terms$F^{(\text{virt})}$ derived in Ref. [38].
  • **For $1S_0^{[8]}$and$3S_1^{[8]}$(CO):** The result (Eq. 4.7) includes an **additional correction term**$\Delta V^{(m,1,8)}_g$(Eq. 4.5) to account for the extra soft divergence specific to octet states. This term involves complex functions including polylogarithms ($Li_2$) and inverse hyperbolic cotangents ($\coth^{-1}$).
  • Quark Channel: The quark-induced impact factors are also provided (Eq. 3.13), showing that the quark channel does not introduce the additional soft divergence seen in the gluon channel for CO states.

• Explicit Expressions: The paper provides the exact squared matrix elements for the real emission processes in Appendix B, which are too lengthy for the main text but are essential for numerical codes.

5. Significance and Future Applications

This work represents a major step forward in high-energy QCD phenomenology:

• Precision Phenomenology: It enables the first NLL-precision studies of forward-backward associated quarkonium production (e.g., quarkonium + jet or quarkonium + quarkonium) at the LHC. This is crucial for testing QCD dynamics in the semi-hard regime.
• Gluon Saturation: The results are vital for the search for gluon saturation (Color Glass Condensate), as quarkonium production is a sensitive probe. Current saturation calculations are mostly LL; this work provides the necessary NLL ingredients.
• Forward Physics: It allows for the computation of cross-sections for single forward quarkonium production with full NLL accuracy, resolving perturbative instabilities found in fixed-order collinear calculations at NLO.
• Future Work: The authors note that this framework paves the way for computing impact factors for P-wave states ($\chi_c, \chi_b$) and extending the analysis to other quarkonium states, completing the NLO NRQCD prediction for hadroproduction in the BFKL regime.

In summary, this paper provides the rigorous theoretical foundation required to move quarkonium production studies from Leading Logarithmic accuracy to Next-to-Leading Logarithmic accuracy within the BFKL formalism, bridging the gap between high-energy resummation and heavy-quark effective theories.