New Recursive Method Simplifies Gluon Production Calculations in Proton-Nucleus Collisions

Calculating Gluon Production Without Gauge Fixing

Physicists studying high-energy particle physics are analyzing a newly published method that calculates gluon production in proton-nucleus collisions without solving the Yang-Mills equations or relying on explicit Feynman diagrams. Published in the Journal of High Energy Physics and covered by Phys.org on November 4, 2026, the work introduces a direct approach to determining scattering amplitudes that bypasses traditional gauge-fixing complications.

Resolving Non-Unique Field Equations

In theoretical models of proton-nucleus collisions, projectiles operate as streams of color charges moving along their respective trajectories. These charges couple to gluons, ultimately generating the observed final-state gluons and quarks, which eventually produce hadrons. Conventionally, researchers calculate the probability of obtaining a specific final state by solving the Yang-Mills equations governing the gluon field relative to these color sources. However, as noted in reports from IPhT, those solutions are not unique. Infinitely many valid solutions describe the exact same physical situation through gauge transformations, forcing researchers to expend significant computational effort calculating terms that do not contribute to physical observables.

To eliminate this redundant computational overhead, François Gélis, a physicist at IPhT, treated solutions of the Yang-Mills equations as sums of tree Feynman graphs. By mapping this analogy onto methods originally developed for scattering amplitudes, Gélis applied an extension of the Britto-Cachazo-Feng-Witten recursion. This mathematical framework computes amplitudes directly without generating Feynman diagrams. The resulting calculations successfully reproduced outcomes previously derived via the traditional Yang-Mills approach, while entirely omitting the need to fix the gauge.

Expanding to Next-to-Leading-Order Corrections

In the color glass condensate framework, the elimination of gauge-fixing constraints opens new routes for higher-order calculations. Another theoretical work submitted to arXiv on April 8, 2026, and revised on June 25, 2026, by Paul Caucal and co-authors, examines forward trijet production in the gluon-initiated channel using covariant perturbation theory within the dilute-dense approximation. That research expresses results as convolutions of multiparton color correlators of light-like Wilson lines and perturbative impact factors, organized in compact expressions in terms of bare topologies, using effective vertices for the quark and gluon propagators.

By contrast, the BCFW recursion route demonstrated by Gélis avoids the explicit diagrammatic expansion entirely. This opens the possibility of applying generalized unitarity and other high-energy physics methods to proton-nucleus collisions to obtain next-to-leading-order corrections in a much simpler manner.

# Conceptual flow of the BCFW recursion approach for gluon amplitudes
def calculate_gluon_amplitude(color_charges, tree_graphs):
    # Bypass Yang-Mills non-unique gauge solutions
    recursive_amplitude = apply_bcfw_extension(tree_graphs)
    return recursive_amplitude

The arXiv work also observes that, in the ggg final state, the four-gluon vertex topology follows a structure similar to the instantaneous contributions, suggesting a simplification in the organization of multigluon color glass condensate amplitudes that may prove useful for future one-loop calculations.