Next-to-Next-to-Leading Order QCD Corrections to Polarized Semi-Inclusive Deep-Inelastic Scattering

Basic Information

• Paper ID: 2404.09959
• Title: Next-to-Next-to-Leading Order QCD Corrections to Polarized Semi-Inclusive Deep-Inelastic Scattering
• Authors: Saurav Goyal, Roman N. Lee, Sven-Olaf Moch, Vaibhav Pathak, Narayan Rana, V. Ravindran
• Classification: hep-ph (High Energy Physics - Phenomenology), hep-ex (High Energy Physics - Experiment)
• Submission Date: April 15, 2024 (arXiv submission), December 31, 2024 (latest revision)
• Paper Link: https://arxiv.org/abs/2404.09959

Abstract

Polarized semi-inclusive deep-inelastic scattering (SIDIS) is a key process for addressing the proton spin puzzle. This paper presents the first complete next-to-next-to-leading order (NNLO) results in perturbative quantum chromodynamics for the polarized SIDIS process. The analysis includes all parton channels for polarized lepton-hadron scattering with spin-averaged hadrons identified in the final state. Numerical analysis demonstrates the importance of NNLO corrections and shows reduced residual scale dependence within the kinematic range accessible to future electron-ion colliders (EIC).

Research Background and Motivation

The Proton Spin Puzzle

The proton spin puzzle is a long-standing unsolved problem in particle physics. Traditionally, the proton's spin was thought to originate primarily from the spin of its constituent quarks; however, experiments revealed that quark spin contributes only approximately 30% to the proton's total spin, far less than the expected 100%. This has motivated scientists to search for alternative sources of contribution, including gluon spin and orbital angular momentum.

Importance of SIDIS

1. Structure Function Measurements: Deep-inelastic scattering (DIS) provides crucial information about the internal structure of hadrons, with structure functions encoding parton distributions.
2. Polarization Information: Polarized SIDIS enables the separation and extraction of helicity distribution functions for different quark flavors.
3. QCD Factorization: The SIDIS process allows separation of short-distance perturbative dynamics from long-distance non-perturbative dynamics.

Existing Theoretical Limitations

• Previous theoretical predictions for polarized SIDIS reached only next-to-leading order (NLO) precision.
• With the construction of new-generation experimental facilities such as the EIC, higher-precision theoretical predictions are required.
• Advances in NNLO calculation techniques have made complete two-loop calculations feasible.

Core Contributions

1. First completion of NNLO QCD calculations for polarized SIDIS, filling an important theoretical gap in the field.
2. Provision of complete coefficient functions for all parton channels, including quark-quark, quark-gluon, gluon-quark, and gluon-gluon scattering.
3. Development of new techniques for handling complex phase-space integrals, particularly multi-loop integrals with additional constraints.
4. Numerical verification demonstrating that NNLO corrections significantly reduce scale dependence, enhancing the reliability of theoretical predictions.
5. Provision of high-precision theoretical predictions for EIC experiments, supporting future investigations of proton spin structure.

Methodology Details

Task Definition

Calculation of NNLO QCD corrections to the polarized semi-inclusive deep-inelastic scattering process: $l(k_l) + H(P) \rightarrow l(k'_l) + H'(P_H) + X$

where the inputs are a polarized lepton and polarized hadron, and the outputs are a scattered lepton, a tagged hadron, and other particles.

Theoretical Framework

QCD Factorization

The polarized structure function $g_1$ can be factorized as: $g_1 = \sum_{a,b} \int_x^1 \frac{dx_1}{x_1} \Delta f_a(x_1, \mu_F^2) \int_z^1 \frac{dz_1}{z_1} D_b(z_1, \mu_F^2) \times G_{1,ab}\left(\frac{x}{x_1}, \frac{z}{z_1}, \mu_F^2, Q^2\right)$

where:

• $\Delta f_a$: polarized parton distribution functions
• $D_b$: parton fragmentation functions
• $G_{1,ab}$: coefficient functions (the primary object of calculation in this work)

Perturbative Expansion

The coefficient functions are expanded in powers of the strong coupling constant: $G_{1,ab}(\mu_F^2) = \sum_{i=0}^{\infty} a_s^i(\mu_R^2) G_{1,ab}^{(i)}(\mu_F^2, \mu_R^2)$

Calculation Techniques

Regularization Scheme

• Dimensional regularization is employed to handle ultraviolet and infrared divergences.
• The Larin scheme is adopted for defining $\gamma_5$ in $D$-dimensional spacetime.
• Conversion to the $\overline{MS}$ scheme is achieved through finite renormalization constants.

Phase-Space Integration Techniques

1. Unitarity Method: Converts phase-space integrals into loop integrals.
2. IBP Identities: Reduces complex integrals to master integrals.
3. Differential Equations Method: Solves systems of differential equations for master integrals.
4. Boundary Conditions: Computed in the threshold limit.

Master Integral Calculations

• Real-virtual (RV) interference process: 7 master integrals
• Pure real (RR) emission process: 20 master integrals
• Results expressed in terms of polylogarithms and Nielsen polylogarithms

Technical Innovations

Constraint Handling

The SIDIS process includes an additional constraint $z/z_1 = p_a \cdot p_b/(p_a \cdot q)$, addressed through: $\delta\left(z' - \frac{p_a \cdot p_b}{p_a \cdot q}\right) \rightarrow -\frac{1}{\pi} \text{Im}\left(\frac{1}{z' - \frac{p_a \cdot p_b}{p_a \cdot q} + i\epsilon}\right)$

Analytic Continuation of Multivalued Functions

Hypergeometric functions encountered require analytic continuation across different kinematic regions: $_2F_1\left(1, \frac{\epsilon}{2}; 1+\frac{\epsilon}{2}; \frac{x'(1-x')}{z'(1-z')}\right)$

Smooth connection across regions is achieved through Pfaff transformations and the Feynman $i\epsilon$ prescription.

Experimental Setup

Numerical Calculation Parameters

• Parton Distribution Functions:
  • LO/NLO: MAPPDF10NLO
  • NNLO: MAPPDF10NNLO
• Fragmentation Functions:
  • LO/NLO: MAPFF10NLO
  • NNLO: MAPFF10NNLO
• Number of Active Quark Flavors: $n_f = 3$
• Scale Choice: $\mu_R^2 = \mu_F^2 = Q^2$

Kinematic Range

• COMPASS Energy: $\sqrt{s} = 17.4$ GeV
• EIC Energy: $\sqrt{s} = 45$ GeV
• z Integration Range: $0.2 < z < 0.85$
• y Constraint: $0.5 \leq y \leq 0.9$

Scale Variation Study

A seven-point scale variation scheme is employed: $(k_1, k_2) \in [1/2, 2]$ with $1/2 \leq k_1/k_2 \leq 2$

Experimental Results

Main Results

Contributions from Various Parton Channels

In $\pi^+$ production, the quark-quark scattering channel dominates, with contributions far exceeding those from quark-gluon and gluon-quark channels.

Numerical Impact of NNLO Corrections

1. Significant Reduction in Scale Dependence: The scale uncertainty band of NNLO predictions is noticeably narrower than NLO.
2. More Pronounced Corrections in Low $Q^2$ Region: At $Q^2 \sim 30$ GeV$^2$, clear differences exist between NNLO and NLO predictions.
3. Convergence in High $Q^2$ Region: As $Q^2$ increases, predictions at various orders gradually converge.

Comparison with Experimental Data

Figure 3 displays comparisons of the $g_1^{\pi^+}/F_1^{\pi^+}$ ratio with COMPASS experimental data:

• NNLO predictions show good agreement with experimental data.
• The reduction in theoretical uncertainty enables more precise comparison with experiments.
• Provides a reliable theoretical benchmark for future EIC experiments.

Scale Dependence Analysis

Scale variation studies at different $Q^2$ values reveal:

• $Q^2 = 30$ GeV$^2$: Strongest scale dependence with most significant NNLO corrections.
• $Q^2 = 800$ GeV$^2$: Weakest scale dependence with predictions from various orders approaching each other.
• Overall trends are consistent with expectations from perturbative QCD.

Related Work

Historical Development

1. NLO Calculations: de Florian, Stratmann, and Vogelsang (1998) first completed NLO calculations for polarized SIDIS.
2. Unpolarized SIDIS: NNLO calculations have been completed in recent years, providing important references for this work.
3. Threshold Logarithm Resummation: Advanced to N³LL precision.

Technical Development

• Loop Integration Techniques: Maturation of IBP methods and differential equation techniques.
• Phase-Space Integration: Application of the unitarity method.
• Multi-Loop Calculations: Development of automated calculation tools.

Relationship to This Work

This paper represents the first application of these advanced techniques to NNLO calculations for polarized SIDIS, achieving an important breakthrough in the field.

Conclusions and Discussion

Main Conclusions

1. Completed the first NNLO calculation for polarized SIDIS, providing complete results for all parton channels.
2. NNLO corrections significantly improve theoretical prediction accuracy, with notably reduced scale dependence.
3. Provides high-precision theoretical support for EIC experiments, facilitating precise measurements of proton spin structure.

Physical Significance

• Enhances the theoretical prediction precision for polarized structure functions.
• Provides a more reliable theoretical framework for extracting polarized parton distribution functions.
• Advances theoretical research on the proton spin puzzle.

Future Directions

1. Electroweak Corrections: Inclusion of electroweak interaction corrections.
2. Higher-Order Calculations: Development toward N³LO precision.
3. Phenomenological Analysis: Global fitting analysis combined with experimental data.

In-Depth Evaluation

Strengths

Technical Innovation

1. Computational Complexity: Successfully handled two-loop phase-space integrals with additional constraints.
2. Methodological Contribution: Developed systematic methods for analytic continuation of multivalued functions.
3. Numerical Stability: Achieved stable numerical calculations across the full kinematic range.

Theoretical Contribution

1. Completeness: Provides complete NNLO results for all parton channels.
2. Accuracy: Results verified through multiple independent methods.
3. Practicality: Results have been packaged into usable numerical programs.

Phenomenological Impact

1. Experimental Relevance: Provides important theoretical support for upcoming EIC experiments.
2. Precision Enhancement: Significantly reduces theoretical uncertainties.
3. Application Prospects: Establishes foundation for precise determination of polarized PDFs.

Limitations

Technical Constraints

1. Computational Complexity: The calculation process is extremely complex, making result verification difficult.
2. Numerical Evaluation: Numerical calculations in certain kinematic regions require further optimization.
3. Higher-Order Effects: Does not include electroweak corrections and other effects.

Methodological Restrictions

1. Regularization Dependence: Results depend on specific regularization scheme choices.
2. Scale Selection: Optimal scale choice requires further investigation.
3. Non-Perturbative Effects: Does not account for higher-twist effects.

Impact Assessment

Academic Impact

1. Pioneering Work: Fills an important theoretical gap in the field.
2. Methodological Value: Provides technical reference for similar complex calculations.
3. Citation Potential: Expected to become an important reference in the field.

Practical Value

1. Experimental Support: Directly serves major experimental projects such as the EIC.
2. Phenomenological Application: Provides important input for global PDF fitting.
3. Standard Establishment: May become the computational standard in the field.

Reproducibility

1. Code Availability: Authors commit to providing complete calculation code.
2. Result Verification: Verified through multiple independent methods.
3. Technical Details: Paper provides sufficient technical details.

Applicable Scenarios

1. High-Energy Physics Experiments: EIC, COMPASS, and other polarized scattering experiments.
2. Theoretical Research: Theoretical analysis of proton spin structure.
3. Phenomenological Research: Global fitting analysis of polarized PDFs.
4. Methodological Research: Technical development in complex multi-loop calculations.

Overall Assessment: This is a high-quality paper of significant theoretical value and practical importance. The authors have successfully completed the first NNLO calculation for the polarized SIDIS process, achieving not only a major technical breakthrough but also providing important theoretical tools for understanding proton spin structure. Despite the extreme complexity of the calculation, the authors have verified the correctness of their results through multiple methods and demonstrated the important physical effects of NNLO corrections. This work will have far-reaching impacts on both theoretical and experimental research in high-energy physics.