{"id":"f84c7e97-da44-4eaf-9e14-3e384861696a","arxiv_id":"2608.10238","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"From STAR proton fluctuation data, a Bayesian extraction constrains chi_2^B at freeze-out, matching lattice QCD below mu_B about 300 MeV and attributing the nonmonotonic C3/C1 to two-baryon correlations plus baryon conservation.","lead":"This paper builds a simulation chain from hydrodynamic expansion to detector acceptance to infer the baryon number susceptibilities of hot QCD matter from measured proton fluctuation data. The extracted second-order susceptibility agrees with lattice QCD up to a baryon chemical potential of about 300 MeV, with a possible enhancement beyond that.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The extraction assumes event-averaged hypersurfaces; unmodeled volume and initial-state fluctuations at low sqrt(s_NN) could bias chi_2^B upward and generate the C3/C1 peak, making the central claim conditional.","rationale":"The paper is methodologically coherent and explicitly discloses the volume-fluctuation limitation, so the appropriate verdict remains CONDITIONAL rather than REJECT. I searched for internal inconsistencies in the MaxEnt mapping, the binomial-thinning reduction, and the two-baryon-dominance derivation in the Supplemental Material and found none: the induced third- and fourth-order susceptibilities in the bDelta chi_3 = bDelta chi_4 = 0 scenario follow algebraically from Eqs. (A6)-(A11). The central vulnerability is the event-averaged hypersurface assumption, which the authors themselves flag as requiring revisit. Because the low-energy points where the data-driven extraction deviates from lattice are precisely where volume fluctuations grow, this assumption is load-bearing for the two headline results: the chi_2^B enhancement at mu_B > 300 MeV and the interpretation of the C3/C1 peak as a chi_2-plus-conservation effect. A concrete volume-fluctuation-corrected extraction would settle whether these results are physical or artifacts. This matches the reader's weakest assumption, so no verdict change is needed.","tokens_in":15792,"tokens_out":6008,"duration_ms":60190,"concrete_test":"Use the STAR measured charged-particle multiplicity distribution in the 0-5% centrality class to construct a volume-fluctuation correction to the proton factorial cumulants, extending the centrality-bin-width formalism [50] to fixed-centrality fluctuations via a gamma-distributed volume with variance matched to the observed multiplicity width. Apply this correction to the data and rerun the Bayesian extraction of chi_2^B at sqrt(s_NN)=7.7 GeV. If the corrected gamma_1 shifts from ~1.06 toward 1.0 or below, beyond the quoted 3.5% uncertainty, the low-energy enhancement and the two-baryon-dominance mechanism do not survive the volume-fluctuation concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The forward model in Eqs. (1)-(3) and the SAM-3.0 correction are built on event-averaged hydrodynamic hypersurfaces. The STAR 0-5% data are used after only the centrality-bin-width correction [50], so event-by-event volume fluctuations and initial-state/baryon-stopping fluctuations are not modeled. These fluctuations contribute directly to the measured factorial cumulants and can mimic or cancel a genuine susceptibility enhancement. The paper itself states they 'are expected to become increasingly important at the lowest collision energies' and that the extraction 'should be revisited' if more complete volume-fluctuation corrections revise the cumulants. The low-energy enhancement of chi_2^B/bar chi_2^B at 7.7 GeV (~1.06-1.07) and the nonmonotonic C3/C1 peak are exactly in this regime. If the neglected volume fluctuations generate the observed enhancement, the claimed agreement with lattice for mu_B <= 300 MeV becomes contingent, and the two-baryon-dominance explanation of C3/C1 (interplay of chi_2^B and conservation) would be replaced by an artifact.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a Bayesian extraction of baryon number susceptibilities from STAR measurements of proton factorial cumulant ratios and net-proton cumulant ratios in 0–5% central Au–Au collisions at √sNN = 7.7–200 GeV. The forward model combines event-averaged MUSIC hydrodynamic hypersurfaces with maximum-entropy freeze-out, binomial thinning to accepted (anti)protons, and the SAM-3.0 exact baryon-number-conservation correction. The authors report tight constraints on χ2B/¯χ2B, quantitative agreement with lattice-QCD-based estimates for μB ≲ 300 MeV, an upward deviation at larger μB, and a demonstration that the nonmonotonic energy dependence of the proton factorial cumulant ratio Ĉ 3/Ĉ 1 can be reproduced without irreducible three- or four-baryon correlations.","tokens_in":16034,"tokens_out":7345,"duration_ms":76347,"significance":"The approach is novel: it converts heavy-ion fluctuation data into a data-driven determination of QCD baryon susceptibilities along the freeze-out line, rather than comparing data with a single assumed equation of state. The forward map is transparent and the supplemental material supplies explicit formulas for the MaxEnt mapping, binomial thinning, and the SAM-3.0 correction. The two independent extractions (proton factorial cumulants and net-proton cumulants) cross-validate each other, and the variation of the switching energy density is a useful robustness check. The two-baryon-dominance scenario makes a concrete, falsifiable statement about the origin of the Ĉ 3/Ĉ 1 peak. However, the current significance is bounded by the paper's own admissions that model systematics are not propagated and that volume fluctuations, initial-state fluctuations, and baryon-stopping fluctuations are neglected; these effects are largest exactly where the most interesting deviations appear.","major_comments":[{"comment":"The quoted uncertainties propagate only experimental errors; the hydrodynamic input, initial state, transport coefficients, switching energy density, and freeze-out prescription are held fixed (Bayesian inference paragraph and footnote [40]). The abstract's claim of 'tight constraints' and the quantitative agreement with lattice-QCD estimates for μB ≲ 300 MeV are therefore conditional on this fixed setup. Because the central claim is a comparison with lattice QCD, the absence of model systematics leaves open the possibility that the agreement, or the low-energy enhancement, is partly a consequence of the chosen hydrodynamic baseline. The paper should either estimate and propagate model uncertainties or state unambiguously in the abstract that the constraints are conditional on the hydrodynamic model.","section":"Bayesian inference and Fig. 3"},{"comment":"The extraction uses event-averaged hypersurfaces and includes only the centrality-bin-width correction for volume fluctuations; the paper itself states that volume-fluctuation, initial-state, and baryon-stopping effects 'are expected to become increasingly important at the lowest collision energies' and that the extraction 'should be revisited' if more complete corrections revise the cumulants. The observed enhancement of χ2B/¯χ2B at 7.7 GeV and the nonmonotonic Ĉ 3/Ĉ 1 peak near 11 GeV lie precisely in this regime. As a result, the claimed deviation from lattice at large μB and the two-baryon-dominance interpretation of Ĉ 3/Ĉ 1 are not yet robust; the manuscript must either include a quantitative estimate of these effects or present the large-μB enhancement as an explicitly model-dependent hint rather than an extracted result.","section":"Comparison with lattice QCD, caveats paragraph"},{"comment":"The Gaussian likelihood neglects correlations among the three measured ratios because the experimental covariance matrices are unavailable. This affects the size of the 68% credible intervals and, in particular, the 2.3σ preference for negative χ4B/χ2B at 19.6 GeV ('Extracted susceptibilities' section). Since the ratios are constructed from the same event samples, their correlations are likely non-negligible; without an estimate, statements of statistical significance for the higher-order susceptibilities are not supported. The authors should either obtain or estimate the covariances, or refrain from quoting significances.","section":"Bayesian inference, footnote [40]"}],"minor_comments":[{"comment":"The phrase 'tight constraints' should be qualified as 'conditional tight constraints' given the fixed hydrodynamic setup; similarly, 'indicate an enhancement' should be flagged as model-dependent.","section":"Abstract"},{"comment":"The upper axis converts √sNN to μB using the freeze-out parametrization of Ref. [43]; this mapping carries model uncertainty and should be mentioned in the caption or in the text.","section":"Fig. 1 caption"},{"comment":"The 200 GeV entry for χ3B/χ1B uses a widened prior and the posterior 'was verified to saturate'; this entry should be clearly marked as unconstrained or removed from the table of extracted values.","section":"Table I"},{"comment":"The statement that Ĉ 3/Ĉ 1 can be described 'without irreducible three- or four-baryon correlations' is correct, but the higher-order susceptibilities in the scenario are nonzero through the induced relations in Eqs. (A7)–(A8); consider wording that makes this distinction clearer, e.g., 'without irreducible multi-baryon correlations beyond second order'.","section":"Two-baryon correlation dominance"}],"recommendation":"major_revision","confidential_remarks":"The paper is methodologically interesting and the forward model is carefully constructed, but the abstract overstates the robustness of the large-μB enhancement and the Ĉ 3/Ĉ 1 explanation given the acknowledged neglect of volume fluctuations and model systematics. The likelihood's neglect of correlations also undermines the 2.3σ significance claim. These issues are fixable by softening the claims and, ideally, adding a quantitative estimate of volume-fluctuation effects, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThis paper is worth your time. It is the first real attempt to extract baryon number susceptibilities from heavy-ion cumulant data using a full forward model: viscous hydro hypersurfaces, MaxEnt freeze-out, binomial thinning to protons, and SAM-3.0 exact baryon conservation. The machinery is serious, and the two extractions (proton factorial cumulants and net-proton cumulants) cross-validate each other. The headline result is that chi_2^B/bar chi_2^B is constrained to 5-10% and agrees with lattice-based estimates for mu_B <= 300 MeV, with an upward deviation at the lowest energies. They also show that the nonmonotonic C3/C1 peak can be reproduced without irreducible three- or four-baryon correlations: chi_2 energy dependence plus conservation does it. That is a clean, useful result, though less surprising than the abstract lets on, since chi_3 and chi_4 are barely constrained.\n\nThe soft spots are the usual ones for this kind of analysis, and the authors mostly own them. The quoted uncertainties propagate only experimental errors; the hydrodynamic setup, epsilon_sw, and freeze-out prescription are fixed. They do check epsilon_sw variation, which is decent. The likelihood neglects correlations among the three measured ratios because the covariance matrices are unavailable; that is a real but probably minor issue. The bigger problem is volume fluctuations. The STAR data use only the centrality-bin-width correction, and the forward model uses event-averaged hypersurfaces. At 7.7 GeV, where they see the enhancement, volume fluctuations and baryon-stopping fluctuations are precisely the ones that could generate a fake signal. The authors say this, but the abstract's 'tight constraints' reads too strong. If a more complete correction changes the cumulants at low energies, the central comparison with lattice shifts.\n\nI also note the self-citations are fine; the SAM-3.0 and hydro references are relevant, and the lattice comparison uses external calculations. No circularity problem.\n\nWho is this for? Anyone working on BES cumulants, freeze-out, or EoS constraints. It deserves a serious referee. My recommendation is to send to peer review, but with the caveat that the paper should either soften the abstract or add a more thorough volume-fluctuation sensitivity test before acceptance.","headline":"A credible, genuinely new extraction of chi_2^B from BES data, with the low-energy enhancement still hostage to unmodeled volume fluctuations.","tokens_in":16602,"tokens_out":1831,"would_cite":true,"duration_ms":18687,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims the first Bayesian extraction of the second-order QCD baryon susceptibility from heavy-ion collision data, matching lattice QCD at baryon chemical potentials up to 300 MeV and showing a fluctuation enhancement at higher…","keywords":["baryon number susceptibility","heavy-ion collisions","net-proton cumulants","maximum-entropy freeze-out","exact baryon-number conservation","Bayesian extraction","QCD equation of state","chemical freeze-out line"],"falsifier":"If centrality-bin-width or more complete volume-fluctuation corrections revise the measured cumulants so that the extracted $\\chi_2^B/\\bar\\chi_2^B$ at the lowest collider energy drops back into the lattice band, or if the nonmonotonic peak in $\\hat{C}_3/\\hat{C}_1$ vanishes when the same framework is applied to fixed-target data, then the central claim of a data-driven low-energy enhancement and of a conservation-driven peak would be falsified.","tokens_in":15583,"feed_emoji":"⚛️","tokens_out":10598,"duration_ms":94568,"temperature":0.7,"pith_summary":"This paper argues that the event-by-event fluctuations of (net-)protons measured in central gold-gold collisions can be turned into a direct measurement of the baryon number susceptibilities of strongly interacting matter at freeze-out. It does this by building a forward model that starts from arbitrary equation-of-state susceptibilities on hydrodynamic freeze-out surfaces, passes them through maximum-entropy freeze-out, maps baryons onto accepted protons, and enforces exact baryon-number conservation. From the data, the second-order susceptibility normalized by the ideal hadron resonance gas value is constrained to 5-10 percent. The extracted values agree with lattice-QCD-based estimates along the chemical freeze-out line up to a baryon chemical potential of about 300 MeV, and rise above them at larger potentials. A striking by-product is that the nonmonotonic peak in the measured third-order proton factorial cumulant ratio can be reproduced without any irreducible three- or four-baryon correlations, as an interplay between the energy dependence of the second-order susceptibility and exact baryon-number conservation.","feed_headline":"Heavy-ion data measure baryon susceptibility along freeze-out","feed_subtitle":"A Bayesian inversion of proton fluctuation data pins the baryon susceptibility to 5-10% and agrees with lattice QCD up to 300 MeV.","key_machinery":"The load-bearing mechanism is a four-step forward map. Maximum-entropy freeze-out, the least-prejudiced assignment of hadronic fluctuations consistent with the local baryon susceptibilities, converts local equation-of-state susceptibilities, parametrized by $\\gamma=(\\chi_2^B/\\bar\\chi_2^B, \\chi_3^B/\\chi_1^B, \\chi_4^B/\\chi_2^B)$, into joint baryon-antibaryon factorial cumulants on every hydrodynamic hypersurface element. Independent binomial thinning with the local proton-to-baryon fraction and the kinematic acceptance maps those cumulants onto accepted (net-)proton cumulants, and the subensemble acceptance method (SAM-3.0) then imposes exact global baryon-number conservation. The whole chain is linear in the three susceptibility deviations before the conservation step, so a Bayesian posterior over the three ratios can be evaluated on a grid and marginalized. This is what lets the extracted $\\chi_2^B$ be compared with lattice QCD along a freeze-out line.","core_discovery":"On its own terms, the paper establishes that the measured cumulants and factorial cumulants of (net-)protons in Au-Au collisions can be inverted into the baryon number susceptibilities $\\chi_2^B$, $\\chi_3^B$, and $\\chi_4^B$ along the chemical freeze-out line, rather than being compared with pre-selected equations of state. The second-order ratio $\\chi_2^B/\\bar\\chi_2^B$ is tightly constrained and shows a non-monotonic energy dependence: near the ideal-gas value at 7.7 GeV, a suppression around 0.73-0.82 at 11.5-19.6 GeV, and a flat suppression at high energies. In absolute terms $\\chi_2^B$ agrees with lattice-based estimates for $\\mu_B \\lesssim 300$ MeV and is enhanced by up to about 50% at $\\sqrt{s_{NN}} = 7.7$ GeV. The third- and fourth-order susceptibilities remain weakly constrained, but a scenario retaining only the irreducible two-baryon correlation describes the measured $\\hat{C}_3/\\hat{C}_1$ peak, meaning the data do not require genuine multi-baryon correlations of order three or four.","pith_inferences":["A natural next step is to apply the same forward map at fixed-target energies; the extraction predicts that if the low-energy enhancement persists once volume fluctuations are controlled, it is a thermodynamic property, whereas if it moves, it was a modeling artifact.","Because the framework treats baryons and antibaryons as a two-species system with only the baryon-density mode fluctuating, an observable that is invariant under baryon-to-proton dilution and conservation corrections would be needed to isolate genuine critical correlations; the current factorial cumulant ratios are not such observables.","The same inversion strategy could be applied to strangeness or electric-charge fluctuations: data on kaon or pion cumulants could yield the corresponding susceptibilities along the freeze-out line, provided the analogous conservation laws are imposed.","If the framework were applied at another centrality or acceptance and the extracted $\\chi_2^B/\\bar\\chi_2^B$ stayed unchanged, that would validate the uniform susceptibility-ratio assumption; if it drifted, the effective single-order-parameter reduction would be too coarse."],"forward_implications":["The second-order baryon susceptibility along the chemical freeze-out line becomes a data-driven observable, and its agreement with lattice QCD at $\\mu_B\\lesssim300$ MeV supports the idea that baryon-number fluctuations equilibrate locally at freeze-out.","At $\\mu_B\\gtrsim300$ MeV the data point to enhanced baryon-number fluctuations, up to about 50% above the lattice/HRG baseline at 7.7 GeV, which, if confirmed, would signal attractive correlations not present in the noncritical baseline.","The measured nonmonotonic peak in $\\hat{C}_3/\\hat{C}_1$ near 11 GeV is reproduced without irreducible three- or four-baryon correlations; it is generated by the energy dependence of $\\chi_2^B$ combined with exact baryon-number conservation.","Third- and fourth-order susceptibilities are too weakly constrained to discriminate critical-point scenarios; definite conclusions about genuine multi-baryon correlations will require reduced experimental uncertainties or new observables.","The two complementary extractions, from proton factorial cumulants and from net-proton cumulants, agree within uncertainties, providing an internal consistency check of the forward model."],"supporting_citations":[{"why":"Maximum-entropy freeze-out: maps local baryon susceptibilities to baryon-antibaryon factorial cumulants.","marker":"[31]"},{"why":"Extends the maximum-entropy mapping to hydrodynamic hypersurfaces, supplying the field-to-hadron fluctuation relations used here.","marker":"[32]"},{"why":"SAM-3.0 imposes exact global baryon-number conservation on the accepted subsample.","marker":"[34]"},{"why":"Defines the hydrodynamic excluded-volume baseline and the acceptance/thinning treatment the framework builds on.","marker":"[22]"},{"why":"Supplies the low-energy beam-energy-scan (net-)proton cumulant measurements used in the fits.","marker":"[20]"},{"why":"Supplies the higher-energy proton fluctuation data included in the extraction.","marker":"[38]"},{"why":"The finite-density lattice-based equation of state used as the primary comparison target.","marker":"[12]"},{"why":"A second lattice-based comparison band for the extracted second-order susceptibility.","marker":"[49]"},{"why":"The chemical freeze-out line used to map collision energies to $(T,\\mu_B)$ points.","marker":"[43]"},{"why":"The hadron list used to set the absolute normalization $\\bar\\chi_2^B$.","marker":"[6]"}],"fun_headline_variants":["Baryon susceptibility measured at freeze-out from heavy-ion data","First Bayesian extraction of baryon susceptibility from proton cumulants","Baryon susceptibility agrees with lattice up to 300 MeV","No triple-baryon correlations needed for RHIC fluctuation peak"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the low-energy rise in the extracted fluctuation strength is a property of equilibrium QCD matter, not an artifact of event-to-event volume or baryon-stopping fluctuations that the model leaves out; the paper itself flags this as the main caveat.","fun_headline_variants_meta":{"raw":{"variants":["Baryon susceptibility measured at freeze-out from heavy-ion data","First Bayesian extraction of baryon susceptibility from proton cumulants","Baryon susceptibility agrees with lattice up to 300 MeV","No triple-baryon correlations needed for RHIC fluctuation peak"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00092,"raw_usage":{"total_tokens":4046,"prompt_tokens":1146,"completion_tokens":2900,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":762,"completion_tokens_details":{"reasoning_tokens":2830}},"tokens_in":762,"tokens_out":2900,"duration_ms":21813,"temperature":1.0,"reasoning_tokens":2830,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:11:11.980930+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If centrality-bin-width or more complete volume-fluctuation corrections revise the measured cumulants so that the extracted $\\chi_2^B/\\bar\\chi_2^B$ at the lowest collider energy drops back into the lattice band, or if the nonmonotonic peak in $\\hat{C}_3/\\hat{C}_1$ vanishes when the same framework is applied to fixed-target data, then the central claim of a data-driven low-energy enhancement and of a conservation-driven peak would be falsified.","supporting_citations":[{"cited_title":"Subensemble Acceptance Method 3.0: General Corrections to Cumulants from Exact Conservation Constraints","cited_arxiv_id":"2607.01783","evidence_quote":"SAM-3.0 imposes exact global baryon-number conservation on the accepted subsample."}],"review_version":1}