AGN Variability with Rubin Observatory in the 2030s

TL;DR

AGN variability provides a direct probe of accretion physics and BLR structure but has been limited by sample size and cadence systematics; this work synthesizes ZTF results with LSST v5.1.0 cadence simulations to forecast continuum and emission-line reverberation mapping at population scale. It emphasizes that optical variability amplitudes and lags track accretion state, with $R_{disk}$ inferred from continuum reverberation sometimes inflated by diffuse BLR continuum emission, and shows that the LSST ocean cadence enhances lag recoverability across a broad $M_{BH}$ and redshift range, aided by standardized variability metrics and alert-broker integration. A two-tier observational strategy—LSST broadband monitoring plus targeted medium-band follow-up—along with multi-wavelength alerts and community-driven data products, is proposed to robustly constrain $R_{disk}$, BLR geometry, and SMBH growth across cosmic time. Fully realizing LSST's potential thus hinges on integrating variability-based physics with scalable data infrastructure, enabling transformative insights into accretion disks and black hole demographics.

Abstract

AGN variability offers a direct probe of accretion physics, disk structure, and black hole growth, but progress has been limited by sample size, cadence heterogeneity, and photometric systematics. The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will deliver multi-band light curves for millions of AGN, enabling variability studies at a true population scale. We synthesize recent results from the Zwicky Transient Facility (ZTF), which demonstrate that optical variability amplitudes and timescales are primarily regulated by accretion state, with secondary dependence on black hole mass and redshift, and establish the feasibility of survey-driven continuum reverberation mapping. ZTF measurements reveal optical continuum-emitting region sizes that often exceed standard thin disk predictions, implicating diffuse continuum emission from the broad line region as a significant contributor to observed inter-band lags. We evaluate the implications of LSST cadence and survey strategy, particularly the deep drilling fields, for continuum and emission line reverberation mapping, changing-look AGN, extreme variability quasars, and periodic variability searches. Key limitations of broadband photometric variability are identified, including variable emission line contamination, diffuse BLR continuum emission, and cadence-dependent lag recoverability. We argue that realizing LSST's full scientific potential requires community-scale, standardized variability metric pipelines, probabilistic classification integrated with alert brokers for follow-up triggering, and complementary medium-band photometric observations to isolate the accretion disk continuum. Together, these elements will enable LSST to convert photometric variability into quantitative constraints on accretion disks, BLR structure, and supermassive black hole growth across cosmic time.

Figures (10)

• Figure 1: Top panel: Survey footprint for the v5.1.0 cadence simulations. The auxiliary axis shows the number of planned visits over the 10-year baseline. Bottom panel: The name and location of the six deep-drilling fields (DDFs) highlighted in the top panel.
• Figure 2: 5$\sigma$ depth as a function of time for a randomly selected location (RA = 0$^{\circ}$, Dec = -30$^{\circ}$) in the Wide-Fast-Deep (WFD) LSST footprint for the v5.1.0 baseline simulations over the 10 years. The six filters (ugrizy) are shown using different colors. The legend highlights the number of visits corresponding to each filter.
• Figure 3: Similar to Figure \ref{['fig:wfd_lightcurve']}, but for COSMOS centered at RA = 150.11$^{\circ}$, Dec = 2.23$^{\circ}$ one of the six deep-drilling fields (DDFs). The legend highlights the number of visits and the standard deviation in the 5$\sigma$-depths corresponding to each filter. Notice the increase in the number of visits in the 2$^{\rm nd}$, 3$^{\rm rd}$, and 4$^{\rm th}$ year of operations in this field because of the adoption of the new "ocean" strategy.
• Figure 4: Cumulative number of visits in each deep-drilling field (DDF) based on the v5.1.0 baseline simulations over the 10 years. Notice the different season(s) of boosting in the number of visits for the DDFs. Careful optimization of science cases requiring dense cadence in each field is suggested.
• Figure 5: The distribution of QSO counts as a function of observed magnitude (AB) in different redshift ranges and filters, from the v5.1.0 simulations based on Li_2025arXiv251208654L. Each column represents a specific filter ($u$, $g$, $r$, $i$, $z$, and $y$ from left to right), and each row corresponds to a redshift range. The blue bins indicate the number of QSOs detected in each magnitude bin. Adapted from Figure 5 in their paper.