{"id":"6f0e270b-a0e0-4e5f-8ca8-bf13337256c7","arxiv_id":"2607.06024","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":5,"one_line_summary":"A review arguing that cosmic-ray escape from sources, near-source regions, and galaxies is a self-regulated nonlinear process that shapes spectra and sets maximum energies, yet is often modeled with ad hoc boundary conditions.","lead":"This paper reviews how cosmic rays escape from their sources, near-source regions, and host galaxies, arguing that escape is a self-regulated nonlinear process often treated with ad hoc assumptions. A generalist might read it to understand why the origin of cosmic rays remains unsettled despite decades of study.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The emergent-halo thesis and its observational support (spectral features at 300 GV and 20 TV) depend on the Lee & Völk (1973) NLLD prescription, which Schroer et al. (2025, 2026) show is incomplete; if the corrected damping rates are used, the halo may not self-consistently emerge.","rationale":"The reader correctly identified the NLLD microphysics as the weakest link. My review confirms this is the most load-bearing concern: the paper's central thesis (emergent halo, self-regulated escape across scales) depends on a specific NLLD formulation that the paper's own cited references show is incomplete. However, the verdict should remain CONDITIONAL rather than moving to REJECT, for three reasons: (1) This is explicitly a review/perspective article for a special issue, not a research paper claiming new quantitative results; the bar for completeness is different. (2) The paper does acknowledge the NLLD uncertainty honestly, even if it does not fully resolve it. (3) The thesis that escape is nonlinear and self-regulated is supported by multiple independent lines of evidence (near-source gamma-ray observations, TeV halos, B/C constraints) that do not all depend equally on the NLLD details. The concern lands on the quantitative claims (spectral features, halo emergence) rather than on the qualitative organizing principle. The CONDITIONAL verdict with MODERATE confidence is appropriate: the perspective is valuable but the quantitative conclusions are not yet secure, and the paper would benefit from more clearly flagging which claims depend on the unsettled NLLD microphysics and which do not.","tokens_in":15446,"tokens_out":2327,"duration_ms":83761,"concrete_test":"Recompute the Galactic CR spectrum at Earth in the Dogiel et al. (2020) / Chernyshov et al. (2022) framework (boundary conditions at infinity, no prescribed halo size) but replacing the Lee & Völk (1973) NLLD rate with the Schroer et al. (2025, 2026) prescription that includes damping contributions from all smaller-k modes and a spectrum of pre-existing turbulence. If the spectral features at 300 GV and 20 TV rigidity disappear or shift by more than a factor of ~2 in rigidity, the emergent-halo picture loses its primary observational support and the paper's central thesis is substantially weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most distinctive claim — that the Galactic halo is an emergent structure determined by the balance of self-generated turbulence and NLLD damping (§4) — rests on calculations (Dogiel et al. 2020; Chernyshov et al. 2022, 2024) that use the Lee & Völk (1973) NLLD rate. Schroer et al. (2025, 2026), cited by the paper itself, demonstrate two things that directly threaten this: (1) the NLLD damping rate at wavenumber k receives contributions from all modes with smaller k, not just the resonant mode as in the Lee & Völk formulation; (2) when pre-existing turbulence is included, it enhances NLLD of self-generated modes so drastically that 'self-generated transport on Galactic scales' is strongly suppressed. The paper acknowledges this but states that the Dogiel/Chernyshov models 'remain viable' without demonstrating why. If the corrected NLLD rates are substantially larger, the critical distance z*(E) where diffusion and advection times equalize shifts inward, the effective halo shrinks, and the spectral features at 300 GV and 20 TV rigidity — the main observational anchor for the emergent-halo picture — may disappear or migrate. This is not a peripheral detail; it is the microphysical foundation of the paper's central organizing thesis. The UHECR confinement claim (§4, Blasi & Amato 2019; Cermenati et al. 2026) is less directly affected since it invokes the non-resonant instability rather than NLLD, but it too depends on the overall self-regulation framework being valid. The paper is a review/perspective and is appropriately cautious in places, but the gap between 'NLLD is unsettled and may drastically change things' and 'the emergent halo is the correct picture' is not bridged.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This review article by Blasi and Amato discusses the role of cosmic-ray (CR) escape in shaping the spectra observed at Earth, covering three scales: escape from SNR shocks (§2), from near-source regions (§3), and from the host Galaxy (§4). The central thesis is that CR escape is a nonlinear, self-regulated process in which escaping particles excite streaming instabilities that in turn confine them, and that this self-regulation—rather than ad hoc boundary conditions—sets maximum energies, shapes source spectra, and determines the structure of the Galactic halo. The paper synthesizes results from the authors' own body of work and that of others, covering the non-resonant Bell instability, the Ptuskin & Zirakashvili source spectrum prescription, near-source grammage and cocoons, and the emergent-halo picture based on self-generated turbulence balanced by nonlinear Landau damping (NLLD).","tokens_in":16292,"tokens_out":1476,"duration_ms":162839,"significance":"The topic is timely and important: the question of how CRs escape their sources and galaxies is central to connecting acceleration theory to observed spectra, and the unifying perspective across scales (shock, near-source, galactic) is valuable. The paper correctly attributes the standard equations (Eqs. 1–5) and honestly acknowledges open problems, including the unsettled NLLD microphysics and the lack of confirmed PeVatron SNRs. The discussion of observational anchors (DAMPE spectral features, LHAASO diffuse emission, LHAASO upper limits around Cas A, Auger UHECR suppression) grounds the theoretical discussion. The UHECR confinement claim (§4, Blasi & Amato 2019; Cermenati et al. 2026) is a falsifiable prediction tied to the non-resonant instability around luminous galaxies.","major_comments":[{"comment":"§4, penultimate paragraph: The paper's most distinctive claim—that the Galactic halo is an emergent structure determined by the balance of self-generated turbulence and NLLD damping (Dogiel et al. 2020; Chernyshov et al. 2022, 2024)—is acknowledged to rest on the Lee & Völk (1973) NLLD prescription. The paper itself cites Schroer et al. (2025, 2026), which show that (1) the NLLD damping rate at wavenumber k receives contributions from all modes with smaller k, not just the resonant mode as in Lee & Völk, and (2) pre-existing turbulence can 'drastically reduce the importance of self-generated transport on Galactic scales.' The paper states that the Dogiel/Chernyshov models 'remain viable' without demonstrating why. Since the spectral features at 300 GV and 20 TV rigidity—the main observational anchor for the emergent-halo picture—depend on the location of z*(E) where diffusion and advecst","section":null},{"comment":"§4, penultimate paragraph (continued): times equalize, and since z*(E) is sensitive to the NLLD rate, a quantitative or even semi-quantitative argument for why the corrected NLLD rates do not destroy the emergent-halo structure is needed. At minimum, the paper should state explicitly whether the Dogiel/Chernyshov calculations use the Lee & Völk rate or the corrected rate, and what the expected shift in z*(E) would be. This is the microphysical foundation of the paper's central organizing thesis and should be addressed before publication.","section":null}],"minor_comments":[{"comment":"§2, Fig. 1 caption: The left panel shows Emax(t) with vertical lines for the Sedov-Taylor transition and the non-resonant-to-resonant transition, but the axis labels and units are not fully described in the caption. Adding explicit labels for the vertical lines and clarifying the x-axis units would make the figure self-contained.","section":null},{"comment":"§2, Eq. (5): The numerical coefficient ~100 is derived assuming Λ ≃ 10 and n ~ 5 e-folds, but the dependence on n is not shown explicitly in the equation. Since n is a free parameter, stating the assumed value in the equation caption or inline would help the reader.","section":null},{"comment":"§3, Fig. 2: The cocoon grammage is described as '0.4 g/cm² (constant in energy with an exponential cutoff at rigidity 20 TV)' but the figure caption says R0 = 20 TV. The notation R0 is introduced only in the caption, not in the main text. Using consistent notation between text and caption would improve clarity.","section":null},{"comment":"§4: The transition from discussing NLLD microphysics to the UHECR confinement claim (Blasi & Amato 2019; Cermenati et al. 2026) is abrupt. A transitional sentence explaining how the emergent-halo picture connects to the circumgalactic non-resonant instability would help the reader.","section":null},{"comment":"§5: The summary states that 'the halo is not a static structure of prescribed size, but rather an emerging one' as a conclusion, but given the unresolved NLLD issue discussed in §4, the summary should acknowledge this caveat more explicitly rather than presenting the emergent-halo picture as established.","section":null},{"comment":"The reference list includes several 2026-dated papers (Cermenati et al. 2026; Capanema et al. 2026; Schroer et al. 2026). The journal should verify that these are properly published or in-press at the time of the manuscript's final version, as DOIs and volume/page numbers may require updating.","section":null},{"comment":"§3: The discussion of turbulent damping cites Farmer & Goldreich (2004) and notes that its importance has been 'questioned by Cerri (2024),' but does not elaborate on the nature of the disagreement. A brief clarification of what Cerri (2024) found would help readers unfamiliar with that work.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a review by leading figures in the field, and much of the cited work is their own. This is not inappropriate for a review, but the NLLD issue is the one place where the authors' own work (Schroer et al. 2025, 2026) potentially undermines their own earlier conclusions (Dogiel/Chernyshov models). The paper's handling of this tension—acknowledging it but asserting viability without argument—is the main weakness. A revised version should either provide a brief quantitative argument or reframe the emergent-halo claim as a hypothesis pending resolution of the NLLD question. The paper is well-suited to the journal's review format."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The referee raises one substantive point, concerning the viability of the emergent-halo picture in light of the corrected NLLD rates from Schroer et al. (2025, 2026). We agree that this deserves a more explicit treatment in the manuscript and will revise accordingly.","responses":[{"response":"The referee is correct that the current text asserts viability of the Dogiel/Chernyshov models without providing the supporting argument. We will revise the manuscript to address this explicitly. To be precise about the physics: the Schroer et al. (2025, 2026) results show that the Lee & Völk (1973) rate underestimates the NLLD damping rate at a given wavenumber k, because modes with smaller k also contribute. This means that self-generated waves are damped more efficiently than in the Lee & Völk prescription, which would push z*(E)—the transition height where diffusive and advective transport times become comparable—closer to the disk. The key question is whether this shift is modest enough that the emergent-halo structure survives, or whether it is so large that self-generation is effectively suppressed on Galactic scales. The Dogiel et al. (2020) and Chernyshov et al. (2022, 2024) calculations use the Lee & Völk rate, not the corrected rate. We will state this explicitly in the revised manuscript. We cannot, at present, provide a full quantitative recalculation of z*(E) with the corrected NLLD rate, as this requires solving the full wave-CR transport problem self-consistently—a task that is beyond the scope of this review. However, we can offer the following semi-quantitative argument for why the models remain viable rather than being ruled out: (1) The Schroer et al. finding that pre-existing turbulence can 'drastically reduce the importance of self-generated transport' applies specifically to the regime where a substantial pre-existing turbulent cascade is present at the relevant wavenumbers. In the Dogiel/Chernyshov models, the situation is different: self-generated waves dominate at the resonant wavenumbers for sub-TeV CRs, and the pre-existing turbulence at thes","revision_made":"partial","referee_comment":"§4, penultimate paragraph: The paper's most distinctive claim—that the Galactic halo is an emergent structure determined by the balance of self-generated turbulence and NLLD damping (Dogiel et al. 2020; Chernyshov et al. 2022, 2024)—is acknowledged to rest on the Lee & Völk (1973) NLLD prescription. The paper itself cites Schroer et al. (2025, 2026), which show that (1) the NLLD damping rate at wavenumber k receives contributions from all modes with smaller k, not just the resonant mode as in Lee & Völk, and (2) pre-existing turbulence can 'drastically reduce the importance of self-generated transport on Galactic scales.' The paper states that the Dogiel/Chernyshov models 'remain viable' without demonstrating why. Since the spectral features at 300 GV and 20 TV rigidity—the main observational anchor for the emergent-halo picture—depend on the location of z*(E) where diffusion and advec[]"},{"response":"We agree with the referee that the manuscript should, at minimum, state explicitly which NLLD prescription the Dogiel/Chernyshov calculations employ, and provide a semi-quantitative discussion of the expected impact of the corrected rates. We will add a dedicated paragraph to §4 addressing the following points: (a) The Dogiel et al. (2020) and Chernyshov et al. (2022, 2024) calculations use the Lee & Völk (1973) NLLD rate. (b) The corrected rate from Schroer et al. (2025, 2026) is larger, because it includes contributions from all modes with k' < k. This would reduce the height z*(E) at which self-generation dominates over damping, i.e., the effective halo size would shrink. (c) The magnitude of this shift depends on the spectrum of pre-existing turbulence, which is itself uncertain. In the limit where pre-existing turbulence is subdominant at the resonant wavenumbers for ~GeV–TeV CRs—a regime that is physically motivated in the inner halo where self-generation is strongest—the additional damping from non-resonant modes is a correction of order unity rather than an order-of-magnitude effect, and the qualitative structure of the emergent halo is preserved. (d) In the opposite limit, where a strong pre-existing turbulent cascade is present, the Schroer et al. results indicate that self-generated transport could be substantially suppressed, and the emergent-halo picture would need to be revisited. We will state clearly that a definitive answer requires repeating the Dogiel/Chernyshov calculations with the corrected NLLD rate, which is an important direction for future work. We acknowledge that our current statement that the models 'remain viable' is stronger than what we can rigorously demonstrate, and we will temper it accordingly, framing it as a plausible but not yet定量地","revision_made":"yes","referee_comment":"§4, penultimate paragraph (continued): times equalize, and since z*(E) is sensitive to the NLLD rate, a quantitative or even semi-quantitative argument for why the corrected NLLD rates do not destroy the emergent-halo structure is needed. At minimum, the paper should state explicitly whether the Dogiel/Chernyshov calculations use the Lee & Völk rate or the corrected rate, and what the expected shift in z*(E) would be. This is the microphysical foundation of the paper's central organizing thesis and should be addressed before publication."}],"tokens_in":15107,"tokens_out":1321,"duration_ms":71146,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"This is a review/perspective, not a new-result paper, and it should be judged as one. The organizing thesis — that CR escape is a nonlinear, self-regulated process from the shock scale up to the circumgalactic medium, and that the Galactic halo is emergent rather than prescribed — is a useful and defensible framing. The paper does a good job connecting the pieces: the non-resonant instability and maximum energies at SNR shocks (§2), the source spectrum shape and its bumps around 5–30 TeV matching DAMPE features, near-source self-confinement and grammage constraints from B/C (§3), and the emergent-halo picture with boundary conditions at infinity (§4). The Cas A / LHAASO discussion is timely and the point that LHAASO upper limits constrain the maximum energy is well taken. The equations quoted for the non-resonant instability are standard and correctly attributed. The synthesis itself has value — I haven't seen these three scales (source, near-source, galactic) discussed together in one place with this level of coherence. Credit where it's earned: the paper is clearly written and the logical chain from shock acceleration through escape to galactic transport is laid out honestly, including the admission that the source spectrum prescription matters and that different recipes give different results. The stress-test concern about NLLD lands, and it's the main soft spot. The paper's most distinctive claim — that the halo self-consistently emerges from the balance of self-generated turbulence and NLLD damping — depends on calculations (Dogiel et al. 2020; Chernyshov et al. 2022, 2024) that use the Lee & Völk (1973) NLLD rate. Schroer et al. (2025, 2026), which the paper itself cites, show two things that threaten this: (1) the NLLD rate at wavenumber k receives contributions from all smaller-k modes, not just the resonant one, and (2) pre-existing turbulence enhances NLLD so drastically that self-generated transport on Galactic scales may be strongly suppressed. The paper acknowledges both points but then states the Dogiel/Chernyshov models 'remain viable' without any argument for why. That's a gap in the load-bearing part of the synthesis. If the corrected NLLD rates are substantially larger, the critical distance z*(E) shifts inward, the effective halo shrinks, and the spectral features at 300 GV and 20 TV — the main observational anchor — may move or disappear. The paper should either sketch why the Dogiel/Chernyshov framework survives or flag this more prominently as an open question. The UHECR confinement claim (§4) is less directly affected since it invokes the non-resonant instability rather than NLLD, but it remains speculative and depends on the overall self-regulation framework holding. A few minor issues: the paper doesn't explicitly label itself as a review/perspective, which could mislead readers into thinking the synthesis is novel. The self-citation rate is high but mostly justified since the authors are central to this program. This paper is for cosmic-ray theorists and transport modelers who want a coherent map of the escape problem across scales. It's a competent and useful review with one important unresolved tension at its core. It deserves a serious referee — mainly to check attribution accuracy, fairness of the NLLD discussion, and whether the 'remain viable' claim should be softened or substantiated.","headline":"Solid review of multi-scale CR escape, but the emergent-halo thesis rests on unsettled NLLD microphysics that the paper acknowledges but doesn't resolve.","tokens_in":16393,"tokens_out":1307,"would_cite":false,"duration_ms":70067,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.85.Ry","98.38.Mz","52.35.Py"],"model":"glm-5.2","headline":"Cosmic ray escape is self-regulated, not passive","keywords":[],"falsifier":"If nonlinear Landau damping rates are substantially different from current models (as recent hybrid-PIC simulations suggest they may be), then self-generated confinement on galactic scales weakens, the halo ceases to be emergent in the way described, and the connection to UHECR suppression breaks.","tokens_in":15472,"feed_emoji":"🧲","tokens_out":1238,"duration_ms":48563,"temperature":0.7,"pith_summary":"This paper argues that the escape of cosmic rays from their acceleration sites, their near-source neighborhoods, and their host galaxies is a nonlinear, self-regulated process rather than a passive leakage. The central mechanism is a feedback loop: escaping particles carry electric currents that excite plasma instabilities (both resonant and non-resonant streaming instabilities), which amplify magnetic fields and in turn scatter and confine those same particles. At supernova remnant shocks, this self-excitation of turbulence sets the maximum achievable energy and shapes the time-integrated spectrum released into the interstellar medium, which is not a clean power law but carries bumps and dips in the 5-30 TeV range. Around sources, the steep cosmic ray gradient excites further instabilities that suppress diffusion, accumulating excess grammage and producing extended gamma-ray emission. On galactic scales, the authors argue that the cosmic ray halo is not a fixed box with prescribed boundary conditions but an emergent structure, self-consistently determined by the balance between self-generated turbulence and its damping (particularly nonlinear Landau damping). They further propose that this self-confinement extends to the circumgalactic medium around luminous galaxies, potentially explaining the observed suppression of ultra-high-energy cosmic ray flux around 10^18 eV.","feed_headline":"Escaping cosmic rays build their own cages","feed_subtitle":"A review argues that cosmic ray escape is a self-regulated process at every scale, from shock to galaxy, shaping spectra and halos.","key_machinery":"The feedback loop connecting escaping cosmic ray currents to plasma instability growth (resonant and non-resonant streaming instabilities), magnetic field amplification, enhanced scattering, and confinement. Nonlinear Landau damping serves as the primary saturation mechanism for self-generated turbulence on galactic scales, and its microphysics determines whether self-confinement dominates over pre-existing turbulence.","core_discovery":"The paper's central claim is that cosmic ray escape at every scale, from the shock to the galaxy to the circumgalactic medium, is governed by the same self-regulating mechanism: escaping particles excite streaming instabilities that confine them, and this feedback, not the acceleration mechanism alone, determines maximum energies, source spectra, transport properties, and the very size of the galactic halo. The galactic halo is recast as an emergent structure arising from the coupling between cosmic ray currents and interstellar plasma, rather than a pre-existing container with fixed boundaries.","pith_inferences":[],"forward_implications":["The source spectrum injected into galactic transport calculations is not a simple power law but carries structure shaped by escape dynamics, meaning transport models that assume clean power-law injection are missing physically motivated spectral features.","If the galactic halo is emergent rather than fixed, standard cosmic ray propagation codes that impose free-escape boundaries at prescribed locations are solving a different problem than nature poses, and their parameter fits may not correspond to physical quantities.","Self-confinement of cosmic rays around luminous galaxies could suppress the flux of particles below 10^18 eV, offering a physical explanation for the spectral features seen in Auger data without requiring exotic injection spectra.","The lack of detected PeVatron candidates among young supernova remnants, combined with theoretical limits on magnetic field amplification, raises the question of whether standard SNRs can produce the knee-region cosmic ray flux at all.","Near-source suppressed diffusivity contributes to diffuse gamma-ray and neutrino emission that is not accounted for in standard diffuse emission models, potentially affecting interpretations of LHAASO and Fermi-LAT data."],"fun_headline_variants":["Cosmic ray escape is self-regulated at every scale","Escaping cosmic rays generate their own confinement","Streaming instabilities set cosmic ray maximum energies","Cosmic ray escape shapes source spectra and galactic halos","The galactic halo emerges from cosmic ray feedback"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The paper's quantitative conclusions depend on how well nonlinear Landau damping is modeled, and the authors acknowledge this microphysics is unsettled: recent simulations show that damping rates depend on all wave modes with smaller wavenumber, not just the resonant mode, and that pre-existing turbulence can drastically reduce the importance of self-generated transport on galactic scales. If the damping physics differs from what is assumed, the conclusions about halo size, 1","fun_headline_variants_meta":{"raw":{"variants":["Cosmic ray escape is self-regulated at every scale","Escaping cosmic rays generate their own confinement","Streaming instabilities set cosmic ray maximum energies","Cosmic ray escape shapes source spectra and galactic halos","The galactic halo emerges from cosmic ray feedback","Cosmic ray currents build the cages that contain them","Self-confinement governs cosmic ray escape from shock to galaxy","Cosmic ray escape is a nonlinear feedback loop at all scales","Streaming instabilities determine cosmic ray transport and spectra"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":550,"prompt_tokens":390,"completion_tokens":160,"prompt_tokens_details":null},"tokens_in":390,"tokens_out":160,"duration_ms":2558,"temperature":1.0,"reasoning_tokens":35,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T18:57:25.820158+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If nonlinear Landau damping rates are substantially different from current models (as recent hybrid-PIC simulations suggest they may be), then self-generated confinement on galactic scales weakens, the halo ceases to be emergent in the way described, and the connection to UHECR suppression breaks.","supporting_citations":[],"review_version":1}