Gordon T Richards

The Legacy Survey of Space and Time (LSST), being conducted by the Vera C. Rubin Observatory, is a wide-field multi-band survey that will revolutionize our understanding of extragalactic sources through its unprecedented combination of area and depth. While the LSST survey strategy is still being finalized, the Rubin Observatory team has generated a series of survey simulations using the LSST Operations Simulator to explore the optimal survey strategy that best accommodates the majority of scientific goals. In this study, we utilize the latest simulated data to predict the number of detectable quasars by LSST in each band and evaluate the impact of different survey strategies. We find that the number of quasars and lower luminosity AGNs detected in the baseline strategy (v4.3.1) in the redshift range z=0.3-6.7 will be highest in the i-band and lowest in the u-band. Over 70% of quasars are expected to be detected within the first year in all bands, as LSST will have already reached the break of the luminosity function at most redshifts. With a limiting magnitude of 25.7 mag, we expect to detect 184 million AGNs in the z-band over the 10-year survey, with quasars constituting only 6% of the total AGNs in each band. This arises because, considering that the luminosities of most low-luminosity AGNs are affected by contamination from their host galaxies, we set a magnitude threshold when predicting the number of quasars. We find that variations in the u-band strategy can impact the number of quasar detections. Specifically, the difference between the baseline strategy and that with the largest total exposure in u is 15%. In contrast, changes in rolling strategies, DDF strategies, weather conditions, and Target of Opportunity observations result in variations below 2%. These results provide valuable insights for optimizing approaches to maximize the scientific output of quasar studies.

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.

The Vera C. Rubin Observatory LSST is expected to discover tens of millions of new Active Galactic Nuclei (AGNs). The survey's exceptional cadence and sensitivity will enable UV/optical/NIR monitoring of a significant fraction of these objects. The unprecedented number of sources makes spectroscopic follow-up for the vast majority of them unfeasible in the near future, so most studies will have to rely on photometric redshifts estimates which are traditionally much less reliable for AGN than for inactive galaxies. This work presents a novel methodology to constrain the photometric redshift of AGNs that leverages the effects of cosmological time dilation, and of the luminosity and wavelength dependence of AGN variability. Specifically, we assume that the variability can be modeled as a damped random walk (DRW) process, and adopt a parametric model to characterize the DRW timescale ($τ$) and asymptotic amplitude of the variability (SF$_\infty$) based on the redshift, the rest-frame wavelength, and the AGN luminosity. We construct variability-based photo-$z$ priors by modeling the observed variability using the expected DRW parameters at a given redshift. These variability-based photometric redshift (VAR-PZ) priors are then combined with traditional SED fitting to improve the redshift estimates from SED fitting.
Validation is performed using observational data from the SDSS, demonstrating significant reduction in catastrophic outliers by more than 10% in comparison with SED fitting techniques and improvements in redshift precision. The simulated light curves with both SDSS and LSST-like cadences and baselines confirm that, VAR-PZ will be able to constrain the photometric redshifts of SDSS-like AGNs by bringing the outlier fractions down to below 7% from 32% (SED-alone) at the end of the survey.

We present the first infrared spectral predictions from a self-consistent simulation of the formation of a quasar in a starburst galaxy, spanning cosmological to innermost stable circular orbit (ISCO) scales. The infrared emission is dominated by a torus-like dust structure composed of the highly magnetized, turbulence-supported outer accretion disk and of accreting gas tidally torn from the interstellar medium (ISM). At these early stages, the AGN is buried and Compton-thick. The near- to mid-IR escaping luminosity varies by almost an order of magnitude across sightlines, largely due to extinction from the inflowing stream of cold dust. Self-absorption within the torus suppresses silicate emission features, and further reprocessing by the ambient ISM leads to prominent silicate absorption and colder IR emission. The sublimation structure is stratified by composition and size, producing sightline-dependent extinction curves that intrinsically vary in shape. However, after repeated scattering in the optically thick dusty medium, these curves emerge substantially grayed. We also demonstrate that bipolar outflows from the central black hole that carve biconical cavities and reveal the central engine in later stages can preserve IR anisotropy and silicate features. These results suggest that dusty starburst quasars can undergo a buried, IR-bright phase early in their evolution.

While C IV is the most common absorption line in Broad Absorption Line Quasar spectra, Balmer absorption lines are among the rarest. We present analysis of Balmer absorption in a sample of fourteen iron low-ionization BAL quasars (FeLoBALQs); eight are new identifications. We measured velocity offset, width, and apparent optical depth. The partial covering ubiquitous in BAL quasar spectra alters the measured Balmer optical depth ratios; taking that into account, we estimated the true H(n= 2) column density. We found the anticipated correlation between Eddington ratio and outflow speed, but it is weak in this sample because nearly all of the objects have the low outflow speeds characterizing loitering outflow FeLoBAL quasars (H. Choi et al. 2022b), objects that are also found to have low accretion rates (K. M. Leighly et al. 2022; H. Choi et al. 2022a). Measures of dN/dv, the differential column density with respect to the outflow speed, are anticorrelated with the luminosity and Eddington ratio: the strongest absorption is observed at the lowest speeds in the lowest luminosity objects. The absorption line width is correlated with αoi, the Fλ point-to-point slope between 5100A and 3 microns. This parameter is strongly correlated with the Eddington ratio among low-redshift quasars (K. M. Leighly et al. 2024). Balmer absorption lines have been recently found in the spectra of Little Red Dots (LRDs), a class of high-redshift objects discovered by JWST. We note suggestive similarities between LRDs and FeLoBAL quasars in the emission line shape, the presence of steep reddening and a scattered blue continuum, the lack of hot dust emission, and X-ray weakness.

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.

We present a multiwavelength spectroscopic survey of 23 luminous mid-infrared-selected Type-2 quasars at redshifts z = 0.88 to 3.49. The targets were selected in the SDSS Stripe 82 field based on their bright WISE W4 detections (flux > 5 mJy) and extremely faint or red optical counterparts (e.g., r > 23 or r - W4 > 8.4), designed to identify heavily obscured quasars. Deep near-infrared (Gemini/GNIRS) and optical (Keck/LRIS and KCWI) spectroscopy confirm 23 out of 24 candidates as Type-2 quasars in this redshift range, including 12 objects at z > 2. The spectra exhibit strong rest-frame UV and optical emission lines (Ly-alpha, C IV, [O III], H-alpha) with a wide range of line widths, indicating significant spectral diversity. Approximately one-third of the sample (8 of 23) shows broad H-alpha emission (FWHM > 2000 km/s) despite their Type-2 classification, while the rest have only narrow lines (FWHM < 2000 km/s) characteristic of classical obscured quasars. Notably, these broad-line Type-2 quasars share similar spectral energy distributions with the JWST-discovered "little red dot" (LRD) AGNs, suggesting that our sample could be lower-redshift analogues of the heavily obscured broad-line AGNs uncovered by JWST. We also find that the [O III] 5007 angstrom emission is relatively weak for their high bolometric luminosities, deviating from trends seen in lower-redshift Type-2 QSOs. A new composite spectrum for Type-2 QSOs is built using our sample. Overall, our results demonstrate that mid-IR selection efficiently uncovers a diverse population of obscured quasars and that spectroscopic follow-up is crucial for revealing their true nature. This study provides new insights into heavily obscured SMBH growth at cosmic noon and bridges the gap to the obscured AGN populations now being revealed by JWST.

We present Chandra observations of 63 sources from the Gemini Near Infrared Spectrograph$-$Distant Quasar Survey (GNIRS-DQS) of which 54 were targeted by snapshot observations in Cycle 24. A total of 55 sources are clearly detected in at least one X-ray band, and we set stringent upper limits on the X-ray fluxes of the remaining eight sources. In combination with rest-frame ultraviolet-optical spectroscopic data for these sources, we assess whether X-rays can provide a robust accretion-rate indicator for quasars, particularly at the highest accessible redshifts. We utilize a recently modified H$β$-based Eddington luminosity ratio estimator, as well as the C IV $λ$1549 emission-line parameter space to investigate trends and correlations with the optical-X-ray spectral slope ($α_{\mathrm{ox}}$) and the effective hard-X-ray power-law photon index ($Γ$). We find that $α_{\mathrm{ox}}$ does not improve current accretion-rate estimates based on H$β$ or C IV. Instead, within the limitations of our sample, we confirm previous findings that the C IV parameter space may be a better indicator of the accretion rate up to $z\sim3.5$. We also find that the average $Γ$ values for a small subset of our sources, as well as the average $Γ$ value in different groupings of our sources, are consistent with their respective relatively high Eddington luminosity ratios. Deeper X-ray observations of our X-ray-detected sources are needed for measuring $Γ$ accurately and testing whether this parameter can serve as a robust, un-biased accretion-rate diagnostic.

We present the multiwavelength spectral energy distributions (SEDs) for 65 luminous broad absorption line (BAL) quasars with redshifts $1.55 \lesssim z \lesssim 3.50$ from the Gemini Near Infrared Spectrograph - Distant Quasar Survey (GNIRS-DQS). We integrate data from a variety of ground- and space-based observatories to construct a comprehensive spectral profile of these objects from radio through X-rays. In addition, we present a mid-infrared to X-ray composite SED of these sources. Our dataset represents the most uniform sample of BAL quasars, providing a statistically robust set of SEDs. Our findings indicate that the BAL quasars in the GNIRS-DQS sample exhibit significant reddening in the ultraviolet-optical continuum relative to their non-BAL counterparts, consistent with previous studies. Notably, our analysis reveals no significant differences in the mid- or near-infrared spectral regime between BAL and non-BAL quasars. In line with previous work, we find no strong evidence that BAL and non-BAL quasars possess fundamentally different SEDs, also consistent with recent findings that both groups display similar rest-frame optical emission-line properties.

We analyze 23 spectroscopically confirmed Type-2 quasars (QSOs) selected from the WISE 22$\rm \mu$m band in the SDSS Stripe 82 region, focusing on their multi-band photometry and spectral energy distributions (SEDs). These objects were selected to be IR-luminous ($\rm flux_{W4} > 5mJy$, i.e., $12.62 < W4 < 14.62 \rm\ AB \, magnitude$), optically faint ($r > 23$) or with red color ($r - W4 >8.38$). Gemini/GNIRS observations were conducted for all 24 candidates, and 18/24 were also observed with Keck/LRIS. The observations confirm 23 to be real Type-2 QSOs in the redshift range $0.88 - 2.99$ (12 are at $z>2$). We collect multi-band photometry and conduct SED fitting. The composite photometry probes the wavelength from 0.1$\rm \mu$m to 10$\rm \mu$m at the rest frame. The IR emission is dominated by dust torus implying an average torus luminosity for the sample of $L_{\rm torus} 10^{46.84} \rm erg/s$. The origin of the rest-UV/optical light is not definitive, but we present three possible scenarios: scattered light, stellar emission, and the reddened accretion disk. Assuming an obscured:unobscured ratio of approximately 1:1, our targets have $L_{\rm bol} = 10^{46.28} \rm erg \,s^{-1} - 10^{47.49} \rm erg \,s^{-1}$ and around SMBH masses $\rm 10^{8.18} M_{\odot} - 10^{9.39} M_{\odot}$, assuming they accreate at the Eddington limit. Compared to previous Type-2 AGN SEDs, our targets have a brighter dust torus and redder optical-IR color. By comparing the SED to the results from JWST `little red dots' (LRDs), we find that these IR-selected Type-2 QSOs have similar SED shapes to the LRDs. This pilot Type-2 QSO survey demonstrates that mid-IR selection is an efficient way to find luminous Type-2 QSOs at $z>2$. Finally, the composite photometry and Type-2 QSOs SED model generated by this sample provide a guide for finding more Type-2 QSOs at higher redshift.