{"id":"db3268a2-c820-4486-980e-9d4ddf612276","arxiv_id":"1908.10367","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"Streaming cosmic-ray heating makes halo-gas thermal instability an overstability with an Alfvén-scale frequency, and yields a Schwarzschild-form convection criterion involving gas plus cosmic-ray pressure.","lead":"This paper works out the linear stability of hot halo gas heated by cosmic rays streaming along magnetic fields. It shows that cosmic-ray pressure changes thermal instability from a pure growth into an oscillation, and gives a new criterion for convection in such gas.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; central results are internally consistent within the stated CR transport model.","rationale":"The reader's ACCEPT verdict with moderate confidence is appropriate. I checked the linearization algebra in Section 4, the high-beta reduction leading to eq. (25), the asymptotic limits in Section 4.5, and the CR-heated background relations in Section 5; no sign error or internal inconsistency surfaced. The only substantive caveat is the CR transport closure, which is also the reader's weakest-assumption. It is a physical assumption, not a mathematical flaw, and the authors explicitly state that super-Alfvenic streaming is not included because its dependence on fluid quantities is uncertain. The central claims are conditional on that closure: if super-Alfvenic streaming dominates, quantitative predictions such as the oscillation frequency and diffusion-stability boundaries would likely need revision, but the eta-controlled growth-rate structure and the qualitative PIE/CIE conclusions are expected to be more robust. This does not undermine the paper as written, so the verdict should remain unchanged.","tokens_in":27533,"tokens_out":19772,"duration_ms":212692,"concrete_test":"Independently re-derive the high-beta quadratic dispersion relation (25) by imposing delta_pc = -delta_pg in eqs. (23)-(24), then compare its roots to the exact eigenvalues of eqs. (20)-(24) at beta = 100 over 10^-2 < eta < 10^2; if any root deviates by more than a few percent, the high-beta reduction hides a term. As a separate transport test, repeat the uniform-medium stability calculation with a super-Alfvenic streaming closure v_st = -v_st0 sign(b_cap dot grad P_c) and v_st0 > v_A, and check whether the overstability frequency and the stability boundaries shift by order unity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"No load-bearing defect found. The linearized equations (20)-(24), the high-beta quadratic dispersion relation (25), and the asymptotic limits (28)-(29) are mutually consistent, and the numerical eigenvalues agree with the analytic limits shown in Figure 2. The strongest physical caveat is the CR transport model: streaming exactly at the Alfven speed with a constant parallel diffusion coefficient kappa (Section 2.2, eqs. 4-5). If super-Alfvenic streaming or non-diffusive transport dominates, the perturbed CR pressure response changes, so the predicted overstability frequency and the diffusion-sensitive stability boundaries in Figures 3-4 would shift. The authors explicitly acknowledge this in Section 2.2 and Section 7, and the paper's headline claims are scoped to the adopted model. Because the eta-dependent isobaric/isochoric growth rates, the CR-heating-induced overstability, and the PIE/CIE stability conclusions are all derived within that stated model, I do not find an unsupported or internally inconsistent step.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a linear thermal stability analysis of dilute plasma heated by streaming cosmic rays, modeled with the standard gas-CR fluid equations (1)-(5). It treats three equilibria: a uniform background with unspecified heating balancing cooling (Section 4), a background in which CR heating balances radiative cooling (Section 5), and a gravitationally stratified CR-heated background (Section 6). The central variables are eta = p_c/p_g and beta, and the characteristic frequencies omega_c, omega_a, omega_d, and omega_ff. The main results are: (i) thermal instability growth rates depend strongly on eta, with perturbations isobaric for eta<<1 and isochoric for eta>>1; (ii) the perturbed CR heating term introduces oscillations at the Alfven frequency, making thermally unstable modes overstable for k dot B != 0; (iii) CR diffusion suppresses a finite band of modes and can introduce a CR Field length; (iv) gas in photoionization equilibrium is thermally stable for any eta, whereas CIE halo gas is generally unstable; and (v) a stratified CR-heated medium obeys a Schwarzschild-like convective instability criterion d s_eff/dz < 0. The analytic asymptotic limits (28)-(29), (40)-(41), and (43)-(44) are checked against the numerical solutions shown in Figures 2-3.","tokens_in":27800,"tokens_out":10698,"duration_ms":107911,"significance":"This is a careful and largely self-contained derivation that corrects and extends earlier heuristic treatments (e.g., Pfrommer 2013) by perturbing the CR energy equation consistently. The explicit linearization (20)-(24), the high-beta quadratic (25), and the asymptotic limits are mutually consistent, and the agreement with the full numerical solution in Figure 2 supports the central claim. The results are physically useful: they give falsifiable stability boundaries in the (eta, omega_d/omega_a) plane, identify the CR Field length, and predict an overstability frequency proportional to the Alfven speed. A notable strength is that the stability boundaries are derived rather than fitted; the only fits, equations (35) and (47), are descriptive broken power laws for the numerically computed Lambda_T,c. The main caveat is the adopted CR transport model: streaming exactly at v_A with a constant parallel diffusion coefficient kappa (Section 2.2). If super-Alfvenic streaming or non-diffusive transport dominates, the perturbed CR heating response changes and the predicted overstability frequencies and diffusion-sensitive boundaries would shift.","major_comments":[],"minor_comments":[{"comment":"In footnote 3, 'equlibrium' should be 'equilibrium', and the notation delta Q/Q should be defined explicitly or by example, since it is not immediately clear which fluid variables are included in Q.","section":"Section 4.1, footnote 3"},{"comment":"In Eq. (49) and the surrounding paragraph, delta appears to denote Lagrangian perturbations, whereas everywhere else in the paper (e.g., Eqs. 20-24 and Appendix C) delta denotes Eulerian perturbations. Please state this explicitly or use a different symbol such as Delta for Lagrangian perturbations; as printed, the triple equality in Eq. (49) is confusing and could be misread as an inconsistency.","section":"Section 6.3.1, Eq. (49)"},{"comment":"The sentence 'CR diffusion does not suppress the overall excitation of thermal instability' is ambiguous because Figure 3 shows a broad region where diffusion does suppress instability. Rephrase to indicate that diffusion does not suppress the instability of all modes, but rather of a range of modes.","section":"Section 4.6 and Appendix B"},{"comment":"There is an extra comma inside the parenthetical factor: '(omega_c/10^-15 s^-1,)^-1/2' should read '(omega_c/10^-15 s^-1)^-1/2'.","section":"Eq. (34)"},{"comment":"The acronym CIE is used without definition at first use; define 'collisional ionization equilibrium (CIE)' when it first appears, and keep the terminology consistent with the abstract.","section":"Abstract and Section 7"},{"comment":"The statement that 'the cosmic-ray pressure equation (5) then implies that p_c proportional to rho^2/3' is not immediately obvious; add a brief derivation or a pointer to the relevant equations so the reader can verify this background relation.","section":"Section 5.1"},{"comment":"The full dispersion relation is not written out, and the numerical eigenvalue solutions depend on MATLAB; for reproducibility, consider providing the code as supplementary material or an appendix.","section":"Sections 4-6"}],"recommendation":"minor_revision","confidential_remarks":"This is a well-executed linear-stability paper well within the scope of MNRAS. The transport-model caveat is real but explicitly scoped, and the internal derivation is sound. I support publication after the minor clarifications above."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth a serious look. The core new result is that if you evolve the CR pressure perturbation instead of assuming an adiabatic p_c–rho relation, thermal instability in CR-heated halos becomes an overstability with an oscillation frequency set by the Alfvén speed, and the growth-rate dependence on eta = p_c/p_g changes substantially. The authors also give a CR-diffusion Field length and a Schwarzschild-form convective criterion for stratified, CR-heated backgrounds. These are genuine extensions beyond Pfrommer 2013 and Wiener et al. 2013b, and the correction to the perturbed CR heating term in earlier work is specific and testable.\n\nThe paper does the linear theory carefully. The perturbed equations are written out in full, the high-beta quadratic dispersion relation is checked against the full numerical matrix solution, and the asymptotic limits reduce properly to the standard isobaric and isochoric TI growth rates. The PIE/CIE distinction is cleanly scoped, and the broken power-law fits to Lambda_T,c are descriptive, not load-bearing. I did not find an internal inconsistency.\n\nThe largest soft spot is not mathematical but physical: the adopted CR transport model, streaming exactly at v_A with a constant parallel diffusion coefficient. The authors flag in Section 2.2 and again in Section 7 that for many damping mechanisms the residual transport is not actually diffusive and can be super-Alfvénic. If super-Alfvénic streaming dominates, the perturbed CR heating term, and therefore the overstability frequency and the diffusion-sensitive stability boundaries, would change. So the headline claims are conditional on that transport model. That is an honest limitation, and they state it, but it is the reason I would not call the PIE/CIE conclusions universal.\n\nA smaller caveat: the uniform-medium calculation uses an ordering with a small background CR pressure gradient (epsilon) and then drops background gradients; the paper notes the breakdown when delta p_c/p_c > (k H_c)^-1. That is a minor technical point, not a fatal one.\n\nWho is this for? Anyone working on multiphase gas in the CGM/ICM, CR feedback, or thermal instability in magnetized plasmas. I would cite it for the corrected perturbed heating term and the CR Field length. It deserves peer review; my own verdict is accept, with the understanding that the transport assumption is the main thing a referee should probe.","headline":"A careful and explicit linear stability analysis that corrects earlier CR-heated TI work and maps out the parameter space; the main caveat is the assumed Alfvén-speed streaming transport.","tokens_in":28315,"tokens_out":1618,"would_cite":true,"duration_ms":19777,"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 shows that thermal instability in halo gas heated by streaming cosmic rays is controlled by the ratio of cosmic-ray pressure to gas pressure, and that the unstable modes oscillate at the Alfvén speed rather than growing…","keywords":["thermal instability","cosmic-ray streaming","Alfvén speed","galaxy halos","multiphase gas","cosmic-ray diffusion","photoionization equilibrium","collisional ionization equilibrium"],"falsifier":"Numerically solve the linearized gas-plus-cosmic-ray equations in a uniform, high-$\\beta$, CR-heated background with $\\kappa=0$, cooling slope $\\Lambda_T=-1$, and $\\eta=1$. The paper predicts the gas entropy mode has frequency $\\omega\\approx-(2/3)\\omega_a-(2/3)i\\Lambda_T\\omega_c$, so the unstable mode oscillates at the Alfvén frequency; if the mode is purely growing with no real oscillatory part at finite $\\eta$, the central overstability claim is wrong.","tokens_in":27346,"feed_emoji":"🌌","tokens_out":13371,"duration_ms":122070,"temperature":0.7,"pith_summary":"The paper asks whether thermal instability can still fragment hot halo gas when that gas is heated by cosmic rays streaming down their pressure gradient. It argues that the answer is controlled by one number, $\\eta=p_c/p_g$: at small $\\eta$ perturbations are isobaric and the usual growth rates apply, while at large $\\eta$ they are isochoric. The perturbed cosmic-ray heating term does not change growth rates much, but it makes every thermally unstable mode oscillate at a frequency of order the Alfvén speed, so the instability is formally an overstability even in a uniform medium. This matters because galaxy halos, groups, and clusters show cold multiphase gas whose origin may be thermal instability; the paper shows when CR heating permits that instability and when it suppresses it.","feed_headline":"Thermal instability in halo gas oscillates at the Alfvén speed","feed_subtitle":"Growth and stability hinge on cosmic-ray pressure, deciding when halo gas forms cold clouds.","key_machinery":"The machinery is the linearized system built from the gas and cosmic-ray pressure equations. In the high-$\\beta$ limit the dispersion relation reduces to a quadratic in the mode frequency $\\omega$, coupling $\\delta p_c$ and $\\delta p_g$ through pressure balance $\\delta p_c\\approx-\\delta p_g$. The organizing parameter is $\\eta=p_c/p_g$, which decides whether the eigenmode is isobaric or isochoric. The frequency $\\omega_a=\\mathbf{k}\\cdot\\mathbf{v}_A$ appears through the perturbed CR heating and supplies the real part of the entropy mode, while $\\omega_d=\\kappa(\\hat{\\mathbf{b}}\\cdot\\mathbf{k})^2$ parametrizes diffusion and can damp a band of unstable wavenumbers, giving a CR Field length $\\lambda_{\\rm CRF}\\sim 2\\pi|\\hat{\\mathbf{b}}\\cdot\\hat{\\mathbf{k}}|\\sqrt{\\eta\\kappa/\\omega_c}$ for $\\eta<1$ when $\\kappa\\omega_c/(\\eta v_A^2)\\lesssim1$.","core_discovery":"The central claim is that in a plasma heated by streaming cosmic rays, linear thermal stability is governed by the ratio of cosmic-ray pressure to gas pressure, $\\eta=p_c/p_g$, through the coupling of the perturbed CR pressure to density. In the high-$\\beta$ limit the entropy mode obeys a quadratic dispersion relation combining the gas and CR energy equations. For $\\eta\\ll1$ the mode is isobaric with growth rate $(2/5)(2-\\Lambda_T)\\omega_c$; for $\\eta\\gg1$ it is isochoric with growth rate $-(2/3)\\Lambda_T\\omega_c$. In both regimes CR streaming introduces an oscillation frequency $\\omega_a=\\mathbf{k}\\cdot\\mathbf{v}_A$ in the gas entropy mode, so thermal instability is an overstability. CR diffusion can suppress instability in an intermediate band of $\\eta$ and wavenumber, acting effectively like thermal conduction with a CR Field length, but only when $\\kappa\\omega_c/(\\eta v_A^2)\\lesssim1$. In gravitationally stratified halos, the convective-instability criterion takes the Schwarzschild form $d s_{\\rm eff}/dz<0$ with $s_{\\rm eff}=\\ln(p_g/\\rho^{5/3})+\\eta\\ln(p_c/\\rho^{4/3})$, and the entropy-mode oscillation frequency can exceed the free-fall frequency.","pith_inferences":["A direct test of the overstability claim is a nonlinear CR-MHD simulation of a cooling parcel with streaming at $v_A$: if the saturated condensation develops coherent oscillations at the local Alfvén frequency, the linear prediction is borne out.","The CR Field length implies a minimum clump size when streaming is subdominant; comparing observed sizes of cold circumgalactic-medium clouds with $\\lambda_{\\rm CRF}$ for adopted $\\kappa$ and $\\eta$ could constrain the diffusion coefficient.","Because photoionization-equilibrium gas is predicted to be always stable while collisional-ionization-equilibrium gas is usually unstable, the mere existence of cold gas in a halo becomes a diagnostic of the dominant ionization and heating balance, apart from CR pressure.","If the correct transport is super-Alfvénic streaming rather than $v_A$ streaming plus constant diffusion, the perturbed heating term changes and the predicted oscillation frequencies and stability maps would need to be recomputed; this is the natural next step the paper leaves open."],"forward_implications":["In collisional-ionization-equilibrium halo gas, which has a cooling slope $\\Lambda_T<0$, cosmic-ray heating does not stabilize thermal instability for realistic $\\eta$; multiphase cold gas can still condense.","Thermally unstable perturbations in CR-heated gas are not stationary growing modes but waves traveling at roughly the Alfvén speed once $\\eta\\sim1$, so the condensation process carries an intrinsic oscillation.","Cosmic-ray diffusion suppresses thermal instability only in a limited range of parameters, and for $\\kappa\\omega_c/(\\eta v_A^2)\\gtrsim1$ even high-wavenumber perturbations remain isobarically unstable; there is no universal CR Field length.","In gravitationally stratified halos the CR-driven entropy-mode frequency can exceed the free-fall frequency, which would shift the critical $t_{\\rm cool}/t_{\\rm ff}$ quoted for multiphase gas formation.","A stratified CR-heated medium is convectively unstable exactly when $d s_{\\rm eff}/dz<0$, giving a Schwarzschild-like criterion that depends only on gas and CR pressure scale heights."],"supporting_citations":[{"why":"provides the baseline thermal-instability dispersion relations and the Field length that the CR-diffusion analogue extends.","marker":"Field 1965"},{"why":"establishes the streaming instability that ties cosmic-ray transport to the Alfvén speed.","marker":"Kulsrud & Pearce 1969"},{"why":"derives the cosmic-ray transport equations with streaming and diffusion used as the paper's starting point.","marker":"Skilling 1971"},{"why":"gives the form of the cosmic-ray heating rate $-v_A\\cdot\\nabla p_c$ that drives the perturbed heating.","marker":"Wentzel 1971"},{"why":"identifies damping regimes in which cosmic-ray transport is non-diffusive and super-Alfvénic, framing the model's applicability.","marker":"Wiener et al. 2013a"},{"why":"shows many damping mechanisms make transport advective rather than diffusive, justifying the caution about the diffusion term.","marker":"Wiener et al. 2018"},{"why":"provides the earlier heuristic cosmic-ray-heated thermal-instability calculation whose perturbed heating this paper corrects.","marker":"Pfrommer 2013"},{"why":"supplies the effective cooling function for photoionization equilibrium used to conclude PIE halo gas is thermally stable.","marker":"Wiersma et al. 2009"},{"why":"first considered thermal instability with cosmic-ray heating in cooling flows, the problem this paper generalizes.","marker":"Loewenstein et al. 1991"}],"fun_headline_variants":["CR pressure ratio tilts halo gas into instability","Halo gas overstability driven by streaming cosmic rays","Alfvén-speed oscillations govern halo gas thermal stability","Cosmic-ray streaming reshapes halo gas multiphase formation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusions assume cosmic rays stream down their pressure gradient at exactly the Alfvén speed, with any residual transport captured by a constant diffusion coefficient along the magnetic field; if super-Alfvénic or non-diffusive transport dominates, the predicted heating perturbations and stability boundaries would change.","fun_headline_variants_meta":{"raw":{"variants":["CR pressure ratio tilts halo gas into instability","Halo gas overstability driven by streaming cosmic rays","Alfvén-speed oscillations govern halo gas thermal stability","Cosmic-ray streaming reshapes halo gas multiphase formation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000238,"raw_usage":{"total_tokens":1607,"prompt_tokens":1135,"completion_tokens":472,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":751,"completion_tokens_details":{"reasoning_tokens":407}},"tokens_in":751,"tokens_out":472,"duration_ms":5283,"temperature":1.0,"reasoning_tokens":407,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:45:51.571832+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the linearized gas-plus-cosmic-ray equations in a uniform, high-$\\beta$, CR-heated background with $\\kappa=0$, cooling slope $\\Lambda_T=-1$, and $\\eta=1$. The paper predicts the gas entropy mode has frequency $\\omega\\approx-(2/3)\\omega_a-(2/3)i\\Lambda_T\\omega_c$, so the unstable mode oscillates at the Alfvén frequency; if the mode is purely growing with no real oscillatory part at finite $\\eta$, the central overstability claim is wrong.","supporting_citations":[{"cited_title":"G., 1971, , 163, 503","cited_arxiv_id":null,"evidence_quote":"gives the form of the cosmic-ray heating rate $-v_A\\cdot\\nabla p_c$ that drives the perturbed heating."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the effective cooling function for photoionization equilibrium used to conclude PIE halo gas is thermally stable."}],"review_version":1}