skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Orbital evolution of binaries in circumbinary discs

We present the to-date largest parameter space exploration of binaries in circumbinary discs (CBDs), deriving orbital evolution prescriptions for eccentric, unequal mass binaries from our suite of hydrodynamic simulations. In all cases, binary eccentricities evolve towards steady state values that increase with mass ratio, and saturate at an equilibrium eccentricity eb,eq ∼ 0.5 in the large mass ratio regime, in line with resonant theory. For binaries accreting at their combined Eddington limit, a steady state eccentricity can be achieved within a few megayears. Once at their steady state eccentricities, binaries with qb ≳ 0.3 evolve towards coalescence, while lower mass ratio systems expand due to CBD torques. We discuss implications for population studies of massive black hole binaries, protostars in binary systems, and post-common envelope binaries observed by ground-based gravitational wave detectors.

more » « less
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Medium: X Size: p. 2707-2717
["p. 2707-2717"]
Sponsoring Org:
National Science Foundation
More Like this

    We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings. Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker–Planck limit of the Boltzmann equation. Our model quantifies the dominant gravitational physics, combining both strong and weak single–binary interactions, which drives the time evolution of binary orbital parameter distributions in dense stellar environments. We focus in particular the distributions of binary orbital energies and eccentricities. We find a novel steady-state distribution of binary eccentricities, featuring strong depletions of both the highest and the lowest eccentricity binaries. In energy space, we compare the predictions of our analytic model to the results of numerical N-body simulations, and find that the agreement is good for the initial conditions considered here. This work is a first step towards the development of a fully self-consistent semi-analytic model for dynamically evolving binary star populations in dense stellar environments due to direct few-body interactions.

    more » « less

    Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs are possible progenitors of black hole mergers detected in gravitational waves by LIGO/VIRGO. To better understand the hydrodynamical evolution of BBHs interacting with the disc gas, we perform a suite of high-resolution 2D simulations of binaries in local disc (shearing-box) models, considering various binary mass ratios, eccentricities and background disc properties. We use the γ-law equation of state and adopt a robust post-processing treatment to evaluate the mass accretion rate, torque and energy transfer rate on the binary to determine its long-term orbital evolution. We find that circular comparable-mass binaries contract, with an orbital decay rate of a few times the mass doubling rate. Eccentric binaries always experience eccentricity damping. Prograde binaries with higher eccentricities or smaller mass ratios generally have slower orbital decay rates, with some extreme cases exhibiting orbital expansion. The averaged binary mass accretion rate depends on the physical size of the accretor. The accretion flows are highly variable, and the dominant variability frequency is the apparent binary orbital frequency (in the rotating frame around the central massive BH) for circular binaries but gradually shifts to the radial epicyclic frequency as the binary eccentricity increases. Our findings demonstrate that the dynamics of BBHs embedded in AGN discs is quite different from that of isolated binaries in their own circumbinary discs. Furthermore, our results suggest that the hardening time-scales of the binaries are much shorter than their migration time-scales in the disc, for all reasonable binary and disc parameters.

    more » « less
  3. We report discovery and characterization of a main-sequence G star orbiting a dark object with mass1.90±0.04M. The system was discovered via Gaia astrometry and has an orbital period of 731 days. We obtained multi-epoch RV follow-up over a period of 639 days, allowing us to refine the Gaia orbital solution and precisely constrain the masses of both components. The luminous star is a12,Gyr-old, low-metallicity halo star near the main-sequence turnoff (,K; ; ;M0.79M) with a highly enhanced lithium abundance. The RV mass function sets a minimum companion mass for an edge-on orbit ofM2>1.67M, well above the Chandrasekhar limit. The Gaia inclination constraint,i=68.7±1.4,deg, then implies a companion mass ofM2=1.90±0.04M. The companion is most likely a massive neutron star: the only viable alternative is two massive white dwarfs in a close binary, but this scenario is disfavored on evolutionary grounds. The system’s low eccentricity (e=0.122±0.002) disfavors dynamical formation channels and implies that the neutron star likely formed with little mass loss (1M) and with a weak natal kick (). Stronger kicks with more mass loss are not fully ruled out but would imply that a larger population of similar systems with higher eccentricities should exist. The current orbit is too small to have accommodated the neutron star progenitor as a red supergiant or super-AGB star. The simplest formation scenario – isolated binary evolution – requires the system to have survived unstable mass transfer and common envelope evolution with a donor-to-accretor mass ratio>10. The system, which we call Gaia NS1, is likely a progenitor of symbiotic X-ray binaries and long-period millisecond pulsars. Its discovery challenges binary evolution models and bodes well for Gaia’s census of compact objects in wide binaries.

    more » « less

    Massive Black Hole (MBH) binaries are considered to be one of the most important sources of Gravitational Waves (GW) that can be detected by GW detectors like LISA. However, there are a lot of uncertainties in the dynamics of MBH binaries in the stages leading up to the GW-emission phase. It has been recently suggested that Nuclear Star Clusters (NSCs) could provide a viable route to overcome the final parsec problem for MBH binaries at the centre of galaxies. NSCs are collisional systems where the dynamics would be altered by the presence of a mass spectrum. In this study, we use a suite of high-resolution N-body simulations with over 1 million particles to understand how collisional relaxation under the presence of a mass spectrum of NSC particles affects the dynamics of the MBH binary under the merger of two NSCs. We consider MBH binaries with different mass ratios and additional non-relaxed models. We find that mass-segregation driven by collisional relaxation can lead to accelerated hardening in lower mass ratio binaries but has the opposite effect in higher mass ratio binaries. Crucially, the relaxed models also demonstrate much lower eccentricities at binary formation and negligible growth during hardening stages leading to longer merger time-scales. The results are robust and highlight the importance of collisional relaxation on changing the dynamics of the binary. Our models are state-of-the-art, use zero softening, and high enough particle numbers to model NSCs realistically.

    more » « less

    Detached eclipsing binaries are a fundamental tool for measuring the physical parameters of stars that are effectively evolving in isolation. Starting from more than 40 000 eclipsing binary candidates identified by the All-Sky Automated Survey for Supernovae (ASAS-SN), we use PHOEBE to determine the sum of the fractional radii, the ratio of effective temperatures, the inclinations, and the eccentricities for 35 576 systems. We visually inspect all the light-curve models to verify the model fits and examine the TESS light curves, when available, to select systems with evidence for additional physics, such as spots, mass transfer, and hierarchical triples. We examine the distributions of the eclipsing binary model parameters and the orbital parameters. We identify two groups in the sum of the fractional radii and effective temperature ratio parameter space that may distinguish systems approaching the semidetached limit. Combining Gaia EDR3 with extinction estimates from three-dimensional dust maps, we examine the properties of the systems as a function of their absolute magnitude and evolutionary state. Finally, we present light curves of selected eclipsing binaries that may be of interest for follow-up studies.

    more » « less