REVIEW 3 major objections 5 minor 1 cited by
Modelling a Transiting Circumbinary Disc in the HD98800 System
T0 review · 3 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [§2.1, Table 3] 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.
- [§2.3, Eq. (1), §3.7] 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.
- [§3.6] 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.
minor comments (5)
- [Abstract, §1, §2.1] 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').
- [§2.2.2] The sentence 'observations ... have also provided estimates of alpha that suggest a that alpha can vary' contains a stray word and should be rewritten.
- [§3.6] 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.
- [Figure 6 caption] The caption states that stars are 'plotted on top' although they pass behind the disc; this is visually confusing and should be clarified.
- [§3.3] 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.
Circularity Check
No significant circularity: synthetic light curves are forward-model predictions from externally calibrated inputs; no fitted parameter is renamed as a prediction.
full rationale
The derivation chain is a forward model. Orbital elements (Zúñiga-Fernández et al. 2021) and disc properties (masses, radial extents; Kennedy et al. 2019) enter as inputs from independent ALMA/astrometric observations; the SPH runs evolve the gas; the dust is added in post-processing by scaling the gas surface density with a fixed gas-to-dust ratio and MRN grain-size distribution (Eq. 1); mcfost then computes optical depths and synthetic light curves (Eq. 2). The transit start/end dates in Table 3 come from applying a 5% flux-drop threshold to these light curves, not from fitting any parameter to the transit itself. The dust masses were calibrated to the system's SED (Fig. 3), but the SED is an integrated observable and does not encode the ingress timing, transit duration, spiral-arm asymmetries, or wavelength dependence that the paper presents; those are independent forward predictions. The cited prior work by the same group is used for observational constraints (ALMA disc sizes/masses, orbital parameters), which are externally checkable and not arguments whose only support is the present paper. The caveats in Sec. 3.7 (no dust evolution, SPH viscosity floor ~0.01, outer edge initialised at 6 AU) mean the no-dimming conclusion is an extrapolation beyond the simulated viscosity range, and the 6 AU edge is an initial condition rather than a converged viscous result; but an extrapolation from an assumed initial condition is a model limitation, not a circular reduction. No equation defines a predicted quantity in terms of itself, and no fitted value is relabelled as a prediction. Therefore there is no significant circularity; score 0.
Assumptions & free parameters
free parameters (8)
- Dust mass M_dust =
0.033, 0.33, 3.3 M_Earth
- Gas mass M_gas =
3.3, 33, 330 M_Earth
- Alpha viscosity =
0.005, 0.01, 0.05
- Initial outer disc edge =
6.0 AU
- Surface density slope p =
1
- Flaring index q =
0.25
- Aspect ratio H/R =
0.05 at 2.5 AU
- 5% flux drop threshold =
0.05 fractional flux decrement
assumptions (7)
- domain assumption Dust grains are small and perfectly coupled to the gas; dust surface density is gas surface density scaled by the gas-to-dust ratio.
- domain assumption An MRN power-law grain size distribution n(a) ~ a^-3.5 over 0.03 to 1000 micron describes the dust.
- domain assumption SPH numerical viscosity floors alpha ~ 0.005 to 0.01; lower viscosities cannot be simulated.
- domain assumption Orbital elements of AaAb, BaBb, and AB from Zuniga-Fernandez et al. (2021) are accurate, especially the AB orbit until 2050.
- domain assumption The disc is initialized in a relaxed state; after 1150 years of evolution, transient features are washed out.
- domain assumption Opacity is dominated by small grains with Stokes number much less than unity, validating the trace-gas approximation.
- standard math Hydrodynamics and radiative transfer equations implemented in phantom and mcfost are correct.
Cite this review
Pith. "Pith review of Modelling a Transiting Circumbinary Disc in the HD98800 System." pith.science (2026). https://pith.science/paper/HMJHVTXF
@misc{pith2026250202376,
author = {Pith},
title = {Pith review of: Modelling a Transiting Circumbinary Disc in the HD98800 System},
year = {2026},
howpublished = {\url{https://pith.science/paper/HMJHVTXF}},
note = {Machine review of arXiv:2502.02376}
}
abstract
We present synthetic optical light curves of the hierarchical HD98800 quadruple system over a decade-long period when the circumbinary disc encircling the system's B binary is expected to eclipse the light from the A binary. We produce and compare light curves of this transit event using hydrodynamical models with different values of the disc's gas mass, dust mass, and $\alpha$-viscosity to determine the observable effect of each parameter. These comparisons provide insight that could aid in the analysis of observational data from the system when the real transit occurs and provide recommendations for how such observations should be made. We find that a higher dust mass or higher value of $\alpha$ correspond to a longer transit, with the gas mass having a more minor effect on the overall shape and duration of the transit. A higher $\alpha$ has an observable effect on the viscous spreading at the outer edge of the disc, though is countered through truncation by the outer binary. It is also shown that long-term interactions between the outer binary and disc can excite spiral arms in the disc, which introduce observable asymmetries to the light curve. Our models suggest that the transit should have begun at the time of writing, but no dimming has yet been observed. It is likely that the disc has a smaller radial extent than our models, due to a lower viscosity than can be simulated with SPH. The transit is expected to last 8-11 years, ending in late 2034 at the latest.
Figures
Figures from the paper (9 more)
Forward citations
Cited by 1 Pith paper
-
Peering through the disc of HD 98800 BaBb. Precise timing predictions for the HD 98800 AaAb occultation
New RVs and multi-wavelength astrometry tighten the HD 98800 outer orbit by ~2×, yielding 5–15-day 1σ windows for the 2025–2031 disc occultation of AaAb.
Reference graph
Works this paper leans on
-
[1]
Akeson R. L., Rice W. K. M., Boden A. F., Sargent A. I., Carpenter J. M., Bryden G., 2007, , 670, 1240
work page 2007
-
[2]
Aly H., Gonzalez J.-F., Nealon R., Longarini C., Lodato G., Price D. J., 2021, , 508, 2743
work page 2021
- [3]
-
[4]
Andrews S. M., Czekala I., Wilner D. J., Espaillat C., Dullemond C. P., Hughes A. M., 2010, , 710, 462
work page 2010
-
[5]
Andrews S. M., Huang J., P \'e rez L. M., Isella A., Dullemond C. P., Kurtovic N. T., Guzm \'a n V. V., Carpenter J. M., Wilner D. J., Zhang S., Zhu Z., Birnstiel T., Bai X.-N., Benisty M., Hughes A. M., \"O berg K. I., Ricci L., 2018, , 869, L41
work page 2018
- [6]
- [7]
-
[8]
Beck T. L., Schaefer G. H., Guilloteau S., Simon M., Dutrey A., Di Folco E., Chapillon E., 2020, , 902, 132
work page 2020
Show all 53 references
-
[9]
C., Savage B
Bohlin R. C., Savage B. D., Drake J. F., 1978, , 224, 132
1978
-
[10]
M., Pani \'c O., Haworth T
Boneberg D. M., Pani \'c O., Haworth T. J., Clarke C. J., Min M., 2016, , 461, 385
2016
-
[11]
E., Henning T., Pavlyuchenkov Y., Dullemond C
Carmona A., van den Ancker M. E., Henning T., Pavlyuchenkov Y., Dullemond C. P., Goto M., Thi W. F., Bouwman J., Waters L. B. F. M., 2008, , 477, 839
2008
-
[12]
M., Jensen E
Czekala I., Chiang E., Andrews S. M., Jensen E. L. N., Torres G., Wilner D. J., Stassun K. G., Macintosh B., 2019, , 883, 22
2019
-
[13]
Dipierro G., Price D., Laibe G., Hirsh K., Cerioli A., Lodato G., 2015, , 453, L73
2015
-
[14]
M., McCabe C., Ceccarelli C., 2005, , 628, 832
Duch \^e ne G., Ghez A. M., McCabe C., Ceccarelli C., 2005, , 628, 832
2005
-
[15]
P., Birnstiel T., Huang J., Kurtovic N
Dullemond C. P., Birnstiel T., Huang J., Kurtovic N. T., Andrews S. M., Guzm \'a n V. V., P \'e rez L. M., Isella A., Zhu Z., Benisty M., Wilner D. J., Bai X.-N., Carpenter J. M., Zhang S., Ricci L., 2018, , 869, L46
2018
-
[16]
H., Martin R
Franchini A., Lubow S. H., Martin R. G., 2019, , 880, L18
2019
-
[17]
M., Bolton C., et al., 2012, , 544, A53
Ga an C., Miko ajewski M., Tomov T., Graczyk D., Apostolovska G., Barzova I., Bellas-Velidis I., Bilkina B., Blake R. M., Bolton C., et al., 2012, , 544, A53
2012
-
[18]
o hler R., Laugier R., Martinache F., Siebenmorgen R., van den Ancker M. E., van Boekel R., Herbst T. M., Pantin E., K \
Kammerer J., Kasper M., Ireland M. J., K \"o hler R., Laugier R., Martinache F., Siebenmorgen R., van den Ancker M. E., van Boekel R., Herbst T. M., Pantin E., K \"a ufl H. U., Petit dit de la Roche D. J. M., Ivanov V. D., 2021, , 646, A36
2021
-
[19]
H., Zuckerman B., Weintraub D
Kastner J. H., Zuckerman B., Weintraub D. A., Forveille T., 1997, Science, 277, 67
1997
-
[20]
M., Matr \`a L., Facchini S., Milli J., Pani \'c O., Price D., Wilner D
Kennedy G. M., Matr \`a L., Facchini S., Milli J., Pani \'c O., Price D., Wilner D. J., Wyatt M. C., Yelverton B. M., 2019, Nature Astronomy, 3, 230
2019
-
[21]
D., Schaefer G., Zhao M., Baron F., McAlister H., Ten Brummelaar T., Che X., Farrington C., et al., 2010, , 464, 870
Kloppenborg B., Stencel R., Monnier J. D., Schaefer G., Zhao M., Baron F., McAlister H., Ten Brummelaar T., Che X., Farrington C., et al., 2010, , 464, 870
2010
-
[22]
K \"o hler R., Kubiak K., 2020, Research Notes of the American Astronomical Society, 4, 73
2020
-
[23]
J., Harsono D., Dipierro G., Pascucci I., Hendler N., Tazzari M., Ragusa E., Salyk C., Edwards S., Lodato G., van de Plas G., Johnstone D., Liu Y., Boehler Y., Cabrit S., Manara C
Long F., Pinilla P., Herczeg G. J., Harsono D., Dipierro G., Pascucci I., Hendler N., Tazzari M., Ragusa E., Salyk C., Edwards S., Lodato G., van de Plas G., Johnstone D., Liu Y., Boehler Y., Cabrit S., Manara C. F., Menard F., Mulders G. D., Nisini B., Fischer W. J., Rigliaco...
2018
-
[24]
S., Rumpl W., Nordsieck K
Mathis J. S., Rumpl W., Nordsieck K. H., 1977, , 217, 425
1977
-
[25]
D., Dominik C., 2012, , 539, A9
Mulders G. D., Dominik C., 2012, , 539, A9
2012
-
[26]
Papaloizou J. C. B., Nelson R. P., 2003, , 339, 983
2003
-
[27]
Pinte C., Dent W. R. F., M \'e nard F., Hales A., Hill T., Cortes P., de Gregorio-Monsalvo I., 2016, , 816, 25
2016
-
[28]
Pinte C., Harries T., Min M., Watson A., Dullemond C., Woitke P., M \'e nard F., Dur \'a n-Rojas M., 2009, , 498, 967
2009
-
[29]
Pinte C., M \'e nard F., Duch \^e ne G., Bastien P., 2006, , 459, 797
2006
-
[30]
L., Skillen I., Collier Cameron A., Christian D
Pollacco D. L., Skillen I., Collier Cameron A., Christian D. J., Hellier C., Irwin J., Lister T. A., Street R. A., West R. G., Anderson D. R., Clarkson W. I., Deeg H., Enoch B., Evans A., Fitzsimmons A., Haswell C. A., Hodgkin S., Horne K., Kane S. R., Keenan F. P., Maxted P. ...
2006
-
[31]
J., 2007, , 24, 159
Price D. J., 2007, , 24, 159
2007
-
[32]
J., Wurster J., Tricco T
Price D. J., Wurster J., Tricco T. S., Nixon C., Toupin S., Pettitt A., Chan C., Mentiplay D., Laibe G., Glover S., et al., 2018, , 35
2018
-
[33]
R., 2017, , 837, 163
Rafikov R. R., 2017, , 837, 163
2017
-
[34]
A., Henry T
Raghavan D., McAlister H. A., Henry T. J., Latham D. W., Marcy G. W., Mason B. D., Gies D. R., White R. J., Theo A., 2010, , 190, 1
2010
-
[35]
A., Leinert C., \'A brah \'a m P., Henning T., Herbst T
Ratzka T., Schegerer A. A., Leinert C., \'A brah \'a m P., Henning T., Herbst T. M., K \"o hler R., Wolf S., Zinnecker H., 2009, , 502, 623
2009
-
[36]
F., 2012, , 537, A128
Rein H., Liu S. F., 2012, , 537, A128
2012
-
[37]
C., Duch \^ e ne G., 2018, , 865, 77
Ribas \' A ., Mac \' as E., Espaillat C. C., Duch \^ e ne G., 2018, , 865, 77
2018
-
[38]
E., Stassun K
Rodriguez J. E., Stassun K. G., Lund M. B., Siverd R. J., Pepper J., Tang S., Kafka S., Gaudi B. S., Conroy K. E., Beatty T. G., Stevens D. J., Shappee B. J., Kochanek C. S., 2016, , 151, 123
2016
-
[39]
P., Guilera O
Ronco M. P., Guilera O. M., Cuadra J., Miller Bertolami M. M., Cuello N., Fontecilla C., Poblete P., Bayo A., 2021, , 916, 113
2021
-
[40]
I., Sunyaev R
Shakura N. I., Sunyaev R. A., 1973, , 24, 337
1973
-
[41]
L., Lubow S
Smallwood J. L., Lubow S. H., Martin R. G., 2022, , 514, 1249
2022
-
[42]
R., King J
Soderblom D. R., King J. R., Siess L., Noll K. S., Gilmore D. M., Henry T. J., Nelan E., Burrows C. J., Brown R. A., Perryman M. A. C., Benedict G. F., McArthur B. J., Franz O. G., Wasserman L. H., Jones B. F., Latham D. W., Torres G., Stefanik R. P., 1998, , 498, 385
1998
-
[43]
P., Barrado y Navascu \'e s D., Stauffer J
Song I., Caillault J. P., Barrado y Navascu \'e s D., Stauffer J. R., Randich S., 2000, , 533, L41
2000
-
[44]
Steinacker A., Papaloizou J. C. B., 2002, , 571, 413
2002
-
[45]
P., Wilner D
Testi L., Birnstiel T., Ricci L., Andrews S., Blum J., Carpenter J., Dominik C., Isella A., Natta A., Williams J. P., Wilner D. J., 2014, in Beuther H., Klessen R. S., Dullemond C. P., Henning T., eds, Protostars and Planets VI Dust Evolution in Protoplanetary Disks . pp 339--361
2014
-
[46]
P., Latham D
Torres G., Stefanik R. P., Latham D. W., Mazeh T., 1995, , 452, 870
1995
-
[47]
Van Leeuwen F., 2007, , 474, 653
2007
-
[48]
van Werkhoven T. I. M., Kenworthy M. A., Mamajek E. E., 2014, , 441, 2845
2014
-
[49]
J., 1977, , 180, 57
Weidenschilling S. J., 1977, , 180, 57
1977
-
[50]
Whipple F., , 1972, From Plasma to Planet ed A. Elvius
1972
-
[51]
L., 1973, in Hemenway C
Whipple F. L., 1973, in Hemenway C. L., Millman P. M., Cook A. F., eds, , Vol. 319, NASA Special Publication. p. 355
1973
-
[52]
Z \'u \ n iga-Fern \'a ndez S., Olofsson J., Bayo A., Haubois X., Corral-Santana J., Lopera-Mej \' a A., Ronco M., Tokovinin A., Gallenne A., Kennedy G., et al., 2021, , 655, A15
2021
-
[53]
write newline
" write newline "" before.all 'output.state := FUNCTION fin.entry write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = 'skip output.state before.all = 'skip after.sentence 'output.stat...
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.