Two-Dimensional Transverse-Momentum Subtraction and Semi-Inclusive Deep-Inelastic Scattering at N$^3$LO in QCD

This paper presents the first next-to-next-to-next-to-leading order (N$^3$LO) calculation of unpolarized semi-inclusive deep-inelastic scattering in perturbative QCD, utilizing a novel two-dimensional transverse-momentum subtraction method to establish a fully differential theoretical framework that enables precision nucleon tomography for upcoming Electron-Ion Collider experiments.

The Gist

Imagine you are a master chef trying to understand exactly how a complex dish is made. You know the ingredients (the basic particles), and you know the recipe (the laws of physics), but when you actually cook it, the kitchen gets messy. Smoke fills the air, ingredients splatter, and it's hard to tell exactly how much of each spice ended up in the final bite.

This is essentially what physicists face when studying Deep-Inelastic Scattering (SIDIS). They smash electrons into protons at incredibly high speeds to see what's inside the proton. The goal is to catch a specific "flavor" of particle (a hadron, like a pion) that flies out of the collision.

For decades, physicists have been able to calculate the recipe with decent accuracy (up to "NNLO" or Next-to-Next-to-Leading Order). But to match the super-precise cameras of future machines like the Electron-Ion Collider (EIC), they needed to go one step further: N3LO (Next-to-Next-to-Next-to-Leading Order).

Here is the problem: Calculating N3LO for this specific type of collision is like trying to count every single grain of sand on a beach while a hurricane is blowing. The math gets overwhelmed by "singularities"—mathematical infinities that happen when particles get too close to each other or move too slowly. Previous methods worked for simple collisions, but they broke down when a specific particle had to be identified in the final mess.

The New "Two-Dimensional" Knife

The authors of this paper, a team of theoretical physicists, have invented a new way to slice through this mathematical mess. They call it Two-Dimensional Transverse-Momentum Subtraction.

Think of the collision debris as a chaotic crowd of people running out of a stadium.

• Old Method: You tried to count everyone by looking at the whole crowd at once. It was too blurry.
• The New Method: The authors realized they could separate the crowd into two distinct groups based on two specific directions:
  1. How far off the main path they are running (Transverse Momentum).
  2. How much they are swaying side-to-side (Azimuthal decorrelation).

They created a "mathematical sieve" with two holes (two dimensions).

• Region A & B (The Quiet Zones): Most of the time, the particles behave predictably. The authors used known formulas (like a pre-written script) to calculate what happens in these calm areas.
• Region C (The Chaotic Zone): This is where the particles are running wild and colliding. Here, they used powerful computer simulations to calculate the exact chaos.

By subtracting the "quiet" predictions from the "chaotic" reality, the infinities cancel each other out, leaving a clean, precise number. It's like using a noise-canceling headphone to hear a single voice in a noisy room.

What Did They Find?

When they applied this new method to calculate the N3LO corrections:

1. The Recipe is Stable: The corrections they found were generally small (moderate). This is great news! It means the previous recipes (NNLO) were already pretty good, and the new math just fine-tuned them.
2. The "Edge" Cases Matter: However, in specific situations (like when the particle flies out at a very sharp angle, near the "threshold"), the corrections were significant. This is crucial for understanding the extreme edges of physics.
3. Less Guesswork: In physics, we often have to guess a "scale" (like the temperature of the oven) to make the math work. Changing this guess usually changes the result. With this new N3LO calculation, the result barely changes even if you tweak the scale. This means the prediction is robust and reliable.

Why Should You Care?

This isn't just about abstract math. This is the foundation for the Electron-Ion Collider (EIC), a massive new machine being built to take "3D movies" of the proton.

• Proton Tomography: Imagine taking an MRI of a human body. The EIC will do the same for a proton, showing how quarks and gluons are arranged inside.
• Precision: To get a clear MRI image, you need a perfect lens. This paper provides the "perfect lens" (the N3LO calculation). Without it, the images from the EIC would be blurry, and we might misinterpret the structure of matter.

The Bottom Line

The authors have successfully built a new, ultra-precise mathematical tool that allows us to predict the outcome of high-energy particle collisions with unprecedented accuracy. They solved a problem that was thought to be too messy to handle by slicing the problem into two manageable dimensions.

This achievement ensures that when the Electron-Ion Collider turns on, we will be ready to interpret its data perfectly, potentially unlocking secrets about how the universe holds itself together and how the proton's spin is generated. It's a giant leap from "guessing the recipe" to "perfectly baking the cake."

Technical Summary

1. Problem Statement

Precision predictions in perturbative Quantum Chromodynamics (QCD) are essential for interpreting collider data, particularly for the upcoming Electron-Ion Collider (EIC). While Next-to-Next-to-Next-to-Leading Order (N3LO) calculations have been achieved for inclusive processes (e.g., Higgs production, Drell-Yan), extending this precision to processes with identified hadrons in the final state has remained a major theoretical challenge.

• The Bottleneck: The primary difficulty lies in handling final-state collinear singularities associated with the observed hadron.
• Current Status: Calculations for identified hadron production have matured at Next-to-Next-to-Leading Order (NNLO), but the N3LO frontier was unexplored, with the sole exception of single-inclusive annihilation ($e^+e^- \to h + X$).
• Motivation: To match the unprecedented precision of future experiments (EIC, future $e^+e^-$ factories) and to enable rigorous extraction of non-perturbative fragmentation functions and nucleon tomography, theoretical predictions must be elevated to N3LO.

2. Methodology

The authors propose a novel Two-Dimensional Transverse-Momentum Subtraction scheme. This method generalizes the established $q_T$-subtraction technique (used for inclusive observables) to handle the specific kinematics of semi-inclusive deep-inelastic scattering (SIDIS).

• Kinematic Framework:

  • The calculation is performed in the Breit frame for the process $\gamma^* + p \to h + X$.
  • At Leading Order (LO), the hadron has zero transverse momentum ($P_{hT} = 0$). Beyond LO, it acquires transverse momentum.
  • The authors define a transverse momentum vector $\vec{P}_{hT}$ decomposed into components parallel ($p_{in}$) and perpendicular ($p_{out}$) to the scattering plane (defined by the beam and the recoiling leading jet).
  • These are parametrized by the beam polar angle ($\delta\theta$) and the azimuthal decorrelation ($\delta\phi$).

• Phase Space Partitioning:
  The cross-section is partitioned into three regions based on two small slicing parameters, $\Delta$ and $\lambda$ (where $\Delta, \lambda \ll 1$):

  1. Region A (Beam Collinear Limit): Defined by $\delta\theta < \Delta$ (i.e., $P_{hT} \ll Q$). The cross-section is governed by Transverse-Momentum-Dependent (TMD) factorization. It involves TMD beam functions, TMD fragmentation functions, and a soft function, all known up to three loops.
  2. Region B (Out-of-Plane Collinear Limit): Defined by $\delta\theta > \Delta$ but $\delta\phi < \lambda\Delta$ (i.e., small out-of-plane momentum $p_{out} \ll P_{hT}$). This region is also governed by a TMD factorization formula involving jet functions and soft functions.
  3. Region C (Resolved Region): Defined by $\delta\theta > \Delta$ and $\delta\phi > \lambda\Delta$. Here, the radiation is fully resolved. The calculation uses fixed-order NLO matrix elements for processes like $e_i \to e_j k l$ and $e_i \to e_j k l m$.

• Subtraction Mechanism:
  The method subtracts the singular contributions from Regions A and B (calculated via TMD factorization) from the full cross-section in the limit of vanishing slicing parameters. The leading-power dependence on $\Delta$ and $\lambda$ cancels out, leaving residual power corrections that vanish as the parameters approach zero.

• Inputs:

  • PDFs: CT18 NNLO (5 active flavors).
  • Fragmentation Functions (FFs): NPC23 NNLO.
  • Jet Functions: The calculation utilizes approximated values for two-loop constant terms of WTA (Winner-Take-All) unpolarized and linearly-polarized gluon jet functions, estimated via scaling relations with Energy-Energy Correlators (EEC).

3. Key Contributions

1. First N3LO Calculation for SIDIS: This is the first-ever N3LO calculation for unpolarized semi-inclusive deep-inelastic scattering ($e^- + p \to e^- + h + X$) involving identified hadrons in both initial and final states.
2. Novel Subtraction Scheme: The introduction of a 2D transverse-momentum subtraction method effectively handles the complex final-state collinear singularities that have hindered previous N3LO attempts.
3. Generalizability: The framework is designed to be applicable to:
   • Single-inclusive annihilation (SIA) in $e^+e^-$ collisions.
   • Identified hadron production within jets.
   • Polarized SIDIS (crucial for proton spin decomposition).

4. Results

The authors present phenomenological studies for charged pion ($\pi^+$) production at the EIC ($\sqrt{s} = 141$ GeV) and the COMPASS experiment ($\sqrt{s} = 17.3$ GeV).

• Numerical Stability:

  • The cross-section was verified to be independent of the arbitrary cutoff $\Delta$ in the limit $\Delta \to 0$.
  • Results for various partonic channels ($uu, ug, gu, gg, u\bar{u}, ud$) showed excellent flatness, confirming the precise cancellation of infrared divergences.

• Perturbative Convergence:

  • Magnitude: N3LO corrections are generally moderate (e.g., $-2\%$ to $+4\%$ for pion production as a function of momentum fraction $z$).
  • Convergence: The expansion shows excellent convergence. For dominant channels ($uu, ug, gu$), N3LO corrections are significantly smaller than NNLO corrections. For channels starting at NNLO ($ud, u\bar{u}, gg$), N3LO corrections can be comparable to or larger than NNLO, indicating the importance of higher orders.
  • Scale Uncertainty: The scale variation bands (hatched regions) are reduced by approximately 50% when moving from NNLO to N3LO, indicating that missing higher-order terms are now under control.

• Threshold Behavior:

  • In the large-$z$ (threshold) region, the N3LO corrections are dominated by threshold logarithms. The results show good agreement with approximated N3LO predictions derived from threshold resummation.

• Comparison with Data:

  • Comparisons with COMPASS data show that while scale uncertainties are reduced to well below experimental errors, the theoretical predictions consistently overshoot the data in the low-$z$ region. This suggests that uncertainties in the input fragmentation functions (FFs) remain a limiting factor.

5. Significance

• Theoretical Foundation for EIC: This work establishes the necessary theoretical precision to match the anticipated experimental accuracy of the Electron-Ion Collider. It enables precision nucleon tomography and the extraction of fragmentation functions with unprecedented rigor.
• Methodological Breakthrough: The 2D subtraction method provides a robust template for calculating N3LO corrections for any process involving identified hadrons, overcoming the long-standing barrier of final-state collinear singularities.
• Spin Physics: The framework is directly extendable to polarized SIDIS, which is critical for resolving the proton spin decomposition puzzle.
• Jet Substructure: The results provide essential perturbative ingredients for calculating energy correlator observables and the internal structure of jets at N3LO.

In summary, this paper represents a major milestone in perturbative QCD, successfully pushing the frontier of identified hadron production calculations to N3LO and providing the tools necessary for the next generation of high-energy physics experiments.