Search for: All records

Abstract Supermassive black hole binaries (SMBHBs) are thought to form in galaxy mergers, possessing the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identification of the host galaxy will be a major challenge due to the ambiguity in possible EM signatures and the poor localization capability of PTAs. To aid EM observations in choosing follow-up sources, we use NANOGrav’s galaxy catalog to quantify the number of plausible hosts in both realistic and idealistic scenarios. We outline a host identification pipeline that injects a single-source GW signal into a simulated PTA data set, recovers the signal using production-level techniques, quantifies the localization region and number of galaxies contained therein, and finally imposes cuts on the galaxies using parameter estimates from the GW search. In an ideal case, the 90% credible areas span 29–241 deg², containing about 14–341 galaxies. After cuts, the number of galaxies remaining ranges from 22 at worst to one true host at best. In a realistic case, these areas range from 287 to 530 deg²and enclose about 285–1238 galaxies. After cuts, the number of galaxies is 397 at worst and 27 at best. While the signal-to-noise ratio is the primary determinant of the localization area of a given source, we find that the area is also influenced by the proximity to nearby pulsars on the sky and the binary chirp mass. 

Abstract Periodic signatures in time-domain observations of quasars have been used to search for binary supermassive black holes (SMBHs). These searches, across existing time-domain surveys, have produced several hundred candidates. The general stochastic variability of quasars, however, can masquerade as a false-positive periodic signal, especially when monitoring cadence and duration are limited. In this work, we predict the detectability of binary SMBHs in the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We apply computationally inexpensive sinusoidal curve fits to millions of simulated LSST Deep Drilling Field light curves of both single, isolated quasars and binary quasars. The period and phase of simulated binary signals can generally be disentangled from quasar variability. Binary amplitude is overestimated and poorly recovered for two-thirds of potential binaries due to quasar accretion variability. Quasars with strong intrinsic variability can obscure a binary signal too much for recovery. We also find that the most luminous quasars mimic current binary candidate light curves and their properties: The false-positive rates are 60% for these quasars. The reliable recovery of binary period and phase for a wide range of input binary LSST light curves is promising for multi-messenger characterization of binary SMBHs. However, pure electromagnetic detections of binaries using photometric periodicity with amplitude greater than 0.1 mag will result in samples that are overwhelmed by false positives. This paper represents an important and computationally inexpensive way forward for understanding the true and false-positive rates for binary candidates identified by Rubin. 

Abstract Supermassive black hole binaries (SMBHBs) are an inevitable consequence of galaxy mergers. At sub-parsec separations, they are practically impossible to resolve, and the most promising technique is to search for quasars with periodic variability. However, searches for quasar periodicity in time-domain data are challenging due to the stochastic variability of quasars. In this paper, we used Bayesian methods to disentangle periodic SMBHB signals from intrinsic damped random walk (DRW) variability in active galactic nuclei light curves. We simulated a wide variety of realistic DRW and DRW+sine light curves. Their observed properties are modeled after the Catalina Real-time Transient Survey (CRTS) and expected properties of the upcoming Legacy Survey of Space and Time (LSST) from the Vera C. Rubin Observatory. Through a careful analysis of parameter estimation and Bayesian model selection, we investigated the range of parameter space for which binary systems can be detected. We also examined which DRW signals can mimic periodicity and be falsely classified as binary candidates. We found that periodic signals are more easily detectable if the period is short or the amplitude of the signal is large compared to the contribution of the DRW noise. We saw similar detection rates both in the CRTS and LSST-like simulations, while the false-detection rate depends on the quality of the data and is minimal in LSST. Our idealized simulations provide an excellent way to uncover the intrinsic limitations in quasar periodicity searches and set the stage for future searches for SMBHBs. 

ABSTRACT Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g. Active Galactic Nuclei (AGNs) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >108M⊙) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multimessenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multimessenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly. 

Supermassive black hole binaries (SMBHBs) are an inevitable consequence of galaxy mergers. At sub-parsec separations, they are practically impossible to resolve and the most promising technique is to search for quasars with periodic variability. However, searches for quasar periodicity in time-domain data are challenging due to the stochastic variability of quasars. In this paper, we used Bayesian methods to disentangle periodic SMBHB signals from intrinsic damped random walk (DRW) variability in AGN light curves. We simulated a wide variety of realistic DRW and DRW+sine light curves. Their observed properties are modeled after the Catalina Real-time Transient Survey (CRTS) and expected properties of the upcoming Legacy Survey of Space and Time (LSST) from the Vera C. Rubin Observatory. Through a careful analysis of parameter estimation and Bayesian model selection, we investigated the range of parameter space for which binary systems can be detected. We also examined which DRW signals can mimic periodicity and be falsely classified as binary candidates. We found that periodic signals are more easily detectable if the period is short or the amplitude of the signal is large compared to the contribution of the DRW noise. We saw similar detection rates both in the CRTS and LSST-like simulations, while the false detection rate depends on the quality of the data and is minimal in LSST. Our idealized simulations provide an excellent way to uncover the intrinsic limitations in quasar periodicity searches and set the stage for future searches for SMBHBs. 

ABSTRACT We present a study of optically selected dual Active Galactic Nuclei (AGN) with projected separations of 3–97 kpc. Using multiwavelength (MWL) information (optical, X-ray, mid-IR), we characterized the intrinsic nuclear properties of this sample and compared them with those of isolated systems. Among the 124 X-ray-detected AGN candidates, 52 appear in pairs and 72 as single X-ray sources. Through MWL analysis, we confirmed the presence of the AGN in >80 per cent of the detected targets in pairs (42 out of 52). X-ray spectral analysis confirms the trend of increasing AGN luminosity with decreasing separation, suggesting that mergers may have contributed to triggering more luminous AGN. Through X-ray/mid-IR ratio versus X-ray colours, we estimated a fraction of Compton-thin AGN (with 1022 cm−2 < NH < 1024 cm−2) of about 80 per cent, while about 16 per cent are Compton-thick sources (with NH > 1024 cm−2). These fractions of obscured sources are larger than those found in samples of isolated AGN, confirming that pairs of AGN show higher obscuration. This trend is further confirmed by comparing the de-reddened [O iii] emission with the observed X-ray luminosity. However, the derived fraction of Compton-thick sources in this sample at the early stages of merging is lower than that reported for late-merging dual-AGN samples. Comparing NH from X-rays with that derived from E(B − V) from narrow-line regions, we found that the absorbing material is likely to be associated with the torus or broad-line regions. We also explored the X-ray detection efficiency of dual-AGN candidates, finding that, when observed properly (at on-axis positions and with long exposures), X-ray data represent a powerful way to confirm and investigate dual-AGN systems. 

Abstract Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15 yr pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese pulsar timing arrays. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational-wave background (GWB). To explore the internal consistency of this result, we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure the correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic GWB in the NANOGrav 15 yr data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data. 

ABSTRACT The variability of quasars across multiple wavelengths is a useful probe of physical conditions in active galactic nuclei. In particular, variable accretion rates, instabilities, and reverberation effects in the accretion disc of a supermassive black hole are expected to produce correlated flux variations in ultraviolet (UV) and optical bands. Recent work has further argued that binary quasars should exhibit strongly correlated UV and optical periodicities. Strong UV–optical correlations have indeed been established in small samples of (N ≲ 30) quasars with well-sampled light curves, and have extended the ‘bluer-when-brighter’ trend previously found within the optical bands. Here, we further test the nature of quasar variability by examining the observed-frame UV–optical correlations among bright quasars extracted from the Half Million Quasars (HMQ) catalogue. We identified a large sample of 1315 quasars in HMQ with overlapping UV and optical light curves from the Galaxy Evolution Explorer and the Catalina Real-time Transient Survey, respectively. We find that strong correlations exist in this much larger sample, but we rule out, at ∼95 per cent confidence, the simple hypothesis that the intrinsic UV and optical variations of all quasars are fully correlated. Our results therefore imply the existence of physical mechanism(s) that can generate uncorrelated optical and UV flux variations.