{"id":"2bc0369f-1e93-45d3-9bcd-e8ba78afb06d","arxiv_id":"2608.06457","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"The 69 au dust ring in LkCa 15 is broad and thick at 0.88 mm and progressively narrower and thinner at 1.34 and 3.08 mm, implying a massive population of small grains alongside a more concentrated large-grain population.","lead":"This paper measures the geometry of dust rings in a protoplanetary disk by modeling how inclined rings appear in ALMA images, then applies the method to the LkCa 15 disk. It finds that the ring's width, height, and optical depth change with wavelength, implying multiple grain populations with an unexpectedly large mass of small grains.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that no single dust population can explain the multi-band data is not tested against the natural null model: one population with common geometry and a wavelength-dependent optical depth. Since high optical depth itself broadens the apparent width (their Fig.","rationale":"The reader's weakest assumption identifies the Gaussian/isothermal/projection model as the fragile premise, and the reader also flags the absence of an explicit single-population fit. My concern is the sharper version of that flag: the paper's uniqueness claim is not established because the null model is never fit. The geometric method itself is reasonable, and the scattering test in Appendix C is a genuine check, but the physical conclusion depends on rejecting a single-population null model. Because high optical depth is known to broaden the apparent ring in the authors' own model, per-band parameter variation is not by itself evidence for multiple populations. This does not require rejecting the paper's geometric measurements; it requires withholding the mass and multi-population interpretation until the null test is run. Since the reader's verdict is already conditional, the verdict stays unchanged, with the condition sharpened to require this specific test.","tokens_in":18772,"tokens_out":8775,"duration_ms":87934,"concrete_test":"Fit a joint single-population model to the extracted azimuthal brightness and width profiles, or to the nine summary values in Table 1, with common sigma_r/r, sigma_h/r, and T, and tau_perp(lambda) = tau_0 (lambda/lambda_0)^-beta, with free tau_0 and beta (five free parameters). Use the same likelihood (Eq. 4) and the same MCMC machinery. If the minimum chi-squared is acceptable, or the Delta-chi-squared relative to the per-band free-geometry fit is not significant, the central claim fails. As a robustness check, repeat with a single population and a non-Gaussian radial density profile (e.g., super-Gaussian or a two-component Gaussian) and apply the paper's Gaussian-extraction pipeline; if the recovered effective sigma_r/r and sigma_h/r track the observed wavelength trends, the Gaussian assumption is driving the conclusion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline physical conclusion is that the LkCa 15 69 au ring cannot be explained by a single dust population (Abstract and Section 3.3). For this to hold, a global single-population model with one spatial geometry and one opacity law must be unable to reproduce all three bands. The paper never performs such a fit. Instead, it fits tau_perp, sigma_r/r, and sigma_h/r independently in each band (Table 1), observes that the inferred Gaussian parameters vary, and treats this as proof. But the inference from brightness and width profiles to intrinsic geometry is already tau-dependent: their own Fig. 2 shows that increasing optical depth beyond unity broadens the apparent width, and Section 4.1 notes that an optically thin, geometrically thin ring can look as asymmetric as a thick, optically thick ring. A single population with fixed geometry and tau_perp(lambda) proportional to lambda^-beta could therefore produce apparent widths and heights that grow toward shorter wavelengths, exactly the observed trend. The fitted per-band sigma_r and sigma_h may simply be absorbing optical-depth broadening. The absence of an explicit null-model fit means the 'can only be explained' claim, and the derived ~100 Earth-mass small-grain population, are not yet established. The Gaussian-profile assumption aggravates this: for a non-Gaussian density, the effective Gaussian width of the emission changes with optical depth even for one population.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a geometric toy model for an inclined, axisymmetric dust ring with Gaussian radial and vertical density profiles (Eq. 1) and a radiative-transfer intensity prescription that allows for large optical depth (Eq. 2). The central methodological claim is that the azimuthal variation of the peak surface brightness and of the apparent radial width of the ring, after beam convolution, jointly constrain the peak optical depth tau_perp, the fractional radial width sigma_r/r, and the fractional vertical thickness sigma_h/r. The authors apply this pipeline to archival ALMA Band 7, 6, and 3 images of the 69 au ring in LkCa 15, all convolved to a common 60 mas beam, and extract azimuthal brightness and width profiles. The per-band and joint MCMC fits (Table 1) show tau_perp decreasing from about 1.6 at 0.88 mm to 0.35 at 3.08 mm, with sigma_r/r decreasing from 0.121 to 0.078 and sigma_h/r from 0.080 to 0.059. On this basis the authors argue that no single dust population can explain the wavelength-dependent geometry and introduce a two-size model in which broadly distributed small grains (less than about 20 microns) contribute roughly 100 Earth masses and dominate the short-wavelength emission, while more concentrated large grains (greater than about 200 microns) dominate at 3 mm. They discuss implications for settling, stirring, and grain growth and compare with DustPy models.","tokens_in":19093,"tokens_out":10588,"duration_ms":91423,"significance":"If the multi-band morphology is interpreted as intrinsic geometry, the paper makes a significant contribution: it extends the DK21 brightness-asymmetry method to the optically thick regime, adds the radial-width asymmetry as a third observable, and shows that multi-band ALMA data can in principle separate optical depth from geometry without assuming a temperature or opacity law. The modeling is transparent, the MCMC implementation is standard, the use of archival data with a common beam is sensible, and the authors include useful robustness checks including a scattering study (Appendix C) and a DustPy comparison (Appendix D). The measured trends are important regardless of the physical interpretation. However, the headline conclusion that the data 'can only be explained' by multiple grain populations is not yet established, because the paper does not fit the natural null model of a single population with common spatial geometry and a wavelength-dependent optical depth, and the inferred 100 Earth mass small-grain population rests on an approximate effective-Gaussian description of the two-population mixture.","major_comments":[{"comment":"The central claim that a single dust population cannot explain the observations is not tested against the natural null model. In Table 1 the three bands are fitted independently, so sigma_r/r and sigma_h/r are free to vary; a decrease in these fitted values with wavelength is then interpreted in §3.3 as evidence for multiple populations. But this interpretation is valid only if the true density profile is Gaussian, because the model's radiative transfer already includes the optical-depth broadening shown in Fig. 2. For a non-Gaussian density, the effective Gaussian width of the emission from one population with fixed intrinsic geometry and tau_perp proportional to lambda^-beta would grow toward shorter wavelengths, exactly reproducing the observed trend. The authors themselves identify the Gaussian assumption as one that 'needs to be refined' (footnote 2) and list it among caveats (§4.3). To support the wording 'can only be explained by the presence of multiple grain populations' (Abstract and §3.3), the paper should present a global three-band fit with common sigma_r, sigma_h, and T and a parametric opacity law, and it should ideally allow a flexible (non-Gaussian) density profile so that the null hypothesis is fairly tested. A joint Gaussian fit that forces equal sigma in all bands would not be a fair null test, because it would already be rejected by the disjoint per-band uncertainties.","section":"§3.3, Eq. (1), Fig. 2"},{"comment":"The width uncertainty is adopted as one third of the beam FWHM without derivation, and this value directly enters the likelihood in Eq. (4) and the reduced chi-squared in Eq. (6) used to constrain the two-size model. Because the quoted confidence intervals in Table 1 and the significance of the cross-band trend depend on these errors, the choice should be justified quantitatively, for example from the azimuthal scatter of the width measurements, from image-plane bootstrap resampling, or from Monte Carlo realizations of the data. The large size of the observed trends suggests that the qualitative conclusions would survive a different choice, but the quantitative claims, in particular the 100 Earth mass small-grain mass and the chi-squared values in Figure 5, need errors that are not set by an unexplained convention.","section":"§3.1, Eq. (4)"},{"comment":"The effective Gaussian width and height of the two-population mixture are computed as the best Gaussian fit to the sum of two Gaussians, weighted by local optical depth, and these effective values are then compared with the single-Gaussian widths extracted from the ALMA images. This is not the same measurement: the data pipeline fits a Gaussian to the radial brightness profile of the full image at each azimuth after beam convolution, while the model approximation ignores the fact that the radial profile of a two-Gaussian mixture with very different widths is non-Gaussian. The authors acknowledge that the approximation 'may fail when the mixture departs significantly from a single Gaussian' (Section 3.3), which is exactly the regime of their best-fit parameters (sigma_r/r for big grains of about 0.03 versus 0.12 for small grains). The derived mass budget and the chi-squared values in Figure 5 should therefore be validated by generating synthetic images of the two-population model and passing them through the same extraction and fitting pipeline used on the real data.","section":"§3.3, Eq. (6), Fig. 5"}],"minor_comments":[{"comment":"The phrasing 'the minor axis are foreshortened' should read 'the minor axes are foreshortened'.","section":"Abstract"},{"comment":"The phrase 'shinning brightly' should be 'shining brightly'.","section":"§3.3"},{"comment":"The SED comparison is a postdiction rather than a validation; because the model was not fitted to the SED and the real disk is inclined and contains multiple rings, panel (d) should be presented as a consistency check rather than as evidence that the model is correct.","section":"§3.3, Fig. 5d"},{"comment":"The radial fitting window of ±1.5 beam FWHM is reasonable, but the sensitivity of the fitted widths and heights to this choice should be tested, since truncation of the profile can bias Gaussian dispersion estimates.","section":"§3.1"},{"comment":"The posterior values of the calibration factors c6 and c7 introduced in the joint fit are not reported; they should appear in Table 1 or in the corner plot in Figure 7.","section":"§3.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well written and the method is a natural and useful extension of the DK21 approach. My main concern is that the paper's strongest physical claim is built on rejecting a null hypothesis that is never actually fitted; this is fixable and should be fixed before publication. I do not think the paper should be rejected, but the revision needs to include the single-population null test and a more careful error treatment for the width measurements. The two-size mass estimate should be presented with its approximation caveats more prominently."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Colleague],\n\nTwo things to know: the paper gives the field a clean, useful extension of the Doi & Kataoka ring-geometry method, and the LkCa 15 multi-band measurements are a genuinely new dataset. But the abstract's headline—that the ring can only be explained by multiple dust populations—goes a step beyond what the analysis actually shows.\n\nWhat's new and good: the authors extend the DK21 model beyond the optically thin limit, add the azimuthal variation of apparent radial width as an observable that breaks the old degeneracy, and correctly emphasize that even a round beam changes the brightness asymmetry when the ring is narrow. The toy model is transparent, and the MCMC procedure is standard. The data work on the 69 au ring in Bands 7/6/3 is careful: the trends in tau_perp, sigma_r/r, and sigma_h/r all decreasing with wavelength are robust and significant, and they agree with existing width measurements where available. The two-size model is a sensible first attempt to interpret those trends, and the small-grain mass estimate is honestly flagged as opacity-dependent.\n\nWhere it gets soft: the central claim that no single grain population can explain the data is not tested against the obvious null model. They fit tau, width, and height independently in each band, see that the fitted values change, and conclude that a single population is excluded. But a single population with fixed intrinsic geometry and a wavelength-dependent optical depth will also appear broader and thicker at shorter wavelengths, because optical depth itself broadens the emission profile—the paper's own Fig. 2 shows that. To make the multiple-population claim stick, they need to fit a global model with common sigma_r and sigma_h and a free tau(lambda) and show it fails. That computation is straightforward and its absence is conspicuous. The per-band fits are suggestive but not a substitute.\n\nWeaker secondary points: the width uncertainty is set to 1/3 of the beam with no derivation; the two-size model is fitted to the per-band parameters rather than to the images, so the ~100 Earth-mass small-grain population is a model output, not a prediction; and the SED comparison is a postdiction with the caveats they acknowledge. None of these are fatal, but they should be addressed. The Gaussian, isothermal, single-opacity assumptions are standard for a first geometric model; the scattering check is a nice extra, and they flag the Gaussian limitation themselves.\n\nBottom line: this is a serious paper with a new method and a new measurement, and it deserves a full referee. The authors should be pushed on the single-population null fit and on direct fitting of the two-size model to the images. I'd want to see that before the physical conclusion is taken at face value.","headline":"Useful methodological extension and solid multi-band measurements, but the multiple-population conclusion needs a single-population null fit before it is secure.","tokens_in":19680,"tokens_out":5620,"would_cite":true,"duration_ms":49350,"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":"The 69 au ring around LkCa 15 changes shape with observing wavelength, ruling out a single dust population and implying a broad swarm of sub-20-micron grains totaling roughly 100 Earth masses.","keywords":["protoplanetary disks","dust rings","LkCa 15","ALMA","optical depth","dust scale height","multi-wavelength observations","grain growth"],"falsifier":"A high-resolution ALMA Band 9 (0.45 mm) image of the 69 au ring, analyzed with the same azimuthal brightness and width extraction, would discriminate between the scenarios: the small-grain model predicts $\\sigma_r/r \\approx 0.12$, $\\sigma_h/r \\approx 0.08$, and $\\tau_\\perp \\gtrsim 1.6$, close to the Band 7 values, whereas a coagulated centimeter-grain disk would appear narrow and thin, with fractional widths near 0.02 and low optical depth.","tokens_in":18490,"feed_emoji":"🪐","tokens_out":7592,"duration_ms":60195,"temperature":0.7,"pith_summary":"This paper tries to show that the geometry of a modestly inclined dust ring—its radial width, vertical thickness, and optical depth—can be read off from how its apparent brightness and width vary around the ellipse. Applying the geometric model to archival ALMA images of the 69 au ring in LkCa 15, the paper finds that the ring's peak optical depth drops from about 1.6 at 0.88 mm to 0.35 at 3.08 mm, while its fractional radial width falls from 0.121 to 0.078 and its vertical thickness from 0.080 to 0.059. Those wavelength-dependent shapes cannot be produced by a single dust population. The paper argues instead for at least two grain populations: broadly distributed grains smaller than about 20 microns carrying roughly 100 Earth masses, and a more concentrated, less massive population of grains larger than about 200 microns. If correct, this reinterprets fluffy rings as small-grain-dominated rather than strongly stirred, and challenges coagulation models that quickly concentrate dust into large grains.","feed_headline":"LkCa 15's main ring holds 100 Earth masses of tiny dust","feed_subtitle":"ALMA shows the 69 au ring narrowing and thinning with wavelength, evidence of multiple dust populations.","key_machinery":"The load-bearing object is a geometric toy model of an inclined, axisymmetric Gaussian ring with dust density $\\rho = (\\Sigma_0/\\sqrt{2\\pi}\\sigma_h)\\exp[-(r-R)^2/(2\\sigma_r^2) - z^2/(2\\sigma_h^2)]$ at constant temperature, observed through a Gaussian beam. From each synthetic image the paper extracts only two azimuthal curves—peak surface brightness and fitted Gaussian radial width as functions of image-plane angle—and compares these with the same extraction from the ALMA images. Projection makes the ansae brighter because the sightline crosses more material and makes the minor axis narrower by foreshortening, so the amplitude and phase of both curves encode three intrinsic parameters: peak perpendicular optical depth $\\tau_\\perp$, fractional radial width $\\sigma_r/r$, and fractional vertical thickness $\\sigma_h/r$. The apparent radial width is the new diagnostic: it breaks the degeneracy of earlier brightness-only analyses, and the beam shape must be included because convolution dilutes the already-narrow minor axis more severely.","core_discovery":"The central discovery is that the 69 au ring of LkCa 15 is geometrically different at different observing wavelengths, and this difference demands multiple dust populations. Using only the azimuthal variation of peak brightness and apparent radial width, with a Gaussian ring seen at inclination 50 degrees, the paper measures a peak optical depth of $\\tau_\\perp \\approx 1.6$ at 0.88 mm, $1.1$ at 1.34 mm, and $0.35$ at 3.08 mm, while the fractional radial width shrinks from about 0.121 to 0.078 and the vertical thickness from 0.080 to 0.059. A single grain species has one opacity law and one spatial distribution, so it cannot reproduce three different apparent widths and heights. The two-size interpretation gives a small-grain population ($\\lesssim 20\\,\\mu$m) that is radially broad and vertically thick, comparable to or slightly above the gas scale height, with a total mass near 100 Earth masses—about ten times the mass of the spatially concentrated large-grain population ($\\gtrsim 200\\,\\mu$m). Standard dust coagulation simulations, in which grains quickly grow to centimeter sizes and settle, predict a much narrower and thinner ring, so the observed geometry conflicts with the usual expectation.","pith_inferences":["We infer that applying the same geometric method to other multi-band ALMA rings could show whether large small-grain reservoirs are common; if they are, disk evolution models may need larger dust masses and different effective opacities.","A direct test is available at shorter wavelengths: at Band 9 (0.45 mm) the small-grain-dominated model predicts an apparent width and height close to the Band 7 values and a still optically thick ring, whereas a centimeter-grain population would look narrow and thin.","If the inferred 100 Earth masses of small grains is correct, the local dust-to-gas ratio inside the ring may exceed the canonical 0.01, which would affect collision timescales, pressure-bump trapping, and planetesimal formation efficiency; this consequence goes beyond what the paper explicitly computes."],"forward_implications":["When a ring is optically thick at short ALMA bands, brightness asymmetry alone underestimates its height; measuring the azimuthal width variation and solving for optical depth is needed for a reliable height.","Multi-band observations can reveal spatial segregation of grain sizes even below the beam width, because different bands see different effective ring widths and heights.","A ring that appears vertically puffy need not be strongly stirred; it may instead be dominated by small grains that are tightly coupled to the gas.","Coagulation models that quickly convert small grains into centimeter-sized bodies predict much narrower, thinner rings than observed, so the dust-mass budget of such disks must include a large small-grain component."],"supporting_citations":[{"why":"introduced the brightness-asymmetry method that this paper extends to optically thick rings with radial-width information.","marker":"Doi & Kataoka 2021"},{"why":"supplies the calibrated multi-band ALMA images and the ring-width measurements used for the comparison.","marker":"Sierra et al. 2025"},{"why":"provides the original Band 6 and Band 7 observations of LkCa 15.","marker":"Long et al. 2022"},{"why":"gives the literature dust-height measurement that the paper's Band 6 result is checked against.","marker":"Villenave et al. 2025"},{"why":"reports a much thinner Band 6 ring, which the paper explains as a consequence of fixing the optical depth too low.","marker":"Jiang et al. 2025"},{"why":"supplies the coagulation simulation showing that grains quickly grow to centimeter sizes, the expectation the observations contradict.","marker":"Stammler & Birnstiel 2022"},{"why":"provides the standard prediction that dust mass quickly concentrates into large grains.","marker":"Birnstiel et al. 2012"},{"why":"defines the DIANA dust composition used for the opacity calculations in the two-size model.","marker":"Woitke et al. 2016"},{"why":"offers the Band 9 image showing a broad ring, independent evidence for abundant small grains.","marker":"Leemker et al. 2022"},{"why":"constrains scattered-light grain sizes and masses, supporting the inferred small-grain population.","marker":"Swastik et al. 2026"}],"fun_headline_variants":["LkCa 15 dust ring hides 100 Earth masses of tiny grains","Wavelength reveals dual dust populations in LkCa 15 ring","ALMA shows LkCa 15 ring morphs with wavelength: multiple dusts","Tiny dust outweighs big grains 10-to-1 in LkCa 15's ring"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the ring is axisymmetric, has a constant temperature, and has a bell-shaped (Gaussian) dust distribution in both radius and height with a single dust opacity per band; non-Gaussian vertical structure, intrinsic brightness asymmetries such as shadows or a hot wall, or a radial temperature gradient would bias the inferred optical depths, widths, and heights.","fun_headline_variants_meta":{"raw":{"variants":["LkCa 15 dust ring hides 100 Earth masses of tiny grains","Wavelength reveals dual dust populations in LkCa 15 ring","ALMA shows LkCa 15 ring morphs with wavelength: multiple dusts","Tiny dust outweighs big grains 10-to-1 in LkCa 15's ring"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000592,"raw_usage":{"total_tokens":2850,"prompt_tokens":1096,"completion_tokens":1754,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":712,"completion_tokens_details":{"reasoning_tokens":1666}},"tokens_in":712,"tokens_out":1754,"duration_ms":11738,"temperature":1.0,"reasoning_tokens":1666,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:34:53.512042+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-resolution ALMA Band 9 (0.45 mm) image of the 69 au ring, analyzed with the same azimuthal brightness and width extraction, would discriminate between the scenarios: the small-grain model predicts $\\sigma_r/r \\approx 0.12$, $\\sigma_h/r \\approx 0.08$, and $\\tau_\\perp \\gtrsim 1.6$, close to the Band 7 values, whereas a coagulated centimeter-grain disk would appear narrow and thin, with fractional widths near 0.02 and low optical depth.","supporting_citations":[],"review_version":1}