One-Dimensional Dipole Cascade Models

• One-dimensional dipole cascade models are analytic frameworks that reduce complex, multi-dimensional evolution to a single directional dynamic for both QCD cascades and lattice dislocation phenomena.
• In QCD, the 1D Mueller dipole and its conformal extensions yield explicit solutions for parton multiplicity and entropy observables, enabling detailed model-to-data comparisons.
• In lattice systems, sequential kink–antikink (dipole) transitions are modeled with closed-form energy solutions, providing insights into nucleation barriers and defect dynamics.

One-dimensional dipole cascade models constitute a tractable framework for describing cascade-type evolution phenomena in both quantum chromodynamics (QCD) and classical lattice systems. Two major axes of research have emerged: QCD-inspired models for parton cascades under rapidity evolution, exemplified by the 1D Mueller dipole model and its conformal-symmetry-motivated generalizations (Kutak et al., 9 Sep 2025); and atomistic models for lattice dislocation or kink–antikink pair cascades, based on discrete nonlinear chains with double-well potentials (Sharma, 2020). These models share the property that the high-dimensional spatial or internal structure is effectively reduced to a single evolution or chain direction, enabling exact solutions, detailed analysis of entropy and multiplicity observables, and rigorous comparison to experiment or simulation.

1. One-dimensional Mueller Dipole Model: Foundations and Solutions

The 1D Mueller dipole model is a stochastic, rapidity-evolution Markov process for parton cascades in QCD, formulated in the limit where all transverse size dependence is neglected and only rapidity $y$ is relevant (Kutak et al., 9 Sep 2025). Each dipole (or gluon) splits independently as $y$ is increased. The strictly linear, non-recombining system is governed by a master equation for $P_n(y)$, the probability of having $n$ dipoles at rapidity $y$: $$\frac{\partial}{\partial y}P_n(y) = -\alpha n P_n(y) + \alpha (n-1) P_{n-1}(y)$$ where $\alpha$ is the splitting kernel, typically set to $4\bar \alpha_s \ln 2$ in BFKL saddle-point approximation. The model assumes no dipole–dipole recombination or saturation effects.

The exact solution is a shifted geometric distribution for $n \geq 1$: $$P_n(y) = \frac{1}{C} e^{-\alpha y} \left(1 - \frac{1}{C} e^{-\alpha y}\right)^{n-1}$$ with normalization $C$. The mean dipole multiplicity evolves exponentially: $\langle n \rangle(y) = C e^{\alpha y}$ In the large-$y$ (high-energy) limit, $C$ is set to unity, yielding a pure geometric distribution.

2. High-Energy Limit and Entropic Observables

Entropy, defined as $S(y) = -\sum_{n=1}^\infty P_n(y)\ln P_n(y)$, characterizes the information content or disorder of the multiplicity distribution. For the geometric distribution,

$$S(y) = \ln\langle n\rangle - (\langle n\rangle-1)\ln(\langle n\rangle-1)$$

In the universal, asymptotic regime ($\langle n\rangle \gg 1$), the leading behavior simplifies to

$$S(\langle n\rangle) \approx \ln\langle n\rangle + 1$$

This form is independent of rapidity window and normalization, motivating $S(\ln\langle n\rangle)$ as a universal observable for final-state hadronic entropy.

3. Conformal-Symmetry Generalization and Negative-Binomial Extension

Extending the 1D dipole cascade model to incorporate conformal symmetry, as motivated by SL(2,ℝ) invariance and Krylov complexity, introduces a conformal weight $h$ (Kutak et al., 9 Sep 2025). The master equation becomes

$$\frac{\partial}{\partial y}P_n(y) = -\alpha (n + 2h) P_n(y) + \alpha (n-1 + 2h) P_{n-1}(y)$$

This generalization leads to a negative-binomial distribution (NBD) for the multiplicity, allowing $n \geq 0$ and parameterized by $k = 2h$: $$P_n(y) = \frac{\Gamma(2h + n)}{n! \Gamma(2h)} p^{2h} (1-p)^n$$ with $p = \frac{1}{C} e^{-\alpha y}$ and mean multiplicity

$$\langle n \rangle(y) = 2h \left( C e^{\alpha y} - 1 \right)$$

The NBD form introduces enhanced variance (“shoulder” in $P_n$), and the vacuum-like term $(n + 2h)$ effectively accounts for subleading, non-BFKL dynamics.

4. Universal Observable $S(\ln\langle n\rangle)$ and Experimental Tests

The physically robust, ambiguity-minimizing observable is $S$ as a function of $x = \ln\langle n\rangle$, where both are defined from any experimental multiplicity distribution $P_n$: $$\langle n\rangle = \sum_n n P_n\quad,\quad S = -\sum_n P_n \ln P_n$$ The conjecture, confirmed empirically, is that the final-state (hadronic) entropy $S_h$ equals the partonic entropy $S$ and that $S(\ln\langle n\rangle)$ is universal across energies and rapidity intervals. This allows model–data comparison that is insensitive to experimental rapidity-window ambiguities.

The practical workflow involves extracting $P_n$ from published $pp$ collision data (ALICE, ATLAS, CMS, LHCb, UA5), computing $\langle n\rangle$ and $S_h$, and plotting $S_h$ versus $\ln\langle n\rangle$. Model predictions (from either distribution) are evaluated for direct comparison.

5. Parameter Fitting and Model Discrimination

Both the original 1D Mueller dipole and the generalized negative-binomial model possess free parameters:

• 1D Mueller: $\alpha$, $C$ (with $h=1/2$ fixed)
• Generalized: $\alpha$ (fixed), $C$ (fixed or fit), $h$

Best-fit parameters to combined $pp$ data across $\sqrt{s} = 0.2$–13 TeV and multiple $\eta$ windows are summarized as:

Parameter	1D Mueller dipole	Generalized model
$h$	0.50 (fixed)	$0.92 \pm 0.05$
$\alpha$	0.32 (fixed)	0.32 (fixed)
$C$	$3.13 \pm 0.48$	3.13 (fixed)
$\chi^2/\mathrm{NDF}$	309/63	24/63

The generalized model achieves tight uncertainty bounds and $\chi^2/\mathrm{NDF}\approx 0.38$, significantly outperforming the $\chi^2/\mathrm{NDF}\approx 4.9$ of the 1D Mueller model. The inclusion of $h$ provides necessary flexibility in the variance structure to match experimental entropy at low and intermediate $\langle n\rangle$.

6. Lattice Dislocation–Dipole (Kink-Antikink) Cascade Models

Beyond QCD, one-dimensional cascade models also arise in the context of discrete lattice systems exhibiting sequential dipole (kink–antikink) formation under applied stress. A prototypical example is an infinite 1D chain with harmonic nearest-neighbor coupling and a two-well (piecewise-quadratic) on-site potential (Sharma, 2020). The equilibrium equations become a discrete nonlinear map: $$u_{n+1} - 2u_n + u_{n-1} - \kappa^2 [w'(u_n) - \sigma] = 0$$ where $\sigma$ is applied stress, $w(u)$ the double-well potential, and $\kappa^2$ a dimensionless parameter.

A “dipole” corresponds to a configuration where sites between $n_l$ and $n_r$ occupy the upper well, and the rest the lower, with separation $\ell = n_r - n_l - 1$. Explicit solutions for $u_n(\ell)$ and total energies $E(\ell;\sigma)$ are available in closed form for the two-well model, as are analytic expressions for energy barriers (Peierls relief), critical nucleation conditions, and the effect of an intermediate “spinodal” region.

Sequential, cascade-like transitions (growth or separation of the dipole) are constructed via path-following procedures that interpolate order parameters and compute corresponding energy landscapes. These results permit analytic study of nucleation barriers, metastability, and the role of lattice trapping, linking directly to phenomena seen in more complex models such as the Frenkel–Kontorova chain.

7. Implications, Extensions, and Outlook

One-dimensional dipole cascade models, in both their QCD and lattice incarnations, provide analytically tractable laboratories for the study of universal cascade phenomena. In the QCD context, the universality of $S(\ln\langle n\rangle)$ supports the maximal entanglement and parton–hadron duality paradigms, and the empirical fit $h \approx 0.9$ points to significant non-BFKL contributions in real cascades (Kutak et al., 9 Sep 2025). The conformal-symmetry generalization via negative-binomial statistics enables accurate reproduction of experimental charged-particle entropy distributions.

Possible directions for further development include the extension to recombination and saturation effects ($\alpha_s^2$ corrections), full 3+1D evolution including transverse dynamics, and applications to heavy-ion and electron–ion collisions for a finer probe of quantum entanglement. In lattice systems, the ability to write closed-form solutions for sequential dipole transitions provides a rigorous framework for analyzing lattice trapping, cascade nucleation, and the effect of spinodal intermediate regimes, with implications for the broader study of nonlinear waves and defect dynamics (Sharma, 2020).