Exclusive photoproduction of a π⁰γ pair in the saturation framework

Basic Information

• Paper ID: 2511.11516
• Title: Exclusive photoproduction of a π⁰γ pair in the saturation framework
• Authors: M. Fucilla, S. Nabeebaccus, L. Szymanowski, S. Wallon, J. Yarwicka
• Institutions: Université Paris-Saclay (France), National Center for Nuclear Research (Poland), University of Manchester (UK)
• Classification: hep-ph (High Energy Physics - Phenomenology), hep-ex, nucl-ex, nucl-th
• Publication: Physics Proceedings Ultra-Peripheral Collisions 2, 010 (2025) 1-10
• Paper Link: https://arxiv.org/abs/2511.11516

Abstract

This paper investigates the exclusive photoproduction of π⁰γ pairs with large invariant mass, which represents a promising channel for studying gluon saturation effects. The process has been shown to be incompatible with collinear factorization approaches based on generalized parton distributions (GPDs) at leading twist. The authors employ a kT-dependent factorization approach in the small-x regime using the shock-wave formalism, which resums multiple gluon exchanges between the projectile and the dense nuclear target. The study reveals that the polarization amplitude undergoes a sign change as a function of the back-to-back transverse momentum |p⃗t| of the pion-photon pair, resulting in a characteristic dip-like structure in the differential cross section as a function of |p⃗t|.

Research Background and Motivation

1. Research Problem

This paper studies the exclusive photoproduction of π⁰γ pairs in ultra-peripheral collisions (UPCs): $\gamma(p_\gamma) + N(p_N) \rightarrow \gamma(k) + N'(p_{N'}) + \pi^0(p_\pi)$

The hard scale of this process is provided by the large invariant mass M²π⁰γ of the final-state photon-pion pair, enabling the application of perturbative QCD (pQCD).

2. Importance of the Problem

• Probing gluon saturation effects: At high-energy collisions and small-x, the gluon density becomes extremely high, making nonlinear effects (saturation) important
• Testing QCD factorization theorems: This process provides a unique perspective for understanding the applicability boundaries of QCD factorization
• UPC physics research: Provides theoretical predictions for future Electron-Ion Collider (EIC) and Large Hadron Collider (LHC) experiments

3. Limitations of Existing Approaches

• Collinear factorization breakdown: References 1,2 have demonstrated that this process violates collinear factorization due to the presence of Glauber pinches
• Double gluon exchange characteristics: This channel allows double gluon exchange with the nucleon, which is absent in other processes respecting charge and charge-parity conservation
• Endpoint singularities: Endpoint singularities appearing in collinear factorization calculations require regularization

4. Research Motivation

• Employ kT-dependent factorization in the small-x regime as an alternative approach
• Use the shock-wave formalism to incorporate gluon saturation effects
• Naturally regularize endpoint singularities through the transverse momentum of t-channel gluons

Core Contributions

1. First calculation: First computation of exclusive π⁰γ pair photoproduction in the saturation framework using the shock-wave formalism
2. Theoretical framework development: Establishes a hybrid factorization scheme combining high-energy factorization (small-x dynamics) and collinear factorization (pion distribution amplitude)
3. Discovery of new physics phenomena: Reveals that the polarization amplitude changes sign at specific transverse momenta, resulting in a unique dip structure in the differential cross section
4. Saturation effects comparison: Systematically compares model predictions with saturation effects (BK evolution) versus without saturation (BFKL evolution)
5. Detailed graphical analysis: Provides comprehensive analysis of the different sign contributions from photon emission before and after the shock wave

Methodology Details

Task Definition

Calculate the amplitude and differential cross section for the exclusive photoproduction reaction γN → π⁰γN', where:

• Input: Incident photon momentum pγ, target nucleon momentum pN
• Output: Final-state photon k, π⁰ meson pπ, recoil nucleon pN'
• Constraints: Semi-hard kinematics sγN ≫ M²π⁰γ ≫ Λ²QCD

Model Architecture

1. Light-Cone Coordinate System

Employs light-cone basis vectors: $n^\mu_1 \equiv \frac{1}{\sqrt{2}}(1,0,0,1), \quad n^\mu_2 \equiv \frac{1}{\sqrt{2}}(1,0,0,-1)$

Any four-momentum is decomposed as: $p^\mu = p^+ n^\mu_1 + p^- n^\mu_2 + p^\mu_\perp$

2. Shock-Wave Formalism

• t-channel gluon separation: Distinguishes internal and external gluons through rapidity cutoff
• Effective Lagrangian: Sums over all external gluon exchange diagrams, producing shock-wave backgrounds represented by Wilson lines
• Dipole operator: Color singlet exchange in diffractive processes represented as dipole operators, evolved according to the B-JIMWLK equation

3. Amplitude Calculation Procedure

Step 1: Projectile collision factor calculation Computes the γ → γqq̄ amplitude Mqq̄ through the shock-wave formalism, including 4 Feynman diagrams (Figure 1 shows 2 representative diagrams): $M^\mu_{qq\bar{q}} = \frac{\delta(p^+_q + p^+_{\bar{q}} + k^+ - p^+_\gamma)}{\sqrt{2p^+_\gamma}\sqrt{2p^+_q}\sqrt{2p^+_{\bar{q}}}\sqrt{2k^+}} \left(\frac{e^2_q\sqrt{N_c}}{(2\pi)}\right) \int d^2p_{1\perp}d^2p_{2\perp} \delta^2(...) \tilde{U}_{12}\Phi^\mu_{qq\bar{q}}$

where the dipole operator is: $\tilde{U}_{12} = \int d^2x_{1\perp}d^2x_{2\perp}e^{-i(x_{1\perp}\cdot p_{1\perp}+x_{2\perp}\cdot p_{2\perp})}\left(1-\frac{1}{N_c}\text{Tr}[U_1U^\dagger_2]\right)$

Step 2: Projection to π⁰ final state Projects the qq̄ pair to neutral pion through Fierz identity: $\Phi^\mu_{qq\bar{q}} \rightarrow \Phi^\mu_{\pi^0} = -\frac{if_{\pi^0}}{4\sqrt{N_c}}\int dz \phi_{\pi^0}(z) \text{Tr}[\hat{H}^\mu\hat{p}_\pi\gamma^5]$

Employs the asymptotic form of the pion distribution amplitude: $\phi_{\pi^0}(z) = 6z(1-z)$

Step 3: Dirac matrix trace calculation In the photoproduction limit (Q²=0), computes traces for transverse polarization configurations: $\text{Tr}[H^i_{1\perp}\hat{p}_\pi\gamma^5] = 4ip^+_\gamma H^i_{1\perp}(z,Q^2,p_{1\perp})$

Includes two contributions from photon emission before and after the shock wave.

Step 4: Complete amplitude construction Convolutes with the target collision factor: $\langle N'|\tilde{U}_{12}|N\rangle = -\frac{(2\pi)^3}{N_c}\delta(p^-_{N'N})\delta^2(p_{N'N\perp}+p_{1\perp}+p_{2\perp})\tilde{F}\left(\frac{p_{12\perp}}{2}\right)$

Final transverse reduced amplitude: $T_T = \left(\frac{-6sf_{\pi^0}e^2_q}{(2\pi)^2N_c}\right)\epsilon_{\gamma\perp\mu}\int d^2p_{1\perp}\tilde{F}(p_{1\perp})\int^1_0 dz\, z(1-z)H^\mu_\perp(z,Q^2,p_{1\perp})$

Technical Innovations

1. Hybrid Factorization Scheme

• High-energy factorization: Handles small-x dynamics of γN → qq̄γN'
• Collinear factorization: Describes qq̄ → π⁰ conversion through pion distribution amplitude
• Natural regularization: Transverse momentum kT of t-channel gluons automatically regularizes endpoint singularities

2. Polarization Structure Analysis

Finds only two non-zero polarization configurations Txy and Tyx (equations 28-29): $T_{xy} = 2C_T\int d^2p'_1\tilde{F}(\vec{p}'_1)\int^1_0 dz\left(\frac{z}{\vec{p}'^2_1}\right)[\text{complex expression}]$

3. Momentum Space Transformation

Performs momentum shift p⃗₁ → p⃗'₁ - zp⃗t to avoid spurious divergences in numerical integration

4. Coordinate Space Representation

Provides coordinate space expressions via inverse Fourier transform for model comparison (equation 32)

Experimental Setup

Data Sets and Models

This paper presents pure theoretical calculations employing the following non-perturbative models:

1. GBW Model (Primary Analysis)

Defines the critical dipole radius: $R_0(x_P) = \frac{1}{Q_s(x_P)} = \left(\frac{x_P}{x_0}\right)^{\lambda/2}\frac{1}{Q_0}$

where xP-Pomeron in the photoproduction limit: $x_P \approx \frac{M^2_{\pi^0\gamma}}{s} = \frac{\vec{p}^2_t}{\alpha_\pi\alpha_k s}$

Parameter values (fitted to proton DIS data):

• Q₀ = 1 GeV
• σ₀ = 23.02 mb
• x₀ = 3.04×10⁻⁴
• λ = 0.288

Momentum space target collision factor: $\tilde{F}\left(\frac{p_{12\perp}}{2}\right) = (2\pi)^2N_c\sigma_0\left(\delta^2\left(\frac{p_{12\perp}}{2}\right) - \frac{R^2_0}{\pi}\exp\left[-\frac{\vec{p}^2_{12}R^2_0}{4}\right]\right)$

2. BK Evolution Model (Saturation)

• Based on MV initial conditions
• Uses BK (Balitsky-Kovchegov) equation evolution
• Parameters determined by fitting HERA DIS data

3. BFKL Evolution Model (No Saturation)

• Obtained by switching off nonlinear terms in the BK equation
• Uses identical initial conditions (for fair comparison)
• Represents linear evolution dynamics

Evaluation Metrics

1. Polarization Amplitudes

• Txy and Tyx as functions of |p⃗t|
• Location of sign-change points
• Absolute values of amplitudes

2. Differential Cross Section

$\frac{d\sigma_{\gamma N\to\gamma\pi^0N'}}{...} = \frac{(2\pi)^4}{4s}\delta^4(p_\gamma+p_N-p_\pi-k-p_{N'})|T_T|^2 dLIPS$

as a function of transverse momentum |p⃗t|, at fixed pion rapidity η and center-of-mass energy √sγN.

Comparison Methods

• GBW vs BK vs BFKL: Compares predictions from different evolution equations
• Saturation vs non-saturation: Investigates the impact of nonlinear effects
• Photon emission timing: Analyzes contributions from emission before and after the shock wave

Implementation Details

Numerical Integration

• t-channel momentum |p⃗₁| and quark momentum fraction z integrals performed numerically
• Polar angle integral over p⃗₁ completed analytically
• Momentum shift employed to avoid numerical instabilities

Kinematic Variables

Three independent kinematic variables:

• Center-of-mass energy sγN
• Back-to-back transverse momentum |p⃗t|
• Pion rapidity η

Computational Framework

• Forward limit: ∆⃗ = 0
• Photoproduction limit: Q² = 0
• Semi-hard kinematics: sγN ≫ M²π⁰γ ≫ Λ²QCD

Experimental Results

Main Results

1. Amplitude Sign-Change Phenomenon (Figure 2)

Key Finding: Polarization amplitude Txy exhibits sign change as a function of |p⃗t|

• Amplitude is negative in low |p⃗t| region
• Amplitude crosses zero at a critical value
• Amplitude becomes positive in high |p⃗t| region

Numerical Characteristics:

• For η = 3, √sγN = 50 GeV, zero crossing occurs at approximately |p⃗t| ≈ 2-3 GeV
• Zero crossing location depends on center-of-mass energy and rapidity
• Amplitude variation near the zero crossing is very steep

2. Graphical Contribution Analysis (Figure 3)

Photon Emission Timing Effects:

• Emission after shock wave (containing D'₁ term): Contribution always negative
• Emission before shock wave (containing D'₂ term): Contribution always positive
• Sign-change mechanism: Competition between two contributions leads to total amplitude sign change

Physical Interpretation: $D'_1 = \frac{\vec{p}^2_t}{\alpha^2_\pi}, \quad D'_2 = \frac{\vec{p}^2_t}{\alpha_\pi\alpha_k} + \frac{(\vec{p}'_1-z\vec{p}_t)^2}{z(1-z)\alpha_\pi}$

Different propagator structures lead to different |p⃗t| dependencies.

3. Saturation Effects Comparison (Figure 4)

Differential Cross Section Behavior:

| Model | Dip Structure | Low |p⃗t| Behavior | High |p⃗t| Behavior | |-------|---------------|-----------------|------------------| | GBW (saturation) | Present | Rapid decrease | Power-law decay | | BK (saturation) | Present | Rapid decrease | Power-law decay | | BFKL (no saturation) | Absent | Monotonic decrease | Faster decay |

Quantitative Differences:

• Saturation models show cross section reduction of 1-2 orders of magnitude at |p⃗t| ≈ 2-4 GeV
• Dip location varies with energy and rapidity
• BK and GBW predictions are similar but differ in details

Important Caveat: The authors emphasize that BFKL predictions are obtained by switching off nonlinear terms in the BK equation using identical initial conditions. More definitive conclusions require refitting initial conditions with the BFKL equation.

Ablation Studies

Polarization Configuration Contributions

• Txy and Tyx are the only two non-zero configurations
• Txx and Tyy configurations vanish through Schouten identity and Lorentz invariance
• Interference between the two non-zero configurations produces the final cross section

Momentum Fraction Integration

Analysis of z-integration contributions shows:

• Endpoint regions (z→0 or z→1) remain finite due to transverse momentum regularization
• Central region (z≈0.5) provides maximum contribution
• Shape of pion distribution amplitude ϕπ⁰(z)=6z(1-z) affects overall normalization

Experimental Findings

1. Universality of Dip Structure

The sign-change phenomenon appears in all tested kinematic configurations, indicating this is an intrinsic feature of the process rather than a numerical artifact.

2. Saturation as a Signature

The presence or absence of dip structure can serve as an experimental signal for saturation effects:

• Saturation models predict pronounced dips
• Linear evolution does not produce this structure
• Provides a clear observational target for EIC experiments

3. Energy Dependence

With increasing center-of-mass energy:

• Overall amplitude increases (consistent with high-energy behavior)
• Zero crossing location shifts slightly
• Dip depth varies

Related Work

1. Collinear Factorization Breakdown Studies

Key References:

• 1,2 Nabeebaccus et al. (2024-2025): First proof that π⁰γ process violates collinear factorization
  • Discovers presence of Glauert pinches
  • Proves factorization breakdown for 2→3 processes at leading twist

Contrast: Other γρ, γπ processes 3-6 preserve collinear factorization 7,8 since charge and parity conservation forbid double gluon exchange

2. Shock-Wave Formalism Applications

Theoretical Foundations:

• 9 McLerran-Venugopalan (1994): Color Glass Condensate theory
• 10 Balitsky (1996): Operator expansion methods

Recent Applications 11-18:

• Diffractive dijet and trijet production
• NLO calculations for exclusive light vector meson production
• Twist corrections to deep virtual meson production
• Present work is first application of this framework to π⁰γ production

3. GPD Approach

Traditional Method:

• 7,8 Qiu-Yu (2022-2023): GPD treatment of single diffractive hard exclusive processes
• Applicable to processes that preserve collinear factorization

Advantages of Present Work:

• kT factorization remains valid when collinear factorization fails
• Naturally incorporates saturation effects
• No endpoint singularity problems

4. Saturation Models

GBW Model 20:

• Simplest phenomenological saturation model
• Parameters fitted to HERA data
• Used for preliminary analysis in this work

BK/BFKL Evolution 21-23:

• More precise small-x evolution equations
• BK includes nonlinear saturation terms
• BFKL represents linear approximation

Conclusions and Discussion

Main Conclusions

1. First complete calculation: First computation of exclusive π⁰γ pair photoproduction in the saturation framework using shock-wave formalism to handle small-x dynamics
2. New physics phenomena: Discovers that polarization amplitude changes sign with transverse momentum |p⃗t|, resulting in a unique dip structure in the differential cross section
3. Saturation effect signature: Dip structure presence is a clear signal of gluon saturation; linear BFKL evolution does not produce this feature
4. Hybrid factorization success: Demonstrates effectiveness of hybrid scheme combining high-energy and collinear factorization for this process
5. Experimental observability: Process shows measurement prospects in future EIC and LHC-UPC experiments

Limitations

1. Theoretical Level

• Leading-order approximation: Only leading-order amplitude computed; NLO corrections not included
• Forward limit: Assumes ∆⃗ = 0, neglecting finite momentum transfer effects
• Asymptotic DA: Uses simple asymptotic pion distribution amplitude without evolution effects

2. Phenomenological Models

• GBW simplification: GBW model is overly simple, not fully capturing saturation dynamics
• Initial condition dependence: BFKL/BK comparison affected by initial condition choice
• Parameter uncertainties: Target model parameters have fitting uncertainties

3. Experimental Challenges

• Cross section magnitude: Absolute cross section predictions require more precise theoretical inputs
• Background processes: Experiments must distinguish signal from background
• Kinematic coverage: Requires sufficient statistics to probe dip structure

Future Directions

Short-term Plans (Explicitly Stated in Paper)

1. Refit initial conditions: Refit DIS data using BFKL evolution equation for fairer comparison of linear versus nonlinear dynamics
2. Finite momentum transfer: Extend to non-forward kinematics, study effects of ∆⃗ ≠ 0
3. More realistic DA: Use evolved pion distribution amplitude, consider higher-order Gegenbauer moments

Medium to Long-term Prospects

1. NLO calculations: Include next-to-leading order corrections for improved theoretical precision
2. Other mesons: Generalize to charged pions, kaons, etc.
3. Nuclear effects: Study corrections in nuclear targets (already in heavy-ion collision context)
4. Experimental comparison: Detailed comparison with future EIC data
5. Machine learning: Explore neural network approaches for target model fitting

In-Depth Evaluation

Strengths

1. Methodological Innovation ⭐⭐⭐⭐⭐

• Breakthrough work: First treatment of this process in saturation framework, opening new research directions
• Hybrid factorization: Cleverly combines advantages of high-energy and collinear factorization
• Shock-wave technique: Successful application of mature framework to new process

2. Physical Insights ⭐⭐⭐⭐⭐

• Sign-change mechanism: Clearly reveals different sign contributions from photon emission timing
• Saturation signature: Dip structure as unique signal of saturation effects has important physical significance
• Factorization breakdown: Deepens understanding of QCD factorization applicability boundaries

3. Computational Completeness ⭐⭐⭐⭐

• Detailed derivations: Complete calculation chain from Feynman diagrams to final amplitude
• Analytical-numerical balance: Appropriate distribution of analytical and numerical integration
• Multi-model comparison: Systematic comparison of GBW, BK, BFKL predictions

4. Writing Clarity ⭐⭐⭐⭐

• Clear logic: Progression from theoretical framework to numerical results
• Complete formulas: Key formulas fully presented for reproducibility
• Effective figures: Figures 2-4 clearly display main findings

Weaknesses

1. Theoretical Completeness ⭐⭐⭐

• Leading-order limitation: NLO corrections may significantly alter numerical predictions
• Forward approximation: Finite-t effects may be important in actual experiments
• Uncertainty estimation: Lacks systematic theoretical uncertainty analysis (scale dependence, DA choice, etc.)

2. Phenomenological Models ⭐⭐⭐

• GBW oversimplification: As primary analysis model, GBW oversimplifies saturation dynamics
• BFKL comparison issue: Authors' own acknowledgment of initial condition inconsistency weakens linear vs. nonlinear comparison
• Parameter space: Insufficient exploration of parameter variation effects on results

3. Experimental Feasibility ⭐⭐⭐

• Absolute normalization: Lacks discussion of absolute cross section magnitude
• Kinematic range: Does not clearly specify optimal measurement window
• Background estimation: Does not discuss major background processes and signal-to-noise ratio

4. Numerical Details ⭐⭐⭐

• Convergence tests: No reported convergence tests for numerical integration
• Grid selection: Justification for integration grid and cutoff choices unclear
• Code availability: No mention of code release plans

Impact Assessment

Academic Impact ⭐⭐⭐⭐⭐

• Opens new direction: Establishes foundation for theoretical research on π⁰γ production
• Methodological contribution: Hybrid factorization scheme generalizable to other processes
• Citation potential: Expected to become important reference in the field

Experimental Value ⭐⭐⭐⭐

• EIC relevance: Directly serves future EIC experiment physics
• Observable predictions: Dip structure is clear experimental observation target
• Complementarity: Provides complementary information to traditional DIS measurements

Reproducibility ⭐⭐⭐

• Formula completeness: Calculations reproducible in principle
• Model parameters explicit: GBW parameters clearly provided
• Detail insufficiency: Some numerical implementation details need supplementation

Applicable Scenarios

1. Direct Applications

• EIC experiments: Electron-proton/nucleus collisions in UPC processes
• LHC-UPC: Proton-proton or proton-lead collisions in ultra-peripheral events
• Fixed-target experiments: High-energy photon beam experiments

2. Method Generalization

• Other meson pairs: π±γ, Kγ, ηγ, etc.
• Heavy quarkonia: J/ψγ, Υγ (requires modification)
• Dimeson production: ππ, KK, etc. (more complex)

3. Theoretical Research

• Saturation dynamics: Universal framework for small-x physics
• Factorization testing: Screening other potentially broken processes
• High-energy QCD: Connection between Regge theory and pQCD

Technical Details Supplement

Light-Cone Projection Techniques

This work cleverly exploits advantages of light-cone quantization:

• Natural hierarchy p⁺≫p⁻ in high-energy limit
• Wilson lines simplify in light-cone gauge A·n₂=0
• Polarization vectors' light-cone components clearly separated

Fierz Identity Application

Projection from qq̄ to π⁰ involves: $\bar{u}H^\mu v = \text{Tr}[H^\mu \hat{p}_\pi\gamma^5]\times\frac{-if_{\pi^0}}{4\sqrt{N_c}}\phi_{\pi^0}(z)$

This step crucially converts fermionic bilinear operators to meson matrix elements.

Numerical Integration Strategy

• z integration: 0,1 interval, endpoints regularized by kT
• |p⃗₁| integration: Segmented treatment, equation (32) divides into 0,z|p⃗t| and [z|p⃗t|,∞)
• Angular integration: Analytically completed using Bessel function properties

Key References

1. S. Nabeebaccus et al., PRD 111 (2025) 034040 - Proves π⁰γ process violates collinear factorization
2. L. McLerran & R. Venugopalan, PRD 50 (1994) 2225 - Color Glass Condensate theory foundation
3. I. Balitsky, Nucl. Phys. B 463 (1996) 99 - Operator expansion and JIMWLK equation
4. K. Golec-Biernat & M. Wüsthoff, PRD 59 (1998) 014017 - GBW saturation model
5. T. Lappi & H. Mäntysaari, PRD 88 (2013) 114020 - HERA data fitting and BK evolution

Overall Assessment

Recommendation: ⭐⭐⭐⭐⭐ (Strongly recommended for high-energy physics, nuclear physics, and small-x physics researchers)