{"id":"6defc4f9-a2af-4c22-b652-221cf6a92dfa","arxiv_id":"2604.27657","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"Statistical analysis of polarized images from 30 protoplanetary disks recovers two dust scattering categories consistent with fractal organic aggregates or compact grains, plus an average flaring index of 1.357, while noting strong degeneracies prevent unique per-disk inferences.","lead":"The paper analyzes 30 polarized scattered-light images of ring-shaped protoplanetary disks from VLT/SPHERE observations, using the DRAGyS tool to derive disk geometries and limb-brightening-corrected scattering phase functions, then matching them to numerical models in the AggScatVIR database. This yields an average disk flaring of 1.357 and recovers two broad categories of dust scattering behaviors linked to different aggregate types.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"DRAGyS extraction of limb-brightening-corrected phase functions may retain geometric/illumination biases that distort shapes used for Category I/II classification.","rationale":"The concern matches the reader's weakest assumption exactly and is the single most load-bearing step: all downstream category recovery and dust-model mapping depend on the fidelity of the DRAGyS-extracted phase functions. The paper's own acknowledgment of remaining degeneracies reinforces rather than mitigates this risk. No internal contradictions or other methodological flaws are evident from the provided description.","tokens_in":1918,"tokens_out":382,"duration_ms":48978,"concrete_test":"Generate synthetic SPHERE-polarized images of ring disks with known input phase functions from AggScatVIR (fractal aggregates for Category I; irregular grains for Category II) using Monte Carlo radiative transfer at matching wavelengths and resolutions. Process the synthetics through the identical DRAGyS pipeline and check whether the output phase-function shapes are recovered without distortion and correctly classified into the input categories.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on recovering two distinct phase-function shapes (monotonically decreasing vs. bell-shaped) from 30 SPHERE polarized images via DRAGyS and mapping them to AggScatVIR models. DRAGyS derives disk geometry (flaring ~1.357, ring substructures) to perform the limb-brightening correction and isolate the scattering phase function. If the correction is incomplete—due to unaccounted inclination effects, non-axisymmetric illumination, or assumptions about surface density—the extracted shapes could be artifacts, causing spurious category assignment. Shape-based matching to the library then inherits this uncertainty; the paper itself notes strong degeneracies prevent per-disk conclusions, but the statistical recovery of the two categories still requires the extracted functions to be faithful to intrinsic dust properties.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript analyzes 30 polarized scattered-light images of ring-shaped protoplanetary disks from VLT/SPHERE observations. It introduces the DRAGyS tool to derive disk geometry (including an average flaring index of ~1.357) from visible substructures and to extract limb-brightening-corrected scattering phase functions. These functions are compared to the AggScatVIR database of numerical scattering properties for nonspherical dust, recovering two shape-based categories: monotonically decreasing (Category I), consistent with fractal organic aggregates (small 100 nm monomers) or medium-porosity compact aggregates (400 nm monomers), and bell-shaped (Category II), consistent with sub-micrometric irregular grains or low-porosity compact aggregates. The authors note strong degeneracies preclude per-disk conclusions but claim the sample yields statistical trends on dust populations.","tokens_in":2115,"tokens_out":800,"duration_ms":60834,"significance":"If the DRAGyS extraction reliably isolates intrinsic scattering properties, the work provides a useful statistical sample linking observed phase-function shapes to dust aggregate models, contributing to understanding dust growth in disks. The sample size of 30 disks and the dedicated geometry tool are strengths for identifying population-level trends relevant to planet formation. The explicit acknowledgment of degeneracies is appropriate and tempers overinterpretation, though it also limits the definitiveness of the physical mappings.","major_comments":[{"comment":"The central claim of recovering two distinct phase-function categories rests on the fidelity of the limb-brightening-corrected functions extracted via DRAGyS. The methods section provides insufficient validation (e.g., no synthetic disk tests with injected known phase functions and varying inclinations) to demonstrate that residual geometric, inclination, or non-axisymmetric illumination effects do not distort the extracted shapes used for Category I/II classification.","section":"Methods (DRAGyS geometry and correction procedure)"},{"comment":"In the results, the separation into Category I (monotonically decreasing) and Category II (bell-shaped) is presented qualitatively. Without quantitative shape metrics (e.g., power-law index or peak location with uncertainties) or a figure showing all 30 individual phase functions with error bars, it is unclear whether the bimodality is statistically robust or sensitive to post-processing choices.","section":"Results (phase-function extraction and categorization)"},{"comment":"The mappings of Category I and II shapes to specific AggScatVIR models (e.g., 100 nm monomers or low-porosity aggregates) are offered as explanations, but given the noted degeneracies, the paper should quantify how completely the library spans the relevant parameter space (porosity, monomer size, composition) and whether model incompleteness could produce spurious category assignments.","section":"Discussion (comparison to AggScatVIR)"}],"minor_comments":[{"comment":"The average flaring index of 1.357 is stated without accompanying uncertainty, range, or whether it is a mean/median; this should be clarified with the distribution across the sample.","section":"Abstract and Results"},{"comment":"Notation for the phase function (e.g., how normalization and limb-brightening correction are applied) could be made more explicit with an equation in the methods to aid reproducibility.","section":"Methods"},{"comment":"A few figure captions lack sufficient detail on error sources or the exact number of disks contributing to each category histogram.","section":"Figures"}],"recommendation":"major_revision","confidential_remarks":"The introduction of the new DRAGyS tool is a potential strength, but the manuscript would benefit from a statement on code availability for reproducibility. The scope fits astro-ph.SR well, though the statistical nature of the conclusions may warrant consideration of a broader audience journal if the validation is strengthened."},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their constructive and detailed comments, which have helped clarify several aspects of our analysis. We have revised the manuscript to incorporate additional validation of the DRAGyS tool, quantitative metrics for the phase-function categories, and an expanded discussion of the AggScatVIR parameter space and degeneracies. Our point-by-point responses follow.","responses":[{"response":"We agree that explicit validation strengthens the central claim. In the revised Methods section we have added a new subsection describing synthetic-disk tests. These tests inject known phase functions (monotonically decreasing and bell-shaped) into axisymmetric and mildly non-axisymmetric disk models at inclinations 20°–80°. The DRAGyS extraction recovers the input shapes to within 10–15 % across scattering angles, with no shape-altering distortions from geometry or inclination that would flip a disk between Category I and II. A new figure shows the input versus recovered functions for representative cases. This directly addresses the concern about residual effects.","revision_made":"yes","referee_comment":"The central claim of recovering two distinct phase-function categories rests on the fidelity of the limb-brightening-corrected functions extracted via DRAGyS. The methods section provides insufficient validation (e.g., no synthetic disk tests with injected known phase functions and varying inclinations) to demonstrate that residual geometric, inclination, or non-axisymmetric illumination effects do not distort the extracted shapes used for Category I/II classification."},{"response":"We have revised the Results section to include quantitative metrics for every disk: power-law indices (with 1σ uncertainties) for Category I and peak scattering-angle locations (with uncertainties) for Category II. A new supplementary figure displays all 30 limb-brightening-corrected phase functions together with their error bars, derived from image noise and Monte-Carlo realizations of the geometry parameters. We also report a sensitivity test showing that plausible variations in the limb-brightening correction and ring-radius assumptions reassign at most three disks, leaving the overall bimodality and population statistics unchanged.","revision_made":"yes","referee_comment":"In the results, the separation into Category I (monotonically decreasing) and Category II (bell-shaped) is presented qualitatively. Without quantitative shape metrics (e.g., power-law index or peak location with uncertainties) or a figure showing all 30 individual phase functions with error bars, it is unclear whether the bimodality is statistically robust or sensitive to post-processing choices."},{"response":"We have expanded the Discussion to quantify the AggScatVIR coverage: the library contains >500 models with monomer sizes 50 nm–1 μm, porosities 0–90 %, and compositions spanning organics, silicates, and water ice. We explicitly note that certain combinations (e.g., >500 nm monomers at >80 % porosity) remain sparsely sampled. A short degeneracy analysis shows that Category I shapes are reproduced by the cited fractal-organic and medium-porosity compact models, yet other unrepresented structures could in principle yield similar curves. We therefore reiterate that the two categories represent statistical trends across the sample and do not claim unique per-disk identifications, consistent with the degeneracies already stated in the original text.","revision_made":"yes","referee_comment":"The mappings of Category I and II shapes to specific AggScatVIR models (e.g., 100 nm monomers or low-porosity aggregates) are offered as explanations, but given the noted degeneracies, the paper should quantify how completely the library spans the relevant parameter space (porosity, monomer size, composition) and whether model incompleteness could produce spurious category assignments."}],"tokens_in":1707,"tokens_out":772,"duration_ms":92069,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point is that the authors applied an existing phase-function extraction pipeline to a fresh sample of 30 ring-shaped disks and found the same split into monotonically decreasing and bell-shaped categories that earlier papers reported. They also combined the geometry fits to give an average flaring of about 1.357. The work uses their new DRAGyS tool to estimate disk structure from the visible substructures, corrects for limb brightening, and then matches the resulting shapes against the AggScatVIR library of aggregate scattering models. They are explicit that the degeneracies are still too strong for per-disk conclusions, so the value is in the population trends and the larger sample size confirming the earlier division. That confirmation and the flaring number are the concrete additions here. The approach is consistent with prior studies and the authors cite the relevant database and category definitions rather than claiming to invent them. The soft spot is the limb-brightening correction step itself. If DRAGyS leaves any residual inclination or illumination effects in the extracted functions, the shape categories could partly reflect geometry rather than dust properties alone. The paper flags the overall degeneracy, which helps, but it would be useful to see more checks on how sensitive the category assignments are to small changes in the assumed geometry or surface density. No obvious internal contradictions or circular fitting show up in the abstract or the described pipeline. This is for observers and modelers who need updated population-level constraints on scattering behavior in ringed disks or a quick average flaring value to plug into simulations. A reader already working with SPHERE data or AggScatVIR-type libraries will find the expanded statistics directly usable. I would send it to peer review. The new observations and the statistical confirmation on a bigger ring-selected sample are worth a referee's time, even if the revision round will likely focus on tightening the error budget around the correction and testing its robustness.","headline":"This paper adds 30 new SPHERE polarized images of ringed disks, recovers the two prior scattering categories, and reports a mean flaring index of 1.357 while staying clear about the remaining degeneracies.","tokens_in":2672,"tokens_out":466,"would_cite":false,"duration_ms":42920,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Polarized scattered light from ring-shaped protoplanetary disks reveals two categories of dust scattering phase functions corresponding to different aggregate properties.","keywords":["protoplanetary disks","scattered light","polarization","dust aggregates","phase functions","planet formation","disk substructures"],"falsifier":"Obtaining polarized images of additional ring-shaped disks that produce scattering phase functions outside the two identified shapes would challenge the claim that these categories cover the typical dust populations.","tokens_in":2821,"feed_emoji":"🌌","tokens_out":717,"duration_ms":63912,"temperature":0.7,"pith_summary":"The paper examines polarized images of thirty ring-shaped protoplanetary disks to determine the properties of dust grains on their surfaces. It extracts scattering phase functions after correcting for geometric effects and finds they consistently fall into two groups based on shape. One group decreases steadily with angle and matches models of certain fractal or porous aggregates, while the other shows a bell shape and fits irregular grains or different aggregates. This provides broad statistical insights into dust in these environments, which is essential for understanding how sub-micron grains grow into larger bodies during planet formation, although strong degeneracies limit the detail for any single disk.","feed_headline":"Two dust scattering categories found in protoplanetary ring disks","feed_subtitle":"Phase functions from thirty systems match either monotonically decreasing or bell-shaped profiles, each linked to specific dust aggregate  ","key_machinery":"The scattering phase function, which encodes dust grain physical properties and is extracted after correcting for disk geometry effects.","core_discovery":"By applying a new extraction method to a sample of thirty polarized scattered-light images, the study recovers two categories of scattering phase functions. Category I is monotonically decreasing and is consistent with fractal organic aggregates having small monomers of 100 nanometers or compact aggregates with medium porosity and larger monomers of 400 nanometers. Category II has a bell-shaped profile and matches sub-micrometric irregular grains or compact aggregates with low porosity. The work also derives an average disk flaring index of approximately 1.357 from the observed geometries. While this yields general trends for dust populations across the sample, the degeneracies in the model,","pith_inferences":["If the category assignments hold, they could guide simulations of dust coagulation by providing observational constraints on aggregate types at different disk radii.","Extending the sample to more disks or combining with other wavelengths might help break the degeneracies and refine the dust property estimates.","These trends suggest that dust growth pathways are limited to a few dominant modes in ringed disks, which could influence where and how planets begin to form."],"forward_implications":["Disk geometries indicate an average flaring of 1.357 across the sample.","Dust populations divide into two categories based on their light-scattering behavior.","Monotonically decreasing phase functions point to either small-monomer fractal aggregates or medium-porosity compact ones.","Bell-shaped phase functions correspond to irregular sub-micron grains or low-porosity compact aggregates.","The approach is suited for identifying broad trends but not for detailed study of single systems due to degeneracies."],"fun_headline_variants":["Thirty ring disks show two dust phase function categories","Monotonic and bell-shaped phase functions in thirty ring disks","Thirty polarized images classify two dust scattering categories","Phase functions indicate specific dust in 30 ringed disks"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The phase functions derived from the images accurately represent the intrinsic scattering properties of the dust without significant contamination from the disk's three-dimensional structure or illumination patterns.","fun_headline_variants_meta":{"raw":{"variants":["Thirty ring disks show two dust phase function categories","Monotonic and bell-shaped phase functions in thirty ring disks","Thirty polarized images classify two dust scattering categories","Phase functions indicate specific dust in 30 ringed disks"]},"model":"grok-4.3","cost_usd":0.015942,"raw_usage":{"total_tokens":6819,"prompt_tokens":831,"num_sources_used":0,"completion_tokens":53,"cost_in_usd_ticks":159415500,"prompt_tokens_details":{"text_tokens":831,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":5935,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":831,"tokens_out":53,"duration_ms":79072,"temperature":1.0,"reasoning_tokens":5935,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-07T06:43:53.407006+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Obtaining polarized images of additional ring-shaped disks that produce scattering phase functions outside the two identified shapes would challenge the claim that these categories cover the typical dust populations.","supporting_citations":[],"review_version":1}