{"id":"831b8802-c4ab-4685-b6de-b6756af0c1cf","arxiv_id":"2505.20070","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"DRAGyS automatically fits ellipses to ring brightness peaks in protoplanetary disk images, estimates disk geometry, and extracts limb brightening corrected scattering phase functions without radiative transfer modeling.","lead":"DRAGyS is a new tool that measures the geometry of ring-shaped protoplanetary disks directly from telescope images and extracts the scattering phase function, the curve showing how dust grains scatter light. It adds a correction for the limb brightening effect that can change the curve by up to about 30 percent.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 3 limb brightening correction vanishes for χ≈1 yet the paper claims 10–30% corrections; the printed formula cannot reproduce the central result.","rationale":"The reader's CONDITIONAL verdict is justified, but the most load-bearing issue is not the circularity/tau=1 assumption (which the paper acknowledges as a limitation) but the limb brightening formula itself. Eq. 3, as printed, is mathematically inconsistent with the stated assumptions and with the reported correction amplitudes. This matters because the LB correction is the paper's main novelty: the claim that 'the limb brightening-corrected SPF is a much better match to the expected SPF' rests entirely on that formula. If the formula is wrong or mis-transcribed, the validation in Fig. 4 is not reproducible from the manuscript. The concrete test—inspecting the GitHub code and comparing with Tazaki et al. 2023—would settle whether the paper text needs a correction or whether the method itself has a bug. This is a more fundamental problem than the missing uncertainty envelopes: it affects the core method, not only its presentation. I therefore maintain the CONDITIONAL verdict (no change to the reader's verdict category) but with an additional, explicit condition: the LB formula in Eq. 3 must be corrected and verified against the code and the original Tazaki et al. derivation.","tokens_in":19442,"tokens_out":13535,"duration_ms":161290,"concrete_test":"Inspect the DRAGyS GitHub repository (https://github.com/mroumesy/DRAGyS) and identify the exact expression used for the limb brightening correction. Compare it with Eq. 3 of the paper and with the derivation in Tazaki et al. (2023). Then evaluate the implemented formula for χ = 1.0, 1.1, 1.2 and for the h/r values of the observed disks in Table 2; check whether it produces correction amplitudes of ≈10–30% claimed in Section 4. If the code differs from Eq. 3 or yields negligible corrections for χ≈1, the paper text must be corrected to match the implementation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.3 defines the limb brightening correction as Eq. 3: LB = cos(γ) sin(γ′−γ) / [cos(γ) sin(γ′−γ) + cos(i) cos(γ′−γ) − sin(γ′) cos(Ψ)], with γ = h/r and γ′ = dh/dr = χ h/r. The text states that a fixed flaring exponent 'close to unity' is adopted, and Fig. 4 shows a substantial LB correction; Section 4 reports corrections of up to ≈30%. However, if χ = 1, then γ′ = γ, sin(γ′−γ) = 0, and Eq. 3 gives LB = 0 identically. For χ = 1.00001 (the value plotted in Fig. 5), sin(γ′−γ) ≈ 10⁻⁶, so LB remains effectively zero unless cos(i) − sin(γ) cos(Ψ) is also of order 10⁻⁶, which is not the case. Even for χ = 1.1 and γ = 0.15, the formula yields LB ≲ 8%, far below the reported ≈30%. Thus the printed equation cannot produce the claimed correction amplitudes. Either the implemented expression differs from Eq. 3, or the adopted χ is much larger than 'close to unity'. In both cases the manuscript is internally inconsistent, and the LB-corrected SPF matches in Fig. 4 are not reproducible from the text alone.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"DRAGyS is a purely geometric pipeline for ring-shaped protoplanetary disks: it detects intensity peaks, fits an ellipse, converts the ellipse parameters into inclination, position angle, and scattering-surface aspect ratio h/r, and extracts total and polarized scattering phase functions. It also applies a limb-brightening correction following Tazaki et al. (2023). The method is validated on noiseless MCFOST images with known geometry and intrinsic SPF, then applied to nine archival SPHERE/IRDIS polarized-intensity images of six disks. The paper claims recovery of inclination within about 2 degrees, position angle within about 1.5 degrees, SPFs consistent with the reference Diskmap method to 2-18 percent, and limb-brightening corrections that can reach about 30 percent.","tokens_in":19684,"tokens_out":15521,"duration_ms":163867,"significance":"The tool addresses a real need: fast SPF extraction for large samples without per-target radiative transfer modeling. The synthetic validation is well designed because the input geometry and intrinsic SPF are known from MCFOST, and the authors make the code publicly available. The sector-based extraction for asymmetric disks is a useful feature. If the limb-brightening correction is correctly implemented, the tool would be a valuable complement to Diskmap. However, the printed limb-brightening formula and the assumed flaring exponent are internally inconsistent, so the quantitative limb-brightening result is not reproducible as written.","major_comments":[{"comment":"As written, the limb-brightening correction is inconsistent with the text. Setting chi=1 gives gamma'=gamma, hence sin(gamma'-gamma)=0 and LB=0 identically; for chi=1.00001 the correction is numerically negligible. The text says a fixed flaring exponent 'close to unity' is adopted, yet Section 4 reports corrections of about 10 percent on average and up to about 30 percent, and Fig. 4 shows a substantial correction. Concretely, for i=70 deg, h/r=0.15, chi=1.1, Eq. (3) gives at most about 8 percent; for the real geometries in Table 2 the correction stays below about 10 percent even for chi=1.2. A 30 percent correction would require chi well above the values shown in Fig. 5, which are 1.00001, 1.1, and 1.2. The adopted chi value must be stated explicitly, or the equation must be corrected, and the corrected SPFs in Figs. 4 and 7 must be reproducible from the text and the code.","section":"Section 3.3, Eq. (3)"},{"comment":"The limb-brightening correction and the scattering-angle mapping assume that the detected ring traces the tau_s=1 surface and that h(r) follows Eq. (1) with a fixed flaring exponent. Appendix D shows that for low dust masses (10^-6 and 10^-7 M_sun per ring) the surface measured by DRAGyS lies below the tau_s=1 surface, and no test in the paper varies the true flaring exponent of the input model. Since the correction depends on chi, the paper should state the range of disk parameters for which the assumption holds, demonstrate that the corrected SPF is not biased by an incorrect chi, or provide a way to estimate chi from the data.","section":"Section 3.3 and Appendix D"}],"minor_comments":[{"comment":"The curve with chi=1.00001 is effectively identical to unity under Eq. (3) and adds no information; it should be replaced by a realistic chi value or removed.","section":"Fig. 5"},{"comment":"The percentage offsets between SPFs are given without error bars; the uncertainty propagation described in Section 2.3 should be applied to these numbers so the reader can judge whether the differences are significant.","section":"Fig. 7"},{"comment":"The position-angle convention is described as measured from the horizontal left axis counter-clockwise; because PA is normally quoted from north to east, please clarify how this convention relates to the values in Table 2 and to the reference values.","section":"Section 2.2"},{"comment":"The text says the sample is six disks observed in J and/or H bands, while the abstract says 'nine images for six disks'; please make the count explicit in the main text to avoid confusion.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of A&A, and the public code and synthetic validation are strengths. The main technical obstacle is the limb-brightening formula/flaring-exponent inconsistency; if the authors can provide the actual formula or the adopted chi and verify the figures, I would support publication. No concerns about citation practices or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things up front. The tool is genuinely useful and the synthetic validation is well done. But the limb brightening correction, which is a headline feature, does not reproduce from the printed Eq. 3. For χ=1 (and the text says they adopt a flaring exponent \"close to unity\"), LB=0 identically, yet the paper reports 10–30% corrections and Fig. 5 shows large effects for χ=1.00001. Even with χ=1.1 and h/r=0.15, Eq. 3 gives at most a few percent. The stress-test note is right: the printed formula cannot generate the figures. This needs to be fixed before the method is published as is.\n\nWhat earns credit: the pipeline automates ellipse-based geometry estimation, SPF extraction, and limb brightening correction in one open-source package, avoiding radiative transfer modeling. The MCFOST validation is solid: geometry recovered within ~2° in inclination, ~1.5° in PA, and the LB-corrected SPF tracks the intrinsic curve. The application to archival SPHERE data shows 2–18% agreement with Diskmap on good cases, and sector extraction is a sensible way to handle asymmetries. The paper is clearly written and honest about multiple scattering and the low-mass caveat in Appendix D.\n\nSoft spots, in proportion: the Eq. 3 issue is the load-bearing one. After that, the synthetic tests are noiseless, so the quoted uncertainties are optimistic. The SPF comparison figures lack uncertainty envelopes, making the 2–18% and the ~30% LB claims hard to evaluate. The geometry assumptions—circular disk, single power-law flaring surface with χ near unity—are reasonable for ringed disks but should be flagged as limitations, which the paper partly does.\n\nThis is a paper for disk observers who want fast, model-free SPF extraction. It deserves a serious referee, but only after the authors correct Eq. 3 or explain how the LB correction was actually computed. As written, someone cannot reproduce the central numerical claim.","headline":"Useful pipeline with solid synthetic validation, but the limb-brightening formula as printed cannot produce the claimed 10–30% corrections; that needs fixing before this is publishable as is.","tokens_in":20323,"tokens_out":4280,"would_cite":false,"duration_ms":43591,"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":"DRAGyS claims that the geometry of ring-shaped protoplanetary disks can be estimated directly from scattered-light images by ellipse-fitting brightness peaks, and that the scattering phase function can then be extracted and corrected for…","keywords":["protoplanetary disks","scattering phase function","polarized intensity","limb brightening","disk geometry","ring fitting","high-contrast imaging","near-infrared polarimetry"],"falsifier":"Run DRAGyS on synthetic images with known intrinsic SPFs while varying the flaring exponent (for example $\\chi$ = 0.5 or 1.5) or using an eccentric ring; if the limb-brightening-corrected SPF departs from the intrinsic one by more than the few-percent recovery claimed for circular, $\\chi\\approx1$ disks, the assumed single circular scattering surface is falsified. The same test applies to low-dust-mass models, where the brightness peak no longer traces the effective scattering surface, a case the paper itself shows in its Appendix D.","tokens_in":19194,"feed_emoji":"🔭","tokens_out":7311,"duration_ms":64890,"temperature":0.7,"pith_summary":"The paper introduces DRAGyS, a tool that claims to recover the scattering phase function (SPF) of ring-shaped protoplanetary disks directly from scattered-light images, without fitting radiative transfer models. The tool fits ellipses to ring brightness peaks to get disk inclination, position angle, and aspect ratio, then uses that geometry to map image pixels to scattering angles and to correct a geometric limb-brightening bias. On synthetic images, the corrected SPF closely matches the intrinsic dust SPF, and on archival polarized-intensity images of six disks the recovered geometry and SPF agree with previously published results to within a few degrees and 2–18 percent. The paper argues this fast geometry-only route makes large-sample SPF surveys feasible and removes a bottleneck when inferring dust grain properties.","feed_headline":"Tool pulls scattering phase function directly from disk images","feed_subtitle":"Fits ring ellipses from images alone, then corrects a near-side brightness bias that can reach 30 percent.","key_machinery":"The mechanical core is the assumption that the brightness peaks of a ring outline the ellipse of the disk scattering surface: least-squares ellipse fitting gives the semi-major axis $M$, semi-minor axis $m$, ellipse-center offset $D$, and position angle, from which the inclination follows as $i=\\cos^{-1}(m/M)$ and the aspect ratio as $h/r=D/(M\\sin i)$ under the flaring law $h(r)=h_{\\rm ref}(r/r_{\\rm ref})^{\\chi}$. The scattering angle at each pixel is then computed with $\\cos(\\Psi)=\\cos(\\gamma)\\cos(\\phi)\\sin(i)+\\sin(\\gamma)\\cos(i)$, where $\\gamma=h/r$, and the limb-brightening correction applies the ratio $\\mathrm{LB} = \\frac{\\cos\\gamma\\,\\sin(\\gamma'-\\gamma)}{\\cos\\gamma\\,\\sin(\\gamma'-\\gamma)+\\cos i\\,\\cos(\\gamma'-\\gamma)-\\sin\\gamma'\\,\\cos\\Psi}$ with $\\gamma'=\\chi h/r$. This formula carries the geometric correction that turns the observed near-side-brightened SPF into the intrinsic one, with uncertainties propagated through the fitted parameters.","core_discovery":"The central claim is that a disk's scattering phase function can be extracted from an image through a purely geometric pipeline that assumes only that the disk is circular and that the brightness maxima of a ring trace the disk scattering surface. Fitting an ellipse to those maxima yields the inclination and position angle directly, while the offset between the ellipse center and the star gives the scattering surface height under the power-law flaring model $h(r)=h_{\\rm ref}(r/r_{\\rm ref})^{\\chi}$. The scattering angle at each pixel follows from the projected geometry, and a published limb-brightening formula, evaluated with the fitted aspect ratio and an assumed flaring exponent near unity, corrects the extracted SPF for the near-side/far-side brightness asymmetry. The paper reports that the limb-brightening-corrected SPF is a much better match to the expected SPF from the model, and that on real observations the differences with reference SPFs are between 2 and 18 percent, while the limb-brightening effect itself can shift the SPF by up to about 30 percent.","pith_inferences":["A natural extension the paper does not quantify is to run the same pipeline on simulated eccentric rings and measure at what eccentricity the extracted SPF departs from the intrinsic curve.","The reported weak dependence of limb brightening on the flaring exponent suggests a two-parameter variant could fit $\\chi$ itself rather than fixing it near unity, which might reduce bias in strongly flared disks.","The same ellipse-plus-offset machinery could plausibly be applied to debris disks or other scattered-light nebulae where brightness peaks trace a single scattering surface, though the paper does not claim this.","Combining sector-by-sector SPF extraction with radiative transfer models could separate genuine dust-asymmetry signals from the geometric limb-brightening effect, which the paper motivates but does not pursue."],"forward_implications":["Large archival samples of ringed disks observed with high-contrast polarimeters can be processed quickly, giving SPFs for many disks without per-disk radiative transfer fits.","SPFs corrected for limb brightening will be systematically shallower than uncorrected ones, so dust grain sizes and compositions inferred from uncorrected curves will need revision.","Because the tool extracts SPFs globally or per sector, brightness asymmetries of the kind highlighted for HD 163296 and RX J1615 can be identified automatically and either excluded or analyzed separately.","The method allows each ring in a multi-ring disk to be treated separately, so radial trends in dust properties can be probed without assuming one SPF for the whole disk."],"supporting_citations":[{"why":"Supplies the previous SPF extraction approach and the extraction-zone formulation that DRAGyS adapts, and serves as the reference tool for comparison on real data.","marker":"Stolker et al. (2016)"},{"why":"Supplies the limb-brightening formula used in Eq. (3), as well as the multiple-scattering caveat for polarized SPFs.","marker":"Tazaki et al. (2023)"},{"why":"Supplies the archival polarized-intensity images, reference geometries, and reference SPFs used for the six-disk validation.","marker":"Ginski et al. (2023)"},{"why":"Supplies the MCFOST radiative transfer simulations that generate the synthetic disk images with known input geometry and intrinsic SPF.","marker":"Pinte et al. (2006)"},{"why":"Supplies the updated MCFOST radiative transfer code used to compute the synthetic images and theoretical SPFs for benchmarking.","marker":"Pinte et al. (2009)"},{"why":"Supplies the least-squares ellipse fitting method used to estimate the geometric parameters of the disk rings.","marker":"Halır & Flusser (1998)"},{"why":"Supplies the azimuthal Stokes $Q_\\phi$ construction used to produce the polarized-intensity images from which SPFs are extracted.","marker":"Schmid et al. (2006)"}],"fun_headline_variants":["Geometric tool extracts scattering phase from disk images","Limb-brightening fix sharpens disk dust scattering reads","DRAGyS fits ring ellipses to pull scattering phase function","Pure-geometry pipeline pulls scattering phase from disks","Fast disk geometry tool corrects limb-brightening bias"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method depends on the assumption that each disk's ring is a perfect circle lying on a single smooth, tilted surface whose height above the midplane rises as a fixed power of radius; if the ring is off-circle, the surface height follows a different law, or the bright arc does not mark the effective scattering surface, the geometry that feeds the correction is biased.","fun_headline_variants_meta":{"raw":{"variants":["Geometric tool extracts scattering phase from disk images","Limb-brightening fix sharpens disk dust scattering reads","DRAGyS fits ring ellipses to pull scattering phase function","Pure-geometry pipeline pulls scattering phase from disks","Fast disk geometry tool corrects limb-brightening bias"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000457,"raw_usage":{"total_tokens":2374,"prompt_tokens":1104,"completion_tokens":1270,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":720,"completion_tokens_details":{"reasoning_tokens":1190}},"tokens_in":720,"tokens_out":1270,"duration_ms":10144,"temperature":1.0,"reasoning_tokens":1190,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:00:21.546189+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run DRAGyS on synthetic images with known intrinsic SPFs while varying the flaring exponent (for example $\\chi$ = 0.5 or 1.5) or using an eccentric ring; if the limb-brightening-corrected SPF departs from the intrinsic one by more than the few-percent recovery claimed for circular, $\\chi\\approx1$ disks, the assumed single circular scattering surface is falsified. The same test applies to low-dust-mass models, where the brightness peak no longer traces the effective scattering surface, a case the paper itself shows in its Appendix D.","supporting_citations":[{"cited_title":"2016, Astronomy & Astrophysics, 596, A70","cited_arxiv_id":null,"evidence_quote":"Supplies the previous SPF extraction approach and the extraction-zone formulation that DRAGyS adapts, and serves as the reference tool for comparison on real data."},{"cited_title":"2023, The Astrophysical Journal, 953, 92","cited_arxiv_id":null,"evidence_quote":"Supplies the archival polarized-intensity images, reference geometries, and reference SPFs used for the six-disk validation."},{"cited_title":"2006, Astronomy & Astrophysics, 459, 797","cited_arxiv_id":null,"evidence_quote":"Supplies the MCFOST radiative transfer simulations that generate the synthetic disk images with known input geometry and intrinsic SPF."},{"cited_title":"J., Min, M., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the updated MCFOST radiative transfer code used to compute the synthetic images and theoretical SPFs for benchmarking."},{"cited_title":"M., Joos, F., & Tschan, D","cited_arxiv_id":null,"evidence_quote":"Supplies the azimuthal Stokes $Q_\\phi$ construction used to produce the polarized-intensity images from which SPFs are extracted."}],"review_version":1}