1D Mueller Dipole Model in High-Energy QCD

• The 1D Mueller dipole model is a simplified QCD cascade that uses rapidity as its sole variable to model gluon emissions and predict multiplicity distributions.
• It derives an exact geometric solution with exponential mean multiplicity growth and logarithmic entropy scaling in the high energy, small-x regime.
• The model’s generalized negative binomial extension, introducing a conformal weight, improves data fits and links the approach to broader saturation physics.

The 1D Mueller dipole model is a simplified QCD parton cascade model in which the splitting dynamics depend solely on rapidity and not on the dipole size; in that sense, “1D” denotes a reduction of the full dipole picture to a rapidity-only stochastic branching process. It provides a minimal description of gluon emissions in the high energy, small-$x$ limit and has been used to study dipole and hadron multiplicity distributions, Shannon entropy, and their asymptotic scaling properties in proton-proton collisions. Recent work has examined both the original model, its high energy limit, and a generalized version motivated by studies of 1D systems with conformal symmetry, with particular emphasis on the observable $S(\ln\langle n\rangle)$ for model–data comparison (Kutak et al., 19 Apr 2026, Kutak et al., 9 Sep 2025).

1. Conceptual setting and physical interpretation

In the formulation used in recent phenomenological studies, the model describes the evolution of partonic cascades in rapidity, focusing on gluon emissions in the high energy limit. The dynamical variable is the probability $P_n(y)$ to find $n$ dipoles at rapidity $y$. The parameter $\alpha$ is the emission kernel characterizing the growth rate, related to BFKL evolution; in the BFKL saddle point approximation it is often identified as $4\bar\alpha_s\ln 2$. The stochastic process is a pure branching process: dipoles split, but there is no merging or recombination in the original 1D construction. In statistical language, the dynamics are similar to a birth process of Yule or geometric type (Kutak et al., 9 Sep 2025).

This reduction has two consequences that define the scope of the model. First, it isolates the rapidity-driven growth of multiplicity. Second, it does not encode the full transverse-coordinate dynamics of the original Mueller dipole formalism. A common misconception is therefore to treat the 1D model as a complete high-energy QCD description; the published analyses instead use it as a minimal QCD-inspired cascade whose chief utility lies in deriving closed-form multiplicity and entropy relations and testing them against data.

2. Evolution equation and exact solution

The probabilistic cascade equation for the dipole multiplicity distribution is

$$\partial_y P_n(y) = -\alpha n P_n(y) + (n-1)\alpha P_{n-1}(y).$$

The loss term describes the depletion of the $n$-dipole sector by splitting, while the gain term accounts for transitions from $n-1$ to $S(\ln\langle n\rangle)$0 dipoles. Introducing an overall normalization parameter $S(\ln\langle n\rangle)$1, the solution is

$S(\ln\langle n\rangle)$2

With

$S(\ln\langle n\rangle)$3

this becomes the geometric distribution

$S(\ln\langle n\rangle)$4

with support on $S(\ln\langle n\rangle)$5 (Kutak et al., 19 Apr 2026).

The mean multiplicity follows immediately: $S(\ln\langle n\rangle)$6 This exponential growth of $S(\ln\langle n\rangle)$7 with rapidity is the characteristic high-energy behavior of the model. Because the solution is geometric, the original 1D Mueller model predicts a relatively narrow multiplicity distribution; this feature becomes central in later comparisons with proton-proton data.

3. Entropy and the high-energy limit

The Shannon entropy of the multiplicity distribution is defined as

$S(\ln\langle n\rangle)$8

For the geometric distribution generated by the 1D Mueller model, the entropy can be written as

$S(\ln\langle n\rangle)$9

and in the large-$P_n(y)$0 limit one obtains

$P_n(y)$1

This limiting form is treated as the universal formula in the high energy, large multiplicity regime (Kutak et al., 9 Sep 2025).

As $P_n(y)$2, the mean multiplicity grows exponentially and the entropy grows logarithmically in the mean multiplicity. The published discussion connects this asymptotic behavior to a maximal-entropy, or maximal-entanglement, picture in which microstates become effectively equiprobable in the high-energy limit. At the same time, the asymptotic formula should not be conflated with a full phenomenological description at finite multiplicity: the data analyses show that the geometric cascade captures the limiting trend more reliably than the detailed shape of measured multiplicity distributions.

4. The observable $P_n(y)$3

To address ambiguities arising from different definitions of rapidity or pseudorapidity intervals in experimental measurements, recent analyses propose the entropy as a function of the logarithm of the average multiplicity,

$P_n(y)$4

as a universal observable. The construction is direct. From a measured multiplicity distribution $P_n(y)$5 one computes

$P_n(y)$6

and then forms the pair $P_n(y)$7 for comparison to model predictions. This observable was built from proton-proton data across $P_n(y)$8 from 200 GeV to 13 TeV and over a wide rapidity range (Kutak et al., 19 Apr 2026).

The significance of this choice is methodological. Rather than comparing observables tied to a single rapidity-window convention, the analysis compares entropy and multiplicity after projecting both onto the same functional relation. This suggests a way to reduce acceptance-dependent ambiguities when confronting cascade models with heterogeneous experimental datasets. The proposal is not that $P_n(y)$9 removes all experimental systematics, but that it provides a more robust cross-experiment comparison than raw multiplicity measures alone.

5. Generalization to a negative binomial cascade

A generalization of the 1D Mueller model was introduced following studies of 1D systems with conformal symmetry. The extended cascade equation is

$n$0

where $n$1 is an additional parameter, identified in the phenomenological discussion as a conformal weight. In the standard Mueller model, $n$2. The solution of the generalized equation is a negative binomial distribution,

$n$3

with mean multiplicity

$n$4

and dispersion parameter

$n$5

This modification allows the width of the distribution to vary independently of the mean, in contrast to the geometric form of the original model (Kutak et al., 19 Apr 2026).

The generalization is phenomenologically consequential because negative binomial distributions are commonly used in fits to proton-proton multiplicity data. In the terms used by the published analyses, the added parameter $n$6 controls the width or dispersion of the cascade. Entropy in the generalized model is more complicated than in the geometric case and does not admit an equally simple explicit closed form; it is instead evaluated numerically or through implicit expressions. A plausible implication is that the generalized model preserves the rapidity-driven branching logic of the 1D cascade while relaxing the overly restrictive fluctuation structure of the pure geometric process.

6. Comparison with proton-proton data

Both the original and generalized cascade models were fitted to $n$7 extracted directly from measured charged-particle multiplicity distributions. In the reported fits, the original Mueller model used $n$8 and $n$9 as fixed parameters and fitted $y$0, while the generalized model fitted $y$1, fixed $y$2, and used $y$3 as fixed. The generalized model provided a significantly better description of the data, with the strongest discrepancy of the original geometric model appearing at lower multiplicity, corresponding in the published discussion to forward rapidity and low-$y$4 regions (Kutak et al., 19 Apr 2026).

Parameter	1D Mueller	Generalized
$y$5	0.50 (fixed)	$y$6
$y$7	0.32 (fixed)	0.32 (fixed)
$y$8	$y$9	3.13 (fixed)
$\alpha$0	309/63	24/63

These results support a clear phenomenological conclusion. The 1D Mueller dipole model gives a qualitative account of mean multiplicity growth and entropy scaling in rapidity at high energies, but it systematically underestimates the width, and therefore the entropy, of measured multiplicity distributions. The data are broader, or more over-dispersed, than the geometric prediction, especially at low and moderate $\alpha$1. The generalized negative-binomial cascade reproduces both the shape and the absolute value of $\alpha$2 much more accurately across the analyzed energy and pseudorapidity range. The reported best-fit value $\alpha$3 is also significantly larger than the standard Mueller value $\alpha$4, reinforcing the conclusion that the original 1D birth process is too restrictive for quantitative hadronic final-state phenomenology.

7. Relation to saturation physics, Pomeron loops, and multiplicity laws

The 1D Mueller dipole model also appears in a broader high-energy framework linking branching dynamics to saturation physics and dipole-dipole scattering. One study connected the exact analytical solution of the Balitsky–Kovchegov equation in the saturation region to the summation of large Pomeron loops within the one-dimensional Mueller dipole model. Using a generating functional

$\alpha$5

the paper derived high-energy parton densities, obtained an $\alpha$6-matrix whose asymptotic suppression matches the rare-fluctuation prediction of Iancu and Mueller, and—via the Abramovsky, Gribov and Kancheli cutting rules—found multiplicity distributions obeying the Koba–Nielsen–Olesen law. In that saturation-region analysis, the entropy was written as

$\alpha$7

with $\alpha$8 taken in the saturation domain (Levin, 2024).

These developments go beyond the original pure-birth formulation used in proton-proton multiplicity fits, but they clarify the model’s role as a scaffold for more elaborate high-energy constructions. The common structure is the same stochastic dipole branching picture. This suggests that the 1D Mueller framework has a dual status in the literature: as a minimal analytic model for entropy and multiplicity scaling, and as a starting point for incorporating unitarization, Pomeron-loop effects, and saturation-region multiplicity laws in dipole-dipole scattering.