Polyakov–KNO Multiplicity Distributions in QCD

• The paper introduces a QCD-motivated expression for Polyakov–KNO multiplicity distributions, where the full energy dependence is absorbed into the mean multiplicity.
• It employs generating functions and Laplace transforms with MDLA resummation to derive moments and the overall shape of the scaling function for jets.
• The findings reveal clear scaling violations at finite energies and distinguish between quark and gluon jet behaviors, matching experimental data.

Polyakov–KNO multiplicity distributions are continuous scaling distributions $\Psi(\nu,Q)$ associated with the discrete multiplicity probabilities $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$, defined through $N(Q)P_n(Q)\simeq \Psi(\nu,Q)$ with $N(Q)\equiv \langle n\rangle(Q)$ and $\nu=n/N(Q)$. In the strict KNO limit, $\Psi$ becomes $Q$-independent; in QCD jets one instead expects approximate P-KNO scaling with slow violations controlled by the running coupling. A universal QCD-motivated expression has been proposed for hadrons in jets, and its moments and overall shape in full events and quark and gluon jets are described with reasonable quantitative precision over a range of energies; in particular, the scaling violation predicted by QCD is seen clearly in the moments and high-multiplicity fluctuations (Dokshitzer et al., 10 Jul 2025).

1. Scaling definition and kinematic variables

The hardness scale $Q$ is the organizing variable for the jet multiplicity problem. For $e^+e^-$ at centre-of-mass $\sqrt{s}$, one has $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$0 for full events or $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$1 per hemisphere; more generally, for a jet of energy $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$2 and opening angle $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$3, $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$4. The Polyakov–KNO scaling variable is

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$5

and the continuous scaling distribution is defined by

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$6

The P-KNO hypothesis, identified in the jet analysis as Polyakov 1971, states that the full energy dependence of $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$7 is absorbed into $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$8, so that $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$9. In this formulation, strict scaling would require $N(Q)P_n(Q)\simeq \Psi(\nu,Q)$0, whereas finite-energy QCD predicts scaling violation through the residual $N(Q)P_n(Q)\simeq \Psi(\nu,Q)$1-dependence of the scaling function.

The physical picture is a cascading QCD parton shower. Hadron multiplicities in jets arise from this cascade, whose universal scaling properties were first anticipated by Polyakov in a conformal field theory context. This places Polyakov–KNO distributions at the intersection of multiparticle phenomenology, parton-cascade dynamics, and asymptotic scaling theory (Dokshitzer et al., 10 Jul 2025).

2. Universal jet distribution in the Dokshitzer–Webber framework

For a gluon-initiated jet, the proposed universal form is

$N(Q)P_n(Q)\simeq \Psi(\nu,Q)$2

Its ingredients are specified by the multiplicity anomalous dimension $N(Q)P_n(Q)\simeq \Psi(\nu,Q)$3, the Polyakov exponent $N(Q)P_n(Q)\simeq \Psi(\nu,Q)$4, the tail-normalisation factor

$N(Q)P_n(Q)\simeq \Psi(\nu,Q)$5

the high-$N(Q)P_n(Q)\simeq \Psi(\nu,Q)$6 variable

$N(Q)P_n(Q)\simeq \Psi(\nu,Q)$7

and the slowly-varying prefactor

$N(Q)P_n(Q)\simeq \Psi(\nu,Q)$8

The anomalous dimension is given to two loops by

$N(Q)P_n(Q)\simeq \Psi(\nu,Q)$9

with $N(Q)\equiv \langle n\rangle(Q)$0 in the CMW (physical) scheme.

For a more general source $N(Q)\equiv \langle n\rangle(Q)$1, such as a quark jet or a full $N(Q)\equiv \langle n\rangle(Q)$2 event, the distribution is written in terms of the hadron-production “power” $N(Q)\equiv \langle n\rangle(Q)$3: $N(Q)\equiv \langle n\rangle(Q)$4 where $N(Q)\equiv \langle n\rangle(Q)$5 is a slowly varying function of $N(Q)\equiv \langle n\rangle(Q)$6 and $N(Q)\equiv \langle n\rangle(Q)$7. Apart from two ad hoc constants $N(Q)\equiv \langle n\rangle(Q)$8, the entire shape and its energy dependence are determined by one function $N(Q)\equiv \langle n\rangle(Q)$9, itself fixed by the two-loop running coupling; quark/gluon differences enter only via a multiplicative $\nu=n/N(Q)$0 factor (Dokshitzer et al., 10 Jul 2025).

3. Resummed derivation from moments to the full distribution

The construction begins with the P-KNO hypothesis and then uses generating functions and Laplace transforms. In the leading double-logarithmic approximation (DLA), one obtains an asymptotic form for the Laplace transform $\nu=n/N(Q)$1 and therefore for the factorial and power moments $\nu=n/N(Q)$2. For $\nu=n/N(Q)$3, the gluon-jet moments behave as

$\nu=n/N(Q)$4

The modified DLA (MDLA), identified with Dokshitzer ’93 in the summary, resums all terms of order $\nu=n/N(Q)$5 in the perturbative series and restores exact energy-momentum balance in the cascade. This yields

$\nu=n/N(Q)$6

Laplace inversion by steepest descent reconstructs $\nu=n/N(Q)$7 from $\nu=n/N(Q)$8, matching the large-$\nu=n/N(Q)$9 behaviour of the moments to the moments of the scaling function. The resulting high-$\Psi$0 tail is

$\Psi$1

A two-parameter completion then introduces subleading shifts $\Psi$2 so as to restore positivity for $\Psi$3 and exactly satisfy $\Psi$4. The minimal choice $\Psi$5 fixes the full-range ansatz. Extension to quark jets and full $\Psi$6 events follows from the DLA relation

$\Psi$7

which is then re-inverted to obtain $\Psi$8. The overall procedure moves from moments to a full scaling distribution while preserving the perturbatively controlled high-multiplicity structure (Dokshitzer et al., 10 Jul 2025).

4. Moments, cumulants, and QCD scaling violation

The mean multiplicity is determined by

$\Psi$9

and the power moments are defined by

$Q$0

In MDLA, the gluon moments retain the form

$Q$1

For a quark jet $Q$2 or a full $Q$3 event $Q$4, the corresponding moments are obtained by repeated differentiation of $Q$5.

Cumulants follow in the usual way from the $Q$6. In particular, the dispersion $Q$7 satisfies

$Q$8

At low $Q$9, DLA fails because $Q$0 and hence $Q$1, whereas MDLA cures this for $Q$2 and predicts a strong KNO-violation for $Q$3 as $Q$4. The energy dependence of the moments is driven entirely by $Q$5: as $Q$6 increases, $Q$7 decreases slowly, the MDLA suppression factor tends to $Q$8, and the asymptotic DLA values are recovered.

Typical moments extracted from the summary are:

$Q$9	$e^+e^-$0: $e^+e^-$1	$e^+e^-$2: $e^+e^-$3
2	$e^+e^-$4	$e^+e^-$5
3	$e^+e^-$6	$e^+e^-$7
4	$e^+e^-$8	$e^+e^-$9

These values were read off Fig. 1 with uncertainties $\sqrt{s}$0. A common misconception is that KNO scaling in QCD should be exact at accessible energies. In the jet calculation, the asymptotic regime is recovered only slowly because $\sqrt{s}$1 only logarithmically with $\sqrt{s}$2; at finite $\sqrt{s}$3, the high-$\sqrt{s}$4 region exhibits noticeable scaling violation, including the narrowing of the distribution as $\sqrt{s}$5 increases (Dokshitzer et al., 10 Jul 2025).

5. Comparison with $\sqrt{s}$6 events and identified jet samples

The distributions are compared in the form $\sqrt{s}$7 versus $\sqrt{s}$8, in both linear and logarithmic scale. For full $\sqrt{s}$9 events at LEP-1, with $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$00, ALEPH, OPAL, DELPHI, and L3 all agree, and MDLA with $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$01 reproduces the shape and high-$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$02 tail to within experimental errors. OPAL data at LEP-2, $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$03 and $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$04, were also compared with the MDLA form. At lower energies, including TASSO at $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$05, HRS at $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$06, ARGUS at $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$07, and TASSO at $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$08 and $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$09, the distributions are broader than MDLA for $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$10, consistent with the proximity to $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$11.

For quark jets, DELPHI at $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$12 and HRS at $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$13 show that the MDLA curves for $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$14 and, optionally, $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$15, bracket the data within uncertainties. For gluon jets, OPAL unbiased gluon jets in three-jet events, with $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$16–$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$17 and mean $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$18, corresponding to $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$19, prefer the $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$20 gluon curve over the $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$21 quark curve, especially in the tail $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$22. OPAL gluon jets recoiling against $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$23, with $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$24 and restricted rapidity $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$25, show reasonable agreement at mid-rapidity.

The high-multiplicity tail is central to the phenomenology. In pure KNO scaling one would have $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$26 as $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$27. In QCD, however, the tail falls faster than exponentially: $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$28 As $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$29, $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$30, so $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$31 and $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$32, giving $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$33 and the pure exponential tail $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$34. At finite $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$35, the super-exponential suppression is therefore stronger, and the narrowing of the distribution in the high-$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$36 region becomes a direct manifestation of scaling violation (Dokshitzer et al., 10 Jul 2025).

6. Other QCD realizations, limitations, and interpretation

Polyakov–KNO-type behaviour also appears in other high-energy QCD limits. In the BFKL approach, the multiplicity generating function for real gluon emissions in a cut ladder takes the form

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$37

with mean multiplicity

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$38

in LLA, and therefore

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$39

This is exactly a Poisson distribution. Defining $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$40 and $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$41, Stirling’s formula gives an asymptotic form in which, up to a slowly-varying prefactor $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$42, the exponent depends only on $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$43. The analysis stresses that this Poisson law describes the number of BFKL “cells” in the cut diagram rather than fully dressed physical gluons, and that realistic cut-offs or restricted rapidity windows produce modest corrections while leaving the leading Poisson behaviour intact (Arakelyan et al., 2020).

A distinct realization arises in deep inelastic scattering in the CGC/saturation framework at large $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$44. In that regime one finds $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$45 with $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$46 at LO BFKL, and for $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$47 the multiplicity distribution almost reproduces KNO scaling with

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$48

For $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$49, however, the distribution is instead

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$50

and these small-$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$51 terms determine the entropy

$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$52

at large $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$53, where the factor $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$54 stems from non-perturbative corrections associated with the large-$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$55 behaviour of the saturation momentum. The derivation relies on LO pQCD with fixed coupling, geometric scaling, AGK cutting rules for CGC Pomeron fan diagrams, and a smooth nonperturbative decay of $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$56; running-coupling effects, NLO corrections, and the model for $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$57 affect both the small-$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$58 region and the entropy coefficient (Levin, 2023).

Taken together, these results indicate that Polyakov–KNO scaling in QCD is not a single universal law with identical microscopic content in all regimes. In jets, the salient feature is a resummed, largely parameter-free analytic form whose deviations from strict scaling are governed by $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$59. In BFKL ladders, the leading multiplicity law is Poissonian. In deep saturation, the large-$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$60 tail is exponential in $P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$61, while the small-$P_n(Q)=\sigma_n(Q)/\sigma_{\rm tot}(Q)$62 sector obeys a different law. This suggests that the Polyakov–KNO concept is best understood as a scaling framework whose concrete realization is determined by the underlying QCD dynamics, the relevant hardness variable, and the treatment of energy conservation, rapidity restrictions, and impact-parameter structure.