⚛️ phenomenology

A simple algorithm for polarized parton evolution

This paper presents a linearly scaling algorithm for polarized parton evolution that incorporates correlations between gluon production and decay planes by identifying charge currents, demonstrating agreement with fixed-order calculations and introducing a new observable for probing beyond-current-current interactions.

Original authors: Stefan Höche, Mareen Hoppe, Daniel Reichelt

Published 2026-03-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Stefan Höche, Mareen Hoppe, Daniel Reichelt

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: The Cosmic Dance of Particles

Imagine you are watching a high-energy particle collision at a place like the Large Hadron Collider (LHC). It's like a cosmic pinball machine where tiny particles smash together and explode into showers of new particles. Physicists use computer programs called Parton Showers to simulate these explosions.

For decades, these simulations have been very good at predicting how many particles come out and where they go. However, they have mostly ignored a subtle but crucial detail: spin and orientation.

Think of a particle like a spinning top or a radio antenna. Just as a radio antenna needs to be aligned correctly with the signal to get a strong reception, particles have a specific "orientation" (polarization) when they are created and when they decay. If your simulation ignores this alignment, it's like trying to tune a radio while holding the antenna upside down—you get the signal, but it's weak and distorted.

The Problem: The Old Way Was Clunky

The current method for handling these orientations (called the Shatz-Collins-Knowles method) is like trying to navigate a city using a massive, detailed paper map that you have to constantly unfold and refold.

It's complicated: It requires complex math for every single step.
It's slow: It gets very heavy on the computer's memory as the number of particles grows.
It's fragile: It breaks down or gets messy when particles are emitted at wide angles (like a car swerving suddenly).

The New Solution: A Simple "Antenna" Algorithm

The authors of this paper propose a new, much simpler way to handle these orientations. They call it a simple algorithm.

Here is the core idea, broken down with analogies:

1. The "Dipole" Antenna Concept

In physics, particles often come in pairs or groups that act like a dipole antenna (think of a classic TV antenna with two prongs).

The Old Way: Tried to calculate the spin of every single particle individually, like trying to track the wind direction for every single leaf in a storm.
The New Way: The authors realized you don't need to track every leaf. You just need to know the direction of the antenna itself. If you know the orientation of the "emitting antenna" (the particle creating the new one) and the "receiving antenna" (the particle absorbing the energy), you can predict how they interact.

2. The "Pass the Baton" Rule

The new algorithm works like a game of passing a baton in a relay race, but the baton is a piece of "orientation information."

Step 1 (Creation): When a particle splits (like a parent having a child), it creates a new particle. The algorithm simply says: "Hey, the new particle inherits the orientation of the parent's antenna." It stores this direction.
Step 2 (Decay): When that new particle eventually splits again, it looks back at the stored orientation. It asks: "Did I inherit the right alignment?" If the "emitting antenna" and the "absorbing antenna" are aligned (co-polarized), the interaction is stronger. If they are misaligned, it's weaker.
Step 3 (Simplicity): The math for this is incredibly simple. It doesn't require heavy lifting. It scales linearly, meaning if you double the number of particles, you only double the work. It's like walking a straight line instead of climbing a mountain.

Why Does This Matter?

1. Speed and Efficiency

The old method was like driving a tank through a forest; it was powerful but slow and clumsy. The new method is like a bicycle. It's fast, efficient, and can handle any terrain (whether particles are moving straight or swerving wide). This allows scientists to run more complex simulations in less time.

2. Seeing the "Hidden" Patterns

The paper introduces a new observable (a way to measure things). Imagine you are watching a dance floor.

Without polarization: You just see people moving around randomly.
With polarization: You see that the dancers are actually moving in synchronized patterns based on the music's rhythm.
The new algorithm allows physicists to see these "dance patterns" (correlations between production and decay planes) that were previously blurred out. This is crucial for the Future Circular Collider (FCC), a planned super-powerful collider that will produce so many particles that we need to understand these tiny details to find new physics.

3. Fixing the "Glitch"

The authors proved that for most situations, this simple "antenna" approach captures 99% of the physics. The tiny bits they ignore are so rare and complex (like quantum interference effects) that they only matter if you are looking for something incredibly specific. For general purposes, this new algorithm is the "Goldilocks" solution: not too simple, not too complex, just right.

The Takeaway

This paper is about simplifying the complex. The authors found a way to describe the orientation of subatomic particles using the simple concept of antenna alignment.

Before: We tried to calculate the spin of every particle individually, which was slow and prone to errors.
Now: We track the "antenna" direction of the particle groups. It's fast, accurate, and easy to implement in computer simulations.

This breakthrough means that in the future, when we smash particles together at record-breaking energies, our computer models will be sharp enough to see the subtle "dance" of the particles, helping us understand the fundamental laws of the universe with unprecedented clarity.

1. Problem Statement

Parton showers are essential tools in collider physics for modeling the evolution of particles from high to low energy scales, bridging the gap between hard scattering processes and the hadronic jets observed in detectors. While standard parton showers accurately model the rates of radiation, they often neglect polarization effects (spin correlations) between the production and decay planes of gluons.

Current Limitations: The standard method for including polarization, the Shatz-Collins-Knowles (SCK) algorithm, relies on spin-density matrices derived from spin-dependent decay matrix elements. A major drawback of the SCK method is that it requires specific, complex modifications to handle soft wide-angle gluon emissions, making it difficult to implement universally across different recoil schemes.
The Need: As experimental precision increases (e.g., at the High-Luminosity LHC and future FCC), the ability to resolve detailed radiation patterns and substructure requires parton showers that correctly account for polarization correlations without sacrificing computational efficiency or generality.

2. Methodology

The authors propose a new, simplified algorithm based on the Lorentz structure of QCD splitting functions, specifically focusing on the decomposition of these functions into charge currents.

A. Theoretical Foundation

The authors re-derive leading-order splitting functions ( $q \to qg$ , $g \to q\bar{q}$ , $g \to gg$ ) using a Sudakov decomposition of momenta. They identify that the polarization dependence of these splittings can be expressed through two fundamental currents:

Scalar Radiation Current ( $S^\mu$ ): Describes the emission of a gluon from a color dipole. It is defined as $S^\mu(p_1, p_2) = (2p_1 + p_2)^\mu / p_{12}^2$ .
Decay Current ( $D^\mu$ ): Describes the absorption of a gluon by a scalar dipole radiator. It is defined as $D^\mu(p_1, p_2) = d^\mu_\nu(p_{12}, \bar{n}) (p_1 - p_2)^\nu / p_{12}^2$ .

Key theoretical insights include:

Factorization: The splitting tensors for $q \to qg$ and $g \to gg$ can be factorized into products of these currents.
Interference Suppression: The authors prove that the "interference" terms in the $g \to gg$ splitting function (which do not factorize into simple current-current products) vanish identically in the limit of zero gluon virtuality or contribute only at $\mathcal{O}(\alpha_s^3)$ . These terms represent true quantum interferences that are sub-leading compared to the dominant classical correlations.
Color Flow: By utilizing the color-flow representation, the algorithm treats gluons as bi-fundamental charges, allowing the scalar current to be written as a sum over charge dipoles.

B. The Algorithm

The proposed Monte Carlo algorithm operates on the principle of storing and transferring polarization vectors associated with color dipoles:

Emission (Production): When a gluon is emitted via the scalar splitting function, the algorithm stores the polarization vector of the emitting color dipole (derived from the normalized scalar current $J^\mu_{i\bar{i}}$ ). If the emission is purely fermionic (non-polarized), no vector is stored.
Decay (Absorption): When a gluon splits (decays), the algorithm constructs a decay polarization vector ( $j^\mu_d$ $j_{d}^{μ}$ ).
- If the incoming gluon has a stored polarization vector, it is contracted with the decay vector. The splitting probability is reweighted by the square of this contraction ( $|j_p \cdot j_d|^2$ ).
- If the incoming gluon has no stored vector (e.g., from a fermionic emission), the decay vector is assigned to its correlated partner.
Propagation: Polarization information and correlation partners are transferred from the incoming particle to the emitter during subsequent splittings.

Complexity: The algorithm scales linearly ( $\mathcal{O}(N)$ ) with the number of particles in the final state regarding both computing time and memory, making it highly efficient compared to methods requiring full matrix element calculations.

3. Key Contributions

Simplification: The algorithm eliminates the need for the complex modifications required by the SCK method for soft wide-angle emissions. It is applicable in both the hard-collinear and soft wide-angle regions.
Current-Based Formalism: It establishes a direct link between the polarization of the vector field and the orientation of the emitting/absorbing QCD color dipoles, providing a physically intuitive interpretation of spin correlations.
Proof of Sufficiency: The authors demonstrate that for matrix elements with up to three additional partons, the current-current correlation is sufficient to capture all polarization effects. Effects from four or more partons are shown to be strongly sub-leading ( $\mathcal{O}(\alpha_s^3)$ ).
New Observable: They propose a current-current correlator ( $C_{i,j}^{1,2}$ ) as a theoretical tool to probe polarization effects beyond standard angular radiation patterns. This observable tests the co-polarization of emitting and absorbing antennae.

4. Results and Validation

The algorithm was implemented in the Alaric parton shower and validated against fixed-order perturbative calculations:

Kinematic Limits: The implementation reproduces the expected behavior in the collinear and soft limits.
- In the double-soft limit, the correlator reduces to the azimuthal angle defined in previous literature.
- In the collinear limit, it matches the squeezed correlator limit.
Validation Cases:
- $e^+e^- \to q\bar{q}q'\bar{q}'$ : The algorithm correctly reproduces the linear dependence on the correlator $C_{i,j}^{1,2}$ for gluon-to-quark splitting, showing a 50% reduction in correlation effects at specific kinematic points, consistent with theory.
- $e^+e^- \to q\bar{q}gg$ : The radiation pattern matches the product of production and decay currents $(J^\mu D_\mu)^2$ .
- Multi-gluon emissions: Tests with 3 and 4 emissions confirmed that the algorithm correctly handles the transfer of polarization information and the suppression of interference terms.
Comparison with SCK: Unlike the SCK algorithm, the new method maintains correct azimuthal modulation even at small rapidities without requiring ad-hoc adjustments.
Full Shower Effects: In full parton-shower evolution, the authors note that polarization effects are suppressed (by a factor of $\sim 50\%$ ) due to the dominance of scalar emissions, a known phenomenon that makes experimental measurement challenging but confirms the algorithm's consistency with standard shower behavior.

5. Significance

Computational Efficiency: The linear scaling makes this the most efficient technique currently available for including polarization in parton showers, enabling its use in large-scale Monte Carlo simulations for the High-Luminosity LHC and Future Circular Collider (FCC).
Precision Physics: By correctly modeling spin correlations, the algorithm improves the accuracy of jet substructure modeling. This is crucial for:
- Calibrating particle detectors.
- Inferring properties of hard interactions.
- Improving hadronization models, which is vital for precision measurements of Higgs couplings.
New Physics Probes: The proposed current-current correlator offers a new observable to search for polarization effects that go beyond standard splitting functions, potentially revealing new quantum interference effects or testing the limits of the factorization hypothesis in QCD.
Implementation: The method is straightforward to implement in existing shower frameworks (dipole, antenna, or traditional showers) that respect color coherence, facilitating immediate adoption by the community.

In summary, this paper presents a robust, efficient, and theoretically sound framework for polarized parton evolution, resolving long-standing implementation difficulties and paving the way for higher-precision QCD simulations in the era of high-luminosity colliders.