An official website of the United States government
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
Dynamical friction modelling of massive black holes in cosmological simulations and effects on merger rate predictions
ABSTRACT In this work, we establish and test methods for implementing dynamical friction (DF) for massive black hole pairs that form in large volume cosmological hydrodynamical simulations that include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of DF from collisionless particles, black holes can effectively sink close to the galaxy centre, provided that the black hole’s dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays a role in assisting the black holes’ orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole’s evolution. DF from gas becomes less than $$1{{\ \rm per\ cent}}$$ of DF from collisionless particles for BH masses >107 M⊙. Using our best DF model, we calculate the merger rate down to z = 1.1 using an Lbox = 35 Mpc h−1 simulation box. We predict ∼2 mergers per year for z > 1.1 peaking at z ∼ 2. These merger rates are within the range obtained in previous work using similar resolution hydrodynamical simulations. We show that the rate is enhanced by factor of ∼2 when DF is taken into account in the simulations compared to the no-DF run. This is due to $${\gt}40{{\ \rm per\ cent}}$$ more black holes reaching the centre of their host halo when DF is added.
more »
« less
- PAR ID:
- 10360860
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 510
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 531-550
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disc galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of $$z$$ ≳ 1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modelled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyse the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a range of galaxy properties from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M⋆ > 1011 M⊙ at $$z$$ = 1). Selecting only snapshots where a disc is present, we calculate the rotational profile $$\bar{v}_\phi (r)$$ of the cool ($$10^{3.5}\,\lt {\it T}\lt 10^{4.5}~\rm {K}$$) gas and compare it to the circular velocity $$v_{\rm c}=\sqrt{GM_{\rm enc}/r}$$. In the simulated galaxies, the gas rotation traces the circular velocity at intermediate radii, but the two quantities diverge significantly in the centre and in the outer disc. Our simulations appear to over-predict observed rotational velocities in the centres of massive galaxies (likely from a lack of black hole feedback), so we focus on larger radii. Gradients in the turbulent pressure at these radii can provide additional radial support and bias dynamical mass measurements low by up to 40 per cent. In both the interior and exterior, the gas’ motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly used analytic models for pressure gradients (or ‘asymmetric drift’) in the ISM of high-redshift galaxies.more » « less
-
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
-
ABSTRACT Accretion discs around supermassive black holes are promising sites for stellar mass black hole mergers detectable with LIGO. Here we present the results of Monte Carlo simulations of black hole mergers within 1-d AGN disc models. For the spin distribution in the disc bulk, key findings are: (1) The distribution of χeff is naturally centred around $$\tilde{\chi }_{\rm eff} \approx 0.0$$, (2) the width of the χeff distribution is narrow for low natal spins. For the mass distribution in the disc bulk, key findings are: (3) mass ratios $$\tilde{q} \sim 0.5\!-\!0.7$$, (4) the maximum merger mass in the bulk is $$\sim 100\!-\!200\, \mathrm{M}_{\odot }$$, (5) $$\sim 1{{\ \rm per\ cent}}$$ of bulk mergers involve BH $$\gt 50\, \mathrm{M}_{\odot }$$ with (6) $$\simeq 80{{\ \rm per\ cent}}$$ of bulk mergers are pairs of first generation BH. Additionally, mergers at a migration trap grow an IMBH with typical merger mass ratios $$\tilde{q}\sim 0.1$$. Ongoing LIGO non-detections of black holes $$\gt 10^{2}\, \mathrm{M}_{\odot }$$ puts strong limits on the presence of migration traps in AGN discs (and therefore AGN disc density and structure) as well as median AGN disc lifetime. The highest merger rate occurs for this channel if AGN discs are relatively short-lived (≤1 Myr) so multiple AGN episodes can happen per Galactic nucleus in a Hubble time.more » « less
-
ABSTRACT JWST has revealed a large population of accreting black holes (BHs) in the early Universe. Recent work has shown that even after accounting for possible systematic biases, the high-z$$M_*{\!-\!}M_{\rm \rm bh}$$ relation can be above the local scaling relation by $$\gt 3\sigma$$. To understand the implications of these overmassive high-z BHs, we study the BH growth at $$z\sim 4{\!-\!}7$$ using the $$[18~\mathrm{Mpc}]^3$$BRAHMA cosmological simulations with systematic variations of heavy seed models that emulate direct collapse black hole (DCBH) formation. In our least restrictive seed model, we place $$\sim 10^5~{\rm M}_{\odot }$$ seeds in haloes with sufficient dense and metal-poor gas. To model conditions for direct collapse, we impose additional criteria based on a minimum Lyman Werner flux (LW flux $$=10~J_{21}$$), maximum gas spin, and an environmental richness criterion. The high-z BH growth in our simulations is merger dominated, with a relatively small contribution from gas accretion. The simulation that includes all the above seeding criteria fails to reproduce an overmassive high-z$$M_*{\!-\!}M_{\rm bh}$$ relation consistent with observations (by factor of $$\sim 10$$ at $$z\sim 4$$). However, more optimistic models that exclude the spin and environment based criteria are able to reproduce the observed relations if we assume $$\lesssim 750~\mathrm{Myr}$$ delay times between host galaxy mergers and subsequent BH mergers. Overall, our results suggest that current JWST observations may be explained with heavy seeding channels if their formation is more efficient than currently assumed DCBH conditions. Alternatively, we may need higher initial seed masses, additional contributions from lighter seeds to BH mergers, and / or more efficient modes for BH accretion.more » « less
