{"id":"a5b76d33-808c-4a02-b100-6cf647d45ca7","arxiv_id":"2502.02376","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Computer models of the HD98800 disc predict an 8 to 11 year optical eclipse beginning around now, with duration and shape sensitive to dust mass and viscosity.","lead":"This paper uses hydrodynamical and radiative transfer simulations to predict how light from the HD98800 quadruple star system will fade as a dusty disc passes in front of one pair of stars. It matters because the upcoming eclipse is a rare chance to measure the disc's mass, size, and viscosity directly.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Predicted ingress dates are set by the assumed 6 AU outer edge, not by viscous evolution: t_nu ~ 10^5 yr >> 1150 yr runtime, so the no-dimming fallback to low alpha is untested and the central timing claim is an initial-condition artefact.","rationale":"The reader's CONDITIONAL verdict already captures the core problem: the quantitative timing prediction is in tension with the authors' own monitoring data. My stress-test sharpens this by identifying why that tension is not resolved by the paper's proposed low-viscosity explanation. The simulations run only 1150 years before the transit, whereas the viscous timescale at 6 AU with the fiducial alpha is about 10^5 years, so the outer edge is essentially fixed by the initial condition at 6.0 AU. Consequently, the start dates in Table 3 are not a robust outcome of viscous spreading; they follow from a chosen initial radius. The paper acknowledges the SPH viscosity floor near alpha~0.01, yet uses alpha=0.005 as fiducial and does not provide any simulation or analytic model for the lower-alpha regime invoked to explain the absence of dimming. This means the central abstract claim is supported only by an extrapolation. The comparative light-curve results and parameter-dependence trends may still be useful, which is why the reader's CONDITIONAL stance is appropriate. No new verdict category is needed, but the condition should explicitly require a sensitivity test on the initial outer edge and a quantitative comparison with the photometric non-detection limits.","tokens_in":16757,"tokens_out":10755,"duration_ms":109310,"concrete_test":"Re-run the fiducial model with the same parameters but initialise the outer disc edge at 4.6 AU instead of 6.0 AU, and recompute the ingress date and 5% decrement epoch in Table 3. If the ingress shifts to after early 2025 or the 5% dip disappears, the central timing claim is an artifact of the assumed initial outer edge and must be reframed as a sensitivity test rather than a prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All seven models in Table 3 place the 5% transit ingress between 03-2023 and 06-2024, yet Sec 3.6 reports no observed dimming as of early 2025. The paper explains this by saying the disc has a smaller radial extent due to lower viscosity, but that explanation is not supported by the simulations. The outer edge is initialised at 6.0 AU (Sec 2.1), and with H/R=0.05 and alpha=0.005, the viscous timescale at 6 AU is t_nu ~ [alpha (H/R)^2 Omega]^-1 ~ 10^5 yr, far longer than the 1150-yr pre-transit integration. The outer-edge radius at transit is therefore effectively an initial condition, not a converged result of viscous spreading, and the alpha=0.005 vs 0.01 comparisons in Sec 3.3 cannot establish physical spreading over the simulated time. The observed non-detection is evidence against the assumed 6 AU optical edge, but the low-alpha fallback is an extrapolation outside the modelled regime (SPH viscosity floor ~0.01, Sec 3.3) and outside the simulated timescale. The central claim that the transit should already have begun is thus load-bearing on an unvalidated initial condition; it also conflates gas extent with the effective optical edge, since CO gas at 6.4 AU can coexist with a smaller dust edge.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents SPH simulations of the HD 98800 circumbinary disc and synthetic optical light curves for the predicted eclipse of the A binary by the disc around B. Seven models vary dust mass, gas mass, and alpha viscosity. The authors find that higher dust mass or alpha lengthen the transit, that gas mass has a smaller effect, and that spiral arms excited by the wide binary introduce asymmetries. They predict the transit began around 2023-2024, should last 8-11 years, and interpret the lack of observed dimming as evidence for a smaller disc due to lower viscosity than SPH can simulate.","tokens_in":16997,"tokens_out":4564,"duration_ms":50053,"significance":"If the predictions hold, this is a rare, testable forecast of a decade-long occultation that could constrain disc mass, viscosity, and structure. The systematic parameter study and the forward-modeling approach are strengths, as are the SED-calibrated dust masses and the explicit discussion of observing cadence and degeneracies. However, the central timing claim rests on assumptions about the initial outer disc edge and on treating dust as a passive gas tracer; the no-dimming interpretation is therefore less secure than the abstract suggests. The paper is valuable as a parameter-dependence study, but the headline prediction needs reframing or additional support.","major_comments":[{"comment":"The predicted ingress dates are set by the assumed 6.0 AU outer edge rather than by viscous evolution. With H/R ~ 0.05 and alpha = 0.005, the viscous timescale at 6 AU is ~1e5 yr, two orders of magnitude longer than the 1150-yr pre-transit integration, so the outer-edge radius during transit is effectively an initial condition. The no-dimming observation therefore constrains the assumed optical edge, not the viscosity, and the inference in §3.6 and the abstract that lower viscosity explains the non-detection is an extrapolation outside both the simulated timescale and, as §3.3 states, the resolved viscosity regime (SPH floor ~0.01). Please present the ingress dates as an initial-condition sensitivity, run models with smaller initial outer edges, or provide an analytic spreading estimate that justifies the adopted edge.","section":"§2.1, Table 3"},{"comment":"The dust is treated as a passive tracer of the gas with a fixed gas-to-dust ratio and an MRN size distribution, so the optical depth profile is proportional to the gas surface density by construction. Because dust evolution is neglected, the observed non-detection can equally be explained by a dust edge smaller than the gas edge (the introduction already notes 4.6 AU dust versus 6.4 AU CO extents) or by a radially varying dust-to-gas ratio. The conclusion that the disc is smaller 'due to lower viscosity' conflates gas extent with the effective optical edge; this degeneracy should be quantified or the claim softened.","section":"§2.3, Eq. (1), §3.7"},{"comment":"The statement that ongoing photometry 'clearly rules out' the higher dust mass and higher viscosity models is not supported by any quantitative upper limit. No cadence, photometric precision, or magnitude of the non-detection is given, and the 5% threshold is an assumption; a gradual ingress below this threshold could be missed. Please provide the actual monitoring limit (or a reference) and propagate relevant uncertainties into the start and end dates in Table 3.","section":"§3.6"}],"minor_comments":[{"comment":"There are several typos, including 'Wepresent' in the abstract, 'To ensure' in §2.1, and 'effect' should be 'affect' in §1 ('may effect the evolution').","section":"Abstract, §1, §2.1"},{"comment":"The sentence 'observations ... have also provided estimates of alpha that suggest a that alpha can vary' contains a stray word and should be rewritten.","section":"§2.2.2"},{"comment":"The sentence 'there is no publication to cite yet' should be replaced with a description of the monitoring program (telescope, cadence, passband, precision) or a citation once available, so readers can assess the non-detection claim.","section":"§3.6"},{"comment":"The caption states that stars are 'plotted on top' although they pass behind the disc; this is visually confusing and should be clarified.","section":"Figure 6 caption"},{"comment":"The statement that the lowest viscosity that can be modelled is ~0.01 should be introduced earlier and explicitly flag that the fiducial alpha = 0.005 model is therefore at the numerical floor, which affects the interpretation of Figures 10 and the fiducial light curve.","section":"§3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper's central prediction is time-sensitive and falsifiable, which is a strength, but the absence of a quantitative monitoring limit is a key gap. The interpretation of no dimming as low viscosity is only one of several possibilities given the dust-tracing assumption; I would encourage the editor to require the authors to either provide the monitoring limit or substantially soften the central claim. The parameter-dependence results are solid enough to warrant revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a useful paper: the first synthetic optical light curves for the HD98800 circumbinary disc transit, and a clean parameter study showing that dust mass and alpha viscosity lengthen the transit, gas mass is minor, and spiral arms introduce asymmetry. The authors are honest about caveats and the modelling is sensible, using public SPH and radiative transfer codes.\n\nThe soft spot is the absolute timing. The predicted ingress dates all come from assuming the disc's outer edge starts at 6 AU (matching CO gas), but the optical edge is set by small dust, which may well be smaller. Over the 1150-year pre-transit run, viscous spreading is negligible — t_nu ~ 1e5 yr at 6 AU for alpha=0.005 — so the outer edge is essentially an initial condition. The alpha=0.005 vs 0.01 comparison therefore cannot demonstrate physical spreading; it mainly shows numerical floor effects. The no-dimming observation is real evidence against a 6 AU optical edge, but the paper's fallback that a lower-than-SPH viscosity explains it is not supported by the simulations. The dust-traces-gas assumption is acknowledged but could change the edge location substantially.\n\nThese are addressable. The light curves and qualitative dependencies are likely robust; the absolute dates should be framed as conditional on the initial disc size, with no error bars. I'd send it to a good referee.","headline":"First synthetic light curves for the HD98800 transit, with a clean parameter study, but the absolute timing predictions are hostage to the assumed 6 AU outer edge and the no-dimming fallback doesn't hold up.","tokens_in":17620,"tokens_out":3307,"would_cite":false,"duration_ms":33791,"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 HD 98800 A binary should be entering a decade-long transit behind the circumbinary disc of the B binary right now, and the paper shows how the light curve will reveal the disc's dust, viscosity, and spiral structure.","keywords":["circumbinary discs","transit light curves","HD 98800","alpha viscosity","dust mass","spiral arms","smoothed particle hydrodynamics","radiative transfer"],"falsifier":"A photometric monitoring campaign at 0.5–0.8 μm that finds no ≥5% decrease in the A-binary flux by the end of 2025 would contradict the models' claim that ingress should be underway, pointing instead to a disc with a smaller radial extent or lower gas-to-dust ratio than simulated; a transit that begins but lasts much longer than 11 years would falsify the assumed outer-edge truncation.","tokens_in":16464,"feed_emoji":"🔭","tokens_out":8815,"duration_ms":80592,"temperature":0.7,"pith_summary":"The paper predicts that the A binary in the HD 98800 quadruple system should currently be entering a transit behind the circumbinary disc around the B binary, an occultation expected to last 8–11 years, though no dimming has yet been observed. Using hydrodynamical simulations that vary dust mass, gas mass, and alpha-viscosity, the authors generate synthetic optical light curves and show that the duration, depth, and asymmetry of the transit map onto disc properties: higher dust mass or higher viscosity lengthens the transit, while gas mass mainly affects the spiral-arm phase. The models also show that the outer binary's periastron passage excites spiral arms that make the second half of the transit longer and asymmetric, providing a direct diagnostic of disc–binary interaction. Because photometric monitoring shows no dimming as of early 2025, the paper concludes that the real disc is probably smaller in radial extent than simulated, likely because the true viscosity is below the floor accessible to smoothed-particle hydrodynamics.","feed_headline":"HD 98800's disc should be dimming its partner star now","feed_subtitle":"Simulations predict an 8–11 year transit that maps the disc's dust, viscosity, and spiral arms.","key_machinery":"The machinery is a suite of three-dimensional smoothed-particle hydrodynamics (SPH) simulations of the quadruple system, run for 1150 years so the disc settles and is tidally truncated, with an additional 25 years covering the transit. Dust is not simulated directly: the gas surface density is scaled by a fixed gas-to-dust ratio and a grain size distribution n(a) ∝ a^−3.5 to obtain the dust opacity, and the stellar fluxes are attenuated by exp(−tau). Light curves are generated by following the A-binary stars along their orbits behind the disc, adding the constant flux of BaBb as a baseline. The parameter space spans dust masses 0.033–3.3 Earth masses, gas masses 3.3–330 Earth masses, and $\\alpha$ values 0.005, 0.01, and 0.05, chosen to bracket the observationally allowed ranges.","core_discovery":"The central claim is a specific, testable prediction: the AaAb binary will pass behind the disc around BaBb, producing a transit that begins around late 2023 to mid-2024 (depending on model) and ends between 2032 and late 2034, with a detectable 5% flux drop defining ingress and egress. The paper establishes a parameter mapping: increasing dust mass or alpha-viscosity widens the transit and shifts ingress earlier, whereas gas mass has a minor effect on timing but makes the spiral-arm region slightly more opaque. A key discovery is that the periastron passage of the outer binary excites spiral arms at the disc's outer edge, which prolongs the second half of the transit and breaks the light curve's symmetry, so asymmetry is a direct diagnostic of outer-binary interaction. The observed absence of dimming as of early 2025 rules out the high-dust and high-viscosity models and points to a disc whose outer radius is smaller than the simulated 6 AU initialisation, consistent with a viscosity below the SPH-simulable floor.","pith_inferences":["The delay itself is a measurement: if the transit starts later than the model range, the lag can be converted into an upper bound on the disc's effective viscosity and a lower bound on its outer radius, independent of assumptions about the dust-to-gas ratio.","Spiral-arm asymmetry could distinguish between viscosity regimes: in a low-viscosity disc the arms persist and the egress asymmetry is sharp, whereas in a high-viscosity disc they are smeared out, offering a probe of turbulent diffusivity.","The same occultation geometry should occur in other hierarchical systems with polar or misaligned discs; predicting transit timing and asymmetry for those systems would turn HD 98800 into a general method rather than a one-off event.","If no dimming is seen by the end of 2025, one should check whether photoevaporation or planet formation has carved the outer gas disc below the ALMA-resolved extent, since the paper's low-viscosity explanation is not the only way to make a compact disc."],"forward_implications":["If the transit begins, the light curve will map the disc's optical depth as a function of radius, giving a direct measure of the dust distribution across the disc.","The transit duration and ingress/egress timing will constrain the disc's radial extent and dust mass, helping break degeneracies with assumed gas-to-dust ratios.","Asymmetry between the first and second halves of the transit will reveal spiral arms excited by the outer binary, and the extra flux at egress indicates those arms are optically thin.","The absence of dimming so far suggests a more compact disc, implying a viscosity below about 0.01, and a delayed ingress would provide a lower bound on the disc's radial extent.","Two-colour observations at 0.5 and 0.8 microns during ingress and egress can separate dust reddening from the stellar colour baseline, potentially probing the grain size distribution."],"supporting_citations":[{"why":"Supplies the disc parameters (dust and gas mass lower bounds, radial extent 2.5–4.6 AU dust and 1.6–6.4 AU gas) and the polar-aligned interpretation that the models take as input.","marker":"Kennedy et al. (2019)"},{"why":"Supplies the orbital parameters for the quadruple system and the predicted ~2026 transit that the paper models.","marker":"Zúñiga-Fernández et al. (2021)"},{"why":"Provides the stellar fluxes used to normalise the synthetic light curves.","marker":"Ribas et al. (2018)"},{"why":"Supplies the MRN grain size distribution n(a) ~ a^-3.5 used to compute dust opacity.","marker":"Mathis et al. (1977)"},{"why":"Prior SPH simulations showing spiral-arm excitation by periastron passage, which this paper connects to light-curve asymmetry.","marker":"Smallwood et al. (2022)"},{"why":"Prior work extracting disc structure from an occultation light curve, the method this paper extends to HD 98800.","marker":"Kloppenborg et al. (2010)"}],"fun_headline_variants":["HD98800: predicted 8-year eclipse hasn't started yet","Disc size tweak explains missing HD98800 dimming","Spiral arms will break symmetry in HD98800's disc transit","HD98800's delayed dimming hints at lower viscosity disc","Transit forecast for HD98800: start delayed, end by 2034"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that dust is a passive tracer of the gas, with the disc's opacity obtained by scaling the simulated gas surface density using a fixed gas-to-dust ratio and a fixed grain size distribution; if the real dust has grown, migrated, or been depleted relative to the gas, the predicted transit timing and shape would be wrong.","fun_headline_variants_meta":{"raw":{"variants":["HD98800: predicted 8-year eclipse hasn't started yet","Disc size tweak explains missing HD98800 dimming","Spiral arms will break symmetry in HD98800's disc transit","HD98800's delayed dimming hints at lower viscosity disc","Transit forecast for HD98800: start delayed, end by 2034"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000264,"raw_usage":{"total_tokens":1649,"prompt_tokens":1039,"completion_tokens":610,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":655,"completion_tokens_details":{"reasoning_tokens":519}},"tokens_in":655,"tokens_out":610,"duration_ms":7252,"temperature":1.0,"reasoning_tokens":519,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T12:22:25.981050+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A photometric monitoring campaign at 0.5–0.8 μm that finds no ≥5% decrease in the A-binary flux by the end of 2025 would contradict the models' claim that ingress should be underway, pointing instead to a disc with a smaller radial extent or lower gas-to-dust ratio than simulated; a transit that begins but lasts much longer than 11 years would falsify the assumed outer-edge truncation.","supporting_citations":[{"cited_title":"M., Matr \\`a L., Facchini S., Milli J., Pani \\'c O., Price D., Wilner D","cited_arxiv_id":null,"evidence_quote":"Supplies the disc parameters (dust and gas mass lower bounds, radial extent 2.5–4.6 AU dust and 1.6–6.4 AU gas) and the polar-aligned interpretation that the models take as input."},{"cited_title":"C., Duch \\^ e ne G., 2018, , 865, 77","cited_arxiv_id":null,"evidence_quote":"Provides the stellar fluxes used to normalise the synthetic light curves."},{"cited_title":"L., Lubow S","cited_arxiv_id":null,"evidence_quote":"Prior SPH simulations showing spiral-arm excitation by periastron passage, which this paper connects to light-curve asymmetry."},{"cited_title":"D., Schaefer G., Zhao M., Baron F., McAlister H., Ten Brummelaar T., Che X., Farrington C., et al., 2010, , 464, 870","cited_arxiv_id":null,"evidence_quote":"Prior work extracting disc structure from an occultation light curve, the method this paper extends to HD 98800."}],"review_version":1}