NNLO QCD corrections to unpolarized and polarized electroweak structure functions in semi-inclusive deep-inelastic scattering

This paper presents next-to-next-to-leading order (NNLO) QCD corrections for both unpolarized and polarized semi-inclusive deep-inelastic scattering structure functions involving neutral and charged current interactions, demonstrating their significant phenomenological impact and reduced scale dependence for future Electron-Ion Collider analyses and global parton distribution function extractions.

The Gist

Imagine the proton, the tiny particle inside an atom's nucleus, not as a solid marble, but as a bustling, chaotic city. Inside this city live tiny citizens: quarks (the workers) and gluons (the glue holding them together). For decades, physicists have been trying to map this city, figuring out exactly where the workers live, how fast they move, and how they spin.

This paper is like a massive, ultra-high-definition upgrade to the map of this city, specifically for a new type of exploration called Semi-Inclusive Deep-Inelastic Scattering (SIDIS).

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

1. The Experiment: The "High-Speed Ping-Pong"

Imagine you are trying to figure out what's inside a sealed, opaque box (the proton). You can't open it, so you throw a very fast ball (an electron) at it.

• The Hit: The ball smashes into the box.
• The Result: The box breaks open, and you see not just the debris flying everywhere, but a specific, identifiable piece of debris (a new particle) flying out in a specific direction.
• The Goal: By watching how that specific piece flies out, you can deduce what the inside of the box looked like before the crash.

2. The Problem: The "Blurry Lens"

In the past, physicists had a "blurry lens." They could see the general shape of the city, but the details were fuzzy.

• Why? The math used to describe these crashes is incredibly complex. For a long time, they only calculated the "main event" (the first hit).
• The Glitch: Just like a photo that is slightly out of focus, their calculations had "fuzziness" (uncertainty) because they ignored the tiny, secondary interactions that happen during the crash. This made it hard to know exactly how the quarks and gluons were behaving.

3. The Solution: The "Next-Next-Next-Generation" Math

This paper presents a massive leap forward. The authors have calculated the effects of the collision up to the NNLO (Next-to-Next-to-Leading Order).

• The Analogy: Think of it like upgrading a video game from 8-bit graphics to 4K Ultra HD.
  • LO (Leading Order): You see the basic shapes of the buildings.
  • NLO (Next-to-Leading Order): You see the windows and doors.
  • NNLO (This Paper): You can see the bricks, the paint texture, and the people walking on the sidewalks.
• What they did: They calculated the "tiny ripples" and "secondary collisions" that happen when the electron hits the proton. They did this for two types of crashes:
  1. Neutral Current (NC): Like a polite handshake (using a photon or Z-boson).
  2. Charged Current (CC): Like a forceful punch (using a W-boson).

4. Why It Matters: The "Future Telescope"

The paper is written in preparation for the Electron-Ion Collider (EIC), a giant new machine being built to act as a super-microscope for the atomic world.

• The Challenge: The EIC will be so powerful that it will take pictures of the proton city in incredible detail. If the scientists' math (the map) isn't as sharp as the machine's lens, they won't be able to interpret the data correctly.
• The Fix: This paper provides the "sharp lens" for the math. By including these complex NNLO corrections, the "fuzziness" (uncertainty) in the predictions drops dramatically.
• The Result: Instead of guessing where the quarks are, scientists can now pinpoint them with high precision. They can also figure out how the "spin" of the proton (which way the city is rotating) is made up of the spinning workers and the glue.

5. The "Flavor" Twist

The paper also looks at "flavor." In the proton city, there are different types of quarks (Up, Down, Strange, etc.).

• The Analogy: Imagine trying to count how many red cars vs. blue cars are in a parking lot.
• The Breakthrough: This new math allows scientists to distinguish between these "flavors" much better. They can tell if a specific piece of debris came from an "Up" quark or a "Down" quark, which helps them understand how these different citizens form the proton.

Summary

In short, this paper is a mathematical masterpiece that cleans up the blurry edges of our understanding of the proton.

• Before: "We think the proton is made of quarks, but our calculations are a bit fuzzy, so we aren't 100% sure of the details."
• After (Thanks to this paper): "We have calculated every tiny interaction up to the highest possible precision. Now, when the new Electron-Ion Collider starts taking pictures, we will have a crystal-clear map to read the results, finally revealing the secrets of how matter is built and how it spins."

This work ensures that when the world's most powerful microscopes start working, the scientists reading the data will have the sharpest possible tools to understand the fundamental building blocks of our universe.

Technical Summary

1. Problem Statement

The primary objective of this research is to achieve a comprehensive understanding of the internal structure of hadrons (specifically protons and nuclei) by accurately determining Parton Distribution Functions (PDFs) and Fragmentation Functions (FFs).

• Context: Semi-Inclusive Deep-Inelastic Scattering (SIDIS) is a critical process for probing these functions simultaneously. It involves a lepton scattering off a hadron, producing a specific final-state hadron.
• The Gap: While Next-to-Leading Order (NLO) QCD corrections for SIDIS have been known for decades, and recent NNLO results exist for specific channels, there is a need for independent, comprehensive calculations that include:
  • Both Neutral Current (NC) ($\gamma, Z$ exchange) and Charged Current (CC) ($W^\pm$ exchange) interactions.
  • Both unpolarized and polarized structure functions (SFs).
  • Full flavor dependence and all partonic channels.
• Motivation: The upcoming Electron-Ion Collider (EIC) will operate with unprecedented precision. To extract PDFs and FFs from EIC data and to study electroweak effects (parity violation, flavor decomposition), theoretical predictions must match experimental precision, requiring Next-to-Next-to-Leading Order (NNLO) accuracy in perturbative QCD.

2. Methodology

The authors employed an independent computational framework to calculate the NNLO QCD corrections. The methodology involves several rigorous steps:

• Theoretical Framework:

  • The process is defined as $l(k_l) + H(P) \to l'(k_{l'}) + H'(P_H) + X$.
  • The differential cross-section is factorized into a leptonic tensor ($L_{\mu\nu}$) and a hadronic tensor ($W_{\mu\nu}$).
  • The hadronic tensor is decomposed into Lorentz structures defining the unpolarized SFs ($F_1, F_2, F_3$) and polarized SFs ($g_1, g_2, g_3, g_4, g_5$).
  • The structure functions are expressed as convolutions of PDFs, FFs, and perturbatively calculable Coefficient Functions (CFs).

• Computational Tools & Techniques:

  • Diagram Generation: Feynman diagrams were generated using QGRAF.
  • Algebraic Manipulation: The output was processed using FORM routines to handle Dirac algebra, Lorentz contractions, and color factors.
  • Handling $\gamma_5$ and Chirality: Since the calculation involves axial-vector couplings and polarized states, the $\gamma_5$ matrix and Levi-Civita tensors (intrinsically 4-dimensional) were handled using Larin's prescription in $d = 4 + \epsilon$ dimensions.
  • Renormalization: A finite renormalization scheme transformation was applied to convert results from Larin's scheme to the standard $\overline{\text{MS}}$ scheme, ensuring the restoration of chiral Ward identities.
  • Phase Space Integration:
    • The calculation categorizes contributions into Virtual-Virtual (VV), Real-Virtual (RV), and Double Real (RR) emissions.
    • Reverse Unitarity was used to map phase-space integrals to loop integrals.
    • Integration-by-Parts (IBP) identities (generated via LiteRed) were applied to reduce integrals to a set of Master Integrals (MIs).
    • The MIs were taken from existing literature, ensuring consistency.

• Mass Factorization:

  • Ultraviolet (UV) divergences were removed via coupling renormalization.
  • Infrared (IR) and collinear divergences were absorbed into the PDFs and FFs using space-like (for PDFs) and time-like (for FFs) Altarelli-Parisi kernels.
  • This process yielded the finite NNLO Coefficient Functions for all relevant channels ($\gamma\gamma, \gamma Z, ZZ, WW$).

3. Key Contributions

• Comprehensive NNLO Results: The paper presents the first independent calculation of NNLO QCD corrections for both unpolarized and polarized SFs in SIDIS, covering all electroweak channels (NC and CC).
• Full Flavor Dependence: The results incorporate contributions from all partonic channels ($q\bar{q}, qg, gg$, etc.) with full flavor dependence, including CKM matrix elements for CC processes.
• Analytical Expressions: All Coefficient Functions are expressed in terms of single-valued polylogarithms, allowing for efficient numerical evaluation across the entire kinematic range.
• Cross-Validation: The results provide a crucial cross-check against recent independent computations (e.g., Bonino et al.), finding complete analytical agreement for distribution terms and numerical agreement for regular parts.

4. Results and Phenomenology

The authors performed a detailed numerical analysis using NNPDF3.1 PDFs and NNFF1.0 FFs, simulating kinematics relevant to the EIC ($\sqrt{s} = 140$ GeV).

• Scale Dependence Reduction:
  • Leading Order (LO): Predictions showed large uncertainties due to strong dependence on the renormalization ($\mu_R$) and factorization ($\mu_F$) scales.
  • NNLO Impact: The inclusion of NNLO corrections significantly reduced the scale variation bands. The perturbative series showed apparent convergence, and the residual scale dependence was drastically minimized, particularly for the polarized structure functions.
• Electroweak Contributions:
  • NC Channel: The photon-mediated contribution dominates at lower $Q^2$. However, $Z$-boson exchange and $\gamma-Z$ interference become relevant at higher $Q^2$, introducing parity-violating asymmetries.
  • CC Channel: $W$-boson mediated contributions are suppressed by the massive propagator ($1/(Q^2 + M_W^2)$) and the chiral nature of the coupling, but they provide unique sensitivity to specific quark flavors (e.g., $u \to d$ transitions).
• Kinematic Sensitivity: The analysis demonstrated that the NNLO corrections are essential for precision predictions in the kinematic ranges ($x, z, Q^2$) expected at the EIC.

5. Significance

• EIC Readiness: These results provide the necessary theoretical infrastructure for the upcoming Electron-Ion Collider. They enable the extraction of PDFs and FFs with precision commensurate with the EIC's experimental capabilities.
• Flavor and Spin Structure: By isolating NC and CC contributions and including polarized SFs, these calculations allow for the disentangling of quark flavor contributions and the spin structure of the proton (how quark/gluon spin and orbital angular momentum contribute to the total nucleon spin).
• Global Fits: The results serve as a critical input for global QCD analyses, reducing theoretical uncertainties in the determination of parton distributions and fragmentation functions.
• Standard Model Tests: The precision of these NNLO predictions allows for stringent tests of the Standard Model, particularly in the electroweak sector (parity violation and flavor decomposition) within the strong interaction environment.

In conclusion, this work establishes a robust, high-precision theoretical framework for SIDIS, bridging the gap between current experimental capabilities and the future precision requirements of the EIC.