2025-11-18T00:28:13.229814

QCD-inspired description of multiplicity distributions in jets

Dokshitzer, Webber
We suggest a universal QCD-motivated expression for the Polyakov-KNO multiplicity distributions of hadrons in jets and compare it with data from $e^+e^-$ annihilation experiments. The moments and overall shape of the distributions in full events and quark and gluon jets, over a range of energies, are described with reasonable quantitative precision. In particular, the scaling violation predicted by QCD is seen clearly in the moments and high-multiplicity fluctuations.
academic

QCD-inspired description of multiplicity distributions in jets

Basic Information

  • Paper ID: 2507.07691
  • Title: QCD-inspired description of multiplicity distributions in jets
  • Authors: Yu.L. Dokshitzer (Riga Technical University), B.R. Webber (University of Cambridge)
  • Classification: hep-ph (High Energy Physics - Phenomenology)
  • Publication Date: October 15, 2025 (v2)
  • Paper Link: https://arxiv.org/abs/2507.07691v2

Abstract

This paper proposes a universal expression based on QCD theory to describe the Polyakov-KNO multiplicity distributions of hadrons in jets, and compares it with experimental data from e+ee^+e^- annihilation. The method successfully describes the moments and overall shape of distributions in complete events as well as in quark and gluon jets across different energy ranges with reasonable precision. Notably, the scaling violation predicted by QCD is clearly manifested in the multiplicity moments and high-multiplicity fluctuations.

Research Background and Motivation

Research Problem

This study addresses the description of multiplicity distributions in multi-particle production in strong interactions, particularly the statistical distribution laws of hadron multiplicities within jets in hard processes.

Importance of the Problem

  1. Theoretical Significance: Multiplicity distributions provide an important window for understanding QCD dynamics, particularly parton cascade processes
  2. Experimental Verification: Provides quantitative tests of QCD theoretical predictions in hadronization processes
  3. Phenomenological Application: Crucial for understanding jet physics in high-energy experiments such as the LHC

Limitations of Existing Methods

  1. DLA Approximation: Traditional double-logarithmic approximation (DLA) exhibits poor convergence at finite coupling constants
  2. Parameter Dependence: Existing models often rely on excessive phenomenological parameters
  3. Energy Scale Dependence: Lack of unified description of scaling violation at different energy scales

Research Motivation

Building on the authors' previous work 1, this paper aims to extend the parameter-independent description of high-multiplicity fluctuations to the complete multiplicity distribution, establishing a universal framework based on perturbative QCD.

Core Contributions

  1. Theoretical Framework: Develops the Modified Double-Logarithmic Approximation (MDLA) method, systematically resumming terms in the perturbative QCD series proportional to the order of multiplicity moments
  2. Universal Formula: Proposes analytical expressions describing the complete P-KNO distribution applicable to different types of jet sources
  3. Experimental Verification: Systematic comparison with e+ee^+e^- annihilation data in the energy range from 29 GeV to 189 GeV
  4. Scaling Violation: First clear observation of QCD-predicted P-KNO scaling violation

Methodology Details

Task Definition

The research objective is to describe the probability distribution of hadron multiplicities in hard processes: Pn(Q)=σn(Q)σtot(Q)P_n(Q) = \frac{\sigma_n(Q)}{\sigma_{tot}(Q)} where QQ is the hardness scale of the process and nn is the number of produced hadrons.

P-KNO Phenomenon

Introducing the scaling variable ν=n/N(Q)\nu = n/N(Q), where N(Q)=n(Q)N(Q) = \langle n \rangle(Q) is the average multiplicity, the P-KNO scaling hypothesis states: N(Q)Pn(Q)Ψ(ν,Q)KNOΨ(ν)N(Q) \cdot P_n(Q) \equiv \Psi(\nu, Q) \xrightarrow{KNO} \Psi(\nu)

MDLA Method Core

Correction of Multiplicity Moments

For gluon jets, MDLA corrects the DLA moments as: gk,MDLA(g)=gk,DLA(g)[Γ(1+γ)]kΓ(1+γk)g_{k,MDLA}^{(g)} = g_{k,DLA}^{(g)} \cdot \frac{[\Gamma(1+\gamma)]^k}{\Gamma(1+\gamma k)}

where γ(αs(Q))\gamma(α_s(Q)) is the anomalous dimension of multiplicity: γ(αs)=2Ncαsπ(β04+10nf3Nc2)αs2π\gamma(α_s) = \sqrt{2N_c}\frac{α_s}{\pi} - \left(\frac{β_0}{4} + \frac{10n_f}{3N_c^2}\right)\frac{α_s}{2\pi}

Complete Form of P-KNO Distribution

For gluon jets, the following expression is proposed: νΨ(ν)2μ2κ2eκ[1+(p1κ)1]p1χ(ksd)\nu\Psi(\nu) \simeq 2\mu^2 \frac{\kappa^2 e^{-\kappa}}{[1+(p_1\kappa)^{-1}]^{p_1}} \chi(k_{sd})

where:

  • κ=[D(γ)ν]μ(γ)\kappa = [D(\gamma)\nu]^{\mu(\gamma)}
  • μ(γ)=11γ\mu(\gamma) = \frac{1}{1-\gamma}
  • D(γ)=Cγγ(1γ)1γΓ(1+γ)D(\gamma) = \frac{C\gamma^\gamma(1-\gamma)^{1-\gamma}}{\Gamma(1+\gamma)}
  • ksd=μκ[1+2(p0κ)1]p0k_{sd} = \frac{\mu\kappa}{[1+2(p_0\kappa)^{-1}]^{p_0}}

Generalization to Arbitrary Sources

For an arbitrary source SS with hadron production capacity ρS\rho_S: νΨ(S)(ν)[νΨ(ν)]ρρGρ1\nu\Psi^{(S)}(\nu) \simeq [\nu\Psi(\nu)]^{\rho} \cdot \sqrt{\rho}G^{\rho-1} where G=2πμ2/[γκ+2(1γ)]G = \sqrt{2\pi\mu^{-2}/[\gamma\kappa + 2(1-\gamma)]}.

Technical Innovations

  1. Energy Balance: MDLA systematically accounts for energy balance in parton cascades, which DLA neglects
  2. All-Order Resummation: Resums all terms of the form (αsk)m(\sqrt{α_s} \cdot k)^m
  3. Positivity Correction: Resolves positivity issues in the low-multiplicity region through parameters p0,p1p_0, p_1
  4. Scale Dependence: All scale dependence is embodied through a single parameter γ(αs(Q))\gamma(α_s(Q))

Experimental Setup

Datasets

Multiple e+ee^+e^- annihilation experimental data are employed:

  1. LEP-1: ALEPH, OPAL, DELPHI, L3 at s=MZ=91.2\sqrt{s} = M_Z = 91.2 GeV
  2. LEP-2: OPAL at s=161,189\sqrt{s} = 161, 189 GeV
  3. Low-energy Data: TASSO (s=44\sqrt{s} = 44 GeV), HRS (s=29\sqrt{s} = 29 GeV)
  4. Jet Data: OPAL quark and gluon jet data

Evaluation Metrics

  1. Multiplicity Moments: gk=nk/[N(Q)]kg_k = \langle n^k \rangle / [N(Q)]^k
  2. Distribution Shape: Complete P-KNO distribution Ψ(ν)\Psi(\nu)
  3. Scaling Ratio: N/D=n/n2n2N/D = \langle n \rangle / \sqrt{\langle n^2 \rangle - \langle n \rangle^2}

Parameter Settings

  • Quark hadron production capacity: ρq=2/3\rho_q = 2/3
  • Model parameters: p0=1,p1=4p_0 = 1, p_1 = 4 (energy and process independent)
  • QCD coupling: Physical scheme with αs(MZ)=0.127α_s(M_Z) = 0.127

Experimental Results

Main Results

Precise Description of Low-Order Moments

MDLA significantly improves the description of low-order multiplicity moments. At LEP energies, the agreement between theoretical predictions and experimental data reaches precision at the level of several percent.

Complete Distribution Fitting

  1. LEP-1 Data: At Q=MZQ = M_Z, the MDLA formula shows good agreement with L3 and other experimental data across the entire multiplicity range
  2. Energy Dependence: Theoretical predictions reasonably describe experimental data across the wide energy range from 29 GeV to 189 GeV
  3. Jet Types: Reasonable precision is achieved for both quark jets (single hemisphere) and gluon jets

Observation of Scaling Violation

Clear observation of QCD-predicted P-KNO scaling violation for the first time:

  • At lower energy scales, larger QCD coupling leads to reduced multiplicity moments
  • Distributions become narrower with faster decline in high-multiplicity tails
  • Comparison between HRS (29 GeV) and LEP-1 (91 GeV) data clearly demonstrates this effect

Ablation Studies

Stability of theoretical predictions is verified by comparing different choices of hardness scales (Q=sQ = \sqrt{s} vs Q=s/2Q = \sqrt{s}/2).

Case Studies

  1. Gluon Jets: OPAL gluon jet data tend to support predictions with ρ=1\rho = 1
  2. Restricted Rapidity: Hadron production in restricted rapidity intervals shows strong interval dependence

Historical Development

  1. KNO Scaling: Koba-Nielsen-Olesen (1972) scaling phenomenon in multiplicity distributions
  2. Polyakov Theory: A.M. Polyakov (1970) predictions within conformal field theory framework
  3. QCD Development: Bassetto, Ciafaloni, Marchesini (1983) theoretical foundations based on QCD cascades

Theoretical Methods

  1. DLA Method: Limitations of traditional double-logarithmic approximation
  2. Monte Carlo Methods: Phenomenological models based on parton shower
  3. MLLA: Modified leading-logarithmic approximation applied to average multiplicity

Conclusions and Discussion

Main Conclusions

  1. Theoretical Success: MDLA method successfully connects perturbative QCD with experimental observations of hadron multiplicity distributions
  2. Universality: The proposed formula applies to different types of jet sources and wide energy ranges
  3. Scaling Violation: Clear experimental observation of QCD-predicted scaling violation for the first time

Limitations

  1. Low-Multiplicity Region: Lack of theoretical control in the ν<1\nu < 1 region
  2. Non-Perturbative Effects: Non-perturbative effects in hadronization require further understanding
  3. Parameter ρ\rho: The hadron production capacity parameter requires better theoretical foundation

Future Directions

  1. LHC Applications: Extension of the method to LHC jet physics
  2. Gluon Jets: More experimental data on gluon jets needed
  3. Restricted Rapidity: Deeper investigation of multiplicity distributions in restricted rapidity intervals
  4. Non-Perturbative Theory: Development of theoretical methods for low-multiplicity fluctuations

In-Depth Evaluation

Strengths

  1. Theoretical Rigor: Systematic method based on perturbative QCD avoiding excessive phenomenological parameters
  2. Experimental Verification: Detailed comparison with extensive experimental data validates the method's effectiveness
  3. Predictive Power: Successfully predicts and observes QCD scaling violation
  4. Universality: Unified description of different jet types and energy ranges

Weaknesses

  1. Theoretical Deficiency: Lack of first-principles theoretical description in the low-multiplicity region
  2. Parameter Dependence: Although fewer parameters, the choice of p0,p1p_0, p_1 retains some arbitrariness
  3. Applicable Range: Theoretical singularities exist near very low energy scales (Q5Q \lesssim 5 GeV)

Impact

  1. Theoretical Contribution: Provides new theoretical tools for QCD applications in soft hadron physics
  2. Experimental Guidance: Establishes theoretical foundation for future LHC jet research
  3. Methodological Value: MDLA method may be applicable to other QCD phenomenology problems

Application Scenarios

  1. High-Energy Physics Experiments: Particularly suited for e+ee^+e^- annihilation and LHC jet physics
  2. Theoretical Research: Provides powerful tools for studying QCD dynamics and parton cascades
  3. Phenomenological Analysis: Applicable to analyzing multi-particle production in various hard processes

References

This paper cites 24 important references, including:

  • Classical KNO scaling and Polyakov theory literature
  • Key papers in QCD theoretical development
  • Multiplicity measurement results from LEP, TASSO, HRS and other experiments
  • Authors' previous work on high-multiplicity fluctuations 1

Overall Assessment: This is a high-quality theoretical physics paper that successfully connects perturbative QCD theory with experimental observations, achieving important progress in the description of multiplicity distributions. Although theoretical deficiencies remain in the low-multiplicity region, its universality and predictive power make it a significant contribution to the field.