{"id":"2ead2ca8-3fb9-47e2-b6a0-b248144f8c9f","arxiv_id":"2506.11726","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":12,"one_line_summary":"ALMA imaging resolves HD 138965's outer debris belt to a radius of 150 au and a width of 49 au, favoring astronomical silicate dust and limiting any inner companion to roughly 2.3 Jupiter masses.","lead":"Using ALMA millimetre-wavelength observations, astronomers spatially resolved the outer debris belt around the young star HD 138965, finding a broad ring at 150 au with a width of 49 au. The new geometry tightens limits on dust composition, scattering albedo, and the mass of any hidden planetary companion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported outer-belt width may be inflated by the unconstrained vertical scale height and the assumed single-Gaussian profile; a visibility refit with fixed h extremes and an alternative radial profile would test whether ΔR/R = 0.77 is robust.","rationale":"The reader's conditional verdict and weakest assumption focus on the fixed inner-belt parameters biasing the dust composition comparison. That is a real and well-identified concern. My stress-test pass targets an earlier link in the chain: the claimed resolved width of the outer belt itself. The ALMA visibility modelling is standard, and the agreement of the fitted radius with the Herschel measurement gives independent support to the radius; no objection is raised to R_peak. However, the width is only marginally resolved and is degenerate with the completely unconstrained vertical scale height and with the assumed Gaussian surface-density profile. Since the subsequent companion and composition constraints are evaluated at fixed radius and width from this fit, the width assumption is at least as load-bearing as the inner-belt assumption, and it is more direct to verify from the same ALMA data. The reader's concern would remain even if the width is robust, so I regard the two concerns as complementary rather than identical, hence 'partial'. The proposed test can be run with existing calibrated visibilities and standard tools, and the result would either retire or substantiate the concern without requiring new observations. The paper honestly flags the scale-height limitation and the marginal resolution, so this is a matter of systematic uncertainty rather than an internal inconsistency; a conditional acceptance with the refit as a check remains the appropriate disposition.","tokens_in":23423,"tokens_out":6871,"duration_ms":80823,"concrete_test":"Refit the calibrated ALMA visibilities with the same RADMC-3D/Galario pipeline, but (a) repeat the fit with h fixed to the lower and upper ends of the prior range (e.g., 0.01 and 0.30), and (b) replace the Gaussian radial profile with a radially sharp-edged power-law belt having inner radius R_in and outer radius R_out as free parameters. Compare R_peak, σ_R (or R_in/R_out), and ΔBIC among the models. If σ_R shifts by more than the quoted 6 au uncertainty between the h-fixed cases, or if the sharp-edged model fits comparably well with a different radial extent, the reported ΔR/R and all downstream quantities derived from the inner edge at 101 au (companion mass and fixed Stardust outer-belt width) must be treated as model-dependent rather than robustly measured.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the outer belt is resolved in width (σ_R = 49+7/−6 au, ΔR/R = 0.77, Table 2) is obtained from a single-Gaussian belt model, while the same fit leaves the vertical scale height h completely unconstrained: the h posterior is flat and uninformative (Appendix B), and Section 3.2 explicitly states that a broad range of scale heights fits the observations. At the fitted inclination of about 50°, a vertically thick belt can broaden the projected emission in ways that a Gaussian radial profile may absorb into σ_R, so the quoted width is not established independently of the vertical and radial profile assumptions. The 2σ NW/SE residual asymmetry in Figure 1 further indicates that the symmetric single-Gaussian model is not uniquely capturing the surface brightness. This matters because the companion constraint uses the inner edge R − σ_R = 101 au (Section 4), and the SED/composition modelling fixes σ_R = 49 au as an input (Section 3.4). If the true radial profile is sharper-edged or vertically different, the propagated inner edge, companion mass limit, and dust composition ranking would each shift. The radius itself is on firmer ground, since it agrees with the Herschel-derived extent, but the width is the part of the geometry on which the subsequent quantitative conclusions lean most heavily.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents ALMA Band 6 millimetre-wavelength imaging of HD 138965, a young A star in the Argus association, and models the debris disc as a single Gaussian belt in the visibility plane. The authors report a resolved outer belt with peak radius 150^{+10}_{-7} au, width sigma_R = 49^{+7}_{-6} au (Delta R/R = 0.77), inclination 49.9 degrees, and position angle 173.3 degrees, with no constraint on the vertical scale height. They combine this geometry with an SED model (Stardust) to infer that astronomical silicate is the preferred outer-belt dust composition, while 10% water-ice inclusions cannot be excluded; from the HST non-detection they derive a scattering albedo limit omega <= 0.09; and from the belt architecture they place a companion mass limit of 2.3 +/- 0.4 M_Jup interior to about 78 au. They also report a non-detection of CO(2-1) emission.","tokens_in":23755,"tokens_out":9883,"duration_ms":93610,"significance":"If the reported geometry holds, this is the first millimetre-wavelength spatial resolution of the outer belt of HD 138965 and it anchors an independent dust-composition ranking, an albedo upper limit, and a companion-mass constraint that improves on direct imaging limits. The paper uses a standard and largely reproducible methodology (MCMC visibility fitting with RADMC-3D and Galario, BIC model comparison) and is appropriately explicit about several limiting assumptions, notably the unconstrained scale height and the degeneracy in the inner-belt parameters. The results are incremental rather than transformative, but they are a solid observational contribution to the debris-disc literature and the strengths of the analysis are the clear separation of the measurement (ALMA geometry) from the derived constraints (composition, albedo, companion mass).","major_comments":[{"comment":"The quoted width sigma_R = 49^{+7}_{-6} au is not robust against the unconstrained vertical scale height and the assumed single-Gaussian radial profile. The scale-height posterior is flat and uninformative (Appendix B), and Section 3.2 states that a broad range of scale heights fits the observations. At i ~= 50 degrees, a vertically thick belt can broaden the projected emission, and this broadening can be absorbed by sigma_R in a Gaussian model. Since sigma_R is subsequently used to set the outer-belt inner edge at R - sigma_R = 101 au for the companion limit (Section 4) and is fixed as an input to the SED modelling (Section 3.4), the width needs a robustness test: repeat the visibility fit with h fixed at the extremes (for example 0.01 and 0.30) and with an alternative radial profile (for example a sharp-edged ring or a power-law surface density) to establish that Delta R/R = 0.77 is not an artefact of the model family. The 2-sigma NW/SE residual asymmetry in Figure 1 reinforces the need for such a test.","section":"3.2, Table 2, Appendix B"},{"comment":"The outer-belt composition ranking is conditional on fixed inner-belt parameters (r_m = 13.7 au, sigma_r = 1.75 au, s_min = 12 micron, q = 3.5, M_s = 0.154 x 10^-3 M_Earth) that were obtained from an initial Stardust fit assuming a pure-silicate outer belt. The text explicitly acknowledges that this may bias the outer-belt comparison toward astronomical silicate, but the inner-belt parameters are not marginalised over or varied in the composition fits. The Delta BIC margins separating S100 from S90WC10 (7.0) and from S90WA10 (4.4) are modest, so a plausible change in the inner-belt parameters could alter the ranking. I request a sensitivity test in which the inner-belt radius, s_min, q, or mass are varied within their posterior ranges, or a joint fit, to show that the silicate preference is not an artefact of the initial assumption.","section":"3.4, Table 4"}],"minor_comments":[{"comment":"The disc flux density is reported as 1.460 +/- 0.230 mJy in Table 2 but as 1.46^{+0.08}_{-0.07} mJy in Figure B1; if the larger uncertainty includes the weather-related calibration systematic described in Section 2.1, this should be stated explicitly, since the SED fit uses the +/- 0.23 mJy value.","section":"3.2, Table 2 vs Figure B1"},{"comment":"The initial MCMC analysis is reported in the text to give an inner-belt mean radius of 15^{+3}_{-2} au, while Table 4 lists r_m = 13.7 au; please state which value (posterior median or maximum likelihood) was adopted for the fixed inner-belt model.","section":"3.4, Table 4"},{"comment":"The sentence in Section 3.4 saying that the best-fit results correspond to the porosity scenarios and to a strong candidate with 90:10 crystalline water ice appears to contradict the BIC table and the abstract, which identify pure astronomical silicate as the best fit; please reconcile the wording.","section":"3.4, Table 5, Abstract"},{"comment":"The companion mass limit is quoted as a <= 78 au in the abstract but as a <= 74^{+8}_{-6} au in Section 4, and the improvement factor is stated as a factor of two in the abstract but as a factor of five over Matthews et al. (2018) and a factor of two over the SHARDDS image in Section 4; unify these numbers and specify the baseline.","section":"Abstract, Section 4"},{"comment":"The SED-derived albedo of 0.56 is computed for 'dirty ice' grains at 1.6 micron, whereas the preferred composition in Table 5 is pure astronomical silicate and the HST limit is at 0.6 micron; recompute the comparison at the HST wavelength for the best-fit composition, or present it explicitly as an illustrative consistency check rather than as a tension.","section":"3.5, Table 5"},{"comment":"The caption says the posterior distributions are based on '10,0000 realisations', which appears to be a typo, and Figure B1 omits the scale-height posterior; since h is unconstrained, please show at least the h-sigma_R covariance or state explicitly that it was not stored.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"This is a competent observational paper whose main conclusions are conditional on two robustness checks: the width measurement under vertical/profile degeneracies and the composition ranking under inner-belt assumptions. Both are fixable with additional modelling and do not appear to be fundamental errors. I would also ask the authors to clean up the inconsistency in the flux-density uncertainties between Table 2 and Figure B1, as the quoted calibration uncertainty is important for the SED constraints. The paper is within the scope of MNRAS."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a solid, incremental ALMA debris disk paper. The genuinely new result is the first spatially resolved mm-wavelength image of HD 138965's outer belt, giving a peak radius of 150 au and a broad width (sigma 49 au, ΔR/R=0.77) at ~50 degrees inclination. The radius is on firm ground—it matches the Herschel extent—and the inclination and position angle are the most precise yet. The paper also derives limits on dust composition, albedo, and a possible sculpting companion. No new methodology, but a useful data point for the sample of radially resolved debris disks.\n\nWhat it does well: the visibility fitting is standard but honestly reported, including the unconstrained vertical scale height and the weather-affected flux calibration. The authors also clearly flag that fixing the unresolved inner belt parameters biases the outer belt composition ranking toward astronomical silicate. That candor is welcome.\n\nSoft spots, in proportion. The width is only marginally resolved and is degenerate with the vertical scale height, which the data do not constrain. At the fitted inclination a vertically thick belt can broaden the projected emission; the quoted ΔR/R should be treated as model-dependent. The companion mass limit uses R − sigma as the inner edge, so this propagates. A refit with fixed h extremes and an alternative radial profile would test how robust the width is. The composition ranking rests on inner belt parameters from an initial fit that assumed a pure-silicate outer belt; the authors acknowledge but do not marginalize this. The BIC values do not rule out 10% water-ice inclusions, so the text calling a 90:10 crystalline water-ice mixture a 'strong candidate' overstates their own numbers. And the 2.3 M_Jup companion 'mass limit' is actually a lower limit from stirring models, not an upper limit; the abstract should say so. The stellar age is uncertain by an order of magnitude, partly mitigated by giving limits at several ages.\n\nBottom line: this deserves a serious referee. The ALMA measurement is credible and the limitations are mostly acknowledged. I would accept it after the width is discussed more cautiously, the companion mass is unambiguously a lower limit, and the 90:10 water-ice overclaim is corrected.","headline":"Credible first resolved mm view of HD 138965's outer belt; radius robust, width and composition model-dependent, but worth refereeing with fixes.","tokens_in":24353,"tokens_out":5964,"would_cite":true,"duration_ms":54126,"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":"First ALMA millimetre imaging of HD 138965 resolves its cool outer debris belt at 150 au and, combined with radiative-transfer modelling, limits any companion to 2.3 Jupiter masses.","keywords":["debris discs","circumstellar matter","radio continuum: planetary systems","planet-disc interactions","stars: individual: HD 138965","millimetre astronomy","dust composition","radiative transfer"],"falsifier":"Spatially resolving the inner belt would settle the question: JWST/MIRI imaging at about $10$ $\\mu$m or longer-baseline ALMA at 1.3 mm could reach the ~0.25 arcsec separation of a 20 au inner belt, and showing its radius, width, or mass differs from the assumed values would directly test whether the outer-belt silicate preference is an artefact of the fixed inner-belt model. A second check is detecting water-ice or carbon spectral features in the outer belt's SED, since the model currently rules out 50% ice mixes and all tested carbon mixes.","tokens_in":2021,"feed_emoji":"🪐","tokens_out":9941,"duration_ms":143007,"temperature":0.7,"pith_summary":"The paper presents the first ALMA 1.3 mm observations to spatially resolve the cool outer debris belt around the young A star HD 138965. The belt peaks at $150^{+10}_{-7}$ au and is unusually broad, with a width of $49^{+7}_{-6}$ au, giving $\\Delta R/R = 0.77$. The resolved architecture is then combined with a radiative-transfer model that moves dust grains under gravity, radiation pressure, and Poynting–Robertson drag; the best-fitting dust composition is astronomical silicate, while mixtures with 10% water ice remain plausible and 50% ice or any carbon mix is excluded. The same geometry yields an optical albedo limit of $\\omega \\le 0.09$ and a companion mass limit of $2.3 \\pm 0.4$ Jupiter masses interior to 78 au. These results matter because a resolved debris belt's width and emission are among the few observable probes of unseen planets and the material composition of planet-forming regions.","feed_headline":"ALMA resolves a 150 au debris belt and caps planet mass","feed_subtitle":"First millimetre-size map of the belt sets a 2.3-Jupiter-mass limit on unseen companions inside 78 au.","key_machinery":"The argument runs through three pieces of machinery. First, a Gaussian-belt model of the ALMA visibilities, generated with a radiative-transfer code and sampled by a Markov chain Monte Carlo ensemble, yields the belt's radius, width, and orientation. Second, the paper's Stardust model redistributes dust grains under gravity, radiation pressure, Poynting–Robertson drag, and a simple collisional lifetime, predicting the spectral energy distribution for each grain composition, with model comparison by the Bayesian Information Criterion. Third, a literature single-planet sculpting relation converts the measured belt inner edge into a companion mass limit, and an established albedo formula turns the HST scattered-light non-detection into $\\omega \\le 0.09$.","core_discovery":"The central claim is that HD 138965's outer debris belt is spatially resolved at millimetre wavelengths for the first time, with a peak radius of $150^{+10}_{-7}$ au and a Gaussian width of $49^{+7}_{-6}$ au (fractional width $\\Delta R/R = 0.77$), at an inclination of $49.9^{+3.3}_{-3.7}$ degrees. From this resolved structure, the authors find that astronomical silicate is the best-fitting dust composition for the outer belt among the tested inclusion mixtures, while scenarios with at least 10% crystalline or amorphous water ice cannot be rejected and 50% ice or any carbon mix is eliminated. Combining the ALMA image with the HST optical non-detection, they derive the scattering albedo upper limit $\\omega \\le 0.09$; and using a single-planet sculpting model on the belt's inner edge at 101 au, they place a companion mass limit of $2.3 \\pm 0.4$ Jupiter masses for separations $a \\le 78$ au.","pith_inferences":["If the star's true age is closer to 350 Myr rather than the assumed Argus age of 31 Myr, the direct-imaging mass limits at 80 au weaken to roughly 15 Jupiter masses, which would make the architecture-based limit of 2.3 Jupiter masses the only meaningful constraint.","If future JWST/MIRI imaging places the inner belt at a radius different from the assumed 13.7 au, then the inner-belt grain temperature and mass would need to change to fit the same mid-infrared excess, which could shift the outer-belt composition preference.","The residual image shows a 2σ brightness asymmetry between the two ansae; if real, this could indicate an eccentric ring or pericentre glow, which deeper or higher-resolution ALMA imaging could confirm.","The discrepancy between the HST albedo limit (≤0.09) and the SED-inferred albedo (~0.56 for the smallest grains) suggests the outer-belt grains may be porous, icy, or strongly forward-scattering, which future polarimetric observations could test."],"forward_implications":["The outer belt's fractional width of 0.77 places HD 138965 among the most radially extended resolved debris discs, comparable to HR 8799's belt, making it a candidate for hidden substructure that current data cannot resolve.","The companion mass limit of $2.3 \\pm 0.4$ Jupiter masses interior to 78 au is about a factor of two tighter than the deepest direct-imaging limit at that separation, showing the power of belt morphology.","The non-detection in scattered light, combined with the resolved ALMA shape, yields an optical albedo upper limit of 0.09, which implies the disc would become detectable with roughly a factor of four deeper imaging.","The CO(2–1) upper limit of $2.7 \\times 10^{-23}$ W m$^{-2}$ is about three times the level predicted for a gas mass of $8.6 \\times 10^{-7}$ Earth masses, so gas is not ruled out but at present is unconstrained."],"supporting_citations":[{"why":"Provides the Herschel-derived disc radius and a previous 50:50 silicate:ice composition that the new modelling refines.","marker":"Morales et al. 2016"},{"why":"Gives the prior high-contrast imaging mass limits that the architecture-based limit improves upon.","marker":"Matthews et al. 2018"},{"why":"Supplies the single-planet sculpting relation used to derive the 2.3 Jupiter-mass limit from the belt's inner edge.","marker":"Pearce et al. 2022"},{"why":"Provides the self-stirring model used to interpret the belt's broad architecture and to estimate the required planetesimal mass.","marker":"Krivov & Booth 2018"},{"why":"Gives the steady-state collisional cascade size distribution index q=3.5 adopted for the inner belt.","marker":"Dohnanyi 1969"},{"why":"Supplies the astronomical silicate optical constants used in all grain composition mixtures.","marker":"Draine 2003"},{"why":"Provides the emcee Markov chain Monte Carlo sampler used for all Bayesian parameter fits.","marker":"Foreman-Mackey et al. 2013"},{"why":"Gives the equation that converts the scattered-light non-detection and disc geometry into an albedo upper limit.","marker":"Marshall et al. 2018"}],"fun_headline_variants":["ALMA resolves 150 au debris belt, caps planet mass","Millimetre imaging of HD 138965: silicate dust, 2.3-Jupiter limit","Debris disc at 150 au: ALMA sets planet and albedo constraints","ALMA map of young star's belt narrows planet and dust answers"],"cache_read_input_tokens":26368,"weakest_assumption_plain":"The ranking of outer-belt dust compositions depends on fixing the unresolved inner belt's radius ($13.7$ au), width ($1.75$ au), minimum grain size ($12$ $\\mu$m), size distribution index ($q = 3.5$), and mass ($0.154 \\times 10^{-3}$ Earth masses) from an initial Stardust fit that assumed the outer belt is pure astronomical silicate; if those inner-belt values are different, the BIC preference for silicate could change.","fun_headline_variants_meta":{"raw":{"variants":["ALMA resolves 150 au debris belt, caps planet mass","Millimetre imaging of HD 138965: silicate dust, 2.3-Jupiter limit","Debris disc at 150 au: ALMA sets planet and albedo constraints","ALMA map of young star's belt narrows planet and dust answers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000207,"raw_usage":{"total_tokens":1476,"prompt_tokens":1098,"completion_tokens":378,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":714,"completion_tokens_details":{"reasoning_tokens":293}},"tokens_in":714,"tokens_out":378,"duration_ms":4882,"temperature":1.0,"reasoning_tokens":293,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:04:47.115209+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Spatially resolving the inner belt would settle the question: JWST/MIRI imaging at about $10$ $\\mu$m or longer-baseline ALMA at 1.3 mm could reach the ~0.25 arcsec separation of a 20 au inner belt, and showing its radius, width, or mass differs from the assumed values would directly test whether the outer-belt silicate preference is an artefact of the fixed inner-belt model. A second check is detecting water-ice or carbon spectral features in the outer belt's SED, since the model currently rules out 50% ice mixes and all tested carbon mixes.","supporting_citations":[{"cited_title":"Y., Bryden G., Werner M","cited_arxiv_id":null,"evidence_quote":"Provides the Herschel-derived disc radius and a previous 50:50 silicate:ice composition that the new modelling refines."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the prior high-contrast imaging mass limits that the architecture-based limit improves upon."},{"cited_title":"P., Milli J., Choquet \\'E ., del Burgo C., Kennedy G","cited_arxiv_id":null,"evidence_quote":"Gives the equation that converts the scattered-light non-detection and disc geometry into an albedo upper limit."}],"review_version":1}