Perturbative high-energy evolution in the IP-Glasma initial state

Bj\"orn Schenke; Chun Shen; Heikki M\"antysaari; Wenbin Zhao

REVIEW 3 major objections 3 minor 15 references

Perturbative high-energy evolution in the IP-Glasma initial state

T0 review · 3 major / 3 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Embedding JIMWLK evolution in IP-Glasma replaces a fitted saturation scale and reproduces the measured RHIC-to-LHC energy dependence of heavy-ion multiplicities.

desk verdict Faithful proceedings summary of a solid PRL; nothing new here, but the JIMWLK-in-IP-Glasma result is worth knowing about. read the letter →

arxiv 2508.20432 v1 pith:2TPSBQEV submitted 2025-08-28 nucl-th

Heikki Mäntysaari , Björn Schenke , Chun Shen , Wenbin Zhao This is my paper

classification nucl-th

keywords IP-GlasmaJIMWLKevolutionColorGlassCondensateheavy-ioncollisionsmultiplicitydistributionmeantransversemomentuminitialstatenucleargeometry

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper argues that the energy dependence of the initial state in heavy-ion collisions can be computed from perturbative QCD rather than fitted by hand. The authors place the JIMWLK equation inside the IP-Glasma model, evolving each nucleus's Wilson lines event-by-event from an initial Bjorken x0=0.01 down to the momentum scale of the collision. The resulting geometry changes with collision energy without tuning a saturation scale to each beam energy. The paper shows the RHIC-to-LHC trend of charged-particle multiplicity is captured by this evolution, and that average transverse momentum data choose JIMWLK over the previous energy-dependent saturation-scale setup. The consequence is that even simple bulk observables are sensitive to how the initial state evolves, so extracting quark-gluon plasma properties requires this theoretical evolution.

What carries the argument

The JIMWLK equation—the perturbative QCD evolution equation for the Wilson lines that represent the color field of a fast-moving nucleus—carries the argument. In this work it is solved numerically for each event, evolving the IP-Glasma initial Wilson-line configurations from x0=0.01 to x = ⟨pT⟩/√sNN, so that the nuclear geometry at a given collision energy is generated rather than parametrized. The non-perturbative starting configurations are fixed from exclusive vector meson data; JIMWLK then supplies all energy dependence.

What would settle it

Measure the centrality dependence of charged-particle multiplicity, or of mean transverse momentum, in O+O collisions at RHIC (200 GeV) and at the LHC (6.8 TeV) with enough precision to distinguish the JIMWLK-evolved curves from the Qs(x) curves in Fig. 1. Because O+O is small, JIMWLK predicts a larger and centrality-dependent difference; data that follow the Qs(x) prediction would falsify the JIMWLK-evolved geometry.

Watch

Extended reading notes

Core claim

On its own terms, the paper's central claim is that the perturbative JIMWLK equation, solved numerically event-by-event on the IP-Glasma Wilson lines, correctly produces the collision-energy dependence of the nuclear initial state. Starting from non-perturbative input at x0=0.01, JIMWLK predicts how the nucleus looks at smaller x; the authors feed those evolved Wilson lines into Yang-Mills and hydrodynamic evolution. They find that the charged-hadron multiplicity distribution in Au+Au, Pb+Pb, and O+O is described across RHIC and LHC energies, and that mean transverse momentum in Pb+Pb and p+Pb is better described by the JIMWLK-evolved setup than by the old Qs(x) parametrization. They also ad

Load-bearing premise

The calculation assumes that the non-perturbative Wilson-line configurations at x0=0.01, from IP-Sat parameters fitted to exclusive vector meson data, are the right starting point; if that input is wrong, the subsequent JIMWLK evolution is not a clean test.

Editorial extensions

If this is right

The RHIC-to-LHC energy dependence of the charged-particle multiplicity distribution is explained by JIMWLK evolution without tuning a saturation scale to each energy.
Average transverse momentum is lower with JIMWLK-evolved, smoother nuclei, and the measured data point to the JIMWLK setup over the Qs(x) setup.
For O+O collisions, which will be measured at both RHIC and LHC, the JIMWLK and Qs(x) setups differ more strongly, so O+O data can discriminate between them.
Predictions for v2 and v3 flow in O+O and Ne+Ne at 5.36 TeV are provided, with ratios expected to be robust against pre-equilibrium uncertainties.
The publicly available code allows future studies to predict energy-dependent initial states rather than parametrize them.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the claim holds, nuclear geometry at LHC energies is smoother than at RHIC; this should show up as an energy-dependent suppression of geometry-driven fluctuations in higher-order flow and event-shape observables, which can be searched for in existing data.
The calculation keeps the x0=0.01 input fixed. Propagating uncertainty in that input, through fits to exclusive vector meson data, would turn the observed agreement into a quantitative test of JIMWLK with error bars.
The same event-by-event evolution could be applied to p+A and light-ion collisions, where the geometry change is larger, making the x-dependence of subnucleonic fluctuations directly measurable.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Read this if you want a one-page advertisement for the authors' PRL (Ref [5]). It says so on page 1: results are based on calculations shown there. All figures are from that PRL. So as a standalone preprint it contains no new derivation, no new numbers, and no new physics. But what it summarizes is genuinely interesting: the authors replace the fitted Qs(x) energy dependence in IP-Glasma with event-by-event perturbative JIMWLK evolution, and show that this captures the RHIC-to-LHC energy dependence of the multiplicity distribution and mean pT. They also give predictions for O-O and Ne-Ne collisions, which is forward-looking. The code is public. That is real work, and the comparison is not circular—JIMWLK is an external evolution equation, so the energy dependence is not fitted to the bulk observables it predicts.

The soft spots are mostly about presentation and robustness. The manuscript gives no equations or numerical details; everything is deferred to Ref [5]. The model curves in Figs 1 and 2 have no uncertainty bands, so the claim that ALICE data "prefers" the JIMWLK setup is based on visual inspection. More substantively, the whole comparison rests on the non-perturbative initial condition at x0=0.01, built from IP-Sat parameters fitted to HERA/LHC data (Ref [11]). The paper notes that a later Bayesian analysis (Ref [12]) provides uncertainty estimates, but those are not propagated here, and there is no scan over x0 or the IR/running-coupling choices. So the JIMWLK-vs-Qs(x) difference could shrink or even reverse under a different equally valid initial-condition fit. This is a real limitation, but it is a limitation of the underlying PRL as much as of this proceedings summary—and the authors are transparent about it.

For a proceedings, the appropriate bar is honesty and fidelity to the cited work, which this meets. If it were submitted as a research article, it should be desk rejected because it is not a new result. Do not cite this; cite Ref [5]. But the science itself is sound and worth engaging with, so if you work on initial-state modeling or QGP property extraction, read the PRL and keep an eye on the Bayesian follow-up for uncertainty quantification.

Referee Report

3 major / 3 minor

Summary. This proceedings paper reports an extension of the IP-Glasma initial-state model in which the Wilson lines describing the colliding nuclei are evolved in Bjorken x by solving the perturbative JIMWLK equation event-by-event, rather than by inserting an energy-dependent saturation scale Q_s(x) as in the standard IP-Sat/IP-Glasma approach. The evolved initial states are coupled to MUSIC hydrodynamics and UrQMD hadronic transport, and the resulting charged-hadron multiplicities and mean transverse momenta are compared with RHIC and LHC data for Pb+Pb, Au+Au, and p+Pb collisions. Predictions are also given for O+O and Ne+Ne elliptic and triangular flow at 5.36 TeV. The central claims are that the JIMWLK-evolved setup captures the RHIC-to-LHC energy dependence of the multiplicity distribution (Sec. 3) and that even simple bulk observables are sensitive to the high-energy evolution, so that a theoretically motivated evolution is needed for precision extractions of QGP properties (Sec. 4). All numerical details, equations, and derivations are deferred to Ref. [5].

Significance. If the results hold, this is a conceptually important step: it replaces a parametrized, fitted energy dependence in the IP-Glasma initial state with a perturbative nonlinear evolution, and it provides publicly available code for the community. The O+O and Ne+Ne flow ratios in Fig. 3 are falsifiable predictions. However, the present manuscript is a proceedings summary: no equations or numerical details are given, and the main evidence in Figs. 1-2 is qualitative, with model curves lacking uncertainty bands and no quantitative goodness-of-fit measure. The initial condition at x0=0.01, constrained by exclusive vector meson data, is the unquantified load-bearing input. These issues limit the strength of the claims as stated, though they are addressable.

major comments (3)

[§3, Figs. 1–2] The central claims that the energy dependence is 'well captured' by JIMWLK and that ALICE data 'prefers' the JIMWLK setup are based on visual inspection. The model curves in Figs. 1 and 2 have no uncertainty bands, and no chi-square or other quantitative comparison is reported. In several centrality bins the difference between the JIMWLK and Q_s(x) curves is comparable to the scatter of the data, so a quantitative test is needed. Please add uncertainty estimates, at least from the parameter posterior of Ref. [12], and report a goodness-of-fit measure for both setups.
[§2, initial condition at x0=0.01] The non-perturbative Wilson-line input at x0=0.01 is the load-bearing premise. The listed parameters (overall saturation scale, nucleon size, hot-spot substructure, saturation fluctuations, running-coupling scale, IR regulators) are constrained by the exclusive vector meson fit of Ref. [11], but their uncertainties are not propagated. The paper explicitly mentions that Ref. [12] provides uncertainty estimates yet does not use them. Without a sensitivity study varying x0, IR cutoff, running-coupling scheme, or sampling the posterior of [12], the difference between JIMWLK and Q_s(x) in Figs. 1-2 could be an artifact of the specific initial condition. Please include such a scan or discuss its impact on the conclusions.
[§3, Fig. 3] The uncertainty bands on the O+O and Ne+Ne flow predictions cover only the pre-equilibrium initialization ambiguity. The paper's own conclusion is that bulk observables are sensitive to the high-energy evolution, so the flow predictions should also reflect the uncertainty in the JIMWLK initial condition and evolution parameters. Without this, the quoted cancellation of systematic uncertainties in the Ne/O ratio is incomplete.

minor comments (3)

[Figure 1 caption] The y-axis label is garbled ('10η/dchdN' should read dN_ch/dη, with appropriate scaling). Please also spell out the centrality definition used.
[§3 and Fig. 2 caption] Fig. 2 is described as p+Pb and Pb+Pb data compared with ALICE data [9], but Ref. [9] is a Xe+Xe measurement. Please check the correct reference for the p+Pb and Pb+Pb <pT> data.
[§1] Typo: 'succesfully' should be 'successfully'. Also, since all equations are deferred to Ref. [5], it would help readers to indicate equation numbers from Ref. [5] at the points where the JIMWLK evolution and the initial conditions are described.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular derivation: the energy evolution is generated by the external JIMWLK equation from independently fitted initial conditions, and the bulk observables are forward predictions.

full rationale

The paper's derivation chain is a forward calculation: Wilson lines at x0=0.01 are built from the IP-Sat/IP-Glasma model, with parameters constrained in Ref. [11] by HERA/LHC exclusive vector meson data (an external, independently published data set), and then evolved to smaller x by numerically solving the JIMWLK equation of Mueller [6], which is an external perturbative QCD evolution equation. The RHIC-to-LHC dependence of dNch/deta and <pT> in Figs. 1-2 is not defined in terms of those observables and is not fitted to them; the comparison is made against ALICE and PHENIX data. The Qs(x) comparison setup is a baseline, not the input to the JIMWLK prediction. The self-citations to Refs. [10, 11, 12] are to data-constrained determinations of the non-perturbative input, so they are real evidence rather than circular support. The paper does honestly flag a limitation: Ref. [12] provides uncertainty estimates but they are not propagated, and no scan over x0 or running-coupling/IR choices is shown. That is an uncertainty shortcoming, not a circularity: it does not make any equation equivalent to its inputs by construction. No step in the manuscript exhibits the reduction pattern (Eq. X = Eq. Y by definition, or a fitted parameter renamed as a prediction).

Assumptions & free parameters 9 free parameters · 6 assumptions · 0 invented entities

The central claim depends on the IP-Sat initial conditions at x0=0.01 with parameters fitted elsewhere, plus the standard hydro/transport setup; none of these are derived in this paper. The JIMWLK evolution itself is an external equation, but its starting point and the full simulation chain are model assumptions.

free parameters (9)

Overall saturation scale of the nucleon = not given in this paper; fitted in Ref [11] to HERA/LHC vector meson data
Sets the normalization of the Wilson line configurations and directly controls the multiplicities.
Nucleon size (impact-parameter width) = not given in this paper; fitted in Ref [11]
Controls the collision geometry and the centrality dependence of observables.
Hot-spot substructure (three hot spots) = from Ref [10]
Determines event-by-event fluctuations that seed flow harmonics.
Saturation scale fluctuation parameter = not given in this paper; fitted in Ref [11]
Controls the magnitude of geometric fluctuations in the Wilson lines.
Coordinate-space running-coupling scale = not given in this paper
Controls the speed of the JIMWLK evolution; different choices change the x-dependence of the result.
IR regulator parameters (Coulomb tail cutoff) = not given in this paper
Suppress unphysical long-range Coulomb tails and affect the evolved density profiles.
Initial Bjorken x0 = 0.01
Starting point for the JIMWLK evolution; chosen because the IP-Sat model is applicable there.
Hydro initial time tau0 = 0.4 fm/c
Proper time at which the energy-momentum tensor is passed from IP-Glasma to MUSIC.
Switching energy esw = 0.18 GeV/fm3
Energy density at which fluid cells are converted into hadrons for UrQMD.

assumptions (6)

domain assumption JIMWLK equation is the correct perturbative QCD evolution for small-x gluon fields
The paper invokes JIMWLK [6] as the high-energy evolution without deriving or testing its validity in this setup.
standard math Eikonal approximation and Wilson lines describe high-energy scattering
The IP-Glasma degrees of freedom are Wilson lines V(x_perp), based on eikonal propagation of quarks through the nuclear color field.
domain assumption IP-Sat dipole model with parameters from Ref [11] describes the nucleus at x0=0.01
The initial non-perturbative input is taken from IP-Sat with parameters fitted to exclusive vector meson data; the paper does not question this input.
domain assumption Classical Yang-Mills evolution from tau=0 to tau0=0.4 fm/c captures the pre-equilibrium stage
The 'glasma' stage is simulated by solving Yang-Mills equations in IP-Glasma for the early time evolution.
domain assumption MUSIC hydrodynamics with the specified equation of state and transport coefficients correctly converts T_mu_nu to hadronic observables
The hydrodynamic stage is assumed to be valid for all centralities and energies studied.
domain assumption UrQMD hadronic cascade correctly models the late-stage hadronic rescattering
The final observable multiplicities and spectra are assumed to be accurately produced by UrQMD after hadronization.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Perturbative high-energy evolution in the IP-Glasma initial state." pith.science (2026). https://pith.science/paper/2TPSBQEV

@misc{pith2026250820432,
  author       = {Pith},
  title        = {Pith review of: Perturbative high-energy evolution in the IP-Glasma initial state},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2TPSBQEV}},
  note         = {Machine review of arXiv:2508.20432}
}

abstract

We include the perturbative JIMWLK energy evolution into the IP-Glasma initial state description used to simulate the early-time dynamics in heavy ion collisions. By numerically solving the JIMWLK equation on an event-by-event basis, we obtain the energy (Bjorken-$x$) dependent structure of the colliding nuclei. Combining the initial state with hydrodynamic simulations, this enables us to predict how observables evolve when moving from RHIC to LHC energies.

Figures

Figures reproduced from arXiv: 2508.20432 by the authors.

**Figure 2.** Average transverse momentum in proton-lead and lead-lead collisions at midrapidity compared to ALICE data [9]. Figure from Ref. [5]. lines V(x⊥), that describe the eikonal propagation of a quark through the color field of the nucleus at a given transverse coordinate x⊥. These Wilson lines are then evolved to smaller x by solving the JIMWLK equations. The free parameters of the model control the overall saturation sc… view at source ↗

**Figure 3.** Elliptic and triangular flow in O+O and Ne+Ne collisions at 5.36 TeV. The uncertainty bands include systematic variation of the pre-equilibrium condition for initializing hydrodynamics, namely matching full T µν, setting initial bulk viscous pressure Π = 0, and depositing only energy density at an early proper time. Those variations largely cancel between the two systems in the vn{2} ratios. will be measured at bo… view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

15 extracted references · 2 canonical work pages

[5]

Mäntysaari, B

H. Mäntysaari, B. Schenke, C. Shen, W. Zhao, Collision-Energy Dependence in Heavy- Ion Collisions from Nonlinear QCD Evolution, Phys. Rev. Lett. 135, 022302 (2025), 2502.05138. 10.1103/gf4y-p5j7

arXiv 2025
[11]

Mäntysaari, F

H. Mäntysaari, F. Salazar, B. Schenke, Nuclear geometry at high energy from ex- clusive vector meson production, Phys. Rev. D 106, 074019 (2022), 2207.03712. 10.1103/PhysRevD.106.074019

arXiv 2022
[12]

Mäntysaari, H

H. Mäntysaari, H. Roch, F. Salazar, B. Schenke, C. Shen, W. Zhao, Global Bayesian Analysis of J/Ψ Photoproduction on Proton and Lead Targets (2025), 2507.14087

arXiv 2025
[1]

Heinz, B

U. Heinz, B. Schenke, Hydrodynamic Description of the Quark-Gluon Plasma (2024), 2412.19393

arXiv 2024
[2]

Schenke, P

B. Schenke, P. Tribedy, R. Venugopalan, Fluctuating Glasma initial conditions and flow in heavy ion collisions, Phys. Rev. Lett. 108, 252301 (2012), 1202.6646. 10.1103/PhysRevLett.108.252301

arXiv 2012
[3]

Schenke, C

B. Schenke, C. Shen, P. Tribedy, Running the gamut of high energy nuclear collisions, Phys. Rev. C 102, 044905 (2020), 2005.14682. 10.1103/PhysRevC.102.044905

arXiv 2020
[4]

Kowalski, D

H. Kowalski, D. Teaney, An Impact parameter dipole saturation model, Phys. Rev. D 68, 114005 (2003), hep-ph/0304189. 10.1103/PhysRevD.68.114005

arXiv 2003
[6]

Mueller, A Simple derivation of the JIMWLK equation, Phys

A.H. Mueller, A Simple derivation of the JIMWLK equation, Phys. Lett. B 523, 243 (2001), hep-ph/0110169. 10.1016/S0370-2693(01)01343-0

arXiv 2001

Show all 15 references

[7]

Adam et al

J. Adam et al. (ALICE), Centrality Dependence of the Charged-Particle Multiplicity Density at Midrapidity in Pb-Pb Collisions at √sNN = 5.02 TeV, Phys. Rev. Lett. 116, 222302 (2016), 1512.06104. 10.1103/PhysRevLett.116.222302

2016 arXiv
[8]

Adler et al

S.S. Adler et al. (PHENIX), Systematic studies of the centrality and √sNN depen- dence of the d ET/dη and d Nch/dη in heavy ion collisions at mid-rapidity, Phys. Rev. C 71, 034908 (2005), [Erratum: Phys.Rev.C 71, 049901 (2005)], nucl-ex/0409015. 10.1103/PhysRevC.71.034908

2005 arXiv
[9]

Acharya et al

S. Acharya et al. (ALICE), Transverse momentum spectra and nuclear modification factors of charged particles in Xe-Xe collisions at√sNN = 5.44 TeV, Phys. Lett. B 788, 166 (2019), 1805.04399. 10.1016/j.physletb.2018.10.052

2019 arXiv
[10]

Mäntysaari, B

H. Mäntysaari, B. Schenke, Evidence of strong proton shape fluctuations from inco- herent diffraction, Phys. Rev. Lett. 117, 052301 (2016), 1603.04349. 10.1103/Phys- RevLett.117.052301

2016 arXiv
[13]

Schenke, S

B. Schenke, S. Jeon, C. Gale, (3+1)D hydrodynamic simulation of relativistic heavy-ion collisions, Phys. Rev. C82, 014903 (2010), 1004.1408. 10.1103/PhysRevC.82.014903

2010 arXiv
[14]

Bass et al., Microscopic models for ultrarelativistic heavy ion collisions, Prog

S.A. Bass et al., Microscopic models for ultrarelativistic heavy ion collisions, Prog. Part. Nucl. Phys. 41, 255 (1998), nucl-th/9803035. 10.1016/S0146-6410(98)00058-1

1998 arXiv
[15]

JIMWLK + IP-Glasma, https: //github.com/schenke/ipglasma/tree/ipglasma_jimwlk (2025)

2025

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.

[1] [5]

Mäntysaari, B

H. Mäntysaari, B. Schenke, C. Shen, W. Zhao, Collision-Energy Dependence in Heavy- Ion Collisions from Nonlinear QCD Evolution, Phys. Rev. Lett. 135, 022302 (2025), 2502.05138. 10.1103/gf4y-p5j7

arXiv 2025

[2] [11]

Mäntysaari, F

H. Mäntysaari, F. Salazar, B. Schenke, Nuclear geometry at high energy from ex- clusive vector meson production, Phys. Rev. D 106, 074019 (2022), 2207.03712. 10.1103/PhysRevD.106.074019

arXiv 2022

[3] [12]

Mäntysaari, H

H. Mäntysaari, H. Roch, F. Salazar, B. Schenke, C. Shen, W. Zhao, Global Bayesian Analysis of J/Ψ Photoproduction on Proton and Lead Targets (2025), 2507.14087

arXiv 2025

[4] [1]

Heinz, B

U. Heinz, B. Schenke, Hydrodynamic Description of the Quark-Gluon Plasma (2024), 2412.19393

arXiv 2024

[5] [2]

Schenke, P

B. Schenke, P. Tribedy, R. Venugopalan, Fluctuating Glasma initial conditions and flow in heavy ion collisions, Phys. Rev. Lett. 108, 252301 (2012), 1202.6646. 10.1103/PhysRevLett.108.252301

arXiv 2012

[6] [3]

Schenke, C

B. Schenke, C. Shen, P. Tribedy, Running the gamut of high energy nuclear collisions, Phys. Rev. C 102, 044905 (2020), 2005.14682. 10.1103/PhysRevC.102.044905

arXiv 2020

[7] [4]

Kowalski, D

H. Kowalski, D. Teaney, An Impact parameter dipole saturation model, Phys. Rev. D 68, 114005 (2003), hep-ph/0304189. 10.1103/PhysRevD.68.114005

arXiv 2003

[8] [6]

Mueller, A Simple derivation of the JIMWLK equation, Phys

A.H. Mueller, A Simple derivation of the JIMWLK equation, Phys. Lett. B 523, 243 (2001), hep-ph/0110169. 10.1016/S0370-2693(01)01343-0

arXiv 2001