Pith. sign in

REVIEW 3 major objections 4 minor 35 references

Accreting Binary Eccentricities follow Predicted Equilibrium Values

T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Actively accreting stellar binaries that are aligned with their circumbinary disks lie on the predicted equilibrium eccentricity curve.

desk verdict First stellar test of the Siwek attractor curve is a good idea, but the claimed fit rests on seven points with no error bars and a pair of systems that already look inconsistent. read the letter →

arxiv 2411.13489 v2 pith:3DFKXV2W submitted 2024-11-20 astro-ph.SR astro-ph.GA

Allen R Murray , Paul C Duffell This is my paper · ORCID
classification astro-ph.SRastro-ph.GA
keywords circumbinarydisksbinaryeccentricityaccretingbinariesequilibriumstellaraccretionbinary-diskinteractionmassratioprotostellar
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Accreting binaries that still have circumbinary disks occupy eccentricity–mass-ratio positions that match the equilibrium curve predicted by disk–binary hydrodynamics. The authors select from a catalog of binary–disk systems those with well-measured orbital parameters, no third body, mutual inclination consistent with zero, and separations too wide for tidal circularization; this aligned subset lies on the predicted curve. These are the first observational points tied to the predicted attractor eccentricity, rather than to the flat eccentricity distribution of the general Milky Way binary population. The paper concludes that eccentricity in such systems is set by disk accretion, and that some later mechanism must act after the disk dissipates to erase that preference.

What carries the argument

The load-bearing object is the equilibrium eccentricity curve: the locus of points in the $(q,e_b)$ plane at which the simulated eccentricity derivative vanishes, $\dot{e}_b(q,e_b)=0$, interpolated from the fixed-binary, two-dimensional hydrodynamic grid of Siwek et al. (2023). Combined with the mass-ratio evolution formula $\dot{q}(q)$ from Duffell et al. (2020), this curve yields a vector map showing binaries migrating onto the curve and then along it toward higher mass ratio. A second criterion, the Martin & Lubow (2019) critical inclination $i_{\min}=\arccos\sqrt{5e_b^2/[4(1+e_b^2)]}$, separates systems that should evolve to planar versus polar alignment, and the paper uses it to define the aligned subsample to which the eccentricity curve applies.

What would settle it

A larger sample of young circumbinary-disk binaries, selected without eccentricity or inclination priors, would refute the claim if the planar-aligned, non-tidally-circularized systems scatter far from the $\dot{e}_b=0$ curve; specifically, if a majority of a dozen new systems fall more than about 0.1 in eccentricity from the Siwek et al. attractor, the equilibrium interpretation would be unsupported.

Watch

Extended reading notes

Core claim

The central discovery is that the observed eccentricities of planar-aligned accreting stellar binaries trace the zero-derivative curve $\dot{e}_b(q,e_b)=0$ computed by Siwek et al. (2023), meaning the binaries have evolved to the eccentricity at which disk torques stop changing it. The aligned systems also sit below the Martin & Lubow (2019) critical inclination curve, as they must if they are to remain aligned, while the two polar-aligned systems sit above the curve and the hierarchical triples scatter outside the simple binary–disk predictions. The authors read the match as evidence that the disks have been present long enough to drive the binaries onto the attractor; the same plot explains why the general field population, whose disks have long since vanished, shows no preferred eccentricity.

Load-bearing premise

The aligned systems have already accreted enough gas for disk torques to have driven their eccentricities to the equilibrium curve, rather than being born near it by chance or selected into the catalog because they lie near it.

Editorial extensions

If this is right

  • A real match would show that circumbinary disk torques reshape stellar binary orbits on observable timescales, providing a local proxy for the same physics predicted for supermassive black hole binaries.
  • Because the merged Milky Way binary population has a flat eccentricity distribution, some mechanism must act after the disk dissipates to erase the preferred eccentricity–mass-ratio pattern.
  • The position of a system like Ak Sco, with $q\approx 1$, $e_b\approx 0.47$, and an age of 18 Myr, is naturally read as a 'twin' binary that has accreted substantially, supporting the predicted mass-ratio evolution.
  • A proposed test using astrometric surveys of young binaries with no third body should recover the same attractor eccentricity if the equilibrium interpretation is correct.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the aligned subsample contains only a handful of systems drawn from a catalog not designed for this test, part of the agreement could reflect selection effects; a blind, volume-limited survey of accreting binaries would sharpen the comparison.
  • The theoretical curve comes from two-dimensional, fixed-binary simulations, while real protoplanetary disks are three-dimensional and the binary's orbit responds to the disk; the level of agreement may be partly forgiving of those simplifications.
  • If the equilibrium is genuine, the same eccentricity–mass-ratio attractor should appear in other circumbinary-disk populations, such as post-asymptotic-giant-branch binaries, where a bimodal eccentricity distribution has already been predicted.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper assembles the well-characterized binary–disk systems from Czekala et al. (2019) and classifies them by mutual binary–disk inclination into tidally circularized, planar aligned, hierarchical trinary, and polar aligned categories, using the Martin–Lubow critical-inclination criterion. For the planar aligned systems, the authors compare the observed eccentricity and mass ratio against the equilibrium (attractor) eccentricity curve of Siwek et al. (2023) and claim that the systems lie on that curve, suggesting that they have accreted enough gas for circumbinary-disk torques to set their eccentricity, and that eccentricity evolution resumes only after disk dispersal.

Significance. If the central comparison is established, this would be a valuable observational test of circumbinary-disk eccentricity theory that was developed largely for supermassive black hole binaries, and it would add an interesting post-disk evolutionary constraint for young stellar binaries. The paper has genuine strengths: the theoretical curve is taken from published simulations rather than fitted to the data, the aligned/polar classification uses an independent published criterion, and the sample is small but carefully selected from a homogeneous catalog. The main significance, however, hinges entirely on the quantitative agreement in Figure 2, and that agreement is currently asserted visually rather than demonstrated.

major comments (3)
  1. [§4, Figure 2] The central claim that the aligned systems lie on the Siwek et al. (2023) attractor curve is not quantified. The paper reports no observational uncertainties for eccentricity or mass ratio, no residual statistic, and no uncertainty on the interpolated theoretical curve. This is not a cosmetic omission: TWA 3A (e=0.628, q=0.84) and HD 200775 (e=0.30, q=0.82) have nearly identical mass ratios but differ in eccentricity by 0.33. Since e_attr(q) is a single-valued, monotonically increasing curve over this q range in Siwek et al., no curve of that shape can pass within a small tolerance of both systems. Either the error bars on e are large enough (about 0.15 or more) that the match carries little information, or at least one of these two systems is inconsistent with the predicted curve as drawn. The authors should provide eccentricity and mass-ratio errors, a measure of scatter or a goodness-of-fit statistic, and the uncertainty on the interpolated attractor curve, so that the claimed agreement can actually be evaluated.
  2. [Table 1 and §4] There is an internal counting inconsistency that suggests the fit was not independently checked: the text in §4 states that 'the eight aligned systems lie approximately upon the attractor eccentricity curve,' but Table 1 lists only seven planar aligned systems (HD 131511, αCrB, TWA 3A, HD 200775, Ak Sco, DQ Tau, UZ Tau E). The two tidally circularized systems (V4046 Sgr, CoRoT 2239) have e=0 and are explicitly placed in a separate category, so they cannot be the eighth aligned system. The paper should correct the count and clarify exactly which systems are used in the Figure 2 fit.
  3. [§2 and §4] The interpretation that the aligned systems 'have accreted enough gas to significantly affect their orbit' is inferred from the match itself, because the sample selection does not include an accretion history or age criterion that would independently identify systems that have reached the equilibrium. As written, the agreement could also arise if these young binaries were simply born near the curve or if the catalog selection favors such systems. To make the evolutionary claim load-bearing, the authors should either add an independent check (for example, comparing accretion indicators or estimated accreted mass against the predicted equilibration timescale) or explicitly reframe the conclusion as a statement about the current e–q distribution only.
minor comments (4)
  1. [§4 and §5] The wording oscillates between 'consistent with' (abstract), 'lie approximately upon' (§4), and 'precisely obeys' (§5); the final wording is too strong given the lack of quantified residuals and should be harmonized with the actual statistical support.
  2. [§3, Equation (1)] Equation (1) is typeset in a garbled way in the manuscript text, and the displayed formula is hard to read; it should be reset so that the argument of the arccosine and the eccentricity dependence are unambiguous.
  3. [§3] The sentence describing the interpolation of points where ė_b=0 and the interpolation for e_b=0.7 should specify whether the interpolation is linear in the (q, e) plane and whether the resulting curve has any systematic uncertainty from the coarse tabulation in Siwek et al. (2023).
  4. [§2, Table 1] The procedure of setting i=0 and using the upper limit as the uncertainty for systems with i<10° is reasonable, but the same treatment is not applied to the eccentricity and mass-ratio columns, which are listed without any uncertainties; providing the original error bars from Czekala et al. (2019) in a supplementary table or in the caption would greatly improve the transparency of the comparison.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the equilibrium curve is an independent simulation prediction (Siwek et al. 2023) compared with external orbital measurements; only minor non-load-bearing self-citations appear.

full rationale

The claimed result is that seven planar-aligned, non-tidally-circularized stellar binaries lie on the equilibrium eccentricity curve e_attr(q) predicted by Siwek et al. (2023). The derivation chain is: (i) read (e, q, i) from the Czekala et al. (2019) catalog; (ii) take the edot_b = 0 locus from the published Siwek et al. (2023) hydrodynamical simulation table and interpolate it; (iii) plot both together and visually compare. The curve is not fitted to the data, no parameter is tuned to the observed points, and the observed eccentricities and mass ratios are independent external measurements. The paper's self-citations (Duffell et al. 2020 for qdot(q); D'Orazio & Duffell 2021 for the attractor) are not load-bearing: the specific attractor curve is from Siwek et al. (2023), and the same attractor is supported by Zrake et al. (2021) and Valli et al. (2024), so the argument does not reduce to a self-citation chain. The sample restrictions (excluding tidally circularized systems and hierarchical triples) are stated scope conditions motivated by the theory's assumptions, not definitions that guarantee the conclusion. No equation in the paper is equal to its input by construction. The absence of a quantitative goodness-of-fit test is a verification weakness, not a circularity. Overall circularity score 2 reflects only the presence of minor, non-load-bearing self-citations.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No free parameters are introduced in this paper; the attractor curve and mass-ratio evolution are imported from published simulations. The load-bearing assumptions are domain assumptions about the applicability of 2D fixed-binary simulations to real disks and about the accretion history of the sample.

assumptions (5)
  • domain assumption The critical inclination criterion of Martin and Lubow (2019), Equation (1), correctly predicts whether a disk evolves to planar or polar alignment for negligible disk mass.
    Used in Figure 1 to classify systems as aligned or polar; if wrong, the sample splits are invalid.
  • domain assumption The 2D fixed-binary simulation attractor curves of Siwek et al. (2023) apply to real 3D circumbinary disks.
    Section 3 explicitly adopts these predictions while noting they rely on fixed binary 2D simulations; Valli et al. (2024) is cited for consistency.
  • domain assumption Kozai-Lidov eccentricity effects are small for mutual inclinations below about 20 degrees.
    Section 3 cites Martin and Lubow (2019) to justify using zero-inclination predictions for the aligned sample.
  • domain assumption Systems without traditional inclination error bars (i less than 10 degrees) are treated as aligned with i equal to 0.
    Section 2 defines this assignment; the central aligned sample depends on it.
  • domain assumption The aligned systems have accreted enough gas to have reached the equilibrium eccentricity.
    Section 4 infers this from the systems' location on the curve; it is not independently measured.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Accreting Binary Eccentricities follow Predicted Equilibrium Values." pith.science (2026). https://pith.science/paper/3DFKXV2W

@misc{pith2026241113489,
  author       = {Pith},
  title        = {Pith review of: Accreting Binary Eccentricities follow Predicted Equilibrium Values},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3DFKXV2W}},
  note         = {Machine review of arXiv:2411.13489}
}
read the original abstract

We investigate observations of circumbinary disks (CBD), to find evidence for an equilibrium eccentricity predicted by current binary accretion theory. Although stellar binary demographics in the Milky Way show no evidence for a preferred eccentricity for binary systems, we show that actively accreting systems lie on a predicted equilibrium eccentricity curve. We constrain our sample to only systems that have well defined orbital parameters (e.g,. eccentricity, mass-ratio, inclination angle). We find observations are consistent with theory for stellar binaries that are aligned with the disk and that are separated enough that tidal circularization is negligible. This suggests that eccentricity in these systems evolves after the dissipation of the CBD, given the flat eccentricity distribution of binary systems in the Milky Way.

Figures

Figures reproduced from arXiv: 2411.13489 by the authors.

Figure 1
Figure 1. Studies that allow live binaries and 3D effects find that for binaries with small inclination, i  20o , eccentricity effects due to Kozai–Lidov oscillations are small. (R. G. Martin & S. H. Lubow 2019). Using the first table in [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. For long-term evolution of binary–disk systems considering semimajor axis, mass ratio, and eccentricity of the binary, see R. Valli et al. (2024). R. Valli et al. (2024) look at simulation results from three separate studies (D. J. D’Orazio & P. C. Duffell 2021; J. Zrake et al. 2021; M. Siwek et al. 2023) and test consistency in CBD dynamics between the different parameter spaces used. They found that the different … view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

35 extracted references · 27 canonical work pages

  1. [1]

    Alencar, S. H. P., Melo, C. H. F., Dullemond, C. P., et al. 2003, A&A, 409, 1037

    work page 2003

  2. [2]

    2015, A&A, 574, A41

    Anthonioz, F., Ménard, F., Pinte, C., et al. 2015, A&A, 574, A41

    work page 2015

  3. [3]

    R., & Bonnell, I

    Bate, M. R., & Bonnell, I. A. 1997, MNRAS, 285, 33

    work page 1997

  4. [4]

    R., Lodato, G., & Pringle, J

    Bate, M. R., Lodato, G., & Pringle, J. E. 2010, MNRAS, 401, 1505

    2010

  5. [5]

    C., Blandford, R

    Begelman, M. C., Blandford, R. D., & Rees, M. J. 1980, Natur, 287, 307

    work page 1980

  6. [6]

    Boss, A. P. 1986, ApJS, 62, 519

    work page 1986

  7. [7]

    Ceppi, S., Longarini, C., Lodato, G., Cuello, N., & Lubow, S. H. 2023, MNRAS, 520, 5817

    work page 2023

  8. [8]

    2013, MNRAS, 429, 3114

    Chapon, D., Mayer, L., & Teyssier, R. 2013, MNRAS, 429, 3114

    work page 2013

Show all 35 references
  1. [9]

    J., Alexander, R

    Cuadra, J., Armitage, P. J., Alexander, R. D., & Begelman, M. C. 2009, MNRAS, 393, 1423

    2009

  2. [10]

    M., Jensen, E

    Czekala, I., Andrews, S. M., Jensen, E. L. N., et al. 2015, ApJ, 806, 154

    2015

  3. [11]

    M., et al

    Czekala, I., Chiang, E., Andrews, S. M., et al. 2019, ApJ, 883, 22 D’Orazio, D. J., & Duffell, P. C. 2021, ApJL, 914, L21

    2019

  4. [12]

    2007, MNRAS, 379, 956

    Dotti, M., Colpi, M., Haardt, F., & Mayer, L. 2007, MNRAS, 379, 956

    2007

  5. [13]

    C., D ’Orazio, D., Derdzinski, A., et al

    Duffell, P. C., D ’Orazio, D., Derdzinski, A., et al. 2020, ApJ, 901, 25

    2020

  6. [14]

    2019, MNRAS, 489, 5822

    El-Badry, K., Rix, H.-W., Tian, H., Duchêne, G., & Moe, M. 2019, MNRAS, 489, 5822

    2019

  7. [15]

    2018, MNRAS, 476, 528

    El-Badry, K., Ting, Y.-S., Rix, H.-W., et al. 2018, MNRAS, 476, 528

    2018

  8. [16]

    B., Coppi, P

    Escala, A., Larson, R. B., Coppi, P. S., & Mardones, D. 2005, ApJ, 630, 152

    2005

  9. [17]

    D., Duffell, P., MacFadyen, A

    Farris, B. D., Duffell, P., MacFadyen, A. I., & Haiman, Z. 2014, ApJ, 783, 134

    2014

  10. [18]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Natur, 585, 357

    2020

  11. [19]

    Hunter, J. D. 2007, CSE, 9, 90

    2007

  12. [20]

    2011, MNRAS, 417, 1466

    Kashi, A., & Soker, N. 2011, MNRAS, 417, 1466

    2011

  13. [21]

    M., Matzner, C

    Kratter, K. M., Matzner, C. D., & Krumholz, M. R. 2008, ApJ, 681, 375

    2008

  14. [22]

    G., & Lubow, S

    Lepp, S., Martin, R. G., & Lubow, S. H. 2023, ApJL, 943, L4

    2023

  15. [23]

    G., & Lubow, S

    Martin, R. G., & Lubow, S. H. 2017, ApJL, 835, L28

    2017

  16. [24]

    G., & Lubow, S

    Martin, R. G., & Lubow, S. H. 2018, MNRAS, 479, 1297

    2018

  17. [25]

    G., & Lubow, S

    Martin, R. G., & Lubow, S. H. 2019, MNRAS, 490, 1332

    2019

  18. [26]

    2007, Sci, 316, 1874 Milosavljević, M., & Merritt, D

    Mayer, L., Kazantzidis, S., Madau, P., et al. 2007, Sci, 316, 1874 Milosavljević, M., & Merritt, D. 2001, ApJ, 563, 34 Milosavljević, M., & Phinney, E. S. 2005, ApJL, 622, L93

    2007

  19. [27]

    L., Clarke, C

    Monin, J. L., Clarke, C. J., Prato, L., & McCabe, C. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 395 (Tucson, AZ: Univ. Arizona Press )

    2007

  20. [28]

    Moody, M. S. L., Shi, J.-M., & Stone, J. M. 2019, ApJ, 875, 66 Muñoz, D. J., Lai, D., Kratter, K., & Miranda, R. 2020, ApJ, 889, 114 Muñoz, D. J., Miranda, R., & Lai, D. 2019, ApJ, 871, 84 The pandas development team 2020, pandas-dev /pandas: Pandas, v2.2.2, Zenodo, doi: 10.52...

    2019 doi

  21. [29]

    M., Hogg, D

    Price-Whelan, A. M., Hogg, D. W., Rix, H.-W., et al. 2020, ApJ, 895, 2

    2020

  22. [30]

    Z., & Hernquist, L

    Siwek, M., Kelley, L. Z., & Hernquist, L. 2024, MNRAS, 534, 2609

    2024

  23. [31]

    2023, MNRAS, 522, 2707

    Siwek, M., Weinberger, R., & Hernquist, L. 2023, MNRAS, 522, 2707

    2023

  24. [32]

    2017, MNRAS, 469, 4258

    Tang, Y., MacFadyen, A., & Haiman, Z. 2017, MNRAS, 469, 4258

    2017

  25. [33]

    Tuna, S., & Metzger, B. D. 2023, ApJ, 955, 125

    2023

  26. [34]

    2024, A&A, 688, A128

    Valli, R., Tiede, C., Vigna-Gómez, A., et al. 2024, A&A, 688, A128

    2024

  27. [35]

    2021, ApJL, 909, L13 4 The Astrophysical Journal, 982:113 (4pp), 2025 April 1 Murray & Duffell

    Zrake, J., Tiede, C., MacFadyen, A., & Haiman, Z. 2021, ApJL, 909, L13 4 The Astrophysical Journal, 982:113 (4pp), 2025 April 1 Murray & Duffell

    2021

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.