⚛️ phenomenology

Non-perturbative effects and soft-gluon dynamics in low- $p_T$ Drell-Yan production

This paper utilizes the novel \pythiaPB\ framework to investigate non-perturbative QCD effects and soft-gluon dynamics in low-transverse-momentum Drell-Yan production, revealing a specific center-of-mass energy dependence of the intrinsic transverse momentum distribution and highlighting the spectrum's sensitivity to the transition between perturbative and non-perturbative regimes at high energies.

Original authors: D. Subotić, H. Jung, A. V. Kotikov, N. Raičević

Published 2026-06-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: D. Subotić, H. Jung, A. V. Kotikov, N. Raičević

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

Imagine you are trying to understand how a massive, invisible cloud of particles behaves when two high-speed trains (protons) smash into each other. Specifically, physicists are interested in a specific "smoke" that comes out of the crash: pairs of particles called Drell-Yan lepton pairs.

The paper focuses on a very specific part of this smoke: the particles that are moving very slowly sideways (low transverse momentum, or low- $p_T$ ). In the world of particle physics, this slow-moving region is tricky because it's where the rules of "standard" math (perturbative QCD) start to break down, and the messy, fuzzy rules of "non-perturbative" physics take over.

Here is a breakdown of what the authors did, using simple analogies:

1. The Two Sources of the "Wobble"

When these particle pairs are created, they don't just sit still; they wiggle sideways. The paper explains that this wobble comes from two distinct sources, like two different people pushing a shopping cart:

The "Intrinsic" Push (The Jitter): Imagine the particles inside the proton are already jittering around before the crash even happens. This is called intrinsic transverse momentum. It's like the particles are shaking in their seats because of the proton's internal structure.
The "Soft Gluon" Push (The Wind): As the particles prepare to collide, they emit a spray of invisible, low-energy particles called soft gluons. These act like a gentle wind blowing the particles sideways. This wind builds up over time through many small emissions.

The paper argues that you cannot explain the sideways wobble by looking at just one of these. You need both the "jitter" and the "wind" to get the right picture.

2. The New Tool: A Better Simulator

The authors used a computer program called PYTHIA8-PDF2ISR. Think of this as a highly sophisticated flight simulator for particle collisions.

Usually, simulators have to make a choice: do we calculate the "wind" (soft gluons) using one set of rules, or do we just guess the "jitter" (intrinsic momentum) with a simple number?
This new approach tries to do both at once, ensuring the "wind" and the "jitter" work together consistently.

3. Fixing the "Recoil" Problem

One of the biggest headaches in these simulations is a concept called recoil.

The Analogy: Imagine you are standing on a skateboard (the proton) and you throw a heavy ball (a particle) forward. To conserve momentum, you and the skateboard must roll backward.
The Problem: In the standard version of the simulator (PYTHIA), when they calculated this backward roll, they accidentally spread the "kick" out among everyone on the skateboard, including the ball thrower. This made the "jitter" look smaller than it actually was, distorting the results.
The Fix: The authors introduced a new rule. They decided that when the ball is thrown, the "kick" should only be absorbed by the empty space behind the thrower (the "beam remnants"), leaving the thrower's original "jitter" untouched.
The Result: With this fix, the simulator's "jitter" setting matched the input setting perfectly. It was like calibrating a scale so that if you put a 1kg weight on it, it actually reads 1kg, not 0.8kg.

4. What They Found

By using this fixed simulator, they compared their predictions against real data from particle colliders (like the LHC at CERN and older experiments).

The "Jitter" Changes with Speed: They found that the amount of "jitter" needed to match the data changes depending on how fast the colliding protons are moving. At higher speeds, the "jitter" needed to be larger. This relationship was different from what other methods predicted.
The "Wind" Matters Most at High Speeds: At the highest energies (like 13 TeV), the "soft gluon wind" becomes the dominant factor. The data becomes so sensitive to this wind that it starts to reveal details about how the "glue" holding particles together (the strong force) behaves when it gets very weak or very strong.
Testing the "Glue" (Strong Coupling): The authors tested three different ways to mathematically describe the "glue" (the strong coupling constant, $\alpha_s$ ) in that messy, low-energy region. They found that at high energies, the shape of the particle spray is sensitive enough to tell the difference between these mathematical descriptions. It's like being able to hear the difference between a whisper and a murmur because the room is so quiet.

5. The Bottom Line

The paper concludes that to understand the messy, low-energy side of particle collisions, you cannot rely on just one trick. You need a simulator that correctly handles:

The natural jitter of the particles.
The accumulation of soft "wind" (gluons).
The correct way to handle the "kick" (recoil) when particles are thrown.

By fixing the "kick" issue in their simulator, they showed that their model works just as well as the most advanced methods currently used, but with a clearer understanding of how the "jitter" and the "wind" interact. This helps physicists use these collisions not just to find new particles, but to map out the invisible, fuzzy rules of the universe at its smallest scales.

Technical Summary: Non-perturbative effects and soft-gluon dynamics in low- $p_T$ Drell–Yan production

Problem Statement
The transverse-momentum ( $p_T$ ) spectrum of Drell–Yan (DY) lepton pairs at small $p_T(\ell\ell)$ serves as a critical probe for the transition between perturbative and non-perturbative Quantum Chromodynamics (QCD). In this kinematic region, the observed transverse momentum arises from two distinct mechanisms: the intrinsic transverse momentum ( $k_T$ ) of partons within the proton (a non-perturbative input) and soft-gluon radiation from initial-state parton evolution. A realistic description requires the simultaneous inclusion of both effects. While Parton Branching (PB) TMD approaches (e.g., in the CASCADE generator) treat intrinsic $k_T$ as a boundary condition and soft gluons dynamically, standard Monte Carlo parton showers (like PYTHIA) often struggle to disentangle these contributions due to recoil treatments and kinematic reconstruction differences. This paper investigates the behavior of the novel PYTHIA8-PDF2ISR framework, which implements a Parton Branching-based initial-state shower consistent with collinear PDFs, to systematically study the interplay of these non-perturbative effects.

Methodology
The authors utilize the PYTHIA8-PDF2ISR framework, which employs the same evolution ingredients as the PB method but operates within a Monte Carlo parton shower consistent with collinear PDFs. The study involves:

Recoil Prescription Modification: The standard PYTHIA 8 recoil-sharing step, which redistributes transverse momentum among all beam-related partons (initiators and remnants), was modified. A new "Remnant-only" recoil scheme was introduced where the intrinsic transverse momentum of hard initiators is preserved, and the full recoil required for momentum conservation is absorbed exclusively by beam remnants. This allows for a direct mapping between the input non-perturbative width and the extracted spectrum.
Intrinsic- $k_T$ Analysis: The width of the intrinsic- $k_T$ distribution ( $\sigma$ , related to $q_s = \sqrt{2}\sigma$ ) was determined by fitting PYTHIA8-PDF2ISR predictions to experimental DY data across a wide range of center-of-mass energies ( $\sqrt{s}$ from 38.8 GeV to 13 TeV).
Strong Coupling ( $\alpha_s$ ) Treatments: The sensitivity of the low- $p_T$ $p_{T}$ spectrum to the infrared behavior of the strong coupling $\alpha_s$ $α_{s}$ was investigated using three prescriptions for the transition from perturbative to non-perturbative scales:
- Option A (Hard Freezing): Freezing $\alpha_s$ below a cutoff $q_{min}$ .
- Option B (Fixed Gluon Mass): Introducing a fixed gluon mass $m_0$ .
- Option C (Running Gluon Mass): Introducing a scale-dependent running gluon mass.
Validation: Predictions were validated against experimental measurements from E605, PHENIX, CDF, D0, and CMS, comparing the PYTHIA8-PDF2ISR results with both PB-TMD calculations and standard PYTHIA recoil schemes.

Key Contributions

Improved Recoil Treatment: The paper introduces and validates a modified recoil prescription in PYTHIA8-PDF2ISR that preserves the intrinsic transverse momentum of hard initiators. This resolves discrepancies between the input and extracted Gaussian widths of the intrinsic- $k_T$ distribution, enabling a transparent comparison with PB-TMD inputs.
Quantitative Discrepancy in $q_s$ : The study demonstrates that the intrinsic- $k_T$ width ( $q_s$ ) required to describe DY data in PYTHIA8-PDF2ISR differs from the PB result. Specifically, the $q_s$ extracted from PYTHIA8-PDF2ISR exhibits a steeper dependence on $\sqrt{s}$ compared to PB, but is flatter than in standard PYTHIA. The authors attribute this to the frame in which transverse momentum is calculated and the specific kinematic reconstruction of the backward-evolution shower.
Sensitivity to $\alpha_s$ Infrared Behavior: The paper identifies that at the highest center-of-mass energies (LHC, $\sqrt{s}=13$ TeV), the low- $p_T(\ell\ell)$ DY spectrum becomes sensitive to the specific functional form of the transition between perturbative and non-perturbative regimes of $\alpha_s$ .

Results

Intrinsic- $k_T$ vs. Soft Gluons: The analysis confirms that intrinsic $k_T$ alone provides only a baseline smearing and cannot account for the observed width and shape of the spectrum at the lowest transverse momenta. The accumulation of transverse momentum through multiple soft-gluon emissions during initial-state evolution is essential.
Energy Dependence: The extracted $q_s$ values from PYTHIA8-PDF2ISR show a steeper $\sqrt{s}$ dependence than those from the PB method. For instance, at $\sqrt{s}=13$ TeV, $q_s \approx 1.7-1.9$ GeV, whereas at $\sqrt{s}=38.8$ GeV, $q_s \approx 0.9-1.0$ GeV.
Recoil Impact: The default PYTHIA recoil scheme leads to an effective reduction of the observed intrinsic- $k_T$ width, causing deviations from PB predictions. The modified "Remnant-only" recoil significantly improves the agreement between the input and extracted parameters.
$\alpha_s$ Sensitivity: At $\sqrt{s}=13$ TeV, variations in the infrared scale ( $m_0$ or $q_{min}$ ) produce visible differences in the lowest $p_T$ bins. Smaller infrared scales (e.g., $m_0 = 0.5$ GeV) enhance soft radiation, leading to an excess in the lowest $p_T$ bins. While different infrared prescriptions (Hard Freezing, Fixed Gluon Mass, Running Gluon Mass) yield similar spectra at lower energies (fixed-target and RHIC), they exhibit distinct shapes in the low- $p_T$ region at Tevatron and LHC energies.

Significance and Claims
The paper claims that the PYTHIA8-PDF2ISR approach, when equipped with the modified recoil treatment, provides a description of the low- $p_T$ DY spectrum comparable in quality to the PB-TMD approach, despite their different conceptual foundations (backward-evolution vs. forward-evolution).

The authors emphasize that the low- $p_T$ spectrum emerges from the interplay of intrinsic parton motion and soft-gluon radiation. They argue that the fitted intrinsic- $k_T$ distribution in PYTHIA8-PDF2ISR effectively absorbs not only genuine non-perturbative transverse motion but also kinematic effects associated with recoil sharing and frame transformations in the backward-evolution algorithm.

Crucially, the paper asserts that at high collision energies, the low- $p_T$ DY spectrum acts as a potential probe for the transition region of the strong coupling $\alpha_s$ . The sensitivity to the infrared behavior of $\alpha_s$ in the first few GeV of the dilepton transverse-momentum spectrum suggests that DY measurements can complement traditional extractions of $\alpha_s$ based on event shapes and jet observables, offering a unique window into non-perturbative QCD dynamics.

1. The Two Sources of the "Wobble"

2. The New Tool: A Better Simulator

3. Fixing the "Recoil" Problem

4. What They Found

5. The Bottom Line

More like this