DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties

Journal Article · · The Astrophysical Journal
ORCiD logo [1]; ORCiD logo [2]; ORCiD logo [3]; ORCiD logo [4]; ORCiD logo [5]; ORCiD logo [6]; ORCiD logo [7];  [8]; ORCiD logo [3]
  1. City Univ. of New York (CUNY), NY (United States); American Museum of Natural History, New York, NY (United States); SLAC
  2. SLAC National Accelerator Laboratory (SLAC), Menlo Park, CA (United States)
  3. City Univ. of New York (CUNY), NY (United States); American Museum of Natural History, New York, NY (United States)
  4. American Museum of Natural History, New York, NY (United States); City Univ. of New York (CUNY), NY (United States)
  5. City Univ. of New York (CUNY), NY (United States); American Museum of Natural History, New York, NY (United States); CUNY Borough of Manhattan Community College, NY (United States); Flatiron Institute, New York, NY (United States)
  6. Flatiron Institute, New York, NY (United States); California Institute of Technology (CalTech), Pasadena, CA (United States)
  7. Harvard-Smithsonian Center for Astrophysics, Cambridge, MA (United States)
  8. Flatiron Institute, New York, NY (United States); Princeton Univ., NJ (United States); New York Univ. (NYU), NY (United States); Carnegie Mellon Univ., Pittsburgh, PA (United States)

Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar’s brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling.

Research Organization:
SLAC National Accelerator Laboratory (SLAC), Menlo Park, CA (United States)
Sponsoring Organization:
USDOE Office of Science (SC)
Grant/Contract Number:
AC02-76SF00515
OSTI ID:
2396968
Journal Information:
The Astrophysical Journal, Journal Name: The Astrophysical Journal Journal Issue: 2 Vol. 965; ISSN 0004-637X
Publisher:
IOP PublishingCopyright Statement
Country of Publication:
United States
Language:
English

References (79)

Conditional Neural Process for nonparametric modeling of active galactic nuclei light curves journal December 2021
Reverberation mapping of active galactic nuclei: From X-ray corona to dusty torus journal June 2021
Investigation of the Radio Source 3C 273 By The Method of Lunar Occultations journal March 1963
3C 273 : A Star-Like Object with Large Red-Shift journal March 1963
Absolute Energy Distribution in the Optical Spectrum of 3C 273 journal March 1963
Red-Shift of the Unusual Radio Source: 3C 48 journal March 1963
A luminous quasar at a redshift of z = 7.085 journal June 2011
An 800-million-solar-mass black hole in a significantly neutral Universe at a redshift of 7.5 journal December 2017
Array programming with NumPy journal September 2020
SciPy 1.0: fundamental algorithms for scientific computing in Python journal February 2020
Long-term optical variability properties of the Palomar-Green quasars journal July 1999
Pan-STARRS1 variability of XMM-COSMOS AGN: I. Impact on photometric redshifts⋆⋆⋆ journal December 2015
The spatially resolved broad line region of IRAS 09149−6206 journal November 2020
Estimating the size of X-ray lamppost coronae in active galactic nuclei journal December 2020
HOLISMOKES: III. Achromatic phase of strongly lensed Type Ia supernovae journal February 2021
Expectations for time-delay measurements in active galactic nuclei with the Vera Rubin Observatory journal July 2023
Optical Identification of 3c 48, 3c 196, and 3c 286 with Stellar Objects. journal July 1963
Accretion of Interstellar Matter by Massive Objects. journal August 1964
Reverberation mapping of the emission line regions of Seyfert galaxies and quasars journal April 1982
Slim accretion disks journal September 1988
Statistics of QSO Broad Emission-Line Profiles. I. The C IV lambda 1549 Line and the lambda 1400 Feature journal October 1993
Magnetic Stress at the Marginally Stable Orbit: Altered Disk Structure, Radiation, and Black Hole Spin Evolution journal January 2000
The Ensemble Photometric Variability of ∼25,000 Quasars in the Sloan Digital Sky Survey journal February 2004
Lag‐Luminosity Relationship for Interband Lags between Variations in B , V , R , and I Bands in Active Galactic Nuclei journal March 2005
Are the Variations in Quasar Optical flux Driven by Thermal Fluctuations? journal May 2009
Modeling the time Variability of sdss Stripe 82 Quasars as a Damped Random walk journal September 2010
The X-Ray Power Spectral Density Function of the Seyfert Active Galactic Nucleus ngc 7469 journal October 2010
Flexible and Scalable Methods for Quantifying Stochastic Variability in the era of Massive Time-Domain Astronomical data sets journal May 2014
SPACE TELESCOPE AND OPTICAL REVERBERATION MAPPING PROJECT. II. SWIFT AND HST REVERBERATION MAPPING OF THE ACCRETION DISK OF NGC 5548 journal June 2015
Is Quasar Optical Variability a Damped Random Walk? journal February 2013
Using quasar X-ray and UV flux measurements to constrain cosmological model parameters journal June 2020
Intensive disc-reverberation mapping of Fairall 9: first year of Swift  and LCO monitoring journal September 2020
Accretion disc sizes from continuum reverberation mapping of AGN selected from the ZTF survey journal January 2022
Optical variability of quasars with 20-yr photometric light curves journal May 2022
Legacy Survey of Space and Time cadence strategy evaluations for active galactic nucleus time-series data in Wide-Fast-Deep field journal March 2022
Survey of gravitationally lensed objects in HSC imaging (SuGOHI) – IX. Discovery of strongly lensed quasar candidates journal September 2023
Accretion disc time lag distributions: applying CREAM to simulated AGN light curves journal December 2015
X-ray/UV/optical variability of NGC 4593 with Swift: reprocessing of X-rays by an extended reprocessor journal August 2018
Winds can ‘blow up’ AGN accretion disc sizes journal October 2018
Reconciling the quasar microlensing disc size problem with a wind model of active galactic nucleus journal November 2018
A systematic study of the condensation of the corona and the application for Γ 2-10 keV-Lbol/LEdd correlation in luminous active galactic nuclei journal March 2018
Optical continuum photometric reverberation mapping of the Seyfert-1 galaxy Mrk509 journal October 2019
Quasar microlensing light-curve analysis using deep machine learning journal April 2019
Constraints on accretion disk size in the massive type 1 quasar PG 2308+098 from optical continuum reverberation lags journal September 2018
Gravitational Field of a Spinning Mass as an Example of Algebraically Special Metrics journal September 1963
Matplotlib: A 2D Graphics Environment journal January 2007
Combined long and short time-scale X-ray variability of NGC 4051 with RXTE and XMM-Newton journal March 2004
The relationship between X-ray variability amplitude and black hole mass in active galactic nuclei journal April 2005
X-ray variability of NGC 3227 and 5506 and the nature of active galactic nucleus 'states' journal October 2005
Testing thermal reprocessing in active galactic nuclei accretion discs journal August 2007
Gaussian Process Regression for Astronomical Time Series journal August 2023
corner.py: Scatterplot matrices in Python journal June 2016
Current and Future Applications of Reverberation-Mapped Quasars in Cosmology journal December 2019
Space Telescope and Optical Reverberation Mapping Project. iii. Optical Continuum Emission and Broadband time Delays in ngc 5548 journal April 2016
Testing wind as an Explanation for the spin Problem in the Continuum-Fitting Method journal April 2016
Fast and Scalable Gaussian Process Modeling with Applications to Astronomical Time Series journal November 2017
The Astropy Project: Building an Open-science Project and Status of the v2.0 Core Package journal August 2018
Modeling Stochastic Variability in Multiband Time-series Data journal November 2020
Searching for Changing-state AGNs in Massive Data Sets. I. Applying Deep Learning and Anomaly-detection Techniques to Find AGNs with Anomalous Variability Behaviors journal October 2021
Space Telescope and Optical Reverberation Mapping Project.Vi. Reverberating disk Models for ngc 5548 journal January 2017
Detection of Time Lags between Quasar Continuum Emission Bands Based On Pan-STARRS Light Curves journal February 2017
Swift Monitoring of NGC 4151: Evidence for a Second X-Ray/UV Reprocessing journal May 2017
Continuum Reverberation Mapping of the Accretion Disks in Two Seyfert 1 Galaxies journal February 2018
Quasar Accretion Disk Sizes from Continuum Reverberation Mapping from the Dark Energy Survey journal July 2018
The First Swift Intensive AGN Accretion Disk Reverberation Mapping Survey journal January 2019
LSST: From Science Drivers to Reference Design and Anticipated Data Products journal March 2019
The Sloan Digital Sky Survey Reverberation Mapping Project: Accretion Disk Sizes from Continuum Lags journal August 2019
Quasar Microlensing Variability Studies Favor Shallow Accretion Disk Temperature Profiles journal June 2020
Deep Modeling of Quasar Variability journal November 2020
Improving Damped Random Walk Parameters for SDSS Stripe 82 Quasars with Pan-STARRS1 journal February 2021
Modeling the Multiwavelength Variability of Mrk 335 Using Gaussian Processes journal June 2021
Examining AGN UV/Optical Variability beyond the Simple Damped Random Walk journal September 2022
Exploring the Link between the X-Ray Power Spectra and Energy Spectra of Active Galactic Nuclei journal August 2022
Continuum Reverberation Mapping of Mrk 876 over Three Years with Remote Robotic Observatories journal August 2023
Quasar Accretion Disk Sizes from Continuum Reverberation Mapping in the DES Standard-star Fields journal January 2020
The 16th Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra journal June 2020
The LSST Cadence Impact on Non-time-critical Eclipsing Binary Science journal July 2023
Estimating AGN Black Hole Masses via Continuum Reverberation Mapping in the Era of LSST journal May 2023
Gaussian Processes for Machine Learning book January 2005