{"id":"0412ee60-259c-4ff7-a3ac-a39c7e9a1898","arxiv_id":"2505.07718","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Simulations find that luminosity outbursts evaporate water ice mantles and reduce pebble mass by about half, with pebble regrowth taking thousands of years.","lead":"This paper simulates how outbursts from young stars heat their surrounding protoplanetary disks and evaporate the icy mantles on dust grains. The simulations show that such outbursts destroy about half of the centimeter-sized pebbles in the disk, with recovery taking thousands of years.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The factor-of-two pebble depletion is governed by an instantaneous dust-destruction assumption that the paper itself flags as one of two competing scenarios; without an erosion run the headline remains conditional.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing point, and the paper itself flags it in Section 2.3, so this concern is not manufactured. It does not justify rejection because the instantaneous-disruption scenario is physically plausible for fragile ice-bound aggregates and is supported by the earlier Molyarova et al. (2021) modeling of mantle-stuck growth; the rest of the simulation—thermal balance, MRI-triggered outbursts, snowline shifts, and spiral sublimation geometry—is internally coherent, and the qualitative sequence of pebble loss and slow regrowth is likely robust. The problem is that the quantitative '≈2×' and 'several thousand years' recovery time are direct outputs of an untested destruction law, and the same paper cites observational and modeling work favoring the slower erosion alternative. The appropriate outcome is the reader's CONDITIONAL verdict, with the first and second conclusions explicitly restricted to the instantaneous-destruction scenario. The lack of released code or data makes independent reproduction costly, which reinforces but does not by itself drive the conditional status.","tokens_in":13212,"tokens_out":5883,"duration_ms":59596,"concrete_test":"Rerun models M1 and M2 with erosion-limited dust destruction as in Stammler & Birnstiel (2022) or with the bi-disperse erosion model of Akimkin et al. (2020), replacing the instantaneous a_max reset in Section 2.1 while keeping all other physics and outburst histories identical, and measure total pebble mass during the 244.5 kyr (M1) and ~340 kyr (M2) outbursts. If the dip remains ≈2× and the minimum occurs during the outburst, the objection is resolved; if the dip is substantially weaker (e.g. <1.3×) or is delayed until after the outburst, the Section 4 first bullet must be reworded as scenario-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative result—pebble mass dropping by ≈2× and recovering over several thousand years—is not a direct output of coagulation physics but of the instantaneous fragmentation-barrier reset described in Section 2.3. When mantles evaporate, a_frag drops by a factor (5/0.5)^2 = 100, and all grown dust above the new barrier is immediately recycled into small dust; in Section 2.1 this is implemented by setting a_max = a_frag and D = 0. The paper itself notes the alternative erosion-limited destruction (Stammler & Birnstiel 2022) operates on coagulation timescales ~10^3 yr, which can exceed the ~200 yr outburst duration, and that Houge et al. (2024) favor prolonged destruction in V883 Ori. Since post-outburst water freeze-out times are days near the snowline, a gradual-erosion scenario would have only the ~200 yr thermal spike in which to act, so the pebble population could survive largely intact. The Molyarova et al. (2021) spiral-structure argument supports mantle-stuck aggregates but does not constrain the disruption rate once mantles are lost. The Section 4 headline is therefore a scenario-dependent quantitative result, not a model-independent one.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the FEOSAD 2D thin-disk hydrodynamics code to simulate the formation and evolution of two protoplanetary disks with different initial core masses, coupling gas dynamics, dust growth and drift, MRI-driven accretion outbursts, and adsorption/desorption of H2O, CO2, CH4, and CO. It focuses on how luminosity outbursts affect snowline positions, ice mantles on grown dust (pebbles), and total pebble mass. The central reported result is that during outbursts the disk-integrated pebble mass drops by roughly a factor of two because water-ice mantles evaporate, the fragmentation velocity drops from 5 to 0.5 m/s, and dust aggregates above the new fragmentation barrier are immediately destroyed into small dust; pebble abundance then recovers on timescales of several thousand years. The paper also reports that CO2, CH4, and CO snowlines shift more than the water snowline, and that ice desorption occurs in non-axisymmetric spiral regions.","tokens_in":13472,"tokens_out":9115,"duration_ms":89516,"significance":"If robust, the half-reduction of the pebble population is significant for planet-formation theory: pebbles are central to streaming-instability planetesimal formation and pebble accretion, so episodic accretion outbursts could modulate the reservoir of pebble-sized solids on thousand-year timescales. The modeling is unusually comprehensive in coupling self-consistent MRI outbursts with multi-species ice chemistry and dust evolution, and the paper is careful to disclose the instantaneous-destruction assumption and the alternative erosion scenario. The non-axisymmetric finding that ice desorbs along spiral arms is a useful, testable prediction. The main caveat is that the quantitative headline depends on a prescribed fragility transition rather than on an emergent collisional outcome; the paper would be strengthened by a conditional framing or an erosion-limited run.","major_comments":[{"comment":"The ≈2× pebble-mass reduction is a direct consequence of the instantaneous-destruction prescription. When ice mantles evaporate, v_frag is reduced from 5 to 0.5 m/s, so a_frag in Eq. (9) drops by a factor of 100, and the implementation sets a_max = a_frag and D = 0, immediately recycling all larger grown dust into small dust. The paper explicitly acknowledges that an erosion-limited destruction scenario (Stammler & Birnstiel 2022) operates on ~10^3 yr timescales and that Houge et al. (2024) favor prolonged destruction for V883 Ori; because the modeled outburst lasts ~200 yr, the erosion scenario would destroy far fewer pebbles. The abstract and the first conclusion bullet state the half-reduction without this caveat, so the headline is scenario-dependent and needs to be qualified there, or supported by an explicit erosion-limited simulation.","section":"Section 2.3"},{"comment":"The paper states that 'similar changes in pebble mass occur during other outbursts' and that the mass decreases 'by a factor of ≈2,' but no statistics over the several dozen outbursts in each model are presented. A median and a measure of scatter (or a histogram of peak-to-trough pebble-mass ratios) would establish that the factor of two is representative, not just a property of the single highlighted event at 244.5 kyr in M1 and the ~340 kyr event in M2.","section":"Section 3.1 (Fig. 4)"}],"minor_comments":[{"comment":"The caption should state explicitly which curve corresponds to the pebble surface density Σ_peb from Eq. (14); the text alternates between 'total mass of refractory components' and 'pebble mass,' and the two are not obviously the same quantity.","section":"Fig. 4"},{"comment":"Both models have high disk-to-star mass ratios (0.55 and 0.60) and are gravitationally unstable; the conclusions should either be restricted to this regime or accompanied by a justification for extrapolating to lower-mass, less gravitationally unstable disks.","section":"Section 2.5"},{"comment":"The adopted v_frag values are quoted to one significant figure without uncertainty; since the fragmentation barrier scales as v_frag^2, a short sensitivity statement would help the reader judge how robust the factor-of-two result is to plausible variations of v_frag for icy versus bare grains.","section":"Section 2.4"},{"comment":"The sentence saying 'the first considered luminosity outburst profile and corresponding masses of refractory components are shown in right panels of Fig. 4' should be corrected to match the figure caption, which places the refractory-component curve in the middle panel.","section":"Section 3.1"},{"comment":"The citation 'Molyarova et al. (2021, see)' in Section 2.4 should point to a specific section or equation if it is intended to justify the adopted v_frag dependence on ice mantles.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for Astronomy Reports and the modeling appears internally consistent. The main issue is that the headline claim is more assertive than the demonstrated scenario: the abstract and conclusions present the factor-of-two pebble reduction as a general outcome, while the paper itself identifies a plausible alternative destruction mechanism that would weaken it. I recommend asking the authors to either add an erosion-limited run or rephrase the abstract and conclusions to make the conditional nature explicit. I see no concerns about originality or misconduct."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper you asked about is a solid, clearly written simulation study from the FEOSAD group. The genuinely new piece is that the outbursts are now self-consistent MRI-driven events rather than manually imposed ones, and the analysis targets pebbles and their ice mantles specifically. The 2D spiral-structure desorption pattern (Fig. 6) and the viscous-versus-irradiation heating dichotomy for snowline sensitivity are the most interesting outcomes. I think those are real contributions.\n\nThe central quantitative claim—that an outburst halves the disk's pebble mass—is an emergent result within the model, but it rests on a specific microphysical choice: when water ice evaporates, the fragmentation barrier drops from 5 to 0.5 m/s and all oversize aggregates are immediately recycled to small dust. The authors say this in Section 2.3 and even cite the competing erosion-limited picture (Stammler & Birnstiel 2022) and Houge et al. (2024), which favors prolonged destruction in V883 Ori. They give a reasonable motivation for their choice (spiral structure keeps aggregates mantle-stuck), but they don't test the alternative. So the factor-of-two is a scenario-dependent number, not a robust prediction. The paper is transparent about this, which is to its credit, but the abstract and Section 4 bullet state the factor-of-two without that qualification. That is the main soft spot.\n\nThe snowline result—H2O snowline moves less because it sits in the viscous-heating-dominated zone—is plausible and well argued. It also explains why their model differs from V883 Ori, where the water snowline shifts a lot. So they're not overclaiming; they put the difference in a useful physical context.\n\nMinor points: only two model masses, no code/data release, and detailed analysis of one outburst per model, though Fig. 4 shows the behavior is recurrent. None of these are problems at this stage.\n\nOverall: a worthwhile paper, honestly framed, with a clear new simulation capability. The recommendation: send it to peer review, but ask the authors to either add an erosion run or explicitly frame the factor-of-two as the instantaneous-destruction scenario. I would bring it to reading group.\n\nBest.","headline":"Self-consistent MRI outbursts in FEOSAD halve pebble mass under the instantaneous-destruction scenario; the 2D spiral desorption pattern is the real new result.","tokens_in":14072,"tokens_out":2493,"would_cite":true,"duration_ms":25282,"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":"Luminosity outbursts can cut a protoplanetary disk's total pebble mass by about half, simulations show.","keywords":["protoplanetary disks","pebbles","luminosity outbursts","ice mantles","snowlines","dust fragmentation","FU Orionis","magnetorotational instability"],"falsifier":"Run the same disk model with the erosion-limited destruction scheme of Stammler & Birnstiel (2022) and compare the total pebble mass through an outburst: if the factor-of-two drop does not appear, the instantaneous-destruction assumption is the deciding ingredient. Observational check: monitor a freshly outbursting FUor object with millimeter-wavelength spectral-index mapping across the snowline region; if the largest grains persist through the outburst, gradual erosion is favored.","tokens_in":12962,"feed_emoji":"🪐","tokens_out":7498,"duration_ms":68507,"temperature":0.7,"pith_summary":"This paper argues that the flare-like luminosity outbursts experienced by young stars in their first few hundred thousand years can cut the total mass of pebbles in the surrounding disk by about a factor of two. The mechanism is not the direct heating itself but the loss of ice mantles: when water ice evaporates from dust aggregates, the grains lose the glue that let them survive collisions, fragment into monomers, and the fragmentation barrier drops. The paper uses global 2D simulations of a self-gravitating disk, including dust growth, drift, and the adsorption/desorption of H2O, CO2, CH4, and CO, to trace pebbles throughout disk formation and early evolution. It also finds that the water snowline is relatively stable because it sits in the viscously heated inner region, while the other snowlines move much more; mantle desorption is a two-dimensional, spiral-shaped phenomenon. The result matters because pebbles are the pipeline to planetesimals and planetary cores, so outbursts may temporarily choke that pipeline.","feed_headline":"Outbursts can halve a young disk's pebble supply","feed_subtitle":"Simulations trace how evaporated ice mantles shatter dust aggregates, leaving disks pebble-starved for millennia.","key_machinery":"The mechanism carrying the argument is the coupling between ice mantle chemistry and the dust fragmentation barrier. The code sets a maximum grown-dust size $a_{\\rm frag}$ that depends on the fragmentation velocity $v_{\\rm frag}$; ice-mantled grains have $v_{\\rm frag}=5\\,\\mathrm{m\\,s^{-1}}$, bare grains have $v_{\\rm frag}=0.5\\,\\mathrm{m\\,s^{-1}}$. When a luminosity outburst heats the disk and evaporates mantles, $v_{\\rm frag}$ drops by an order of magnitude, the fragmentation barrier $a_{\\rm frag}$ changes abruptly, and the model recycles all dust exceeding the new barrier into small dust. The outbursts themselves are produced self-consistently in the FEOSAD thin-disk hydrodynamics code through a layered MRI-effective $\\alpha$ parameterization that creates episodic accretion luminosity. Pebbles are defined as grown dust with Stokes number $St>0.01$ and size above 0.05 cm, and their ice masses are tracked via adsorption and desorption of H2O, CO2, CH4, and CO.","core_discovery":"The central claim is that a typical MRI-driven luminosity outburst in a young self-gravitating disk reduces the total mass of pebbles by roughly a factor of two, by thermally desorbing water ice mantles and thereby destroying the aggregates that the mantles held together. In the FEOSAD simulations, the disk loses about half its pebble mass during an outburst of roughly two hundred years; pebble growth then rebuilds the population over several thousand years, much longer than the days-to-decades freeze-out timescales of the volatiles. The paper also finds that snowline shifts are asymmetric: the H2O snowline barely moves because it sits in the viscously heated region inside about 7 au, while CO2, CH4, and CO snowlines, which lie in irradiation-dominated regions, move outward strongly. Mantle desorption takes place in a non-axisymmetric spiral pattern, so the pebble mass loss is not simply a radial band but a patterned 2D depletion. These conclusions hold in both a low-mass (0.66 solar mass) and a higher-mass (1.0 solar mass) disk model.","pith_inferences":["If real dust destruction is gradual (erosion) rather than instantaneous, the factor-of-two pebble loss is likely an upper bound; the same outburst could leave more pebbles if aggregates shed mass slowly instead of shattering.","With outbursts repeating on roughly 10^5-year timescales and recovery taking roughly 10^3 years, early disks may spend a nontrivial fraction of their lives in a pebble-depleted state, possibly modulating the efficiency of streaming instability and pebble accretion.","The paper's own planned bidisperse erosion model would directly test whether the factor-of-two result survives; the instantaneous-destruction assumption is the main uncertainty it flags.","The distinction between viscous-heating and irradiation-heating snowline regions implies that young, massive disks should show weak water-snowline shifts during outbursts, whereas older, less massive disks like V883 Ori should show strong shifts; this could be used to age- or mass-date outbursting disks."],"forward_implications":["During an outburst, pebble mass drops by about half and recovers only after several thousand years of coagulation, so the pebble population is depressed for much longer than the roughly 200-year outburst itself.","Because water ice is the glue, any process that removes water ice from grains, whether radial drift across the snowline or outburst heating, will suppress pebble numbers in the affected region.","Snowline responses to outbursts are not universal: molecules whose snowlines sit in irradiation-heated zones (CO2, CH4, CO) move far more than water's, which sits in the viscously heated zone; interpreting observed snowline shifts therefore requires knowing the local heating regime.","The non-axisymmetric, spiral-shaped desorption pattern means that spatially unresolved disk averaging can hide strong local pebble destruction; resolved observations of an outbursting disk should see azimuthal structure in ice tracers and dust size."],"supporting_citations":[{"why":"Supplies the FEOSAD 2D thin-disk hydrodynamic model that produces the simulated disks and their self-consistent outbursts.","marker":"Vorobyov et al. (2018)"},{"why":"Provides the volatile adsorption/desorption and ice-dependent fragmentation velocity scheme used to track H2O, CO2, CH4, and CO on dust.","marker":"Molyarova et al. (2021)"},{"why":"Establishes the pebble identification criterion and the initial pebble/ice properties on which this work builds.","marker":"Topchieva et al. (2024)"},{"why":"Gives the fragmentation-barrier formula $a_{\\rm frag}$ that sets the maximum grain size and drives grain destruction when the barrier drops.","marker":"Birnstiel et al. (2012)"},{"why":"Supplies the 'many seeds' picture of aggregates held together by ice mantles, the physical justification for instantaneous breakup into monomers.","marker":"Schoonenberg et al. (2017)"},{"why":"Shows idealized-outburst snowline shifts and post-outburst dust-growth timescales, the basis for the recovery timescales reported here.","marker":"Vorobyov et al. (2022)"},{"why":"Provides the layered MRI effective-alpha parameterization that produces the episodic accretion luminosity outbursts in the model.","marker":"Kadam et al. (2022)"},{"why":"Supplies the alternative erosion treatment of dust destruction, used as the contrasting scenario that would delay the pebble loss.","marker":"Stammler & Birnstiel (2022)"},{"why":"Presents the V883 Ori spectral-index comparison that favors prolonged rather than instantaneous dust destruction, motivating the paper's discussion.","marker":"Houge et al. (2024)"}],"fun_headline_variants":["Outbursts reduce pebble abundance by half in disks","Ice mantle loss during flares shatters dust aggregates","Simulations show pebble mass drops 50% after outburst","Outbursts trigger 2D spiral pebble depletion in disks","Pebble recovery after outbursts takes millennia, not years"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on the assumption that dust aggregates shatter into monomers the moment their ice mantles evaporate; if destruction instead proceeds gradually through collisions, the pebble-mass drop would be smaller or delayed.","fun_headline_variants_meta":{"raw":{"variants":["Outbursts reduce pebble abundance by half in disks","Ice mantle loss during flares shatters dust aggregates","Simulations show pebble mass drops 50% after outburst","Outbursts trigger 2D spiral pebble depletion in disks","Pebble recovery after outbursts takes millennia, not years"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000391,"raw_usage":{"total_tokens":2145,"prompt_tokens":1123,"completion_tokens":1022,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":739,"completion_tokens_details":{"reasoning_tokens":939}},"tokens_in":739,"tokens_out":1022,"duration_ms":9825,"temperature":1.0,"reasoning_tokens":939,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:09:14.198938+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same disk model with the erosion-limited destruction scheme of Stammler & Birnstiel (2022) and compare the total pebble mass through an outburst: if the factor-of-two drop does not appear, the instantaneous-destruction assumption is the deciding ingredient. Observational check: monitor a freshly outbursting FUor object with millimeter-wavelength spectral-index mapping across the snowline region; if the largest grains persist through the outburst, gradual erosion is favored.","supporting_citations":[],"review_version":1}