Scalable and Robust Multiband Modeling of AGN Light Curves in Rubin-LSST

Abstract

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor tens of millions of active galactic nuclei (AGNs) for a period of 10 years with an average cadence of 3 days in six broad photometric bands. This unprecedented dataset will enable robust characterizations of AGN UV/optical variability across a wide range of AGN physical properties. However, existing tools for modeling AGN light curves are not yet capable of fully leveraging the volume, cadence, and multiband nature of LSST data. We present EzTaoX, a scalable light curve modeling tool designed to take advantage of LSST's multiband observations to simultaneously characterize AGN UV/optical stochastic variability and measure interband time delays. EzTaoX achieves a speed increase of $\sim 10^2-10^4\times$ on CPUs over current tools with similar capabilities, while maintaining equal or better accuracy in recovering simulated variability properties. This performance gain enables continuum time-delay measurements for all AGNs discovered by LSST -- both in the Wide Fast Deep survey and the Deep Drilling Fields -- thereby opening new opportunities to probe AGN accretion-flow geometries. In addition, EzTaoX's multiband capability allows robust characterization of AGN stochastic variability down to hourly timescales, facilitating the identification of accreting low-mass AGNs -- such as those residing in dwarf galaxies -- through their distinctive variability signatures.

arXiv commentAccepted for publication in ApJ, 24 pages, 15 figures

References in corpus (104)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• emcee: The MCMC Hammer
• The Sloan Digital Sky Survey: Technical Summary
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• The Zwicky Transient Facility: System Overview, Performance, and First Results
• Rank-normalization, folding, and localization: An improved $\hat{R}$ for assessing convergence of MCMC
• The Zwicky Transient Facility: Data Processing, Products, and Archive
• On characterising the variability properties of X-ray light curves from active galaxies
• Fast and scalable Gaussian process modeling with applications to astronomical time series
• The SDSS Quasar Survey: Quasar Luminosity Function from Data Release Three
• The Zwicky Transient Facility: Science Objectives
• Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations?
• Intermediate-Mass Black Holes
• Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk
• On Uncertainties in Cross-Correlation Lags and the Reality of Wavelength-Dependent Continuum Lags in Active Galactic Nuclei
• The Ensemble Photometric Variability of ~25000 Quasars in the Sloan Digital Sky Survey
• The Quasar Accretion Disk Size - Black Hole Mass Relation
• The Quasar Accretion Disk Size -- Black Hole Mass Relation
• Quantifying Quasar Variability As Part of a General Approach To Classifying Continuously Varying Sources
• An Alternative Approach To Measuring Reverberation Lags in Active Galactic Nuclei
• Is Quasar Optical Variability a Damped Random Walk?
• The flux-dependent amplitude of broadband noise variability in X-ray binaries and active galaxies
• Space Telescope and Optical Reverberation Mapping Project. III. Optical Continuum Emission and Broad-Band Time Delays in NGC 5548
• A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging
• 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
• Explaining changing-look AGN with state transition triggered by rapid mass accretion rate drop
• Kepler Observations of Rapid Optical Variability in Active Galactic Nuclei
• The Sizes of the X-ray and Optical Emission Regions of RXJ1131-1231
• Colors of 2625 Quasars at 0<z<5 Measured in the Sloan Digital Sky Survey Photometric System
• Reverberation mapping of Active Galactic Nuclei: from X-ray corona to dusty torus
• Quasar Accretion Disks Are Strongly Inhomogeneous
• X-Ray and Optical Flux Ratio Anomalies in Quadruply Lensed Quasars: I. Zooming in on Quasar Emission Regions
• Space Telescope and Optical Reverberation Mapping Project. I. Ultraviolet Observations of the Seyfert 1 Galaxy NGC 5548 with the Cosmic Origins Spectrograph on Hubble Space Telescope
• A characteristic optical variability timescale in astrophysical accretion disks
• Accretion disk reverberation with Hubble Space Telescope observations of NGC 4593: evidence for diffuse continuum lags
• Towards an Understanding of Changing-Look Quasars: An Archival Spectroscopic Search in SDSS
• The Variable Diffuse Continuum Emission of Broad Line Clouds
• Optimal Time-Series Selection of Quasars
• Revisiting Stochastic Variability of AGNs with Structure Functions
• Are the Variability Properties of the Kepler AGN Light Curves Consistent with a Damped Random Walk?
• Alert Classification for the ALeRCE Broker System: The Light Curve Classifier
• X-ray/UV/optical variability of NGC 4593 with Swift: Reprocessing of X-rays by an extended reprocessor
• The Kepler Light Curves of AGN: A Detailed Analysis
• Understanding extreme quasar optical variability with CRTS: II. Changing-state quasars
• A new direct method for measuring the Hubble constant from reverberating accretion discs in active galaxies
• Identifying AGNs in low-mass galaxies via long-term optical variability
• Limitations on the recovery of the true AGN variability parameters using Damped Random Walk modeling
• X-ray and Ultraviolet Properties of AGN in Nearby Dwarf Galaxies
• Quasar Selection Based on Photometric Variability
• Quantifying the Diffuse Continuum Contribution of BLR Clouds to AGN Continuum Inter-band Delays
• Improving Damped Random Walk parameters for SDSS Stripe 82 Quasars with Pan-STARRS1
• The Color Variability of Quasars
• Quasar Accretion Disk Sizes From Continuum Reverberation Mapping From the Dark Energy Survey
• Modeling the UV/optical continuum time-lags in AGN
• The Average Size and Temperature Profile of Quasar Accretion Disks
• 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
• AGN Variability Analysis Handbook
• Modelling type 1 quasar colours in the era of Rubin and Euclid
• Examining AGN UV/optical Variability Beyond the Simple Damped Random Walk
• Lognormal variability in BL Lacertae
• Quasar Accretion Disk Sizes from Continuum Reverberation Mapping in the DES Standard Star Fields
• Extracting Information from AGN Variability
• The QUEST-La Silla AGN variability survey: connection between AGN variability and black hole physical properties
• Accretion Disk Sizes from Continuum Reverberation Mapping of AGN Selected from the ZTF Survey
• AGN Continuum Reverberation Mapping Based on Zwicky Transient Facility Light Curves
• The dual nature of blazar fast variability. Space and ground observations of S5 0716+714
• A Deep Learning Approach for Active Anomaly Detection of Extragalactic Transients
• On Reverberation Mapping Lag Uncertainties
• Variability-selected intermediate mass black hole candidates in dwarf galaxies from ZTF and WISE
• The complex variability of blazars: Time-scales and periodicity analysis in S4 0954+65
• Dwarf AGNs from Optical Variability for the Origins of Seeds (DAVOS): Insights from the Dark Energy Survey Deep Fields
• A Non-parametric Approach to Constrain the Transfer Function in Reverberation Mapping
• Bayesian Analysis of Quasar Lightcurves with a Running Optimal Average: New Time Delay Measurements of COSMOGRAIL Gravitationally Lensed Quasars
• A structure function analysis of VST-COSMOS AGN
• Low frequency X-ray timing with Gaussian processes and reverberation in the radio-loud AGN 3C 120
• A Microlensing Accretion Disk Size Measurement in the Lensed Quasar WFI 2026-4536
• The LSST era of supermassive black holes accretion-disk reverberation mapping
• A Fast, 2D Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization
• A random forest-based selection of optically variable AGN in the VST-COSMOS field
• Gaussian Process Modeling Blazar Multiwavelength Variability: Indirectly Resolving Jet Structure
• Modelling the Multiwavelength Variability of Mrk 335 using Gaussian Processes
• Persistent and occasional: searching for the variable population of the ZTF/4MOST sky using ZTF data release 11
• Testing Super-Eddington Accretion onto a Supermassive Black Hole: Reverberation Mapping of PG 1119+120
• Twisted quasar light curves: implications for continuum reverberation mapping of accretion disks
• Ensemble power spectral density of SDSS quasars in UV/optical bands
• How can the optical variation properties of active galactic nuclei be unbiasedly measured?
• Continuum optical-UV and X-ray variability of AGN: current results and future challenges
• Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties
• Dwarf AGNs from Variability for the Origins of Seeds (DAVOS): Optical Variability of Broad-line Dwarf AGNs from the Zwicky Transient Facility
• Using Gaussian Processes to detect AGN flares
• Examining AGN UV/Optical Variability Beyond the Simple Damped Random Walk. II. Insights from 22 Years Observations of SDSS, PS1 and ZTF
• Modelling variability power spectra of active galaxies from irregular time series
• The Success of Optical Variability in Uncovering AGNs in Low-stellar Mass Galaxies
• Application of the FRADO model of BLR formation to the Seyfert galaxy NGC 5548 and the first step toward determining the Hubble constant

Cited by in corpus (1)

• Seeing through the dust: Unraveling near-infrared variability in type 2 active galactic nuclei