{"id":"0db0494b-b316-41e2-9649-d8b697168fef","arxiv_id":"2604.00302","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In AMPT and UrQMD simulations, fluid acceleration in heavy-ion collisions reaches a few hundred MeV, peaks at fireball boundaries, and evolves from stopping-dominated deceleration at low energies to boundary/tilt pulses at high energies.","lead":"This paper uses two transport models (AMPT and UrQMD) to map how fast pieces of the hot fireball created in heavy-ion collisions accelerate, in space and time. The acceleration can reach hundreds of MeV and peaks at the fireball edge, which matters for proposed Unruh-heating, chiral-transition, and spin-polarization effects.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Peak acceleration values are derivative quantities whose magnitude may be set by the Gaussian smearing width and energy-density threshold; the claimed insensitivity is not demonstrated.","rationale":"The reader flagged the smeared velocity field and early-time reliability as the weakest assumption. My concern is closely related but more specific: even at nominal later times, the derivative nature of acceleration makes the quantitative peak values sensitive to the smearing resolution and threshold, which are not systematically tested. This is load-bearing because the abstract's headline number and the boundary-localization result feed directly into the Unruh-effect discussion. The paper earns credit for using established transport models, multiple energies, and honest caveats; the qualitative directions (outward transverse acceleration, stopping-induced longitudinal deceleration) are plausibly robust across models. However, the quantitative range 'several hundred MeV' should be conditional on a resolution-convergence check. Since the reader's verdict is already CONDITIONAL and explicitly requests sensitivity checks, my assessment does not change the verdict.","tokens_in":15676,"tokens_out":6417,"duration_ms":76334,"concrete_test":"Recompute the b=0, √s=200 GeV AMPT volume-averaged proper acceleration ⟨a⟩(t) for a 3×3 grid of smearing widths σ∈{0.3,0.6,1.2} fm and thresholds ε_c∈{20,50,100} MeV/fm³, all other settings fixed; if the peak ⟨a⟩ shifts by more than 50% across the grid, the reported few-hundred-MeV values are resolution artifacts and should be reported only as range estimates with sensitivity caveats. Optionally, check consistency with Eq. (2.4) using an assumed EoS.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that proper acceleration peaks at several hundred MeV and is boundary-enhanced. These numbers are obtained as derivatives of a velocity field extracted from a Gaussian-smeared energy-momentum tensor (Eq. 3.6, σ⊥=σz=0.6 fm). The gradient scale of any such smeared field is controlled by σ: larger σ suppresses gradients, smaller σ amplifies particle-noise fluctuations. The paper states 'varying these width parameters by up to 50% leads only to minor changes' but shows no quantitative sensitivity scan. Similarly, the fireball boundary and the volume average are defined by an energy-density threshold ε_c=50 MeV/fm³ (Sec. IV). Since the largest accelerations live at that boundary, the threshold choice directly controls which regions contribute to the peak and average values. A shift in ε_c (or in σ) could move the quoted 100–500 MeV range substantially. The early-time caveat (Sec. IV C, t<0.5 fm/c) does not address this: the 200 GeV acceleration pulse peaks at t≈0.5–1.0 fm/c, just outside the shaded region, and its amplitude is still a smearing-scale-dependent derivative. If the 'several hundred MeV' are not stable under σ and ε_c variation, the Unruh-temperature and spin-polarization implications lose their quantitative foundation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses AMPT and UrQMD transport models, with a Gaussian smearing procedure (Eq. 3.6) to define the local energy-momentum tensor and fluid 4-velocity, and then computes the proper acceleration field of the produced medium in heavy-ion collisions. It reports spatial maps of the acceleration components, their impact-parameter dependence, and their time evolution across beam energies from 3.5 GeV to 2.76 TeV. The central claims are that the peak proper acceleration reaches a few hundred MeV, that transverse acceleration is outward and strongest at the fireball boundary, that longitudinal acceleration is strongly collision-energy dependent (nuclear stopping at low energies, passing-nucleus drag pulses at high energies), and that the volume-averaged acceleration is nearly centrality independent. The authors connect these results to possible Unruh temperatures, chiral restoration, and spin-polarization effects beyond vorticity.","tokens_in":15972,"tokens_out":3179,"duration_ms":39416,"significance":"If the reported acceleration fields are robust, the paper provides a useful first systematic map of a quantity that is usually only discussed qualitatively: the local 4-acceleration of the QGP fireball. The use of two independent transport models over a wide energy range is a genuine strength, as is the explicit recognition, in Sec. IV C, that very early times are model-dependent. The qualitative patterns—outward transverse acceleration, boundary enhancement via steeper pressure gradients and lower enthalpy density, and the longitudinal stopping-versus-drag crossover—are plausible and physically motivated. The quantitative claims, however, are derivative quantities whose magnitude may be controlled by the smearing width and the energy-density threshold. The paper's strongest quantitative conclusion, 'several hundred MeV,' is not yet supported without a systematic sensitivity study and without statistical uncertainties on the event-averaged results. The paper does not provide machine-checkable proofs or code, but it is transparent about the extraction method; the main weakness is the absence of demonstrated robustness of the central numbers.","major_comments":[{"comment":"The peak acceleration values are spatial derivatives of a velocity field obtained from Gaussian smearing with σ⊥=σz=0.6 fm. The paper states that varying these widths by up to 50% leads to only minor changes, but no quantitative sensitivity scan is provided. This is load-bearing because the gradient scale of the smeared field is set by σ; larger σ suppresses gradients and smaller σ amplifies particle-noise fluctuations. The authors should show ⟨a⟩, peak a, and the boundary emphasis for at least σ=0.3, 0.6, 0.9 fm. Without this, the quoted hundreds-of-MeV magnitudes and the boundary-enhancement claim are not established.","section":"Sec. IV A and Eq. (3.6)"},{"comment":"All volume averages and the definition of the fireball boundary use a fixed energy-density threshold εc=50 MeV/fm³. Since the largest accelerations live at the boundary, the threshold choice directly controls which regions enter the peak and volume-averaged values. A scan over εc (e.g., 25, 50, 100 MeV/fm³) is necessary to show that the centrality independence and the 'several hundred MeV' numbers are not threshold artifacts. The early-time caveat in Sec. IV C does not address this, because the high-energy pulse peaks at t≈0.5–1.0 fm/c, just outside the shaded unreliable region.","section":"Sec. IV A and Sec. IV C (εc=50 MeV/fm³)"},{"comment":"No statistical uncertainties are reported on any averaged curve. For UrQMD only 200 events per setting are used, and the acceleration is sensitive to rare early-time fluctuations. The finite-event-number effect should be quantified, either as error bars on ⟨a⟩ or as an event-by-event spread. This is particularly important for the small ⟨a_z⟩ signal in central collisions at 7.7 GeV, which the authors themselves suspect may be a numerical artifact (Sec. IV C).","section":"Figs. 2–4 and Sec. IV"},{"comment":"The impact-parameter comparison uses different time slices for low energies (t at maximum energy density) and high energies (fixed t=2 fm/c), and the paper explicitly says absolute magnitudes cannot be directly compared. Yet the text then uses this figure to conclude that proper acceleration increases with collision energy within each regime and that the volume average is nearly centrality independent. The claim of a smooth energy dependence is therefore not directly supported by the presented comparison. A consistent time selection, or a clear statement that only within-regime trends are being inferred, is needed.","section":"Sec. IV B, Fig. 2"},{"comment":"The boundary enhancement is interpreted using the ideal-fluid Euler equation a^μ=∇^μP/(ε+P), but the acceleration is not computed from a separately constructed pressure and enthalpy; it is computed from velocity gradients. The interpretation is plausible, but it is not a verification. If the authors wish to claim that pressure gradients and low enthalpy density are the driving mechanism, they should demonstrate this directly by computing ∇P and ε+P from the same smeared T^μν, or at least check the Euler relation approximately for the late-time, near-equilibrium stage.","section":"Sec. II, Eqs. (2.4) and Sec. IV A"}],"minor_comments":[{"comment":"There are several typos and stylistic issues: 'borotropic' should be 'barotropic,' 'inivertablly' should be 'inevitably,' 'anologue' should be 'analogue,' and 'The is a relativistic effect' needs rewording. Equation (2.10) writes q^μ_C = -κ(T a^μ - ∇^μT), which conflicts with the usual sign convention; please check the sign convention consistency.","section":"Sec. II"},{"comment":"The phrase 'the same applies to the figures below' is informal and should be replaced with an explicit statement in each caption about the collision system and energy labels.","section":"Fig. 2 caption"},{"comment":"The scaling factor K is said to be fit to multiplicities, but since it multiplies all components of T^μν uniformly, it does not affect the eigenvector velocity u^μ. The paper should state this explicitly to avoid any impression that the K fits influence the acceleration results beyond the εc cut.","section":"Sec. III"},{"comment":"The shaded region t<0.5 fm/c is useful, but some curves for low energies peak near or inside this region. It would help to report the peak time for each energy and clarify how much of the 'up to ~500 MeV' claim rests on the shaded, less reliable interval.","section":"Sec. IV C, Fig. 3"},{"comment":"The discussion of Unruh temperatures and their implications in Sec. V is somewhat speculative relative to the numerical results. A quantitative estimate of T_U from the reported acceleration values and a comparison with the thermodynamic temperature at the same spacetime points would strengthen the connection.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"This is a phenomenological transport-model study in an area that is currently active. The qualitative findings are likely correct and interesting, but the quantitative headline numbers ('several hundred MeV', 'centrality independence') are not yet supported because they depend on the smearing width, the energy-density threshold, and the statistical precision of the event sample. A sensitivity analysis and error bars are essential before the Unruh-temperature and spin-polarization implications can be taken seriously. I do not see a load-bearing logical inconsistency or circularity; the issues are about robustness and quantification, not about the method being invalid. The paper fits the journal if these revisions are made."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, systematic map of fluid 4-acceleration from AMPT and UrQMD, with honest caveats, but the headline peak values are derivative quantities and the paper does not demonstrate they are stable against the smearing width and density threshold. The qualitative patterns are probably right; the quantitative numbers are not yet backed up.\n\nWhat's new: it is the first AMPT/UrQMD survey of acceleration across 3.5 GeV to 2.76 TeV, with component-resolved spatial, temporal and centrality dependence. The velocity extraction follows the established Gaussian smearing of Refs. [40,47]; the physics questions follow the PHSD work Ref. [70], but the energy range and two-model comparison are new. The transverse outward acceleration, boundary enhancement, and energy-dependent longitudinal structure are plausible and supported by both models. Using the Euler equation to interpret the boundary enhancement is legitimate, not circular.\n\nSoft spots, in rough order of severity:\n\n- The magnitude of a is a gradient of a smeared velocity field. The paper says varying sigma by up to 50% \"leads only to minor changes\" but shows no scan. The boundary is defined by epsilon_c = 50 MeV/fm^3, and the largest accelerations live at that boundary; no threshold scan either. The stress-test concern is not refuted. For the 200 GeV case the peak sits at t ~ 0.5-1 fm/c, just outside the shaded unreliable region, so the caveat about early times does not cover the headline number. If the several-hundred-MeV values move when sigma or epsilon_c changes, the Unruh-temperature and spin-polarization implications lose their quantitative foundation.\n\n- No statistical errors on any averaged curve. 200 UrQMD events per setting is enough for qualitative shapes, not for claiming a 500 MeV peak to three significant figures.\n\n- Low- and high-energy panels are compared at different time slices (t(<epsilon>_max) vs t = 2 fm/c). The paper says so honestly, but it makes the energy dependence claim weaker than it looks.\n\n- No code or data released; reproducibility is limited.\n\nWhat holds: the qualitative story — boundary enhancement, outward transverse acceleration, stopping vs drag — is consistent across two models and is physically reasonable. The paper is honest about early-time model dependence. This is not a circular or fabricated result.\n\nWho should read it: heavy-ion phenomenologists working on acceleration, Unruh effect, chiral transition, or spin polarization. It deserves a serious referee. My recommendation: send it out, and require a quantitative sensitivity study for sigma and epsilon_c, plus error bars on the averaged curves, before the numbers are used.","headline":"Useful systematic map of fluid acceleration from two transport models, but the headline peak values are derivative quantities whose stability against smearing width and density threshold is not demonstrated.","tokens_in":16490,"tokens_out":2688,"would_cite":true,"duration_ms":31049,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q"],"model":"deepseek-v4-flash","headline":"The proper acceleration of quark–gluon plasma in heavy-ion collisions peaks at a few hundred MeV, is strongest at the fireball boundary, and can make Unruh temperatures exceed the local thermodynamic temperature early in the collision.","keywords":["fluid acceleration","heavy-ion collisions","quark-gluon plasma","relativistic hydrodynamics","Unruh effect","transport models","pressure gradients","spin polarization"],"falsifier":"Extract the same acceleration field from the same transport events with Gaussian widths doubled and halved (σ=1.2 fm and σ=0.3 fm); if the peak proper accelerations change by more than the model-to-model spread, the headline values are smearing artifacts rather than physical fields. A complementary check is to evolve the same initial profile with an ideal hydrodynamic code and see whether a^mu = ∇^mu P/(ε+P) reproduces the extracted boundary field.","tokens_in":15547,"feed_emoji":"⚛️","tokens_out":8407,"duration_ms":87551,"temperature":0.7,"pith_summary":"This paper sets out to map the fluid acceleration field inside the hot, dense matter created in heavy-ion collisions, a quantity studied far less than vorticity. Using transport simulations and Gaussian smearing to define a local fluid velocity, it finds that the proper acceleration reaches a few hundred MeV, with the strongest values at the fireball's peripheral boundary rather than in the bulk. The boundary enhancement is explained by the relativistic Euler equation: steep pressure gradients combine with small enthalpy density at the rim. The longitudinal component behaves very differently at low and high collision energies—nuclear stopping decelerates matter early at low energies, while passing nuclei drag the midrapidity plasma at high energies—and the volume-averaged acceleration is nearly independent of centrality. If these fields are physical, the associated Unruh temperature can exceed the thermodynamic temperature in early stages, affecting chiral symmetry restoration, deconfinement, and spin polarization.","feed_headline":"Peak fireball accelerations hit hundreds of MeV","feed_subtitle":"Boundary pressure spikes and low enthalpy make the QGP fastest at its rim; early Unruh temperatures may then exceed the critical QCD tempera","key_machinery":"The load-bearing identity is the relativistic Euler equation a^mu = (1/(ε+P)) ∇^mu P, which isolates acceleration as the pressure gradient divided by enthalpy density. The numerical machinery is Gaussian smearing of particle four-momenta from transport events into an energy-momentum tensor T^{mu nu}; the local fluid four-velocity is the timelike eigenvector of T^{mu nu}, and the 4-acceleration is u^nu ∂_nu u^mu, with proper acceleration a = sqrt(-a_mu a^mu). The boundary enhancement follows from the two-fold effect of a steep pressure drop and a small ε+P at the fireball rim.","core_discovery":"The paper claims that proper acceleration in heavy-ion fireballs reaches several hundred MeV and follows a^mu = ∇^mu P/(ε+P). Transverse acceleration points outward, peaking at the boundary where pressure gradients are steep and enthalpy density low; this persists at low energies and early times. Longitudinal acceleration is energy-dependent: nuclear stopping decelerates early at low energies, passing nuclei drag midrapidity plasma at high energies, and a dip-peak marks the crossover to transverse expansion. Volume-averaged acceleration is nearly centrality independent, as extremes localize at boundaries. These fields can make Unruh temperatures exceed the thermodynamic temperature early, af","pith_inferences":["Beyond the paper: because the peak acceleration sits at the boundary, any Unruh-induced signal should track the fireball surface-to-volume ratio across centralities; testing that scaling would separate acceleration effects from bulk thermodynamic ones.","Beyond the paper: comparing the extracted field with a full ideal-hydrodynamic evolution from the same initial profiles would show whether the boundary peaks reflect genuine pressure gradients or smearing artifacts.","Beyond the paper: the low-energy longitudinal deceleration from nuclear stopping should correlate with rapidity-dependent directed flow; measuring that correlation could give an independent handle on early-time acceleration magnitudes."],"forward_implications":["If the peak proper acceleration is indeed a few hundred MeV, the Unruh temperature a/(2π) in the early fireball can exceed the local thermodynamic temperature, making acceleration a control parameter for early-time QCD matter.","The boundary-localized structure implies that acceleration-driven effects—such as Unruh-like radiation or acceleration contributions to spin polarization—should be strongest near the fireball surface, not in the bulk.","The near-centrality independence of the volume-averaged acceleration would make bulk-averaged acceleration-driven observables roughly flat across centrality, distinct from energy-density-driven signals.","Strong longitudinal deceleration at low energy and drag pulses at high energy provide a direct kinematic readout of nuclear stopping and passing-nucleus dynamics, which future flow measurements could probe.","If acceleration contributes to spin polarization as in the mean-spin formula, the acceleration-induced p×a term would add to the vorticity term and could help explain local polarization measurements."],"fun_headline_variants":["Fireball acceleration peaks at hundreds of MeV","Boundary acceleration spikes in heavy-ion collisions","Unruh temperatures may exceed QCD critical in fireballs","Acceleration fields point outward, strongest at rim","Nuclear stopping shapes early-time acceleration"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The headline peak values depend on treating the Gaussian-smeared particle output of transport models as a genuine local fluid velocity whose gradients are physical accelerations, even during the earliest times when the matter is not thermalized and the fields carry large, possibly unphysical fluctuations.","fun_headline_variants_meta":{"raw":{"variants":["Fireball acceleration peaks at hundreds of MeV","Boundary acceleration spikes in heavy-ion collisions","Unruh temperatures may exceed QCD critical in fireballs","Acceleration fields point outward, strongest at rim","Nuclear stopping shapes early-time acceleration"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000486,"raw_usage":{"total_tokens":2197,"prompt_tokens":669,"completion_tokens":1528,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":413,"completion_tokens_details":{"reasoning_tokens":1459}},"tokens_in":413,"tokens_out":1528,"duration_ms":13432,"temperature":1.0,"reasoning_tokens":1459,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T05:33:18.552827+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Extract the same acceleration field from the same transport events with Gaussian widths doubled and halved (σ=1.2 fm and σ=0.3 fm); if the peak proper accelerations change by more than the model-to-model spread, the headline values are smearing artifacts rather than physical fields. A complementary check is to evolve the same initial profile with an ideal hydrodynamic code and see whether a^mu = ∇^mu P/(ε+P) reproduces the extracted boundary field.","supporting_citations":[],"review_version":1}