An official website of the United States government
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
Massive black hole binary inspiral and spin evolution in a cosmological framework
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
- PAR ID:
- 10285763
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 501
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2531 to 2546
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We assess the possibility of detecting both eccentricity and gas effects (migration and accretion) in the gravitational wave (GW) signal from LISA massive black hole binaries at redshift $z=1$. Gas induces a phase correction to the GW signal with an effective amplitude ($$C_{\rm g}$$) and a semimajor axis dependence (assumed to follow a power-law with slope $$n_{\rm g}$$). We use a complete model of the LISA response and employ a gas-corrected post-Newtonian inspiral-only waveform model TaylorF2Ecc. By using the Fisher formalism and Bayesian inference, we constrain $$C_{\rm g}$$ together with the initial eccentricity $$e_0$$, the total redshifted mass $$M_z$$, the primary-to-secondary mass ratio q, the dimensionless spins $$\chi _{1,2}$$ of both component BHs, and the time of coalescence $$t_c$$. We find that simultaneously constraining $$C_{\rm g}$$ and $$e_0$$ leads to worse constraints on both parameters with respect to when considered individually. For a standard thin viscous accretion disc around $$M_z=10^5~{\rm M}_{\odot }$$, $q=8$, $$\chi _{1,2}=0.9$$, and $$t_c=4$$ years MBHB, we can confidently measure (with a relative error of $$\lt 50$$ per cent) an Eddington ratio $${\rm f}_{\rm Edd}\sim 0.1$$ for a circular binary and $${\rm f}_{\rm Edd}\sim 1$$ for an eccentric system assuming $$\mathcal {O}(10)$$ stronger gas torque near-merger than at the currently explored much-wider binary separations. The minimum measurable eccentricity is $$e_0\gtrsim 10^{-2.75}$$ in vacuum and $$e_0\gtrsim 10^{-2}$$ in gas. A weak environmental perturbation ($${\rm f}_{\rm Edd}\lesssim 1$$) to a circular binary can be mimicked by an orbital eccentricity during inspiral, implying that an electromagnetic counterpart would be required to confirm the presence of an accretion disc.more » « less
-
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
-
ABSTRACT Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrounding geometrically thin disc on this precession, as would be the case in the truncated disc picture of the hard state. To address this, we perform a set of general relativistic magnetohydrodynamics simulations of truncated discs misaligned with the spin axes of their central black holes. Our results confirm that the inner-hot flow still undergoes precession, though at a rate that is only 5 per cent of what is predicted for an isolated precessing torus. We find that the exchange of angular momentum between the outer thin and the inner thick disc causes this slow-down in the precession rate and discuss its relevance to type-C QPOs.more » « less
-
Binary black holes with misaligned spins will generically induce both precession and nutation of the orbital angular momentum 𝐋 about the total angular momentum 𝐉. These phenomena modulate the phase and amplitude of the gravitational waves emitted as the binary inspirals to merger. We introduce a “taxonomy” of binary black-hole spin precession that encompasses all the known phenomenology, then present five new phenomenological parameters that describe generic precession and constitute potential building blocks for future gravitational waveform models. These are the precession amplitude ⟨𝜃𝐿⟩, the precession frequency ⟨Ω𝐿⟩, the nutation amplitude Δ𝜃𝐿, the nutation frequency 𝜔, and the precession-frequency variation ΔΩ𝐿. We investigate the evolution of these five parameters during the inspiral and explore their statistical properties for sources with isotropic spins. In particular, we find that nutation of 𝐋 is most prominent for binaries with high spins (𝜒≳0.5) and moderate mass ratios (𝑞∼0.6).more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa3826
