Heavy-quark contributions to the polarized DIS structure functions at NLO in the ACOT scheme

This paper presents analytical and numerical calculations of heavy-quark contributions to the polarized deep-inelastic scattering structure functions $g_1$, $g_4$, $g_5$, $g_6$, and $g_7$ at next-to-leading order within the ACOT scheme, offering new insights into spin-dependent QCD dynamics while excluding higher-twist suppressed functions $g_2$ and $g_3$.

The Gist

Imagine the proton (the core of an atom) not as a solid marble, but as a bustling, chaotic city. Inside this city, there are tiny citizens called quarks and messengers called gluons that zip around, holding everything together.

For decades, physicists have been trying to understand how these citizens spin. This "spin" is what gives the proton its magnetic personality. To study this, scientists smash high-speed electrons into protons (a process called Deep-Inelastic Scattering or DIS) and watch how the debris flies off.

However, there's a catch. Inside the proton, there are also "heavy" citizens: Charm and Bottom quarks. These are like the VIPs of the city—they are much heavier and harder to create than the regular citizens.

This paper is a new, highly detailed map of how these heavy VIPs behave when the proton is spinning, specifically looking at what happens when we smash it with electrons at the next level of precision (called NLO).

Here is the breakdown using simple analogies:

1. The Problem: The "Massless" Approximation vs. Reality

In the past, when physicists calculated how these heavy quarks behave, they often used a shortcut. They pretended the heavy quarks had no mass (like they were weightless ghosts).

• The Analogy: Imagine trying to calculate the traffic flow in a city where you pretend all the trucks are actually bicycles. It works fine when the road is empty and fast, but as soon as you get to a busy intersection or a steep hill (low energy), your bicycle math breaks down. The trucks (heavy quarks) actually have weight, and that weight changes how they move and interact.

2. The Solution: The ACOT Scheme

The authors used a specific mathematical framework called the ACOT scheme.

• The Analogy: Think of the ACOT scheme as a smart traffic controller.
  • When the energy is low (near the "threshold"), the controller treats the heavy quarks as heavy, slow-moving trucks that need extra space to turn.
  • When the energy is very high, the controller realizes the trucks are moving so fast they act like bicycles, and it switches to the simpler math.
  • The magic of this paper is that the authors figured out exactly how to switch between these two modes smoothly, without creating a traffic jam (mathematical errors) in the middle.

3. The Two Main Events

The paper calculates two main ways these heavy quarks get involved in the crash:

• Quark Scattering (QS): Imagine a heavy quark already inside the proton gets hit by the electron and kicks out a gluon (a messenger).
  • The Paper's Job: They calculated the exact "bounce" of this heavy quark, keeping its weight in mind the whole time.
• Boson-Gluon Fusion (GF): Imagine a gluon (messenger) inside the proton splits apart to create a new pair of heavy quarks (a heavy quark and an anti-quark).
  • The Paper's Job: They calculated the energy cost and the spin dynamics of creating these new heavy twins from scratch.

4. Why "Subtraction" is Necessary

When you combine the "heavy truck" math with the "bicycle" math, you might accidentally count the same traffic twice.

• The Analogy: If you count the trucks as trucks and then count them again as bicycles because they are moving fast, you have double-counted the traffic.
• The Fix: The authors developed a "correction formula" (subtraction terms) to remove the double-counting. It's like a referee blowing a whistle to say, "Stop, we already counted that!" ensuring the final number is perfectly accurate.

5. The Results: Why It Matters

The authors ran these calculations on a computer to see how much the "heavy truck" math changes the results compared to the old "bicycle" math.

• The Finding: Near the energy levels where heavy quarks are just starting to be created (the "threshold"), the difference is huge—up to 10%.
• The Metaphor: If you were trying to predict the winner of a race, and you ignored the fact that one runner is carrying a heavy backpack, you might be wrong by a significant margin. In the world of subatomic physics, a 10% error is massive.

6. The Future: The Electron-Ion Collider (EIC)

The paper is written in preparation for a new giant machine called the Electron-Ion Collider (EIC), which is being built to study these exact questions.

• The Takeaway: The EIC will be a high-precision camera. To interpret the photos it takes, the scientists need a perfect map. This paper provides that map. It tells the EIC scientists: "When you see a heavy quark spinning, here is exactly how to calculate its behavior, accounting for its weight, so you don't get confused."

Summary

This paper is a mathematical instruction manual for understanding how heavy, slow-moving particles behave inside a spinning proton. It fixes old shortcuts that ignored the particles' weight, ensuring that when the next generation of particle colliders (like the EIC) starts taking data, physicists will have the precise tools needed to decode the secrets of the proton's spin.

Technical Summary

Here is a detailed technical summary of the paper "Heavy-quark contributions to the polarized DIS structure functions at NLO in the ACOT scheme" by Spezzano et al.

1. Problem Statement

The spin structure of the nucleon is a fundamental open question in hadronic physics. While unpolarized heavy-quark contributions to Deep-Inelastic Scattering (DIS) are well-established at Next-to-Leading Order (NLO) and Next-to-Next-to-Leading Order (NNLO), the polarized sector has historically lacked a consistent, full-mass treatment at NLO.

• Context: Upcoming facilities like the Electron-Ion Collider (EIC) and the proposed EicC in China will provide high-precision data on spin-dependent observables.
• Challenge: In the kinematic regimes targeted by these colliders (near heavy-flavor thresholds and at high $x$), the zero-mass approximation (where heavy quarks are treated as massless partons) fails. It cannot account for the restricted phase space of massive final states or the transition from the massive threshold region to the asymptotic limit where heavy quarks act as active partons.
• Goal: To provide a rigorous, NLO calculation of heavy-quark contributions to polarized structure functions ($g_1, g_4, g_5, g_6, g_7$) that retains full dependence on heavy-quark masses ($m_c, m_b, m_t$) within a General-Mass Variable-Flavor-Number Scheme (GM-VFNS).

2. Methodology

The authors employ the ACOT (Aivazis-Collins-Olness-Tung) renormalization scheme, which ensures theoretical consistency across all kinematic regions by treating heavy quarks as massive near the threshold and as active partons at high scales.

A. Kinematics and Formalism

• Structure Functions: The study focuses on $g_1, g_4, g_5, g_6,$ and $g_7$. Functions $g_2$ and $g_3$ are excluded as they are higher-twist suppressed or determined by Wandzura-Wilczek relations.
• Processes: The calculation includes two primary partonic processes at NLO:
  1. Quark Scattering (QS): $B^* + Q_1 \to Q_2 + g$ (Real emission and Virtual corrections).
  2. Boson-Gluon Fusion (GF): $B^* + g \to Q_1 + Q_2$.
• Helicity Handling: Unlike massless schemes where $\gamma_5$ requires specific $D$-dimensional prescriptions (like HVBM), this calculation treats $\gamma_5$ strictly in $D=4$. The heavy-quark mass acts as a physical cutoff for collinear divergences, justifying the use of 4-dimensional projectors. Phase-space integration is performed in $D$ dimensions to regulate soft singularities.

B. Analytical Calculation

• Real Contributions: The authors derived explicit, closed-form analytic expressions for the polarized partonic structure functions ($\hat{g}_i$) for both QS and GF channels. These expressions depend on Mandelstam variables and the full heavy-quark masses ($m_1, m_2$).
• Virtual Corrections: One-loop vertex corrections were calculated. The self-energy insertions on external legs were canceled by mass and wave-function renormalization counterterms in the on-shell scheme.
• Subtraction Terms: To avoid double-counting (resummed logarithms in PDFs vs. massive coefficient functions), explicit subtraction terms were derived. These terms remove the zero-mass limit contributions from the massive NLO results. The authors adopted the CWZ decoupling prescription, ensuring subtractions are active only when the factorization scale $Q > m_Q$.

C. Numerical Implementation

• A dedicated Mathematica code was developed.
• The code interfaces with the ManeParse package to utilize the NNPDFpol1.1 polarized Parton Distribution Functions (PDFs).
• Electroweak parameters (masses, couplings, CKM matrix) were fixed to Particle Data Group values.
• The study covers Neutral Current (NC) and Charged Current (CC) channels.

3. Key Contributions

1. First Full-Mass NLO Polarized Calculation: This work provides the first complete NLO calculation of heavy-quark contributions to polarized structure functions ($g_1, g_4, g_5, g_6, g_7$) retaining full mass dependence for both initial and final state quarks.
2. Analytic Coefficients: The paper presents explicit analytic formulas for the massive Wilson coefficients in the ACOT scheme, including the necessary subtraction kernels.
3. Resolution of $\gamma_5$ Artifacts: By using 4-dimensional projectors and a specific subtraction term ($-2(1-\xi')$), the authors address the helicity non-conservation artifacts often associated with dimensional regularization in time-like quark scattering, ensuring consistency with standard $\overline{\text{MS}}$ results in the massless limit.
4. Violation of Dicus Relation: The authors explicitly demonstrate that the relation linking $g_4$ to $g_5$ and $g_6$ (analogous to the Callan-Gross relation in unpolarized DIS) is violated at NLO in the quark-scattering channel, even in the massless limit, due to the breakdown of the Dicus relation.

4. Results and Phenomenological Discussion

• K-Factors: The NLO corrections induce significant modifications to Leading Order (LO) predictions. Large K-factors (NLO/LO ratios) are observed, particularly as $x \to 1$, driven by Sudakov logarithms from soft-gluon emissions.
• Mass Effects near Thresholds:
  • Near the heavy-quark production thresholds (e.g., $Q^2 \approx m_b^2$), the massive ACOT results show a substantial suppression (up to 10% deviation) compared to the Zero-Mass (ZM) approximation.
  • This suppression reflects the physical constraint of the reduced phase space for massive final states, which the ZM approximation fails to capture.
• Asymptotic Behavior: At high scales ($Q^2 \gg m_b^2$), the ACOT results smoothly converge to the ZM limit, validating the scheme's consistency.
• Channel Specifics:
  • NC Channels: Results agree with previous literature for $g_1, g_4, g_5$.
  • CC Channels: The boson-gluon fusion process contributes significantly to charged-current structure functions, a feature absent in neutral-current interactions due to charge-conjugation symmetry.
• $g_6$ and $g_7$: These functions are suppressed by lepton mass squared ($m_\ell^2/Q^2$) and are negligible for electron beams (EIC) but could be relevant for $\tau$-lepton scattering.

5. Significance

• Precision Physics for EIC/EicC: The results provide a critical theoretical baseline for the upcoming Electron-Ion Collider. Accurate extraction of polarized gluon distributions and heavy-quark helicity densities requires a massive treatment near flavor thresholds, which this work supplies.
• Global QCD Analyses: The derived coefficients and subtraction terms can be directly integrated into global PDF fitting frameworks (like NNPDF) to reduce theoretical uncertainties in spin-dependent analyses.
• Future Directions: The framework is flexible and sets the stage for future NNLO calculations, threshold resummations, and the investigation of higher-twist effects relevant for $g_2$ and $g_3$.

In summary, this paper bridges a critical gap in perturbative QCD by delivering a theoretically consistent, numerically robust, and fully massive NLO description of heavy-flavor effects in polarized DIS, essential for the next generation of spin-physics experiments.