The launch of spacebased gravitationalwave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays will extend the GW window to low frequencies, opening new investigations into dynamical processes involving massive black hole binaries (MBHBs) and their mergers across cosmic time. MBHBs are expected to be among the primary sources for the upcoming lowfrequency (10−4–10−1 Hz) window probed by LISA. It is important to investigate the expected supermassive BH merger rates and associated signals, to determine how potential LISA events are affected by physics included in current models. To study this, we postprocess the large population of MBHBs in the Illustris simulation to account for dynamical friction time delays associated with BH infall/inspiral. We show that merger delays associated with binary evolution have the potential to decrease the expected merger rates, with $M_{\rm {BH}}\ \gt\ 10^6\ \mathrm{M}_\odot$ MBHBs (the lowest mass in Illustris) decreasing from ∼3 to ∼0.1 yr−1, and shifting the merger peak from z ∼2 to ∼1.25. During this time, we estimate that accretion grows the total merging mass by as much as 7x the original mass. Importantly, however, dynamical frictionassociated delays (which shift the mergers toward lower redshift and higher masses) lead to a stronger signal/strain for the emitted GWs in the LISA band, increasing mean frequency from 10−3.1 to 10−3.4–10−4.0 Hz, and mean strain from 10−17.2 to 10−16.3–10−15.3. Finally, we show that after including a merger delay and associated MBH growth, mergers still tend to lie on the typical MBH–M* relation, but with an increased likelihood of an undermassive BH.
more » « less NSFPAR ID:
 10406817
 Publisher / Repository:
 Oxford University Press
 Date Published:
 Journal Name:
 Monthly Notices of the Royal Astronomical Society
 Volume:
 512
 Issue:
 4
 ISSN:
 00358711
 Page Range / eLocation ID:
 p. 60076020
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

null (Ed.)ABSTRACT Massive black hole (MBH) binary inspiral timescales 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 subresolution 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 unequalmass 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

ABSTRACT Massive black holes in the centres of galaxies today must have grown by several orders of magnitude from seed black holes formed at early times. Detecting a population of intermediate mass black holes (IMBHs) can provide constraints on these elusive BH seeds. Here, we use the large volume cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and dynamical friction to investigate the population of IMBH seeds. Dynamical friction is largely inefficient at sinking and merging seed IMBHs at highz. This leads to an extensive population (several hundred per galaxy) of wandering IMBHs in large haloes at $z\sim 2$. A small fraction of these IMBHs are detectable as HLXs, Hyper Luminous Xray sources. Importantly, at $z\sim 2$, IMBHs mergers produce the peak of GW events. We find close to a million GW events in Astrid between $z=\rm{2\!\!3}$ involving seed IMBH mergers. These GW events (almost all detectable by LISA) at cosmic noon should provide strong constraints on IMBH seed models and their formation mechanisms. At the centre of massive galaxies, where the number of IMBHs can be as high as 10–100, SMBHIMBH pairs can form. These Intermediate mass ratio inspirals (IMRIs) and extreme mass ratio inspirals (EMRIs), will require the next generation of milli$\mu$Hz spacebased GW interferometers to be detected. Large populations of IMBHs around massive black holes will probe their environments and MBH causal structure.

ABSTRACT We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the largevolume 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 subgrid 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 kpcscales, 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 insimulation 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 seedmass 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 seedmass 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}}$ seedseed mergers, $\sim 30\,{\rm{per\,cent}}$ involving only one seedmass MBH, and $\sim 10\,{\rm{per\,cent}}$ mergers of nonseed MBHs.

Abstract We evaluate the cosmological coalescence and detection rates for massive black hole (MBH) binaries targeted by the gravitational wave observatory Laser Interferometer Space Antenna (LISA). Our calculation starts with a population of gravitationally unbound MBH pairs, drawn from the TNG503 cosmological simulation, and follows their orbital evolution from kiloparsec scales all the way to coalescence using a semianalytic model developed in our previous work. We find that for the majority of MBH pairs that coalesce within a Hubble time dynamical friction is the most important mechanism that determines their coalescence rate. Our model predicts an MBH coalescence rate ≲0.45 yr −1 and a LISA detection rate ≲0.34 yr −1 . Most LISA detections should originate from 10 6 to 10 6.8 M ⊙ MBHs in gasrich galaxies at redshifts 1.6 ≤ z ≤ 2.4 and have a characteristic signaltonoise ratio S/N ∼100. We however find a dramatic reduction in the coalescence and detection rates, as well as the average S/N, if the effects of radiative feedback from accreting MBHs are taken into account. In this case, the MBH coalescence rate is reduced by 78% (to ≲0.1 yr −1 ), and the LISA detection rate is reduced by 94% (to 0.02 yr −1 ), whereas the average S/N is ∼10. We emphasize that our model provides a conservative estimate of the LISA detection rates, due to the limited MBH mass range in TNG503, consistent with other works in the literature that draw their MBH pairs from cosmological simulations.more » « less

null (Ed.)Abstract Among the potential milliHz gravitational wave (GW) sources for the upcoming spacebased interferometer LISA are extreme or intermediatemass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range 105M⊙ ≲ M ≲ 107M⊙ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGN), where torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform. Here we use a suite of twodimensional hydrodynamical simulations with the movingmesh code DISCO to present a systematic study of disc torques. We measure torques on an inspiraling BH and compute the corresponding waveform deviations as a function of the binary mass ratio q ≡ M2/M1, the disc viscosity (α), and gas temperature (or equivalently Mach number; $\mathcal {M}$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends nontrivially on the combination of these parameters. The gas imprint is detectable by LISA for binaries embedded in AGN discs with surface densities above $\Sigma _0\ge 10^{46} \rm \, g cm^{2}$, depending on q, α and $\mathcal {M}$. Deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GWdriven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of discsatellite interaction.more » « less