{"id":"05996d2d-49cd-4279-8049-cca13f9499bc","arxiv_id":"2607.04191","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"An expanding spherical fireball model with Cooper–Frye freeze-out describes STAR pT spectra of π±, K±, p, p̄ at √sNN=7.7–39 GeV and predicts Gaussian rapidity distributions.","lead":"A spherical expanding fireball with blast-wave flow fits STAR mid-rapidity pT spectra of light hadrons in Au+Au collisions at RHIC BES energies (7.7–39 GeV) using shared freeze-out parameters. It supplies a low-parameter alternative to full hydrodynamics for spectra and Gaussian rapidity shapes.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Spherical approximation for integrated spectra is least secure for non-central collisions, where initial ellipticity is large and may affect fit quality and parameter reliability.","rationale":"The reader correctly isolated the spherical-geometry assumption as the weakest link supporting the five-parameter claim. The manuscript itself flags the initial ellipticity and then asserts (without quantitative check) that integrated observables are insensitive; the elevated peripheral χ² and the existence of the authors’ own eight-parameter elliptic alternative make this the single most load-bearing vulnerability. No stronger internal inconsistency or data-description failure appears once that assumption is granted. Consequently the CONDITIONAL verdict and the rest of the reader’s assessment stand; no adjustment is required.","tokens_in":20751,"tokens_out":531,"duration_ms":33617,"concrete_test":"For the 50–60 % and 70–80 % bins at √sNN=19.6 GeV, recompute the azimuthally integrated mid-rapidity pT spectra with the expanding elliptic fire-cylinder of Ref. [59], fixing its volume and mean radial flow to the spherical values used here; if the resulting spectra differ from the spherical results by more than experimental uncertainties or raise χ² by >20 %, the insensitivity claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (good pT description + Gaussian rap predictions with five parameters) rests on the assertion in Sec. II (paragraph after Eq. 5) that azimuthally integrated pT spectra and rapidity distributions are largely insensitive to initial elliptic geometry and anisotropic flow. This lets the authors replace the true elliptic-cylinder overlap (Eqs. 2–4) by an effective spherical r0 (Eq. 5) and a purely radial blast-wave profile vr=(r/rB)ṙB (Eqs. 7–8). The assumption is weakest for mid-central to peripheral bins (40–80 %), where a0/b0 is large; the higher set-1 χ²/NDF values (up to 7.3) in those bins may signal residual geometry dependence rather than only hard-process contamination. If the integrated spectra do feel the anisotropy, the common (Tkin,A,v∞,tf) extracted from pions, the species-dependent µkin, and the untested Gaussian rapidity shapes all become unreliable, and the claimed five-parameter economy relative to the eight-parameter elliptic model of Ref. [59] collapses.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript models the medium in Au+Au collisions at RHIC BES energies (√sNN = 7.7–39 GeV) as an expanding spherical fireball whose radius evolves as rB(t) = r0 + v∞[t − (1 − e−At)/A], with a linear radial flow profile vr = (r/rB)ṙB. Particle spectra for π±, K±, p and p̄ are obtained via the Cooper–Frye formula on an instantaneous freeze-out hypersurface. Common kinetic freeze-out parameters (Tkin, A, v∞, tf) are fixed from mid-rapidity pion pT spectra; only a species-dependent µkin is then adjusted for the remaining hadrons. The model is shown to describe STAR pT spectra across six centralities and is used to predict approximately Gaussian rapidity distributions.","tokens_in":21109,"tokens_out":1155,"duration_ms":26061,"significance":"If the description holds, the work supplies a compact, five-parameter alternative to more elaborate blast-wave or hydrodynamic parametrizations for azimuthally integrated spectra and rapidity distributions in the BES regime. The self-consistent link between surface velocity and radius evolution, the systematic centrality trends of the extracted parameters, and the explicit comparison of two fitting windows (full-range versus STAR-like restricted intervals) are useful for phenomenological surveys and as a baseline for more differential studies. The reduction relative to the eight-parameter elliptic fire-cylinder of Ref. [59] is a concrete practical advantage when only integrated observables are required.","major_comments":[{"comment":"Sec. II (paragraph following Eq. (5) and Eqs. (1)–(8)): The central claim of five-parameter economy rests on the assertion that azimuthally integrated pT spectra and rapidity distributions are largely insensitive to the initial elliptic geometry. This is least secure for the 40–80 % bins, where a0/b0 is large and Set-1 χ²/NDF reaches 4–7.3 (Table I). A quantitative estimate of residual anisotropy effects, or a direct side-by-side comparison with the elliptic model of Ref. [59] for at least one mid-central and one peripheral bin, is needed to substantiate that the spherical reduction does not degrade the description or bias the common flow parameters.","section":"Sec. II, Eqs. (1)–(8), Table I"},{"comment":"Sec. III and Figs. 5–7: The rapidity distributions are presented as model predictions, yet only the mid-rapidity STAR points are shown. Because dN/dy|y=0 is simply the pT integral of the already-fitted spectrum, agreement at that single point is largely by construction once the pT description is acceptable. The Gaussian shape itself therefore remains untested; either a comparison to existing full-rapidity data (SPS/GSI or RHIC) or a clear statement that no such data exist for these energies/centralities is required to support the predictive claim.","section":"Sec. III, Figs. 5–7"},{"comment":"Table I and the accompanying discussion of fitting ranges: The main rapidity results are generated with Parameter Set-1 (full pT range), whose χ²/NDF values are systematically higher than both Set-2 and the STAR blast-wave fits. While the text attributes the excess to resonance and hard-process contributions, it is not demonstrated that the extracted (A, v∞, tf) remain stable enough under this contamination to justify their use for the unmeasured rapidity shapes. A short robustness check (e.g., propagating Set-2 parameters into the rapidity distributions) would clarify which set underpins the central claims.","section":"Table I, Sec. III"}],"minor_comments":[{"comment":"Fig. 1 caption and the paragraph defining angles: the dual conventions for ϕ (position space in the reaction plane versus momentum-space ϕp in the XY plane) are easy to misread; a single clarifying sentence or a small inset would help.","section":"Fig. 1, Sec. II"},{"comment":"Table II: several µ entries for peripheral bins at low energy appear with large magnitude and opposite sign for particles and antiparticles; a brief remark on whether these remain within the expected range of effective chemical potentials (after feed-down absorption) would be useful.","section":"Table II"},{"comment":"Throughout: the notation switches between Tkin/µkin and T/µ; consistent subscripts would improve readability.","section":null},{"comment":"References: the recent spherical/spheroidal works [29–32] are cited, but a one-sentence contrast with the present radial-velocity prescription would better locate the novelty.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid but incremental phenomenology paper. Its main selling point—the spherical reduction relative to the authors’ own elliptic model [59]—needs the quantitative checks requested above before the “fewer free parameters” claim can be regarded as demonstrated. Scope is appropriate for a nuclear-theory journal that publishes blast-wave and fireball studies; I see no citation or novelty-disclosure issues."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This paper takes a spherical expanding fireball with a simple radius-evolution law (r_B(t) driven by A and v_∞) and applies Cooper–Frye freeze-out to the full RHIC BES Au+Au data set for π, K, p. What is actually new is the systematic extraction of the shared flow parameters from mid-rapidity pion p_T spectra across six centralities and five energies, then the reuse of those parameters (plus one µ_kin per species) for the other hadrons, plus the resulting Gaussian rapidity predictions and the two tables of numbers. That is a legitimate, transparent extension of the Siemens–Rasmussen / Bondorf line and of the recent few-GeV spherical work they cite.\n\nThey do the technical job carefully. The master integral (Eq. 10) is standard, the initial r_0 is obtained from the usual elliptic-cylinder volume, the flow profile is linear in r, and the figures show that the p_T spectra are reproduced once the pion parameters are fixed. Mid-rapidity yields come out right by construction, and the centrality trends (longer t_f, cooler T_kin, faster flow development in central collisions) make physical sense. Citation pattern is honest; they flag the eight-parameter elliptic alternative and the resonance-feed-down issue.\n\nSoft spots are real but proportionate. Set-1 χ²/NDF reaches 5–7 in the more peripheral bins when the full p_T range is used; restricting the range (Set-2) cleans it up, as expected. µ_kin is an effective parameter that swallows feed-down, which they acknowledge. The rapidity distributions are only confronted with the single mid-rapidity point that STAR published, so the Gaussian shape remains a prediction. The spherical approximation itself is weakest for 40–80 % centrality where the initial a_0/b_0 ratio is large; the higher χ² there may partly reflect residual geometry dependence rather than pure hard-process contamination. Still, for azimuthally integrated spectra the assumption is the usual one and does not collapse the five-parameter economy they claim relative to the elliptic model.\n\nMath and data handling look solid; nothing is incoherent. This is useful for anyone who needs a quick, analytic freeze-out background at BES energies or who wants tabulated kinetic parameters under a controlled spherical ansatz. I would send it to peer review; a referee can push on the peripheral χ² and the untested rapidity shapes, but the work is already clear enough to deserve that scrutiny. Worth engaging if you care about simple dynamical models or BES freeze-out systematics.","headline":"Clean, low-parameter spherical blast-wave fit to STAR BES spectra that delivers usable freeze-out tables, with the expected limitations of the geometry assumption and unrestricted χ².","tokens_in":21684,"tokens_out":637,"would_cite":true,"duration_ms":13745,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","25.75.Ld","24.10.Nz"],"model":"grok-4.5","headline":"A spherical expanding fireball with five parameters reproduces light-hadron pT spectra at RHIC BES energies and yields Gaussian rapidity distributions.","keywords":["spherical fireball","blast-wave","Cooper-Frye","kinetic freeze-out","RHIC Beam Energy Scan","pT spectra","rapidity distributions","radial flow"],"falsifier":"Measure the full rapidity distributions of identified light hadrons at the same BES energies and centralities; if they deviate systematically from the predicted Gaussians, or if the same five parameters fail to describe the pT spectra once resonance feed-down is removed, the model is ruled out.","tokens_in":21618,"feed_emoji":"🔥","tokens_out":820,"duration_ms":7366,"temperature":0.7,"pith_summary":"The paper shows that Au+Au collisions from 7.7 to 39 GeV can be described by a spherically expanding fireball whose radius grows according to a simple two-parameter law, with a blast-wave-like radial velocity profile fixed by that growth. Particle spectra are obtained from the Cooper–Frye formula on an instantaneous freeze-out surface. Kinetic freeze-out temperature, expansion parameters and freeze-out time are fixed once and for all by the mid-rapidity pion spectra; every other light hadron then requires only its own chemical potential. The resulting pT spectra match STAR data across centralities, while the same parameters automatically produce approximately Gaussian rapidity distributions. The practical payoff is a compact, hydrodynamics-inspired description that needs far fewer free parameters than full anisotropic blast-wave models yet still captures the bulk observables that do not depend on azimuthal anisotropies.","feed_headline":"Five-parameter spherical fireball fits RHIC light-hadron spectra","feed_subtitle":"Same parameters yield Gaussian rapidity distributions from 7.7 to 39 GeV across centralities","key_machinery":"The radius law rB(t)=r0+v∞[t-(1-e-At)/A] that sets both the surface velocity and the blast-wave-like radial rapidity profile vr=(r/rB)ṛB, inserted into the Cooper–Frye integral over a constant-time spherical freeze-out surface.","core_discovery":"An expanding spherical fireball whose surface velocity is identified with the time derivative of its radius, combined with a linear interior flow profile and Cooper–Frye freeze-out, simultaneously describes the measured mid-rapidity pT spectra of π±, K±, p and p-bar at RHIC Beam Energy Scan energies and predicts Gaussian-like rapidity distributions, using only five parameters fixed primarily by the pion spectra.","pith_inferences":["Because the model already works with only azimuthally integrated data, it supplies a cheap baseline against which the necessity of elliptic or viscous corrections can be judged.","The two fitting windows (full versus STAR-restricted pT ranges) give a direct handle on how much the extracted freeze-out parameters are contaminated by hard or resonance contributions.","The same radius law can be exported to lower-energy fixed-target experiments where full hydrodynamics is still expensive."],"forward_implications":["Kinetic freeze-out temperature falls and lifetime rises from peripheral to central collisions, giving a quantitative map of cooling and expansion.","Radial flow develops faster in more central events, consistent with stronger pressure gradients.","Effective chemical potentials absorb missing resonance feed-down and can be read off species by species.","The same parameter set can be reused for electromagnetic or heavy-flavor probes that need an analytic medium evolution."],"fun_headline_variants":["Five-parameter spherical fireball fits RHIC light-hadron pT spectra","Expanding fireball describes π K p spectra from 7.7–39 GeV Au+Au","Cooper-Frye fireball model yields Gaussian rapidity distributions at RHIC","Blast-wave spherical fireball unifies light hadron data across BES energies","Radial flow from fireball radius growth matches STAR midrapidity spectra"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"Azimuthally integrated spectra and rapidity distributions stay essentially unchanged when the real initial elliptic geometry and anisotropic flow are replaced by an effective spherical fireball.","fun_headline_variants_meta":{"raw":{"variants":["Five-parameter spherical fireball fits RHIC light-hadron pT spectra","Expanding fireball describes π K p spectra from 7.7–39 GeV Au+Au","Cooper-Frye fireball model yields Gaussian rapidity distributions at RHIC","Blast-wave spherical fireball unifies light hadron data across BES energies","Radial flow from fireball radius growth matches STAR midrapidity spectra"]},"model":"grok-4.5","effort":"low","cost_usd":0.004424,"raw_usage":{"total_tokens":1324,"prompt_tokens":784,"num_sources_used":0,"completion_tokens":108,"cost_in_usd_ticks":44240000,"prompt_tokens_details":{"text_tokens":784,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":432,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":784,"tokens_out":108,"duration_ms":4160,"temperature":1.0,"reasoning_tokens":432,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T21:03:05.975038+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Measure the full rapidity distributions of identified light hadrons at the same BES energies and centralities; if they deviate systematically from the predicted Gaussians, or if the same five parameters fail to describe the pT spectra once resonance feed-down is removed, the model is ruled out.","supporting_citations":[],"review_version":1}