{"id":"efe9edc1-90a9-44e5-9b2f-c82545a244b0","arxiv_id":"1909.00643","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Kinetic freeze-out in central Au+Au collisions is a continuous process spanning broad ranges of temperature and baryon chemical potential, with averages nearly flat in transverse momentum and rapidity.","lead":"Using a transport model with coarse-graining, the authors map where and when hadrons stop interacting in central gold-gold collisions at five beam energies, extracting the temperature and baryon chemical potential at kinetic freeze-out. The work shows freeze-out is a gradual, spread-out process rather than a single point, which matters for how experimental data are interpreted with Blast-Wave fits.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quantitative T/µB values hinge on mapping non-equilibrium UrQMD cells onto a full-equilibrium HRG EoS; the paper concedes this may bias results, and no test quantifies the shift.","rationale":"I read the paper as a transparent model study whose main qualitative message, that kinetic freeze-out is a continuous process with a wide decoupling-time distribution, is well supported by the UrQMD last-interaction samples. That part does not depend on the EoS mapping. The load-bearing weakness is the conversion of coarse-grained energy density and net-baryon density into thermodynamic variables via a full-equilibrium HRG EoS. The paper explicitly states that the EoS may introduce a bias when the system is out of chemical equilibrium, so the concern is not hidden, but it is also not quantified. Because the paper's most impactful quantitative comparison, the claimed tension with Blast-Wave temperatures, uses these mapped values, the absence of a sensitivity test means the comparison is conditional. The reader's weakest_assumption identifies the same local-equilibrium step, and the reader's CONDITIONAL verdict is appropriate. A dedicated re-extraction with a chemically non-equilibrium EoS would settle whether the reported averages and the Blast-Wave discrepancy are robust or a modeling artifact.","tokens_in":14628,"tokens_out":10449,"duration_ms":103923,"concrete_test":"Re-extract ⟨T⟩ and ⟨µB⟩ at √sNN=200 GeV from the same coarse-grained cells using a partial-chemical-equilibrium HRG EoS, e.g. with a pion chemical potential or species fugacities adjusted cell-by-cell to reproduce UrQMD's hadron abundances, instead of the full-equilibrium EoS of Ref. [63]. Compare the resulting ⟨T⟩ with the Blast-Wave point in Fig. 11; if the shift exceeds about 10 MeV, the claimed tension is an artifact of the equilibrium assumption rather than a robust physical difference. As a second arm, rerun the extraction including the roughly 15% of freeze-out hadrons dropped at 200 GeV by the 100-particle cell cut and check whether the quoted averages move outside the claimed 10% band.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The values that drive the comparison with Blast-Wave fits (Fig. 11) are produced by the two-step procedure of Sec. II: coarse-grained (ε, ρB) from UrQMD cells are mapped to (T, µB) through a full-equilibrium Hadron Resonance Gas EoS. UrQMD is a non-equilibrium cascade, and nothing in the extraction checks that the cell-level momentum or species distributions are thermally and chemically equilibrated. The paper itself concedes in Sec. IV that the chosen EoS is inadequate for a system out of chemical equilibrium and 'might introduce a bias'. If the cell abundances deviate from HRG equilibrium, the same (ε, ρB) does not single out a unique T/µB, so the reported ⟨T⟩, the claim of no more than 10% variation with pT, and the high-energy tension with Blast-Wave fits are all conditional on an unquantified bias. The qualitative central claim survives because the broad freeze-out time distribution follows directly from the last-interaction samples, but the quantitative comparison to Blast-Wave fits is not secure without testing the equilibrium mapping.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extracts kinetic freeze-out temperatures and baryon chemical potentials from UrQMD cascade simulations combined with a coarse-graining procedure. For central Au+Au collisions at five beam energies from 2.4 to 200 GeV, the authors define kinetic freeze-out as the last interaction point of each hadron, compute coarse-grained energy density and net-baryon density around those points, and map the resulting (ε, ρB) values to (T, μB) using a tabulated hadron resonance gas equation of state. They present freeze-out time, temperature, and chemical-potential distributions, study the transverse-momentum and rapidity dependence of the averaged freeze-out parameters, and compare their average (T, μB) points with Blast-Wave fit results and statistical hadronization model results. The central qualitative claim is that kinetic freeze-out is a continuous process occurring over 15–20 fm/c and over broad ranges of T and μB, so that representing it by a single point in the QCD phase diagram is an oversimplification. The paper also claims a tension with Blast-Wave fits, particularly at high collision energy.","tokens_in":14803,"tokens_out":4173,"duration_ms":42963,"significance":"If the quantitative extraction is accepted, the paper provides a useful phenomenological benchmark showing that single-point freeze-out parameters hide a substantial spread in decoupling conditions, and it offers a direct microscopic definition of the kinetic freeze-out hypersurface in a transport approach. Strengths of the paper are its large event statistics, the transparent step-by-step presentation of the coarse-graining and mapping procedure, and the explicit listing of limitations in Section IV, including the equilibrium-EoS assumption and the absence of a particlization stage. The qualitative claim about the continuous nature of freeze-out is robust because it follows directly from the last-interaction samples, whereas the quantitative claims about average T and μB, the 10% flatness statement, and the comparison with Blast-Wave fits depend on the unvalidated equilibrium mapping of non-equilibrium cascade cells.","major_comments":[{"comment":"The central quantitative step is the mapping of coarse-grained (ε, ρB) from non-equilibrium UrQMD cells to unique (T, μB) via the hadron resonance gas EoS. This assumes local thermal and chemical equilibrium in every cell, but UrQMD is a non-equilibrium cascade and no test in the paper checks whether cell-level momentum or species distributions are thermal. The manuscript itself concedes in Section IV that the EoS 'might introduce a bias' for a system out of chemical equilibrium. Because this mapping feeds directly into the average ⟨T⟩, the pT- and y-dependence claims, and the Blast-Wave comparison in Fig. 11, the authors should quantify the bias, for example by comparing coarse-grained cell momentum distributions with thermal fits, by repeating the extraction with a lattice-based EoS at high T, or by restricting the analysis to cells where an equilibrium criterion is satisfied. Without such a test, the quantitative values are conditional on an unquantified assumption.","section":"Section II, Eqs. (2)–(5) and interpolation step"},{"comment":"The claimed tension with Blast-Wave fits at high collision energy is interpreted as weaker transverse expansion in the cascade. However, Blast-Wave Tkin is a fit parameter strongly correlated with the fitted transverse velocity and with the selected pT ranges, so a direct comparison of the present average freeze-out temperature with the Blast-Wave parameter is not apples-to-apples. The paper notes the pT cuts and resonance feed-down as possible sources of difference but does not quantify them. A quantitative robustness test is needed, for instance by applying the same Blast-Wave fitting procedure to UrQMD-generated spectra and comparing the resulting Tkin with the coarse-grained ⟨T⟩, or by varying the pT ranges used in the data fits. This is load-bearing for the paper's comparison and for the interpretation that the cascade produces weaker radial flow than the data.","section":"Section III.C, Fig. 11"},{"comment":"The sentence 'except at √sNN = 200 GeV, where we drop at ≈ 85%' is ambiguous and potentially serious. If 85% of the kinetically frozen-out hadrons are discarded because their cells have fewer than 100 particles, the 200 GeV results in Figs. 4–11 are based on only about 15% of the freeze-out sample, which could bias the averages and distributions toward high-density central cells. The authors should clarify the retention fraction and, if the drop is indeed 85%, quantify the selection bias by comparing the accepted and rejected samples or by lowering the acceptance threshold with a statistical correction.","section":"Section II, coarse-graining acceptance at √sNN=200 GeV"},{"comment":"The claim that the average freeze-out temperature and baryon chemical potential are 'essentially independent of rapidity and transverse momentum' is too strong as stated. The paper itself notes that at √sNN = 2.4 GeV the pT dependence is not flat (Figs. 4–6 show a pronounced rise of ⟨T⟩ and ⟨μB⟩ with pT at the lowest energy), and the conclusion is already qualified in the body of Section III.B. The abstract and conclusions should either explicitly restrict the flatness claim to √sNN ≥ 4.5 GeV or quantify the size of the deviations at each energy, since the 10% statement is one of the paper's quantitative results.","section":"Section III.B and Section IV"}],"minor_comments":[{"comment":"There are several typos and inconsistencies: 'for for high collision energies' in Section III.A, 'UrMQD' in Section IV, the comma in '√sNN = 2.4, GeV' in Section III.A, and the axis label 'T [Mev]' in Fig. 10. These should be corrected.","section":"Throughout"},{"comment":"The text refers to √sNN = 7 GeV when the simulated energies are 7.7 GeV and 19.6 GeV; please make the numbers consistent.","section":"Section III.C"},{"comment":"The caption of Fig. 12 describes a density profile in the (T, μB) plane but does not state the axis labels or the meaning of the color scale. Please clarify what is plotted and how the density is normalized.","section":"Fig. 12 caption"},{"comment":"Reference [49] appears incomplete (no year or journal information is given), and Refs. [45] and [49] share the same authors; please provide full bibliographic information and check for duplication.","section":"References"},{"comment":"The statement that the peak emission temperature 'does not rise above a certain threshold of approx. 150 MeV' is based on the HRG EoS used in the mapping; this should be stated together with the caveat that the EoS itself restricts the accessible temperature range, as is already partially acknowledged later in the same paragraph.","section":"Section III.A"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest and within the scope of J. Phys. G; the qualitative freeze-out-continuity message is well supported by the last-interaction distributions. The main issue is the unquantified equilibrium-mapping bias in the quantitative T/μB extraction and the Blast-Wave comparison, plus the ambiguous 200 GeV acceptance statement. These are fixable with additional validation and a clearer presentation, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful thing here is the systematic five-energy extraction of kinetic freeze-out T and muB from the UrQMD coarse-graining framework, including pT and rapidity dependence. The qualitative central claim—that freeze-out is a continuous, distributed process rather than a single phase-diagram point—is well supported by the last-interaction samples: decoupling time distributions are broad with FWHM 15-20 fm/c, and the T/muB distributions are broad too. I buy that. The paper is honest about what it is; it does not oversell.\n\nWhere it gets softer: the quantitative T and muB values rest on mapping non-equilibrium cascade cells onto a full-equilibrium HRG EoS. The authors themselves concede in Sec. IV that the EoS may be inadequate out of chemical equilibrium and might bias results. That is a real limitation, and it matters for the comparison to Blast-Wave fits in Fig. 11. They present that comparison without error bars on the model points, and they do not test sensitivity to the EoS or cell acceptance. The high-energy tension with Blast-Wave—their most impactful quantitative statement—is therefore conditional on an unquantified mapping bias. The stress-test concern lands here. It is not fatal to the main message, because the continuous nature of freeze-out is evident directly from the last-interaction coordinates and broad distributions; but if the highlighted result is the tension with Blast-Wave, that piece is not secure as it stands.\n\nAlso note the claim of no more than 10% variation with pT and rapidity is stated without statistical uncertainties. The energy-ordering of the curves is clear, so this is a moderate, not fatal, weakness.\n\nCitation pattern looks fine; the authors build on their own coarse-graining refs and cite the prior Bravina work. The paper is aimed at heavy-ion phenomenologists working on freeze-out. It deserves a serious referee. Conditional acceptance with requests for uncertainty estimation and an EoS sensitivity check would be the right path.","headline":"Useful five-energy map of kinetic freeze-out from UrQMD coarse-graining; the continuous-freeze-out claim holds up, but the quantitative tension with Blast-Wave is model-dependent and not yet error-barred.","tokens_in":15367,"tokens_out":2153,"would_cite":true,"duration_ms":79226,"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":"Kinetic freeze-out is not a single point in the phase diagram","keywords":["kinetic freeze-out","coarse graining","heavy ion collisions","baryon chemical potential","QCD phase diagram","hadron resonance gas","UrQMD","freeze-out hypersurface"],"falsifier":"Measure the momentum distribution of two hadron species (say pions and protons) at the same coarse-grained freeze-out cell and check whether both are described by the same temperature and chemical potential; if the species-dependent temperatures differ by more than the claimed 10% at any cell, the equilibrium mapping fails. Alternatively, if rerunning the analysis with a lattice-based equation of state or with a cascade that includes a hydrodynamically evolving quark-gluon plasma stage changes the average kinetic freeze-out temperature by more than about 10%, the tension with Blast-Wave fits would not be robust.","tokens_in":1966,"feed_emoji":"💥","tokens_out":5205,"duration_ms":332678,"temperature":0.7,"pith_summary":"This paper sets out to show that kinetic freeze-out in relativistic heavy-ion collisions cannot be summarized as a single point on the QCD phase diagram. Using a hadron cascade model together with a coarse-graining procedure, the authors track the last interaction of each hadron and assign the local temperature and baryon chemical potential of the surrounding medium at that space-time point. For central Au+Au collisions from $\\sqrt{s_{NN}}=2.4$ GeV to 200 GeV, they find that hadrons decouple continuously over a time window of 15-20 fm/c (full width at half maximum) and over broad ranges of temperature and baryon chemical potential, even though the average decoupling temperature and chemical potential vary by less than about 10% across transverse momentum and show only weak rapidity dependence. The paper therefore argues that quoting a single freeze-out point from Blast-Wave fits is an oversimplification, and it highlights a tension: at high collision energy the kinetic freeze-out temperature extracted here is higher than in Blast-Wave fits.","feed_headline":"Kinetic freeze-out is a continuous process, not a single point","feed_subtitle":"In central Au+Au collisions, hadrons decouple over 15-20 fm/c while average T and muB vary less than 10% with pT.","key_machinery":"The argument rests on a coarse-graining procedure that converts the microscopic output of a hadron cascade into continuum thermodynamic fields. At fixed times, the net-baryon four-current and energy-momentum tensor are computed in cells of size roughly 0.8 fm, the Eckart frame is used to define the fluid velocity, and the local rest-frame energy density $\\varepsilon$ and net-baryon density $\\rho_B$ are mapped to temperature and baryon chemical potential through a tabulated Hadron Resonance Gas equation of state with the same degrees of freedom as the cascade. These continuum values are then assigned to the space-time point of each hadron's last interaction, defining the kinetic freeze-out hypersurface. This machinery allows the paper to attach thermodynamic meaning to the decoupling distribution and to compare the resulting averages with Blast-Wave fits.","core_discovery":"The central discovery is that kinetic freeze-out is a continuous, dynamical process rather than a sharp moment. On the decoupling hypersurface defined by the last interaction of each hadron (including strong decays), the temperature and baryon chemical potential are spread over wide ranges whose widths reflect the expansion dynamics and energy-dependent cross sections: for example, at intermediate beam energies the baryon chemical potential distribution is particularly broad because the system transitions from baryon-dominated to meson-dominated matter. Averaged over the hypersurface, however, the temperatures and chemical potentials order smoothly with collision energy, with the average temperature rising and the average baryon chemical potential falling as $\\sqrt{s_{NN}}$ increases from 2.4 to 200 GeV, and the averaged values vary by less than 10% with transverse momentum and only mildly with rapidity. The paper further claims that the average kinetic freeze-out temperature at high energy exceeds the Blast-Wave fitted value, which it attributes to the weaker transverse expansion generated by the pure hadron cascade without a hydrodynamic or quark-gluon plasma stage.","pith_inferences":["The same coarse-grained freeze-out extraction could be applied species-by-species (pions, kaons, protons) to test whether different hadrons decouple at systematically different temperatures and to search for the species-splitting that a single freeze-out temperature cannot capture.","The broad spread in $(T,\\mu_B)$ at freeze-out implies that inclusive momentum-space observables average over many different thermodynamic conditions; this may contribute to apparent temperatures extracted from slope ratios and may be connected to transverse-momentum-dependent fluctuations.","If the tension with Blast-Wave fits stems from the absence of a quark-gluon plasma stage, extending the analysis to hybrid transport-hydrodynamic models at RHIC and LHC energies would be a direct test; the method could also be applied to smaller systems such as pp or pA collisions.","The claimed time-extended decoupling of 15-20 fm/c suggests that two-particle correlations and interferometry measurements carry information about the freeze-out duration; comparing the coarse-grained emission function with measured correlations would provide an independent check."],"forward_implications":["Freeze-out parameters should be reported as distributions on a decoupling hypersurface rather than as a single point in the $(T,\\mu_B)$ plane if one wants to capture the actual dynamics.","The weak transverse-momentum and rapidity dependence of the averages means that a single average $(\\langle T\\rangle, \\langle\\mu_B\\rangle)$ per collision energy remains a useful summary despite the underlying spread.","At high collision energy the kinetic freeze-out temperature is higher than Blast-Wave fits suggest, indicating that the hadron cascade produces weaker transverse expansion than the data; this points to the need for a hydrodynamic or particlization stage or a reassessment of Blast-Wave fit assumptions.","The separation between chemical and kinetic freeze-out temperatures grows from roughly 5-10 MeV at low energies to more than 40-50 MeV at high energies, quantifying the strength of the expansion between the two stages."],"supporting_citations":[{"why":"Supplies the hadron cascade model used both to generate the event ensemble for coarse graining and to determine the space-time points of last interaction.","marker":"[52, 53]"},{"why":"Develops the coarse-graining method that reconstructs thermal fields from the cascade output, the technical basis of the extraction.","marker":"[55]"},{"why":"Provides the tabulated Hadron Resonance Gas equation of state used to convert local energy and baryon densities into temperature and chemical potential.","marker":"[63]"},{"why":"Gives the Blast-Wave fit results at 200 GeV that the paper compares against its own kinetic freeze-out temperature.","marker":"[46]"},{"why":"Provides Blast-Wave fit kinetic freeze-out temperatures at 7.7 and 19.6 GeV and chemical freeze-out temperatures used in the comparison.","marker":"[48]"},{"why":"Supplies Blast-Wave fit kinetic freeze-out temperatures at low collision energies (2.7-4.3 GeV) used for comparison.","marker":"[50]"},{"why":"Gives a prior Kadanoff-Baym based analysis whose freeze-out time distribution the paper finds consistent with its own continuous decoupling picture.","marker":"[21]"},{"why":"Provides chemical freeze-out temperatures from statistical hadronization model fits at low energies used for comparing with kinetic freeze-out.","marker":"[74]"}],"fun_headline_variants":["Freeze-out is a process, not a single point","Kinetic freeze-out is gradual, spanning 15-20 fm/c","Average T and muB shift <10% with pT at freeze-out","Au+Au freeze-out: continuous decoupling, not a snapshot"],"cache_read_input_tokens":17536,"weakest_assumption_plain":"The whole extraction assumes that the local mix of hadrons at a freeze-out cell is close enough to thermal equilibrium that a single temperature and baryon chemical potential, read off from an equilibrium hadron-resonance-gas equation of state, can meaningfully describe the cell's energy and baryon densities.","fun_headline_variants_meta":{"raw":{"variants":["Freeze-out is a process, not a single point","Kinetic freeze-out is gradual, spanning 15-20 fm/c","Average T and muB shift <10% with pT at freeze-out","Au+Au freeze-out: continuous decoupling, not a snapshot"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000296,"raw_usage":{"total_tokens":1735,"prompt_tokens":976,"completion_tokens":759,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":683}},"tokens_in":592,"tokens_out":759,"duration_ms":13896,"temperature":1.0,"reasoning_tokens":683,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:40:44.915205+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the momentum distribution of two hadron species (say pions and protons) at the same coarse-grained freeze-out cell and check whether both are described by the same temperature and chemical potential; if the species-dependent temperatures differ by more than the claimed 10% at any cell, the equilibrium mapping fails. Alternatively, if rerunning the analysis with a lattice-based equation of state or with a cascade that includes a hydrodynamically evolving quark-gluon plasma stage changes the average kinetic freeze-out temperature by more than about 10%, the tension with Blast-Wave fits would not be robust.","supporting_citations":[{"cited_title":"Zschiesche, S","cited_arxiv_id":null,"evidence_quote":"Provides the tabulated Hadron Resonance Gas equation of state used to convert local energy and baryon densities into temperature and chemical potential."}],"review_version":1}