{"id":"fec7caa6-93e3-4f90-b524-b91cefcc8768","arxiv_id":"2502.01753","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Cosmic rays from past black hole growth can rival gas pressure at the outskirts of group-mass halos, potentially driving baryons out and reshaping halo outskirts.","lead":"This paper argues that cosmic rays produced by supermassive black holes during their growth can build up enough pressure in the outer parts of galaxy groups to push gas out of the halos. The idea matters because it could explain why matter in the universe appears less clumped than expected, the so-called S8 tension.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim depends on a point-source delta-function injection model (Eq. 3) that may double-count the CR energy budget when paired with the assumed efficiency; the cumulative CR energy inside R200 likely exceeds the halo's thermal energy in groups, so an energetic consistency check is needed.","rationale":"The paper is transparent about its order-of-magnitude nature and identifies CR transport as the major uncertainty; the reader's weakest_assumption focuses on the diffusion coefficient window. I agree that the transport window is the primary parametric uncertainty, but the most load-bearing structural issue is the normalization of Eq. 3 relative to Eq. 2. The paper's own Figure 2 shows E_CR,BH can exceed E_vir for groups, so if Eq. 3 yields p_CR ~ p_th at R200, the implied CR energy inside R200 must be checked for consistency with the assumed total budget. If the integrated energy is much smaller than E_CR,BH, the conclusion is actually more robust (only a fraction of the CRs are needed); if it is of order E_CR,BH, then the model is self-consistent but the CRs are mostly inside R200, somewhat undermining the 'outskirts' narrative. This is a mechanical normalization check, not a fundamental flaw, so the verdict stays CONDITIONAL. The concern is distinct from but partially overlapping with the reader's transport-window worry: even within the quoted kappa window, the energetic integral must be verified.","tokens_in":23,"tokens_out":1965,"duration_ms":136311,"concrete_test":"Recompute Eq. 3 numerically: integrate p_CR(r,t) over a sphere of radius R200 for M200 = 1e13 Msun, kappa = 1e30 cm2/s, t = 10 Gyr, and E_CR,BH = 1e60 erg. Compare the resulting integrated CR energy inside R200 to E_CR,BH and to the thermal energy within R200 from the Arnaud et al. (2010) profile. If the integrated CR energy exceeds E_CR,BH or is not within ~30% of it, the normalization is inconsistent; if it is much less than E_CR,BH, then the CR energy that produces p_CR ~ p_th at R200 is only a small fraction of the available budget and the claim is strengthened, not weakened. Also check the same integral for the cluster case M200 = 8e14 Msun.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim is that CRs from BH growth can have p_CR ~ p_th near R200 in ~10^13-14 Msun halos. The quantitative basis is Eq. 3, a point-source, instantaneous delta-function diffusion solution with E_CR,BH from Eq. 2. Within the window kappa ~ 3e29-1e31 cm2/s highlighted in Figure 3, the CR energy remains largely interior to R200 at t ~ 10 Gyr (diffusion radius ~ (kappa t)^1/2 ranges from ~0.3-1.7 R200), so p_CR ~ p_th near R200 implies a total CR energy within R200 of order the thermal energy of the halo. However, Eq. 2 gives E_CR,BH ~ 1e60 erg for a 1e13 Msun halo, which is comparable to or larger than the total thermal energy E_vir of the halo (as the paper itself shows in Figure 2). If p_CR ~ p_th at R200, the volume-integrated CR energy inside R200 would be of order E_th, and the model would put essentially all of E_CR,BH inside the halo, meaning the CRs are not actually being transported to and beyond R200 in the way the narrative suggests. The issue is not that E_CR,BH is too small, but that Eq. 3 with the quoted normalization is not self-consistently normalized to E_CR,BH unless the Gaussian is integrated over all space; when evaluated at R200 with the quoted parameters, the resulting p_CR must be checked against the integral constraint. A related concern is that the empirical thermal pressure profile used (Arnaud et al. 2010) is extrapolated from cluster-mass systems to groups and to R200; if the true group thermal pressure at R200 is lower or higher than this extrapolation, the required E_CR shifts accordingly. Nonetheless, the dominant load-bearing step is the energetic consistency of Eq. 3 with Eq. 2.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents analytic order-of-magnitude estimates for the pressure exerted by cosmic rays (CRs) produced by supermassive black hole accretion in massive halos. It argues that, for CR diffusion coefficients kappa ~ 3e29-1e31 cm^2/s and CR energy injection efficiencies eps_BH ~ 1e-3-1e-2 of the black hole rest energy, the CR pressure at R200 in ~1e13-1e14 M_sun halos can be comparable to the thermal gas pressure. The authors discuss CR transport via streaming, compare with Fermi gamma-ray upper limits in clusters and groups, estimate non-thermal emission detectability, and discuss implications for the baryon distribution, the S8 tension, and kSZ observations. The paper is explicitly framed as a plausibility argument rather than a definitive calculation, and all central results are presented with transparent scaling relations.","tokens_in":23744,"tokens_out":33538,"duration_ms":293178,"significance":"If the mechanism operates, it introduces a long-lived, spatially extended feedback channel in group-mass halos that is absent from most cosmological simulations, and it would provide a concrete physical route to strong feedback at large halo radii. The paper's strengths are its transparent analytic framework, explicit scaling relations (e.g., p_CR proportional to eps_BH kappa^{-3/2}), and use of observational constraints such as Fermi gamma-ray upper limits. The authors consistently flag the dominant uncertainties: eps_BH, the CR transport coefficient, the plasma beta, and the extrapolation of the Arnaud et al. pressure profile. I checked the energetic-consistency concern raised in the stress test: Eq. (3) is correctly normalized to E_CR,BH over all space, and for the fiducial parameters (M200=1e13 M_sun, eps_BH=3e-3, kappa=3e30 cm^2/s, t=10 Gyr) only about 20% of E_CR,BH lies inside R200, with a corresponding CR energy that is a modest fraction of the halo thermal energy. That specific concern therefore does not materialize.","major_comments":[],"minor_comments":[{"comment":"The delta-function diffusion solution is properly normalized to E_CR,BH, but the paper should state explicitly what fraction of E_CR,BH lies within R200 for the fiducial parameters. For M200=1e13 M_sun, eps_BH=3e-3, kappa=3e30 cm^2/s, and t=10 Gyr, this fraction is roughly 20%, and the CR energy inside R200 is a modest fraction of the thermal energy. Adding this one-line bookkeeping check would preempt concerns about apparent over-injection of CR energy within the halo.","section":"2.2, Eq. (3)"},{"comment":"In the sentence introducing the CR injection efficiency, the ordering 'eps_BH ~ 10^-2 - 10^-3' should read '10^-3 - 10^-2' to agree with Eq. (2) and Figure 3, and with the subsequent discussion that treats 10^-3 as the lower end and 10^-2 as the upper end of the range.","section":"2.1"},{"comment":"The paper notes that the effective isotropic diffusivity is likely a factor of a few smaller than the field-aligned streaming value, but it does not quantify how this factor shifts the viable range of kappa in Figure 3. Since the factor is degenerate with the uncertain plasma beta, a sentence stating the resulting shift would help the reader connect Eq. (10) to the window highlighted in Figure 3.","section":"3, after Eq. (10)"},{"comment":"Equation (13) normalizes the predicted gamma-ray intensity to the Arnaud et al. (2010) self-similar pressure at R200. Given the paper's own caveat that this profile is extrapolated from cluster-mass systems to groups, please add a sentence quantifying how the detectability estimate changes if the true group thermal pressure at R200 differs from the adopted normalization by a factor of two or three.","section":"4, Eq. (13)"},{"comment":"In the Appendix, the sentence 'Equation A3 gives the usual result' should refer to Eq. (A2); there is no separately labeled Eq. (A3). In addition, the notation for the CR injection efficiency should be harmonized between the text (eps_BH) and the Figure 3 caption (eps_CR,BH), and a few typographical errors (e.g., 'signficantly' in Section 5) should be corrected.","section":"Appendix and notation"}],"recommendation":"minor_revision","confidential_remarks":"For the editor: this is a speculative but well-argued order-of-magnitude paper. The central claim is conditional on uncertain parameters, but the paper is transparent about these assumptions and the analytic estimates are internally consistent. The stress-test concern about normalization and energy bookkeeping does not, on closer reading, land: Eq. (3) is correctly normalized, and the CR energy inside R200 in the fiducial cases is not anomalously large. I see no grounds for rejection on technical merit; the remaining issues are presentation and clarity, so minor revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead Quataert & Hopkins. The new thing here is the specific estimate that cosmic rays produced by black hole growth over a Hubble time can yield p_CR ~ p_th at the virial radius of 10^13-14 Msun halos, but not in more massive clusters. That specific claim, with the scaling plots, is not in the earlier literature. The paper is honest about being an order-of-magnitude motivation, and it does a good job checking against Fermi gamma-ray limits on clusters and noting that groups are less constrained.\n\nThe stress-test note about Eq. 3 being mis-normalized does not hold up. The solution is the standard delta-function diffusion solution and is correctly normalized: integrate 3p over all space and you recover E_CR,BH. If p_CR ~ p_th at R200, then indeed a large fraction of E_CR,BH is inside the halo at that time; that is the point, and it is consistent with Fig. 2's energetics. So no double-counting there. The note's energetic consistency check is actually satisfied by construction.\n\nThe real soft spots are the ones the authors flag. The ratio scales linearly with eps_BH, the fraction of BH rest mass going into CRs, which is uncertain by an order of magnitude. The thermal pressure profile is extrapolated from clusters. And the diffusion coefficient or streaming speed window that gives p_CR ~ p_th at R200 is exactly the uncertain regime. But the paper says all of this. It also gives a plausible transport estimate (streaming at Alfven speed) that lands in the needed range, which is the main reason to take the idea seriously. The S8 and kSZ connection is speculation, and they label it as such.\n\nWho is this for? People working on baryon feedback in groups, CR transport, and weak lensing systematics. It is a useful paper to motivate simulations. I would send it to a good referee; it's a well-posed proposal with no internal errors I could find. I'd cite it if I wrote about baryon distributions in groups.\n\nBest.","headline":"A transparent, well-scoped order-of-magnitude case that BH-generated CRs can matter at R200 in groups; the central claim is parameter-dependent but the paper says so.","tokens_in":24309,"tokens_out":2293,"would_cite":true,"duration_ms":22650,"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":"The paper argues that cosmic rays produced during the growth of a massive central black hole accumulate for billions of years and end up exerting pressure comparable to the hot gas pressure in the outskirts of galaxy groups (halos of…","keywords":["cosmic rays","black hole feedback","galaxy groups","circumgalactic medium","baryon distribution","virial radius","kinetic Sunyaev-Zeldovich effect","S8 tension"],"falsifier":"A stacked gamma-ray search of ~$10^{13}$–$10^{14}$ solar-mass galaxy groups that measures the cosmic-ray pressure fraction at the virial radius at less than about 10 percent of the thermal pressure would rule out the mechanism; today's strongest limits, from Fermi observations of nearby groups, already start to exclude the largest diffusion coefficients. A cheaper check is a cosmological simulation with streaming transport: if the effective diffusion coefficient realized at R200 in group-mass halos falls outside the $10^{30}$–$10^{31}$ $cm^{2}$/s window, the predicted pressure ratio collapses.","tokens_in":23123,"feed_emoji":"🌌","tokens_out":10669,"duration_ms":91855,"temperature":0.7,"pith_summary":"The authors argue that cosmic rays generated by the growth of a central massive black hole accumulate over cosmic time and, by today, exert a pressure comparable to the thermal gas pressure in the outskirts of galaxy groups, halos of roughly $10^{13}$–$10^{14}$ solar masses. If correct, this pressure would push gas from near the virial radius out past it, redistributing baryons around groups and changing what future surveys infer about structure growth. Because the cosmic rays were produced mostly at redshifts 1–3 and take billions of years to travel outward, this feedback is set by the whole history of black hole activity in a halo rather than by its current cooling rate or current accretion state. That makes it the kind of mechanism that standard cosmological simulations, which couple feedback to the instantaneous mechanical output of black holes, would systematically miss. The authors connect it to long-standing anomalies in weak lensing and kinetic Sunyaev-Zeldovich measurements that hint at unusually strong feedback at large halo radii.","feed_headline":"Black-hole cosmic rays can push gas out of galaxy groups","feed_subtitle":"Pressure at the virial radius of 10^13–10^14 solar-mass halos can rival thermal gas and push baryons out.","key_machinery":"The argument runs on one ratio, p_CR/p_th at r = R200, with the cosmic-ray pressure given by the Green's-function solution for isotropic diffusion from a point burst, p_CR(r,t) ≃ E_CR/[3(4πκt)^{3/2}] exp(−r²/4κt), and the thermal pressure from the empirical universal pressure profile of Arnaud et al. (2010). The other load-bearing piece is a transport estimate: if the cosmic rays stream down their pressure gradient at the Alfvén speed, the effective isotropic diffusion coefficient is κ ~ 6×$10^{31}$ $β^{{−1/2}}$(r/R200)(M200/$10^{13}$ M☉)^{2/3} cm²/s, which for plausible plasma $\\beta$ values lands inside the window that makes p_CR/p_th ~ 1 at R200 in groups while keeping clusters untouched.","core_discovery":"During its growth a massive black hole releases of order $10^{60}$ erg of energy into cosmic rays, about ten times the supernova cosmic-ray budget of the same halo and comparable to the total virialized thermal energy of a $10^{12}$–$10^{13}$.5 solar-mass halo. Modeled as a burst of cosmic-ray energy diffusing outward from the halo center for about 10 Gyr, this reservoir produces a cosmic-ray pressure that is within a factor of order unity of the thermal pressure at the virial radius in $10^{13}$–$10^{14}$ solar-mass halos, provided the effective diffusion coefficient is near $10^{30}$–$10^{31}$ $cm^{2}$/s. The same calculation applied to $10^{14}$.5 solar-mass clusters keeps the ratio far below unity, consistent with Fermi gamma-ray limits of p_CR/p_th ≲ 0.02–0.04 in clusters. The paper's central claim is that this accumulated cosmic-ray pressure is a previously underappreciated, dynamically important component at large radii in group-mass halos, acting even after the jets and winds that made the cosmic rays have shut off.","pith_inferences":["A concrete prediction not worked out in the paper: the baryon fraction of groups should show a sharp radial deficit that grows toward and beyond R200 and saturates at cluster masses, a signature stackable with X-ray and Sunyaev-Zeldovich profiles.","If groups later merge into clusters, their cosmic-ray-loaded gas would deliver a pre-heated, non-thermal pressure component to cluster outskirts, so evidence of this mechanism might appear in clusters as a small pressure excess at large radii rather than at the center.","The fate of the mechanism hinges on the composition of jets; if relativistic jets are primarily electron-positron pairs, only the ~0.1–1 GeV leptons survive to large radii, so the effective energy budget may shrink, an uncertainty the paper explicitly leaves open.","A simulation-based test: run a 10^13–10^14 solar-mass halo with streaming cosmic-ray transport and read the effective diffusion coefficient at R200; values of κ outside 10^30–10^31 cm^2/s would move the predicted pressure ratio far from unity."],"forward_implications":["In group-mass halos cosmic-ray pressure becomes a long-lived feedback channel that stays active for gigayears after the black hole activity that produced the cosmic rays has ended.","Cosmological simulations that model only mechanical black hole feedback will systematically overpredict the gas retained inside 10^13–10^14 solar-mass halos and underpredict baryons beyond R200.","The baryon content of a group's outskirts depends on the integrated growth history of its black hole, decoupling large-scale gas properties from the current radiative cooling rate of the inner halo.","Removing baryons from near the virial radius suppresses the matter power spectrum at k ~ 1 Mpc^{-1}, offering a single mechanism that can address both the weak lensing S8 tension and the kSZ evidence for diffuse, extended gas in halo outskirts.","Clusters of about 10^14.5 solar masses and above remain essentially unaffected, so the mechanism naturally breaks the degeneracy between groups and clusters and survives existing Fermi constraints."],"supporting_citations":[{"why":"Supplies the empirical universal thermal pressure profile p_th(M200, r) that the cosmic-ray pressure is compared against.","marker":"Arnaud et al. (2010)"},{"why":"Fermi stacked-cluster gamma-ray upper limits that constrain p_CR/p_th to about 0.02–0.04 in massive clusters.","marker":"Ackermann et al. (2014)"},{"why":"Stacked gamma-ray analysis providing the cluster pressure limits shown in Figure 4 and the nearby-group limit p_CR < 1e-13 erg/cm^3.","marker":"Huber et al. (2013)"},{"why":"Gives the stellar-mass–halo-mass relation used to set central stellar and black hole masses in the energy estimates.","marker":"Moster et al. (2013)"},{"why":"Black-hole–stellar-mass relation used to estimate M_BH and hence the total cosmic-ray energy from black hole growth.","marker":"Kormendy & Ho (2013)"},{"why":"Identifies collisionless damping as the dominant wave damping in hot halo plasmas, yielding the streaming speed used in the effective diffusion coefficient estimate.","marker":"Wiener et al. (2018)"},{"why":"Zoom-in simulations showing that a cosmic-ray energy fraction near 10^-3 can quench galaxies in ~10^13 solar-mass halos, motivating the fiducial energy fraction.","marker":"Wellons et al. (2023)"},{"why":"Provides the X-ray luminosity–halo mass relation used to compare radiated energy and maintenance-mode cosmic-ray pressure in Figure 2 and equation (5).","marker":"Anderson et al. (2015)"}],"fun_headline_variants":["Cosmic rays from black holes expel gas from galaxy groups","Black hole cosmic rays push gas out of group halos","Stored cosmic rays from black holes drive baryons out","Black-hole cosmic rays act long after jets switch off","Cosmic ray pressure from black hole rivals thermal in groups"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise, one the paper itself flags in §3 and §5, is that cosmic rays travel through the outer halo at an effective diffusion coefficient near $10^{30}$–$10^{31}$ $cm^{2}$/s (roughly streaming at the Alfvén speed); faster transport spreads them out before they reach the virial radius, slower transport keeps them in the inner halo, and only the middle window gives p_CR comparable to p_th there.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic rays from black holes expel gas from galaxy groups","Black hole cosmic rays push gas out of group halos","Stored cosmic rays from black holes drive baryons out","Black-hole cosmic rays act long after jets switch off","Cosmic ray pressure from black hole rivals thermal in groups"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000311,"raw_usage":{"total_tokens":1845,"prompt_tokens":1090,"completion_tokens":755,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":706,"completion_tokens_details":{"reasoning_tokens":674}},"tokens_in":706,"tokens_out":755,"duration_ms":6360,"temperature":1.0,"reasoning_tokens":674,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T14:35:36.757828+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A stacked gamma-ray search of ~$10^{13}$–$10^{14}$ solar-mass galaxy groups that measures the cosmic-ray pressure fraction at the virial radius at less than about 10 percent of the thermal pressure would rule out the mechanism; today's strongest limits, from Fermi observations of nearby groups, already start to exclude the largest diffusion coefficients. A cheaper check is a cosmological simulation with streaming transport: if the effective diffusion coefficient realized at R200 in group-mass halos falls outside the $10^{30}$–$10^{31}$ $cm^{2}$/s window, the predicted pressure ratio collapses.","supporting_citations":[],"review_version":1}