Supermassive Black Hole Mergers Research Articles

Noise test of an electromagnetic antenna for low frequency scattered gravitational waves

Abstract Gravitational waves (GW) proved hard to detect experimentally and yet carry high energy fluxes. Conversely the technology for electromagnetic waves (EMW) is sophisticated, easily available and highly sensitive. Various proposals have been made to use the conversion of a gravitational wave into an electromagnetic wave to generate signals that could be detected. We present test results from a novel electromagnetic antenna, efficient at the very low frequencies needed to explore the signals from supermassive black holes mergers. These will eventually be studied by the Laser Interferometric Space Antenna (LISA) at wavelengths between 108 and 1011 m. Common forms of electromagnetic antennae are only efficient when the ratio of wavelength, λ, to antenna length, L, is roughly between 0.1 and 10. This paper proposes a scheme for a realistic ground-based electromagnetic antenna , ”the resistive antenna”, which operates efficiently at a ratio $\frac{\lambda }{L}$ of more than 106. The first estimates of ambient noise at these low frequencies from field trials of such a device are presented. The measured levels of ambient noise result in this detection scheme having a sensitivity much worse than LISA. The antenna design described in this communication may also have application in other fields.

Dual Active Galactic Nuclei: Precursors of Binary Supermassive Black Hole Formation and Mergers

The presence of dual active galactic nuclei (AGNs) on scales of a few tens of kiloparsecs can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. This study uses the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs (L bol > 1043 erg s−1) and their neighboring SMBHs (any SMBHs closer than 30 pkpc to an AGN) at z ≤ 2. This is our underlying population. We then applied the luminosity threshold of L bol > 1043 erg s−1 to the neighboring SMBHs thereby identifying dual and multiple AGNs. Our findings indicate an increase in the number of both single and dual AGNs from lower to higher redshifts. We also find that the number of dual AGNs with separations of 0.5–4 kpc is twice the number of duals with separations of 4–30 kpc. All dual AGNs in our sample resulted from major mergers. Compared to single AGNs, duals have a lower black hole-to-halo mass ratio. We found that the properties of dual AGN host halos, including halo mass, stellar mass, star formation rate, and gas mass, are generally consistent with those of single AGN halos, albeit tending toward the higher end of their respective property ranges. Our analysis uncovered a diverse array of evolutionary patterns among dual AGNs, including rapidly evolving systems, slower ones, and instances where SMBH mergers are ineffective.

Counter-rotation and slow precession in aligned eccentric nuclear discs due to gravitational wave recoil kicks

ABSTRACT The M31 nucleus contains a supermassive black hole embedded in a massive stellar disc of apsidally aligned eccentric orbits. It has recently been shown that this disc is slowly precessing at a rate consistent with zero. Here, we demonstrate using N-body methods that apsidally aligned eccentric discs can form with a significant ($\sim$0.5) fraction of orbits counter-rotating as the result of a gravitational wave recoil kick of merging supermassive black holes. Higher amplitude kicks map to a larger retrograde fraction in the surrounding stellar population, which in turn results in slow precession. We furthermore show that discs with significant counter-rotation are more stable (i.e. apsidal alignment is most pronounced and long lasting), more eccentric, and have the highest rates of stars entering the black hole’s tidal disruption radius.

Big Galaxies and Big Black Holes: The Massive Ends of the Local Stellar and Black Hole Mass Functions and the Implications for Nanohertz Gravitational Waves

We construct the z = 0 galaxy stellar mass function (GSMF) by combining the GSMF at stellar masses M * ≲ 1011.3 M ⊙ from the census study of Leja et al. and the GSMF of massive galaxies at M * ≳ 1011.5 M ⊙ from the volume-limited MASSIVE galaxy survey. To obtain a robust estimate of M * for local massive galaxies, we use MASSIVE galaxies with M * measured from detailed dynamical modeling or stellar population synthesis modeling (incorporating a bottom-heavy initial mass function) with high-quality spatially resolved spectroscopy. These two independent sets of M * agree to within ∼7%. Our new z = 0 GSMF has a higher amplitude at M * ≳ 1011.5 M ⊙ than previous studies, alleviating prior concerns of a lack of mass growth in massive galaxies between z ∼ 1 and 0. We derive a local black hole mass function (BHMF) from this GSMF and the scaling relation of supermassive black holes (SMBHs) and galaxy masses. The inferred abundance of local SMBHs above ∼1010 M ⊙ is consistent with the number of currently known systems. The predicted amplitude of the nanohertz stochastic gravitational-wave background is also consistent with the levels reported by Pulsar Timing Array teams. Our z = 0 GSMF therefore leads to concordant results in the high-mass regime of the local galaxy and SMBH populations and the gravitational-wave amplitude from merging SMBHs. An exception is that our BHMF yields a z = 0 SMBH mass density that is notably higher than the value estimated from quasars at higher redshifts.

Noise analysis of six pulsars and a limit on the gravitational wave background

The direct detection of ultra-low frequency gravitational waves (GWs, 10−9 − 10−6 Hz) using pulsar timing arrays (PTAs) is a significant application of millisecond pulsars. The initial detection of GWs in the nano-Hertz range will probably be an isotropic stochastic GW background (GWB) emitted by the combination of merging super-massive black hole binaries (SMBHBs). In this study, we utilize data from six pulsars collected with four European PTA (EPTA) telescopes to place an upper limit on the GWB amplitude A=2.34−1.32+0.75×10−15 (95%) as a function of the spectral slope α = − 2/3, corresponding to a GWB generated by inspiraling SMBHBs. The posterior for the amplitude that we obtained from the six brightest EPTA pulsars is strongly consistent with previously published upper limits set by the three PTAs. We also compare the sensitivity of our reprocessed pulsar pulse times of arrival (ToAs) determined by the superior ToA creation method, which utilizes a combination of the best template generation algorithm and template-matching algorithm, with the ToAs determined by the EPTA groups. Our analysis shows that the reprocessed ToAs obtained with the superior ToA creation method yield narrower posteriors, indicating an overall increased sensitivity to GW.

Self-Interacting Dark Matter Solves the Final Parsec Problem of Supermassive Black Hole Mergers.

Evidence for a stochastic gravitational wave (GW) background, plausibly originating from the merger of supermassive black holes (SMBHs), is accumulating with observations from pulsar timing arrays. An outstanding question is how inspiraling SMBHs get past the "final parsec" of separation, where they have a tendency to stall before GW emission alone can make the binary coalesce. We argue that dynamical friction from the dark matter (DM) spike surrounding the black holes is sufficient to resolve this puzzle, if the DM has a self-interaction cross section of order cm^{2}/g. The same effect leads to a softening of the GW spectrum at low frequencies as suggested by the current data. For collisionless cold DM, the friction deposits so much energy that the spike is disrupted and cannot bridge the final parsec, while for self-interacting DM, the isothermal core of the halo can act as a reservoir for the energy liberated from the SMBH orbits. A realistic velocity dependence, such as generated by the exchange of a massive mediator like a dark photon, is favored to give a good fit to the GW spectrum while providing a large enough core. A similar velocity dependence has been advocated for solving the small-scale structure problems of cold DM.

Viscous torque in turbulent magnetized active galactic nucleus accretion disks and its effects on the gravitational waves of extreme mass ratio inspirals

The merger of supermassive black holes produces millihertz gravitational waves (GWs), which are potentially detectable by the future Laser Interferometer Space Antenna (LISA). Such binary systems are usually embedded in an accretion disk environment at the center of an active galactic nucleus (AGN). Recent studies suggest the plasma environment imposes measurable imprints on the GW signal if the mass ratio of the binary is around q ∼ 10−4 − 10−3. The effect of the gaseous environment on the GW signal is strongly dependent on the disk’s parameters; therefore, it is believed that future low-frequency GW detections will provide us with precious information about the physics of AGN accretion disks. We investigated this effect by measuring the viscous torque via modeling of the evolution of magnetized tori around the primary massive black hole. Using the general relativistic magnetohydrodynamic HARM-COOL code, we performed 2D and 3D simulations of weakly magnetized, thin accretion disks, with a possible truncation and transition to advection-dominated accretion flow. We studied the angular momentum transport and turbulence generated by the magnetorotational instability. We quantified the disk’s effective alpha viscosity and its evolution over time. We applied our numerical results to quantify the relativistic viscous torque on a hypothetical low-mass secondary black hole via a 1D analytical approach, and we estimated the GW phase shift due to the gas environment.

From Seeds to Supermassive Black Holes: Capture, Growth, Migration, and Pairing in Dense Protobulge Environments

The origins and mergers of supermassive black holes (SMBHs) remain a mystery. We describe a scenario from a novel multiphysics simulation featuring rapid (≲1 Myr) hyper-Eddington gas capture by a ∼1000 M ⊙ “seed” black hole (BH) up to supermassive (≳106 M ⊙) masses in a massive, dense molecular cloud complex typical of high-redshift starbursts. Due to the high cloud density, stellar feedback is inefficient, and most of the gas turns into stars in star clusters that rapidly merge hierarchically, creating deep potential wells. Relatively low-mass BH seeds at random positions can be “captured” by merging subclusters and migrate to the center in ∼1 freefall time (vastly faster than dynamical friction). This also efficiently produces a paired BH binary with ∼0.1 pc separation. The centrally concentrated stellar density profile (akin to a “protobulge”) allows the cluster as a whole to capture and retain gas and build up a large (parsec-scale) circumbinary accretion disk with gas coherently funneled to the central BH (even when the BH radius of influence is small). The disk is “hypermagnetized” and “flux-frozen”: dominated by a toroidal magnetic field with plasma β ∼ 10−3, with the fields amplified by flux-freezing. This drives hyper-Eddington inflow rates ≳1 M ⊙ yr−1, which also drive the two BHs to nearly equal masses. The late-stage system appears remarkably similar to recently observed high-redshift “little red dots.” This scenario can provide an explanation for rapid SMBH formation, growth, and mergers in high-redshift galaxies.

Can Supercooled Phase Transitions Explain the Gravitational Wave Background Observed by Pulsar Timing Arrays?

Several pulsar timing array collaborations recently reported evidence of a stochastic gravitational wave background (SGWB) at nHz frequencies. While the SGWB could originate from the merger of supermassive black holes, it could be a signature of new physics near the 100MeV scale. Supercooled first-order phase transitions (FOPTs) that end at the 100MeV scale are intriguing explanations, because they could connect the nHz signal to new physics at the electroweak scale or beyond. Here, however, we provide a clear demonstration that it is not simple to create a nHz signal from a supercooled phase transition, due to two crucial issues that could rule out many proposed supercooled explanations and should be checked. As an example, we use a model based on nonlinearly realized electroweak symmetry that has been cited as evidence for a supercooled explanation. First, we show that a FOPT cannot complete for the required transition temperature of around 100MeV. Such supercooling implies a period of vacuum domination that hinders bubble percolation and transition completion. Second, we show that even if completion is not required or if this constraint is evaded, the Universe typically reheats to the scale of any physics driving the FOPT. The hierarchy between the transition and reheating temperature makes it challenging to compute the spectrum of the SGWB.

Tracking Supermassive Black Hole Mergers from kpc to sub-pc Scales with AXIS

We present an analysis showcasing how the Advanced X-ray Imaging Satellite (AXIS), a proposed NASA Probe-class mission, will significantly increase our understanding of supermassive black holes undergoing mergers—from kpc to sub-pc scales. In particular, the AXIS point spread function, field of view, and effective area are expected to result in (1) the detection of hundreds to thousands of new dual AGNs across the redshift range 0<z<5 and (2) blind searches for binary AGNs that are exhibiting merger signatures in their light curves and spectra. AXIS will detect some of the highest-redshift dual AGNs to date, over a large range of physical separations. The large sample of AGN pairs detected by AXIS (over a magnitude more than currently known) will result in the first X-ray study that quantifies the frequency of dual AGNs as a function of redshift up to z=4.

Supermassive primordial black holes from inflation

There is controversy surrounding the origin and evolution of our universe's largest supermassive black holes (SMBHs). In this study, we consider the possibility that some of these black holes formed from the direct collapse of primordial density perturbations. Since the mass of a primordial black hole is limited by the size of the cosmological horizon at the time of collapse, these SMBHs must form rather late, and are naively in conflict with constraints from CMB spectral distortions. These limits can be avoided, however, if the distribution of primordial curvature perturbations is highly non-Gaussian. After quantifying the departure from Gaussianity needed to evade these bounds, we explore a model of multi-field inflation — a non-minimal, self-interacting curvaton model — which has all the necessary ingredients to yield such dramatic non-Gaussianities. We leave the detailed model building and numerics to a future study, however, as our goal is to highlight the challenges associated with forming SMBHs from direct collapse and to identify features that a successful model would need to have. This study is particularly timely in light of recent observations of high-redshift massive galaxy candidates by the James Webb Space Telescope as well as evidence from the NANOGrav experiment for a stochastic gravitational wave background consistent with SMBH mergers.

Primordial Black Holes from Spatially Varying Cosmological Constant Induced by Field Fluctuations in Extra Dimensions

The origin and evolution of supermassive black holes (SMBHs) in our universe have sparked controversy. In this study, we explore the hypothesis that some of these black holes may have seeded from the direct collapse of dark energy domains with density significantly higher than the surrounding regions. The mechanism of the origin of such domains relies on the inflationary evolution of a scalar field acting in D dimensions, which is associated with the cosmological constant in our four-dimensional spacetime manifold. Inner space quantum fluctuations of the field during inflation are responsible for the spatial variations of the dark energy density in our space. This finding holds particular significance, especially considering recent evidence from pulsar timing array observations, which supports the existence of a stochastic gravitational wave background consisting of SMBH mergers.

RABBITS – I. The crucial role of nuclear star formation in driving the coalescence of supermassive black hole binaries

ABSTRACT In this study of the ‘Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations’ (RABBITS) series, we focus on the hardening and coalescing process of supermassive black hole (SMBH) binaries in galaxy mergers. For simulations including different galaxy formation processes (i.e. gas cooling, star formation, SMBH accretion, stellar, and AGN feedback), we systematically control the effect of stochastic eccentricity by fixing it to similar values during the SMBH hardening phase. We find a strong correlation between the SMBH merger time-scales and the presence of nuclear star formation. Throughout the galaxy merging process, gas condenses at the centre due to cooling and tidal torques, leading to nuclear star formation. These recently formed stars, which inherit low angular momenta from the gas, contribute to the loss cone and assist in the SMBH hardening via three-body interactions. Compared to non-radiative hydrodynamical runs, the SMBH merger time-scales measured from the runs including cooling, stellar, and SMBH physical processes tend to be shortened by a factor of ∼1.7. After fixing the eccentricity to the range of e ∼ 0.6–0.8 during the hardening phase, the simulations with AGN feedback reveal merger time-scales of ∼100–500 Myr for disc mergers and ∼1–2 Gyr for elliptical mergers. With a semi-analytical approach, we find that the torque interaction between the binary and its circumbinary disc has minimal impact on the shrinking of the binary orbit in our retrograde galaxy merger. Our results are useful in improving the modelling of SMBH merger time-scales and gravitational-wave event rates.

Analysis of Kozai cycles in equal-mass hierarchical triple supermassive black hole mergers in the presence of a stellar cluster

ABSTRACT Supermassive black holes (SMBHs) play an important role in galaxy evolution. Binary and triple SMBHs can form after galaxy mergers. A third SMBH may accelerate the SMBH merging process, possibly through the Kozai mechanism. We use N-body simulations to analyse oscillations in the orbital elements of hierarchical triple SMBHs with surrounding star clusters in galaxy centres. We find that SMBH triples spend only a small fraction of time in the hierarchical merger phase (i.e. a binary SMBH with a distant third SMBH perturber). Most of the time, the enclosed stellar mass within the orbits of the innermost or the outermost SMBH is comparable to the SMBH masses, indicating that the influence of the surrounding stellar population cannot be ignored. We search for Eccentric Kozai–Lidov (EKL) oscillations for which (i) the eccentricity of the inner binary and inclination are both oscillate and are antiphase or in-phase and (ii) the oscillation period is consistent with EKL time-scale. We find that EKL oscillations are short-lived and rare: the triple SMBH spends around 3 per cent of its time in this phase over the ensemble of simulations, reaching around 8 per cent in the best-case scenario. This suggests that the role of the EKL mechanism in accelerating the SMBH merger process may have been overestimated in previous studies. We follow-up with three-body simulations, using initial conditions extracted from the simulation, and the result can to some extent repeat the observed EKL-like oscillations. This comparison provides clues about why those EKL oscillations with perturbing stars are short-lived.

Anisotropic Star Clusters around Recoiling Supermassive Black Holes

Gravitational-wave recoil kicks from merging supermassive black hole binaries can have a profound effect on the surrounding stellar population. In this work, we study the dynamic and kinematic properties of nuclear star clusters following a recoil kick. We show that these postkick structures present unique signatures that can provide key insight to observational searches for recoiling supermassive black holes. In our previous paper, we showed that an in-plane recoil kick turns a circular disk into a lopsided, eccentric disk such as the one we observe in the Andromeda nucleus. Building on this work, here we explore many recoil kick angles as well as initial stellar configurations. For a circular disk of stars, an in-plane kick causes strong apsidal alignment with a significant fraction of the disk becoming retrograde at large radii. If the initial orbits are highly eccentric, an in-plane kick forms a bar-like structure made up of two antialigned lopsided disks. An out-of-plane kick causes clustering in the argument of periapsis, ω, regardless of the initial eccentricity distribution. Initially, isotropic configurations form anisotropies in the form of a torus of eccentric orbits oriented perpendicular to the recoil kick. Postkick surface density and velocity maps are presented in each case to highlight the distinct, observable structures of these systems.

High-redshift supermassive black hole mergers in simulations with dynamical friction modelling

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.

A Candidate Dual QSO at Cosmic Noon

We report the discovery of a candidate dual QSO at z = 1.889, a redshift that is in the era known as “cosmic noon” where most of the universe’s black hole and stellar mass growth occurred. The source was identified in Hubble Space Telescope WFC3/IR images of a dust-reddened QSO that showed two closely separated point sources at a projected distance of 0.″26, or 2.2 kpc. This red QSO was targeted for imaging to explore whether red QSOs are hosted by merging galaxies. We subsequently obtained a spatially resolved Space Telescope Imaging Spectrograph spectrum of the system, covering the visible spectral range, and verifying the presence of two distinct QSO components. We also obtained high-resolution radio continuum observations with the Very Long Baseline Array at 1.4 GHz (21 cm L band) and found two sources coincident with the optical positions. The sources have similar black hole masses, bolometric luminosities, and radio-loudness parameters. However, their colors and reddenings differ significantly. The redder QSO has a higher Eddington ratio, consistent with previous findings. We consider the possibility of gravitational lensing and find that it would require extreme and unlikely conditions. If confirmed as a bona fide dual QSO, this system would link dust reddening to galaxy and supermassive black hole mergers, opening up a new population in which to search for samples of dual active galactic nuclei.

Aligning Retrograde Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disc

ABSTRACT Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination i < 90°) decrease inclination, as well as semimajor axis (a) and eccentricity (e) until orbital alignment with the gas disc (‘disc capture’). Captured stellar-origin black holes (sBH) add to the embedded AGN population that drives sBH–sBH mergers detectable in gravitational waves using LIGO–Virgo–KAGRA or sBH–SMBH mergers detectable with Laser Interferometer Space Antenna. Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive ‘immortal’ stars. Here, we investigate the behaviour of polar and retrograde orbiters (i ≥ 90°) interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation (i < 90°) during capture, and decrease a dramatically towards the SMBH. For sBH, we find a critical angle iret ∼ 113°, below which retrograde sBH decay towards embedded prograde orbits (i → 0°), while for io > iret sBH decay towards embedded retrograde orbits (i → 180°). sBH near polar orbits (i ∼ 90°) and stars on nearly embedded retrograde orbits (i ∼ 180°) show the greatest decreases in a. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial a. For sBH, disc capture time is longest for polar orbits, low-mass sBH, and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the spin distribution of embedded sBH.

Searching for clues of past binary supermassive black hole mergers in nuclear star clusters

ABSTRACTGalaxy mergers are common processes in the Universe. As a large fraction of galaxies hosts at their centres a central supermassive black hole (SMBH), mergers can lead to the formation of a supermassive black hole binary (SMBHB). The formation of such a binary is more efficient when the SMBHs are embedded in a nuclear star cluster (NSC). NSCs are dense and massive stellar clusters present in the majority of the observed galaxies. Their central densities can reach up to $10^7\, {\rm M_{\odot }}\,{\rm pc^{-3}}$ and their masses can be as large as a few $10^7\, {\rm M_{\odot }}$. The direct detection of an SMBHB is observationally challenging. In this work, we illustrate how the large-scale structural and dynamical properties of an NSC can help to identify nucleated galaxies that recently went through a merger that possibly led to the formation of a central SMBHB. Our models show that the merger can imprint signatures on the shape, density profile, rotation, and velocity structure of the NSC. The strength of the signatures depends on the mass ratio between the SMBHs and on the orbital initial conditions of the merger. In addition, the number of hypervelocity stars produced in the mergers is linked to the SMBHB properties. The merger can also contribute to the formation of the nuclear stellar disc of the galaxy.

Chasing supermassive black hole merging events withAthenaandLISA

ABSTRACTThe European Space Agency is studying two large-class missions bound to operate in the decade of the 30s, and aiming at investigating the most energetic and violent phenomena in the Universe. Athena is poised to study the physical conditions of baryons locked in large-scale structures from the epoch of their formation, as well as to yield an accurate census of accreting supermassive black holes down to the epoch of reionization; LISA will extend the hunt for Gravitational Wave (GW) events to the hitherto unexplored mHz regime. We discuss in this paper the science that their concurrent operation could yield, and present possible Athena observational strategies. We focus on Supermassive (M$\lesssim 10^7\, \rm {M_\odot }$) Black Hole Mergers (SMBHMs), potentially accessible to Athena up to z ∼ 2. The simultaneous measurement of their electromagnetic (EM) and GW signals may enable unique experiments in the domains of astrophysics, fundamental physics, and cosmography, such as the magnetohydrodynamics of fluid flows in a rapidly variable space–time, the formation of coronae and jets in Active Galactic Nuclei, and the measurement of the speed of GW, among others. Key to achieve these breakthrough results will be the LISA capability of locating a SMBHM event with an error box comparable to, or better than the field-of-view of the Athena Wide Field Imager ($\simeq 0.4\,$ deg2) and Athena capability to slew fast to detect the source during the inspiral phase and the post-merger phase. Together, the two observatories will open in principle the exciting possibility of truly concurrent EM and GW studies of the SMBHMs