First Light And Reionization Epoch Simulations (FLARES) -- XIX: Supermassive black hole mergers in the early Universe and their environmental dependence

TL;DR

This study uses cosmological simulations to analyze supermassive black hole mergers at high redshifts, predicting detection rates for LISA and emphasizing the environmental dependence of these events.

Contribution

It provides new insights into SMBH merger rates, mass ratios, and environmental effects at early cosmic times using the FLARES simulation suite.

Findings

Merger rates increase at lower redshifts.

Denser regions host more massive and earlier mergers.

LISA can detect SMBH mergers with specific mass and mass ratio ranges.

Abstract

The upcoming space-based gravitational wave (GW) observatory, LISA, is expected to detect GW signals from supermassive black hole (SMBH) mergers occurring at high redshifts. However, understanding the origin and growth of SMBHs in the early Universe remains an open problem in astrophysics. In this work, we utilize the First Light And Reionization Epoch Simulations (FLARES), a suite of cosmological hydrodynamical zoom-in simulations, to study SMBH mergers at $5 \lesssim z \lesssim 10$ across a wide range of environments. Most mergers in FLARES involve secondary SMBHs near the seed mass ($m_{\rm seed} \approx 1.5 \times 10^{5}~{\rm M}_{\odot}$) while primary SMBHs span up to $10^{9}~{\rm M}_{\odot}$, resulting in mass ratios from $q \sim 10^{-4}$ to $1$, with a peak at $q \sim 1$. The number of mergers increases rapidly towards lower redshifts, and the comoving total number density scales with overdensity as $n_{\rm merger} = 10^{-3.81} (1 + \delta)^{4.78}$. Denser regions host more massive mergers, with higher merger redshifts and lower mass ratios. Within the FLARES redshift range, LISA is expected to detect mergers with $10^{5} \lesssim M_{\rm tot}/{\rm M}_{\odot} \lesssim 10^{8}$ and $q \gtrsim 10^{-2}$, corresponding to a detection rate of $0.030~{\rm yr}^{-1}$ for events with signal-to-noise ratio ${\rm SNR} \geq 10$. Our study demonstrates the sensitivity of GW predictions at high redshifts to SMBH seed models and merger time delays, highlighting the need for improved modeling in future cosmological simulations to maximize LISA's scientific return.