{"id":"2ff41912-ccd7-4ccb-9bd2-916c17823a1f","arxiv_id":"2511.23267","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In simulations of a Perseus-like cluster with observed-level pre-existing turbulence, turbulent heating is insufficient to balance cooling in the core, and the XRISM central velocity rise is transient jet-driven bulk flow.","lead":"This paper simulates a Perseus-like galaxy cluster with both observed-level pre-existing turbulence and an AGN jet, and finds that turbulence dissipates too little heat to offset radiative cooling in the cluster core. The result pushes back against the idea that AGN-powered turbulence can solve the cooling-flow problem, and it reinterprets XRISM's central velocity broadening as transient bulk motion rather than sustained turbulence.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim rests on an undefined Q_cool: §3.5/Fig. 10 compare Q_turb to radiative cooling, but the manuscript never specifies the cooling function, metallicity, or radial bins.","rationale":"The reader's verdict is CONDITIONAL, and the reader noted in the rationale that 'the central comparison to radiative cooling lacks a defined baseline (Q_cool is never specified)', but the formal 'weakest_assumption' field points to the Kolmogorov-cascade assumption. My read overlaps with the rationale but elevates the missing Q_cool definition to the primary load-bearing concern. The paper's abstract and conclusions depend directly on Fig. 10, which plots Q_turb against an undefined Q_cool. Because radiative cooling is excluded from the simulations, the comparison requires an externally supplied cooling baseline; without specifying it, the central result is untestable. This is an internal omission, not a disagreement with consensus, so it is a correctness/reproducibility risk. Even if the Kolmogorov constants and inertial-range selection are correct, the conclusion only follows if Q_cool is known. The concern does not move the verdict away from CONDITIONAL — the paper may be correct, but the quantitative claim cannot be assessed or replicated until Q_cool is specified.","tokens_in":17210,"tokens_out":5003,"duration_ms":51971,"concrete_test":"Request from the authors the exact definition of Q_cool used in Fig. 10: the cooling function (e.g., Gnat & Ferland 2012, Raymond-Smith), metallicity, gas density/temperature profiles, and shell radii. Independently recompute the five shell cooling rates using a standard ICM cooling function applied to the same temperature/density data. If the recomputed Q_cool differs from the plotted curves by more than the shaded uncertainty in the inner shells (r < 50 kpc), the central shortfall claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim — that turbulent heating cannot offset radiative cooling, especially in the core — is established by comparing Q_turb estimates from Eqs. (9)/(10) with 'radiative cooling rates' in §3.5/Fig. 10. However, Q_cool is never defined. The simulations explicitly exclude radiative cooling (Eqs. 1–3), so Q_cool must be an external input, yet no cooling function, metallicity, density/temperature profile, or clustercentric radial bin edges are given, and no equation or reference supports the plotted cooling curves. Without this baseline, the central inequality Q_turb < Q_cool cannot be independently verified, even if the Kolmogorov-based ε estimates are accurate. The reader's weakest_assumption focuses on cascade constants; while relevant, it is secondary. The comparison itself is unconstrained: a plausible standard cooling function might or might not reproduce the plot, but as written the quantitative result is not reproducible.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The authors present 3D hydrodynamic simulations of a Perseus-like galaxy cluster using the FLASH code, comparing three controlled runs: pre-existing turbulence only, AGN jet only, and both. The turbulence is stirred by spectral Ornstein–Uhlenbeck forcing calibrated to the Hitomi line-of-sight velocity dispersion, and the AGN is a single 10 Myr kinetic bipolar jet. The central claim is that when the simulated velocity field is analyzed with the second-order velocity structure function and the kinetic-energy power spectrum under Kolmogorov scalings, the inferred turbulent heating rate Q_turb is smaller than the radiative cooling rate Q_cool in the cluster core, so turbulent heating alone cannot offset radiative cooling in cool-core clusters. A secondary claim is that the central σ_LOS enhancement recently reported by XRISM is reproduced by the combined run but is better interpreted as transient, coherent bulk motion than as sustained turbulence.","tokens_in":17487,"tokens_out":6799,"duration_ms":75368,"significance":"If the central comparison is made reproducible, the paper would strengthen the simulation-based case that turbulent heating is subdominant in cool-core cluster cores, while explicitly including pre-existing turbulence stirred to the observationally constrained Perseus level. Strengths of the paper are the controlled run design, the use of two independent estimators of the dissipation rate, the explicit interpretation of the estimates as upper limits, the convergence check, and the direct engagement with the new XRISM velocity-dispersion measurement. The main quantitative result, however, currently rests on an undefined and unreferenced radiative-cooling baseline, and the Kolmogorov-based estimates would benefit from a quantitative justification of the underlying assumptions. These issues are local and fixable; they do not by themselves invalidate the qualitative conclusion, but they block verification as written.","major_comments":[{"comment":"The central inequality Q_turb < Q_cool is not testable as written because Q_cool is never defined. Radiative cooling is excluded from the simulations (Eqs. 1–3), so Q_cool must be an external input, but the manuscript supplies no cooling function, metallicity, redshift, density/temperature profile, radial shell edges, or reference. Figure 10 also lacks explicit axis labels and units. Please specify how Q_cool was computed or from which observational estimate it was taken, and give the numerical values in the five radial shells. Without this, the abstract's claim cannot be verified or reproduced.","section":"§3.5 / Fig. 10"},{"comment":"The dissipation-rate estimates assume Kolmogorov incompressible turbulence with fixed constants (C2 ≈ 2.0, CK ≈ 1.65) applied over a 20–100 kpc inertial range. The justification that this is appropriate—solenoidal dominance and Fr > 0.1—is qualitative. Because the absolute magnitude of ε carries the whole argument, please quantify: (i) the compressive kinetic-energy fraction in the inertial range; (ii) the sensitivity of ε to the chosen fitting interval and to constants within published ranges; (iii) the variation of D_LL/E(k) across radial shells if shell-resolved estimates are used. This matters most in the outer bins, where Q_turb and Q_cool become comparable and a factor-of-order-unity error could alter the conclusion.","section":"§2.4, §3.5"},{"comment":"The method used to assign Q_turb to the five clustercentric shells is underspecified. The text says Q_turb = ρε with ρ averaged within each shell, but does not state whether ε is evaluated separately per shell or taken from a global VSF/E(k), nor the shell boundaries or weighting (emissivity, volume). The claim that the two methods bracket the true heating rate also requires that the same radial binning and fitting procedure be applied to both. Please specify the radial bins and the shell-wise computation, and give the resulting ε values (or Q_turb values) in a table so the comparison can be reproduced.","section":"§3.5 / Fig. 10"},{"comment":"The secondary conclusion that the central σ_LOS enhancement is \"not turbulence\" is based on the visual similarity of E(k) at t = 30 Myr and the short-lived bump in the combined run. This is not quantitative enough to support the strong wording. Please provide a quantitative measure: the fraction of injected jet energy that appears as solenoidal vs irrotational power on scales below the driving range, the decay time of excess large-scale power relative to the local eddy turnover time, or a comparison of the shell-resolved VSF slopes in the core with and without the jet. As written, the statement that XRISM's central enhancement is primarily coherent bulk motion is suggestive but not established.","section":"§4.2 / Fig. 11"}],"minor_comments":[{"comment":"The abstract and conclusions should soften \"CC clusters\" to \"the Perseus-like cluster modeled here\"; the simulations contain one cluster model with a single episodic jet and no radiative cooling.","section":"Abstract / §5"},{"comment":"Please state the units of Q_turb and Q_cool and explicitly define the radial bins. The text refers to \"radiative cooling rates\" but does not say whether these are per unit volume, per unit mass, or integrated luminosities in shells.","section":"Fig. 10 / §3.5"},{"comment":"The σ_LOS calibration is described as domain-averaged, while Hitomi/XRISM measure inner-core values. Clarify the region used for matching and verify that the simulated central σ_LOS profile is consistent with the observational aperture.","section":"§2.2"},{"comment":"The Froude number is quoted as Fr > 0.1 but never defined. Give the formula and the calculated value (with the integral-scale quantities used) so the stratification argument can be checked.","section":"§4.1"},{"comment":"The data availability statement says data will be shared upon reasonable request. Consider depositing analysis scripts and derived profiles in a permanent repository with a DOI; this would make the central figures auditable.","section":"General / Data availability"},{"comment":"Typographical and formatting issues: the Table 1 header appears as \"T able 1\"; \"velocity disperson\" in §4.1; a repeated sentence after Fig. 11; inconsistent use of spaces in \"T urb+Jet\" in captions and text.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The blocking issue is the undefined Q_cool baseline in §3.5; this is fixable and should be addressed before the paper can be accepted. The other major comments are also within the scope of a revision. I see no evidence of circularity or of the authors fitting parameters to force Q_turb < Q_cool; the numerical experiment itself is reasonable and timely. If the authors can supply a proper, referenced cooling baseline and make the shell-by-shell numbers available, the paper could be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nKnow this paper as a solid, controlled simulation study that adds a genuinely new ingredient — pre-existing turbulence stirred to the observed Perseus level — to the AGN-feedback simulations that previously found turbulent heating subdominant. The authors show that even with this strong background turbulence, a single jet episode heats the gas via shocks and bubbles, not by feeding a sustained cascade. That's a clean result, and it's supported by two independent estimators of the dissipation rate, a convergence check, and an explicit statement that the heating rates are upper limits. The reinterpretation of the XRISM central sigma_LOS enhancement as transient bulk flow rather than sustained turbulence is plausible and worth taking seriously.\n\nThe soft spot is real and it's in the paper's central quantitative claim. In Section 3.5 and Figure 10 they compare Q_turb to \"radiative cooling rates,\" but Q_cool is never defined. The simulations exclude cooling entirely, so this is an external input, yet no cooling function, metallicity, or radial bin edges are provided. A reader cannot reproduce the plotted cooling curve or verify the inequality. That's a load-bearing omission for a paper whose abstract says turbulent heating cannot offset radiative cooling. It's fixable — one paragraph and a reference — but as written the quantitative result is not reproducible.\n\nThe tangent-point method for choosing the inertial-range scale is also a bit ad hoc. They pick the point tangent to the Kolmogorov power law and call the result an upper limit. That's defensible, but the rationale is thin, and the Kolmogorov constants carry their own uncertainty in a compressible, stratified flow. The reader's concern there is secondary, though; the bigger issue is the missing baseline.\n\nCredit where due: the paper is honest about its limitations, the controlled comparison across the three runs is methodologically clean, and it engages fairly with the observational literature. The convergence check is a plus. Nothing about the citation pattern or the setup sets off alarms.\n\nWho gets value: anyone working on cluster feedback or interpreting XRISM velocity data. The XRISM comparison alone makes it worth a serious referee. I'd accept it for peer review and ask for the Q_cool definition and a more transparent description of the tangent selection. A revision should handle that; no need to reject on these grounds.","headline":"A solid controlled simulation study that makes a real addition to the turbulent-heating debate, but the central Q_turb < Q_cool comparison is missing its baseline because Q_cool is never defined.","tokens_in":17926,"tokens_out":1898,"would_cite":true,"duration_ms":21276,"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":"This paper argues that turbulent heating alone cannot balance radiative cooling in cool-core clusters, even when pre-existing turbulence is included at observed levels, and that the XRISM central velocity rise is transient bulk motion.","keywords":["cool-core clusters","AGN feedback","intracluster medium turbulence","turbulent heating","Perseus cluster","velocity structure function","energy power spectrum","XRISM"],"falsifier":"Measure the two-point velocity statistics in the Perseus core over scales 20–100 kpc with a future X-ray spectrometer that has both high spectral resolution and sufficient spatial resolution (or use a large sample of emission-line velocity maps) and check whether D_LL(ℓ) follows ℓ^{2/3} with the assumed constant; if the measured dissipation rate, integrated over the core, equals or exceeds the radiative cooling rate, the paper's central claim fails. A complementary check: if the central σ_LOS enhancement observed by XRISM is resolved spatially and shows a Kolmogorov power spectrum extending to","tokens_in":17122,"feed_emoji":"🌌","tokens_out":5751,"duration_ms":52336,"temperature":0.7,"pith_summary":"Cool-core clusters radiate away their central gas energy faster than expected, yet the gas does not collapse into a cooling flow; something must heat it. A leading candidate has been turbulence, stirred either by AGN jets or by pre-existing motions, and some X-ray brightness fluctuation analyses have suggested turbulence alone could offset the cooling. This paper tests that idea in three-dimensional hydrodynamic simulations of a Perseus-like cluster, including a strong pre-existing turbulent field calibrated to the observed ~185 km/s line-of-sight velocity dispersion, plus a single AGN jet outburst. Using two independent Kolmogorov-based estimators—the velocity structure function and the kinetic energy power spectrum—the authors find that the turbulent heating rate falls below the radiative cooling rate in the cluster core (r < 50 kpc), implying turbulence alone cannot prevent catastrophic cooling. They further show that the jet's contribution to the velocity field is a short-lived coherent flow, not a self-sustaining cascade, which they argue explains the central line-broadening enhancement seen by XRISM without invoking AGN-powered turbulence.","feed_headline":"Turbulent heating cannot offset cooling in cluster cores","feed_subtitle":"Perseus-like simulations show a core heating deficit even with strong turbulence; the XRISM signal is likely transient bulk flow.","key_machinery":"The key machinery is the Kolmogorov-cascade estimator of turbulent dissipation. The velocity structure function D_LL(ℓ) gives the variance of longitudinal velocity differences and is fitted to C_2 ε^{2/3} ℓ^{2/3} in the inertial range; the energy power spectrum E(k) is fitted to C_K ε^{2/3} k^{-5/3}. The inertial range is identified at 20–100 kpc, ε is read off at the tangent points (ℓ_VSF ≈ 53.1 kpc, ℓ_Ek ≈ 52.6 kpc), and the volumetric heating rate Q_turb = ρ ε is compared with cooling in radial shells. A second element is the time-resolved E(k) comparison between the run with and without the jet, which separates a self-sustaining cascade from transient coherent large-scale motions.","core_discovery":"Central discovery: in a Perseus-like cool-core cluster with pre-existing turbulence stirred to the observed σ_LOS ≈ 185 km/s, two independent Kolmogorov estimators—the longitudinal velocity structure function (ℓ_VSF ≈ 53.1 kpc) and the energy power spectrum (ℓ_Ek ≈ 52.6 kpc)—give turbulent heating rates below the radiative cooling rate in the core (r < 50 kpc), and the rates are treated as upper limits. Early-time energy spectra show the AGN jet's injected energy decays within ~20–30 Myr without cascading, so the jet drives transient coherent bulk flows, not sustained turbulence; the authors argue the XRISM central σ_LOS enhancement in Perseus is this bulk motion, not AGN-powered turbulence.","pith_inferences":["If the central velocity enhancement in Perseus is indeed a transient bulk flow, then X-ray line measurements with spatial resolution sufficient to resolve 20–100 kpc scales might see a coherent shearing or dipole pattern in the velocity centroid map rather than random small-scale eddies; this is a testable distinction from sustained turbulence.","The simulation's conclusion depends on the assumption that numerical dissipation on the grid mimics the physical dissipation at the true ICM Reynolds number; if the real ICM has a much longer inertial range or the cascade is modified by magnetic fields, the inferred Q_turb could shift—an extension the authors note but do not quantify.","Because the simulation excludes radiative cooling, magnetic fields, thermal conduction, and cosmic rays, a natural next step is to repeat the comparison in self-regulated feedback runs with cooling; if cooling changes the density and stratification, the inertial-range scalings and the inferred heating rates may change, potentially moving the conclusion in either direction.","The paper's reinterpretation of XRISM data suggests a caution for the broader practice of inferring turbulent velocities from unresolved line broadening or surface brightness fluctuations: attributing all velocity variance to a turbulent cascade may systematically overestimate the dissipation available to heat the gas."],"forward_implications":["Turbulent heating is not the primary mechanism solving the cooling-flow problem in cool-core clusters; other AGN channels (bubble mixing, weak shocks, sound waves, cosmic-ray heating) must supply the missing heat.","Even when pre-existing turbulence is driven to the maximum level consistent with Hitomi/XRISM velocity measurements, the core heating rate remains below the cooling rate, so boosting turbulence alone will not close the energy budget.","The central σ_LOS enhancement observed by XRISM in Perseus should not be interpreted as direct evidence of AGN-powered turbulence; if the simulated picture is right, it is a transient coherent flow, and turbulent heating estimates based on that line broadening would be overestimates.","Both VSF and E(k) methods give the same radial trend (heating shortfall in the core, comparability at large radii), with the E(k) estimate systematically higher, bracketing the systematic uncertainty in the upper limit.","The jet's contribution to the global velocity statistics is minor, yet it still heats the ICM through pressure work, so AGN feedback remains effective even when the velocity field is turbulence-dominated."],"fun_headline_variants":["Turbulent heating falls short in cluster cores","Simulations: turbulence can't stop cooling in Perseus-like clusters","AGN jets create transient flows, not sustained turbulence","XRISM signal in Perseus likely bulk gas motion, not turbulence","Preexisting turbulence doesn't rescue cluster cores from cooling"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the 20–100 kpc velocity fluctuations in the cluster core obey a classical incompressible Kolmogorov cascade, so that the dissipation rate can be read off from the structure function and power spectrum with standard constants; if compressibility, stratification, or numerical viscosity distort the inertial-range scaling, the inferred turbulent heating rate could be off by a factor and the central shortfall might shrink or grow.","fun_headline_variants_meta":{"raw":{"variants":["Turbulent heating falls short in cluster cores","Simulations: turbulence can't stop cooling in Perseus-like clusters","AGN jets create transient flows, not sustained turbulence","XRISM signal in Perseus likely bulk gas motion, not turbulence","Preexisting turbulence doesn't rescue cluster cores from cooling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00046,"raw_usage":{"total_tokens":2139,"prompt_tokens":741,"completion_tokens":1398,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":1313}},"tokens_in":485,"tokens_out":1398,"duration_ms":12491,"temperature":1.0,"reasoning_tokens":1313,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T19:32:46.094732+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the two-point velocity statistics in the Perseus core over scales 20–100 kpc with a future X-ray spectrometer that has both high spectral resolution and sufficient spatial resolution (or use a large sample of emission-line velocity maps) and check whether D_LL(ℓ) follows ℓ^{2/3} with the assumed constant; if the measured dissipation rate, integrated over the core, equals or exceeds the radiative cooling rate, the paper's central claim fails. A complementary check: if the central σ_LOS enhancement observed by XRISM is resolved spatially and shows a Kolmogorov power spectrum extending to","supporting_citations":[],"review_version":1}