{"id":"3acbf822-829d-4db6-9b27-3cd7f9eb3496","arxiv_id":"2502.04452","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Fragile dust destroyed beyond a dead zone limits mm-dust emission to the dead zone radius, naturally producing the compact discs seen around young stars.","lead":"This paper proposes that compact protoplanetary discs form when fragile dust is shattered by turbulence in the outer disc, so the disc size is set by the quiet dead zone. It uses dust evolution simulations and synthetic ALMA images to argue this mechanism matches observed compact discs.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pre-relaxing gas without dust removes the dead-zone-edge pressure bump that would otherwise trap dust, so the compact-disc result may be an artifact of the initial conditions.","rationale":"The reader's weakest assumption (sharp alpha transition and fragile dust) is related to my concern, but the more specific issue is the pre-relaxation step that removes the pressure bump at the dead-zone edge. This is a modeling choice that may predetermine the outcome: without a dust trap, mm dust cannot survive in the MRI-active region, so the disc size naturally becomes the dead-zone radius. The paper's own Fig. 1 shows that without relaxation a positive pressure gradient develops at the transition, which is known to trap dust. Since dust is present in real discs during the viscous adjustment, the initial condition is physically questionable. This is a concrete, testable concern, not an unsupported objection. The reader's verdict of conditional acceptance remains appropriate, but the condition should include a co-evolution test from a realistic initial state, not just variation of alpha-transition shape and vfrag.","tokens_in":27720,"tokens_out":2718,"duration_ms":30656,"concrete_test":"Re-run Model 2 (sharp transition, Rt=30 au) without the 6-Myr gas-only relaxation: initialize the self-similar gas disc and the MRN dust distribution simultaneously, and evolve both together. At t=1 Myr, compare the mm-dust surface density profile and the measured observed radius Rd,90 at 1.3 mm. If a bright ring at ~Rt forms and Rd,90 is significantly larger than the 26.67 au reported in Table 2, the compact-disc mechanism is not supported; if the disc remains compact despite the initial pressure bump, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that fragmentation in the MRI-active zone sets the disc size to the dead-zone radius. This requires that no long-lived dust trap exists at the dead-zone edge; otherwise mm dust would be retained there and the disc would not be compact. The paper intentionally removes this trap: Section 4.1 evolves the gas-only disc for 6 Myr so that the pressure bump at the sharp alpha-transition edge relaxes before dust is introduced. Previous work (e.g., Pinilla et al. 2016) shows this bump traps dust and can persist for several Myr. In a real disc, dust is present throughout the viscous relaxation and would encounter the bump. By initializing dust only after the bump has dissipated, the model selects the regime where it is absent. The paper does not test co-evolution of gas and dust from the initially self-similar state. Therefore 'fragmentation beyond the dead zone produces compact discs' is not robust to the choice of initial condition; it may be an artifact of pre-relaxing the gas without dust.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that fragmentation of fragile dust (vfrag ~ 1 m/s) in moderately turbulent MRI-active regions outside a dead zone can explain radially compact protoplanetary discs. Using 1D DustPy models with a radially varying alpha_SS profile, the authors first evolve the gas alone for 6 Myr to remove the dead-zone-edge pressure bump, then introduce dust and follow its coagulation, fragmentation, drift, and emission. They find that the mm-dust disc is truncated near the dead-zone transition radius, yielding Rd,90 values of roughly 19--38 au, wavelength-independent sizes for sharp transitions, and a size--luminosity relation broadly consistent with observed compact discs. The paper also explores planet-carved pressure bumps, porous dust, higher fragmentation velocities, and smooth versus sharp alpha transitions, and compares synthetic ALMA images with objects such as Sz 66 and MP Mus.","tokens_in":27952,"tokens_out":10037,"duration_ms":106790,"significance":"If the mechanism survives the model assumptions, it is significant: it gives a physical interpretation of compact discs as tracing the dead-zone radius, and it produces testable predictions (similar 1.3/3 mm sizes for sharp transitions, a weak outer-edge 'shoulder' or ring, and a degeneracy between alpha_SS and vfrag). The study is transparent and well-documented: the parameter space is clearly tabulated, no parameters are fitted to the target observations, and the synthetic observations use standard, reproducible tools (RADMC-3D and CASA). The main reservation is that the central result depends on three conditions --- a sharp alpha transition, vfrag ~ 1 m/s, and dust being absent during the gas relaxation phase --- and the paper does not fully stress-test the third condition, which is the most directly load-bearing for the title claim.","major_comments":[{"comment":"The central claim that the mm-dust disc size is set by the dead-zone radius depends on the absence of a dust trap at the dead-zone outer edge, but that absence is imposed by construction: the gas is evolved alone for 6 Myr before dust is introduced, so the pressure bump at Rt has relaxed before any dust can encounter it. In a real disc, dust is present throughout the viscous relaxation, and Pinilla et al. (2016) --- cited in this section --- show that the bump can trap dust for several Myr; the text itself says the trap 'remains as a boost to the pressure gradient ... for another a few million years'. The model therefore selects the regime in which the trap has already dissipated, and the reported Rd,90 values at 1 Myr would likely be larger in a co-evolution run that includes dust from t = 0. I recommend adding such a run, or providing a physical justification for why late dust introduction is appropriate. The 6-Myr gas-only preconditioning also means that the quoted 1 Myr and 3 Myr times are time since dust insertion, not the physical age of the disc, which complicates the comparison to young observed discs.","section":"Section 4.1 (initial condition)"},{"comment":"The smooth-transition model is presented as only mildly different, but Table 2 shows a substantial quantitative change: at 1 Myr, Model 4 has Rd,90 = 37.68 au and F1.3mm = 97.71 mJy, versus 26.67 au and 19.71 mJy for the sharp-transition Model 2; even at 3 Myr the sizes are 29.93 au versus 23.77 au. Because the slow-transition profile is the case motivated by non-ideal MHD simulations (Bai et al. 2016, cited in Section 4.2.1), the size-limiting mechanism is not established for realistic dead-zone edges. The text should either quantify this sensitivity in the conclusions or explicitly restrict the claim to sharp transitions and explain why that regime is physical.","section":"Section 4.2.1 and Table 2 (Model 4)"},{"comment":"The robustness test for the fragmentation velocity shows that mm-size dust begins to form beyond the dead zone for vfrag >= 5 m/s (Fig. 14), and the quoted M2',5 case gives Rd,90 = 29.46 au at 0.5 Myr, compared with 26.67 au for Model 2 at 1 Myr. The comparison uses different epochs, so the conclusion that 'less fragile dust does not significantly increase disc sizes' is not cleanly demonstrated. Since vfrag is a free parameter with considerable experimental uncertainty, the abstract and conclusions should state more carefully that the mechanism is contingent on vfrag ~ 1 m/s, and the analysis should specify how much dust outside Rt can be present before the observed size criterion changes.","section":"Section 5.3.2 and Fig. 14"}],"minor_comments":[{"comment":"There is a typo in the text: 'Sefan-Boltzmann constant' should read 'Stefan-Boltzmann constant'.","section":"Section 2.1"},{"comment":"The software list contains 'DustyPyLib'; the correct name used elsewhere in the paper is 'DustPyLib'.","section":"Data Availability"},{"comment":"The citation '(Anania et al, submitted)' appears in the text but has no entry in the reference list; please add the reference or remove the citation.","section":"Section 5.1"},{"comment":"Models M2',3 through M2',5 are listed in Table 4 without any resulting quantities; since the note explains why, it would be helpful to point the reader explicitly to Fig. 14 as the quantitative output for these runs.","section":"Table 4"},{"comment":"The M2' runs with high vfrag are quoted at 0.5 Myr, while all other size measurements in Table 2 are at 1 and 3 Myr; please state the reason for this different evaluation time in the text.","section":"Section 5.3.2"},{"comment":"The caption for Figure B1 says 'orange and green lines', but the text and legend identify the two profiles as 'noise-free' and 'thermal noise'; please make the colour and label terminology consistent.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"I agree with the conditional assessment. The paper is a competent and clearly written parameter study with a nice observational framing, and the stress-test concern about the 6-Myr gas-only relaxation is real and load-bearing. I would like to see a co-evolution run (dust present during the relaxation phase) before acceptance, together with a more honest statement about the sensitivity to smooth transitions and higher vfrag. The topic fits MNRAS well, and I see no novelty-disclosure concerns: the prior dead-zone-trap literature is properly cited."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here's my take on Tong & Alexander. The new idea is simple and worth knowing: if dust is as fragile as laboratory experiments suggest (v_frag ~ 1 m/s), the turbulent MRI-active zone beyond the dead zone prevents mm grains from growing, so the observed mm-dust disc truncates at the dead zone radius. That is a genuinely new formation channel for compact discs, distinct from drift-dominated or born-small scenarios. The 'hitchhiking' effect—a sharp drop in gas surface density at the dead zone edge locally accelerating the drift of ~mm dust and producing a low-contrast shoulder—is a nice piece of physics with a concrete, testable observational signature. The paper is also honest: Section 5.5 lists the missing MHD winds, back-reaction, and porosity issues. And no parameters are fit to observations; the size-luminosity comparison is post-hoc.\n\nBut the central claim is narrower than the abstract suggests. The mechanism requires fragile dust and a sharp alpha transition. The paper's own Model 4 (smooth transition) gives R_d,90 = 37.68 au versus 26.67 au for the sharp case, and v_frag >= 5 m/s allows mm dust beyond the dead zone. The robustness section shows the effect weakens, yet the conclusion still states that disc size is determined by the dead zone radius. That overstates the explored parameter space. Moderate issue, not fatal—the paper is transparent and the relevant regime (fragile ice grains, sharp transition) is plausible.\n\nThe bigger soft spot is the initial condition. The gas is pre-relaxed for 6 Myr without dust, which removes the pressure bump at the dead zone edge. Pinilla et al. 2016 shows that bump can trap dust for several Myr. The paper frames the relaxation as mimicking gradual dead zone formation, but it never tests co-evolution from the initially self-similar state. In a real disc, dust is present while the dead zone edge forms; if the bump persists, the outer disc retains mm dust and the compact-disc outcome is not realized. This is addressable—run one model with dust from the start—but it directly affects the core mechanism and should be acknowledged explicitly.\n\nMinor: data and code are not public. For a study whose main product is synthetic observations, that is a reproducibility weakness.\n\nOverall, a serious and well-executed model study with a testable new mechanism. I would bring it to reading group and would cite it in a compact-disc paper. It deserves peer review. Under review, I would ask for the co-evolution test and a more carefully worded conclusion that states the regime of validity.","headline":"Plausible new compact-disc mechanism via dead-zone fragmentation, but the central claim depends on a narrow parameter regime and the pre-relaxed gas initial condition needs testing.","tokens_in":28465,"tokens_out":2941,"would_cite":true,"duration_ms":30044,"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 can explain the prevalence of compact protoplanetary discs: when fragile dust fragments in the moderately turbulent region beyond the dead zone, the observed mm-sized dust disc is cut off at the dead zone edge, so disc size…","keywords":["protoplanetary discs","dead zones","dust fragmentation","dust evolution","magnetorotational instability","compact discs","synthetic observations","size-luminosity relation"],"falsifier":"A direct test would be to measure the 1.3-mm and 3-mm disc sizes of a sample of compact discs at high resolution: the model predicts that sharp-transition discs have nearly identical sizes at both wavelengths, whereas drift-dominated discs should appear smaller at longer wavelengths. A second, equally direct check is to measure the fragmentation velocity of realistic icy aggregates in the laboratory; if v_frag ≥ 5 m/s, the models themselves show that mm dust forms across the MRI-active region and the dead-zone-size correspondence breaks down.","tokens_in":27522,"feed_emoji":"🪐","tokens_out":3867,"duration_ms":37483,"temperature":0.7,"pith_summary":"The paper proposes that radially compact protoplanetary discs, common in nearby star-forming regions, can form when fragile dust grains fragment in moderately turbulent regions beyond the dead zone. By running one-dimensional dust transport and collision models and converting them into synthetic ALMA observations, it shows that the mm-sized dust disc is truncated at the dead zone edge, so the observed disc size directly traces the dead zone extent. This offers an alternative to the usual explanations (drift-dominated or born-small discs) and predicts specific observables: similar disc sizes at 1.3 and 3 mm, a sharp outer edge that can appear as a low-contrast ring, and a size–luminosity relation consistent with current samples. The mechanism works with dust porosity and with fragmentation velocities up to a few m/s, but requires both fragile dust and a sharp transition in turbulence at the dead zone edge.","feed_headline":"Dead zones decide how big compact discs look","feed_subtitle":"Fragile dust shatters beyond the dead zone, so the mm-dust disc size marks the zone radius - a new way to read ALMA images.","key_machinery":"The load-bearing mechanism is the fragmentation-limited dust growth regime, defined by the balance between turbulent relative velocities and the fragmentation velocity through the Stokes number (St_frag = (1/3 δ_turb) (v_frag/c_s)^2). Beyond the dead zone, α_SS jumps from about $10^{-4}$ to $10^{-3}$ or higher, so the maximum grain size collapses to sub-mm sizes, preventing the growth of detectable mm dust. The dead zone edge acts as a radial gate: a sharp (but continuous) transition in α_SS produces a local steepening of the pressure gradient that accelerates radial drift of ≲ mm dust ('hitchhiking'), which further inhibits growth and creates a slight dust pile-up at the outer edge, visible as a low-contrast ring in synthetic images. The models also track the full dust size distribution with DustPy and post-process it with RADMC-3D and CASA to generate realistic synthetic observations, which are then compared to actual ALMA data. This machinery connects an analytically known growth limit (from Birnstiel et al. 2009, 2012) to a specific disc-structure feature (the dead zone) and uses it to predict an observable property (disc size).","core_discovery":"The central discovery is that the fragmentation of fragile dust (fragmentation velocity ≈ 1 m/s) in moderate turbulence (alpha_SS ≈ $10^{-3}$), as expected in MRI-active regions beyond the dead zone, naturally produces radially compact millimetre-dust discs whose observed size is set by the radial extent of the dead zone. In the models, mm-sized dust cannot grow in the MRI-active region because turbulent relative velocities exceed the fragmentation threshold, so the mm-dust disc is bounded by the transition radius R_t; model discs with R_t = 20, 30, and 50 au have observed radii R_d,90% ≈ 19, 27, and 34 au at 1 Myr. The sharp drop in gas surface density at the dead zone edge accelerates inward drift of sub-mm grains ('hitchhiking'), creating a subtle outer-edge bump that appears as a low-contrast ring or shoulder at high resolution. The paper also shows that pressure bumps inside the dead zone are inefficient dust traps at radii ≲ 10 au because fragile dust fragments near the bump and the small fragments replenish the optically thick inner disc, hiding the trap; only more resilient dust (v_frag ≈ 3 m/s or above) produces a visible ring under those conditions. These synthetic morphologies resemble recent high-resolution observations of compact discs such as Sz 66 and MP Mus.","pith_inferences":["A natural extension the paper leaves implicit is that if dead-zone radius sets disc size, then in a population of compact discs the measured sizes should correlate with stellar properties that control ionization (e.g., X-ray luminosity, cosmic-ray rate, magnetic field geometry) - a test that could be done with existing surveys.","The hitchhiking effect is not unique to dead zones: it should operate at any sharp outward drop in gas surface density, such as at the outer edge of a deep planetary gap or a photoevaporation front; this could be tested in substructured discs showing a weak outer ring at similar contrast.","The degeneracy between turbulence α_SS and fragmentation velocity v_frag highlighted in Section 5.2 implies that dust-size measurements alone cannot separately constrain these parameters; breaking the degeneracy requires independent turbulence estimates from molecular-line widths or dust scale-heights, combined with laboratory and microphysical constraints on v_frag.","The models indicate that dust grains in the outer disc are systematically limited to sub-mm sizes in the MRI-active region even at α_SS = 5 × 10^-4; if real discs have such low turbulence, then planetesimal formation by coagulation in the outer disc appears difficult, which may push the field toward alternative growth mechanisms such as streaming instability or gravitational collapse."],"forward_implications":["If correct, the observed mm-dust size of compact discs becomes a direct tracer of the dead zone radius, turning ALMA size measurements into a probe of dead zone location and hence of disc ionization and magnetic-field structure.","The model predicts nearly identical disc sizes at 1.3 mm and 3 mm for sharp dead zone edges, which can be tested against multi-wavelength ALMA observations; a strong size difference would favour drift-dominated or alternative mechanisms.","High-resolution (≈0.02 arcsec) imaging of seemingly smooth compact discs should reveal a low-contrast ring or a change in slope at the outer edge, as seen in Sz 66 and MP Mus, providing a concrete observational signature of the hitchhiking effect.","Compact discs formed by this mechanism follow the observed size–luminosity relation (R_eff ∝ L_mm^0.5) and naturally yield a shallower slope at later stages, matching the trend seen in older star-forming regions such as Upper Sco.","The model predicts that inner pressure bumps (≲ 10 au) are hidden by the optically thick inner disc, so a faint ring in scattered light or CO may still be present even when the mm continuum looks smooth, informing the interpretation of substructure surveys."],"supporting_citations":[{"why":"Provides the analytic fragmentation-limited Stokes number and particle size (Eq. 11-12) that underpin the prediction that mm dust cannot form in the MRI-active region.","marker":"Birnstiel et al. 2009"},{"why":"Derives the drift-fragmentation limited sizes used to interpret the hitchhiking regime near the dead zone edge.","marker":"Birnstiel et al. 2012"},{"why":"Supplies the dust transport and collision code used to run the one-dimensional models.","marker":"Stammler & Birnstiel 2022 (DustPy)"},{"why":"Provides laboratory evidence that ice grains fragment at velocities as low as ≈1 m/s, the canonical fragility assumed in all main models.","marker":"Musiolik & Wurm 2019"},{"why":"Establishes the dead-zone concept (MRI-quenched region) whose radial extent the models vary as R_t.","marker":"Gammie 1996"},{"why":"Non-ideal MHD simulations used to motivate the slow-transition profile (Model 4) that tests the sharpness assumption.","marker":"Bai et al. 2016"},{"why":"One of the high-resolution observations of compact discs to which the synthetic radial intensity profiles are compared (e.g., Sz 66).","marker":"Miley et al. 2024"},{"why":"Another compact-disc observation (MP Mus) whose outer-edge slope change and weak ring are reproduced by the hitchhiking feature.","marker":"Ribas et al. 2023"},{"why":"Provides the size–luminosity relation predicted for drift-dominated discs, which the paper compares against its own model trend.","marker":"Rosotti et al. 2019"}],"fun_headline_variants":["Dead zone extent sets compact disc size","Dust shattering locks disc radius to dead zone","Compact discs: size follows dead zone edge","Fragile dust beyond dead zone caps disc size"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mechanism requires both fragile dust (fragmentation velocity around 1 m/s) and a sharp jump in turbulence at the dead zone edge; if real dead zone edges are smooth or if grains can survive collisions at speeds above a few m/s, the dead zone size would no longer set the observed disc size.","fun_headline_variants_meta":{"raw":{"variants":["Dead zone extent sets compact disc size","Dust shattering locks disc radius to dead zone","Compact discs: size follows dead zone edge","Fragile dust beyond dead zone caps disc size"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1523,"prompt_tokens":1025,"completion_tokens":498,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":641,"completion_tokens_details":{"reasoning_tokens":440}},"tokens_in":641,"tokens_out":498,"duration_ms":5685,"temperature":1.0,"reasoning_tokens":440,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T22:39:20.620490+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to measure the 1.3-mm and 3-mm disc sizes of a sample of compact discs at high resolution: the model predicts that sharp-transition discs have nearly identical sizes at both wavelengths, whereas drift-dominated discs should appear smaller at longer wavelengths. A second, equally direct check is to measure the fragmentation velocity of realistic icy aggregates in the laboratory; if v_frag ≥ 5 m/s, the models themselves show that mm dust forms across the MRI-active region and the dead-zone-size correspondence breaks down.","supporting_citations":[],"review_version":1}