Michael S Vogeley

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.

Observations and theoretical simulations suggest that the large-scale environment plays a significant role in how galaxies form and evolve and, in particular, whether and when galaxies host an actively accreting supermassive black hole in their center (i.e., an active galactic nucleus; AGN). One signature of AGN activity is luminosity variability, which appears in the mid-IR when circumnuclear dust reprocesses UV and optical photons from the AGN accretion disk. We present here a suite of constraints on the fraction of AGN activity in the most underdense regions of the Universe (cosmic voids) relative to the rest of the Universe (cosmic walls) by using ∼12 yr of combined multiepoch data from AllWISE and NEOWISE to quantify mid-IR variability. We find clear evidence for a larger mid-IR variability−AGN fraction among high- and moderate-luminosity void galaxies compared to their wall counterparts. We also show that mid-IR variability identifies a rather large and unique population of AGNs, the majority of which have eluded detection using more traditional AGN selection methods such as single-epoch mid-IR color selection. The fraction of these newly recovered AGNs is larger among galaxies in voids, suggesting once again more prolific AGN activity in the most underdense large-scale structures of the Universe.

A damped random walk (DRW) process is often used to describe the temporal UV/optical continuum variability of active galactic nuclei (AGN). However, recent investigations have shown that this model fails to capture the full spectrum of AGN variability. In this work, we model the 22-year-long light curves of $21,767$ quasars, spanning the redshift range $0.28 < z < 2.71$, as a noise-driven damped harmonic oscillator (DHO) process. The light curves, in the optical $g$ and $r$ bands, are collected and combined from the Sloan Digital Sky Survey, the Panoramic Survey Telescope and Rapid Response System, and the Zwicky Transient Facility. A DHO process can be defined using four parameters, two for describing its long-term behavior/variability, and the other two for describing its short-term behavior/variability. We find that the best-fit DHO model describes the observed variability of our quasar light curves better than the best-fit DRW model. Furthermore, the best-fit DHO parameters exhibit correlations with the rest-frame wavelength, the Eddington ratio, and the black hole mass of our quasars. Based on the power spectral density shape of the best-fit DHOs and these correlations, we suggest that the observed long-term variability of our quasars can be best explained by accretion rate or thermal fluctuations originating from the accretion disk, and the observed short-term variability can be best explained by reprocessing of X-ray variability originating from the corona. The additional information revealed by DHO modeling emphasizes the need to go beyond DRW when analyzing AGN light curves delivered by next-generation wide-field time-domain surveys.

Observations and theoretical simulations suggest that the large scale environment plays a significant role in how galaxies form and evolve and, in particular, whether and when galaxies host an actively accreting supermassive black hole in their center (i.e., an Active Galactic Nucleus, or AGN). One signature of AGN activity is luminosity variability, which appears in the mid-infrared (mid-IR) when circumnuclear dust reprocesses UV and optical photons from the AGN accretion disk. We present here a suite of constraints on the fraction of AGN activity in the most underdense regions of the universe (cosmic voids) relative to the rest of the universe (cosmic walls) by using ~12 years of combined multi-epoch data from AllWISE and NEOWISE to quantify mid-IR variability. We find clear evidence for a larger mid-IR variability-AGN fraction among high and moderate-luminosity void galaxies compared to their wall counterparts. We also show that mid-IR variability identifies a rather large and unique population of AGNs, the majority of which have eluded detection using more traditional AGN-selection methods such as single-epoch mid-IR color selection. The fraction of these newly-recovered AGNs is larger among galaxies in voids, suggesting once again more prolific AGN activity in the most underdense large scale structures of the universe.

A damped random walk (DRW) process is often used to describe the temporal UV/optical continuum variability of active galactic nuclei (AGN). However, recent investigations have shown that this model fails to capture the full spectrum of AGN variability. In this work, we model the 22 yr long light curves of 21,767 quasars, spanning the redshift range 0.28 <  z  < 2.71, as a noise-driven damped harmonic oscillator (DHO) process. The light curves, in the optical g and r bands, are collected and combined from the Sloan Digital Sky Survey, the Panoramic Survey Telescope and Rapid Response System, and the Zwicky Transient Facility. A DHO process can be defined using four parameters, two for describing its long-term behavior/variability, and the other two for describing its short-term behavior/variability. We find that the best-fit DHO model describes the observed variability of our quasar light curves better than the best-fit DRW model. Furthermore, the best-fit DHO parameters exhibit correlations with the rest-frame wavelength, the Eddington ratio, and the black hole mass of our quasars. Based on the power spectral density shape of the best-fit DHOs and these correlations, we suggest that the observed long-term variability of our quasars can be best explained by accretion rate or thermal fluctuations originating from the accretion disk, and the observed short-term variability can be best explained by reprocessing of X-ray variability originating from the corona. The additional information revealed by DHO modeling emphasizes the need to go beyond DRW when analyzing AGN light curves delivered by next-generation wide-field time-domain surveys.

We investigate the spatial distribution of Lyman-$\alpha$ (Ly $\alpha$) absorbers within cosmic voids. We create a catalogue of cosmic voids in Sloan Digital Sky Survey Data Release 7 (SDSS DR7) with the VoidFinder algorithm of the Void Analysis Software Toolkit (VAST). Using the largest catalogue of low-redshift (z $\leq$ 0.75) IGM absorbers to date, we identify 392 Ly $\alpha$ absorbers inside voids. The fraction of Ly $\alpha$ absorbers inside voids (65 per cent) is comparable to the volume filling fraction of voids (68 per cent), and significantly greater than the fraction of galaxies inside voids (21 per cent). Inside voids, the spatial distribution of Ly $\alpha$ absorbers differs markedly from that of galaxies. Galaxy density rises sharply near void edges, while Ly $\alpha$ absorber density is relatively uniform. The radial distribution of Ly $\alpha$ absorbers inside voids differs marginally from a random distribution. We find that lower column density Ly $\alpha$ absorbers are more centrally concentrated inside voids than higher column density Ly $\alpha$ absorbers. These results suggest the presence of two populations of Ly $\alpha$ absorbers: low column density systems that are nearly uniformly distributed in the interiors of voids and systems associated with galaxies at the edges of voids.

The Kepler satellite potentially provides the highest precision photometry of active galactic nuclei (AGNs) available to investigate short-timescale optical variability. We targeted quasars from the Sloan Digital Sky Survey that lie in the fields of view of the Kepler/K2 campaigns. Based on those observations, we report the discovery and properties of a previously unidentified instrumental signature in K2. Systematic errors in K2, beyond those due to the motion of the detector, plague our AGNs and other faint-target, guest observer science proposals. Weakly illuminated pixels are dominated by low-frequency trends that are both nonastrophysical and correlated from object to object. The instrumental signature lags in time as a function of radius from the center of the detector, crossing channel boundaries. Thus, systematics documented in this investigation are unlikely to be due to Moire noise, rolling band, or pointing jitter. A critical clue to understanding this instrumental systematic is that different targets observed in the same channels of Campaign 8 (rear facing) and Campaign 16 (forward facing) have nearly identical light curves after time reversal of one of the campaigns. We find evidence of temperature trends that also reverse according to the Sun-spacecraft field orientation and that may dominate the systematics. These temperature variations are larger in K2 than in the nominal Kepler mission and strongly support our hypothesis of temperature-driven focus changes. Further characterization of this signature is crucial for rehabilitating K2 data for use in investigations of AGN light curves.

The $\textit{Kepler}$ satellite potentially provides the highest precision
photometry of active galactic nuclei (AGN) available to investigate
short-timescale optical variability. We targeted quasars from the Sloan Digital
Sky Survey that lie in the fields of view of the $\textit{Kepler/K2}$
campaigns. Based on those observations, we report the discovery and properties
of a previously unidentified instrumental signature in K2. Systematic errors in
K2, beyond those due to the motion of the detector, plague our AGN and other
faint-target, guest-observer science proposals. Weakly illuminated pixels are
dominated by low frequency trends that are both non-astrophysical and
correlated from object to object. A critical clue to understanding this
instrumental noise is that different targets observed in the same channels of
Campaign 8 (rear facing) and Campaign 16 (forward facing) had nearly identical
light curves after time reversal of one of the campaigns. This observation
strongly suggests that the underlying problem relates to the relative
Sun-spacecraft-field orientation, which was approximately the same on day 1 of
Campaign 8 as the last day of Campaign 16. Furthermore, we measure that the
instrumental signature lags in time as a function of radius from the center of
the detector, crossing channel boundaries. Systematics documented in this
investigation are unlikely to be due to Moir\'{e} noise, rolling band, or
pointing jitter. Instead this work strongly suggests temperature-dependent
focus changes that are further subject to channel variations. Further
characterization of this signature is crucial for rehabilitating K2 data for
use in investigations of AGN light curves.

The advent of new time domain surveys and the imminent increase in astronomical data expose the shortcomings of traditional time series analysis (such as power spectra analysis) in characterizing the abundantly varied, complex, and stochastic light curves of Active Galactic Nuclei (AGNs). Recent applications of novel methods from non-linear dynamics have shown promise in characterizing higher modes of variability and time-scales in AGN. Recurrence analysis in particular can provide complementary information about characteristic time-scales revealed by other methods, as well as probe the nature of the underlying physics in these objects. Recurrence analysis was developed to study dynamical trajectories in phase space, which can be constructed from 1D time series such as light curves. We apply the methods of recurrence analysis to two optical light curves of Kepler-monitored AGN. We confirm the detection and period of an optical quasi-periodic oscillation in one AGN, and confirm multiple other time-scales recovered from other methods ranging from 5 to 60 d in both objects. We detect regions in the light curves that deviate from regularity, provide evidence of determinism and non-linearity in the mechanisms underlying one light curve (KIC 9650712), and determine realizations of a linear stochastic process describe the dominant variability in the other light curve (Zwicky 229-015). We discuss possible underlying processes driving the dynamics of the light curves and their diverse classes of variability.