{"id":"dabda6a0-fb36-4afd-aa8a-81f07a538295","arxiv_id":"1908.08061","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Cosmic rays accelerated at a T Tauri star's accretion shock can raise disk ionization far above the Galactic background and extend the region where the magneto-rotational instability operates, possibly creating a self-regulating accretion cycle.","lead":"This paper models how fast-moving charged particles (cosmic rays) made by a young star's feeding shock ionize the disk of gas and dust around it. The extra ionization can switch on a magnetic instability that speeds up the star's feeding, possibly explaining why young stars flicker in brightness.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central MRI result rests on applying the ISM-calibrated CR attenuation law (Eq. 9) to a shock-accelerated E^-2 spectrum; a direct transport calculation with the post-shock spectrum is needed.","rationale":"The reader identifies Eq. 9 as the weakest assumption, and I agree. In good faith, the paper is transparent about the extrapolation, uses a public chemistry code, and explores several propagation and dust parameters, so there is no internal inconsistency or overreach in framing. The central claim, however, is a quantitative bridge from the shock-acceleration model to the MRI-active boundary, and Eq. 9 is the only element in that bridge that is explicitly acknowledged to be assumed rather than computed for the relevant spectrum. The MRI boundary at Sigma ~ 1 g/cm2 is near the column where the attenuation exponent has the largest leverage, so even a moderate error in the exponent changes the claimed extension toward the midplane by a meaningful factor. Since the test could go in either direction, this concern does not justify rejecting the paper, but it does justify keeping the reader's CONDITIONAL verdict. The proposed direct transport calculation is concrete, uses existing tools, and would settle whether the central MRI claim stands at its stated strength.","tokens_in":21151,"tokens_out":12404,"duration_ms":135062,"concrete_test":"Run the Padovani et al. (2018) CR transport code, or an equivalent continuous-slowing-down calculation with secondary production, using the post-region-II shock spectrum from Figure 2 as the input spectrum. Tabulate zeta_CR(N_H) for N_H = 1e18-1e27 cm^-2 and fit the effective power-law index over the disk-relevant range N_H = 6e22-6e25 cm^-2. If the fitted index differs from -0.34 by more than about 0.05, recompute the MRI-active boundaries in Figures 6-8 with the new attenuation curve and report whether the boundary at Sigma ~ 1 g/cm2 inside 20 au shifts by more than a factor of roughly 3 in surface density.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.4.1 adopts the Padovani et al. (2018) attenuation law zeta_CR = zeta_0 (N_H / 1e18 cm^-2)^-0.34 (Eq. 9) inside the disk. That law was calibrated for interstellar CR background spectra, not for the E^-2 accretion-shock spectrum of Gaches & Offner (2018b) used here. The shock spectrum is also reshaped in regions I and II by energy losses in the accretion flow, inverse funneling, and secondary production (Figure 2), so its low-energy cutoff and slope at the disk surface are not the ISM case. The authors explicitly state 'we assume the attenuation will be qualitatively similar in our case,' which is a testable assumption rather than a derived result. The MRI-active boundary at Sigma ~ 1 g/cm2 (N_H ~ 6e23 cm^-2) sits in the column range where the exponent -0.34 matters: changing the exponent from -0.34 to -0.6 lowers zeta by roughly a factor of 30 at that column, moving the CR-dominated active region to substantially lower surface density and weakening the 'extends toward the midplane' claim. The fit is also applied or extrapolated beyond N_H ~ 1e25 cm^-2 for the MMSN midplane, outside the stated range. A direct transport calculation could strengthen or weaken the headline result; without it, the quantitative MRI boundary remains conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the gas-grain chemistry code UCLCHEM to study whether cosmic rays accelerated at the accretion shock of a T Tauri star can ionize its circumstellar disk. Starting from the shock-accelerated E^-2 spectrum of Gaches & Offner (2018b), the authors propagate the spectrum through the accretion flow, the magnetosphere, and the disk, and compute the local ionization rate from cosmic rays together with FUV, X-rays, and radionuclides. They then run the chemistry at fixed density and temperature structure and map the electron fraction and ion abundances, using the magnetic Reynolds number and the ambipolar diffusion number to define the MRI-active region. The main claims are that local cosmic rays produce zeta >= 1e-15 to 1e-14 s^-1 at the disk surface within about 10 au, extend the MRI-active region to Sigma ~ 1 g cm^-2 inside about 20 au, and leave the MMSN midplane MRI-dead except for diffusive cosmic-ray propagation. The paper also proposes a feedback loop connecting accretion rate, cosmic-ray attenuation, and MRI activity.","tokens_in":21444,"tokens_out":13781,"duration_ms":143341,"significance":"If the results hold, this paper would establish a local, accretion-powered cosmic-ray source as a major term in the ionization budget of protoplanetary disks, with consequences for dead-zone structure, disk chemistry, and accretion variability. The modeling is transparent and forward-looking: the parameters are listed in Table 2, the chemistry code is public, the cosmic-ray spectrum is taken from an independent acceleration calculation rather than fitted to disk observables, and the paper states its main limitations explicitly. The predictions for dominant ions, H3+ abundances, CH3CN/HCN ratios, and the location of the MRI-active boundary are falsifiable and should stimulate observations. The quantitative MRI boundary is the key deliverable, and its robustness is the main thing the revision needs to establish.","major_comments":[{"comment":"The attenuation law zeta_CR = zeta_0 (N_H / 1e18 cm^-2)^-0.34 is adopted from Padovani et al. (2018) for interstellar cosmic-ray spectra and applied to the accretion-shock E^-2 spectrum, with the justification 'we assume the attenuation will be qualitatively similar in our case.' This assumption is load-bearing: the MRI-active boundary at Sigma ~ 1 g cm^-2 corresponds to N_H ~ 6e23 cm^-2, precisely the column range where the exponent matters. Changing the exponent from -0.34 to -0.6 lowers zeta by roughly a factor of 30 at that column and moves the CR-dominated active region to substantially lower surface density. The law is also applied beyond the stated validity range N_H < 1e25 cm^-2 when the MMSN midplane is evaluated. I ask for either a direct transport calculation of the post-shock spectrum through region III using the same loss-function machinery as regions I and II, or an explicit sensitivity study over the attenuation exponent and column range, with the resulting MRI boundaries reported. Without this, the quantitative MRI-active boundary and the midplane dead-zone conclusion remain conditional.","section":"Section 2.4.1, Eq. (9)"},{"comment":"The model assumes that photons and cosmic rays enter the disk only vertically and neglects side-entry at the inner rim. The authors acknowledge in Section 4.6 that this underestimates ionization at the midplane within the inner ~10 au and may extend the MRI-active region. This is not a minor caveat: the paper's negative conclusion that the MMSN midplane remains MRI-dead except for diffusive propagation is derived in the vertical-only geometry, and the proposed feedback loop depends on where the dead zone actually sits. Please quantify the side-entry contribution, at least with an approximate treatment of the inner-edge column, or explicitly restrict the midplane-dead-zone claim to the vertical-entry geometry.","section":"Section 2.1 and Section 4.6"},{"comment":"The abstract and conclusions present a cosmic-ray feedback loop that 'mediates accretion and may produce luminosity variability,' but no time-dependent or even steady-state feedback model is presented. The paper demonstrates monotonic relationships between accretion rate, cosmic-ray production, attenuation, and MRI activity, but a negative-feedback loop that regulates accretion or produces variability requires an additional stability or timescale argument. Please add a simple estimate, such as the equilibrium accretion rate or the loop response time, or soften the wording to describe a possible feedback pathway rather than an established accretion-regulation mechanism.","section":"Section 4.1 and Conclusions"}],"minor_comments":[{"comment":"The radial dependence in Eq. (11) appears inverted: as written G0 grows as r^2, whereas Eq. (10) gives G0 proportional to r^-2. This should presumably be (3 au / r)^2, and the implementation should be checked against the intended scaling.","section":"Eq. (11)"},{"comment":"The normalization column should be written with consistent units, N_H in cm^-2, and the meaning of N_H as the vertical hydrogen column from the disk surface should be stated at first use.","section":"Eq. (9)"},{"comment":"The phrase 'a minor affect on the disk ionization' should read 'a minor effect.'","section":"Section 3.1"},{"comment":"The text 'grains colagulate and sediment' should read 'coagulate,' and 'ionization fraction of up to 10^-8 s^-1' should have dimensionless units for a fraction, not s^-1.","section":"Section 4.6"},{"comment":"In the sentence about large grains, 'Sigma ~ 10 g cm^-3' should be 'g cm^-2.'","section":"Section 3.4.2"},{"comment":"The two Gaches & Offner references appear with identical bibliographic data (ApJ 861, 87); if 2018a and 2018b are distinct papers, the page or journal data should be corrected, and if they are the same paper, one citation should be removed.","section":"References"},{"comment":"The time unit in the radionuclide ionization expression should be specified explicitly, since the exponent 1.04 t is only sensible with t in a stated unit such as Myr.","section":"Eq. (15)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of the journal and the authors are unusually transparent about their assumptions. My main concern is technical rather than conceptual: the spectrum-dependent attenuation law in region III is the load-bearing step for the quantitative MRI boundaries, and the revision should either compute it or bracket it with a sensitivity study. I have no concerns about citation practices beyond the duplicate Gaches & Offner reference noted in the minor comments."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—the headline is that this paper makes a solid, well-caveated modeling case that accretion-shock cosmic rays can raise disk ionization well above the Galactic background and shift the MRI-active boundary inward. The genuinely new thing is coupling the authors' own shock-accelerated E^-2 spectrum to a full gas-grain chemistry code (UCLCHEM) and mapping the MRI boundary, plus the feedback loop idea linking accretion rate, CR attenuation, and MRI activity. The code is public, the parameter exploration is thorough (accretion rate, hot-spot covering fraction, dust size, CR coupling, free-streaming vs diffusive propagation), and the paper is unusually transparent about its assumptions. The chemistry predictions (C+, S+, Mg+ at surface, H3+ at depth) are testable.\n\nThe soft spots are real but proportionate. The load-bearing one is Eq. 9: they apply the Padovani et al. (2018) ISM attenuation law, calibrated for the interstellar CR background, to their harder E^-2 shock spectrum, stating they assume the attenuation is 'qualitatively similar.' That is not derived, and the stress-test note quantifies the effect: at the column where the MRI boundary sits (~1 g/cm2), changing the exponent from -0.34 to -0.6 lowers the ionization rate by a factor of ~30, which would substantially move the boundary. They also extrapolate beyond the law's stated column range. This is testable: a direct transport calculation with the post-shock spectrum would settle it. The neglect of horizontal CR/X-ray entry, which they explicitly admit, likely makes their results conservative for the inner 10 au. The feedback loop is qualitative and the midplane-dead-zone conclusion depends on the assumed propagation regime. None of this sinks the paper; it is a forward model with clearly identified uncertainty.\n\nI'd send this out. The right referee will ask for the attenuation calculation, but the framework and chemistry maps are worth publishing and the caveats are honestly stated.","headline":"A useful, transparent forward model of accretion-shock CR ionization in disks, whose main quantitative claim rests on an unvalidated attenuation-law extrapolation.","tokens_in":22011,"tokens_out":2095,"would_cite":true,"duration_ms":20385,"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":"Accreting T Tauri stars accelerate cosmic rays at their accretion shocks, producing disk ionization rates at least an order of magnitude above the Galactic background.","keywords":["cosmic rays","circumstellar disks","T Tauri stars","accretion shocks","magneto-rotational instability","disk chemistry","ionization rate","protoplanetary disks"],"falsifier":"Measure the cosmic-ray ionization rate in the inner ~10 au of a T Tauri disk, for example through H3+ line emission or absorption; if the rate is consistent with the shielded Galactic background (~1e-16 $s^{-1}$) rather than the predicted zeta >= 1e-14 $s^{-1}$, the accretion-shock CR source is not effective at the claimed level.","tokens_in":1940,"feed_emoji":"🌌","tokens_out":4025,"duration_ms":96428,"temperature":0.7,"pith_summary":"An accreting young star's own accretion shock may be the dominant source of cosmic rays in its circumstellar disk. The paper argues that for T Tauri accretion rates, these locally accelerated cosmic rays produce disk-surface ionization rates at least ten times the Galactic cosmic-ray background, and up to $\\zeta \\sim 10^{-14}$ s$^{-1}$ inside 10 au. That extra ionization changes which ions dominate the disk chemistry and pushes the magneto-rotational instability (MRI) active region toward the midplane, where the MRI can drive accretion. The coupling between accretion rate, cosmic-ray attenuation by the accretion flow, and MRI activity forms a feedback loop that could naturally produce accretion variability. A sympathetic reader would care because ionization sets disk chemistry, gas-field coupling, and planet formation conditions.","feed_headline":"Accretion shocks ionize disks 10x the cosmic-ray background","feed_subtitle":"Local cosmic rays reshape disk chemistry and widen the MRI-active zone, possibly regulating accretion.","key_machinery":"The load-bearing machinery is a three-stage propagation model for an accretion-shock cosmic-ray spectrum with energy dependence $E^{-2}$ and energies up to a few GeV. It is attenuated by inverse funneling and energy losses in the ionized accretion flow between the star and the inner disk, by geometric spreading and energy losses above the disk surface, and by a power-law column-density attenuation within the disk, $\\zeta_{\\rm CR} = \\zeta_0 \\, (N_{\\rm H}/10^{18}\\,\\mathrm{cm}^{-2})^{-0.34}$, taken from interstellar CR studies. The propagated spectrum feeds a gas-grain chemistry calculation that yields electron and ion abundances, which are then converted into magnetic Reynolds and Ambipolar numbers to locate MRI-active regions. The same machinery produces the proposed feedback loop: accretion rate sets CR luminosity, accretion-column column density sets CR attenuation, and MRI activity sets accretion rate.","core_discovery":"The central discovery is that cosmic rays accelerated by the stellar accretion shock of a T Tauri star can dominate disk ionization in the inner ~10 au. For accretion rates $\\dot{M}_* \\sim 10^{-9}$ to $10^{-6}\\,M_\\odot$ yr$^{-1}$, the paper computes surface ionization rates $\\zeta \\geq 10^{-15}$ s$^{-1}$, exceeding the Galactic background by at least an order of magnitude, with $\\zeta \\geq 10^{-14}$ s$^{-1}$ inside 10 au. This CR flux raises the ionization at intermediate to high column densities ($\\Sigma > 10$ g cm$^{-2}$), makes C$^{+}$, S$^{+}$, and Mg$^{+}$ the surface ions and H$_3^+$ the ion above 1 g cm$^{-2}$, and extends the MRI-active region toward the midplane. Only with diffusive propagation ($\\zeta \\propto r^{-1}$) does the minimum-mass solar nebula midplane become MRI-active; otherwise it remains a dead zone. The paper concludes that the same accretion that accelerates CRs, the accretion column that attenuates them, and the MRI they enable form a feedback loop that can mediate accretion and drive luminosity variability.","pith_inferences":["A testable extension: molecular-line observations of H$_3^+$ in the inner ~10 au of an accreting T Tauri disk could distinguish shock-accelerated CR ionization from X-ray ionization, since the paper predicts H$_3^+$ dominates above 1 g cm$^{-2}$ with a CR-driven floor.","If the CR feedback loop operates, young-star accretion should be self-limiting on timescales of a few years at the MRI-active boundary; this is an inference, as the paper only sketches the loop.","Applying a steeper or shallower attenuation law for the $E^{-2}$ shock spectrum would shift the MRI-active boundary; a particle-transport calculation through disk gas would directly test the assumed $-0.34$ power law.","The paper's discussion implies that CR-enhanced gas-phase CO could make CO observations of accreting sources overestimate the luminosity of past accretion bursts; this follows from its discussion but is not a central claim."],"forward_implications":["For accretion rates $10^{-9}$ to $10^{-6}\\,M_\\odot$ yr$^{-1}$, shock-accelerated CRs give disk-surface ionization rates $\\zeta \\geq 10^{-15}$ s$^{-1}$, at least ten times the Galactic CR background, reaching about $10^{-14}$ s$^{-1}$ within 10 au.","CR ionization dominates over X-rays and FUV at surface densities above $10$ g cm$^{-2}$ inside roughly 10 au, shifting the ion-neutral transition and changing the dominant ions: C$^{+}$, S$^{+}$, and Mg$^{+}$ at the surface, H$_3^+$ above 1 g cm$^{-2}$.","The MRI-active region extends toward the midplane; in the fiducial model the disk is MRI-active at column densities near 1 g cm$^{-2}$ inside about 20 au, but the minimum-mass solar nebula midplane remains a dead zone.","If CRs propagate diffusively ($\\zeta \\propto r^{-1}$), the midplane can become MRI-active, enabling an accretion self-regulation loop.","At very high accretion rates the dense accretion flow attenuates the CRs, so CR feedback is strongest in T Tauri disks ($10^{-9}$ to $10^{-7}\\,M_\\odot$ yr$^{-1}$) and weaker in protostellar disks."],"supporting_citations":[{"why":"Supplies the initial accretion-shock cosmic-ray spectrum (E^-2 up to a few GeV) and shock acceleration framework.","marker":"Gaches & Offner 2018b"},{"why":"Provides the interstellar CR attenuation law zeta = zeta0 (N_H/1e18)^-0.34 used for disk penetration and the claim that attenuation is weak below N_H ~1e25 cm^-2.","marker":"Padovani et al. 2018"},{"why":"Supplies the gas-grain chemistry code that computes species abundances from the local ionization inputs.","marker":"Holdship et al. 2017"},{"why":"Supplies the 522-species chemical network and elemental abundances used in the disk chemistry calculations.","marker":"Quénard et al. 2018"},{"why":"Supplies the disk vertical structure, grazing-angle geometry, and magnetic Reynolds number formulation used to map MRI-active regions.","marker":"Perez-Becker & Chiang 2011"},{"why":"Defines the T-Tauriosphere shielding baseline and the comparison model for X-ray, FUV, radionuclide, and CR ionization and MRI activity.","marker":"Cleeves et al. 2013a"},{"why":"Sets the Galactic cosmic-ray background ionization rate (~1e-16 s^-1) against which the enhanced rates are compared.","marker":"Indriolo & McCall 2012"},{"why":"Provides the X-ray ionization prescription and the standard MRI-active region calculation that the paper's results are compared with.","marker":"Bai & Goodman 2009"},{"why":"Motivates the free-streaming versus diffusive cosmic-ray propagation limits (r^-2 versus r^-1) that determine whether the midplane can become MRI-active.","marker":"Rodgers-Lee et al. 2017"}],"fun_headline_variants":["Accretion shocks power inner-disk ionization 10x cosmic background","Stellar cosmic rays extend MRI-active zone in protoplanetary disks","Cosmic-ray feedback loop may regulate accretion and drive variability","Young star's own cosmic rays reshape disk chemistry and MRI"],"cache_read_input_tokens":24064,"weakest_assumption_plain":"The calculation assumes the interstellar cosmic-ray attenuation law (ionization scales as column density to the -0.34 power) applies to the steeper, lower-energy spectrum accelerated at the accretion shock; if the true attenuation is materially different, the computed midplane ionization and MRI boundaries shift.","fun_headline_variants_meta":{"raw":{"variants":["Accretion shocks power inner-disk ionization 10x cosmic background","Stellar cosmic rays extend MRI-active zone in protoplanetary disks","Cosmic-ray feedback loop may regulate accretion and drive variability","Young star's own cosmic rays reshape disk chemistry and MRI"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000257,"raw_usage":{"total_tokens":1677,"prompt_tokens":1139,"completion_tokens":538,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":755,"completion_tokens_details":{"reasoning_tokens":466}},"tokens_in":755,"tokens_out":538,"duration_ms":5831,"temperature":1.0,"reasoning_tokens":466,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:51:25.639117+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the cosmic-ray ionization rate in the inner ~10 au of a T Tauri disk, for example through H3+ line emission or absorption; if the rate is consistent with the shielded Galactic background (~1e-16 $s^{-1}$) rather than the predicted zeta >= 1e-14 $s^{-1}$, the accretion-shock CR source is not effective at the claimed level.","supporting_citations":[{"cited_title":"V., Galli, D., & Caselli, P","cited_arxiv_id":null,"evidence_quote":"Provides the interstellar CR attenuation law zeta = zeta0 (N_H/1e18)^-0.34 used for disk penetration and the claim that attenuation is weak below N_H ~1e25 cm^-2."},{"cited_title":"2011, ApJ, 735, 8 Cosmic-Ray Mediated Disks 19","cited_arxiv_id":null,"evidence_quote":"Supplies the disk vertical structure, grazing-angle geometry, and magnetic Reynolds number formulation used to map MRI-active regions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets the Galactic cosmic-ray background ionization rate (~1e-16 s^-1) against which the enhanced rates are compared."},{"cited_title":"2009, ApJ, 701, 737","cited_arxiv_id":null,"evidence_quote":"Provides the X-ray ionization prescription and the standard MRI-active region calculation that the paper's results are compared with."},{"cited_title":"M., Ray, T","cited_arxiv_id":null,"evidence_quote":"Motivates the free-streaming versus diffusive cosmic-ray propagation limits (r^-2 versus r^-1) that determine whether the midplane can become MRI-active."}],"review_version":1}