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.
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- NSF-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
-
more » « less
ABSTRACT -
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
-
more » « less
ABSTRACT In the near future, projects like Laser Interferometer Space Antenna (LISA) and pulsar timing arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here, we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological simulation, which uses a range of seed masses to probe down to low-mass black holes (BHs), and directly incorporates dynamical friction so as to accurately model the dynamical processes that bring black holes to the galaxy centre where binary formation and coalescence will occur. The black hole populations in Astrid include black holes down to $\sim 10^{4.5} \, \mathrm{M}_\odot$, and remain broadly consistent with the TNG simulations at scales $\gt 10^6 \, \mathrm{M}_\odot$ (the seed mass used in TNG). By resolving lower mass black holes, the overall merger rate is ∼5× higher than in TNG. However, incorporating dynamical friction delays mergers compared to a recentring scheme, reducing the high-z merger rate mass-matched mergers by a factor of ∼2×. We also calculate the expected LISA signal-to-noise values, and show that the distribution peaks at high SNR (>100), emphasizing the importance of implementing a seed mass well below LISA’s peak sensitivity ($\sim 10^6 \, \mathrm{M}_\odot$) to resolve the majority of LISA’s gravitational wave detections.
-
more » « less
ABSTRACT We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\lesssim$4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $\gt t_{\rm H}^{-1}$ to $\lt 0.2t_{\rm H}^{-1}$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\lesssim 1{{\ \rm per\ cent}}$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the H2 content, but for $M_*\gtrsim 10^{10.5} \,\mathrm{ M}_{\odot }$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift, ∝(1 + z)3.5, but the quenching rate evolves much more slowly, ∝(1 + z)0.9; there are insufficient mergers to explain the quenching rate at $z\lesssim 1.5$. simba first quenches galaxies at $z\gtrsim 3$, with a number density in good agreement with observations. The quenching time-scales τq are strongly bimodal, with ‘slow’ quenchings (τq ∼ 0.1tH) dominating overall, but ‘fast’ quenchings (τq ∼ 0.01tH) dominating in M* ∼ 1010–1010.5 M$\odot$ galaxies, likely induced by simba’s jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.
-
more » « less
ABSTRACT We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple massive black hole (MBH) systems at z = 2–3. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find 4 per cent of the $M_{\rm BH}\gt 10^7\, M_\odot$ are in tuples with $\Delta r_{\rm max} \lt 200\, {\rm kpc}$. The tuple systems span a range of separations with the majority of the observable AGN systems at Δr ∼ 50–100 kpc. They include some of the most massive BHs (up to $10^{10} \, M_\odot$) but with at least one of the components of $M_{\rm BH} \sim 10^7 \, {\rm M}_{\odot }$. Tuples’ host galaxies are typically massive with $M_* \sim 10^{10-11} \, M_\odot$. We find that $\gt 10~{{\ \rm per\ cent}}$ massive haloes with Mhalo > 1013 M⊙ host MBH tuples. Following the subsequent interactions between MBHs in tuples, we found that in $\sim 5~{{\ \rm per\ cent}}$ of the triplets all three MBHs merge within a Gyr, and 15 per cent go through one merger. As a by-product of the complex multigalaxy interaction of these systems, we also find that up to $\sim 5~{{\ \rm per\ cent}}$ of tuples lead to runaway MBHs. In ASTRID, virtually all of the ultramassive black holes ($\gt 10^{10} \, M_\odot$) have undergone a triple quasar phase, while for BHs with $M_{\rm BH} \sim 10^9 \, M_\odot$, this fraction drops to 50 per cent.
