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

Abstract

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.

22 pages, 12 figures

References in corpus (58)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• A luminous quasar at a redshift of z = 7.085
• The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra
• An 800-million-solar-mass black hole in a significantly neutral Universe at redshift 7.5
• Fast and scalable Gaussian process modeling with applications to astronomical time series
• Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations?
• Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk
• The Ensemble Photometric Variability of ~25000 Quasars in the Sloan Digital Sky Survey
• Combined Long and Short Timescale X-ray Variability of NGC4051 with RXTE and XMM-Newton
• Is Quasar Optical Variability a Damped Random Walk?
• Space Telescope and Optical Reverberation Mapping Project. III. Optical Continuum Emission and Broad-Band Time Delays in NGC 5548
• Space Telescope and Optical Reverberation Mapping Project. II. Swift and HST Reverberation Mapping of the Accretion Disk of NGC 5548
• Testing thermal reprocessing in AGN accretion discs
• Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets
• Reverberation mapping of Active Galactic Nuclei: from X-ray corona to dusty torus
• The First Swift Intensive AGN Accretion Disk Reverberation Mapping Survey
• Swift monitoring of NGC 4151: Evidence for a Second X-ray/UV Reprocessing
• X-ray/UV/optical variability of NGC 4593 with Swift: Reprocessing of X-rays by an extended reprocessor
• The relationship between X-ray variability amplitude and black hole mass in active galactic nuclei
• Intensive disc-reverberation mapping of Fairall 9: 1st year of Swift & LCO monitoring
• Space Telescope and Optical Reverberation Mapping Project VI: reverberating Disk Models for NGC 5548
• Using quasar X-ray and UV flux measurements to constrain cosmological model parameters
• Improving Damped Random Walk parameters for SDSS Stripe 82 Quasars with Pan-STARRS1
• Quasar Accretion Disk Sizes From Continuum Reverberation Mapping From the Dark Energy Survey
• The spatially resolved broad line region of IRAS 09149-6206
• Accretion Disc Time Lag Distributions: Applying CREAM to Simulated AGN Light Curves
• Detection of Time Lags Between Quasar Continuum Emission Bands based on Pan-STARRS Light-curves
• The Sloan Digital Sky Survey Reverberation Mapping Project: Accretion-Disk Sizes from Continuum Lags
• Continuum Reverberation Mapping of the Accretion Disks in Two Seyfert 1 Galaxies
• Optical Variability of Quasars with 20-Year Photometric Light Curves
• Examining AGN UV/optical Variability Beyond the Simple Damped Random Walk
• Quasar Accretion Disk Sizes from Continuum Reverberation Mapping in the DES Standard Star Fields
• Accretion Disk Sizes from Continuum Reverberation Mapping of AGN Selected from the ZTF Survey
• Searching for changing-state AGNs in massive datasets -- I: applying deep learning and anomaly detection techniques to find AGNs with anomalous variability behaviours
• Winds Can "Blow Up" AGN Accretion Disk Sizes
• Current and future applications of Reverberation-mapped quasars in Cosmology
• Deep modeling of quasar variability
• Quasar Microlensing Variability Studies Favor Shallow Accretion Disk Temperature Profiles
• HOLISMOKES -- III. Achromatic Phase of Strongly Lensed Type Ia Supernovae
• Pan-STARRS1 variability of XMM-COSMOS AGN. I. Impact on photometric redshifts
• Reconciling the Quasar Microlensing Disc Size Problem with a Wind Model of Active Galactic Nucleus
• Quasar microlensing light curve analysis using deep machine learning
• Estimating the size of X-ray lamppost coronae in active galactic nuclei
• Testing wind as an explanation for the spin problem in the continuum-fitting method
• Modeling Stochastic Variability in Multi-Band Time Series Data
• Optical continuum photometric reverberation mapping of the Seyfert-1 galaxy Mrk509
• Modelling the Multiwavelength Variability of Mrk 335 using Gaussian Processes
• Estimating AGN Black Hole Masses via Continuum Reverberation Mapping in the Era of LSST
• A systematic study of the condensation of the corona and the application for ${\Gamma }_{2-10keV}-L_{\mathrm{bol}}/L_{\mathrm{Edd}}$ correlation in luminous active galactic nuclei
• Expectations for time-delay measurements in active galactic nuclei with the Vera Rubin Observatory
• The X-ray Power Spectral Density Function of the Seyfert Active Galactic Nucleus NGC 7469
• LSST Cadence Strategy Evaluations for AGN Time-series Data in Wide-Fast-Deep Field
• Constraints on accretion disk size in the massive type 1 quasar PG 2308+098 from optical continuum reverberation lags
• Conditional Neural Process for non-parametric modeling of AGN light curve
• Exploring the Link between the X-ray Power Spectra and Energy Spectra of Active Galactic Nuclei

Cited by in corpus (8)

• Insights from leptohadronic modelling of the brightest blazar flare
• Joint Modeling of Quasar Variability and Accretion Disk Reprocessing using Latent Stochastic Differential Equations
• Scalable and Robust Multiband Modeling of AGN Light Curves in Rubin-LSST
• Predicting High-magnification Events in Microlensed Quasars in the Era of LSST using Recurrent Neural Networks
• Amoeba: An AGN Model of Optical Emissions Beyond steady-state Accretion discs
• Modeling quasar variability through self-organizing map-based process
• Random walks around black holes and low-frequency X-ray variability
• Modeling Irregular Astronomical Time Series with Neural Stochastic Delay Differential Equations