Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Thursday, July 16 until 12:00 AM ET on Friday 17 due to maintenance. We apologize for the inconvenience.


Title: Reliable Identification of Binary Supermassive Black Holes from Rubin Observatory Time-domain Monitoring
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.  more » « less
Award ID(s):
1945546 2108668 2106990 1909933 2020265
PAR ID:
10512857
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
965
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
34
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Binary supermassive black holes (SMBHs) are consequences of galaxy mergers and dominate the low-frequency gravitational-wave background. Finding binary SMBHs in existing time-domain observations has proven difficult, as their periodic, electromagnetic signals can be confused with the natural variability of single quasars. In this work, we investigate the effects of host-galaxy contamination and survey design (cadence and duration) on the detectability of binary SMBHs with the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We simulate millions of LSST light curves of single and binary quasars, with a distribution of quasar and host-galaxy properties motivated by empirical observations and the anticipated LSST detection space. We then apply simple sinusoidal curve fits as a potential computationally inexpensive detection method. We find that host-galaxy contamination will increase false-positive rates and decrease binary parameter recovery rates. Lower-mass, lower-luminosity binary systems are most likely to be negatively affected by host-galaxy contamination. We also find that monitoring duration affects binary detection more than survey effective cadence for this detection method. As the light-curve duration increases, false-positive rates are suppressed and binary parameter recovery rates, especially for binary periods, are improved. Increasing the light-curve duration from 5 to 10 yr shows the most dramatic improvement for successful binary detection and false-positive rejection, with additional improvement from extending the light-curve duration to 20 yr. The observation duration increase is especially critical for recovering binary periods that are longer than a decade. 
    more » « less
  2. 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. 
    more » « less
  3. 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. 
    more » « less
  4. ABSTRACT The Legacy Survey of Space and Time (LSST) by the Vera C. Rubin Observatory is expected to discover tens of millions of quasars. A significant fraction of these could be powered by coalescing massive black hole (MBH) binaries, since many quasars are believed to be triggered by mergers. We show that under plausible assumptions about the luminosity functions, lifetimes, and binary fractions of quasars, we expect the full LSST quasar catalogue to contain between 20 and 100 million compact MBH binaries with masses M = 105–9M⊙, redshifts z = 0–6, and orbital periods P = 1–70 d. Their light-curves are expected to be distinctly periodic, which can be confidently distinguished from stochastic red-noise variability, because LSST will cover dozens, or even hundreds of cycles. A very small subset of 10–150 ultracompact (P ≲ 1 d) binary quasars among these will, over ∼5–15 yr, evolve into the mHz gravitational-wave frequency band and can be detected by LISA. They can therefore be regarded as ‘LISA verification binaries’, analogous to short-period Galactic compact-object binaries. The practical question is how to find these handful of ‘needles in the haystack’ among the large number of quasars: this will likely require a tailored co-adding analysis optimized for this purpose. 
    more » « less
  5. Abstract Supermassive black hole binary (SMBHB) systems are expected to form as a consequence of galaxy mergers. At subparsec separations, SMBHBs can be identified as quasars with periodic variability, with previous periodicity searches uncovering significant candidates. However, these searches focused primarily on sinusoidal signals, while theoretical models and hydrodynamical simulations predict that binaries produce more complex non-sinusoidal pulse shapes. Here we examine the efficacy of the Lomb–Scargle periodogram (LSP; one of the most popular tools for periodicity searches in unevenly sampled lightcurves) to detect periodicities with a sawtooth shape mimicking results of hydrodynamical simulations. We simulate idealized well-sampled lightcurves, lightcurves that mimic the data in the Palomar Transient Factory (PTF) analyzed in M. Charisi et al. (2016), and lightcurves that resemble our expectations for single-band data in the upcoming Legacy Survey of Space and Time (LSST) of the Rubin Observatory. We approximate quasar variability with a damped random walk (DRW) model, inject sinusoidal and sawtooth pulse shapes, and assess their statistical significance. We find that in the presence of red noise, the LSP detects a relatively low fraction of the sinusoidal signals (∼45%, ∼24%, and ∼23%, in the PTF-like, idealized, and LSST-like lightcurves, respectively). The fraction is significantly reduced for sawtooth periodicity (with only ∼9% in PTF-like and ∼1% in idealized and LSST-like lightcurves). These low recovery rates imply that previous searches have missed the large majority of binaries. They also have significant implications for the detection of SMBHBs in upcoming LSST necessitating the development of advanced tools that go beyond the simple LSP. 
    more » « less