{"id":"48e36a3c-17f6-40c7-a734-63a17382164f","arxiv_id":"1908.09728","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The authors extract the temperature (222 MeV, about 2.6 trillion kelvin), entropy density (20 fm^-3), and speed of sound (c_s^2 = 0.24) of the quark-gluon plasma from LHC lead-lead data, consistent with lattice QCD.","lead":"This paper extracts the temperature, entropy density, and speed of sound of the quark-gluon plasma created in lead-lead collisions at the LHC, using particle data plus hydrodynamic simulations. It offers a new way to turn heavy-ion collision measurements into direct thermodynamic properties of the hottest fluid ever made in a laboratory.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The s/T^3 agreement is not an independent check: Eq. (1) defines V_eff through the input EOS, so Eq. (4) mostly returns that EOS; only c_s^2 from Eq. (7) is non-circular, and its assumed scaling laws are tested only with the soft EOS.","rationale":"The paper proposes a clever, physically motivated way to convert measured and multiplicity into an effective temperature and thermodynamic observables. The central claim requires that the calibration relation = 3.07 T_eff and the effective-volume behavior are robust properties of the hydrodynamic description, not artifacts of one EOS choice. The reader's weakest assumption points at this calibration. My stress-test sharpens it: for the entropy density, the circularity is not merely partial. Because Eq. (1) defines V_eff through the same EOS, Eq. (4) can reduce to the input s_EOS(T_eff) when the hydrodynamic entropy is normalized to the measured multiplicity with the same entropy-per-charged-particle input. The paper does not report whether its simulation's entropy per charged particle equals the adopted S/N_ch = 6.7 from Ref. [25], so the extent of tautology is undetermined. The speed-of-sound result is more independent, since Eq. (7) uses only relative variations and cancels the absolute calibration factors, but it still relies on the energy-independence of V_eff and of the /T_eff ratio, which are checked only within the soft-EOS framework. The stiff-EOS cross-check is too coarse to set the claimed error bars. A concrete, sufficient test is to rerun the extraction with an independent lattice EOS, e.g., the HotQCD parametrization, retuning initial conditions to the same data. If T_eff and s/T^3 move beyond the quoted errors, the confirmation claim must be downgraded. If they do not, the concern is resolved. I also note a secondary, easily checked issue: the text cites Ref. [13] for the 681 MeV value, but Ref. [13] is a Xe-Xe 5.44 TeV paper, while Eq. (3) applies the value to Pb-Pb at 5.02 TeV; this citation should be corrected or the data source verified, although it does not change the conceptual concern. Overall, the reader's CONDITIONAL verdict is appropriate: the method and c_s^2 result are valuable, but the entropy-density agreement and the absolute temperature calibration need an explicit EOS-sensitivity study before the 'agreement with ab-initio calculations' claim can be taken at face value.","tokens_in":10994,"tokens_out":13064,"duration_ms":137485,"concrete_test":"Repeat the central calibration using a second modern lattice EOS (e.g., HotQCD) in MUSIC, retuning initial conditions to reproduce the same dN_ch/deta and = 681 MeV, and recompute T_eff (Eq. 3), s/T^3 (Eqs. 4-5), and c_s^2 (Eq. 8). If these move outside the quoted error bars, the s/T^3 agreement and the +/- 9 MeV and +/- 0.04 uncertainties are calibration-driven; if they stay within errors, the circularity charge is largely answered.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Reader's circularity concern lands, and it is more specific than stated. In Eq. (1), V_eff is defined by S = s(T_eff)V_eff, so if the hydrodynamic S is used, V_eff = S/s_EOS(T_eff). Eq. (4) then gives s(T_eff) = (1/V_eff)(S/N_ch)(dN_ch/dy) = S/V_eff = s_EOS(T_eff). If the default run's entropy per charged particle coincides with the assumed S/N_ch = 6.7, the 'measured' s(T_eff)/T^3 in the top panel of Fig. 4 is a restatement of the s95p-v1 input EOS, not a confirmation of lattice QCD. The paper does not report this consistency check. The genuinely non-circular result is c_s^2 from Eq. (7), because relative variations cancel V_eff and the 3.07 calibration; but Eq. (7) assumes V_eff is approximately constant in energy and = 3.07 T_eff from 2.76 to 5.02 TeV, validated only with the soft EOS. The single stiff-EOS test changes /T_eff by 6% and by 40%, so it cannot bound the sensitivity at the quoted +/- 9 MeV and +/- 0.04 level. This leaves the claimed 'agreement with ab-initio calculations' dependent on unquantified EOS sensitivity in the calibration.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a method to extract the effective temperature, entropy density, and speed of sound of the quark-gluon plasma produced in ultrarelativistic heavy-ion collisions. Using hydrodynamic simulations with the soft lattice-QCD equation of state s95p-v1, the authors establish that the mean transverse momentum of charged hadrons is proportional to an effective freeze-out temperature, = 3.07 Teff, across centralities, collision energies, and different transport coefficients. They then combine this calibration with the measured = 681 MeV to obtain Teff = 222 ± 9 MeV, use the charged-particle multiplicity and the hydrodynamic effective volume to obtain s(Teff) = 20 ± 5 fm^-3 (i.e., s/T^3 = 14 ± 3.5), and use the relative variation of and dNch/dη between 2.76 and 5.02 TeV to obtain c_s^2 = 0.24 ± 0.04. The paper concludes that these results agree with lattice QCD and confirm the production of a deconfined phase.","tokens_in":11317,"tokens_out":7333,"duration_ms":72665,"significance":"If the extraction is robust, the paper offers an unusually direct experimental route to the thermodynamics of the quark-gluon plasma, with the attractive feature that many model-dependent normalization factors cancel in ratios. The systematic exploration of ideal, shear-viscous, and bulk-viscous hydrodynamics, as well as the use of both soft and stiff equations of state, is a real strength, and the authors make their hydrodynamic code publicly available. The speed-of-sound extraction from the energy dependence of and multiplicity is elegant and is the most model-independent part of the paper. However, the entropy-density result is largely circular as presented, and the quoted uncertainties do not yet include the equation-of-state sensitivity of the central calibration. These issues affect the main quantitative claims and therefore need to be addressed before the paper can be accepted.","major_comments":[{"comment":"The reported s(Teff)/T^3 = 14 ± 3.5 is not an independent confirmation of lattice QCD as it stands. In the simulation, Veff is defined through S = s_EOS(Teff) Veff using the s95p-v1 equation of state, and the initial entropy is normalized centrality-by-centrality to the measured dNch/dη. When Eq. (4) is evaluated with this hydrodynamic Veff and the measured multiplicity, it approximately returns the input s95p-v1 entropy density, provided the hydrodynamic entropy per charged particle agrees with the assumed S/Nch = 6.7. The manuscript should report this consistency check explicitly (e.g., compare the hydrodynamic S/Nch with 6.7) and state which part of the extracted s is genuinely independent of the input equation of state. As written, the top panel of Fig. 4 demonstrates consistency of the simulation with its own equation of state rather than an ab-initio confirmation.","section":"Eq. (1) and Eq. (4), Fig. 4 (top)"},{"comment":"The quoted uncertainty Teff = 222 ± 9 MeV includes the 4% freeze-out temperature variation but excludes the 6% change in /Teff obtained with the stiff equation of state. The argument that the stiff equation of state overpredicts by 40% disfavors that particular equation of state, but it does not bound the calibration ratio for other equations of state that still reproduce the measured . Because Teff enters all subsequent thermodynamic extractions, the paper should either include an equation-of-state modeling contribution in Eq. (3) or provide a calibration test over a range of equations of state that remain compatible with the data. Without this, the uncertainty on Teff is incomplete.","section":"Eq. (3) and Fig. 2 (dot-dashed lines)"},{"comment":"The speed-of-sound extraction c_s^2 = 0.24 ± 0.04 assumes that the proportionality coefficients in Eq. (6) are energy-independent between 2.76 and 5.02 TeV. The authors estimate the Veff variation (at most 3%) but do not propagate the equation-of-state sensitivity of the -to-Teff normalization across this energy interval. Since the stiff-equation-of-state test changes this normalization by 6% in absolute value, a relative change of similar magnitude between the two energies could shift the derived c_s^2 well beyond the quoted ±0.04. Please show the stiff-equation-of-state predictions for /Teff at both 2.76 and 5.02 TeV, or otherwise justify that this effect is negligible.","section":"Eqs. (6)-(8)"}],"minor_comments":[{"comment":"The black and red curves are plotted on different vertical scales that differ by the factor 3.07, so the statement in the text that 'black lines and red lines overlap' is only true after rescaling; please clarify this in the caption to avoid confusion.","section":"Fig. 2 caption"},{"comment":"The approximation dN/dy ≃ 1.15 dN/dη is used without derivation or an error estimate; a short justification or reference would improve the reproducibility of Eq. (5).","section":"Methods"},{"comment":"The grey box is shifted 5 MeV lower than the Teff of Eq. (3) assuming a linear interpolation of the temperature between 2.76 and 5.02 TeV; please state this assumption explicitly.","section":"Fig. 4 bottom"},{"comment":"The phrase 'the first such determination' should be tempered or carefully justified, since thermal fits (Ref. [24]) already provide freeze-out temperature estimates; the novelty is better described as a simultaneous determination of the temperature with a second thermodynamic quantity from the same data.","section":"Abstract and introduction"},{"comment":"The code availability statement links to the generic MUSIC repository but not to the specific configuration files, initial conditions, or scripts used for this paper; providing these would make the results fully reproducible.","section":"Code Availability"}],"recommendation":"major_revision","confidential_remarks":"The main issue is the circularity of the s/T^3 comparison: if the authors can show, by reporting the hydrodynamic S/Nch consistency check and by quantifying the equation-of-state uncertainty in the /Teff calibration, that the extracted entropy density is not just the input s95p-v1 equation of state, the paper would be substantially strengthened. The speed-of-sound extraction is the most novel and least circular result, but its sensitivity to the equation of state across the energy interval needs to be quantified. I recommend major revision rather than rejection because the method is promising and the issues appear addressable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper before the next heavy-ion seminar. The authors take the old Van Hove idea that tracks temperature and give it a modern, quantitative footing: they calibrate = 3.07 T_eff in full hydrodynamic simulations and use it to turn ALICE data into T_eff = 222 ± 9 MeV, s/T^3, and c_s^2 = 0.24 ± 0.04. The speed-of-sound extraction is the real new physics. The entropy-density result is not an independent check.\n\nWhat is genuinely good: the method is clean, the simulations are standard and public, and the observation that /T_eff is stable across centralities, collision energies, and different viscosity settings is a useful empirical fact. The paper is honest about freeze-out uncertainties and does not oversell transport conclusions. The c_s^2 measurement, obtained by taking ratios of and dN/dy between 2.76 and 5.02 TeV, is a new handle on the QGP equation of state that is largely free of the usual model dependence.\n\nNow the soft spot, and it is load-bearing. The s/T^3 agreement with lattice QCD is close to circular. In the simulations, V_eff is defined through S = s(T_eff)V_eff, and the entropy is tuned to the measured multiplicity per centrality. When you then insert the experimental dN/dy into Eq. (4), you are mostly recovering the s95p-v1 input equation of state. The paper does not flag this. The c_s^2 result is more independent because relative variations cancel V_eff and the calibration constant, but the assumption that V_eff is energy-independent is tested only within the same soft EOS framework. The single stiff-EOS run changes by 40% and /T_eff by only 6%, which is reassuring but does not justify the ±9 MeV uncertainty without a more systematic scan. The back-of-envelope g~30 estimate is fine as an illustration but adds nothing.\n\nWho should read this: heavy-ion phenomenologists, lattice practitioners who want to see their EOS confronted with data, and anyone doing Bayesian constraints on transport or EOS. The paper deserves serious peer review. The referee should ask for an explicit EOS-sensitivity study of the /T_eff calibration and a reframing of the entropy density as a consistency check rather than a confirmation.\n\nBottom line: accept conditional on those revisions. The c_s result alone justifies the paper.","headline":"A genuinely new method for extracting the QGP effective temperature and speed of sound from LHC data, but the entropy-density 'confirmation' of lattice QCD is largely circular.","tokens_in":11927,"tokens_out":1863,"would_cite":true,"duration_ms":21643,"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":"A single momentum observable fixes the temperature of the quark-gluon plasma created in heavy-ion collisions.","keywords":["quark-gluon plasma","heavy-ion collisions","effective temperature","lattice QCD equation of state","speed of sound","entropy density","mean transverse momentum"],"falsifier":"Compare the speed-of-sound value obtained from the ratio of percent changes in mean transverse momentum and multiplicity across a scan of collision energies to the lattice QCD prediction: a mismatch larger than the stated 0.04 uncertainty at any one energy would invalidate the universal calibration.","tokens_in":10744,"feed_emoji":"🔥","tokens_out":7330,"duration_ms":68415,"temperature":0.7,"pith_summary":"The paper claims the first determination of a temperature and a second thermodynamic quantity from heavy-ion collision data, using the measured mean transverse momentum to read off an effective temperature. It derives $T_{\\rm eff} = 222 \\pm 9$ MeV in central lead-lead collisions at 5.02 TeV, and from the same data an entropy density $s(T_{\\rm eff}) = 20 \\pm 5$ fm$^{-3}$ and a squared speed of sound $c_s^2 = 0.24 \\pm 0.04$. These values agree with lattice QCD calculations, which is why a careful reader would care: a direct, almost parameter-free experimental thermometer for the deconfined phase would establish that the matter created in these collisions is a quark-gluon plasma with the expected number of degrees of freedom.","feed_headline":"One momentum observable now fixes the quark-gluon plasma's temperature","feed_subtitle":"Charged-particle momentum alone gives 2.58 trillion kelvin, plus entropy density and sound speed that match lattice QCD.","key_machinery":"The central object is the effective temperature $T_{\\rm eff}$ and effective volume $V_{\\rm eff}$, defined by $E = \\epsilon(T_{\\rm eff}) V_{\\rm eff}$ and $S = s(T_{\\rm eff}) V_{\\rm eff}$ from the integrals of the stress-energy tensor and entropy current over the freeze-out hypersurface. Taking the ratio $E/S$ eliminates $V_{\\rm eff}$, so $T_{\\rm eff}$ is fixed by the equation of state alone. The load-bearing identity is the hydrodynamic proportionality $\\langle p_t\\rangle = 3.07\\, T_{\\rm eff}$ with the soft lattice equation of state, which converts a measured momentum into a temperature, and the thermodynamic identity $c_s^2 = d\\ln\\langle p_t\\rangle/d\\ln(dN_{\\rm ch}/d\\eta)$ for the speed of sound.","core_discovery":"The paper establishes that the average transverse momentum of charged particles, $\\langle p_t\\rangle$, is proportional to an effective temperature $T_{\\rm eff}$ defined by equating the total energy and entropy at freeze-out to those of a uniform fluid at rest, using the same lattice-QCD equation of state as in the hydrodynamic simulation. In simulations with a soft equation of state, the proportionality constant is $\\langle p_t\\rangle/T_{\\rm eff} = 3.07$ and is nearly insensitive to shear or bulk viscosity, centrality, and collision energy. Using the measured $\\langle p_t\\rangle = 681$ MeV therefore yields $T_{\\rm eff} = 222 \\pm 9$ MeV, and combining the charged multiplicity with the effective volume gives $s(T_{\\rm eff}) = 20 \\pm 5$ fm$^{-3}$; forming the logarithmic derivative of $\\langle p_t\\rangle$ with respect to multiplicity between two collision energies gives $c_s^2 = 0.24 \\pm 0.04$. The paper presents these results as agreeing with ab initio lattice QCD and as evidence that a deconfined phase with active colour degrees of freedom is produced.","pith_inferences":["If the 3.07 calibration is universal, this turns a single high-statistics observable into a thermometer, and the same procedure could be applied to smaller systems such as proton-nucleus collisions where hydrodynamic analyses are less settled.","The energy-derivative route to $c_s^2$ suggests a direct experimental search for a softest point in the equation of state: scan $\\langle p_t\\rangle$ as a function of collision energy and look for a minimum in its slope with respect to multiplicity.","A test worth checking is whether the constant 3.07 survives a different hydrodynamic code or a different tuned initial-condition model; if it shifts, the quoted temperature would move with it, even though the paper's internal checks cover shear viscosity, bulk viscosity, and the stiff equation of state."],"forward_implications":["The measured temperature $T_{\\rm eff}=222\\pm9$ MeV ($2.58\\times10^{12}$ K) gives an experimental anchor for the quark-gluon plasma equation of state.","The entropy density $s/T^3 = 14\\pm3.5$, compared with lattice QCD, indicates roughly 30 active degrees of freedom, i.e., deconfined quarks and gluons.","The speed of sound $c_s^2=0.24\\pm0.04$, about half the speed of light squared, constrains the equation of state in the crossover region.","Because $T_{\\rm eff}$ is nearly independent of centrality, all centralities at 5.02 TeV probe essentially the same temperature, simplifying comparisons across collision systems.","If $\\langle p_t\\rangle$ were measured at lower beam energies, the same derivative method would extend the thermodynamic curve to lower temperatures and nonzero baryon chemical potential."],"supporting_citations":[{"why":"Supplies the experimental $\\langle p_t\\rangle = 681$ MeV that anchors $T_{\\rm eff}$.","marker":"[13]"},{"why":"Provides the ab initio lattice QCD equation-of-state results used for comparison.","marker":"[3]"},{"why":"Supplies the soft lattice-based equation of state used in the simulations that yield the 3.07 factor.","marker":"[17]"},{"why":"Supplies the charged multiplicity at 5.02 TeV used for entropy density and speed-of-sound estimates.","marker":"[15]"},{"why":"Supplies the entropy per charged particle, $S/N_{\\rm ch}=6.7\\pm0.8$, converting multiplicity to entropy density.","marker":"[25]"},{"why":"Supplies the charged multiplicity at 2.76 TeV needed for the energy derivative.","marker":"[26]"},{"why":"Supplies the transverse momentum spectra used to obtain the relative change of $\\langle p_t\\rangle$ between collision energies.","marker":"[27]"},{"why":"Supplies the tuned initial-condition parametrization whose transverse radius fixes the effective volume.","marker":"[21]"},{"why":"Supplies the hydrodynamic code in which the $\\langle p_t\\rangle$-to-$T_{\\rm eff}$ proportionality is established.","marker":"[42-44]"},{"why":"Supplies the bulk-viscosity parametrization used to test whether the proportionality survives transport effects.","marker":"[16]"}],"fun_headline_variants":["First temperature and entropy of quark-gluon plasma measured","Transverse momentum alone sets quark-gluon plasma temperature","Experiments pin down QGP temperature, entropy, and sound speed","Heavy-ion collisions yield first QGP temperature and entropy","First measurement of QGP temperature and entropy matches lattice QCD"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The measured average sideways momentum of the produced particles is assumed to be locked to the true temperature by a factor of 3.07, a lock calibrated in one computer model of the collision; if that model's equation of state or freeze-out description is wrong for the real plasma, every extracted number shifts.","fun_headline_variants_meta":{"raw":{"variants":["First temperature and entropy of quark-gluon plasma measured","Transverse momentum alone sets quark-gluon plasma temperature","Experiments pin down QGP temperature, entropy, and sound speed","Heavy-ion collisions yield first QGP temperature and entropy","First measurement of QGP temperature and entropy matches lattice QCD"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00074,"raw_usage":{"total_tokens":3327,"prompt_tokens":989,"completion_tokens":2338,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":2254}},"tokens_in":605,"tokens_out":2338,"duration_ms":15983,"temperature":1.0,"reasoning_tokens":2254,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:03:53.025320+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the speed-of-sound value obtained from the ratio of percent changes in mean transverse momentum and multiplicity across a scan of collision energies to the lattice QCD prediction: a mismatch larger than the stated 0.04 uncertainty at any one energy would invalidate the universal calibration.","supporting_citations":[{"cited_title":"Transverse momentum spectra and nuclear modification factors of charged particles in Xe-Xe collisions at √sNN = 5.44 TeV,","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental $\\langle p_t\\rangle = 681$ MeV that anchors $T_{\\rm eff}$."},{"cited_title":"Full result for the QCD equation of state with 2+1 flavors,","cited_arxiv_id":null,"evidence_quote":"Provides the ab initio lattice QCD equation-of-state results used for comparison."},{"cited_title":"QCD Equation of State and Hadron Resonance Gas,","cited_arxiv_id":null,"evidence_quote":"Supplies the soft lattice-based equation of state used in the simulations that yield the 3.07 factor."},{"cited_title":"Centrality de- pendence of the charged-particle multiplicity density at midrapidity in Pb-Pb collisions at √sNN = 5.02 TeV,","cited_arxiv_id":null,"evidence_quote":"Supplies the charged multiplicity at 5.02 TeV used for entropy density and speed-of-sound estimates."},{"cited_title":"Entropy production in pp and Pb-Pb collisions at energies avail- able at the CERN Large Hadron Collider,","cited_arxiv_id":null,"evidence_quote":"Supplies the entropy per charged particle, $S/N_{\\rm ch}=6.7\\pm0.8$, converting multiplicity to entropy density."},{"cited_title":"Centrality de- pendence of the charged-particle multiplicity density at mid-rapidity in Pb-Pb collisions at √sN N = 2.76 TeV,","cited_arxiv_id":null,"evidence_quote":"Supplies the charged multiplicity at 2.76 TeV needed for the energy derivative."},{"cited_title":"Transverse momentum spectra and nuclear modification factors of charged particles in pp, p-Pb and Pb-Pb collisions at the LHC,","cited_arxiv_id":null,"evidence_quote":"Supplies the transverse momentum spectra used to obtain the relative change of $\\langle p_t\\rangle$ between collision energies."},{"cited_title":"Hydrodynamic predictions for 5.44 TeV Xe+Xe collisions,","cited_arxiv_id":null,"evidence_quote":"Supplies the tuned initial-condition parametrization whose transverse radius fixes the effective volume."},{"cited_title":"Applying Bayesian parameter estimation to relativistic heavy-ion collisions: simultaneous char- acterization of the initial state and quark-gluon plasma medium,","cited_arxiv_id":null,"evidence_quote":"Supplies the bulk-viscosity parametrization used to test whether the proportionality survives transport effects."}],"review_version":1}