{"id":"4af1e14c-787e-4086-a664-25556e2f7450","arxiv_id":"2506.13422","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"AGN jet feedback is predicted to add only 4 to 6 percent non-thermal pressure in typical cluster cores, and the peak value can be used to infer the AGN jet duty cycle.","lead":"This paper predicts how much non-thermal pressure (turbulence-like pressure) AGN jets add to the hot gas in galaxy cluster cores, finding only a few percent. It then uses this prediction to estimate how often the central black hole's jets switch on and off in the Perseus cluster.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Perseus duty-cycle inference rests on an untested equivalence between the ensemble-mean and time-averaged NTP; if turbulence decays between outbursts, the single-epoch Hitomi value cannot be reliably mapped to a duty cycle.","rationale":"After reading the manuscript in good faith, the forward model is coherent and the comparison to observations is encouraging. The strongest part is the parameter-free kinetic-energy coupling (Section 3.2.3), which avoids tuning the NTP fraction. However, the headline claim—that the peak NTP fraction can be used to infer the AGN jet duty cycle of an individual cluster—hinges on the ensemble mean <F> being the right descriptor for a single-epoch observation of one cluster. The paper's own time-averaging (sqrt(delta) factor) is derived for the mean over the on/off cycle, but Hitomi provides an instantaneous snapshot. If turbulent NTP decays on timescales comparable to or shorter than the off-phase, the observed value is phase-dependent and the inferred duty cycle is meaningless. This is the same assumption the reader flagged as weakest. I agree with the reader's conditional verdict. A hydrodynamic test with repeated outbursts would settle the issue. No other concern—such as the post hoc epsilon cut or the equilibrium-time cut—is as damaging to the central claim, since those would merely shift the quantitative coefficients rather than invalidate the mapping from NTP to duty cycle.","tokens_in":44664,"tokens_out":13544,"duration_ms":131602,"concrete_test":"Run a 3D hydrodynamic simulation of a Perseus-like cluster core with a sequence of jet outbursts whose powers and ages are drawn from the same distributions as in Section 2.2, with a prescribed duty cycle delta (e.g., 0.2), and measure the volume-weighted NTP fraction in the central ~10 kpc as a function of time over at least 10 outburst cycles. Compare the time-averaged NTP fraction with the model's prediction from Eq. (62) using the Perseus-like coefficients in Table 4. If the time-average differs by more than ~30%, or if the instantaneous value at a random epoch deviates from the time-average by more than a factor of 2, the ensemble-to-time-average assumption is violated and the Perseus duty-cycle inference is not secure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central inference in Section 5.2 — mapping Hitomi's single-epoch NTP measurement to Perseus' AGN jet duty cycle via the relation in Eq. (62) — requires that the observed NTP fraction equals the time-averaged NTP over the AGN lifecycle. The model computes <F> as an ensemble average over a Monte Carlo population of outbursts (Eq. 57), and then assumes this equals the time average for a single cluster. This ergodicity assumption is stated rather than tested in Section 5.2.1 ('Perseus is known to have experienced frequent outbursts'). However, the model's own time-averaging in Eq. (41) (the sqrt(delta) factor) presupposes that NTP responds to the square of the instantaneous heating rate and that the time-average is the relevant observable. A snapshot observation like Hitomi is not a time average: if the turbulent velocity dispersion decays on a timescale shorter than the off-phase t_off = (1-delta)/delta * t_on, the observed NTP at a random epoch will be systematically below the time-averaged value unless the cluster is caught during an outburst. For Perseus, the ripple period is t_on+t_off ~ 9.2 Myr and the bubble inflation time ~2.5 Myr, so t_off ~ 6.7 Myr; the sound crossing time of the 3 kpc core is ~3 Myr. The paper cites Bîrzan et al. (2012) for long-lived bubble heating, but that concerns enthalpy injection from bubbles, not the decay of turbulence. No decay timescale is modelled or measured. If the NTP decays between outbursts, the inferred duty cycle of ~13% (or 3-48%) is biased low, and the central claim that NTP can be used to infer the duty cycle of a single cluster is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an analytical model for the non-thermal pressure (NTP) fraction induced by AGN jet outbursts in the cores of galaxy clusters. The authors sample jet powers and ages from observationally motivated probability distributions, evolve each source with a two-phase (ballistic and flaring) FR-I jet model, and couple the injected energy to spherical gas shells using the first law of thermodynamics. The model predicts mean radial NTP fraction profiles for different cluster environments and AGN jet duty cycles, finding peak values of roughly 4-6% for duty cycles of 10-30% in a typical cluster. The central new claim is a relation between the peak mean NTP fraction and the AGN jet duty cycle, which the authors apply to the Hitomi measurements in the Perseus cluster to infer a duty cycle of about 13%, with a broader range of 3-48% from the 2018 Hitomi data; this is claimed to be consistent with independent estimates of Perseus' recent jet activity. The paper also combines the core predictions with the large-scale NTP profiles from previous work, and presents a perturbation analysis that suggests the pristine-equilibrium assumption is self-consistent.","tokens_in":45149,"tokens_out":7027,"duration_ms":68450,"significance":"If the central assumptions hold, the paper offers a novel and potentially practical framework: inferring past AGN jet duty cycles from single-epoch X-ray measurements of core gas motions. The first-law-based derivation in Section 3.2 is internally consistent, the perturbation scheme in Section 6.2.2 converges, the predicted NTP fractions are of the order seen by Hitomi and the first XRISM results, and the code is publicly available. However, the duty-cycle inference depends on two load-bearing modeling choices that are not adequately tested: the equivalence between an ensemble-averaged population mean and a single cluster's time-averaged NTP, and the identification of the thermodynamic shell expansion rate with the turbulent velocity that produces non-thermal pressure. Because the practical headline result (the Perseus duty-cycle estimate) rests on these choices, the paper needs substantial revision before its central claim can be regarded as robust.","major_comments":[{"comment":"The inference of Perseus's duty cycle rests on treating the Monte Carlo ensemble-averaged mean <F> (Eq. (57)) as equal to the time-averaged NTP of a single cluster over its AGN lifecycle. The sqrt(delta) factor introduced in Eq. (41) already assumes that NTP responds as the square of the instantaneous heating rate and that the time average is the relevant quantity. A Hitomi snapshot is not a time average: if the turbulent velocity dispersion decays on timescales shorter than the off-phase t_off = (1-delta)/delta * t_on, the observed NTP at a random epoch will be systematically below the time-averaged value unless the cluster happens to be caught in an outburst. For Perseus, the ripple period is about 9.2 Myr and the outburst time about 2.5 Myr, giving t_off about 6.7 Myr, comparable to the roughly 3 Myr sound-crossing time of the 3 kpc core. The paper cites Bîrzan et al. (2012) for long-lived bubble heating, but that work concerns enthalpy injection from bubbles rather than the decay of turbulence. No NTP decay timescale is modelled or observationally constrained. Please either model the time dependence explicitly and test the ergodicity assumption, or reframe the inferred duty cycle as an upper or lower limit that depends on the assumed NTP lifetime.","section":"Section 5.2.1 and Eq. (41)"},{"comment":"The identification of dr/dt as the gas velocity that generates non-thermal pressure is a modeling leap. The first law of thermodynamics, Eq. (43), determines the rate of change of shell volume under a given energy injection rate; the 'velocity kick' in Eq. (47) is defined as that radial expansion rate. The paper then uses v_gas^2 in Eqs. (55)-(57) as the kinetic energy density of gas motions contributing to NTP. In a real cluster, the injected energy can go into thermal heating, bulk expansion, turbulence, waves, and other forms, and there is no derivation given for why the adiabatic shell expansion rate equals (or maps to) the observable velocity dispersion producing non-thermal pressure. The perturbation analysis in Section 6.2.2 validates the pristine-equilibrium assumption, not this identification. Because the NTP fraction scales as v_gas^2, the predicted profiles and the peak-vs-duty-cycle relation are directly sensitive to this assumption. Please validate the mapping against simulations (for example, the Perseus-like TNG-Cluster runs of Truong et al. 2024 or the jet simulations of Bourne & Sijacki 2017), or clearly state that the model predicts the NTP only under this identification.","section":"Section 3.2.2, Eqs. (47)-(48)"},{"comment":"All jet sources with t_outburst < t_eq are removed from the population, with the justification that they would deposit all their energy in a very small volume and produce unphysically high heating rates. Figure 2 indicates that this is a substantial cut: roughly 80% of the ~60% of sources that escape the BCG, i.e., about 48% of the original sample, are removed. The paper does not quantify how the results change if these sources are retained with an alternative prescription (for example, a minimum deposition volume, or a purely thermal coupling model). The mean NTP amplitude and the shape of the peak-vs-duty-cycle relation could be sensitive to how this large subpopulation is handled. Please add a sensitivity test, or at least quantify the fraction removed and justify that their exclusion does not bias the central results.","section":"Section 2.4.4"},{"comment":"The claimed relationship between the peak NTP fraction and the AGN jet duty cycle is partly constructed rather than independently predicted. Equation (61) reduces algebraically to F approximately equal to delta/(A + delta), and the coefficients k1 and k2 in Eq. (62) are fit to the model's own peak values at only three duty cycles (10%, 20%, 30%). The resulting fit is therefore, to a large extent, a restatement of the sqrt(delta) time-averaging ansatz introduced in Eq. (41). The application to Perseus in Section 5.2.2 inverts this relation, so the inferred duty cycle inherits the assumptions of the model's time-averaging prescription. Please make this dependence explicit, and ideally test the relation against a simulated cluster whose duty cycle is known from the simulation, such as the TNG-Cluster runs.","section":"Section 5.1.1, Eqs. (61)-(62) and Table 4"},{"comment":"The manuscript contains unresolved author notes that state unfinished work: 'Will find a scaling relation to convert [delta] to a radio-loud fraction? Suggestions welcome.', 'This feels like the most significant finding, and worth discussing the correlation we find...', and 'Other things to discuss: — Gentle heating assumption in the methods. — Relationship between [delta] and halo mass...'. These passages show that central parts of the analysis, including the interpretation of the duty-cycle relation and the gentle-heating assumption, are still under development. The paper is not in a publishable state until these notes are removed and the promised analyses are either completed in the text or explicitly deferred to future work.","section":"Multiple passages (e.g., the interleaved author notes after Section 3.2.2 and before Section 5, and in the vicinity of…"}],"minor_comments":[{"comment":"The submitted text has severe typesetting problems: large duplicated blocks, repeated figure captions, interleaved page headers such as '24 Andrew Sullivan et al.', and equations numbered inconsistently in the inserted summary blocks. This must be cleaned up before any resubmission.","section":"Whole manuscript"},{"comment":"The time-averaging argument states that since NTP is proportional to the square of the heating rate, the effective heating rate carries a sqrt(delta) factor. This proportionality is asserted rather than derived from microphysics or simulation; a brief justification or a reference to supporting simulations would strengthen the presentation.","section":"Section 3.1.3"},{"comment":"The abstract states that the mean NTP fraction peaks at ~4-6% for duty cycles of 10-30%, but the text also notes that these peaks occur at only ~3 kpc and may be obscured by the BCG. This caveat should appear in the abstract or at least in the conclusions, since it affects observational applicability.","section":"Abstract and Section 4.1.2"},{"comment":"The sentence excluding the upper Hitomi bound (F ~ 11-12%) says 'if the gas is undergoing sloshing motion', but does not explain why sloshing invalidates that particular bound. Please clarify the reasoning or cite the relevant discussion.","section":"Section 5.2.2"},{"comment":"The Perseus-like cluster uses a different virial mass and radius (r500 = 1.3 Mpc, M500 = 5.8e14 Msun) than the 'typical cluster' used elsewhere (Table 2), so comparisons with the other environments are not directly on the same normalization. This is acceptable, but it should be stated explicitly when the fits are compared in Table 4 and Figure 7.","section":"Section 5.2.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears to be an early draft: it contains embedded author notes, repeated text blocks, and unfinished statements about the central duty-cycle relation. The square-root-delta time-averaging is an ad hoc assumption that is not independently tested, and the Perseus inference hinges on it. I would encourage the editor to require a revision that addresses the ergodicity issue and the unfinished-manuscript problem before sending the paper back to referees."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Genuinely new: population-based forward model of AGN-induced NTP in cluster cores, with FR-I jets sampled from observed power and age distributions and evolved through a two-phase lobe model. The NTP-duty cycle relation applied to Perseus is a new application. The derivation from the first law to the velocity kick is transparent, the Monte Carlo population is well described, and the code is on GitHub. The sensitivity runs on jet power distribution, BCG escape calibration, and the converging perturbation of hydrostatic equilibrium add real value.\n\nThe soft spots are mostly at the inference stage. The Perseus duty-cycle estimate assumes the ensemble-mean NTP equals the time-averaged NTP for a single cluster. That is stated in Section 5.2.1 but not tested. The sqrt(delta) time-averaging in Eq. (41) already presumes time-average is the observable. Hitomi gives a snapshot; if turbulence decays in the ~7 Myr off-phase, the observed NTP is biased low and so is the inferred duty cycle. This is the load-bearing concern for the paper's central claim.\n\nAlso, identifying dr/dt from the first law as the gas velocity is a modelling leap. The perturbation analysis shows internal consistency of the NTP fraction, not that this velocity equals the turbulent dispersion. Second, they remove t_outburst < t_eq sources without a sensitivity test, so the impact of that cut is unknown. Third, the restriction to epsilon in [0,0.5] is made after seeing the divergent epsilon=1 profiles; they argue those cores are inside the BCG, but it is still a post hoc choice. Fourth, the manuscript contains author notes ('suggestions welcome', 'this feels like the most significant finding') scattered through the text. That is a hygiene issue rather than a scientific one, but it should be cleaned.\n\nThe NTP-duty cycle relation itself is largely the assumed sqrt(delta) scaling rearranged: Eq. (61) reduces to F ~ delta/(A+delta), and k1 and k2 are fit to the model's own peak values. So it is not an independent prediction, though the forward comparison with Hitomi is still a useful consistency check.\n\nAudience: cluster feedback theorists and XRISM observers. This deserves a serious referee, but I would want the ergodicity issue addressed before acceptance.","headline":"Useful and transparent population model for AGN-driven NTP, but the Perseus duty-cycle inference rests on an untested ensemble-to-time-average equivalence that should be addressed before publication.","tokens_in":45624,"tokens_out":2931,"would_cite":false,"duration_ms":29703,"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":"AGN jet outbursts create a small, duty-cycle-dependent non-thermal pressure in cluster cores, so its peak fraction can reveal past jet activity.","keywords":["AGN feedback","non-thermal pressure","galaxy clusters","jet duty cycle","FR-I jets","intracluster medium","Perseus cluster","kinetic AGN feedback"],"falsifier":"Measure the NTP profile in several clusters whose jet duty cycles are independently known from radio bubble counts, as done for Perseus, and check whether the inferred duty cycles match the bubble-counting values; systematic disagreement would falsify the claimed relation. A sharper test would be a hydrodynamic simulation with a prescribed on/off jet cycle, checking whether the time-averaged core NTP fraction matches the paper's predicted peak for that duty cycle.","tokens_in":44515,"feed_emoji":"🌌","tokens_out":7402,"duration_ms":67713,"temperature":0.7,"pith_summary":"This paper argues that the non-thermal pressure (NTP) left in the core of a galaxy cluster by AGN jet outbursts is small but measurable: for typical cluster gas and dark matter profiles, the mean NTP fraction peaks at roughly 4-6% of the total pressure when the AGN jets are active for 10-30% of the AGN lifecycle. This prediction matches X-ray observations, which find 2-7% in Perseus, and supports the picture that kinetic AGN feedback does not dominate the core pressure balance. The paper then shows that the height of this peak rises with the AGN jet duty cycle in a systematic way, so a measured NTP fraction can be converted into the fraction of time the jets have been active. Applied to the Perseus cluster, the model infers a jet duty cycle of about 13% (or 3-48% given the observed range), consistent with independent estimates from the cluster's X-ray bubbles. If correct, this turns X-ray NTP measurements into a new observational probe of past AGN activity in clusters.","feed_headline":"Jet outbursts leave cluster cores only 4-6% non-thermal pressure","feed_subtitle":"The peak pressure fraction tracks the jet duty cycle, turning X-ray pressure maps into a history of AGN activity.","key_machinery":"The machinery is a Monte Carlo population of Fanaroff-Riley type I jets, sampled from a broken power-law in jet power and a pink-noise power-law in jet age, each evolved with a two-phase analytical jet-lobe model: a ballistic phase with $R(t) \\propto t^{2/(4-\\varepsilon)}$ until the jet reaches pressure equilibrium with the ambient cluster gas, followed by a flaring phase with $R(t) \\propto t^{1/(3-\\varepsilon)}$. Each outburst's energy is spread through the jet lobe volume, and the first law of thermodynamics applied to spherical gas shells yields a 'velocity kick', $v_{\\rm gas}$; the NTP fraction is then the ratio of the induced kinetic energy density to the total energy density. A calibration based on hydrodynamic simulations removes jets that cannot escape the central galaxy, and a population of compact sources is added to keep the duty cycle definition consistent with the observed radio-loud fraction.","core_discovery":"The central claim is that the peak of the mean non-thermal pressure fraction, defined as $F \\equiv p_{\\rm nt}/p$, in a cluster core is set by the AGN jet duty cycle $\\delta \\equiv t_{\\rm on}/(t_{\\rm on}+t_{\\rm off})$, and that this peak is small for realistic conditions. For a typical cluster with an NFW dark matter halo and a weakly cusped gas profile, the predicted peak values are $\\langle F\\rangle \\simeq 4.1\\%$, $5.4\\%$, and $6.3\\%$ for $\\delta = 10\\%$, $20\\%$, and $30\\%$, respectively. The paper also claims that this relationship is environment-dependent but robust across a range of non-cool core and cool core clusters, and that the combined predictions agree with existing constraints from the Hitomi satellite, Dupourqué et al., and early XRISM results. Applying the relation to the observed NTP in Perseus yields a jet duty cycle of about 13% (or 3-48% within the observed range), which is consistent with independent evidence of Perseus' recent jet activity.","pith_inferences":["Inference: if the relation holds, high-resolution X-ray spectroscopy of many clusters could produce a statistical map of jet duty cycles as a function of halo mass and cool-core status, effectively turning NTP maps into an AGN activity census.","Inference: clusters with similar NTP but no detectable jet activity, such as merger-dominated systems, should map to near-zero duty cycle, so the method could help separate jet-driven turbulence from merger-driven turbulence.","Inference: the duty cycle inferred from NTP could be cross-calibrated with the independently measured radio-loud fraction, providing a new check on the relationship between jet activity and the observable radio-loud population.","Inference: a direct hydrodynamic test would be to prescribe a known on/off jet cycle in a simulation and compare the time-averaged core NTP fraction against the paper's predicted peak for that duty cycle; agreement would strengthen the interpretation of NTP as a duty-cycle probe."],"forward_implications":["In typical clusters, the mean NTP fraction peaks at roughly 4-6% for jet duty cycles of 10-30%, matching the Hitomi and early XRISM observations and implying that AGN jet feedback is not the dominant source of non-thermal pressure in cluster cores.","Because the peak NTP fraction increases with duty cycle in every non-cool core and cool core environment considered, cooler and cuspier gas profiles are more susceptible to kinetic AGN feedback.","Combining the new core predictions with earlier large-scale NTP profiles gives a complete radial NTP fraction profile, with a pronounced dip around 100-200 kpc where the core contribution fades and the outer contribution grows.","For the Perseus cluster, the observed NTP implies a jet duty cycle of about 13%, with a range of 3-48% allowed by the measurements; this is consistent with independent estimates from the cluster's X-ray bubbles.","The predicted NTP fraction changes by only a few percent when the cluster's hydrostatic equilibrium is perturbed to account for the NTP itself, indicating a convergent steady-state description of the core pressure balance."],"supporting_citations":[{"why":"provides the single NTP fraction measurement (about 4%) used to infer Perseus' duty cycle.","marker":"Hitomi Collaboration et al. (2016)"},{"why":"provides the 2-7% (up to 11-13%) NTP range against which the predicted peak profiles are compared.","marker":"Hitomi Collaboration et al. (2018)"},{"why":"constrains the low-power slope ($s_Q = 1$) and the pink-noise age slope ($s_t = 1$) of the sampled jet population.","marker":"Shabala et al. (2020)"},{"why":"constrains the high-power slope ($s_Q = 1.5$) of the broken power-law jet power distribution.","marker":"Quici et al. (2025)"},{"why":"motivates the broken power-law shape of the jet power distribution from the radio luminosity function.","marker":"Kaiser & Best (2007)"},{"why":"supplies the observationally calibrated energetics and lifetimes of local radio AGN used in the population model.","marker":"Turner & Shabala (2015)"},{"why":"provides the two-phase FR-I jet expansion dynamics used to evolve each outburst.","marker":"Turner & Shabala (2023)"},{"why":"supplies the fitted gas density and temperature profiles that define the Perseus-like cluster.","marker":"Churazov et al. (2003)"},{"why":"calibrates how long a jet takes to escape the brightest cluster galaxy, setting the minimum jet power.","marker":"Young et al. (2024)"},{"why":"provides the large-scale NTP profiles that are combined with the new core predictions to make a complete radial profile.","marker":"Sullivan et al. (2024b)"}],"fun_headline_variants":["AGN jets add only 4-6% non-thermal pressure to cluster cores","Cluster cores see just 4-6% pressure from AGN jets","Non-thermal pressure peak tracks AGN jet duty cycle","Jet duty cycle read from non-thermal pressure in clusters","AGN jet outbursts leave a 4-6% pressure imprint"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The population average over many modelled outbursts is treated as the time average in a single real cluster, and the measured extra pressure is assumed to come from jets rather than from gas sloshing or mergers.","fun_headline_variants_meta":{"raw":{"variants":["AGN jets add only 4-6% non-thermal pressure to cluster cores","Cluster cores see just 4-6% pressure from AGN jets","Non-thermal pressure peak tracks AGN jet duty cycle","Jet duty cycle read from non-thermal pressure in clusters","AGN jet outbursts leave a 4-6% pressure imprint"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000323,"raw_usage":{"total_tokens":1870,"prompt_tokens":1060,"completion_tokens":810,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":676,"completion_tokens_details":{"reasoning_tokens":717}},"tokens_in":676,"tokens_out":810,"duration_ms":7462,"temperature":1.0,"reasoning_tokens":717,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:02:02.892472+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the NTP profile in several clusters whose jet duty cycles are independently known from radio bubble counts, as done for Perseus, and check whether the inferred duty cycles match the bubble-counting values; systematic disagreement would falsify the claimed relation. A sharper test would be a hydrodynamic simulation with a prescribed on/off jet cycle, checking whether the time-averaged core NTP fraction matches the paper's predicted peak for that duty cycle.","supporting_citations":[],"review_version":2}