Researcher profile

Gordon T. Richards

Gordon T. Richards contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2025arXiv

Early Stages of Dusty Tori: The First Infrared Spectra from a Highly Multiscale Quasar Simulation

We present the first infrared spectral predictions from a self-consistent simulation of the formation of a quasar in a starburst galaxy, spanning the cosmological environment to scales well below the dust sublimation region. The infrared (IR) 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 active galactic nuclei (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.

preprint2022arXiv

Can X-ray Observations Improve Optical-UV-based Accretion-Rate Estimates for Quasars?

Current estimates of the normalized accretion rates of quasars (L/L_Edd), rely on measuring the velocity widths of broad optical-UV emission lines (e.g., H$β$ and Mg II $\lambda2800$). However, such lines tend to be weak or inaccessible in the most distant quasars, leading to increasing uncertainty in L/L_Edd estimates at $z > 6$. Utilizing a carefully selected sample of 53 radio-quiet quasars that have H$β$ and C IV $\lambda1549$ spectroscopy as well as {\sl Chandra} coverage, we searched for a robust accretion-rate indicator for quasars, particularly at the highest-accessible redshifts ($z \sim 6-7$). Our analysis explored relationships between the H$β$-based L/L_Edd, the equivalent width (EW) of C IV, and the optical-to-X-ray spectral slope (a_ox). Our results show that EW(C IV) is the strongest indicator of the H$β$-based L/L_Edd parameter, consistent with previous studies, although significant scatter persists particularly for sources with weak C IV lines. We do not find evidence for the a_ox parameter improving this relation, and we do not find a significant correlation between a_ox and H$β$-based L/L_Edd. This absence of an improved relationship may reveal a limitation in our sample. X-ray observations of additional luminous sources, found at $z \gtrsim 1$, may allow us to mitigate the biases inherent in our archival sample and test whether X-ray data could improve L/L_Edd estimates. Furthermore, deeper X-ray observations of our sources may provide accurate measurements of the hard-X-ray power-law photon index ($Γ$), which is considered an unbiased L/L_Edd indicator. Correlations between EW(C IV) and a_ox with $Γ$-based L/L_Edd may yield a more robust prediction of a quasar normalized accretion rate.

preprint2022arXiv

Connecting Low- and High-Redshift Weak Emission-Line Quasars via HST Spectroscopy of Ly$α$ Emission

We present ultraviolet spectroscopy covering the Ly$α$ + N V complex of six candidate low-redshift ($0.9 < z < 1.5$) weak emission-line quasars (WLQs) based on observations with the Hubble Space Telescope. The original systematic searches for these puzzling Type 1 quasars with intrinsically weak broad emission lines revealed an $N \approx 100$ WLQ population from optical spectroscopy of high-redshift ($z > 3$) quasars, defined by a Ly$α$ + N V rest-frame equivalent width (EW) threshold $< 15.4$ Å. Identification of lower-redshift ($z < 3$) WLQ candidates, however, has relied primarily on optical spectroscopy of weak broad emission lines at longer rest-frame wavelengths. With these new observations expanding existing optical coverage into the ultraviolet, we explore unifying the low- and high-$z$ WLQ populations via EW[Ly$α$+NV]. Two objects in the sample unify with high-$z$ WLQs, three others appear consistent with the intermediate portion of the population connecting WLQs and normal quasars, and the final object is consistent with typical quasars. The expanded wavelength coverage improves the number of available line diagnostics for our individual targets, allowing a better understanding of the shapes of their ionizing continua. The ratio of EW[Ly$α$+NV] to EW[MgII] in our sample is generally small but varied, favoring a soft ionizing continuum scenario for WLQs, and we find a lack of correlation between EW[Ly$α$+NV] and the X-ray properties of our targets, consistent with a &#34;slim-disk&#34; shielding gas model. We also find indications that weak absorption may be a more significant contaminant in low-$z$ WLQ populations than previously thought.

preprint2022arXiv

Examining AGN UV/optical Variability Beyond the Simple Damped Random Walk

We present damped harmonic oscillator (DHO) light-curve modeling for a sample of 12,714 spectroscopically confirmed quasars in the Sloan Digital Sky Survey Stripe 82 region. DHO is a second-order continuous-time autoregressive moving-average (CARMA) process, which can be fully described using four independent parameters: a natural oscillation frequency ($ω_{0}$), a damping ratio ($ξ$), a characteristic perturbation timescale ($τ_{\mathrm{perturb}}$), and an amplitude for the perturbing white noise ($σ_{\mathrmε}$). The asymptotic variability amplitude of a DHO process is quantified by $σ_{\mathrm{DHO}}$ -- a function of $ω_{0}$, $ξ$, $τ_{\mathrm{perturb}}$, and $σ_{\mathrmε}$. We find that both $τ_{\mathrm{perturb}}$ and $σ_{\mathrmε}$ follow different dependencies with rest-frame wavelength ($λ_{\mathrm{RF}}$) on either side of 2500 Å, whereas $σ_{\mathrm{DHO}}$ follows a single power-law relation with $λ_{\mathrm{RF}}$. After correcting for wavelength dependence, $σ_{\mathrm{DHO}}$ exhibits anti-correlations with both the Eddington ratio and the black hole mass, while $τ_{\mathrm{perturb}}$ -- with a typical value of days in the rest-frame -- shows an anti-correlation with the bolometric luminosity. Modeling AGN variability as a DHO offers more insight into the workings of accretion disks close to the supermassive black holes (SMBHs) at the center of AGN. The newly discovered short-term variability (characterized by $τ_{\mathrm{perturb}}$ and $σ_{\mathrmε}$) and its correlation with bolometric luminosity pave the way for new algorithms that will derive fundamental properties (e.g., Eddington ratio) of AGN using photometric data alone.

preprint2022arXiv

The Physical Properties of Low-Redshift FeLoBAL Quasars: II. The Rest-Frame Optical Emission Line Properties

We report the results of analysis of the Hbeta emission-line region of a sample of thirty low-redshift (z<1) iron low-ionization broad absorption line quasars (FeLoBALQs). Eleven of these objects are newly classified as FeLoBALQs. A matched sample of 132 unabsorbed quasars was analyzed in parallel. The emission lines showed the well known anticorrelation between the [OIII] and FeII emission (Boroson & Green 1992). Using a summary statistic called E1 to quantify this anticorrelation, we found that while the distribution of E1 for the unabsorbed quasars has a single peak, the FeLoBALQs have a bimodal shape in this parameter. Previous studies have shown that the line emission properties of BAL and non-BAL quasars are consistent, and therefore the difference in the Hbeta region emission between FeLoBAL quasars and unabsorbed quasars is a new result. The two populations of FeLoBAL quasars are characterized by low and high bolometric luminosities and Eddington ratios. Some previous studies have suggested that BAL quasars are high accretion-rate objects, and therefore the discovery of the low accretion-rate branch of FeLoBAL quasars was unexpected. We also found that the Hbeta FWHM is systematically broader among the FeLoBALQs compared with the non-BAL quasars implying a higher inclination viewing angle or a dearth of low-velocity line-emitting gas.

preprint2021arXiv

Blazar variability with the Vera C. Rubin Legacy Survey of Space and Time (LSST)

With their emission mainly coming from a relativistic jet pointing towards us, blazars are fundamental sources to study extragalactic jets and their central engines, consisting of supermassive black holes (SMBHs) fed by accretion discs. They are also candidate sources of high-energy neutrinos and cosmic rays. Because of the jet orientation, the non-thermal blazar emission is Doppler beamed; its variability is unpredictable and occurs on time-scales from less than one hour to years. The comprehension of the diverse mechanisms producing the flux and spectral changes requires well-sampled multiband light curves on long time periods. In particular, outbursts are the best test bench to shed light on the underlying physics, especially when studied in a multiwavelength context. The Vera C. Rubin Legacy Survey of Space and Time (Rubin-LSST) will monitor the southern sky for ten years in six photometric bands, offering a formidable tool to study blazar variability features in a statistical way. The alert system will allow us to trigger follow-up observations of outstanding events, especially at high (keV-to-GeV) and very high (TeV) energies. We here examine the simulated Rubin-LSST survey strategies with the aim of understanding which cadences are more suitable for the blazar variability science. Our metrics include light curve and colour sampling. We also investigate the problem of saturation, which will affect the brightest and many flaring sources, and will have a detrimental impact on follow-up observations.

preprint2021arXiv

Exploring the link between C IV outflow kinematics and sublimation-temperature dust in quasars

Using data from SDSS, UKIDSS and WISE, we investigate the properties of the high-frequency cutoff to the infrared emission in $\simeq$5000 carefully selected luminous ($L_{bol} \simeq 10^{47}$) type 1 quasars. The strength of $\simeq$2 $μ$m emission, corresponding to emission from the hottest (T>1200K) dust in the sublimation zone surrounding the central continuum source, is observed to correlate with the blueshift of the C IV $λ$1550 emission line. We therefore find that objects with stronger signatures of nuclear outflows tend to have a larger covering fraction of sublimation-temperature dust. When controlling for the observed outflow strength, the hot dust covering fraction does not vary significantly across our sample as a function of luminosity, black hole mass or Eddington fraction. The correlation between the hot dust and the C IV line blueshifts, together with the lack of correlation between the hot dust and other parameters, therefore provides evidence of a link between the properties of the broad emission line region and the infrared-emitting dusty regions in quasars.

preprint2021arXiv

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey&#39;s massive data throughput will be transformational for many other astrophysics domains and Rubin&#39;s data access policy sets the stage for a huge potential users&#39; community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

preprint2021arXiv

Placing LOFAR-detected quasars in CIV emission space: implications for winds, jets and star formation

We present an investigation of the low-frequency radio and ultraviolet properties of a sample of $\simeq$10,500 quasars from the Sloan Digital Sky Survey Data Release 14, observed as part of the first data release of the Low-Frequency-Array Two-metre Sky Survey. The quasars have redshifts $1.5 < z < 3.5$ and luminosities $44.6 < \log(L_{\text{bol}}/\text{erg s}^{-1}) < 47.2$. We employ ultraviolet spectral reconstructions based on an independent component analysis to parametrize the CIV$λ$1549-emission line that is used to infer the strength of accretion disc winds, and the HeII$λ$1640 line, an indicator of the soft X-ray flux. We find that radio-detected quasars are found in the same region of CIV blueshift versus equivalent-width space as radio-undetected quasars, but that the loudest, most luminous and largest radio sources exist preferentially at low CIV blueshifts. Additionally, the radio-detection fraction increases with blueshift whereas the radio-loud fraction decreases. In the radio-quiet population, we observe a range of HeII equivalent widths as well as a Baldwin effect with bolometric luminosity, whilst the radio-loud population has mostly strong HeII, consistent with a stronger soft X-ray flux. The presence of strong HeII is a necessary but not sufficient condition to detect radio-loud emission suggesting some degree of stochasticity in jet formation. Using energetic arguments and Monte Carlo simulations, we explore the plausibility of winds, compact jets and star formation as sources of the radio quiet emission, ruling out none. The existence of quasars with similar ultraviolet properties but differing radio properties suggests, perhaps, that the radio and ultraviolet emission is tracing activity occurring on different timescales.

preprint2020arXiv

BAL and non-BAL quasars: continuum, emission and absorption properties establish a common parent sample

Using a sample of $\simeq$144,000 quasars from the Sloan Digital Sky Survey data release 14 we investigate the outflow properties, evident both in absorption and emission, of high-ionization Broad Absorption Line (BAL) and non-BAL quasars with redshifts 1.6 $\lesssim z \leq$ 3.5 and luminosities 45.3 $< \log_{10}(L_{bol}) < $ 48.2 erg s$^{-1}$. Key to the investigation is a continuum and emission-line reconstruction scheme, based on mean-field independent component analysis, that allows the kinematic properties of the CIV$λ$1550 emission line to be compared directly for both non-BAL and BAL quasars. CIV-emission blueshift and equivalent-width (EW) measurements are thus available for both populations. Comparisons of the emission-line and BAL-trough properties reveal strong systematic correlations between the emission and absorption properties. The dependence of quantitative outflow indicators on physical properties such as quasar luminosity and luminosity relative to Eddington-luminosity are also shown to be essentially identical for the BAL and non-BAL populations. There is an absence of BALs in quasars with the hardest spectral energy distributions (SEDs), revealed by the presence of strong HeII$λ$1640 emission, large CIV$λ$1550-emission EW and no measurable blueshift. In the remainder of the CIV-emission blueshift versus EW space, BAL and non-BAL quasars are present at all locations; for every BAL-quasar it is possible to identify non-BAL quasars with the same emission-line outflow properties and SED-hardness. The co-location of BAL and non-BAL quasars as a function of emission-line outflow and physical properties is the key result of our investigation, demonstrating that (high-ionization) BALs and non-BALs represent different views of the same underlying quasar population.

preprint2020arXiv

Discovery of a Remarkably Powerful Broad Absorption Line Quasar Outflow in SDSS J135246.37+423923.5

Broad absorption line (BAL) features in quasar spectra reveal an unambiguous signature of energetic outflows from central supermassive black holes, and thus BAL quasars are prime targets for investigating the potential process of luminous quasar feedback on galaxies. We analyzed the rest-UV spectrum of an &#34;overlapping trough&#34; iron low-ionization broad absorption line quasar (FeLoBAL) SDSS J135246.37+423923.5 using the novel spectral synthesis code SimBAL (Leighly et al. 2018) and discovered an extraordinarily fast and energetic BAL outflow. Our analysis revealed outflow velocities reaching $\sim -38000\rm \, km\, s^{-1}$ with a velocity width of $\sim 10000\rm \, km\, s^{-1}$ which is the largest FeLoBAL outflow velocity measured to date. The column density of the outflow gas is log$N_H\sim23.2\,[\rm cm^{-1}]$ with the log kinetic luminosity $\log L_{KE}\sim48.1$ [erg $\rm s^{-1}$] which exceeds the bolometric luminosity of the quasar and is energetic enough to effectively drive quasar feedback. The energy estimate for the outflow is far greater than the estimates from any BAL object previously reported. The object also shows &#34;anomalous reddening&#34; and a significant scattered component that we were able to model with SimBAL. We found the first definitive case for radiation filtering in an additional zero-velocity absorption component that required an absorbed continuum to produce the particular absorption lines observed (MgII, AlIII and AlII) without also producing the high ionization lines such as CIV.

preprint2020arXiv

Placing High-Redshift Quasars in Perspective: a Catalog of Spectroscopic Properties from the Gemini Near Infrared Spectrograph -- Distant Quasar Survey

We present spectroscopic measurements for 226 sources from the Gemini Near Infrared Spectrograph - Distant Quasar Survey (GNIRS-DQS). Being the largest uniform, homogeneous survey of its kind, it represents a flux-limited sample ($m_{i}$ ${\lesssim}$ 19.0 mag, $H$ ${\lesssim}$ 16.5 mag) of Sloan Digital Sky Survey (SDSS) quasars at 1.5 ${\lesssim}$ $z$ ${\lesssim}$ 3.5 with a monochromatic luminosity ($λL_λ$) at 5100 ${\unicode{xC5}}$ in the range of $10^{44} - 10^{46}$ erg $\rm{s}^{-1}$. A combination of the GNIRS and SDSS spectra covers principal quasar diagnostic features, chiefly the C IV $λ$1549, Mg II $λλ$2798, 2803, H$β$ $λ$4861, and [O III] $λλ$4959, 5007 emission lines, in each source. The spectral inventory will be utilized primarily to develop prescriptions for obtaining more accurate and precise redshifts, black hole masses, and accretion rates for all quasars. Additionally, the measurements will facilitate an understanding of the dependence of rest-frame ultraviolet-optical spectral properties of quasars on redshift, luminosity, and Eddington ratio, and test whether the physical properties of the quasar central engine evolve over cosmic time.

preprint2020arXiv

The Bolometric Quasar Luminosity Function at z = 0-7

In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new best-fit bolometric quasar luminosity function behaves qualitatively different from the Hopkins et al. 2007 model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at $z\gtrsim 2$. Due to the paucity of measurements at the faint end, the faint end slope at $z\gtrsim 5$ is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at $z\gtrsim 5$. Further multi-band observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at $z\lesssim 3$. The predicted CXB spectrum, cosmic SMBH mass density and SMBH mass function are generally consistent with existing observations.