pith. sign in

arxiv: 2510.24127 · v2 · pith:KR6AZXB2new · submitted 2025-10-28 · 🌌 astro-ph.SR · astro-ph.HE

Radiatively Cooled Binary Mass Transfer: Flow Structure, Luminosities, and L2 Outflows Across Mass Transfer Rates

Pith reviewed 2026-05-21 20:53 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.HE
keywords binary mass transferL2 outflowshydrodynamical simulationsradiative coolingcircumbinary outflowsstellar transientsmass transfer ratesangular momentum
0
0 comments X

The pith

In binary mass transfer, significant outflows through the L2 point occur only for rates above about 0.001 solar masses per year, carrying the specific angular momentum of L2.

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

Hydrodynamical simulations of a 10 solar mass donor and 5 solar mass accretor include approximate radiative cooling to model mass transfer across a range of rates. The results show that mass flows into a disk around the accretor. Significant equatorially concentrated outflows through L2 happen at high mass transfer rates of 10^{-3} solar masses per year and above, while lower rates result in mostly conservative transfer. Outflowing gas always approximately carries the angular momentum of L2. Cooling luminosity and temperature rise with increasing mass transfer rate, reaching 10^5 solar luminosities and 10^4 K in strong outflow cases, contributing to luminous transients.

Core claim

Mass flows from the donor into a disk around the accretor, with significant equatorially concentrated outflows through the outer Lagrange point L2 occurring for MT rates ≳ 10^{-3} M_⊙/yr, while the MT remains mostly conservative for lower MT rates. In all cases, any outflowing gas approximately carries the specific angular momentum of L2. The gas cooling luminosity L and temperature increases with MT rate, with L ∼ 10^5 L_⊙ and T ∼ 10^4 K for simulations featuring the strongest outflows, with contributions from both the CBO and the accretor's disk.

What carries the argument

Hydrodynamical simulations incorporating approximate radiative cooling, run at varying orbital separations to achieve mass transfer rates from 10^{-5} to 10^{-1} solar masses per year and track the resulting stream, disk, and outflow structures.

If this is right

  • Mass transfer remains mostly conservative for rates below 10^{-3} solar masses per year.
  • Significant equatorially concentrated outflows through L2 develop at rates of 10^{-3} solar masses per year and higher.
  • Any outflowing gas carries approximately the specific angular momentum of L2 regardless of rate.
  • Cooling luminosity and temperature increase with mass transfer rate, reaching roughly 10^5 solar luminosities and 10^4 K in the strongest outflow cases.
  • Contributions to the total luminosity come from both the circumbinary outflow and the disk around the accretor.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Observed circumbinary material in some stellar transients may trace back to these L2 outflows at high mass transfer rates.
  • The angular momentum removed by L2 outflows could speed up orbital shrinkage or merger in the binary.
  • Optical light curves from transients could test the predicted rise in luminosity and temperature with mass transfer rate.
  • Similar flow patterns and luminosity scaling may hold in binaries with other mass ratios if the cooling model remains valid.

Load-bearing premise

The approximate radiative cooling prescription used in the hydrodynamical simulations sufficiently captures the thermal and dynamical evolution of the mass-transfer stream and disk across the explored range of orbital separations.

What would settle it

A direct measurement or comparison showing whether equatorially concentrated outflows through L2 become significant precisely at mass transfer rates around 10^{-3} solar masses per year and whether the specific angular momentum of those outflows matches L2.

Figures

Figures reproduced from arXiv: 2510.24127 by Jim Fuller, Peter Scherbak, Wenbin Lu.

Figure 1
Figure 1. Figure 1: The density ρ in the equatorial plane, for simulations of varying MT rate. The orbital separation a and Lagrange points are labeled. The snapshots, representative of the simulations once they have reached a quasi-steady state, are labeled with t in both code units and physical units (days). The full domain extends to r = 5a, but these plots extend to r = 3a to zoom in near the accretion disk. Note the diff… view at source ↗
Figure 3
Figure 3. Figure 3: The predicted fraction of material escaping through L2, fL2, for various a and MT rates. This follows the figures of Lu et al. (2023), but adapted for a 10 M⊙ primary and 5 M⊙ primary. Solid lines show contours of constant fL2, and the cyan dotted line shows where M˙ donor = M˙ edd. Stars show where our simulations lie in this parameter space, with the internal color corresponding to their fL2 value. The h… view at source ↗
Figure 2
Figure 2. Figure 2: Top panel: The fraction of material retained in the accretion disk, β, versus the MT rate M˙ donor for our simulations. Bottom panel: the same quantity, plotted ver￾sus M˙ donor in units of the Eddington MT rate (Eq. 11). The vertical dashed line shows where M˙ donor = M˙ edd. The values shown are time-averaged once the simulations reach a quasi￾steady state, at t greater than ∼ 300 (code units). opacity v… view at source ↗
Figure 4
Figure 4. Figure 4: The cooling luminosity calculated within the ac￾cretion disk (dashed lines, Eq. 12, and the entire computa￾tional domain (solid lines, Eq. 13). Curves are smoothed using a rolling average, and values are not shown at t < 200 because the MT rate has not ramped up to its steady-state value. Lcool,tot = Z entire domain E˙ cooldV. (13) Note that there is no cooling inside the Roche lobe of the donor M1 [PITH_… view at source ↗
Figure 5
Figure 5. Figure 5: The gas temperature T in the equatorial plane for all simulations. The snapshots, representative of the simulations once they have reached a quasi-steady state, are labeled with t in both code units and physical units (days). The T floor is reached in the outer regions of the mid mdot simulation and almost everywhere in the low mdot simulation (see text for details). orbit. Unsurprisingly, the luminosity i… view at source ↗
Figure 6
Figure 6. Figure 6: The average gas effective temperature Teff in the accretion disk near M2 and in the circumbinary disk. The curves shown are calculated over cylinders of constant cylindrical radius Rdisk about M2 and constant Rout about the origin (see text for details). Teff values are averaged starting at t = 300 because the MT rate has not reached a steady-state value and the accretion disk has not yet formed at early t… view at source ↗
Figure 7
Figure 7. Figure 7: The specific AM of outflowing material hloss, in units of hL2, versus time in code units, where 1 orbit = 2π. hloss is calculated over spheres centered at the COM, of radius 3, 4, and 4.65 (dashed, dotted, solid lines) and all curves are smoothed using a rolling average. L2 and the COM. In units of √ GMtota, hL2 = 1.56 for the q = 0.5 binaries that we simulate. The specific AM ⃗h of an arbitrary fluid elem… view at source ↗
Figure 8
Figure 8. Figure 8: Properties related to the energy and velocity of outflowing gas. Curves are time-smoothed using a rolling average. Left: The Bernoulii parameter B, in units of GMtot/a, for COM-centered spheres of |⃗r − ⃗rcom| = 3, 4 and rmax ≈ 4.65, in units of a = 1 (dashed, dotted, solid lines)). Right: The volume-averaged radial velocity of the outflow, v ′ r,inertial (Eq. 25), in units of the orbital velocity for each… view at source ↗
Figure 9
Figure 9. Figure 9: The ratio of gas pressure to radiation pressure (see Eq. A1) in the equatorial plane of the simulations, for snapshots representative of the quasi-steady state behavior. both the accretion disk and in the circumbinary disk, meaning that material rapidly cools to near the inter￾nal energy floor (which is set via a temperature floor of 2000 K in the higher density regions such as the accre￾tion disk). C. ADO… view at source ↗
Figure 10
Figure 10. Figure 10: Plots of the opacity used in this work. Colored dots correspond to a sample of grid cells, corresponding to the labeled t, and where they lie in ρ, T space. Cyan points are located in the Roche lobe surrounding M2. Red points are in the equatorial plane, excluding any points falling into M2’s Roche lobe. Orange points are sampled everywhere else in the domain, i.e. above the equatorial plane. The highest … view at source ↗
read the original abstract

High rates of stable mass transfer (MT) occur for some binary star systems, resulting in luminous transients and circumbinary outflows (CBOs). We perform hydrodynamical simulations of a $10 \ M_\odot$ donor star and a $5\ M_\odot$ point mass accretor, incorporating approximate effects of radiative cooling. By varying the orbital separation of the system, we probe MT rates between $10^{-5}$ and $10^{-1} M_\odot$/yr. Mass flows from the donor into a disk around the accretor, with significant equatorially concentrated outflows through the outer Lagrange point L2 occurring for MT rates $\gtrsim 10^{-3} M_\odot$/yr, while the MT remains mostly conservative for lower MT rates. In all cases, any outflowing gas approximately carries the specific angular momentum of L2. The gas cooling luminosity $L$ and temperature increases with MT rate, with $L \sim 10^{5} L_\odot$ and $T \sim 10^{4}$ K for simulations featuring the strongest outflows, with contributions from both the CBO and the accretor's disk. The most luminous transients associated with mass outflows will be rare due to the high MT rate requirement, but generate significant optical emission from both near the accretor and the CBO.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The paper reports hydrodynamical simulations of a 10 M_⊙ donor and 5 M_⊙ accretor binary incorporating approximate radiative cooling. Varying orbital separation probes mass-transfer rates from 10^{-5} to 10^{-1} M_⊙ yr^{-1}. Mass flows into a disk around the accretor; significant equatorially concentrated L2 outflows appear for rates ≳ 10^{-3} M_⊙ yr^{-1} while lower rates remain mostly conservative. Outflowing gas carries approximately the specific angular momentum of L2. Cooling luminosity and temperature rise with MT rate, reaching L ∼ 10^5 L_⊙ and T ∼ 10^4 K in the strongest-outflow cases, with contributions from both the circumbinary outflow and the disk.

Significance. If the central results hold, the work supplies concrete flow structures, outflow fractions, and luminosities for high-rate mass transfer, directly relevant to luminous transients and circumbinary outflows. A clear strength is the direct numerical integration of the hydrodynamical equations with cooling across a controlled range of orbital separations, yielding internally consistent trends in outflow onset and angular-momentum transport without fitted parameters.

major comments (1)
  1. [Numerical setup / cooling implementation] The approximate radiative cooling prescription that governs the thermal evolution of the mass-transfer stream and disk is load-bearing for the reported transition at MT rates ≳ 10^{-3} M_⊙ yr^{-1}. Because outflow dynamics and the fraction of mass reaching L2 depend on local sound speed and temperature, the lack of any resolution study, comparison run with an alternative cooling function, or cross-check against gray or multi-group radiative transfer leaves the quantitative location of the threshold insecure (see abstract and numerical-setup description).
minor comments (1)
  1. [Abstract] The abstract states that 'any outflowing gas approximately carries the specific angular momentum of L2' in all cases; a brief quantitative statement of the measured deviation (e.g., percentage scatter) would strengthen this claim.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful review and constructive feedback on our manuscript. We address the major comment below and describe the revisions we will incorporate.

read point-by-point responses
  1. Referee: The approximate radiative cooling prescription that governs the thermal evolution of the mass-transfer stream and disk is load-bearing for the reported transition at MT rates ≳ 10^{-3} M_⊙ yr^{-1}. Because outflow dynamics and the fraction of mass reaching L2 depend on local sound speed and temperature, the lack of any resolution study, comparison run with an alternative cooling function, or cross-check against gray or multi-group radiative transfer leaves the quantitative location of the threshold insecure (see abstract and numerical-setup description).

    Authors: We agree that the approximate cooling prescription is central to the reported transition and that the absence of dedicated resolution studies or comparisons with alternative cooling functions or full radiative transfer leaves the precise numerical value of the threshold somewhat uncertain. Full radiative transfer calculations remain computationally prohibitive for the broad parameter space of orbital separations and mass-transfer rates we explore. Within our controlled suite, however, the trend is robust: higher mass-transfer rates produce denser, more efficiently cooled streams and disks that drive equatorially concentrated L2 outflows, while lower rates remain largely conservative. The threshold is already phrased as approximate (≳ 10^{-3} M_⊙ yr^{-1}) to reflect the discrete sampling of our models. We have revised the numerical-setup section to include an explicit discussion of the cooling-function limitations, the expected sensitivity of the exact threshold to radiative-transfer details, and the physical robustness of the qualitative outflow onset. These changes clarify the scope of the quantitative claims without altering the core hydrodynamical results. revision: yes

Circularity Check

0 steps flagged

No circularity: results from direct hydrodynamical integration

full rationale

The paper obtains its claims on mass flow structure, L2 outflows at MT rates ≳ 10^{-3} M_⊙/yr, conservative MT at lower rates, and associated luminosities by performing hydrodynamical simulations of a 10 M_⊙ donor and 5 M_⊙ accretor while varying orbital separation to span MT rates from 10^{-5} to 10^{-1} M_⊙/yr. The governing equations are integrated numerically with an approximate radiative cooling term; no central quantity (e.g., outflow angular momentum or transition threshold) is defined in terms of itself, fitted to a subset of the same data and then relabeled as a prediction, or reduced to a load-bearing self-citation whose validity depends on the present work. The derivation chain therefore remains independent of the enumerated circularity patterns and is self-contained within the simulation setup.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the hydrodynamical code accurately evolving the flow under an approximate cooling term and on the point-mass accretor plus donor model capturing the essential Roche-lobe overflow dynamics.

axioms (1)
  • domain assumption Approximate radiative cooling effects are adequate to model the thermal evolution of the gas in the mass-transfer stream and disk.
    Invoked to justify the inclusion of cooling while varying orbital separation to reach different MT rates.

pith-pipeline@v0.9.0 · 5782 in / 1301 out tokens · 47144 ms · 2026-05-21T20:53:37.953420+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Formation of extremely low-mass white dwarf binaries undergoing enhanced angular momentum loss

    astro-ph.SR 2026-05 unverdicted novelty 5.0

    Enhanced AML via L2-point mass loss in the RLOF channel alters ELM WD internal structure and mass-radius relation, reproducing observed shorter orbital periods.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages · cited by 1 Pith paper · 4 internal anchors

  1. [1]

    Podsiadlowski, P., & Rupen, M. P. 2001, ApJ, 562, L79, doi: 10.1086/324573

  2. [2]

    B., & Church, R

    Bobrick, A., Davies, M. B., & Church, R. P. 2017, Mon Not R Astron Soc, 467, 3556, doi: 10.1093/mnras/stx312

  3. [3]

    A., Mohamed, S., & Podsiadlowski, P

    Booth, R. A., Mohamed, S., & Podsiadlowski, P. 2016, Monthly Notices of the Royal Astronomical Society, 457, 822, doi: 10.1093/mnras/stw001

  4. [4]

    Bowler, M. G. 2010, A&A, 521, A81, doi: 10.1051/0004-6361/201014711 —. 2011a, A&A, 531, A107, doi: 10.1051/0004-6361/201016381 —. 2011b, A&A, 534, A112, doi: 10.1051/0004-6361/201117756

  5. [5]

    M., Postnov, K

    Cherepashchuk, A. M., Postnov, K. A., & Belinski, A. A. 2018, Monthly Notices of the Royal Astronomical Society, 479, 4844, doi: 10.1093/mnras/sty1853

  6. [6]

    2020, Monthly Notices of the Royal Astronomical Society, 492, 2208, doi: 10.1093/mnras/stz3598

    Clark, P., Maguire, K., Inserra, C., et al. 2020, Monthly Notices of the Royal Astronomical Society, 492, 2208, doi: 10.1093/mnras/stz3598

  7. [7]

    Davidge, T. J. 2023, The Astronomical Journal, 165, 189, doi: 10.3847/1538-3881/acc580

  8. [8]

    Cuadra and P

    Forgan, D., Rice, K., Stamatellos, D., & Whitworth, A. 2009, Monthly Notices of the Royal Astronomical Society, 394, 882, doi: 10.1111/j.1365-2966.2008.14373.x

  9. [9]

    2021, The Astrophysical Journal, 922, 110, doi: 10.3847/1538-4357/ac2610

    Kalogera, V. 2021, The Astrophysical Journal, 922, 110, doi: 10.3847/1538-4357/ac2610

  10. [10]

    Harris, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2

  11. [11]

    1963, The Astrophysical Journal, 138, 471, doi: 10.1086/147659

    Huang, S.-S. 1963, The Astrophysical Journal, 138, 471, doi: 10.1086/147659

  12. [12]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55

  13. [13]

    S., Bauer, E

    Jermyn, A. S., Bauer, E. B., Schwab, J., et al. 2023, ApJS, 265, 15, doi: 10.3847/1538-4365/acae8d

  14. [14]

    2022, Monthly Notices of the Royal Astronomical Society, 511, 5936, doi: 10.1093/mnras/stac415

    Korol, V., Hallakoun, N., Toonen, S., & Karnesis, N. 2022, Monthly Notices of the Royal Astronomical Society, 511, 5936, doi: 10.1093/mnras/stac415

  15. [15]

    A., & Leahy, J

    Leahy, D. A., & Leahy, J. C. 2015, Computational Astrophysics and Cosmology, 2, 4, doi: 10.1186/s40668-015-0008-8 LIGO Scientific Collaboration, and KAGRA Collaboration, V. C., Abbott, R., Abbott, T., et al. 2023, Phys. Rev. X, 13, 041039, doi: 10.1103/PhysRevX.13.041039

  16. [16]

    C., McInally, W

    Lombardi, Jr, J. C., McInally, W. G., & Faber, J. A. 2015, Monthly Notices of the Royal Astronomical Society, 447, 25, doi: 10.1093/mnras/stu2432

  17. [17]

    2025, Tidal equilibrium of a star in Roche potential, Zenodo, doi: 10.5281/zenodo.15499473

    Lu, W. 2025, Tidal equilibrium of a star in Roche potential, Zenodo, doi: 10.5281/zenodo.15499473

  18. [18]

    Physical model

    Lu, W., Fuller, J., Quataert, E., & Bonnerot, C. 2023, Monthly Notices of the Royal Astronomical Society, 519, 1409, doi: 10.1093/mnras/stac3621

  19. [19]

    C., & Stone, J

    MacLeod, M., Ostriker, E. C., & Stone, J. M. 2018a, The Astrophysical Journal, 863, 5, doi: 10.3847/1538-4357/aacf08 —. 2018b, The Astrophysical Journal, 868, 136, doi: 10.3847/1538-4357/aae9eb

  20. [20]

    Marchant, P., Pappas, K. M. W., Gallegos-Garcia, M., et al. 2021, Astronomy and Astrophysics, 650, A107, doi: 10.1051/0004-6361/202039992

  21. [21]

    2007, The Astrophysical Journal Supplement Series, 170, 228, doi: 10.1086/513316

    Mignone, A., Bodo, G., Massaglia, S., et al. 2007, The Astrophysical Journal Supplement Series, 170, 228, doi: 10.1086/513316

  22. [22]

    2012, ApJS, 198, 7, doi: 10.1088/0067-0049/198/1/7

    Mignone, A., Zanni, C., Tzeferacos, P., et al. 2012, ApJS, 198, 7, doi: 10.1088/0067-0049/198/1/7

  23. [23]

    Mink, S. E. d., Pols, O. R., & Hilditch, R. W. 2007, A&A, 467, 1181, doi: 10.1051/0004-6361:20067007

  24. [24]

    2012, Baltic Astronomy, 21, 88, doi: 10.1515/astro-2017-0362

    Mohamed, S., & Podsiadlowski, P. 2012, Baltic Astronomy, 21, 88, doi: 10.1515/astro-2017-0362

  25. [25]

    V., Glazunova, L

    Nazarenko, V. V., Glazunova, L. V., & Shakun, L. S. 2005, Astronomy Reports, 49, 284, doi: 10.1134/1.1898406

  26. [26]

    2011, The Astrophysical Journal Supplement Series, 192, 3, doi: 10.1088/0067-0049/192/1/3

    Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3, doi: 10.1088/0067-0049/192/1/3

  27. [27]

    2013, The Astrophysical Journal Supplement Series, 208, 4, doi: 10.1088/0067-0049/208/1/4

    Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4, doi: 10.1088/0067-0049/208/1/4

  28. [28]

    2015, The Astrophysical Journal Supplement Series, 220, 15, doi: 10.1088/0067-0049/220/1/15

    Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJS, 220, 15, doi: 10.1088/0067-0049/220/1/15

  29. [29]

    B., et al

    Paxton, B., Schwab, J., Bauer, E. B., et al. 2018, ApJS, 234, 34, doi: 10.3847/1538-4365/aaa5a8

  30. [30]

    2019, The Astrophysical Journal Supplement Series, 243, 10, doi: 10.3847/1538-4365/ab2241

    Paxton, B., Smolec, R., Schwab, J., et al. 2019, ApJS, 243, 10, doi: 10.3847/1538-4365/ab2241

  31. [31]

    2014, The Astrophysical Journal, 788, 22, doi: 10.1088/0004-637X/788/1/22

    Pejcha, O. 2014, The Astrophysical Journal, 788, 22, doi: 10.1088/0004-637X/788/1/22

  32. [32]

    D., & Tomida, K

    Pejcha, O., Metzger, B. D., & Tomida, K. 2016, Monthly Notices of the Royal Astronomical Society, 455, 4351, doi: 10.1093/mnras/stv2592 Perez M., S., & Blundell, K. M. 2010, Mon Not R Astron Soc, 408, 2, doi: 10.1111/j.1365-2966.2010.16638.x

  33. [33]

    2024, Astronomy and Astrophysics, 681, A31, doi: 10.1051/0004-6361/202347090

    Picco, A., Marchant, P., Sana, H., & Nelemans, G. 2024, Astronomy and Astrophysics, 681, A31, doi: 10.1051/0004-6361/202347090

  34. [34]

    C., & Hsu, J

    Podsiadlowski, P., Joss, P. C., & Hsu, J. J. L. 1992, The Astrophysical Journal, 391, 246, doi: 10.1086/171341

  35. [35]

    A., & Yungelson, L

    Postnov, K. A., & Yungelson, L. R. 2014, Living Rev. Relativ., 17, 3, doi: 10.12942/lrr-2014-3 18Scherbak, Lu, and Fuller

  36. [36]

    E., de Koter, A., et al

    Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444, doi: 10.1126/science.1223344

  37. [37]

    2025, Rapid binary mass transfer: Circumbinary outflows and angular momentum losses

    Scherbak, P., Lu, W., & Fuller, J. 2025, Rapid binary mass transfer: Circumbinary outflows and angular momentum losses. https://arxiv.org/abs/2505.21264

  38. [38]

    Schneider, F. R. N., Podsiadlowski, P., & M¨ uller, B. 2021, A&A, 645, A5, doi: 10.1051/0004-6361/202039219

  39. [39]

    R., Schaefer, G

    Shepard, K., Gies, D. R., Schaefer, G. H., et al. 2024, ApJ, 977, 236, doi: 10.3847/1538-4357/ad82e7

  40. [40]

    Stability Criteria for Mass Transfer in Binary Stellar Evolution

    Soberman, G. E., Phinney, E. S., & van den Heuvel, E. P. J. 1997, Stability criteria for mass transfer in binary stellar evolution., arXiv, doi: 10.48550/arXiv.astro-ph/9703016

  41. [41]

    P., Bisbas, T., & Goodwin, S

    Stamatellos, D., Whitworth, A. P., Bisbas, T., & Goodwin, S. 2007, A&A, 475, 37, doi: 10.1051/0004-6361:20077373

  42. [42]

    D., Sollerman , J., et al

    Taddia, F., Stritzinger, M. D., Sollerman, J., et al. 2013, A&A, 555, A10, doi: 10.1051/0004-6361/201321180

  43. [43]

    Oliphant, Matt Haberland, Tyler Reddy, David Cournapeau, Evgeni Burovski, Pearu Peterson, Warren Weckesser, Jonathan Bright, St´ efan J

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nat Methods, 17, 261, doi: 10.1038/s41592-019-0686-2

  44. [44]

    arXiv:1104.3868 [astro-ph.CO] 284 Roedig C, Sesana A, Dotti M, Cuadra J, Amaro-Seoane P, Haardt F (2012) Evolution of binary black holes in self gravitating discs

    Wilkins, D. R., & Clarke, C. J. 2012, Monthly Notices of the Royal Astronomical Society, 419, 3368, doi: 10.1111/j.1365-2966.2011.19976.x

  45. [45]

    C., & Fuller, J

    Wu, S. C., & Fuller, J. 2022, ApJL, 940, L27, doi: 10.3847/2041-8213/ac9b3d