REVIEW
Small-x TMD distributions: JIMWLK evolution and the Gaussian truncation
High Energy Physics - Phenomenology · Submitted 31 Aug 2026 · 71 references
Pith's verdict: First systematic JIMWLK-vs-GT evolution benchmark for small-x TMDs; the representation-dependent slowdown is believable, though the finite-step extraction deserves a convergence check.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
Gaussian truncation fails to reproduce full JIMWLK evolution for most small-x TMD distributions
The reading
What carries the argument
The central object is the two-point correlation function γ^t_r extracted from the dipole evolution, which under the Gaussian truncation determines all higher correlators via Wick contractions. The paper also uses the representation-theoretic decomposition of the relevant operators into irreducible components Ω^ω_ag (based on SU(3) birdtracks), which allows systematic comparison of the truncation's performance across different color channels. The full JIMWLK evolution is implemented numerically via the Langevin formulation on a 6144² lattice, with the initial condition given by the MV model.
Extended reading notes
Core claim
The Gaussian truncation of the Balitsky hierarchy, when applied to the JIMWLK evolution of small-x TMD distributions, fails to reproduce the full numerical evolution for all operators except those directly related to the Laplacian of the dipole ($Ω^{3}$_qg, $Ω^{{8a}}$_gg, $Ω^{{8s}}$_gg). For all other operators in the quark-gluon and gluon-gluon sectors, the truncation predicts a slower evolution than the full JIMWLK equation, with the mismatch growing for higher-dimensional representations. Moreover, for the TMD distributions that are linear combinations of the building-block operators, the failure can become qualitative: the truncation predicts extra nodes or misses nodes that appear in the full evolution. The ordering of the discrepancy by representation dimension suggests that the truncation misses multi-reggeon exchanges that accelerate the evolution in larger irreps.
Load-bearing premise
The extraction of the two-point function γ^(evo)_r from the dipole using Eq. (3.13) assumes an infinitesimal rapidity step, but the actual step is finite (α_s δY = 0.05), and the fit ansatz (3.14), b_n $r^{2}$ + c_n $r^{{δ_n}}$, is purely phenomenological with no theoretical justification for its functional form.
Editorial extensions
If this is right
- Phenomenological studies that use the Gaussian truncation to compute small-x TMD distributions (e.g., for di-jet or heavy-boson production) likely underestimate the evolution speed, leading to incorrect predictions for momentum spectra and rapidity dependencies.
- An effective rescaling of rapidity by an irrep-dependent factor (e.g., Y → Y_ω) may partially correct the truncation's speed, but a deeper theoretical understanding is needed before such a patch can be trusted.
- The large-Nc limit (BK equation) cannot fully replace the JIMWLK evolution for TMD distributions, since the truncation's failure is linked to subleading-Nc effects that grow with representation dimension.
- Future approximations must incorporate at least a three-point function or a more sophisticated treatment of multi-gluon correlations to capture the correct evolution of all TMD operators.
Reading between the lines
- The systematic under-evolution in the Gaussian truncation likely arises from neglected multi-Reggeon exchanges (higher-order correlators) that become more important in higher color representations, consistent with the observed violation of Casimir scaling in the adjoint dipole evolution.
- The fact that the singlet operator Ω^1_gg in the gluon-gluon sector also shows slower evolution, despite having the smallest Casimir, suggests that the failure is not solely governed by color factors but also by the operator's structure (derivatives of Wilson lines mixing with the dipole evolution).
- A natural test of the inference that higher-point correlations accelerate the evolution would be to compare the GT and JIMWLK for the hexapole (six-point) operator, which should show even larger discrepancies than the quadrupole.
- The method of extracting γ^t_r solely from the dipole is insufficient; a simultaneous fit including other simple operators (like Ω^15_qg) would improve the truncation's accuracy but would then lose the Gaussian truncation's main advantage.
Original abstract
We extend the systematic study of small-$x$ Transverse Momentum Dependent (TMD) distributions, initiated in [Cougoulic:2026drr], by including the energy evolution. We use the JIMWLK evolution equation to study the energy dependence of three distributions in the quark-gluon sector, $\mathcal{F}_{qg}$, and seven distributions in the gluon-gluon sector, $\mathcal{F}_{gg}$. The initial conditions are provided by the McLerran-Venugopalan model. We study these TMD distributions in position space and with $N_c = 3$ fixed. We then compare these results to the Gaussian truncation of the Balitsky hierarchy. Although the Gaussian truncation worked very well for the TMD distributions at the initial condition, now we identify some cases where the Gaussian truncation provides a reasonably good description, but in the majority of cases it fails to fully reproduce numerical data. Using the fixed representation building blocks, $\Omega_{ag}^\omega$, we are able to partially explain the outcomes. This work constitutes the next step with the goal of understanding the interplays between the large-$N_c$ limit and Gaussian truncation as well as the BK/JIMWLK evolution equations for phenomenologically relevant observables.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We extend the systematic study of small-$x$ Transverse Momentum Dependent (TMD) distributions, initiated in [Cougoulic:2026drr], by including the energy evolution. We use the JIMWLK evolution equation to study the energy dependence of three distributions in the quark-gluon sector, $\mathcal{F}_{qg}$, and seven distributions in the gluon-gluon sector, $\mathcal{F}_{gg}$. The initial conditions are provided by the McLerran-Venugopalan model. We study these TMD distributions in position space and with $N_c = 3$ fixed. We then compare these results to the Gaussian truncation of the Balitsky hierarchy. Although the Gaussian truncation worked very well for the TMD distributions at the initial condition, now we identify some cases where the Gaussian truncation provides a reasonably good description, but in the majority of cases it fails to fully reproduce numerical data. Using the fixed representation building blocks, $\Omega_{ag}^\omega$, we are able to partially explain the outcomes. This work constitutes the next step with the goal of understanding the interplays between the large-$N_c$ limit and Gaussian truncation as well as the BK/JIMWLK evolution equations for phenomenologically relevant observables.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We extend the systematic study of small-$x$ Transverse Momentum Dependent (TMD) distributions, initiated in [Cougoulic:2026drr], by including the energy evolution. | contribution 0.87 | 0.87 | 0.04 | asserted 0.99 |
| s1 | We use the JIMWLK evolution equation to study the energy dependence of three distributions in the quark-gluon sector, $\mathcal{F}_{qg}$, and seven distributions in the gluon-gluon sector, $\mathcal{F}_{gg}$. | method 0.56 | 0.93 | 0.04 | asserted 0.99 |
| s2 | The initial conditions are provided by the McLerran-Venugopalan model. | method 1.00 | 0.22 | 0.02 | asserted 1.00 |
| s3 | We study these TMD distributions in position space and with $N_c = 3$ fixed. | method 0.97 | 0.91 | 0.02 | asserted 0.98 |
| s4 | We then compare these results to the Gaussian truncation of the Balitsky hierarchy. | method 0.67 | 0.89 | 0.05 | asserted 0.73 |
| s5 | Although the Gaussian truncation worked very well for the TMD distributions at the initial condition, now we identify some cases where the Gaussian truncation provides a reasonably good description, but in the majority of cases it fails to fully reproduce numerical data. | result 0.99 | 0.94 | 0.76 | asserted 0.98 |
| s6 | Using the fixed representation building blocks, $\Omega_{ag}^\omega$, we are able to partially explain the outcomes. | interpretation 0.68 | 0.92 | 0.07 | asserted 0.90 |
| s7 | This work constitutes the next step with the goal of understanding the interplays between the large-$N_c$ limit and Gaussian truncation as well as the BK/JIMWLK evolution equations for phenomenologically relevant observables. | future work 0.95 | 0.35 | 0.11 | asserted 0.93 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
References
How the literature uses it
What this paper does with the arXiv papers it cites
What it cites · 66 sentences without a confident label
Builds on 1 paper
-
In this study, we built upon previous research on small-$x$ TMD distributions by incorporating evolution into our framework in [ref] .
Small-x TMD distributions initial condition: Nc-dependence and Gaussian approximations · Discussion and prospects · 0.99
Uses the method or tool of 8 papers
-
The initial condition is generated according to the MV-model, and we refer to [ref] for further details.
Small-x TMD distributions initial condition: Nc-dependence and Gaussian approximations · Implementation of the JIMWLK evolution · 0.87
-
Gluon - Gluon: according to [ref] , we write TMD operators as:
TMD gluon distributions for multiparton processes · Definition of the operators · 0.93
-
Gluon - Gluon: according to [ref] , we write TMD operators as:
Transverse-momentum-dependent gluon distributions from JIMWLK evolution · Definition of the operators · 0.93
Compares with 1 paper
-
This situation is distinct from the cases explored in [ref] , where the dependence on the rapidity parametrization $y$ disappeared in favor of $Γ_$ being the potential integrated over the longitudinal profile.
When JIMWLK evolution really matters: the example of incoherent diffraction · Analytic results for the Gaussian truncation · 0.99
Cites as support 2 papers
-
While this program has obvious advantages, it has been shown, e.g. [ref] , that the GT does not automatically capture all relevant features of correlation functions for the quadrupole and the sextupole.
Forward dijets in high-energy collisions: evolution of QCD n-point functions beyond the dipole approximation · Introduction · 0.97
-
While this program has obvious advantages, it has been shown, e.g. [ref] , that the GT does not automatically capture all relevant features of correlation functions for the quadrupole and the sextupole.
When JIMWLK evolution really matters: the example of incoherent diffraction · Introduction · 0.97
Cites as background 26 papers
-
This approach has been widely used over the years, see e.g. [ref] and references therein.
Evolution at small x_bj: The Color Glass Condensate · Introduction · 0.99
-
This approach has been widely used over the years, see e.g. [ref] and references therein.
Forward dijets in high-energy collisions: evolution of QCD n-point functions beyond the dipole approximation · Introduction · 0.99
-
This approach has been widely used over the years, see e.g. [ref] and references therein.
JIMWLK evolution in the Gaussian approximation · Introduction · 0.99
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.