{"id":"56628a61-c5b1-41fe-9905-fa7b42b3a225","arxiv_id":"1908.02515","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Simulations of layered protoplanetary disks find that rings form in the inner dead zone, migrate inward, and trigger FUor-like accretion events.","lead":"A set of computer simulations shows that gaseous rings, dense and low-viscosity bands, naturally form in the inner dead zone of forming protoplanetary disks and migrate inward before breaking in bright outbursts. The result suggests new ways to trap dust and build planets close to young stars.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'ubiquitous ring' result is tested only at α_rd=10^-5, the lowest end of the residual dead-zone viscosity range cited in the paper; higher residual viscosity (10^-4–10^-3) is never explored.","rationale":"The central claim—layered disks form dense, long-lived gaseous rings—is well supported within the model: the fiducial comparison is clean, the spacetime diagrams show persistent axisymmetric structures, and the viscous-torque sign changes across the rings are consistent with the proposed mechanism. The most load-bearing concern is external rather than internal: ring strength is set by the residual dead-zone viscosity α_rd, and the paper's own text allows this parameter to be as large as 10^-3–10^-5 while all simulations use 10^-5. Since α_eff in the ring is approximately α_rd when Σ ≫ Σ_a, a residual viscosity of 10^-4 or 10^-3 would reduce the viscosity contrast that drives the rings from roughly 10^3 to 10^2 or 10, potentially erasing the reported two-order-of-magnitude ring contrast. The parameter study varies Tcrit and Σ_a but not α_rd, so the 'ubiquitous' claim is not yet bracketed by the model's own uncertainty range. This is a missing sensitivity test, not an internal contradiction. The reader's weakest_assumption already flagged smoother transitions or higher residual viscosity; my emphasis is specifically on the untested residual-viscosity value, which is the most direct control on ring amplitude. The conditional verdict remains appropriate, pending a rerun with α_rd = 10^-4 and 10^-3.","tokens_in":24729,"tokens_out":13539,"duration_ms":149551,"concrete_test":"Rerun the fiducial model1 T1300 S100 with α_rd = 10^-4 and α_rd = 10^-3, keeping all other settings identical, and measure the main ring's surface-density contrast Σ_ring/Σ_av and the number of rings in spacetime diagrams over 0.3–0.5 Myr. If rings persist with contrast ≳2 in both runs, the central claim is robust to the residual-viscosity uncertainty; if rings weaken or disappear at α_rd = 10^-3, the paper's conclusion is conditional on the lowest end of its own adopted parameter range.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that magnetically layered disks ubiquitously form high-contrast gaseous rings rests on the viscosity contrast between the MRI-active surface (α_a=0.01) and the dead midplane. In Equations (5)–(8), the dead-zone residual viscosity is set through α_rd = min(10^-5, α_a Σ_a/Σ_d), and every layered simulation fixes α_rd = 10^-5. The paper itself states in Section 2.1 that residual dead-layer viscosity can be as large as 10^-3–10^-5, citing Okuzumi & Hirose (2011). Choosing 10^-5 selects the most favorable end of that range. Because α_eff in the ring is approximately α_rd when Σ ≫ Σ_a, the viscosity contrast that drives the ring-forming torque and the reported two-order-of-magnitude surface-density excess shrinks from about 10^3 to about 10^2 or 10 if α_rd = 10^-4 or 10^-3 are adopted. Section 3.4 varies Tcrit and Σ_a but never α_rd, so the 'ubiquitous rings' conclusion is not bracketed by the paper's own stated uncertainty range. This is not an internal inconsistency; it is a missing sensitivity test on the parameter most directly controlling ring amplitude, and it leaves the astrophysical claim conditional on an extreme value.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents global thin-disk hydrodynamic simulations of protoplanetary disk formation starting from collapsing cloud cores, with a magnetically layered disk model implemented as an adaptive α_eff (Equations 5–8). The authors compare a fiducial layered model (model1 T1300 S100) with a constant-α fully MRI-active model, and explore variations in Tcrit, Σ_a, and core mass (10 models total). The main reported result is that the layered models develop a dead zone (inner ~15 au in the fiducial case) inside which long-lived, axisymmetric, high-surface-density gaseous rings form via viscous torques at sharp viscosity transitions. The rings migrate inward, occasionally become gravitationally unstable, and their MRI-triggered destruction bursts resemble FUor outbursts. The paper also argues that pressure maxima in the rings can trap dust and increase the fragmentation barrier.","tokens_in":25019,"tokens_out":5482,"duration_ms":53423,"significance":"If the result holds, it provides a plausible and observationally relevant mechanism for gaseous ring formation in the innermost few au of protoplanetary disks, linking layered accretion, episodic outbursts, and dust trapping. The paper's strengths are its self-consistent collapse initial conditions, the smallest inner sink cell (0.4 au) used in such global simulations, the explicit comparison between layered and constant-α disks, and the candid discussion of model limitations (dust evolution, vortex resolution, sparse event statistics). The ring formation mechanism is not imposed by hand; it emerges from the viscous evolution of the layered-α prescription, and no parameter is fitted to reproduce rings. However, the central claim of ubiquitous rings rests on the assumed low residual dead-zone viscosity (α_rd = 10^-5), which is the most favorable value in the range cited by the paper, and the model dependence on that parameter is not explored.","major_comments":[{"comment":"All layered-disk models set α_rd = 10^-5, which is the lowest (most favorable) end of the residual-viscosity range cited in Section 2.1 (10^-3 to 10^-5, based on Okuzumi & Hirose 2011). Since α_eff ≈ α_rd in the ring interiors where Σ ≫ Σ_a, and since the ring-forming torque in Section 3.2 is said to scale with the gradient of kinematic viscosity, the reported two-order-of-magnitude contrast in Σ and α_eff is expected to shrink to about one order of magnitude or less if α_rd = 10^-4 or 10^-3, respectively. Section 3.4 varies Tcrit and Σ_a but never α_rd, so the claim that rings form ubiquitously in layered disks is not bracketed by the paper's own stated uncertainty range. I request either simulations with α_rd = 10^-4 and 10^-3, or a quantitative argument (e.g., a scaling analysis of the ring amplitude with α_a/α_rd) explaining why the conclusion survives over the stated range.","section":"§2.1, Eq. (8); §3.4"},{"comment":"The dead zone is defined as the region where α_eff falls below 80% of α_a = 0.01, with the justification that 'below this threshold ... surface density started to diverge.' This threshold is arbitrary, and the quantitative dead-zone extents quoted in the paper (≈16 au for the fiducial model, ≈78 au for Σ_a = 10 g cm^-2, ≈3.3 au for the low-mass model; Section 3.4 and Figures 9–12) will shift if the threshold is changed. Since these extents are used as the basis for parameter-dependence conclusions (e.g., 'five times larger' and 'five times smaller'), the definition should be made more robust—for example, by quoting the radius where the surface-density profile first diverges from the constant-α model, or by showing the sensitivity of the extent to the chosen threshold.","section":"§3.1, Figure 3"},{"comment":"The migration rate of about −25 au/Myr and the mean outburst interval of about 38,000 yr are derived from only two ring-discontinuity events (at approximately 0.305 and 0.375 Myr), one of which was not directly sampled in the output (the text notes that the 0.305 Myr event was not caught). The first panel of Figure 7 also shows that the ring's trajectory is highly nonlinear. The reported 'average' quantities therefore carry a large, unquantified uncertainty; the text should either provide an error estimate, extend the time baseline with more events, or explicitly soften the quantitative claims in the abstract and conclusion.","section":"§3.2, Figure 7"}],"minor_comments":[{"comment":"The model shown in green for the low-mass case is called 'model2 T1500 S100' in the text and figure, but Table 1 and the surrounding text refer to 'model2 T1300 S100' (which is also the subject of Figure 12); please harmonize the model naming.","section":"§3.4, Figure 9"},{"comment":"The abstract states 'the smallest possible inner computational boundary,' but the text more precisely says 'the smallest sink cell used in global collapse simulations'; please rephrase the abstract to avoid overstatement.","section":"§2.2"},{"comment":"The viscous stress tensor uses 'e' for the unit tensor, which can be confused with the internal energy density e in Eq. (3); please use a distinct symbol or explicitly note the difference.","section":"§2.1, Eq. (4)"},{"comment":"The vertical lines marking the two rings at 1.2 and 3.2 au are described in the text but are not labeled in the figure captions; adding explicit labels or a legend would improve readability.","section":"§3.1, Figures 3 and 8"},{"comment":"The explanation that the inner dead-zone boundary undergoes a sharp transition while the outer boundary is smooth is intuitive but would benefit from an explicit expression for the ∝ −dν/dr component of the viscous torque, which is invoked but not written out.","section":"§3.2"},{"comment":"The dust-trapping analysis uses a single snapshot, and the text already acknowledges that the rings have short lifetimes; the abstract's statement that dust 'could be trapped' should be qualified as 'while the rings exist,' since the ring lifetime may be much shorter than the dust growth timescale.","section":"§3.3"},{"comment":"The word 'ubiquitously' is used for the ring-formation claim, but the conclusion rests on 8 layered-disk models with a specific, narrow parameter choice; consider replacing 'ubiquitously' with 'in all models explored here' to match the actual scope.","section":"Abstract and Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest and technically solid within its assumptions, and its main limitation—the absence of a sensitivity study on α_rd—is a natural extension rather than a fundamental flaw. However, because that parameter directly controls the ring amplitude and the 'ubiquitous' claim, I would request the additional simulations or a quantitative scaling argument before publication. The 'ubiquitous' language in the abstract is somewhat strong for a parameter space of 8 layered runs, but this is fixable by rewording."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth reading and worth refereeing, but the central claim comes with a caveat the paper does not test. The simulations are a real step up from earlier dead-zone ring studies: they follow collapse from a prestellar core to the T Tauri phase in a thin-disk global model with a 0.4 au inner boundary, and they show that a layered viscosity prescription produces long-lived, axisymmetric rings in the inner few au while a constant-alpha run does not. The torque analysis makes the ring formation mechanism concrete: positive feedback between viscous torques and the sharp viscosity transition piles gas into rings. The parameter study of Tcrit and Sigma_a is sensible, and the dust-trapping implications are presented as plausible rather than proven, with the missing dust evolution explicitly acknowledged.\n\nThe soft spot is alpha_rd. Equation (8) sets the dead-zone residual viscosity to 1e-5, the lower bound of the 1e-3–1e-5 range the paper itself cites from Okuzumi & Hirose. The ring's effective alpha is essentially alpha_rd once the column density exceeds Sigma_a, so the viscosity contrast that builds the rings is about 1e3 with alpha_rd = 1e-5. At 1e-4 or 1e-3 that contrast drops to 1e2 or 10, and the distinct rings would likely weaken or disappear. The authors never run a sensitivity test on this parameter, which directly controls the amplitude of the phenomenon they are claiming. That is the difference between \"ubiquitous rings\" and \"rings under an optimistic residual viscosity.\" The migration rate and outburst period also come from only two events, one of which was missed in the output cadence, so the quantitative statistics are thin. The vortices are admittedly marginally resolved. None of these are fatal; the main simulation result is clear and the qualitative behavior is probably robust. But the alpha_rd point is a genuine gap.\n\nFor a colleague working on dead zones, episodic accretion, or inner-disk planetesimal formation, this is a useful paper to know. It deserves a serious referee. I would recommend sending it out with the request that the authors either run an alpha_rd sensitivity sequence covering 1e-4 and 1e-3 or provide a stronger physical justification for fixing it at the extreme value. Without that, the paper is a good conditional result, not a settled one.","headline":"Solid global simulations of dead-zone ring formation with a clean fiducial comparison, but the headline result rides on the least favorable residual viscosity setting and the authors never bracket it.","tokens_in":25638,"tokens_out":2067,"would_cite":true,"duration_ms":25226,"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":"Dead zones in layered protoplanetary disks are not passive bottlenecks but dynamical factories of dense gaseous rings that trap dust, migrate inward, and erupt in FUor-like accretion events.","keywords":["protoplanetary disks","star formation","T Tauri stars","dead zones","magnetorotational instability","gaseous rings","episodic accretion","numerical hydrodynamics"],"falsifier":"Run the same simulation with the MRI-on switch smoothed over a finite temperature or density range and with the residual dead-zone viscosity raised from $10^{-5}$ to $10^{-3}$: if the dense axisymmetric rings no longer form, the sharp viscosity transition is the load-bearing element. Observational complement: sub-au millimeter continuum imaging of a nearby T Tauri disk that resolves inside a few au should show a compact bright ring or gap with the predicted dust-trap contrast if the rings are generic; a smooth, featureless inner disk would argue against the mechanism.","tokens_in":24519,"feed_emoji":"🪐","tokens_out":8936,"duration_ms":90813,"temperature":0.7,"pith_summary":"This paper argues that the magnetically layered structure long assumed for protoplanetary disks is not a quiet bottleneck but a persistent source of dense gaseous rings. Global thin-disk simulations that start from collapsing cloud cores show that a dead zone develops inside roughly 15 au, within which long-lived, axisymmetric rings form wherever viscous torques pile gas against a sharp viscosity transition. The rings are dense, low-viscosity, marginally gravitationally unstable structures that migrate inward and end in MRI-triggered accretion events resembling FUor outbursts. Because the rings create pressure maxima with dust fragmentation barriers up to about 100 m, the paper concludes they can trap dust and provide rapid sites for planetesimal formation, linking disk ionization structure to planet formation.","feed_headline":"Layered disks grow dense gas rings that trap dust and erupt","feed_subtitle":"These rings trap dust at centimeter-to-meter sizes, giving planet formation a head start.","key_machinery":"The central object is the effective viscosity parameter $\\alpha_{\\rm eff} = (\\Sigma_a \\alpha_a + \\Sigma_d \\alpha_d)/\\Sigma$, with $\\alpha_a = 0.01$, residual dead-zone viscosity $\\alpha_d = \\min(10^{-5},\\ \\alpha_a \\Sigma_a/\\Sigma_d)$, and thermal MRI activation when the midplane temperature exceeds $T_{\\rm crit}$. The key mechanism is the viscous torque, which has a component proportional to the negative gradient of the kinematic viscosity; at the sharp inner edge of the dead zone this torque piles up gas instead of letting it flow inward. The pileup increases surface density and lowers $\\alpha_{\\rm eff}$, which steepens the viscosity gradient and strengthens the pileup, producing a positive feedback that builds and maintains the rings.","core_discovery":"The paper's central claim is that the dead zone of a magnetically layered protoplanetary disk is not a uniform, quiescent region but a self-organizing system of gaseous rings. With canonical parameters (MRI-active surface layer column density $\\Sigma_a = 100$ g cm$^{-2}$ and activation temperature $T_{\\rm crit}=1300$ K), the dead zone extends to about 15 au, and the inner few au host multiple axisymmetric rings with surface densities up to two orders of magnitude higher and effective viscosities about two orders of magnitude lower than a fully MRI-active disk. The rings form through a positive feedback: viscous torques with a component proportional to the negative gradient of kinematic viscosity push gas into low-viscosity regions, which raises surface density and lowers $\\alpha_{\\rm eff}$, steepening the gradient further. The rings migrate inward, become gravitationally unstable, generate spiral waves and short-lived vortices, and terminate when the midplane temperature crosses $T_{\\rm crit}$, triggering an MRI event that rapidly accretes the ring onto the star in an FUor-like outburst. Within the rings, dust faces fragmentation barriers of centimeters to meters, occasionally approaching 100 m, in contrast to a few millimeters in a fully MRI-active disk.","pith_inferences":["The paper does not say this, but the same positive-feedback mechanism should operate at any sharp radial viscosity drop in a Keplerian disk, so opacity or dust-sublimation transitions could produce similar rings even without MRI layering.","An outside inference from the ring lifetimes: gas-only runs give ring lifetimes of only a few percent of the roughly 1 Myr viscous timescale before MRI erupts, so whether dust actually reaches planetesimal sizes depends on whether growth beats about $10^4$ yr, a question the paper leaves for future dust evolution models.","If the ring phase shortens with decreasing stellar mass, as the lower-core-mass run suggests, then low-mass stars should show both weaker inner-disk dust trapping and less frequent FUor-like accretion, a statistical prediction that young stellar cluster surveys could test.","A direct numerical check of the mechanism would replace the step-function MRI activation with a smooth, microphysically motivated transition and raise the residual dead-zone viscosity from $10^{-5}$ to $10^{-3}$; if the dense rings disappear, the sharp viscosity jump is what carries the claim."],"forward_implications":["Layered protoplanetary disks generically develop long-lived, axisymmetric gaseous rings inside a few au, whereas fully MRI-active disks with constant $\\alpha$ do not.","The rings migrate inward at tens of au per Myr and end in MRI-triggered accretion episodes, providing a mechanistic route to FUor-like outbursts with intervals of roughly tens of thousands of years.","Ring pressure maxima raise the local dust fragmentation barrier from millimeter sizes to centimeter-meter sizes, occasionally reaching about 100 m, making the inner few au a favorable site for dust trapping and rapid planetesimal growth.","The inner-disk structure is highly sensitive to active-layer thickness and initial cloud core mass: a tenfold thinner active layer widens the dead zone to roughly 78 au and slows ring evolution, while a lower-mass core yields a short-lived ring phase without MRI-triggered outbursts.","Because large bodies decouple from the gas, planetesimals that form in the rings may survive the MRI-triggered dispersal of the gaseous ring, so the rings can assist planet formation despite their short gas lifetime."],"supporting_citations":[{"why":"Introduces the magnetically layered disk model with a dead zone, the physical setup the simulations implement.","marker":"Gammie (1996)"},{"why":"Supplies the alpha-viscosity prescription that the effective alpha parameter adapts to layered structure.","marker":"Shakura, & Sunyaev (1973)"},{"why":"Provides the adaptive $\\alpha_{\\rm eff}$ formula and MRI activation temperature used to emulate the layered disk.","marker":"Bae et al. (2014)"},{"why":"Justifies the nonzero residual dead-zone viscosity driven by active-layer turbulence, used in the alpha_d prescription.","marker":"Okuzumi, & Hirose (2011)"},{"why":"Gives the canonical cosmic-ray ionized column density $\\Sigma_a = 100$ g cm$^{-2}$ that sets the active layer thickness.","marker":"Umebayashi, & Nakano (1981)"},{"why":"Provides the thin-disk hydrodynamics framework and viscous torque formalism in which the rings are found.","marker":"Vorobyov, & Basu (2009)"},{"why":"Establishes magnetorotational instability as the turbulence mechanism whose suppression defines the dead zone.","marker":"Hawley et al. (1995)"},{"why":"Supplies the layered-disk episodic accretion framework and $T_{\\rm crit}$ values against which the ring-driven outbursts are compared.","marker":"Zhu et al. (2010)"}],"fun_headline_variants":["Dead zones spawn migrating gas rings that trap dust","Gas rings in disk dead zones funnel dust to planet birth","Simulations show dead zones build dust-trapping ring arrays","Layered disks churn out rings that seed planet formation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The rings grow from a sharp jump in gas viscosity between the magnetically active surface layers and the dead midplane, so if real dead zones have smoother transitions or a higher residual viscosity, the rings would be weaker or absent.","fun_headline_variants_meta":{"raw":{"variants":["Dead zones spawn migrating gas rings that trap dust","Gas rings in disk dead zones funnel dust to planet birth","Simulations show dead zones build dust-trapping ring arrays","Layered disks churn out rings that seed planet formation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000414,"raw_usage":{"total_tokens":2195,"prompt_tokens":1058,"completion_tokens":1137,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":674,"completion_tokens_details":{"reasoning_tokens":1072}},"tokens_in":674,"tokens_out":1137,"duration_ms":11120,"temperature":1.0,"reasoning_tokens":1072,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:41:41.438347+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same simulation with the MRI-on switch smoothed over a finite temperature or density range and with the residual dead-zone viscosity raised from $10^{-5}$ to $10^{-3}$: if the dense axisymmetric rings no longer form, the sharp viscosity transition is the load-bearing element. Observational complement: sub-au millimeter continuum imaging of a nearby T Tauri disk that resolves inside a few au should show a compact bright ring or gap with the predicted dust-trap contrast if the rings are generic; a smooth, featureless inner disk would argue against the mechanism.","supporting_citations":[{"cited_title":"I., & Sunyaev, R","cited_arxiv_id":null,"evidence_quote":"Supplies the alpha-viscosity prescription that the effective alpha parameter adapts to layered structure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the adaptive $\\alpha_{\\rm eff}$ formula and MRI activation temperature used to emulate the layered disk."},{"cited_title":"2011, ApJ, 742, 65 Padovani, M., & Galli, D","cited_arxiv_id":null,"evidence_quote":"Justifies the nonzero residual dead-zone viscosity driven by active-layer turbulence, used in the alpha_d prescription."},{"cited_title":"1981, Publications of the Astronomical Society of Japan, 33,","cited_arxiv_id":null,"evidence_quote":"Gives the canonical cosmic-ray ionized column density $\\Sigma_a = 100$ g cm$^{-2}$ that sets the active layer thickness."},{"cited_title":"I., & Basu, S","cited_arxiv_id":null,"evidence_quote":"Provides the thin-disk hydrodynamics framework and viscous torque formalism in which the rings are found."},{"cited_title":"F., Gammie, C","cited_arxiv_id":null,"evidence_quote":"Establishes magnetorotational instability as the turbulence mechanism whose suppression defines the dead zone."},{"cited_title":"F., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the layered-disk episodic accretion framework and $T_{\\rm crit}$ values against which the ring-driven outbursts are compared."}],"review_version":1}