Reliable Identification of Binary Supermassive Black Holes from Rubin Observatory Time-Domain Monitoring

Abstract

Periodic signatures in time-domain observations of quasars have been used to search for binary supermassive black holes. 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 supermassive black holes 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. 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: 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 supermassive black holes. However, pure electromagnetic detections of binaries using photometric periodicity with amplitude greater than 0.1 magnitude 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.

21 pages, 14 figures, 3 tables

References in corpus (25)

• The NumPy array: a structure for efficient numerical computation
• The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background
• A Catalog of Quasar Properties from SDSS DR7
• Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array
• An Observational Determination of the Bolometric Quasar Luminosity Function
• Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk
• Measuring the rotation period distribution of field M-dwarfs with Kepler
• A Compact Supermassive Binary Black Hole System
• The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
• Is Quasar Optical Variability a Damped Random Walk?
• The second data release from the European Pulsar Timing Array I. The dataset and timing analysis
• Viscous Hydrodynamics Simulations of Circumbinary Accretion Discs: Variability, Quasi-Steady State, and Angular Momentum Transfer
• A Catalog of Quasar Properties from Sloan Digital Sky Survey Data Release 16
• The Local Nanohertz Gravitational-Wave Landscape From Supermassive Black Hole Binaries
• On the variability of quasars: a link between Eddington ratio and optical variability?
• An all-sky search for continuous gravitational waves in the Parkes Pulsar Timing Array data set
• Limitations on the recovery of the true AGN variability parameters using Damped Random Walk modeling
• Constraining the Orbit of Supermassive Black Hole Binary 0402+379
• Improving Damped Random Walk parameters for SDSS Stripe 82 Quasars with Pan-STARRS1
• The ANTARES Astronomical Time-Domain Event Broker
• UGC 4211: A Confirmed Dual Active Galactic Nucleus in the Local Universe at 230 pc Nuclear Separation
• Gravitational Self-Lensing in Populations of Massive Black Hole Binaries
• A Method to Measure the Unbiased Decorrelation Timescale of the AGN Variable Signal from Structure Functions
• ANTARES: Progress towards building a `Broker' of time-domain alerts
• Eccentric Minidisks in Accreting Binaries

Cited by in corpus (5)

• Time Evolution of Mg II in SDSS J2320+0024: Implications for a Subparsec Binary Supermassive Black Hole System
• Identifying massive black hole binaries via light curve variability in optical time-domain surveys
• Search for Quasar Pairs with $Gaia$ Astrometric Data. I. Method and Candidates
• Primordial Black Holes from Spatially Varying Cosmological Constant Induced by Field Fluctuations in Extra Dimensions
• Identifying Changing-Look AGN Transitions in Light Curve Data with the Zwicky Transient Facility