{"id":"d1d8efaf-034f-43b4-923c-ea6698e9d71d","arxiv_id":"2411.17057","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Optical polarimetric imaging of the post-AGB binary IRAS 08544-4431 reveals a bright, grey-coloured circumbinary disc with forward scattering, consistent with submicron porous dust aggregates.","lead":"Researchers obtained new optical polarimetric images of the dusty disc around the evolved binary star IRAS 08544-4431 using the SPHERE instrument on the VLT. The images show the disc is bright in polarized light, grey in colour from optical to near-infrared, and likely made of small porous dust clumps, similar to discs around young stars.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Grey-colour and porous-aggregate conclusions depend on band-to-band geometry that the paper itself shows is not constant; the V/I' PSF-smearing correction even used the H-band orientation, so the polarimetric colour may be partly a geometric artifact.","rationale":"The paper is a careful pilot study: the PDI reduction, frame selection, telescope-polarization correction, unresolved-central-polarization subtraction, PSF-smearing correction, deconvolution, and SNR masking are described with enough detail that the resolved detection and ~1.5% polarized brightness are plausible. I do not challenge the observational reality of the disc or the qualitative complexity of its morphology. The central inference that the surface dust consists of large porous aggregates of submicron monomers is, however, a chain: grey polarimetric colour + forward scattering + high polarization efficiency leads to porous aggregates. The weakest link is the colour measurement. Equation (3) is only a dust-scattering diagnostic if the scattering geometry is the same at all wavelengths; Section 4.4.1 states this explicitly, and Table 2 shows it is not true. Moreover, the PSF-smearing correction applied to the V/I' brightness was computed using the H-band orientation rather than the V/I' orientations, so part of the wavelength dependence in Q/I may be introduced by the correction itself. The quoted errors already make the grey classification weak on statistical grounds: eta_V/I = 0.24 ± 0.72 extends well outside the grey window. A rerun with orientation-matched corrections and matched masks is the decisive check. This supports the existing CONDITIONAL verdict rather than a rejection, because the data are new and the reduction is transparent; the condition should explicitly include the orientation-matched colour recomputation.","tokens_in":22102,"tokens_out":6421,"duration_ms":61408,"concrete_test":"Recompute the PSF-smearing correction maps of §3.2.5 for V and I' using the orientations from Table 2 (i = 55°, PA = 102°; i = 40°, PA = 102°) instead of the H-band orientation (i = 23°, PA = 121°), re-derive Q_phi/I_tot in V and I', and propagate the new values through Eq. (3). If either eta_V/I or eta_V/H moves by more than its quoted 1σ uncertainty (0.72 or 0.34), the grey-colour and porous-aggregate conclusions are not robust. In the same rerun, restrict the colour estimate to the intersection of the SNR≥5 masks deprojected with each band's own inclination; a significant shift in eta would confirm that the 'similar scattering geometry' assumption is the source.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The grey polarimetric colour (Eq. 3), on which the porous-aggregate dust conclusion rests, is measured under the explicit assumption in §4.4.1 of 'similar scattering geometries' across bands. That assumption is inconsistent with the paper's own Table 2: the fitted ring inclination changes from 55° (V) to 40° (I') to 23° (H). If different bands sample different surface elements, or the same elements at different scattering angles, then eta_V/I = 0.24 ± 0.72 and eta_V/H = -0.24 ± 0.34 mix dust scattering with geometric selection. In addition, the PSF-smearing correction maps in §3.2.5 used for V and I' were computed with the H-band orientation (i = 23°, PA = 121°), not with the V/I' orientations from Table 2; the authors checked consistency with IR-interferometric orientation, but not with their own optical fits. Since the correction factor is large close to the binary (Fig. 4), an orientation mismatch could shift V and I' polarized fluxes by more than the quoted uncertainties and change both colour gradients. The forward-scattering asymmetry in §4.3.2 is interpreted from a one-sided brightness rise without a quantitative phase-function fit, so it provides only weak independent support for the grain-size conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents V- and I'-band ZIMPOL polarimetric imaging of the post-AGB binary IRAS 08544-4431, combined with previously published H-band IRDIS data, to characterize the circumbinary disc in scattered polarized light. The authors describe a detailed reduction: polarimetric differential imaging, frame selection, telescope polarization correction, subtraction of unresolved central polarization, PSF-smearing correction, SNR masking, and deconvolution. They measure polarized disc brightness of about 1.3-1.5% of the stellar intensity, fit an elliptical ring whose fitted inclination changes from 55 deg (V) to 40 deg (I') to 23 deg (H), find a forward-scattering asymmetry, and compute polarimetric colour gradients consistent with a grey disc. They conclude that the surface dust is dominated by large porous aggregates of submicron monomers and emphasize the similarity of this post-AGB disc to protoplanetary discs.","tokens_in":22388,"tokens_out":7992,"duration_ms":76061,"significance":"If the conclusions hold, this is the first optical multi-wavelength polarimetric study of a post-AGB circumbinary disc and provides quantitative constraints on dust scattering properties that can be compared directly with protoplanetary-disc studies. The data reduction is careful: the paper includes telescope-polarization correction, unresolved central polarization subtraction, PSF-smearing correction following Ma et al. (2024), SNR-based masking, and Richardson-Lucy deconvolution with a reference PSF. The central quantities (polarized brightness ratios, colour gradients, ellipse parameters) are falsifiable and presented with uncertainties. The comparison with external theoretical models (Tazaki & Dominik 2022) and with the compilation of Ma et al. (2024) is appropriate. However, as detailed below, the key 'grey colour and porous aggregates' conclusion currently depends on an assumption about constant scattering geometry that is inconsistent with the paper's own orientation measurements, and on PSF-smearing corrections whose orientation dependence has not been quantified. These issues are addressable with additional analysis, so the manuscript merits revision rather than rejection.","major_comments":[{"comment":"The grey polarimetric colour, which is the basis for the porous-aggregate conclusion, is not statistically secure and is derived under the explicit assumption of 'similar scattering geometries across these wavelengths'. The measured gradients are eta_V/I = 0.24 ± 0.72 and eta_V/H = -0.24 ± 0.34; with these 1-sigma uncertainties the data are also fully consistent with significantly blue or red colours, so the classification 'grey' is only marginal. More importantly, Table 2 shows that the fitted inclination changes from 55 deg (V) to 40 deg (I') to 23 deg (H), so the geometry is demonstrably not the same across bands. If the resolved ring traces different surface elements or scattering angles in different bands, the measured gradients mix true dust scattering with geometric selection. The V-band ellipse in particular (a = 22 mas) is only marginally larger than the 30-mas pre-deconvolution resolution, and no validation with injected disc models is provided for the inclination fits. The authors should quantify the systematic effect on the colour gradients of using the actual V/I'/H orientations (e.g., by re-projecting to a common inclination or restricting the comparison to a common radial range), or should soften the grey-colour and dust-composition claims to reflect this degeneracy.","section":"§4.4.1, Eq. (3), Table 2"},{"comment":"The PSF-smearing correction maps used to derive the V- and I'-band polarized brightness were computed with the H-band orientation (i = 23 deg, PA = 121 deg), and consistency was checked only against IR interferometric orientation, not against the authors' own optical ellipse fits (i = 55 deg in V, i = 40 deg in I'; Table 2). Since Fig. 4 shows that the correction factor is large close to the binary, where a significant part of the measured Q_phi signal originates, an orientation mismatch could shift the V and I' corrected fluxes by more than the quoted uncertainties and thereby change both eta_V/I and eta_V/H. The authors should quantify this sensitivity, e.g., by recomputing the correction maps with the V and I' orientations and reporting the resulting brightness ratios.","section":"§3.2.5"},{"comment":"The forward-scattering asymmetry is inferred from a qualitative one-sided rise in the azimuthal brightness profile, without fitting a scattering phase function or accounting for the varying scattering angle along the inclined ring. In its current form this provides only weak, non-quantitative support for the micron-size grain interpretation. The authors should either fit a simple phase function (or at minimum report the contrast ratio between the two sides with uncertainties) or explicitly state that the porous-aggregate conclusion does not rely on the asymmetry.","section":"§4.3.2 and §5.2"}],"minor_comments":[{"comment":"The text contains a typo, 'SPHERE/ZMPOL', which should read 'SPHERE/ZIMPOL'.","section":"§3.2.5"},{"comment":"The caption says the white circles represent the size of the resolution element '(see Section 5)', but the resolution is defined in Section 3.2.6, not Section 5.","section":"Fig. 3 caption"},{"comment":"The notation for the colour gradient can be made clearer: define exactly which intensity is used (Q_phi/I_tot) and specify that the wavelengths lambda are the central filter wavelengths; currently the relation between the labels V, I', H and the numerical values is only implicit.","section":"Eq. (3) and §4.4.1"},{"comment":"The definition of position angle as 'rising counterclockwise from the vertical axis (North)' is ambiguous; standard astronomical usage is 'east of north', which should be stated explicitly to avoid confusion with the usual mathematical convention.","section":"§4.2"},{"comment":"When comparing with protoplanetary discs, the phrase 'IRAS08544-4431 circumbinary disc among the brightest protoplanetary discs' should be phrased as 'among the brightest discs in terms of polarized intensity' to avoid implying this evolved system is itself a protoplanetary disc.","section":"§5.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents a substantial new dataset and the reduction is generally careful. The principal concern is that the grey-colour and porous-aggregate conclusion relies on an assumption of similar scattering geometry across bands that the paper's own measurements contradict, and on a PSF-smearing correction whose orientation sensitivity is not tested. These are fixable with additional analysis, so I recommend revision rather than rejection: the authors should either demonstrate robustness of the colour gradients to the orientation choices or substantially soften the dust-composition claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a real observational first — SPHERE/ZIMPOL V and I' polarimetric imaging of a post-AGB circumbinary disc, with the IRDIS H data re-reduced for PSF smearing. The data reduction is careful and transparent: PDI, frame selection, telescope polarisation correction, subtraction of unresolved central polarisation, PSF-smearing correction, SNR masking. That earns credit. The basic detection — resolved elliptical ring, high azimuthal polarisation up to ~1.5%, forward-scattering asymmetry — is solid and new for this target in the optical.\n\nThe grey-colour claim is the soft spot, and it is soft in exactly the way the stress test says. The paper needs 'similar scattering geometries across these wavelengths' for Eq. 3 to be a clean dust diagnostic, but Table 2 shows inclination changing from 55° in V to 40° in I' to 23° in H. If different bands are seeing different surface elements at different scattering angles, the colour gradient mixes geometry with dust properties. On top of that, the V/I' PSF-smearing correction maps were computed with the H-band orientation (i = 23°, PA = 121°), not with the V/I' orientations from Table 2. The authors checked against IR interferometric orientation but not against their own optical ellipse fits; near the binary where the correction factor is large, that could shift the V/I' fluxes enough to matter. The uncertainties on eta are large enough that 'grey' is not strongly constrained: eta_V/I = 0.24 ± 0.72, eta_V/H = -0.24 ± 0.34. The porous-aggregate conclusion is an interpretation consistent with the data, not a tight measurement.\n\nThe forward-scattering asymmetry is qualitatively presented without a phase-function fit, so it supports the grain-size argument but does not carry it alone. The warping discussion is speculation, properly flagged as such. The interstellar polarisation check is reasonable.\n\nNet: new data and a careful reduction that will be useful to the evolved-binary community. The colour and dust-composition conclusions need revision — either a quantitative model that handles wavelength-dependent geometry, or a re-analysis of the PSF correction with the optical orientations — before they can be taken as established. This deserves peer review; I would not desk-reject it.","headline":"A careful, genuinely new optical polarimetric dataset for a post-AGB circumbinary disc; the grey-colour and porous-aggregate conclusions rest on assumptions the paper's own morphology table undermines.","tokens_in":22957,"tokens_out":1831,"would_cite":true,"duration_ms":17658,"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 circumbinary disc of IRAS 08544-4431 is resolved in optical scattered light with polarized brightness up to 1.5% of the stellar intensity, a grey polarimetric colour across V, I', and H, and a forward-scattering asymmetry.","keywords":["post-AGB binary","circumbinary disc","polarimetric imaging","SPHERE/ZIMPOL","scattered light","dust aggregates","polarimetric colour","second-generation disc"],"falsifier":"A radiative-transfer model that adopts a single disc inclination and scattering surface across the V, I', and H bands and predicts the polarized intensities; if the predicted colour is non-grey, then the observed grey colour is a geometric artifact rather than an intrinsic dust property.","tokens_in":18,"feed_emoji":"🔭","tokens_out":3098,"duration_ms":91617,"temperature":0.7,"pith_summary":"This paper reports the first optical multi-wavelength polarimetric imaging of a post-AGB circumbinary disc, resolving the disc around IRAS 08544-4431 in the V and I' bands with SPHERE/ZIMPOL and combining the result with previously published H-band IRDIS data. The disc shows azimuthally polarized scattered light reaching about 1.5% of the stellar intensity, a grey polarimetric colour across the three bands, and a forward-scattering asymmetry that places the near side of the disc to the north. The authors argue that these observed properties match the predicted scattering behaviour of large porous aggregates of submicron-sized monomers, which are also invoked for protoplanetary discs. If correct, this strengthens the analogy between second-generation discs around evolved binaries and planet-forming discs around young stars.","feed_headline":"Dust disc around post-AGB binary resolved at three wavelengths","feed_subtitle":"Optical polarimetry reveals grey colour and forward scattering, tying evolved-star discs to planet-forming discs.","key_machinery":"The analysis is carried by azimuthally polarized intensity Q_phi, which isolates single scattering of stellar light on the disc surface and is used to measure disc brightness, morphology, and substructures. The photometric calibration relies on correction maps for PSF smearing and polarimetric cancellation applied to the V, I', and H bands, and the wavelength dependence is quantified with the logarithmic polarimetric colour index eta. The forward-scattering asymmetry and the grey colour are then compared with theoretical dust-scattering models for aggregates to infer the surface dust composition.","core_discovery":"The central claim is that the surface layers of the circumbinary disc around IRAS 08544-4431 scatter starlight with the same efficiency in the optical and near-infrared, producing a grey polarimetric colour (eta_V/I = 0.24 +/- 0.72 and eta_V/H = -0.24 +/- 0.34), while the resolved disc brightness is high at about 1% of the stellar intensity and is dominated by forward scattering. The disc is resolved in the V and I' bands as a bright elliptical ring with substructures, and the resolved surface appears more elliptical at shorter wavelengths, which the paper interprets as a possible sign of disc warping or shadowing. Combining these multi-wavelength polarimetric measurements, the paper concludes that the dust on the disc surface is consistent with large porous aggregates composed of submicron-sized monomers, matching the properties inferred for many protoplanetary discs.","pith_inferences":["A testable extension is to fit the V, I', and H images with a single disc geometry; if the grey colour persists under a common inclination and scattering surface, it is a true dust property rather than a geometric artifact.","The wavelength-dependent ring ellipticity could be caused by the scattering surface height increasing with wavelength, which would be testable with high-resolution ALMA continuum observations of the disc midplane.","The close similarity to TW Hya suggests that large porous aggregates may be common in both post-AGB circumbinary discs and young protoplanetary discs, implying a shared dust-growth pathway despite very different formation histories."],"forward_implications":["If the grey polarimetric colour is real, the disc surface dust is not dominated by small Rayleigh scatterers but by larger porous aggregates.","If the forward-scattering asymmetry is real, the northern part of the circumbinary disc is tilted toward the observer.","If the wavelength-dependent ring morphology is real, the disc may be warped or shadowed at shorter wavelengths.","The observed polarized brightness of about 1% of the stellar intensity places this post-AGB disc among the brightest protoplanetary discs, supporting the disc-evolution analogy."],"supporting_citations":[{"why":"Supplies the H-band IRDIS polarimetric data and the initial disc orientation and morphology that this paper extends to optical wavelengths.","marker":"Andrych et al. 2023"},{"why":"Provides the PSF-smearing correction method used to recover intrinsic polarized disc brightness in all three bands.","marker":"Ma et al. 2024"},{"why":"Provides the theoretical dust-scattering models for porous aggregates used to interpret the grey colour and forward scattering.","marker":"Tazaki & Dominik 2022"},{"why":"Supplies the protoplanetary disc polarized-brightness statistics and the deprojection methodology used for radial brightness profiles.","marker":"Avenhaus et al. 2018"},{"why":"Defines the polarized-contrast methodology used to place IRAS 08544-4431 among the brightest scattered-light discs.","marker":"Garufi et al. 2014"},{"why":"Documents the SPHERE/ZIMPOL instrument and its polarimetric calibration, establishing the accuracy of the optical measurements.","marker":"Schmid et al. 2018"}],"fun_headline_variants":["Grey disc around evolved star hints at planet-forming analog","Polarimetry reveals scattering uniformity in post-AGB disc","Evolved binary's disc mimics protoplanetary dust properties","Optical and NIR polarimetry find grey scattering disc"],"cache_read_input_tokens":25088,"weakest_assumption_plain":"The colour comparison assumes the disc scatters light from the same geometry in all three bands, but the paper's own fits give different inclinations for the resolved ring in each band, so a geometric difference could produce the apparent grey colour.","fun_headline_variants_meta":{"raw":{"variants":["Grey disc around evolved star hints at planet-forming analog","Polarimetry reveals scattering uniformity in post-AGB disc","Evolved binary's disc mimics protoplanetary dust properties","Optical and NIR polarimetry find grey scattering disc"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000197,"raw_usage":{"total_tokens":1423,"prompt_tokens":1062,"completion_tokens":361,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":678,"completion_tokens_details":{"reasoning_tokens":293}},"tokens_in":678,"tokens_out":361,"duration_ms":4446,"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-12T12:34:13.842824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A radiative-transfer model that adopts a single disc inclination and scattering surface across the V, I', and H bands and predicts the polarized intensities; if the predicted colour is non-grey, then the observed grey colour is a geometric artifact rather than an intrinsic dust property.","supporting_citations":[],"review_version":1}