Supermassive Black Hole Binaries Research Articles

On the anisotropies of the cosmological gravitational-wave background from pulsar timing array observations

Significant evidence for a stochastic gravitational-wave background has recently been reported by several Pulsar Timing Array observations. These studies have shown that, in addition to astrophysical explanations based on supermassive black hole binaries (SMBHBs), cosmological origins are considered equally important sources for these signals. To further explore these cosmological sources, in this study, we discuss the anisotropies in the cosmological gravitational wave background (CGWB) in a model-independent way. Taking the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 15-year dataset as a benchmark, we estimate the angular power spectra of the CGWB and their cross-correlations with cosmic microwave background (CMB) fluctuations and weak gravitational lensing. We find that the NANOGrav 15-year data implies suppressed Sachs-Wolf (SW) effects in the CGBW spectrum, leading to a marginally negative cross-correlation with the CMB at large scales. This procedure is applicable to signals introduced by different early universe processes and is potentially useful for identifying unique features about anisotropies of CGWB from future space-based interferometers and astrometric measurements.

Gravitational torque in circumbinary discs: global radial oscillations

ABSTRACT Circumbinary discs (CBDs) arise in many astrophysical settings, including young stellar binaries and supermassive black hole binaries. Their structure is mediated by gravitational torques exerted on the disc by the central binary. The spatial distribution of the binary torque density (so-called excitation torque density) in CBDs is known to feature global large-amplitude, quasi-periodic oscillations, which are often interpreted in terms of the local resonant Lindblad torques. Here, we investigate the nature of these torque oscillations using 2D, inviscid hydrodynamic simulations and theoretical calculations. We show that torque oscillations arise due to the gravitational coupling of the binary potential to the density waves launched near the inner cavity and freely propagating out in the disc. We provide analytical predictions for the radial periodicity of the torque density oscillations and verify them with simulations, showing that disc sound speed and the multiplicity of the density wave spiral arms are the key factors setting the radial structure of the oscillations. Resonant Lindblad torques play no direct role in determining the radial structure and periodicity of the torque oscillations and manifest themselves only by driving the density waves in the disc. We also observe the formation of vortices at the inner edge of the disc, which can provide a non-trivial contribution to the angular momentum transport in the CBD and may be involved in the development of a non-axisymmetric central cavity.

Gravitational waves from supermassive black hole binaries in light of the NANOGrav 15-year data

Gravitational waves from supermassive black hole binaries in light of the NANOGrav 15-year data

Deep learning for parameter estimation of supermassive binary black holes with simulated LISA data

Previous studies have demonstrated that deep learning is an effective approach for processing gravitational wave (GW) data obtained from ground-based detectors, as it can enhance the efficiency of data processing and offer great potential for real-time parameter estimation. In this paper, we explore three different deep learning architectures (MCD, TFP, and VI) for inferring the chirp mass (Mc) and the luminosity distance (DL) of supermassive binary black holes (SMBBHs) using simulated data from the Laser Interferometer Space Antenna (LISA). We train the neural networks with the simulated data and evaluate their performance on predicting the parameters (Mc,DL). The results show that more than 97.5% of the true values fall within the 3σ confidence interval of the predicted values with the optimal network. To verify the accuracy of the network on parameter estimation, we also calculate the estimation error of the parameters using the Fisher Information Matrix (FIM) with the same simulated data. By comparing the root-mean-square error (RMSE) between deep learning and FIM, we find that the two methods are comparable, which implies that deep learning can be reliable for GW parameter estimation.

On the evolution of a twisted thin accretion disc in eccentric inclined supermassive binary black holes

ABSTRACT We propose a model of a twisted accretion disc around a Kerr black hole interacting with a secondary black hole of a smaller mass on an inclined eccentric orbit. We use parameters of the system, which may be appropriate for the so-called precessing massive model of OJ 287. We calculate expressions for torque exerted on the disc by the secondary and a contribution of the secondary to the apsidal precession of disc elements by a double averaging procedure over the periods of the secondary and the disc elements. These expressions are used at all scales of interest, including the ones inside the binary orbit. We calculate numerically the evolution of the disc tilt and twist assuming a flat initial configuration. We consider the disc aspect ratio h/r = 10−3, a rather large viscosity parameter α = 0.1 and several values of the primary rotational parameter, χ. We find that, after a few periods of Lense–Thirring precession of the orbit, the disc relaxes to a quasi-stationary configuration in the precessing frame with a non-trivial distribution of the disc inclination angle, β, over the radial scale. We propose an analytic model for this configuration. We show that the presence of the twisted disc leads to multiple crossings of the disc by the secondary per one orbital period, with time periods between the crossings being different from the flat disc model. Our results should be taken into account in the modelling of OJ 287. They can also be applied to similar sources.

How Many Supermassive Black Hole Binaries Are Detectable through Tracking Relative Motions by (Sub)millimeter Very Long Baseline Interferometry?

The (sub)millimeter wavelengths (86–690 GHz) very long baseline interferometry will provide ∼5–40 μas angular resolution, ∼10 mJy baseline sensitivity, and ∼1 μas yr−1 proper-motion precision, which can directly detect supermassive black hole binary (SMBHB) systems by imaging the two visible sources and tracking their relative motions. Such a way exhibits an advantage compared to indirect detect methods of observing periodic signals in motion and light curves, which are difficult to confirm from competing models. Moreover, tracking relative motion at (sub)millimeter wavelengths is more reliable, as there is a negligible offset between the emission region and the black hole center. In this way, it is unnecessary to correct the black hole location by a prior of jet morphology as it would be required at longer wavelengths. We extend the formalism developed in D’Orazio & Loeb (2018) to link the observations with the orbital evolution of SMBHBs from the ≲10 kpc dynamical friction stages to the ≲0.01 pc gravitational radiation stages, and estimate the detectable numbers of SMBHBs. By assuming 5% of active galactic nuclei holding SMBHBs, we find that the number of detectable SMBHBs with redshift z ≤ 0.5 and mass M ≤ 1011 M ⊙ is about 20. Such a detection relies heavily on proper-motion precision and sensitivity. Furthermore, we propose that the simultaneous multifrequency technique plays a key role in meeting the observational requirements.

The Central Kinematics and Black Hole Mass of 4C+37.11

We report on integral field unit (IFU) measurements of the host of the radio source 4C+37.11. This massive elliptical contains the only resolved double compact nucleus at parsec-scale separation, likely a bound supermassive black hole binary (SMBHB). i-band photometry and GMOS-N IFU spectroscopy show that the galaxy has a large r b = 1.″5 core and that the stellar velocity dispersion increases inside of a radius of influence r SOI ≈ 1.″3. Jeans Anisotropic Modeling analysis of the core infers a total SMBHB mass of , making this one of the most massive black hole systems known. Our data indicate that there has been significant scouring of the central kiloparsec of the host galaxy.

Tidal Disruption Events from the Combined Effects of Two-body Relaxation and the Eccentric Kozai–Lidov Mechanism

Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH’s sphere of influence, commonly referred to as two-body relaxation. This process can alter the angular momentum of stars at large distances and place them into nearly radial orbits, thus driving them to disruption. On the other hand, gravitational perturbations from an SMBH companion via the eccentric Kozai–Lidov (EKL) mechanism have also been proposed as a promising stellar disruption channel. Here we demonstrate that the combination of EKL and two-body relaxation in SMBH binaries is imperative for building a comprehensive picture of the rates of TDEs. Here we explore how the density profile of the surrounding stellar distribution and the binary orbital parameters of an SMBH companion influence the rate of TDEs. We show that this combined channel naturally produces disruptions at a rate that is consistent with observations and also naturally forms repeated TDEs, where a bound star is partially disrupted over multiple orbits. Recent observations show stars being disrupted in short-period orbits, which is challenging to explain when these mechanisms are considered independently. However, the diffusive effect of two-body relaxation, combined with the secular nature of the eccentricity excitations from EKL, is found to drive stars on short eccentric orbits at a much higher rate. Finally, we predict that rTDEs are more likely to take place in the presence of a steep stellar density distribution.

Disentangling the primordial nature of stochastic gravitational wave backgrounds with CMB spectral distortions

ABSTRACT The recent detection of a stochastic gravitational wave background (SGWB) at nanohertz frequencies by pulsar timing arrays (PTAs) has sparked a flurry of interest. Beyond the standard interpretation that the progenitor is a network of supermassive black hole binaries, many exotic models have also been proposed, some of which can potentially offer a better fit to the data. We explore how the various connections between gravitational waves (GWs) and cosmic microwave background (CMB) spectral distortions (SDs) can be leveraged to help determine whether an SGWB was generated primordially or astrophysically. To this end, we present updated k-space window functions that can be used for distortion parameter estimation on enhancements to the primordial scalar power spectrum. These same enhancements can also source GWs directly at second order in perturbation theory, so-called scalar-induced GWs (SIGWs), and indirectly through the formation of primordial black holes (PBHs). We perform a mapping of scalar power spectrum constraints into limits on the GW parameter space of SIGWs for δ-function features. We highlight that broader features in the scalar spectrum can explain the PTA results while simultaneously producing an SD within reach of future experiments. We additionally update PBH constraints from μ- and y-type SDs. Refined treatments of the distortion window functions widen existing SD constraints, and we find that a future CMB spectrometer could play a pivotal role in unravelling the origin of GWs imprinted at or below CMB anisotropy scales.

Probing the Mass Relation between Supermassive Black Holes and Dark Matter Halos at High Redshifts by Gravitational Wave Experiments

Numerous observations have shown that almost all galaxies in our Universe host supermassive black holes (SMBHs), but there is still much debate about their formation and evolutionary processes. Recently, gravitational waves (GWs) have been expected to be a new and important informative observation, in particular, in the low-frequency region by making use of the Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays (PTAs). As an evolutionary process of the SMBHs, we revisit a dark matter (DM) halo–SMBH coevolution model based on the halo merger tree employing an ansatz for the mass relation between the DM halos and the SMBHs at z = 6. In this model, the mass of SMBHs grows through their mergers associated with the halo mergers, and hence, the evolutionary information must be stored in the GWs emitted at the mergers. We investigate the stochastic gravitational background from the coalescing SMBH binaries, which the PTAs can detect, and also the GW bursts emitted at the mergers, which can be detected by the mHz band observations such as LISA. We also discuss the possibility of probing the mass relation between the DM halos and the SMBHs at high redshift by future GW observations.

Uncovering Hidden Massive Black Hole Companions with Tidal Disruption Events

Dynamical perturbations from supermassive black hole (SMBH) binaries can increase the rates of tidal disruption events (TDEs). However, most previous work focuses on TDEs from the heavier black hole in the SMBH binary (SMBHB) system. In this work, we focus on the lighter black holes in SMBHB systems and show that they can experience a similarly dramatic increase in their TDE rate due to perturbations from a more massive companion. While the increase in TDEs around the more massive black hole is mostly due to chaotic orbital perturbations, we find that, around the smaller black hole, the eccentric Kozai–Lidov mechanism is dominant and capable of producing a comparably large number of TDEs. In this scenario, the mass derived from the light curve and spectra of TDEs caused by the lighter SMBH companion is expected to be significantly smaller than the SMBH mass estimated from galaxy scaling relations, which are dominated by the more massive companion. This apparent inconsistency can help find SMBHB candidates that are not currently accreting as active galactic nuclei and that are at separations too small for them to be resolved as two distinct sources. In the most extreme cases, these TDEs provide us with the exciting opportunity to study SMBHBs in galaxies where the primary SMBH is too massive to disrupt Sun-like stars.

Observational Evidence of a Centi-parsec Supermassive Black Hole Binary Existing in the Nearby Galaxy M81

Studying a centi-parsec supermassive black hole binary (SMBHB) would allow us to explore a new parameter space in active galactic nuclei, and these objects are also potential sources of gravitational waves. We report evidence that an SMBHB with an orbital period of ∼30 yr may be resident in the nearby galactic nucleus M81. This orbital period and the known mass of M81 imply an orbital separation of ∼0.02 pc. The jet emanating from the primary black hole showed a short period of jet wobbling at ∼16.7 yr, superposing a long-term precession at a timescale of several hundred years. Periodic radio and X-ray outbursts were also found two times per orbital period, which could be explained by a double-peaked mass accretion rate variation per binary orbit. If confirmed, M81 would be one of the closest SMBHB candidates, providing a rare opportunity to study the final parsec problem.

How to Detect an Astrophysical Nanohertz Gravitational Wave Background

Analyses of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nanohertz frequency band. The most plausible source of this background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for this background and assess its significance make several simplifying assumptions, namely (i) Gaussianity, (ii) isotropy, and most often, (iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated data sets. The data-set length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15 yr data set. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated data sets, even though their fundamental assumptions are not strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them.

Axionic domain walls at Pulsar Timing Arrays: QCD bias and particle friction

The recent results from the Pulsar Timing Array (PTA) collaborations show the first evidence for the detection of a stochastic background of gravitational waves at the nHz frequencies. This discovery has profound implications for the physics of both the late and the early Universe. In fact, together with the interpretation in terms of supermassive black hole binaries, many sources in the early Universe can provide viable explanations as well. In this paper, we study the gravitational wave background sourced by a network of axion-like-particle (ALP) domain walls at temperatures around the QCD crossover, where the QCD-induced potential provides the necessary bias to annihilate the network. Remarkably, this implies a peak amplitude at frequencies around the sensitivity range of PTAs. We extend previous analysis by taking into account the unavoidable friction on the network stemming from the topological coupling of the ALP to QCD in terms of gluon and pion reflection off the domain walls at high and low temperatures, respectively. We identify the regions of parameter space where the network annihilates in the scaling regime ensuring compatibility with the PTA results, as well as those where friction can be important and a more detailed study around the QCD crossover is required.

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.

Implications for the supermassive black hole binaries from the NANOGrav 15-year data set

Implications for the supermassive black hole binaries from the NANOGrav 15-year data set

Does NANOGrav observe a dark sector phase transition?

Gravitational waves from a first-order cosmological phase transition, at temperatures at the MeV-scale, would arguably be the most exciting explanation of the common red spectrum reported by the NANOGrav collaboration, not the least because this would be direct evidence of physics beyond the standard model. Here we perform a detailed analysis of whether such an interpretation is consistent with constraints on the released energy deriving from big bang nucleosynthesis and the cosmic microwave background. We find that a phase transition in a completely secluded dark sector with sub-horizon sized bubbles is strongly disfavoured with respect to the more conventional astrophysical explanation of the putative gravitational wave signal in terms of supermassive black hole binaries. On the other hand, a phase transition in a dark sector that subsequently decays, before the time of neutrino decoupling, remains an intriguing possibility to explain the data. From the model-building perspective, such an option is easily satisfied for couplings with the visible sector that are small enough to evade current collider and astrophysical constraints. The first indication that could eventually corroborate such an interpretation, once the observed common red spectrum is confirmed as a nHz gravitational wave background, could be the spectral tilt of the signal. In fact, the current data already show a very slight preference for a spectrum that is softer than what is expected from the leading astrophysical explanation.

Shock-driven periodic variability in a low-mass-ratio supermassive black hole binary

ABSTRACT We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disc, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of an SMBHB with mass ratio q = 0.01 and separation a = 100 Rg, using a physically self-consistent steady disc model. We estimate the electromagnetic variability from the system by monitoring accretion on to the secondary and using an artificial viscosity scheme to capture shocks and monitor the energy dissipated. The SMBHB produces a wide, eccentric gap in the disc, previously only observed for larger mass ratios, which we attribute to our disc model being much thinner (H/R ≈ 0.01 near the secondary) than is typical of previous works. The eccentric gap drives periodic accretion on to the secondary SMBH on a time-scale matching the orbital period of the binary, $t_{\rm {bin}}\approx 0.1\,\,\rm {yr}$, implying that the variable accretion regime of the SMBHB parameter space extends to lower mass ratios than previously established. Shocks driven by the binary are periodic, with a period matching the orbital period, and the shocks are correlated with the accretion rate, with peaks in the shock luminosity lagging peaks in the accretion rate by 0.43 tbin. We propose that the correlation of these quantities represents a useful identifier of SMBHB candidates, via observations of correlated variability in X-ray and ultraviolet monitoring of active galactic nuclei, rather than single-waveband periodicity alone.

Quasi-simultaneous Optical Flux and Polarization Variability of the Binary Super Massive Black Hole Blazar OJ 287 from 2015 to 2023: Detection of an Anticorrelation in Flux and Polarization Variability

We study the optical flux and polarization variability of the binary black hole blazar OJ 287 using quasi-simultaneous observations from 2015 to 2023 carried out using telescopes in the USA, Japan, Russia, Crimea, and Bulgaria. This is one of the most extensive quasi-simultaneous optical flux and polarization variability studies of OJ 287. OJ 287 showed large amplitude, ∼3.0 mag flux variability, large changes of ∼37% in degree of polarization, and a large swing of ∼215° in the angle of the electric vector of polarization. During the period of observation, several flares in flux were detected. Those flares are correlated with a rapid increase in the degree of polarization and swings in electric vector of polarization angle. A peculiar behavior of anticorrelation between flux and polarization degree, accompanied by a nearly constant polarization angle, was detected from JD 2,458,156 to JD 2,458,292. We briefly discuss some explanations for the flux and polarization variations observed in OJ 287.

Limits on scalar-induced gravitational waves from the stochastic background by pulsar timing array observations

Recently, the NANOGrav, PPTA, EPTA, and CPTA Collaborations independently reported their evidence of the Stochastic Gravitational Waves Background (SGWB). While the inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from the population of supermassive black-hole binaries (SMBHBs), the search for new physics remains plausible in this observational window. In this work, we explore the possibility of explaining such a signal by the scalar-induced gravitational waves (IGWs) in the very early universe. We use a parameterized broken power-law function as a general description of the energy spectrum of the SGWB and fit it to the new results of NANOGrav, PPTA and EPTA. We find that this approach can put constraints on the parameters of IGW energy spectrum and further yield restrictions on various inflation models that may produce primordial black holes (PBHs) in the early universe, which is also expected to be examined by the forthcoming space-based GW experiments.