skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: MBH binary intruders: triple systems from cosmological simulations
ABSTRACT Massive black hole (MBH) binaries can form following a galaxy merger, but this may not always lead to a MBH binary merger within a Hubble time. The merger time-scale depends on how efficiently the MBHs lose orbital energy to the gas and stellar background, and to gravitational waves (GWs). In systems where these mechanisms are inefficient, the binary inspiral time can be long enough for a subsequent galaxy merger to bring a third MBH into the system. In this work, we identify and characterize the population of triple MBH systems in the Illustris cosmological hydrodynamic simulation. We find a substantial occurrence rate of triple MBH systems: in our fiducial model, 22 per cent of all binary systems form triples, and $$\gt 70{{\ \rm per\ cent}}$$ of these involve binaries that would not otherwise merge by z = 0. Furthermore, a significant subset of triples (6 per cent of all binaries, or more than a quarter of all triples) form a triple system at parsec scales, where the three BHs are most likely to undergo a strong three-body interaction. Crucially, we find that the rate of triple occurrence has only a weak dependence on key parameters of the binary inspiral model (binary eccentricity and stellar loss-cone refilling rate). We also do not observe strong trends in the host galaxy properties for binary versus triple MBH populations. Our results demonstrate the potential for triple systems to increase MBH merger rates, thereby enhancing the low-frequency GW signals detectable with pulsar timing arrays and with LISA.  more » « less
Award ID(s):
1909933
PAR ID:
10479241
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
527
Issue:
3
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 7424-7437
Size(s):
p. 7424-7437
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT Massive black hole (MBH) binary inspiral time-scales are uncertain, and their spins are even more poorly constrained. Spin misalignment introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with Laser Interferometer Space Antenna, and their retention rate in galaxies. Here, we introduce a sub-resolution model for gas- and gravitational wave (GW)-driven MBH binary spin evolution using accreting MBHs from the Illustris cosmological hydrodynamic simulations. We also model binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes that the circumbinary disc always removes angular momentum from the binary. It also assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47 per cent of the MBHs in our population merge by z = 0. Of these, 19 per cent have misaligned primaries and 10 per cent have misaligned secondaries at the time of merger in our fiducial model with initial eccentricity of 0.6 and accretion rates from Illustris. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor of 100, in a thicker disc, yields 79 and 42 per cent misalignment for primaries and secondaries, respectively. Even in the more conservative fiducial model, more than 12 per cent of binaries experience recoils of >500 km s−1, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession. 
    more » « less
  2. ABSTRACT We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid . Astrid includes galaxy formation and black hole models recently updated with an MBH seed population between 3 × 104h−1M⊙ and 3 × 105h−1M⊙ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to 1.5 ckpc h−1. We calculate the initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above 0.7 for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below $${\sim 200\, \text{pc}}$$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($$\sim 500\, \text{Myr}$$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries (MBH < 2Mseed). As a result, only $$\lesssim 20{{\rm per \,cent}}$$ of seed MBH pairs merge at z > 3 after considering both unresolved DF evolution and binary hardening. These z > 3 seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $$\gt 10^9\, {\rm M}_\odot$$. With the higher initial eccentricity prediction from Astrid , we estimate an expected merger rate of 0.3−0.7 per year from the z > 3 MBH population. This is a factor of ∼7 higher than the prediction using the circular orbit assumption. The Laser Interferometer Space Antenna events are expected at a similar rate, and comprise $$\gtrsim 60\,{\rm{per\,cent}}$$ seed-seed mergers, $$\sim 30\,{\rm{per\,cent}}$$ involving only one seed-mass MBH, and $$\sim 10\,{\rm{per\,cent}}$$ mergers of non-seed MBHs. 
    more » « less
  3. ABSTRACT The fraction of stars that are in binaries or triples at the time of stellar death and the fraction of these systems that survive the supernova explosion are crucial constraints for evolution models and predictions for gravitational wave source populations. These fractions are also subject to direct observational determination. Here, we search 10 supernova remnants containing compact objects with proper motions for unbound binaries or triples using Gaia EDR3 and new statistical methods and tests for false positives. We confirm the one known example of an unbound binary, HD 37424 in G180.0−01.7, and find no other examples. Combining this with our previous searches for bound and unbound binaries, and assuming no bias in favour of finding interacting binaries, we find that 72.0 per cent (52.2–86.4 per cent, 90 per cent confidence) of supernova producing neutron stars are not binaries at the time of explosion, 13.9 per cent (5.4–27.2 per cent) produce bound binaries, and 12.5 per cent (2.8–31.3 per cent) produce unbound binaries. With a strong bias in favour of finding interacting binaries, the medians shift to 76.0 per cent were not binaries at death, 9.5 per cent leave bound binaries, and 13.2 per cent leave unbound binaries. Of explosions that do not leave binaries, $${\lt}18.9{{\ \rm per\ cent}}$$ can be fully unbound triples. These limits are conservatively for $$M\gt 5\, \mathrm{M}_\odot$$ companions, although the mass limits for some individual systems are significantly stronger. At birth, the progenitor of PSR J0538+2817 was probably a 13–$$19\, \mathrm{M}_\odot$$ star, and at the time of explosion, it was probably a Roche limited, partially stripped star transferring mass to HD 37424 and then producing a Type IIL or IIb supernova. 
    more » « less
  4. Abstract The Gaia mission has detected many white dwarfs (WDs) in binary and triple configurations, and while observations suggest that triple-stellar systems are common in our Galaxy, not much attention was devoted to WDs in triples. For stability reasons, these triples must have hierarchical configurations, i.e., two stars are on a tight orbit (the inner binary), with the third companion on a wider orbit about the inner binary. In such a system, the two orbits torque each other via the eccentric Kozai–Lidov mechanism, which can alter the orbital configuration of the inner binary. We simulate thousands of triple-stellar systems for over 10 Gyr, tracking gravitational interactions, tides, general relativity, and stellar evolution up to their WD fate. As demonstrated here, three-body dynamics coupled with stellar evolution is a critical channel to form tight WD binaries or merge a WD binary. Among these triples, we explore their manifestations as cataclysmic variables, Type Ia supernovae, and gravitational-wave events. The simulated systems are then compared to a sample of WD triples selected from the Gaia catalog. We find that including the effect of mass-loss-induced kicks is crucial for producing a distribution of the inner binary–tertiary separations that is consistent with Gaia observations. Lastly, we leverage this consistency to estimate that, at minimum, 30% of solar-type stars in the local 200 pc were born in triples. 
    more » « less
  5. Abstract More than half of all main-sequence (MS) stars have one or more companions, and many of those with initial masses <8Mare born in hierarchical triples. These systems feature two stars in a close orbit (the inner binary) while a tertiary star orbits them on a wider orbit (the outer binary). In hierarchical triples, three-body dynamics combined with stellar evolution drives interactions and, in many cases, merges the inner binary entirely to create a renovated “post-merger binary” (PMB). By leveraging dynamical simulations and tracking binary interactions, we explore the outcomes of merged triples and investigate whether PMBs preserve signatures of their three-body history. Our findings indicate that in 26%–54% of wide double white dwarf (DWD) binaries (s≳ 100 au), the more massive white dwarf (WD) is a merger product, implying that these DWD binaries were previously triples. Overall, we estimate that 44% ± 14% of observed wide DWDs originated in triple star systems and thereby have rich dynamical histories. We also examine MS+MS and MS+red giant mergers manifesting as blue straggler stars (BSSs). These PMBs have orbital configurations and ages similar to most observed BSS binaries. While the triple+merger formation channel can explain the observed chemical abundances, moderate eccentricities, and companion masses in BSS binaries, it likely only accounts for ∼20%–25% of BSSs. Meanwhile, we predict that the majority of observed single BSSs formed as collisions in triples and harbor long-period (>10 yr) companions. Furthermore, both BSS binaries and DWDs exhibit signatures of WD birth kicks. 
    more » « less