Analytical solution for QCD ⊗ QED evolution

Basic Information

• Paper ID: 2505.03520
• Title: Analytical solution for QCD ⊗ QED evolution
• Authors: Daniel de Florian, Lucas Palma Conte (ICAS, ICIFI & ECyT-UNSAM, Argentina)
• Classification: hep-ph (High Energy Physics - Phenomenology), hep-th (High Energy Physics - Theory)
• Publication Date: arXiv:2505.03520v2 hep-ph 15 Oct 2025
• Paper Link: https://arxiv.org/abs/2505.03520v2

Abstract

This paper presents analytical solutions for parton distribution function evolution incorporating mixed-order QCD ⊗ QED corrections, covering both polarized and unpolarized cases. By extending to mixed-order Altarelli-Parisi kernels, the DGLAP equations are solved exactly in Mellin N-space, and the corresponding Wilson coefficients for the polarized structure function g₁ are derived. The analytical approach not only improves computational efficiency but also enhances the theoretical prediction accuracy relevant to current and future phenomenological applications.

Research Background and Motivation

Problem Definition

1. Precision Requirements: QCD calculations have reached NNLO and even N³LO precision levels, at which point previously neglected QED corrections begin to play an important role in theoretical predictions
2. Complete Description Needs: To achieve comprehensive descriptions of hadronic processes, matching precision is required on the non-perturbative side (PDFs and polarized PDFs)
3. Mixed-Order Corrections: Inclusion of O(αₛα) QCD ⊗ QED mixed-order corrections is necessary to enhance the accuracy of theoretical predictions

Significance

• Phenomenological Importance: As QCD calculation precision improves, the contribution of QED corrections becomes non-negligible
• Experimental Requirements: Modern particle physics experiments demand increasingly higher precision in theoretical predictions
• Methodological Value: Analytical methods offer higher computational efficiency and better stability compared to numerical approaches

Limitations of Existing Methods

• Existing QED corrections remain primarily at LO order
• Systematic treatment of mixed-order O(αₛα) corrections lacks analytical solutions
• Computational efficiency requires improvement, particularly in repeated calculation scenarios

Core Contributions

1. First Complete Solution: Analytical solution of the DGLAP equations for mixed-order QCD ⊗ QED
2. Methodological Innovation: Development of two methods for solving mixed-order evolution equations (U-matrix method and Magnus expansion method)
3. Kernel Function Calculation: Computation of mixed-order AP splitting functions in the polarized case using the "Abelianization" algorithm
4. Wilson Coefficients: Derivation of mixed-order Wilson coefficients for the polarized structure function g₁
5. Efficiency Enhancement: Analytical methods significantly improve computational efficiency, particularly suitable for phenomenological studies

Detailed Methodology

Task Definition

Solving the DGLAP evolution equations with mixed-order QCD ⊗ QED corrections:

• Input: Parton distribution functions at initial scale Q₀
• Output: Parton distribution functions at arbitrary scale Q
• Constraints: Maintain O(αₛα) precision, covering both polarized and unpolarized cases

Core Theoretical Framework

1. Evolution Equation System

General form of the DGLAP equation:

df/dt = Σᵢⱼ Pᵢⱼ ⊗ fⱼ

where t = ln(Q²), Pᵢⱼ are AP splitting functions containing QCD, QED, and mixed-order contributions:

Pᵢⱼ = Σ_{a,b} aₛᵃ aᵇ P^(a,b)ᵢⱼ

2. Parton Basis Selection

To minimize mixing of different parton distributions during evolution, singlet/non-singlet basis is adopted:

• Non-singlet: fᴺˢ = {uᵥ, dᵥ, sᵥ, cᵥ, bᵥ, Δuc, Δds, Δsb, Δct}
• Singlet: fˢ = {ΔUD, Σ, g, γ}

3. Mixed-Order Splitting Functions

Mixed-order contributions are computed from NLO QCD results using the "Abelianization" algorithm. Main results include:

Photon-initiated processes:

ΔP^(1,1)_qγ = 2CᶠCₐe²_q{-22 + 27x - 9ln(x) + 8(1-x)ln(1-x) + ...}

Gluon-initiated processes:

ΔP^(1,1)_qg = 2Tᵣe²_q{-22 + 27x - 9ln(x) + 8(1-x)ln(1-x) + ...}

Solution Methods

Method 1: Extended U-Matrix Method

The evolution operator is decomposed as:

E(Q,Q₀) = E^QCD(Q,Q₀) E^QED(Q,Q₀) E^MIX(Q,Q₀)

where the mixed-order operator satisfies:

dE^MIX/dt = aₛa R^(1,1) E^MIX

The solution is:

E^MIX = exp(R^(1,1) ∫ aₛa dt)

Method 2: Magnus Expansion Method

Based on the Magnus expansion exponential form solution:

E^Mag = exp(Ω(t))

where:

Ω₁(t) = ∫ A(t₁)dt₁
Ω₂(t) = (1/2)∫∫ [A(t₁),A(t₂)]dt₁dt₂

Technical Innovations

1. Coupling Constant Evolution: Systematic treatment of RGE equations including mixed-order corrections
2. Non-commutativity Handling: Correction of non-commutative behavior of matrix kernels through Aᵢ operators
3. Mellin Space Techniques: Exploitation of convolution-to-product transformation advantages for improved computational efficiency
4. Analytical Integration: Development of analytical calculation methods for mixed-order integrals ∫aₛa dt

Experimental Setup

Datasets

• Polarized PDFs: DSSV polarized PDF set at initial scale Q₀² = 1 GeV²
• Polarized Photon PDFs: Calculations based on the LuxQED method
• Parameter Settings: Variable flavor number scheme with minimum 3 flavors; charm and bottom quark thresholds at m²_c = 2 GeV² and m²_b = 20.25 GeV² respectively

Evaluation Metrics

Relative corrections defined as:

δf = (f^QCD+QED - f^QCD+QED+MIX)/f^QCD+QED

Implementation Details

• Initial coupling constants: aₛ(M²_Z) ≃ 0.00946, a(M²_Z) ≃ 0.000610
• Implementation in FORTRAN combined with MT-1.0 Mathematica package for Mellin transforms
• Factorization scale setting: μ²_F = Q²

Experimental Results

Main Results

1. Mixed-Order Corrections to Coupling Constants

In the range Q² = 1000-10000 GeV², mixed-order corrections are of order O(10⁻⁴), with QED coupling corrections slightly larger than QCD coupling corrections.

2. Parton Distribution Function Corrections

At Q² = 1000 GeV²:

• Most polarized PDFs: Relative corrections of order O(10⁻⁴)
• Polarized Photon PDF: Corrections reach percentage level, reaching ~12% in high-x regions
• Non-singlet Combinations: Complete consistency between the two methods

3. Structure Function g₁ Corrections

At Q² = 1000 GeV²:

• At x ≃ 0.1: Relative corrections of order O(10⁻⁴)
• As x → 1: Corrections increase to O(10⁻³)
• Primary contributions arise from mixed-order Wilson coefficients, particularly quark channel contributions

Method Comparison

The U-matrix method and Magnus expansion method show basic consistency in mixed-order correction predictions, with differences primarily arising from different treatments of higher-order QCD terms.

Important Findings

1. Photon PDF Sensitivity: Mixed-order corrections have the most significant impact on photon PDFs
2. Wilson Coefficient Dominance: Mixed-order corrections to structure function g₁ primarily originate from Wilson coefficients rather than PDF evolution
3. Computational Efficiency: Analytical methods significantly enhance computational speed compared to numerical approaches

Related Work

Historical Development

1. QCD Side: NNLO splitting function calculations Moch et al. and modern PDF analyses NNPDF, CT, MSHT, etc.
2. QED Corrections: LO QED kernels Roth, Weinzierl and "Abelianization" algorithm de Florian et al.
3. Polarized Case: Polarized QCD splitting functions LO to NNLO and global analyses DSSV, NNPDF, JAM, etc.

Contributions of This Work

• First systematic treatment of mixed-order QCD ⊗ QED evolution
• Development of two complementary analytical solution methods
• Provision of complete mixed-order corrections in the polarized case

Conclusions and Discussion

Main Conclusions

1. Successfully solved the DGLAP equations for mixed-order QCD ⊗ QED
2. Mixed-order corrections have significant impact on photon PDFs, reaching percentage levels
3. Analytical methods substantially enhance computational efficiency while maintaining precision
4. Mixed-order corrections to polarized structure function g₁ primarily originate from Wilson coefficients

Limitations

1. Perturbative Expansion: Only considered to O(αₛα) order; higher-order corrections not included
2. Lepton Distributions: Coupling of lepton distributions not included (factorizable at current precision)
3. Numerical Approximations: Some integrals still require numerical calculation, particularly higher-order terms in the Magnus method

Future Directions

1. Extension to higher-order mixed corrections
2. Complete treatment including lepton distributions
3. Application to specific phenomenological process calculations
4. Optimization of numerical implementation for further efficiency enhancement

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Systematic treatment of mixed-order corrections with rigorous mathematical derivations
2. Methodological Innovation: Development of two complementary analytical solution methods
3. Practical Value: Significant computational efficiency improvement suitable for phenomenological applications
4. Completeness: Coverage of both polarized and unpolarized cases with complete theoretical framework
5. Sufficient Verification: Consistency between two methods enhances result reliability

Weaknesses

1. Complexity: Magnus method becomes complex and difficult to handle at higher orders
2. Approximate Treatment: Some non-commutative terms handled with approximations
3. Phenomenological Verification: Lacks direct comparison with experimental data
4. Applicable Range: Primarily suitable for intermediate Q² ranges; extreme regions may require additional considerations

Impact

1. Theoretical Contribution: Provides important tools for high-precision QCD calculations
2. Methodological Value: Analytical methods can be generalized to similar problems
3. Phenomenological Applications: Important implications for single-photon production in electron-proton collisions and similar processes
4. Future Development: Establishes foundation for higher-order correction calculations

Applicable Scenarios

1. Precision Measurements: Experimental analyses requiring high-precision theoretical predictions
2. Phenomenological Studies: Research requiring repeated PDF evolution calculations
3. Global Fitting: PDF global fitting including QED corrections
4. New Physics Searches: Research requiring precise Standard Model predictions as background

References

The paper includes 61 references, providing a comprehensive literature system spanning from classical Altarelli-Parisi work to the latest NNLO calculations and modern PDF analyses, establishing a solid bibliographic foundation for related research.