The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe.
Gravitational-wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. The ngVLA offers unparalleled capabilities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales, opening a new avenue for multi-messenger studies in the dawn of low-frequency GW astronomy.
The nature of the first seeds of supermassive black holes (SMBHs) is currently unknown, with postulated initial masses ranging from $\sim10^5~M_{\odot}$ to as low as $\sim10^2~M_{\odot}$. However, most existing cosmological simulations resolve BHs only down to $\sim10^5-10^6~M_{\odot}$. In this work, we introduce a novel sub-grid BH seed model that is directly calibrated from high resolution zoom simulations that can trace the formation and growth of $\sim 10^3~M_{\odot}$ seeds forming in halos with pristine, star-forming gas. We trace the BH growth along merger trees until their descendants reach masses of $\sim10^4$ or $10^5~M_{\odot}$. The descendants assemble in galaxies with a broad range of properties (e.g., halo masses $\sim10^7-10^9~M_{\odot}$) that evolve with redshift and are sensitive to seed parameters. The results are used to build a new stochastic seeding model that directly seeds these descendants in lower resolution versions of our zoom region. Remarkably, we find that by seeding the descendants simply based on total galaxy mass, redshift and an environmental richness parameter, we can reproduce the results of the detailed gas based seeding model. The baryonic properties of the host galaxies are well reproduced by the mass-based seeding criterion. The redshift-dependence of the mass-based criterion captures the influence of halo growth, star formation and metal enrichment on seed formation. The environment based seeding criterion seeds the descendants in rich environments with higher numbers of neighboring galaxies. This accounts for the impact of unresolved merger dominated growth of BHs, which produces faster growth of descendants in richer environments with more extensive BH merger history. Our new seed model will be useful for representing a variety of low mass seeding channels within next generation larger volume uniform cosmological simulations.
We model continuous gravitational wave (CW) sources detectable by pulsar timing arrays (PTAs), using the ASTRID cosmological simulation. The most detectable single sources are in the low-frequency bins and are produced by the most supermassive black hole (SMBH) mergers (with M_BH ~ 10^10 - 10^11 solarmass) in the most massive galaxies M_star ~ 10^12 solarmass. Remarkably, these mergers (mostly at z< 0.3) occur within massive galaxies residing at the center of galaxy clusters. Particularly striking in ASTRID is a triple merger event, where two subsequent mergers in the same cluster core generate high-detection-probability CW signals at ~ 2nHz and ~ 10nHz. We expect electromagnetic signatures from these events: either single or dual active galactic nuclei (AGN) and massive host galaxies that are undergoing significant star formation. We establish a novel connection between these high-mass, central cluster mergers, occurring in galaxies with active star formation and AGN, and the formation process of the central cluster galaxy in ASTRID. This provides new insights into the low-frequency gravitational wave sky and informs future multi-messenger searches for PTA CW sources.
Abstract The coevolution of galaxies and the supermassive black holes (SMBHs) at their centers via hierarchical galaxy mergers is a key prediction of ΛCDM cosmology. As gas and dust are funneled to the SMBHs during the merger, the SMBHs light up as active galactic nuclei (AGNs). In some cases, a merger of two galaxies can encounter a third galaxy, leading to a triple merger, which would manifest as a triple AGN if all three SMBHs are simultaneously accreting. Using high spatial resolution X-ray, near-IR, and optical spectroscopic diagnostics, we report here a compelling case of an AGN triplet with mutual separations <10 kpc in the advanced merger SDSS J084905.51+111447.2 at z = 0.077. The system exhibits three nuclear X-ray sources, optical spectroscopic line ratios consistent with AGN in each nucleus, a high excitation near-IR coronal line in one nucleus, and broad Pa α detections in two nuclei. Hard X-ray spectral fitting reveals a high column density along the line of sight, consistent with the picture of late-stage mergers hosting heavily absorbed AGNs. Our multiwavelength diagnostics support a triple AGN scenario, and we rule out alternative explanations such as star formation activity, shock-driven emission, and emission from fewer than three AGN. The dynamics of gravitationally bound triple SMBH systems can dramatically reduce binary SMBH inspiral timescales, providing a possible means to surmount the “Final Parsec Problem.” AGN triplets in advanced mergers are the only observational forerunner to bound triple SMBH systems and thus offer a glimpse of the accretion activity and environments of the AGNs prior to the gravitationally bound triple phase.
Abstract Recently we found compelling evidence for a gravitational-wave background with Hellings and Downs (HD) correlations in our 15 yr data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes, which produce different interpulsar correlations. In this work, we search the NANOGrav 15 yr data set for evidence of a gravitational-wave background with quadrupolar HD and scalar-transverse (ST) correlations. We find that HD correlations are the best fit to the data and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors ∼2 when comparing HD to ST correlations, and ∼1 for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.
Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.
Abstract Evidence for a low-frequency stochastic gravitational-wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these data sets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15 yr data set. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of ∼4 nHz, and another at ∼170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings–Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems. At our most sensitive frequency of 6 nHz, we place a sky-averaged 95% upper limit of 8 × 10 −15 on the strain amplitude. We also calculate an exclusion volume and a corresponding effective radius, within which we can rule out the presence of black hole binaries emitting at a given frequency.
