{"id":"cd10d2a3-3ce8-464d-afcc-61cb82ee27b9","arxiv_id":"2507.14443","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Modeling six wavelengths of HL Tau's disk shows small (about 100 micron) grains beyond 40 au and suggests planet formation by pebble accretion is possible only if the outer dust is porous and rich in organic material.","lead":"We mapped the dusty disk around the young star HL Tau at six radio wavelengths and fitted a computer model to find what the dust grains are made of, how porous they are, and how big they grow. The results suggest the outer disk holds fairly small grains and, depending on the dust composition, could either feed planet formation quickly or not at all.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Polarization prior in Eq. 15, with its global C≈0.02 threshold, is the load-bearing link to the porous organics-rich branch; if C or the radius-independent assumption is off, the high pebble-accretion conclusion loses its anchor.","rationale":"The paper is a careful, honest MCMC analysis of an excellent multi-wavelength dataset. The new Band 3 data and the corrected Band 4 data genuinely improve on previous work, and the derived outer spectral index ~3.7 and the corresponding amax_ffill~100 µm for compact/AC-rich dust are likely robust; I would not contest the observational characterization. The load-bearing weakness is the step from Stokes I to the favored microphysical model. Section 4.2 shows that four combinations of (fAC,ffill) all fit the intensity profiles, so the Stokes I data alone cannot distinguish them. The abstract's exclusion of compact amorphous-carbon-rich dust and preference for ffill~0.03–0.3 comes from the polarization prior in Eq. (15), which the authors explicitly flag as an approximation (P≈C·P90ω_eff with C≈0.02, threshold profile constant across the disk). If C or the effective threshold changes by even a factor of two—well within the uncertainty of the scattering-to-polarization conversion—the posterior over ffill and fAC shifts, and the pebble accretion rate '10 M⊕ Myr−1 within 1 Myr' in §5.4 is no longer the favorable branch. The paper is conditional in wording, so I would not reject it, but the CONDITIONAL verdict is correct. My proposed check—rerunning the MCMC with a radius-dependent polarization likelihood and sampling C over its plausible range—would settle whether the planet-formation claim survives. Independent literature (Zhang et al. 2023; Lin et al. 2024) supports moderate porosity, but this paper's internal argument must stand on its own prior.","tokens_in":25740,"tokens_out":9998,"duration_ms":119000,"concrete_test":"Re-run the MCMC at r=40–120 au replacing the global p2 in Eq. (15) with a radius-dependent likelihood based on the observed 0.87 mm polarization fraction profile (e.g., from Stephens et al. 2023), sampling C over 0.01–0.04 and allowing the ring/gap decomposition uncertainty; check whether the 68% credible interval for log ffill still lies within roughly [−1.5, −0.5] and whether p(ffill<0.03)+p(ffill>0.3) remains small, and whether the fAC≈0.8 branch is excluded. If the posterior shifts so that compact/AC-rich or very porous dust is allowed, the strongest pebble-accretion claim in §5.4 is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the polarization prior p2 in Eq. (15), applied uniformly across all radii. The paper itself states in §4.3 that the actual P90ω_eff 'is expected to vary with radius' and that Eq. (15) is adopted only because 'a definitive profile of the scattering polarization degree has not yet been obtained.' This prior is what converts the observed 0.87 mm polarization (scattering component ≳0.5% in rings) into a model probability via P ≈ C·P90ω_eff with a single constant C≈0.02 and a global threshold (p2≈1 for P90ω_eff≳0.25, ≈0 for ≲0.05). Section 4.2 shows that all four fixed dust models (compact/organics-rich, porous/organics-rich, compact/AC-rich, porous/AC-rich) reproduce the Stokes I profiles, so the abstract's claim that 'compact, amorphous-carbon-rich dust is unlikely, and moderately porous dust is favored'—and hence the pebble-accretion rate of ~10 M⊕/Myr—rests almost entirely on this approximate prior. If C is actually 0.01 or 0.04, or if the scattering polarization fraction in the outer rings is lower or higher than 0.5%, the threshold moves and the 68% credible interval for ffill (0.03–0.3) and the exclusion of fAC≈0.8 could shift, changing the planet formation verdict. The concern is not that the prior is absurd; it is that the quantitative conclusion is not robust to its stated approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a resolved multi-wavelength study of the HL Tau disk, modeling azimuthally averaged radial intensity profiles at 0.45, 0.87, 1.3, 2.1, 3.1, and 7.9 mm with a common 0.05 arcsec resolution. The authors fit, independently at each radius, six dust parameters (temperature, surface density, maximum grain size, amorphous-carbon fraction, filling factor, and size-distribution power-law index) using a standard scattering-aware intensity model and an MCMC implementation. They introduce a corrected 2.1 mm image that includes short-baseline ALMA data and add new 3.1 mm data. The principal claims are that the outer disk (r ≳ 40 au) has a long-wavelength spectral index of about 3.7; that amax ffill is about 100 μm for compact or amorphous-carbon-rich dust; that the 0.87 mm scattering polarization, encoded in a prior p2, favors moderately porous (ffill ~ 0.03-0.3), organics-rich (fAC ≲ 0.4) dust; and that, under that branch, the inferred surface density and grain size imply pebble accretion rates near 10 M⊕/Myr at r > 40 au, enabling giant-planet cores to form within 1 Myr. The paper also predicts that 1.3 and 2.0 cm ngVLA observations can distinguish the organic-rich and amorphous-carbon-rich models.","tokens_in":26079,"tokens_out":6705,"duration_ms":67251,"significance":"If the dust-property inference is robust, this is a valuable contribution: it combines the broadest wavelength coverage (0.45-7.9 mm) yet used for a resolved HL Tau analysis, corrects a known short-baseline problem in the 2.1 mm data, and provides quantitative posterior distributions rather than single best-fit models. The new 3.1 mm data and the Band 4 correction are concrete observational improvements, and the ngVLA predictions in Section 5.5 are a useful falsifiable output. The main significance is conditional, however: Section 4.2 shows that the Stokes I data alone do not distinguish compact, porous, organic-rich, and amorphous-carbon-rich models, so the paper's headline planet-formation conclusion rests largely on the approximate polarization prior in Eq. (15). The pebble accretion rates in Section 5.4 are also a remapping of the same MCMC posteriors rather than an independent test. With sensitivity analyses for the prior and a more careful framing of the planet-formation implications, the paper would be a solid advance; in its current form, the central claims are plausible but not fully supported.","major_comments":[{"comment":"The polarization prior p2 is load-bearing but is applied as a single global function of P90 ω_eff at 0.87 mm, with C ≈ 0.02 and one threshold profile for the entire disk, even though the paper itself states that 'the actual value of P90ωeff in the disk is expected to vary with radius' and that Eq. (15) is adopted because a definitive profile has not been obtained. Section 4.2 shows that all four fixed composition/porosity models reproduce the observed Stokes I profiles at similar quality, so the exclusion of compact amorphous-carbon-rich dust and the preference for 0.03 ≲ ffill ≲ 0.3 are effectively determined by this approximate prior. Because the Section 5.4 pebble accretion conclusion is then evaluated on the porous organics-rich branch selected by this prior, the headline claim is not robust as presented. Please add sensitivity tests that vary C over a plausible range (e.g., 0.01 and 0.04), vary the location and width of the threshold in Eq. (15), and adopt a radially varying prior informed by the observed ring/gap polarization fractions; report how the 68% credible intervals for ffill, fAC, and the pebble accretion rate shift under these changes.","section":"Section 4.3, Eq. (15)"},{"comment":"At each radius, six free parameters are constrained by six observed intensities through Eq. (6), so the posterior is inevitably shaped by the priors in Eqs. (9) and (15). The paper acknowledges this in part, and Section 4.2 shows that fixing fAC and ffill leaves the intensity fits essentially unchanged, but the reported 68% intervals for ffill, fAC, and pd (Figures 3 and 6) are presented as empirical constraints even though Figures 3, 5, and 11 show broad, often multimodal posteriors. Please quantify the effective information added by the data, for example by comparing the full six-parameter model with reduced models using the Bayesian information criterion or nested-sampling evidence, and report prior-sensitivity runs (e.g., uniform versus log-uniform priors on T, Σd, and amax) so the reader can distinguish data-driven constraints from prior-driven constraints.","section":"Section 3.3 and Section 4.1"},{"comment":"The pebble accretion rate and the statement that 'a giant planet core can form within 1 Myr' are computed by feeding the MCMC-fitted Σd and amax into Eqs. (16)-(24); they are therefore not an independent test of planet formation but a remapping of the same posterior that already selected the porous organics-rich branch. The paper is partly careful in using conditional language and noting the 3 M⊕ seed assumption, but the abstract and Section 6 state the 10 M⊕/Myr rate and the 1 Myr core formation as general conclusions. Please restate these as model-dependent consequences of the polarization-prior-selected branch, propagate the full posterior including model uncertainty among the four branches rather than only the favored branch, and test sensitivity to the assumed gas surface density from 13C17O, since the Stokes numbers and drift velocities in Eqs. (14), (17), and (18) scale with Σg.","section":"Section 5.4, Eqs. (16)-(24)"}],"minor_comments":[{"comment":"The rendered title contains typographical artifacts: 'Multi-W avelength' and 'T au' should be corrected to 'Multi-Wavelength' and 'Tau'.","section":"Title"},{"comment":"The 3.1 mm beam is described as 0.061 × 0.041 arcsec and then said to have a beam area equivalent to 0.05 × 0.05 arcsec; please state the effective circular beam size explicitly in the caption to avoid confusion.","section":"Section 2, Figure 1"},{"comment":"The sentence introducing Eq. (15) derives the threshold from P ≈ C P90ω_eff with C ~ 0.02 in a single clause; please give the uncertainty in C and in the observed scattering polarization fraction (≳0.5%) in the main text, since these directly set the prior in Eq. (15).","section":"Section 4.3"},{"comment":"The ngVLA integration-time estimates assume ffill = 0.1 and only compare fAC = 0.3 and 0.8; the quoted times are therefore conditional on the porous branch and do not include the full posterior spread in ffill and fAC. Please clarify this in the text.","section":"Section 5.5"},{"comment":"The statement that the inferred surface density is 'consistent with the gravitational instability threshold' should explicitly note the adopted fd2g = 0.01 and Q = 1.4 values, since Eq. (11) scales with these assumptions.","section":"Section 6, summary item 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is candid about its limitations, and the data-processing improvements are genuine. I am not recommending rejection because the core observational finding—the outer-disk spectral index of about 3.7 and the associated small amax ffill for compact/amorphous-carbon-rich dust—is likely robust and interesting. The problem is the quantitative overreach from the approximate polarization prior to the pebble-accretion conclusion. If the authors add the requested sensitivity analyses and reframe the planet-formation statements as conditional, the paper would be publishable; without those changes, the headline claim is not supported to the standard implied by the abstract."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real contribution here is the data. The authors fixed the missing short-baseline problem in the 2.1 mm ALMA image and added a new 3.1 mm map, and that alone moves the outer-disk spectral index to ~3.7 and the preferred amax*ffill down to ~100 microns. That is a concrete revision of Carrasco-González et al. (2019), and it is well documented. The MCMC treatment of composition, porosity, and size distribution as free parameters is a natural next step, and the paper is honest about degeneracies: it shows compact dust and amorphous-carbon-rich dust both fit Stokes I, and it clearly flags when priors are doing the work.\n\nThe soft spot is exactly where the stress-test puts it. The polarization prior in Eq. 15 is the step that excludes compact AC-rich dust and prefers ffill = 0.03–0.3, but it rests on a single global constant C ~ 0.02 and a radius-independent threshold that the authors concede is a placeholder. I don't think that is fatal—the prior is grounded in existing scattering polarization models—but it does mean the compositional and porosity distinction is not as secure as the abstract's wording suggests. The pebble accretion rates in Section 5.4 are propagated from those fitted parameters, so they are scenarios, not predictions. The paper itself notes the inner-disk inconsistency (t_life < 1 Myr) rather than hiding it, which I take as a good sign.\n\nI would send this to review. The data correction alone is worth publishing, and the analysis gives other groups a clear template plus testable long-wavelength predictions (Figure 8). I would push the authors to either get more constraints on the polarization conversion or soften the claims accordingly, but the core observational result—small grains in the outer disk, spectral index ~3.7—should survive scrutiny. No code or products are shipped, but the recipe is detailed enough to reproduce. For anyone working on HL Tau or disk dust evolution, this is a useful reference.","headline":"The paper delivers a genuine data improvement—corrected 2.1 mm and new 3.1 mm images—that revises the outer-disk grain size downward, but its stronger claims about porosity and pebble accretion rest on a polarization prior the authors themselves flag as approximate.","tokens_in":26717,"tokens_out":1703,"would_cite":true,"duration_ms":21554,"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":"Six-wavelength maps pin HL Tau's outer dust to ~100 micron grains, making pebble accretion fast enough to build giant-planet cores within 1 Myr.","keywords":["HL Tau","protoplanetary disks","dust opacity","dust porosity","pebble accretion","spectral index","millimeter interferometry","scattering polarization"],"falsifier":"Measure the radial profile of 0.87 mm scattering polarization in HL Tau at higher signal-to-noise: if the polarization efficiency in the rings falls below the threshold where the adopted prior drops to zero, or varies strongly with radius, the prior in Eq. (15) is invalidated and the porous-dust conclusion must be rederived. Alternatively, a 1.3 cm intensity measurement beyond 40 au that falls outside the factor-of-two separation predicted between the $f_{\\rm AC}=0.3$ and $f_{\\rm AC}=0.8$ models would falsify one of the two compositional branches.","tokens_in":25512,"feed_emoji":"🪐","tokens_out":7756,"duration_ms":76369,"temperature":0.7,"pith_summary":"Using six wavelengths from 0.45 to 7.9 mm at 0.05 arcsecond resolution, this paper tries to pin down what the dust in the HL Tau disk is made of, how porous it is, how big the grains are, and how much mass sits in the disk. The key new observational inputs are a corrected 2.1 mm image that recovers flux lost when only long baselines were used, and a new 3.1 mm image; together they show that the outer disk ($r\\gtrsim40$ au) has a steep spectral slope $\\alpha\\sim3.7$ at long wavelengths. The authors argue that this slope, combined with the scattering polarization already observed at 0.87 mm, rules out compact grains made mostly of amorphous carbon and favors moderately porous, organics-rich grains with an effective maximum size near $100\\,\\mu$m beyond 40 au. If that dust model is correct, the inferred surface density and grain size imply pebble accretion rates near $10\\,M_\\oplus\\,{\\rm Myr}^{-1}$, enough to grow a giant-planet core within about a million years. If the dust is instead amorphous-carbon-rich, pebble accretion is too slow, and planet formation in this disk would need another route such as gravitational instability.","feed_headline":"HL Tau's outer dust can build giant planet cores in 1 Myr","feed_subtitle":"Six-wavelength maps plus polarization favor porous organic grains, making pebble accretion fast enough to form a giant core within 1 Myr.","key_machinery":"The engine of the analysis is a one-dimensional, scattering-inclusive radiative transfer model for the radial intensity, with dust temperature, surface density, maximum grain size, amorphous-carbon fraction, filling factor, and size-distribution index treated as free parameters in Markov chain Monte Carlo fits at each radius. Opacities come from a standard dust mixture of water ice, silicates, troilite, and refractory organics, in which part of the carbonaceous material is replaced by amorphous carbon (parameter $f_{\\rm AC}$); porosity is included as vacuum mixed into the grains, and the key grain-scale quantity is $a_{\\rm max}f_{\\rm fill}$, the mass-to-area ratio that controls both opacity and aerodynamic drift. The decisive additional input is a polarization-based prior that converts the modeled polarization efficiency $P_{90}\\omega^{\\rm eff}_{345\\,{\\rm GHz}}$ into a probability, using the observed ring polarization to exclude parameter combinations that scatter too little. The corrected 2.1 mm image built from concatenated short- and long-baseline data is also part of this machinery, because the old long-baseline-only image artificially suppressed outer-disk flux and steepened the spectral index.","core_discovery":"The central claim is that the HL Tau disk's outer region is populated by small, fluffy, organic-rich grains rather than the millimeter-sized compact grains suggested by earlier four-wavelength fits. Modeling six azimuthally averaged intensity profiles (0.45, 0.87, 1.3, 2.1, 3.1, and 7.9 mm) with a scattering-aware radiative transfer formula, the paper finds a spectral index $\\alpha\\sim3.7$ between 0.87 and 3.1 mm beyond 40 au, which drives the product of maximum grain radius and filling factor, $a_{\\rm max} f_{\\rm fill}$, down to roughly $100\\,\\mu$m for compact or amorphous-carbon-rich dust. Adding a prior constructed from the observed 0.87 mm scattering polarization, the compact amorphous-carbon-rich branch becomes unlikely, and the filling factor is bracketed between 0.03 and 0.3. With moderately porous organic-rich dust, the same intensity profiles imply a high enough dust surface density and small enough drift speeds that the radial pebble flux reaches up to $\\sim10^3\\,M_\\oplus\\,{\\rm Myr}^{-1}$, and the authors conclude that pebble accretion could build a giant-planet core within 1 Myr. The amorphous-carbon-rich scenario instead leaves gravitational instability as the more plausible formation mechanism.","pith_inferences":["The polarization prior carries much of the weight: if future imaging resolves a radially varying scattering-polarization profile or finds ring polarization below the threshold used here, the inferred porosity range and the fast-pebble-accretion conclusion would need to be revisited.","The short-baseline correction may matter beyond HL Tau: other disks analyzed with long-baseline-only ALMA data could have artificially depressed outer-disk fluxes and overstated spectral indices, so re-reducing those datasets might change dust-size estimates elsewhere.","A natural testable extension is to apply the same six-parameter fitting with a polarization prior to other disks that have measured ring polarization, checking whether the porous organic-rich solution is generic or specific to HL Tau.","The leaky-gap result implies that dust surface density inside gap regions is not strongly depleted; combining these opacity-based densities with gas tracers could test whether the gaps are carved by planets or by dust evolution effects."],"forward_implications":["The steep spectral slope at long wavelengths in the outer disk rules out mm-sized compact grains as the dominant population; the remaining solutions are small grains with $a_{\\rm max}f_{\\rm fill}\\sim100\\,\\mu$m or moderately porous dust with sizes up to roughly 1 mm.","If the porous organics-rich interpretation holds, the HL Tau gaps are leaky: neither dust surface density nor grain size drops sharply across the gaps, so the gaps do not efficiently trap pebbles and the inner disk is not starved of dust.","Pebble accretion at roughly $10\\,M_\\oplus\\,{\\rm Myr}^{-1}$ allows a few-Earth-mass seed to grow into a giant-planet core in under 1 Myr, making pebble accretion a viable formation path in this disk.","If the dust is instead amorphous-carbon-rich, the low surface density and small grain size keep pebble accretion far slower, favoring disk fragmentation by gravitational instability as the origin of the observed substructures.","Centimeter-wavelength observations near 1.3 and 2.0 cm should separate the organics-rich and amorphous-carbon-rich models by roughly a factor of two in brightness beyond 40 au, so the compositional ambiguity is testable."],"supporting_citations":[{"why":"Supplies the 0.87, 1.3, and 7.9 mm data and the earlier four-wavelength analysis whose dust-size conclusion this paper revises.","marker":"Carrasco-González et al. (2019)"},{"why":"Supplies the 0.45 mm Band 9 data and the previous multi-wavelength constraints that the six-wavelength fit extends.","marker":"Guerra-Alvarado et al. (2024)"},{"why":"Provides the 3.1 mm Science Verification data used as one of the six wavelengths.","marker":"ALMA Partnership et al. (2015)"},{"why":"Reports the 0.87 mm scattering polarization in rings and gaps that anchors the polarization-based prior.","marker":"Stephens et al. (2023)"},{"why":"Interprets Stokes I and polarized emission with a porous-dust model, supporting the porosity range favored here.","marker":"Zhang et al. (2023)"},{"why":"Provides multi-wavelength polarimetric modeling of HL Tau that also favors moderately porous dust.","marker":"Lin et al. (2024)"},{"why":"Supplies the approximation that converts modeled polarization efficiency to an observed polarization degree, used in the prior.","marker":"Kataoka et al. (2016)"},{"why":"Defines the base dust composition with refractory organics that the amorphous-carbon fraction parameter modifies.","marker":"Birnstiel et al. (2018)"},{"why":"Provides the 13C17O-based gas surface density used to compare dust surface densities and assess gravitational instability.","marker":"Booth & Ilee (2020)"},{"why":"Provides the analytical pebble flux and pebble accretion framework used for the planet-formation implications.","marker":"Lambrechts & Johansen (2014)"}],"fun_headline_variants":["Porous grains let HL Tau build giant cores in 1 Myr","HL Tau's fluffy dust speeds planet core formation","Pebble accretion wins in HL Tau with porous grains","Organic-rich dust makes HL Tau planet formation fast","Six-wavelength fits favor porous grains for HL Tau cores"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the observed 0.87 mm scattering polarization maps onto dust properties through one global conversion formula applied across the entire disk, even though the true polarization profile likely varies with radius; if that conversion is wrong, the preference for porous organics-rich dust and the fast-pebble-accretion scenario lose their main independent observational anchor.","fun_headline_variants_meta":{"raw":{"variants":["Porous grains let HL Tau build giant cores in 1 Myr","HL Tau's fluffy dust speeds planet core formation","Pebble accretion wins in HL Tau with porous grains","Organic-rich dust makes HL Tau planet formation fast","Six-wavelength fits favor porous grains for HL Tau cores"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000239,"raw_usage":{"total_tokens":1614,"prompt_tokens":1144,"completion_tokens":470,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":760,"completion_tokens_details":{"reasoning_tokens":400}},"tokens_in":760,"tokens_out":470,"duration_ms":5644,"temperature":1.0,"reasoning_tokens":400,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:57:19.579582+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the radial profile of 0.87 mm scattering polarization in HL Tau at higher signal-to-noise: if the polarization efficiency in the rings falls below the threshold where the adopted prior drops to zero, or varies strongly with radius, the prior in Eq. (15) is invalidated and the porous-dust conclusion must be rederived. Alternatively, a 1.3 cm intensity measurement beyond 40 au that falls outside the factor-of-two separation predicted between the $f_{\\rm AC}=0.3$ and $f_{\\rm AC}=0.8$ models would falsify one of the two compositional branches.","supporting_citations":[],"review_version":1}