Multiplicity distributions in QCD jets and jet topics

✨

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 at a massive, chaotic fireworks display. You throw a single firework into the sky, and it explodes. That explosion isn't just one big bang; it's a cascade. The first explosion shoots out sparks, which hit other sparks, which hit more sparks, creating a shower of light and debris.

In the world of particle physics, this is exactly what happens when we smash protons together at the Large Hadron Collider (LHC). We create a Jet: a high-speed spray of particles (hadrons) shooting out from a single, high-energy collision.

This paper is about trying to predict how many particles will be in that spray.

The Problem: Counting the Sparks

Physicists have a theory called Quantum Chromodynamics (QCD) that explains how these particles interact. They wanted to know: If I shoot a jet with a certain energy, can I predict the probability of it having 10 particles, 20 particles, or 100 particles?

Decades ago, scientists proposed a neat rule called KNO Scaling. Think of it like a "universal recipe."

The Idea: If you have a small firework and a giant firework, the pattern of how the sparks spread out is the same. You just have to scale the numbers up or down.
The Analogy: Imagine a small crowd of people and a huge stadium crowd. If you look at the density of people relative to the total size, the "crowd shape" looks similar. The theory said all jets, no matter how energetic, should follow this same universal shape when you adjust for their size.

The Old Theory vs. Reality

For a long time, physicists used a simplified version of the math (called DLA) to predict this shape. It was like using a basic sketch to draw a masterpiece.

The Flaw: The sketch was okay, but when they compared it to real data from the ATLAS experiment at CERN, the picture was blurry. The real jets had more "spikes" and "valleys" in their particle counts than the simple math predicted. The old theory was missing a crucial ingredient: Energy Conservation.

The Missing Ingredient:
Imagine you have a fixed amount of money (energy) to spend on fireworks. If you buy a huge, expensive firework, you can't buy as many small ones. The old math didn't strictly enforce this budget. It let the particles "borrow" energy in ways that don't happen in real life.

The New Solution: The "Modified" Recipe

The authors of this paper decided to fix the math. They took the old sketch and added a strict budget constraint (Energy Conservation). They call this new method MDLA (Modified Double Logarithmic Approximation).

The Result: When they added the budget rule, the predicted shape of the particle distribution changed dramatically. It stopped looking like a blurry sketch and started looking like a sharp, detailed photograph.
The Analogy: It's like switching from a blurry, low-resolution photo of a crowd to a high-definition 4K image. The new math finally matches what the universe actually does.

The "Jet Topics" Trick

Here is the tricky part: In a real collision, you don't just get one type of firework. You get a mix of Quark Jets (like a specific type of firework) and Gluon Jets (a different, more chaotic type). They are mixed together in the detector, making it hard to tell which is which.

To solve this, the authors used a clever statistical trick called "Jet Topics."

The Analogy: Imagine you have a bag of red and blue marbles mixed together. You can't see the colors individually. But, if you take two different bags of mixed marbles (one from a factory that makes mostly red, one from a factory that makes mostly blue) and compare the ratios of colors in specific spots, you can mathematically "unmix" them.
The Application: They used this method on real data from ATLAS to separate the "Quark" jets from the "Gluon" jets without needing to look at every single particle.

The Big Discovery

When they compared their new, budget-conscious math (MDLA) against the real data:

It Worked: The new theory matched the experimental data almost perfectly across a huge range of energies.
The "Universal" Shape is Real: They confirmed that yes, Quark Jets and Gluon Jets do follow their own unique, universal shapes (KNO scaling), but only if you use the right math that respects the energy budget.
Simulation Check: They also checked their math against computer simulations (PYTHIA), and it matched there too.

Why Does This Matter?

This isn't just about counting particles. It's about understanding the fundamental rules of the universe.

Validation: It proves that our understanding of how energy turns into matter is solid, provided we do the math correctly.
Future Tools: By understanding these "shapes" of jets, physicists can build better detectors and algorithms to spot new, rare particles in the future. It's like learning the exact sound of a violin so you can hear it clearly even in a noisy orchestra.

In short: The authors fixed a broken math formula by adding a "budget rule," used a clever statistical trick to separate mixed-up particle types, and proved that the universe follows a beautiful, predictable pattern when you look at it with the right tools.

1. Problem Statement

The paper addresses the theoretical description of particle multiplicity distributions ( $P(n)$ ) in Quantum Chromodynamics (QCD) jets. While the Koba–Nielsen–Olesen (KNO) scaling hypothesis suggests that at asymptotically high energies, the distribution becomes universal when rescaled by the mean multiplicity ( $\bar{n}P(n) = \Psi(n/\bar{n})$ ), standard perturbative QCD calculations face significant challenges:

Failure of DLA: The standard Double Logarithmic Approximation (DLA) predicts KNO scaling functions that differ substantially from experimental data and Monte Carlo simulations (e.g., PYTHIA), particularly failing to describe the shape and peak position of the distributions.
Inclusive vs. Exclusive: In proton-proton ($pp$) collisions, inclusive measurements mix quark- and gluon-initiated jets, obscuring the specific scaling behaviors of each flavor.
Energy Conservation: The standard DLA neglects strict energy conservation in the parton cascade, which is crucial for accurate phenomenological predictions at LHC energies.

The authors aim to derive improved KNO scaling functions by incorporating energy conservation into the DLA framework (Modified DLA or MDLA) and to test these against LHC data (ATLAS) and advanced jet substructure techniques (jet topics).

2. Methodology

A. Theoretical Framework: Modified DLA (MDLA)

Generating Functions (GF): The authors utilize the GF method to describe the evolution of parton multiplicities. They start with the evolution equation for the GF, $Z_a(u, Q)$ , which governs the probability of finding $n$ particles.
Incorporating Energy Conservation: Unlike standard DLA, which approximates splitting functions as $\hat{P} \approx 2C/z$ , the MDLA enforces energy conservation in the NLL evolution equation. This modifies the recursive relations for multiplicity moments.
Derivation of Scaling Functions:
- The authors calculate the moments $f^{(k)}_a$ of the KNO scaling function $\Psi_a(x)$ iteratively.
- They express the scaling functions as an expansion in Laguerre polynomials: $\Psi_a(x) \approx e^{-x} \sum \hat{f}^{(m)}_a L_m(x)$ .
- Parameters: The gluon scaling function depends on the anomalous dimension $\gamma_0$ . The quark scaling function depends on $\gamma_0$ and the ratio of mean multiplicities $c_q = \bar{n}_q / \bar{n}_g$ in the asymptotic limit.
- Fitting: Instead of fixing $c_q$ to the DLA value ( $C_F/C_A$ ), the authors treat it as a phenomenological parameter fitted to data/simulations.

B. Phenomenological Application

Inclusive Distributions: The authors construct the inclusive multiplicity distribution for $pp$ collisions as a weighted sum of quark and gluon distributions: $P(n) = r_q P_q(n) + r_g P_g(n)$ , where $r_{q,g}$ are the fractions derived from Leading Order (LO) cross-sections.
Mean Multiplicity Inputs: To compare with charged-particle data, they combine the MDLA scaling functions with mean multiplicities calculated at N3LO (Next-to-Next-to-Next-to-Leading Order) and the Local Parton–Hadron Duality (LPHD) hypothesis to bridge the parton-hadron gap.
Jet Topics: To isolate quark and gluon jets without relying on theoretical cross-section fractions, they apply the jet topics method (a statistical unmixing technique) to ATLAS data. This separates the "forward" and "central" jet samples into two pure "topics" (approximating quark and gluon jets).

C. Data and Simulation

Experimental Data: ATLAS measurements at $\sqrt{s} = 13$ TeV for the two leading jets with $p_T \in [0.1, 2.5]$ TeV.
Simulation: PYTHIA 8 (A14 tune) with MPI and hadronization, used to validate the MDLA parameters and as a benchmark.

3. Key Contributions

Derivation of MDLA KNO Functions: The paper provides explicit analytical expressions for KNO scaling functions for both quark and gluon jets that incorporate energy conservation. These functions differ qualitatively from standard DLA predictions, shifting the peak toward higher multiplicities.
Phenomenological Success of MDLA: By fixing two parameters ( $\gamma_0 \approx 0.43$ and $c_q \approx 0.8$ ) based on PYTHIA, the MDLA framework successfully describes the inclusive charged-particle multiplicity distributions measured by ATLAS across a wide $p_T$ range.
Validation via Jet Topics: The authors demonstrate that the MDLA scaling functions are consistent with quark- and gluon-dominated distributions extracted from ATLAS data using the jet topics method. This provides a direct test of KNO universality for specific jet flavors without relying on LO cross-section fractions.
Comparison of Mean Multiplicities: The study highlights that while DLA mean multiplicities deviate from data (underestimating at low $p_T$ , overestimating at high $p_T$ ), N3LO calculations combined with MDLA scaling functions provide a precise description of the data.

4. Results

KNO Scaling Verification:
- DLA Failure: Standard DLA scaling functions fail to describe the shape of the distributions, particularly the peak position and width.
- MDLA Success: The MDLA scaling functions (solid lines in Fig. 1 and 2) align well with PYTHIA simulations and the QCD-inspired expressions from recent literature, though with slight quantitative differences in the quark sector.
Inclusive Distributions (Fig. 5 & 6):
- Theoretical predictions using MDLA + N3LO mean multiplicities match ATLAS data ( $\bar{n}_{ch}P(n_{ch})$ ) across all $p_T$ bins (0.1–2.5 TeV) within experimental uncertainties.
- Predictions using DLA mean multiplicities show significant discrepancies, especially at low and high $p_T$ .
Jet Topics Analysis (Fig. 7 & 9):
- The multiplicity distributions extracted for "Topic 1" (quark-like) and "Topic 2" (gluon-like) from ATLAS data are consistent with the MDLA predictions within propagated uncertainties.
- While experimental uncertainties (statistical and systematic) are sizable, the central values support the universality of the KNO scaling functions for distinct jet flavors.
- A slight discrepancy is noted for gluon jets at high $p_T$ ($1.6-2.0$ TeV), attributed to differences in mean multiplicity between PYTHIA and experimental extractions.

5. Significance and Outlook

Theoretical Validation: The work confirms that energy conservation is a critical ingredient for describing multiplicity distributions in QCD jets, validating the MDLA approach over pure DLA.
Experimental Utility: The study demonstrates that jet topic modeling is a viable experimental tool for isolating quark and gluon jet properties, allowing for direct tests of perturbative QCD predictions without relying on theoretical jet fractions.
Future Directions:
- The authors note that a full calculation beyond DLA (including running coupling effects and higher-order corrections) is still needed.
- Improving the precision of experimental measurements (finer binning, reduced uncertainties) is necessary to resolve discrepancies in the high-multiplicity tails and high- $p_T$ regions.
- The work sets a foundation for using infrared-safe multiplicity definitions and more sophisticated quark-gluon discrimination methods in future high-precision QCD studies.

In conclusion, this paper establishes a robust theoretical framework (MDLA) that, when combined with modern jet substructure techniques (jet topics) and high-order mean multiplicity calculations (N3LO), successfully describes the complex multiplicity fluctuations observed in high-energy $pp$ collisions at the LHC.