{"id":"585b69fe-a640-4855-989c-f3ef707dd792","arxiv_id":"2508.20432","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"JIMWLK evolution is added to IP-Glasma initial conditions, and the model is shown to capture the RHIC-to-LHC energy dependence of bulk observables.","lead":"This proceedings paper incorporates the perturbative JIMWLK energy evolution into the IP-Glasma initial state model for heavy ion collisions, replacing a parametrized energy dependence. It reports that the resulting setup describes RHIC and LHC multiplicity and mean transverse momentum data, and offers new predictions for oxygen and neon collisions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Initial condition at x0=0.01 is the unquantified load-bearing premise; without propagating its uncertainty, the JIMWLK-vs-Qs(x) comparison is not a clean test.","rationale":"The reader's weakest assumption correctly identifies the non-perturbative input at x0=0.01 as the load-bearing premise. My stress test concurs and sharpens it: the manuscript explicitly mentions a Bayesian analysis with uncertainty estimates (Ref. [12]) but does not use it, and it reports no variation of x0 or the IR/alpha_s regulators. Because the entire evidence for the central claim is the comparison between JIMWLK-evolved and Qs(x)-parametrized initial states, an unquantified dependence of that comparison on the initial-condition choices would undermine the claim. The proposed sensitivity scan, using the existing posterior from Ref. [12] and varying x0, would directly test whether the conclusion is robust. Since the reader already assigned CONDITIONAL on essentially these grounds, the verdict should stand unchanged; the stress test reinforces rather than alters the assessment.","tokens_in":4216,"tokens_out":12250,"duration_ms":143779,"concrete_test":"Run a sensitivity scan for Pb-Pb 5.02 TeV and Au-Au 200 GeV: (i) choose x0 = 0.005, 0.01, 0.02, re-deriving the IP-Sat parameters for each x0 from the exclusive-vector-meson data used in Ref. [11] (or from the posterior of [12]), and (ii) draw several initial-condition parameter sets from the Bayesian posterior of Ref. [12]. For each case, solve JIMWLK to the same final x, run the same MUSIC+UrQMD chain, and recompute dNch/deta and <pT>. If the spread across these predictions is comparable to or larger than the difference between the JIMWLK and Qs(x) curves, or is larger than the experimental uncertainties, the preference for JIMWLK is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central comparison in Figs. 1–2 attributes the RHIC-to-LHC energy dependence to JIMWLK, but this conclusion rests entirely on the non-perturbative initial condition at x0=0.01. Section 2 lists the free parameters (overall saturation scale, nucleon size, three-hot-spot substructure, saturation fluctuations, running-coupling scale, IR regulators) said to be constrained by the exclusive-vector-meson fit in Ref. [11]. The paper explicitly notes that a later Bayesian analysis [12] provides uncertainty estimates, but these are not propagated, and no scan over x0 or the IR/running-coupling choices is shown. If a different equally valid fit (e.g., from the posterior of [12]) or a different x0 changes the evolved nuclear geometry enough, the difference between the JIMWLK and Qs(x) setups could shrink or reverse. Then the claim that the energy dependence is 'well captured' would not be a robust prediction of perturbative evolution but an artifact of a particular choice of initial conditions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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].","tokens_in":4532,"tokens_out":3757,"duration_ms":46140,"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":[{"comment":"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.","section":"§3, Figs. 1–2"},{"comment":"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.","section":"§2, initial condition at x0=0.01"},{"comment":"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.","section":"§3, Fig. 3"}],"minor_comments":[{"comment":"The y-axis label is garbled ('10η/dchdN' should read dN_ch/dη, with appropriate scaling). Please also spell out the centrality definition used.","section":"Figure 1 caption"},{"comment":"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.","section":"§3 and Fig. 2 caption"},{"comment":"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.","section":"§1"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings summary of a result already published in Ref. [5]. The main issue is not novelty but the strength of the evidence presented: the comparison is qualitative, and the initial-condition uncertainty is not propagated. The authors can likely fix these issues by adding a quantitative comparison and a sensitivity study; if they prefer to keep the proceedings format short, the central claims should be softened or explicitly refer to the quantitative analysis in Ref. [5]."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nThe 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.\n\nFor 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.","headline":"Faithful proceedings summary of a solid PRL; nothing new here, but the JIMWLK-in-IP-Glasma result is worth knowing about.","tokens_in":5011,"tokens_out":2370,"would_cite":false,"duration_ms":27672,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Embedding JIMWLK evolution in IP-Glasma replaces a fitted saturation scale and reproduces the measured RHIC-to-LHC energy dependence of heavy-ion multiplicities.","keywords":["IP-Glasma","JIMWLK evolution","Color Glass Condensate","heavy-ion collisions","multiplicity distribution","mean transverse momentum","initial state","nuclear geometry"],"falsifier":"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.","tokens_in":4118,"feed_emoji":"⚛️","tokens_out":7456,"duration_ms":79456,"temperature":0.7,"pith_summary":"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.","feed_headline":"JIMWLK evolution matches RHIC-to-LHC multiplicity data","feed_subtitle":"Putting perturbative QCD evolution into IP-Glasma replaces a fitted saturation scale and is favored by mean-pT data.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduces the IP-Glasma initial state model whose Wilson lines are evolved here.","marker":"[2]"},{"why":"Defines the baseline IP-Glasma plus hydrodynamics setup whose energy dependence used the parametrized Qs(x); JIMWLK is compared against it.","marker":"[3]"},{"why":"Companion publication containing the full calculations this proceeding reports.","marker":"[5]"},{"why":"Derives the JIMWLK equation used for the high-energy evolution.","marker":"[6]"},{"why":"Supplies the Pb-Pb multiplicity data at 5.02 TeV used for comparison.","marker":"[7]"},{"why":"Supplies the Au-Au multiplicity data at 200 GeV used for the RHIC comparison.","marker":"[8]"},{"why":"Supplies the mean-transverse-momentum data that prefers the JIMWLK setup.","marker":"[9]"},{"why":"Provides the three-hot-spot proton substructure that seeds the initial Wilson-line configurations.","marker":"[10]"},{"why":"Constrains the non-perturbative parameters at x0=0.01 through exclusive vector meson production data.","marker":"[11]"},{"why":"Hydrodynamic evolution code used after the Yang-Mills stage to produce final-state observables.","marker":"[13]"}],"fun_headline_variants":["JIMWLK evolution predicts RHIC-to-LHC multiplicity and mean-pT","JIMWLK evolution replaces fitted Qs in IP-Glasma and matches data","Event-by-event JIMWLK evolution describes RHIC-to-LHC energy dependence","JIMWLK evolution yields RHIC-to-LHC multiplicity without fitted Qs","JIMWLK evolution replaces fitted scale to match RHIC-to-LHC mean-pT"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["JIMWLK evolution predicts RHIC-to-LHC multiplicity and mean-pT","JIMWLK evolution replaces fitted Qs in IP-Glasma and matches data","Event-by-event JIMWLK evolution describes RHIC-to-LHC energy dependence","JIMWLK evolution yields RHIC-to-LHC multiplicity without fitted Qs","JIMWLK evolution replaces fitted scale to match RHIC-to-LHC mean-pT"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001087,"raw_usage":{"total_tokens":4319,"prompt_tokens":623,"completion_tokens":3696,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":367,"completion_tokens_details":{"reasoning_tokens":3587}},"tokens_in":367,"tokens_out":3696,"duration_ms":27748,"temperature":1.0,"reasoning_tokens":3587,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:05:02.225265+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Systematic Studies of the Centrality and sqrt(s_NN) Dependence of dE_T/deta and dN_ch/deta in Heavy Ion Collisions at Mid-rapidity","cited_arxiv_id":"nucl-ex/0409015","evidence_quote":"Supplies the Au-Au multiplicity data at 200 GeV used for the RHIC comparison."}],"review_version":1}