{"id":"073a4fae-3254-412e-b70b-85486524b146","arxiv_id":"1908.02866","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A hydrodynamic-attractor relation between initial-state energy and final charged-particle multiplicity explains the centrality dependence of yields and lets the authors estimate pre-equilibrium energy densities.","lead":"This paper derives a formula connecting the energy of the initial quark-gluon plasma state to the number of particles eventually measured, using hydrodynamic attractor theory. It then uses the formula to estimate the initial energy density and to constrain the plasma's shear viscosity.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The promotion of the homogeneous attractor relation Eq. (6) to a pointwise local map in Eq. (12) is the least-secure step; the cited support only checks the longitudinal pressure relation, not the pointwise entropy map itself.","rationale":"The reader's weakest-assumption analysis identified the same load-bearing concern: the promotion of the homogeneous, boost-invariant result to a transverse local map. This is indeed the pivot of the paper because both the centrality dependence and the subsequent extraction of dE0/dηs and the η/s constraint rely on Eq. (12). I agree with the reader's assessment that the paper's justification via Ref. [43] is incomplete: that reference tests the validity of the 1D longitudinal-pressure relation, not the pointwise equality of entropy production to the homogeneous attractor formula. However, this concern does not by itself invalidate the central relation Eq. (6), which is well-supported by multiple independent microscopic theories (QCD kinetic theory, RTA, AdS/CFT) and collapses onto a universal attractor. The paper's honest discussion of limitations, including the assumption of complete equilibration, the fitted normalization, and the omission of propagated systematic uncertainties, is already reflected in a CONDITIONAL verdict. My proposed test would settle whether the local promotion is quantitatively safe; if it passes, the central claims stand, and if it fails, the paper would need to incorporate transverse corrections, but that would be a modification rather than a complete rejection. Therefore the reader's CONDITIONAL verdict remains appropriate.","tokens_in":62,"tokens_out":5524,"duration_ms":263557,"concrete_test":"Run a full (3+1)-dimensional kinetic-theory simulation of the pre-equilibrium phase (e.g., the code of Kurkela et al., Phys. Rev. C99, 034910) for a Pb-Pb collision at √sNN = 2.76 TeV, using the same initial (eτ)_0 profile from Eq. (17) and a fluctuating Glauber geometry. At τ_hydro ≈ 1 fm/c, compute the transverse entropy profile (sτ)_hydro(x⊥) and the resulting ⟨dNch/dη⟩ in several centrality bins (e.g., 0–5%, 20–30%, 40–50%). Compare to the pointwise map (sτ)_hydro(x⊥) ∝ [(eτ)_0(x⊥)]^{2/3} integrated as in Eq. (12). If the ratio of the full simulation to the map deviates by more than ~10% in any centrality bin, the local promotion in Eq. (12) is not justified and the centrality shape prediction needs to be revised to include transverse corrections.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (6) is rigorously derived for a boost-invariant, transversely homogeneous conformal system. The step that makes the paper's phenomenological claims possible is the promotion of Eq. (6) to a pointwise local map at each transverse coordinate, leading to Eq. (12). This requires that the entropy per rapidity produced at each x⊥ depends only on the local initial (eτ)_0(x⊥), with no mixing from transverse gradients, pre-flow, or nonlocal attractor dynamics. The paper justifies this by stating that transverse gradients are negligible over τ_hydro and cites Ref. [43]. However, Ref. [43] validates the one-dimensional constitutive relation for PL/e in the presence of transverse gradients, not the stronger claim that the total entropy production at each transverse point equals the homogeneous attractor result evaluated at the local energy density. The linearized fluctuation relation Eq. (13) covers only long-wavelength perturbations, whereas Eq. (12) is a fully nonlinear local map. If transverse gradients cause energy flux between fluid cells before hydrodynamization, the exponent 2/3 and the centrality dependence could be modified. This effect is expected to be largest in peripheral collisions, where the model without fluctuations overpredicts data; the paper attributes this to geometrical fluctuations, but breakdown of the local map is a competing explanation. Because the extraction of dE0/dηs in Eq. (21) and the η/s bound also rely on this local promotion via the effective area A⊥, the concern affects all main phenomenological conclusions, not just the shape comparison.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript derives a closed relation between the initial-state energy per unit rapidity and the hydrodynamized entropy per unit rapidity, Eq. (6), using the existence of a hydrodynamic attractor for boost-invariant, transversely homogeneous conformal systems. It then promotes this relation to a pointwise local map, Eq. (12), to compute the centrality dependence of charged-particle multiplicities from CGC/Glauber initial conditions, compares the resulting shapes to data at RHIC and LHC, and inverts the relation to estimate the initial energy density dE0/dηs and to discuss constraints on η/s. The central derivation in Eq. (6) is clean and is benchmarked against QCD kinetic theory, Yang-Mills kinetic theory, Boltzmann RTA, and AdS/CFT, with the constant C∞ varying at the 10–20% level across theories. The phenomenological applications, however, rely on two additional steps that are less secure: the local promotion of Eq. (6) and the normalization of the theoretical curves to one centrality bin per system.","tokens_in":12338,"tokens_out":14112,"duration_ms":143470,"significance":"If the central relation Eq. (6) can be justified beyond the idealized homogeneous setup, it is a valuable and falsifiable result: it turns measured multiplicities into a direct probe of the pre-equilibrium energy density and of η/s, with all prefactors explicit and with remarkable insensitivity to the microscopic equilibration mechanism. The paper also provides a transparent derivation that connects attractor theory to a standard heavy-ion observable, and the comparison in Fig. 2 is a clean shape test of the centrality dependence. The main value is in the explicit analytic relation and in the proposal that the early-time entropy production can be quantified by the exponent 2/3 rather than by model-dependent parametric estimates. The significance is reduced, but not eliminated, by the fact that the data comparison uses one fitted normalization per system and that the local-map step is supported only indirectly.","major_comments":[{"comment":"The promotion of Eq. (6) to a pointwise local map is the key step that produces the centrality prediction Eq. (12), but the cited justification is insufficient. Ref. [43] validates that the one-dimensional constitutive relation for P_L/e remains approximately valid in the presence of transverse gradients; it does not establish that the entropy per rapidity produced at each transverse point equals the homogeneous-attractor result evaluated at the local energy density. Eq. (13) is a linearized long-wavelength relation, whereas Eq. (12) is a fully nonlinear local map. In peripheral collisions, where transverse gradients are largest, the overprediction of dNch/dη could equally be a symptom of the breakdown of this local map rather than of the neglect of geometrical fluctuations. The authors should test the local map against a 3D kinetic-theory or KøMPøST evolution with transverse structure and show that the local entropy profile, not only P_L/e, follows the attractor relation.","section":"Centrality dependence of particle multiplicity, Eq. (12)"},{"comment":"The data comparison in Fig. 2 uses one fitted normalization per collision system: the combination in Eq. (20) is matched to data, with values 10.7 and 4.8 for Pb-Pb and Au-Au, and the other systems' data are scaled by arbitrary factors shown in the caption. Consequently, the agreement tests the shape of the centrality curve, not an ab initio prediction of the absolute multiplicity. The abstract's claim of an \"ab initio model of energy deposition\" is therefore stronger than what is demonstrated, because the normalization constant N_e (or the combination in Eq. (20)) is constrained by experiment. In addition, no systematic uncertainties from the Glauber parameters, the S/Nch range, ν_eff, or the normalization are propagated into the theoretical curves; the authors should either propagate them or explicitly state that the comparison is shape-only.","section":"Fig. 2 and Eq. (20) in the supplemental material"},{"comment":"The claimed upper limit η/s ≲ 0.4 is conditional on the assumption of (nearly-)complete equilibration, and the paper itself states that for η/s ≳ 0.4 the estimates need to be revised because the QGP would not equilibrate in peripheral collisions. The bound is therefore derived using the same assumption that it is supposed to constrain, and no quantitative model of the incomplete-equilibration regime is provided to close the loop. Also, the uncertainty bands in Fig. 3 vary only η/s and C∞; the substantial uncertainties in S/Nch (6.7–8.5) and ν_eff enter Eq. (21) with powers 3/2 and −1/2 and are not shown. The η/s constraint should be presented as an illustration of potential sensitivity rather than as a firm constraint.","section":"Discussion and Eq. (21), Fig. 3"}],"minor_comments":[{"comment":"Eq. (18) appears to have the Jacobian J in the numerator, whereas Eq. (7) has 1/J; since J≈1.1, this is a 10% discrepancy that is absorbed by the fitted normalization N_e, but the two expressions should be made consistent.","section":"Eq. (18), supplemental material"},{"comment":"The denominator in the (η/s) factor is printed as 2/(4π), which is inconsistent with the values η/s=0.08 and 0.16 quoted in Fig. 3; if the intended reference value is the KSS bound, it should read 1/(4π).","section":"Eq. (14)"},{"comment":"The text states that the variation of C∞ is only at the ~10% level, but the values shown in Fig. 1 span 0.87–1.06, a spread of about 20% around the central value; since C∞ enters Eq. (21) with the power −9/8, this variation induces a ~20–25% uncertainty in dE0/dηs and should be described accordingly.","section":"Fig. 1 and surrounding text"},{"comment":"The multiplicative factors attached to the Xe-Xe, U-U, Au-Au, and Cu-Cu data are not defined in the caption; it should be stated explicitly that these are overall normalization factors used to collapse the data onto a universal curve and whether the same factors are applied to the theory curves.","section":"Fig. 2 caption"},{"comment":"There is a typo in \"Relativstic Heavy-Ion Collider\" that should be corrected to \"Relativistic\".","section":"Introduction (page 1)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well-written and the central analytic result Eq. (6) is solid under its stated assumptions. The main issue is that the phenomenological claims in Figs. 2 and 3 and the η/s constraint rely on the local promotion of Eq. (6), and the current justification is weaker than the claim. The authors should at least quantify the validity of the local map with 3D kinetic-theory or KøMPøST evolutions, and they should revise the abstract and text to avoid overstating the 'ab initio' nature of the comparison given the per-system normalization. This is fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The central formula, Eq. (6), is the real thing: a closed-form attractor relation between (e tau)_0 and (s tau)_hydro, with the 2/3 power and (eta/s)^{1/3} dependence made explicit. It goes beyond the earlier parametric estimates, and the benchmarking in Fig. 1 makes the case that C_infinity is theory-insensitive at the ~10% level. That part is solid.\n\nThe centrality derivation, Eq. (12), is also genuinely new, and the comparison in Fig. 2 is interesting: after one overall normalization per system, the shape tracks the data much better than the initial-gluon estimate in Eq. (11). The paper is honest that the normalization is fitted, that fluctuations are needed in peripheral bins, and that the eta/s bound assumes complete equilibration. Those are real caveats, but they are named.\n\nThe soft spot that worries me is the step from Eq. (6) to Eq. (12). The paper promotes a boost-invariant, transversely homogeneous result to a pointwise local map, justified by saying transverse gradients are negligible over tau_hydro and citing Ref. [43]. But [43] validates the one-dimensional constitutive relation for P_L/e in the presence of transverse gradients, not the stronger claim that the entropy produced at each x_perp equals the homogeneous attractor value evaluated at the local (e tau)_0. The linearized relation Eq. (13) covers only long wavelengths; Eq. (12) is fully nonlinear. If transverse energy flux or memory of the 3D structure matters before hydrodynamization, the 2/3 exponent and the centrality curve could shift. This mainly bites in peripheral collisions, the same place the model overpredicts data; the paper blames geometrical fluctuations, but local-map breakdown is a competing explanation. That deserves a direct test—for instance, KøMPøST or kinetic theory with transverse structure—or at least a softer claim about Eq. (12).\n\nNone of this undermines the central attractor relation. It is a useful and citable result; the phenomenological reach is real but conditional. I would send this to a serious referee, asking for a sharpened discussion of the local promotion and, if possible, a check of Eq. (12) against transversely inhomogeneous kinetic theory. For a heavy-ion or attractor audience, it deserves the time.\n\nWould I cite it? Yes, for Eq. (6) and the benchmark curves.","headline":"A genuinely useful closed-form attractor relation (Eq. 6) with a clever centrality prediction; the main weak point is the local-map promotion in Eq. (12), which needs stronger backing.","tokens_in":12925,"tokens_out":2968,"would_cite":true,"duration_ms":32264,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper establishes a two-thirds-power relation between the initial energy of a relativistic nuclear collision and its final charged-particle multiplicity, and uses it to extract the energy density and shear viscosity of the…","keywords":["hydrodynamic attractor","heavy-ion collisions","quark-gluon plasma","entropy production","initial-state energy density","charged particle multiplicity","shear viscosity","pre-equilibrium dynamics"],"falsifier":"Compute the initial energy per unit rapidity $(e\\tau)_0$ from a first-principles calculation at one collision energy and compare the multiplicity predicted by Eqs. (6)-(7) with measured $dN_{\\rm ch}/d\\eta$ across centrality; if the implied $\\eta/s$ varies with centrality or energy, or if the reconstructed $dE_0/d\\eta_s$ falls below the measured final-state energy $dE_{\\rm final}/dy$ in any centrality bin, the central relation is falsified.","tokens_in":11825,"feed_emoji":"⚛️","tokens_out":12996,"duration_ms":123625,"temperature":0.7,"pith_summary":"This paper claims that the number of charged particles eventually produced in a high-energy nuclear collision is fixed, to the two-thirds power, by the energy deposited in the initial state, through the universal early-time behavior known as a hydrodynamic attractor. The central formula, Eq. (6), gives the entropy per unit rapidity after equilibration as a known constant times $(e\\tau)_0^{2/3}$, where $(e\\tau)_0$ is the initial energy per unit rapidity. If correct, this turns the accurately measured particle multiplicities into a direct measurement of the energy density of the far-from-equilibrium quark-gluon plasma before it thermalizes, and a constraint on its shear viscosity. The same relation explains the measured universal centrality dependence of multiplicities across collision systems using only nuclear geometry and pre-equilibrium entropy production.","feed_headline":"Particle yields reveal pre-collision energy via 2/3 law","feed_subtitle":"Under this attractor relation, particle counts become a direct probe of the plasma before it thermalizes.","key_machinery":"The load-bearing object is the hydrodynamic attractor for the energy density, the curve\n$$E(\\tilde w)=\\frac{e(\\tau)\\$tau^{{4/3}}$}{e_{\\rm hydro}\\tau_{\\rm hydro}^{4/3}},\\qquad \\tilde w=\\frac{T_{\\rm eff}(\\tau)\\tau}{4\\pi\\eta/s},$$\nwhich interpolates between $C_\\infty^{-1}\\tilde w^{4/9}$ at early free-streaming times and $1-2/(3\\pi\\tilde w)$ near viscous hydrodynamics. Its role is to supply the far-from-equilibrium constitutive relation $P_L/e=f(\\tilde w)$, so that Eq. (1) can be integrated without committing to one microscopic theory. The constant $C_\\infty$, which encodes how much entropy is produced between free streaming and thermalization, is nearly the same across very different microscopic descriptions, which is why the resulting 2/3-power relation carries over across theories.","core_discovery":"The paper's central discovery is a quantitative one-to-one relation between the initial-state energy and the final-state entropy in relativistic nuclear collisions. Starting from energy-momentum conservation for a boost-invariant, transversely homogeneous conformal system and using the hydrodynamic-attractor curve $E(\\tilde w)$ shared by different microscopic equilibration theories, it derives\n$$(s\\tau)_{\\rm hydro} = \\frac{4}{3} C_\\$infty^{{3/4}}$ \\left(4\\pi\\eta/s\\right)^{1/3} \\left(\\frac{\\$pi^{2}$\\nu_{\\rm eff}}{30}\\right)^{1/3} (e\\tau)$_0^{{2/3}}$.$$\nBecause the entropy at freeze-out is proportional to the measured charged-particle multiplicity, this relation lets one read the initial energy per unit rapidity $(e\\tau)_0$ backwards from $dN_{\\rm ch}/d\\eta$. For central Pb-Pb collisions at 2.76 TeV the paper quotes $e_0\\approx 270$ GeV/fm$^3$ at $\\tau_0=0.1$ fm/c, nearly three orders of magnitude above the QCD crossover energy density. Applied to centrality dependence, the same formula predicts a nontrivial dependence on the nuclear thickness through $\\int d^2 x_\\perp (T_< \\sqrt{T_>})^{2/3}$, which reproduces the measured universal curve once pre-equilibrium entropy production is included.","pith_inferences":["If the 2/3 relation holds event by event, multiplicity fluctuations should carry the imprint of the pre-equilibrium phase; comparing the variance of measured $dN_{\\rm ch}/d\\eta$ with the variance of modeled $(e\\tau)_0$ in a Glauber Monte Carlo would test the exponent without invoking final-state hydrodynamics.","Because the extracted initial energy and $\\eta/s$ enter only through $(e\\tau)_0^{2/3}$ and $(\\eta/s)^{1/3}$, an independent first-principles computation of the initial energy would convert measured multiplicities into a direct measurement of $\\eta/s$, bypassing flow modeling entirely.","Deviations from the predicted centrality curve in smaller systems, such as oxygen-oxygen or proton-nucleus collisions, would delimit where attractor-based entropy production stops being the whole story and where three-dimensional pre-equilibrium dynamics must be included.","The difference between reconstructed initial energy and measured final energy could be cross-checked in existing hydrodynamic simulations: the integrated longitudinal work should match the energy lost between the initial and final states, making the work estimate internally testable."],"forward_implications":["Measured charged-particle multiplicities become an indirect calorimeter for the energy density of the pre-equilibrium quark-gluon plasma; for central Pb-Pb at 2.76 TeV the paper quotes $e_0\\approx 270$ GeV/fm$^3$ at $\\tau_0=0.1$ fm/c.","The centrality dependence of $dN_{\\rm ch}/d\\eta$ across Au-Au, Cu-Cu, Pb-Pb, Xe-Xe and U-U collisions follows from pre-equilibrium entropy production combined with color-glass-condensate energy deposition, whereas counting initial gluons alone does not reproduce the data.","The reconstructed initial-state energy per unit rapidity exceeds the measured final-state energy by a factor of two to three in central collisions, quantifying the work done against the longitudinal expansion.","Large values of $\\eta/s$ and $C_\\infty$ are excluded because they would make the reconstructed initial energy fall below the measured final energy; assuming full equilibration, the paper obtains an upper limit $\\eta/s\\lesssim 0.4$.","Fluctuations of the initial energy are transmitted to the entropy with a 2/3 exponent, $\\delta s_{\\rm hydro}/s_{\\rm hydro}=\\frac{2}{3}\\delta e_0/e_0$, which alters the statistics of initial-state fluctuations and helps describe peripheral collisions."],"supporting_citations":[{"why":"Supplies the boost-invariant energy-momentum conservation law, Eq. (1), and the Bjorken scaling used to factor out the late-time evolution.","marker":"[12]"},{"why":"Establishes the hydrodynamic attractor concept that the paper uses to describe early-time dynamics.","marker":"[14]"},{"why":"Provides the AdS/CFT attractor curve and the late-time viscous asymptotic form used in Eq. (4) and Fig. 1.","marker":"[15]"},{"why":"Supplies the QCD kinetic theory attractor curve with $C_\\infty=0.87$ and the linearized kinetic evolution that underlies the quantitative prefactors.","marker":"[26]"},{"why":"Gives the entropy per charged particle, in the range 6.7 to 8.5, needed to convert the hydrodynamized entropy in Eq. (6) into a measured multiplicity in Eq. (7).","marker":"[32]"},{"why":"Shows that transverse gradients remain negligible during the pre-equilibrium phase, justifying the local promotion of Eq. (6) to Eq. (12).","marker":"[43]"},{"why":"Provides the relation between initial energy density and saturation scales, $(e\\tau)_0\\propto (Q_<^s)^2 Q_>^s$, used in Eq. (8).","marker":"[45]"},{"why":"Supports the dilute-dense $k_T$-factorization estimate of initial energy deposition used to connect saturation scales to nuclear thickness.","marker":"[49]"},{"why":"Supplies the optical and Monte Carlo Glauber model of nuclear thickness and centrality classes used for the comparisons in Figs. 2 and 3.","marker":"[59]"}],"fun_headline_variants":["Attractor scaling lets particle yields reveal initial energy","Particle counts become a direct probe of pre-collision energy","The 2/3 law: how multiplicity reads the plasma's starting energy","Hydrodynamic attractors expose the plasma's initial energy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the early-time plasma in a real collision is close enough to a boost-invariant, transversely homogeneous conformal system that the one-dimensional attractor relation can be applied point by point across the transverse plane; if transverse gradients or three-dimensional memory of the initial state matter before hydrodynamization, the 2/3 exponent and the centrality predictions would change.","fun_headline_variants_meta":{"raw":{"variants":["Attractor scaling lets particle yields reveal initial energy","Particle counts become a direct probe of pre-collision energy","The 2/3 law: how multiplicity reads the plasma's starting energy","Hydrodynamic attractors expose the plasma's initial energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000871,"raw_usage":{"total_tokens":3758,"prompt_tokens":919,"completion_tokens":2839,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":535,"completion_tokens_details":{"reasoning_tokens":2770}},"tokens_in":535,"tokens_out":2839,"duration_ms":22614,"temperature":1.0,"reasoning_tokens":2770,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:31:47.773777+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the initial energy per unit rapidity $(e\\tau)_0$ from a first-principles calculation at one collision energy and compare the multiplicity predicted by Eqs. (6)-(7) with measured $dN_{\\rm ch}/d\\eta$ across centrality; if the implied $\\eta/s$ varies with centrality or energy, or if the reconstructed $dE_0/d\\eta_s$ falls below the measured final-state energy $dE_{\\rm final}/dy$ in any centrality bin, the central relation is falsified.","supporting_citations":[],"review_version":1}