{"id":"3806d958-fa6a-42d9-8d7a-c8b348cb9bdd","arxiv_id":"2507.08246","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"New JVLA observations of the Class I protostar WL 17 indicate millimeter-sized dust grains in its ring, with a best-fit maximum grain size of about 4.2 mm.","lead":"Astronomers used the JVLA radio telescope to observe the young star WL 17 at seven frequencies and measured how its dust ring emits at centimeter wavelengths. The data suggest that some dust in the ring has grown to millimeter sizes, which may mean planet formation can start very early.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 4.2 mm amax is not uniquely established: the unresolved 18-48 GHz SED can be reproduced by the authors' own free-free components, so the quoted value is an upper limit, not a measurement of millimeter grains.","rationale":"The reader's conditional verdict is appropriate. My concern sharpens it: free-free-II is not a negligible contaminant, because the Table 5 parameters give order-unity optical depth at 20 GHz and non-negligible optical depth at 40 GHz, affecting precisely the frequency range that drives the amax fit. The paper is nevertheless a solid observational study with new multi-epoch JVLA data, transparent uvmodelfit flux measurement, and a two-component model whose limitations are disclosed in Section 3.2.2. The central claim is best read as an upper limit on amax; as a positive detection of 4.2 mm grains in the ring it overreaches, but this does not change the conditional verdict already assigned.","tokens_in":18653,"tokens_out":7412,"duration_ms":93821,"concrete_test":"Re-run the Section 3.2 MCMC with amax of the Grown component fixed to 0.1, 0.3, and 1.0 mm, leaving free-free-I/II EM and Omega free, and compare BIC/delta-chi2 against the amax=4.2 mm model using the same data. If a small-amax plus free-free model fits within delta-chi2 < ~2 per added freedom, the millimeter-grain claim is not required. A spatial test would be decisive: JVLA A-configuration at 30-50 GHz (<=50 mas) to see whether the 18-48 GHz emission is extended on dust-ring scales or dominated by a compact ionized knot.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2.2 constrains the Grown dust component with the 18-48 GHz JVLA bands and ALMA Band 3 at 100 GHz. The source is unresolved at >0.5\" (Section 2), so the decomposition into dust plus free-free-I/II is purely spectral. The paper's own best-fit free-free-II (Table 5) has EM ~ 5.8e8 to 1.1e9 cm^-6 pc and Omega ~ 0.06-0.14e-14 sr; with T=8000 K and Eq. A1 this component has tau ~ 1 at 20 GHz and tau ~ 0.3 at 40 GHz, so it contributes substantially across exactly the bands used to measure amax. A mixture of partially optically thick free-free (alpha < 2) with dust can flatten the observed 14-48 GHz spectral indices, and a model with much smaller amax plus stronger free-free may fit equally well. The authors acknowledge this in Section 3.2.2, stating that amax and Sigma_dust should be regarded as upper limits if spinning dust or optically thick free-free contributes. The quoted uncertainty (+1.8/-2.1 mm, Table 4) comes from the MCMC sampler scatter only and does not include this systematic. Therefore the abstract's claim that 'structures internal to the ring achieved a maximum grain size of approximately 4.2 mm' goes beyond what the unresolved SED fit can establish.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents JVLA observations of the Class I protostar WL 17 at 2–48 GHz with five epochs, detecting unresolved emission at 4–48 GHz. The authors construct a spectral energy distribution combining their JVLA fluxes with ALMA data from the literature and fit it with a radiative transfer model consisting of two dust components (a 'Grown' high-column-density component and a 'Small' low-column-density halo) plus one or two free-free emission components. The best fit yields a maximum grain size of amax = 4.2 +1.8/-2.1 mm for the Grown dust component (Table 4). On this basis the authors argue that millimeter-sized grains are present in the WL 17 ring and discuss the ring's gravitational stability and the possibility of forming a planetary core by pebble accretion.","tokens_in":19055,"tokens_out":4064,"duration_ms":47762,"significance":"If the inferred grain size is robust, this would be one of the few measurements of millimeter-sized grains in a Class I disk and would be relevant to early planet formation. The paper reports a careful, multi-epoch JVLA campaign, uses visibility-domain flux fitting, and specifies the calibration and error treatment in detail. The radiative transfer and MCMC fitting are described transparently, and the authors correctly identify several degeneracies (T_dust, Omega_dust, free-free contamination). The main weakness is that the headline claim is presented as a measurement in the abstract while the body text repeatedly characterizes it as an upper limit; the unresolved nature of the source and the spectral decomposition into dust and free-free make the '4.2 mm' value model-dependent rather than uniquely constrained. The data set is useful and the caveats are partly acknowledged, but the framing needs substantial revision before the paper can be accepted.","major_comments":[{"comment":"The pebble-accretion core mass calculation uses the central values from Table 4 (Mdust ~ 1007 M_earth and amax = 4.2 mm) without propagating the stated upper-limit caveat. While the text notes the result 'should also be considered an upper limit,' the range of Mcore corresponding to the allowed ranges of Mdust and amax (or to the alternative free-free-dominated decomposition) is not given. Please provide a numerical range for Mcore under the systematic variations discussed in §3.2.2, since the qualitative statement is insufficient for the quantitative claim that a gas giant could form if the gas-to-dust ratio is below ~10.","section":"§4.2, Eq. (2)"}],"minor_comments":[{"comment":"The program ID is given as 23A-124 in the text but 24A-001 in Table 1; please correct the inconsistency.","section":"Section 2 vs Table 1"},{"comment":"The phrase 'spatial resolution exceeding 0.5 arcsec' is ambiguous; it should read 'spatial resolution coarser than 0.5 arcsec' or 'with a beam size of ≳0.5 arcsec.'","section":"Abstract"},{"comment":"Equation (2) uses 'fr s' and 'frp' while the text defines 'frs' and 'frp'; please unify the notation. Also, 'Mdust = 1007 M_earth' has an unwarranted precision; prefer 3.1 M_Jup or 1000 M_earth.","section":"§4.2, near Eq. (2)"},{"comment":"The column header 'Adapted error' appears to mean 'Adopted error'; please correct the typo.","section":"Table 2"},{"comment":"The statement 'Some symbols are larger than their error bars' is informal and underspecified; please clarify which epochs/frequencies are meant or adjust the plotting so the error bars are visible.","section":"Figure 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a valuable multi-epoch JVLA data set for an interesting Class I disk, but the abstract and conclusion currently overstate the robustness of the amax measurement. The authors' own text in §3.2.2 correctly labels the result as an upper limit, so the revision path is clear: reframe the claims and add the requested sensitivity analysis. I would not recommend rejection, as the observational data and modeling effort are solid and the caveated result is still of interest. However, the paper needs a substantive reframing before publication, not just cosmetic edits."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful data paper. The team observed the Class I protostar WL 17 with the JVLA across 2–48 GHz in five epochs, got clean detections at 4–48 GHz, and fit the centimeter SED with a two-dust-component plus one-or-two-free-free-component model. The data handling looks careful; uvmodelfit and peak-flux measurements agree within 10%, and the error budget includes both thermal noise and absolute calibration. The genuinely new result is the centimeter-wavelength SED for this source and the constraint, conditional on DSHARP opacities, that the high-column-density dust has amax near 4 mm. That adds a data point to the small set of cm grain-size constraints in Class I disks.\n\nThe soft spot is the gap between the abstract and what the fit can support. The ring is unresolved in these data (beam >0.5\"), so the split into 'grown' and 'small' dust and free-free components is purely spectral. The authors acknowledge in §3.2.2 that amax and Sigma_dust should be read as upper limits if spinning dust or optically thick free-free contributes. That caveat belongs in the abstract. The stress-test math checks out: their own best-fit free-free-II has tau ~1 at 20 GHz and ~0.3 at 40 GHz, so it contributes significantly in exactly the bands that constrain amax, and a smaller amax with stronger free-free could plausibly fit the same points. The quoted +1.8/-2.1 mm is MCMC scatter, not a systematic error budget over the free-free/opacity degeneracy. So the defensible claim is an upper limit of a few mm, not a point measurement of 4.2 mm.\n\nMinor issues: T_dust fixed at 35 K and small-grain amax fixed at 45 µm deserve a sensitivity test; the pebble accretion discussion is speculative but clearly framed as such; the gravitational stability argument rests on a gas-to-dust ratio they admit is unknown. None of this sinks the paper.\n\nMy take: it deserves a serious referee. The data are real, the analysis is transparent, and the caveats are mostly present. The revision should state the upper-limit interpretation in the abstract and lead with that. I'd cite it as an upper limit and bring it to a reading group.","headline":"Solid new JVLA data on WL 17, but the 4.2 mm grain size is an upper limit, and the abstract should say so.","tokens_in":19578,"tokens_out":2888,"would_cite":true,"duration_ms":29687,"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":"The paper claims that WL 17's dust ring contains grains grown to about 4.2 mm, based on JVLA 2–48 GHz observations combined with ALMA data and radiative transfer modeling.","keywords":["Circumstellar disks","Dust continuum emission","Planet formation","Protoplanetary disks","Protostars","Spectral energy distribution","Very Large Array","Maximum grain size"],"falsifier":"Take the JVLA in A-configuration at 30–50 GHz to image the WL 17 ring at ≲50 mas resolution. If the 18–48 GHz emission turns out to be spatially smooth and extended rather than concentrated in a compact, high-column-density substructure, the two-component SED decomposition that yields $a_{\\rm max} = 4.2$ mm would be falsified.","tokens_in":18466,"feed_emoji":"🪐","tokens_out":13537,"duration_ms":131137,"temperature":0.7,"pith_summary":"The paper sets out to measure how large dust grains have grown in the compact ring around WL 17, a Class I protostar only a few hundred thousand years old, using JVLA observations from 2 to 48 GHz. By fitting the observed spectral energy distribution with radiative transfer models that assume the DSHARP dust opacity, the authors conclude that the ring's high-column-density component contains grains with a maximum size of about 4.2 mm. Because grain size controls how efficiently streaming instability and pebble accretion build planets, the result would mean that millimeter-sized particles — and possibly the first steps of planet formation — already exist in a very young disk. The authors also use the inferred dust mass and grain size to argue that, if the gas-to-dust ratio is around 10, pebble accretion in the ring could form a core massive enough to trigger runaway gas accretion and produce a gas giant.","feed_headline":"Dust grains up to 4.2 mm found in young star WL 17's ring","feed_subtitle":"If correct, it means planet-building dust can appear before the star finishes forming.","key_machinery":"The central mechanism is the interpretation of the centimeter-to-millimeter spectral energy distribution as the emission of geometrically flat, isothermal, optically thick dust slabs, using the analytic radiative-transfer solution of Birnstiel et al. (2018) (following Miyake & Nakagawa 1993) with the DSHARP dust opacity model and a power-law grain-size distribution $n(a)\\propto a^{-3.5}$. The observed SED is decomposed into a low-column-density 'Small' component ($a_{\\rm max} = 45\\,\\mu$m) and a high-column-density 'Grown' component whose $a_{\\rm max}$ is a free parameter, plus one or two free-free emission components whose emission measure and solid angle absorb the time-varying low-frequency signal. The curvature of the 18–48 GHz spectrum — where optically thick dust gives a spectral index near 2 that rolls off as the opacity drops — is what constrains $a_{\\rm max}$, and the same model, fit with an MCMC routine, simultaneously yields the column density, solid angle, and dust mass of each component.","core_discovery":"The paper's central claim is that the spectral energy distribution of WL 17's ring from 2 to 345 GHz is best reproduced by a model with two dust components: a low-column-density 'Small' component with $a_{\\rm max} = 45\\,\\mu$m that dominates the ALMA (sub)millimeter emission, and a high-column-density 'Grown' component that is optically thick up to ~30 GHz and yields a maximum grain size of $a_{\\rm max} = 4.2^{+1.8}_{-2.1}$ mm under the DSHARP dust opacity model (a standard prescription for the emission and absorption opacities of compact dust grains). The 18–48 GHz flux is attributed to this optically thick dust slab, whose spectral index of about 2 flattens at lower frequencies; the curvature of the SED across the JVLA bands is what pins down the grain size. Time-variable free-free emission accounts for the <14 GHz behavior and for part of the 14–18 GHz excess, but the authors argue that without a dust component of roughly millimeter-sized grains it would be difficult to keep the emission optically thick at Q band (40–48 GHz). They therefore conclude that grain growth to millimeter sizes has already occurred in the WL 17 ring, while noting that both free-free and spinning dust could raise the non-dust fraction of the 18–48 GHz emission, making the inferred $a_{\\rm max}$ an upper limit.","pith_inferences":["If the 4.2 mm grain size survives higher-resolution imaging, WL 17 would join HL Tau and a few other Class 0/I disks where substructures and large grains appear within the first ~0.5 Myr, strengthening the case that planet formation is not confined to the Class II stage.","The two-component structure of the model is not spatially verified; a natural test is to image the ring with the JVLA in A-configuration at 30–50 GHz and check whether the high-column-density 'Grown' component coincides with a compact substructure, or whether the unresolved SED fit has artificially split a single component.","A broader extension would be to apply the same two-component SED fit to other Class I disks with cm-wavelength data, using the time variability of the free-free components as a tag to separate dust from ionized gas; if mm-sized grains are common, early pebble accretion could be the default pathway to gas giants."],"forward_implications":["Grain growth to ~4 mm can occur within the Class I stage (age ≲1 Myr), before the disk reaches the Class II phase, so millimeter-sized pebbles are available for planetesimal formation earlier than the standard core-accretion timeline assumes.","If the gas-to-dust ratio in the ring is ~10, pebble accretion around an already-formed planetesimal could build a core of ~16 $M_\\oplus$, exceeding the critical core mass of ~6 $M_\\oplus$ needed for runaway gas accretion and potentially producing a gas giant.","The ring may be gravitationally unstable if the disk-to-star mass ratio is above ~0.1, a condition the authors flag as an upper limit because the gas mass is poorly constrained.","The existence of a compact, optically thick component with mm-sized grains would make the ring's substructures directly testable: future observations at ≲50 mas resolution in the 30–50 GHz bands could image the narrow rings or vortices that harbor the grown dust.","Because free-free and spinning dust could contribute to the 18–48 GHz emission, the inferred $a_{\\rm max}$ of 4.2 mm is an upper limit; if either mechanism is significant, the true maximum grain size could be much smaller."],"supporting_citations":[{"why":"Supplies the DSHARP dust opacity model and the analytic optically thick dust-slab radiative transfer solution used to compute model SEDs; the central assumption for converting observed SED curvature into grain size.","marker":"Birnstiel et al. (2018)"},{"why":"The SED modeling approach (two dust components plus free-free) that the paper follows directly; provides the formulation for combined dust and free-free emission.","marker":"Liu et al. (2019b)"},{"why":"Same modeling approach; provides the fitting framework and the treatment of optically thick dust components.","marker":"Liu et al. (2021)"},{"why":"Established the WL 17 dust ring, its 0.1 arcsec scale, low inclination, and ring properties; supplies the ALMA Band 3 (100 GHz) flux used in the SED fit and the ring location r0 = 16 au.","marker":"Sheehan & Eisner (2017)"},{"why":"Provides the ALMA fluxes at 100–345 GHz and the spectral index α = 2.30 ± 0.01 at >90 GHz that anchor the optically thin 'Small' dust component in the fit.","marker":"Gulick et al. (2021)"},{"why":"Gives the stellar mass of 0.3 M⊙ and the spatially resolved spectral-index analysis that supports the two-component interpretation with different dust distributions.","marker":"Han et al. (2023)"},{"why":"Recent imaging of the WL 17 ring providing inclination, temperature, and azimuthal asymmetry constraints used in the model setup and discussion.","marker":"Shoshi et al. (2024)"},{"why":"Supplies the free-free optical depth formula used to separate the time-variable ionized-gas component from the dust emission, which is essential for isolating the grain-size signal.","marker":"Keto (2003)"},{"why":"Provides the semi-analytic pebble-accretion formula used to estimate the planetary core mass that could form in the ring.","marker":"Jiang & Ormel (2023)"}],"fun_headline_variants":["Millimeter grains found in young star's ring","Early planet seeds: 4.2 mm grains in WL 17 disk","Class I star's ring hosts grains up to 4.2 mm","Dust growth to millimeter sizes in star-forming disk","4.2 mm grains in young star's ring hint at early planets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The grain-size result depends on the assumption that the 18–48 GHz radio emission is dominated by thermally radiating dust that is optically thick at those frequencies, with free-free emission and spinning dust making only minor contributions; if either of those contributes more than modeled, the inferred 4.2 mm maximum grain size would be an upper limit and could be much smaller.","fun_headline_variants_meta":{"raw":{"variants":["Millimeter grains found in young star's ring","Early planet seeds: 4.2 mm grains in WL 17 disk","Class I star's ring hosts grains up to 4.2 mm","Dust growth to millimeter sizes in star-forming disk","4.2 mm grains in young star's ring hint at early planets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001524,"raw_usage":{"total_tokens":6193,"prompt_tokens":1126,"completion_tokens":5067,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":742,"completion_tokens_details":{"reasoning_tokens":4979}},"tokens_in":742,"tokens_out":5067,"duration_ms":40343,"temperature":1.0,"reasoning_tokens":4979,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:23:06.359816+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the JVLA in A-configuration at 30–50 GHz to image the WL 17 ring at ≲50 mas resolution. If the 18–48 GHz emission turns out to be spatially smooth and extended rather than concentrated in a compact, high-column-density substructure, the two-component SED decomposition that yields $a_{\\rm max} = 4.2$ mm would be falsified.","supporting_citations":[{"cited_title":"B., Tsai, A.-L., Chen, W","cited_arxiv_id":null,"evidence_quote":"Same modeling approach; provides the fitting framework and the treatment of optically thick dust components."},{"cited_title":"D., & Eisner, J","cited_arxiv_id":null,"evidence_quote":"Established the WL 17 dust ring, its 0.1 arcsec scale, low inclination, and ring properties; supplies the ALMA Band 3 (100 GHz) flux used in the SED fit and the ring location r0 = 16 au."},{"cited_title":"C., Sadavoy, S., Matrà, L., Sheehan, P., & van der Marel, N","cited_arxiv_id":null,"evidence_quote":"Provides the ALMA fluxes at 100–345 GHz and the spectral index α = 2.30 ± 0.01 at >90 GHz that anchor the optically thin 'Small' dust component in the fit."}],"review_version":1}