Source author record

Gordon T. Richards

Gordon T. Richards appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

94works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

94 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 "slim-disk" 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's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' 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 "overlapping trough" 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 "anomalous reddening" 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.

preprint2016arXiv

CIV Emission Line Properties and Systematic Trends in Quasar Black Hole Mass Estimates

Black-hole masses are crucial to understanding the physics of the connection between quasars and their host galaxies and measuring cosmic black hole-growth. At high redshift, z > 2.1, black hole masses are normally derived using the velocity-width of the CIV broad emission line, based on the assumption that the observed velocity-widths arise from virial-induced motions. In many quasars, the CIV-emission line exhibits significant blue asymmetries (`blueshifts') with the line centroid displaced by up to thousands of km/s to the blue. These blueshifts almost certainly signal the presence of strong outflows, most likely originating in a disc wind. We have obtained near-infrared spectra, including the H$α$ emission line, for 19 luminous ($L_{Bol}$ = 46.5-47.5 erg/s) Sloan Digital Sky Survey quasars, at redshifts 2 < z < 2.7, with CIV emission lines spanning the full-range of blueshifts present in the population. A strong correlation between CIV-velocity width and blueshift is found and, at large blueshifts, > 2000 km/s, the velocity-widths appear to be dominated by non-virial motions. Black-hole masses, based on the full width at half maximum of the CIV-emission line, can be overestimated by a factor of five at large blueshifts. A larger sample of quasar spectra with both CIV and H$β$, or H$α$, emission lines will allow quantitative corrections to CIV-based black-hole masses as a function of blueshift to be derived. We find that quasars with large CIV blueshifts possess high Eddington luminosity ratios and that the fraction of high-blueshift quasars in a flux-limited sample is enhanced by a factor of approximately four relative to a sample limited by black hole mass.

preprint2016arXiv

Extremely Red Quasars from SDSS, BOSS and WISE: Classification of Optical Spectra

Quasars with extremely red infrared-to-optical colours are an interesting population that can test ideas about quasar evolution as well as orientation, obscuration and geometric effects in the so-called AGN unified model. To identify such a population we match the quasar catalogues of the Sloan Digital Sky Survey (SDSS), the Baryon Oscillation Spectroscopic Survey (BOSS) to the Wide-Field Infrared Survey Explorer (WISE) to identify quasars with extremely high infrared-to-optical ratios. We identify 65 objects with r(AB)-W4(Vega)>14 mag (i.e., F_nu(22um)/F_nu(r) > ~1000). This sample spans a redshift range of 0.28<z<4.36 and has a bimodal distribution, with peaks at z~0.8 and z~2.5. It includes three z>2.6 objects that are detected in the W4-band but not W1 or W2 (i.e., W1W2-dropouts). The SDSS/BOSS spectra show that the majority of the objects are reddened Type 1 quasars, Type 2 quasars (both at low and high redshift) or objects with deep low-ionization broad absorption lines (BALs) that suppress the observed r-band flux. In addition, we identify a class of Type 1 permitted broad-emission line objects at z~2-3 which are characterized by emission line rest-frame equivalent widths (REWs) of >~150Ang , much larger than those of typical quasars. In particular, 55% (45%) of the non-BAL Type 1s with measurable CIV in our sample have REW(CIV) > 100 (150)Ang, compared to only 5.8% (1.3%) for non-BAL quasars in BOSS. These objects often also have unusual line ratios, such as very high NV/Ly-alpha ratios. These large REWs might be caused by suppressed continuum emission analogous to Type 2 quasars; however, there is no obvious mechanism in standard Unified Models to suppress the continuum without also obscuring the broad emission lines.

preprint2016arXiv

Insights into Quasar UV Spectra Using Unsupervised Clustering Analysis

Machine learning can provide powerful tools to detect patterns in multi-dimensional parameter space. We use K-means -a simple yet powerful unsupervised clustering algorithm which picks out structure in unlabeled data- to study a sample of quasar UV spectra from the Quasar Catalog of the 10th Data Release of the Sloan Digital Sky Survey of Paris et al. (2014). Detecting patterns in large datasets helps us gain insights into the physical conditions and processes giving rise to the observed properties of quasars. We use K-means to find clusters in the parameter space of the equivalent width (EW), the blue- and red-half-width at half-maximum (HWHM) of the Mg II 2800 A line, the C IV 1549 A line, and the C III] 1908 A blend in samples of Broad Absorption-Line (BAL) and non-BAL quasars at redshift 1.6-2.1. Using this method, we successfully recover correlations well-known in the UV regime such as the anti-correlation between the EW and blueshift of the C IV emission line and the shape of the ionizing Spectra Energy distribution (SED) probed by the strength of He II and the Si III]/C III] ratio. We find this to be particularly evident when the properties of C III] are used to find the clusters, while those of Mg II proved to be less strongly correlated with the properties of the other lines in the spectra such as the width of C IV or the Si III]/C III] ratio. We conclude that unsupervised clustering methods (such as K-means) are powerful methods for finding "natural" binning boundaries in multidimensional datasets and discuss caveats and future work.

preprint2016arXiv

SpIES: The Spitzer IRAC Equatorial Survey

We describe the first data release from the Spitzer-IRAC Equatorial Survey (SpIES); a large-area survey of 115 deg^2 in the Equatorial SDSS Stripe 82 field using Spitzer during its 'warm' mission phase. SpIES was designed to probe sufficient volume to perform measurements of quasar clustering and the luminosity function at z > 3 to test various models for "feedback" from active galactic nuclei (AGN). Additionally, the wide range of available multi-wavelength, multi-epoch ancillary data enables SpIES to identify both high-redshift (z > 5) quasars as well as obscured quasars missed by optical surveys. SpIES achieves 5σ depths of 6.13 μJy (21.93 AB magnitude) and 5.75 μJy (22.0 AB magnitude) at 3.6 and 4.5 microns, respectively - depths significantly fainter than WISE. We show that the SpIES survey recovers a much larger fraction of spectroscopically-confirmed quasars (98%) in Stripe 82 than are recovered by WISE (55%). This depth is especially powerful at high-redshift (z > 3.5), where SpIES recovers 94% of confirmed quasars, whereas WISE only recovers 25%. Here we define the SpIES survey parameters and describe the image processing, source extraction, and catalog production methods used to analyze the SpIES data. In addition to this survey paper, we release 234 images created by the SpIES team and three detection catalogs: a 3.6 μm-only detection catalog containing 6.1 million sources, a 4.5 μm-only detection catalog containing 6.5 million sources, and a dual-band detection catalog containing 5.4 million sources.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: An Investigation of Biases in CIV Emission-Line Properties

We investigate the dependence on data quality of quasar properties measured from the CIV emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey (SDSS) Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch and coadded spectra, focusing on the CIV emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ~2-4), we find increasing systematic offsets with decreasing S/N. The systematic difference (measurement uncertainty) in our lowest S/N (<5) subsample between the single-epoch and coadded spectrum (i) CIV equivalent width is 17A (31A), (ii) centroid wavelength is <1A (2A), and fractional velocity widths, Delta V/V, characterized by (iii) the line dispersion, sigma_l, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1-sigma measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the CIV full-width at half maximum (FWHM) velocity width and the CIV profile characterized by FWHM/sigma_l are only smaller than the statistical uncertainties when S/N>10, although offsets in lower S/N spectra exceed the statistical uncertainties by only a factor of ~1.5. Characterizing the CIV line profile by the kurtosis is the least robust property investigated, as the median systematic coadded--single-epoch measurement differences are larger than the statistical uncertainties for all S/N subsamples.

preprint2016arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Velocity Shifts of Quasar Emission Lines

Quasar emission lines are often shifted from the systemic velocity due to various dynamical and radiative processes in the line-emitting region. The level of these velocity shifts depends both on the line species and on quasar properties. We study velocity shifts for the line peaks of various narrow and broad quasar emission lines relative to systemic using a sample of 849 quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. The coadded (from 32 epochs) spectra of individual quasars have sufficient signal-to-noise ratio (SNR) to measure stellar absorption lines to provide reliable systemic velocity estimates, as well as weak narrow emission lines. The sample also covers a large dynamic range in quasar luminosity (~2 dex), allowing us to explore potential luminosity dependence of the velocity shifts. We derive average line peak velocity shifts as a function of quasar luminosity for different lines, and quantify their intrinsic scatter. We further quantify how well the peak velocity can be measured for various lines as a function of continuum SNR, and demonstrate there is no systematic bias in the line peak measurements when the spectral quality is degraded to as low as SNR~3 per SDSS pixel. Based on the observed line shifts, we provide empirical guidelines on redshift estimation from [OII]3728, [OIII]5008, [NeV]3426, MgII, CIII], HeII1640, broad Hbeta, CIV, and SiIV, which are calibrated to provide unbiased systemic redshifts in the mean, but with increasing intrinsic uncertainties of 46, 56, 119, 205, 233, 242, 400, 415, and 477 km/s, in addition to the measurement uncertainties. These more realistic redshift uncertainties are generally much larger than the formal uncertainties reported by the redshift pipelines for spectroscopic quasar surveys, and demonstrate the infeasibility of measuring quasar redshifts to better than ~200 km/s with only broad lines.

preprint2015arXiv

Are the Variability Properties of the Kepler AGN Light Curves Consistent with a Damped Random Walk?

We test the consistency of active galactic nuclei (AGN) optical flux variability with the $\textit{damped random walk}$ (DRW) model. Our sample consists of 20 multi-quarter $\textit{Kepler}$ AGN light curves including both Type 1 and 2 Seyferts, radio-loud and -quiet AGN, quasars, and blazars. $\textit{Kepler}$ observations of AGN light curves offer a unique insight into the variability properties of AGN light curves because of the very rapid ($11.6-28.6$ min) and highly uniform rest-frame sampling combined with a photometric precision of $1$ part in $10^{5}$ over a period of 3.5 yr. We categorize the light curves of all 20 objects based on visual similarities and find that the light curves fall into 5 broad categories. We measure the first order structure function of these light curves and model the observed light curve with a general broken power-law PSD characterized by a short-timescale power-law index $γ$ and turnover timescale $τ$. We find that less than half the objects are consistent with a DRW and observe variability on short timescales ($\sim 2$ h). The turnover timescale $τ$ ranges from $\sim 10-135$ d. Interesting structure function features include pronounced dips on rest-frame timescales ranging from $10-100$ d and varying slopes on different timescales. The range of observed short-timescale PSD slopes and the presence of dip and varying slope features suggests that the DRW model may not be appropriate for all AGN. We conclude that AGN variability is a complex phenomenon that requires a more sophisticated statistical treatment.

preprint2015arXiv

Bayesian High-Redshift Quasar Classification from Optical and Mid-IR Photometry

We identify 885,503 type 1 quasar candidates to i<22 using the combination of optical and mid-IR photometry. Optical photometry is taken from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III/BOSS), while mid-IR photometry comes from a combination of data from the Wide-Field Infrared Survey Explorer (WISE) "ALLWISE" data release and several large-area Spitzer Space Telescope fields. Selection is based on a Bayesian kernel density algorithm with a training sample of 157,701 spectroscopically-confirmed type-1 quasars with both optical and mid-IR data. Of the quasar candidates, 733,713 lack spectroscopic confirmation (and 305,623 are objects that we have not previously classified as photometric quasar candidates). These candidates include 7874 objects targeted as high probability potential quasars with 3.5<z<5 (of which 6779 are new photometric candidates). Our algorithm is more complete to z>3.5 than the traditional mid-IR selection "wedges" and to 2.2<z<3.5 quasars than the SDSS-III/BOSS project. Number counts and luminosity function analysis suggests that the resulting catalog is relatively complete to known quasars and is identifying new high-z quasars at z>3. This catalog paves the way for luminosity-dependent clustering investigations of large numbers of faint, high-redshift quasars and for further machine learning quasar selection using Spitzer and WISE data combined with other large-area optical imaging surveys.

preprint2015arXiv

Detection of Rest-frame Optical Lines from X-shooter Spectroscopy of Weak Emission Line Quasars

Over the past 15 years, examples of exotic radio-quiet quasars with intrinsically weak or absent broad emission line regions (BELRs) have emerged from large-scale spectroscopic sky surveys. Here, we present spectroscopy of seven such weak emission line quasars (WLQs) at moderate redshifts (z=1.4-1.7) using the X-shooter spectrograph, which provides simultaneous optical and near-infrared spectroscopy covering the rest-frame ultraviolet through optical. These new observations effectively double the number of WLQs with spectroscopy in the optical rest-frame, and they allow us to compare the strengths of (weak) high-ionization emission lines (e.g., CIV) to low-ionization lines (e.g., MgII, Hb, Ha) in individual objects. We detect broad Hb and Ha emission in all objects, and these lines are generally toward the weaker end of the distribution expected for typical quasars (e.g., Hb has rest-frame equivalent widths ranging from 15-40 Ang.). However, these low-ionization lines are not exceptionally weak, as is the case for high-ionization lines in WLQs. The X-shooter spectra also display relatively strong optical FeII emission, Hb FWHM <4000 km/s, and significant CIV blueshifts (1000-5500 km/s) relative to the systemic redshift; two spectra also show elevated ultraviolet FeII emission, and an outflowing component to their (weak) MgII emission lines. These properties suggest that WLQs are exotic versions of "wind-dominated" quasars. Their BELRs either have unusual high-ionization components, or their BELRs are in an atypical photoionization state because of an unusually soft continuum.

preprint2015arXiv

Do the Kepler AGN Light Curves Need Re-processing?

We gauge the impact of spacecraft-induced effects on the inferred variability properties of the light curve of the Seyfert 1 AGN Zw 229-15 observed by \Kepler. We compare the light curve of Zw 229-15 obtained from the Kepler MAST database with a re-processed light curve constructed from raw pixel data (Williams & Carini, 2015). We use the first-order structure function, $SF(δt)$, to fit both light curves to the damped power-law PSD of Kasliwal, Vogeley & Richards, 2015. On short timescales, we find a steeper log-PSD slope ($γ= 2.90$ to within $10$ percent) for the re-processed light curve as compared to the light curve found on MAST ($γ= 2.65$ to within $10$ percent)---both inconsistent with a damped random walk which requires $γ= 2$. The log-PSD slope inferred for the re-processed light curve is consistent with previous results (Carini & Ryle, 2012, Williams & Carini, 2015) that study the same re-processed light curve. The turnover timescale is almost identical for both light curves ($27.1$ and $27.5$~d for the reprocessed and MAST database light curves). Based on the obvious visual difference between the two versions of the light curve and on the PSD model fits, we conclude that there remain significant levels of spacecraft-induced effects in the standard pipeline reduction of the Kepler data. Re-processing the light curves will change the model inferenced from the data but is unlikely to change the overall scientific conclusion reached by Kasliwal et al. 2015---not all AGN light curves are consistent with the DRW.

preprint2015arXiv

Mean and Extreme Radio Properties of Quasars and the Origin of Radio Emission

We investigate the evolution of both the radio-loud fraction (RLF) and (using stacking analysis) the mean radio-loudness of quasars. We consider how these values evolve as a function of redshift and luminosity, black hole (BH) mass and accretion rate, and parameters related to the dominance of a wind in the broad emission line region. We match the FIRST source catalog to samples of luminous quasars (both spectroscopic and photometric), primarily from the Sloan Digital Sky Survey. After accounting for catastrophic errors in BH mass estimates at high-redshift, we find that both the RLF and the mean radio luminosity increase for increasing BH mass and decreasing accretion rate. Similarly both the RLF and mean radio loudness increase for quasars which are argued to have weaker radiation line driven wind components of the broad emission line region. In agreement with past work, we find that the RLF increases with increasing luminosity and decreasing redshift while the mean radio-loudness evolves in the exact opposite manner. This difference in behavior between the mean radio-loudness and the RLF in L-z may indicate selection effects that bias our understanding of the evolution of the RLF; deeper surveys in the optical and radio are needed to resolve this discrepancy. Finally, we argue that radio-loud (RL) and radio-quiet (RQ) quasars may be parallel sequences but where only RQ quasars at one extreme of the distribution are likely to become RL, possibly through slight differences in spin and/or merger history.

preprint2015arXiv

Mining for Dust in Type 1 Quasars

We explore the extinction/reddening of ~35,000 uniformly selected quasars with 0<z<5.3 in order to better understand their intrinsic optical/ultraviolet spectral energy distributions. Using rest-frame optical-UV photometry taken from the Sloan Digital Sky Survey's (SDSS) 7th data release, cross-matched to WISE in the mid-infrared, 2MASS and UKIDSS in the near-infrared, and GALEX in the UV, we isolate outliers in the color distribution and find them well described by an SMC-like reddening law. A hierarchical Bayesian model with a Markov Chain Monte Carlo sampling method was used to find distributions of powerlaw indices and E(B-V) consistent with both the broad absorption line (BAL) and non-BAL samples. We find that, of the ugriz color-selected type 1 quasars in SDSS, 2.5% (13%) of the non-BAL (BAL) sample are consistent with E(B-V)>0.1 and 0.1% (1.3%) with E(B-V)>0.2. Simulations show both populations of quasars are intrinsically bluer than the mean composite, with a mean spectral index ($α_λ$) of -1.79 (-1.83). The emission and absorption-line properties of both samples reveal that quasars with intrinsically red continua have narrower Balmer lines and stronger ionizing spectral lines, the latter indicating a harder continuum in the extreme-UV and the former pointing to differences in black hole mass and/or orientation.

preprint2015arXiv

Quasar Classification Using Color and Variability

We conduct a pilot investigation to determine the optimal combination of color and variability information to identify quasars in current and future multi-epoch optical surveys. We use a Bayesian quasar selection algorithm (Richards et al. 2004) to identify 35,820 type 1 quasar candidates in a 239 square degree field of the Sloan Digital Sky Survey (SDSS) Stripe 82, using a combination of optical photometry and variability. Color analysis is performed on 5-band single- and multi-epoch SDSS optical photometry to a depth of r ~22.4. From these data, variability parameters are calculated by fitting the structure function of each object in each band with a power law model using 10 to >100 observations over timescales from ~1 day to ~8 years. Selection was based on a training sample of 13,221 spectroscopically-confirmed type-1 quasars, largely from the SDSS. Using variability alone, colors alone, and combining variability and colors we achieve 91%, 93%, and 97% quasar completeness and 98%, 98%, and 97% efficiency respectively, with particular improvement in the selection of quasars at 2.7<z<3.5 where quasars and stars have similar optical colors. The 22,867 quasar candidates that are not spectroscopically confirmed reach a depth of i ~22.0; 21,876 (95.7%) are dimmer than coadded i-band magnitude of 19.9, the cut off for spectroscopic follow-up for SDSS on Stripe 82. Brighter than 19.9, we find 5.7% more quasar candidates without confirming spectra in sky regions otherwise considered complete. The resulting quasar sample has sufficient purity (and statistically correctable incompleteness) to produce a luminosity function comparable to those determined by spectroscopic investigations. We discuss improvements that can be made to the process in preparation for performing similar photometric selection and science on data from post-SDSS sky surveys.

preprint2015arXiv

The Discovery of the First "Changing Look" Quasar: New Insights into the Physics & Phenomenology of AGN

SDSS J015957.64+003310.5 is an X-ray selected, $z=0.31$ AGN from the Stripe 82X survey that transitioned from a Type 1 quasar to a Type 1.9 AGN between 2000 and 2010. This is the most distant AGN, and first quasar, yet observed to have undergone such a dramatic change. We re-observed the source with the double spectrograph on the Palomar 5m telescope in July 2014 and found that the spectrum is unchanged since 2010. From fitting the optical spectra, we find that the AGN flux dropped by a factor of 6 between 2000 and 2010 while the broad H$α$ emission faded and broadened. Serendipitous X-ray observations caught the source in both the bright and dim state, showing a similar 2-10 keV flux diminution as the optical while lacking signatures of obscuration. The optical and X-ray changes coincide with $g$-band magnitude variations over multiple epochs of Stripe 82 observations. We demonstrate that variable absorption, as might be expected from the simplest AGN unification paradigm, does not explain the observed photometric or spectral properties. We interpret the changing state of J0159+0033 to be caused by dimming of the AGN continuum, reducing the supply of ionizing photons available to excite gas in the immediate vicinity around the black hole. J0159+0033 provides insight into the intermittency of black hole growth in quasars, as well as an unprecedented opportunity to study quasar physics (in the bright state) and the host galaxy (in the dim state), which has been impossible to do in a single sources until now.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Ensemble Spectroscopic Variability of Quasar Broad Emission Lines

We explore the variability of quasars in the MgII and Hbeta broad emission lines and UV/optical continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic epochs from SDSS-RM over $6$ months, containing 357 quasars with MgII and 41 quasars with Hbeta . On longer timescales, the study is also supplemented with two-epoch data from SDSS-I/II. The SDSS-I/II data include an additional $2854$ quasars with MgII and 572 quasars with Hbeta. The MgII emission line is significantly variable ($Δf/f$ 10% on 100-day timescales), a necessary prerequisite for its use for reverberation mapping studies. The data also confirm that continuum variability increases with timescale and decreases with luminosity, and the continuum light curves are consistent with a damped random-walk model on rest-frame timescales of $\gtrsim 5$ days. We compare the emission-line and continuum variability to investigate the structure of the broad-line region. Broad-line variability shows a shallower increase with timescale compared to the continuum emission, demonstrating that the broad-line transfer function is not a $δ$-function. Hbeta is more variable than MgII (roughly by a factor of $1.5$), suggesting different excitation mechanisms, optical depths and/or geometrical configuration for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM project have important consequences for future studies of reverberation mapping and black hole mass estimation of $1<z<2$ quasars.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: First Broad-line Hbeta and MgII Lags at z>~0.3 from six-Month Spectroscopy

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

preprint2014arXiv

Close companions to two high-redshift quasars

We report the serendipitous discoveries of companion galaxies to two high-redshift quasars. SDSS J025617.7+001904 is a z=4.79 quasar included in our recent survey of faint quasars in the SDSS Stripe 82 region. The initial MMT slit spectroscopy shows excess Lyman alpha emission extending well beyond the quasar's light profile. Further imaging and spectroscopy with LBT/MODS1 confirms the presence of a bright galaxy (i_AB = 23.6) located 2arcsec (12 kpc projected) from the quasar with strong Lyman alpha emission (EW_0 ~ 100Ang) at the redshift of the quasar, as well as faint continuum. The second quasar, CFHQS J005006.6+344522 (z=6.25), is included in our recent HST SNAP survey of z~6 quasars searching for evidence of gravitational lensing. Deep imaging with ACS and WFC3 confirms an optical dropout ~4.5 mag fainter than the quasar (Y_AB=25) at a separation of 0.9 arcsec. The red i_775-Y_105 color of the galaxy and its proximity to the quasar (5 kpc projected if at the quasar redshift) strongly favor an association with the quasar. Although it is much fainter than the quasar it is remarkably bright when compared to field galaxies at this redshift, while showing no evidence for lensing. Both systems may represent late-stage mergers of two massive galaxies, with the observed light for one dominated by powerful ongoing star formation and for the other by rapid black hole growth. Observations of close companions are rare; if major mergers are primarily responsible for high-redshift quasar fueling then the phase when progenitor galaxies can be observed as bright companions is relatively short.

preprint2014arXiv

High Signal-to-Noise Ratio Mid-Infrared Quasar Spectral Templates

Mid-infrared (MIR) quasar spectra exhibit a suite of emission features including high ionization coronal lines from the narrow line region (NLR) illuminated by the ionizing continuum, and hot dust features from grains, as well as polycyclic aromatic hydrocarbons (PAH) features from star formation in the host galaxy. Few features are detected in most spectra because of typically low signal-to-noise ratio (S/N) data. By generating spectral composites in three different luminosity bins from over 180 Spitzer Infrared Spectrograph (IRS) observations, we boost the S/N and reveal important features in the complex spectra. We detect high-ionization, forbidden emission lines in all templates, PAH features in all but the most luminous objects, and broad silicate and graphite features in emission whose strength increases relative to the continuum with luminosity. We find that the intrinsic quasar spectrum for all luminosity templates is consistent, and the differences in the spectra can be explained by host galaxy contamination in the lower luminosity templates. We also posit that star formation may be active in most quasar host galaxies, but the spectral features of star formation are only detectable if the quasar is faint.

preprint2014arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science.

preprint2013arXiv

Characterizing Quasars in the Mid-infrared: High Signal-to-Noise Spectral Templates

Mid-infrared (MIR) quasar spectra exhibit a suite of emission features including high ionization coronal lines from the narrow line region illuminated by the ionizing continuum, broad dust bumps from silicates and graphites, and polycyclic aromatic hydrocarbon (PAH) features from star formation in the host galaxy. However, in Spitzer Infared Spectrograph (IRS) data, few features are detected in most individual spectra because of typically low signal-to-noise ratios (S/N). By generating spectral composites from over 180 IRS observations of Sloan Digital Sky Survey broad-line quasars, we boost the S/N and reveal features in the complex spectra that are otherwise lost in the noise. In addition to an overall composite, we generate composites in three different luminosity bins that span the range of 5.6 micron luminosities of 10$^{40}$--10$^{46}$ (erg~s$^{-1}$). We detect the high-ionization, forbidden emission lines of [SIV], [OIV], and [NeV] 14 micron in all templates and PAH features in all but the most luminous template. Ratios of lines with a range of ionization potentials show no evidence for a strong difference in the shape of the 41--97 eV ionizing continuum over this range of luminosities. The scaling of the emission-line luminosities as a function of continuum luminosity is consistent with what is expected, and shows no indication of a ``disappearing narrow-line region.'' The broad 10 and 18 micron silicate features in emission increase in strength with increasing luminosity, and a broad 3--5 micron black body consistent with graphite emission at 750 K is evident in the highest luminosity template. We find that the intrinsic quasar continua for all luminosity templates are consistent; apparent differences arise primarily from host galaxy contamination most evident at low luminosity.

preprint2013arXiv

Mean Spectral Energy Distributions and Bolometric Corrections for Luminous Quasars

We explore the mid-infrared (mid-IR) through ultraviolet (UV) spectral energy distributions (SEDs) of 119,652 luminous broad-lined quasars with 0.064<z<5.46 using mid-IR data from Spitzer and WISE, near-infrared data from Two Micron All Sky Survey and UKIDSS, optical data from Sloan Digital Sky Survey, and UV data from Galaxy Evolution Explorer. The mean SED requires a bolometric correction (relative to 2500A) of BC=2.75+-0.40 using the integrated light from 1um-2keV, and we further explore the range of bolometric corrections exhibited by individual objects. In addition, we investigate the dependence of the mean SED on various parameters, particularly the UV luminosity for quasars with 0.5<z<3 and the properties of the UV emission lines for quasars with z>1.6; the latter is a possible indicator of the strength of the accretion disk wind, which is expected to be SED dependent. Luminosity-dependent mean SEDs show that, relative to the high-luminosity SED, low-luminosity SEDs exhibit a harder (bluer) far-UV spectral slope, a redder optical continuum, and less hot dust. Mean SEDs constructed instead as a function of UV emission line properties reveal changes that are consistent with known Principal Component Analysis (PCA) trends. A potentially important contribution to the bolometric correction is the unseen extream-UV (EUV) continuum. Our work suggests that lower-luminosity quasars and/or quasars with disk-dominated broad emission lines may require an extra continuum component in the EUV that is not present (or much weaker) in high-luminosity quasars with strong accretion disk winds. As such, we consider four possible models and explore the resulting bolometric corrections. Understanding these various SED-dependent effects will be important for accurate determination of quasar accretion rates.

preprint2013arXiv

Millijansky Radio Variability in SDSS Stripe 82

We report on a blind survey for extragalactic radio variability that was carried out by comparing two epochs of data from the FIRST survey with a third epoch from a new 1.4 GHz survey of SDSS Stripe 82. The three epochs are spaced seven years apart and have an overlapping area of 60 deg^2. We uncover 89 variable sources down to the millijansky level, 75 of which are newly-identified, and we find no evidence for transient phenomena. This new sample of variable sources allows us to infer an upper limit to the mean characteristic timescale of AGN radio variability of 14 years. We find that only 1% of extragalactic sources have fractional variability f_var >3, while 44% of Galactic sources vary by this much. The variable sample contains a larger fraction of quasars than a comparable non-variable control sample, though the majority of the variable sources appear to be extended galaxies in the optical. This implies that either quasars are not the dominant contributor to the variability of the sample, or that the deep optical data allow us to detect the host galaxies of some low-z quasars. We use the new, higher resolution data to report on the morphology of the variable sources. Finally, we show that the fraction of sources that are variable remains constant or increases at low flux densities. This may imply that next generation radio surveys with telescopes like the Australian Square Kilometer Array Pathfinder and MeerKAT will see a constant or even increasing fraction of variable sources down into the submillijansky regime.

preprint2013arXiv

Optical selection of quasars: SDSS and LSST

Over the last decade, quasar sample sizes have increased from several thousand to several hundred thousand, thanks mostly to SDSS imaging and spectroscopic surveys. LSST, the next-generation optical imaging survey, will provide hundreds of detections per object for a sample of more than ten million quasars with redshifts of up to about seven. We briefly review optical quasar selection techniques, with emphasis on methods based on colors, variability properties and astrometric behavior.

preprint2013arXiv

The z=5 Quasar Luminosity Function from SDSS Stripe 82

We present a measurement of the Type I quasar luminosity function at z=5 using a large sample of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright end (M_1450<-26) with Sloan Digital Sky Survey (SDSS) data covering ~6000 deg^2, then extend to lower luminosities (M_1450<-24) with newly discovered, faint z~5 quasars selected from 235 deg^2 of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT and Magellan telescopes. We construct a well-defined sample of 4.7<z<5.1 quasars that is highly complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations reach below the break in the luminosity function (M_1450* ~ -27). The bright end slope is steep (beta <~ -4), with a constraint of beta < -3.1 at 95% confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation for the flattening of the bright end slope reported previously. We find a factor of ~2 greater decrease in the number density of luminous quasars (M_1450<-26) from z=5 to z=6 than from z=4 to z=5, suggesting a more rapid decline in quasar activity at high redshift than found in previous surveys. Our model for the quasar luminosity function predicts that quasars generate ~30% of the ionizing photons required to keep the universe ionized at z=5.

preprint2012arXiv

A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging

We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST.

preprint2012arXiv

CIV Emission as a Probe of Accretion Disk Winds

We present a brief description of a model for the broad emission line region (BELR) in quasars, which is supported by analysis of CIV and other emission lines in the spectra of high-z SDSS quasars. Specifically we consider a two-component BELR with a disk and wind where the relative strength of each component is a function of luminosity. The implications of such a model for our understanding of quasar outflows and estimates of their black hole masses and accretion rates are discussed.

preprint2012arXiv

PHL 1811 Analogs: A Population of X-ray Weak Quasars

We report on a population of X-ray weak quasars with similar UV emission-line properties to those of the remarkable quasar PHL 1811. All radio-quiet PHL 1811 analogs are notably X-ray weak by a mean factor of ~13, with hints of heavy X-ray absorption. Correlations between the X-ray weakness and UV emission-line properties suggest that PHL 1811 analogs may have extreme wind-dominated broad emission-line regions (BELRs). We propose an AGN geometry that can potentially unify the PHL 1811 analogs and the general population of weak-line quasars.

preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

preprint2012arXiv

The Chandra view of the Largest Quasar Lens SDSS J1029+2623

We present results from Chandra observations of the cluster lens SDSS J1029+2623 at z_l=0.58, which is a gravitationally lensed quasar with the largest known image separation. We clearly detect X-ray emission both from the lensing cluster and the three lensed quasar images. The cluster has an X-ray temperature of kT = 8.1 (+2.0, -1.2) keV and bolometric luminosity of L_X = 9.6e44 erg s^-1. Its surface brightness is centered near one of the brightest cluster galaxies, and it is elongated East-West. We identify a subpeak North-West of the main peak, which is suggestive of an ongoing merger. Even so, the X-ray mass inferred from the hydrostatic equilibrium assumption appears to be consistent with the lensing mass from the Einstein radius of the system. We find significant absorption in the soft X-ray spectrum of the faintest quasar image, which can be caused by an intervening material at either the lens or source redshift. The X-ray flux ratios between the quasar images (after correcting for absorption) are in reasonable agreement with those at optical and radio wavelengths, and all the flux ratios are inconsistent with those predicted by simple mass models. This implies that microlensing effect is not significant for this system and dark matter substructure is mainly responsible for the anomalous flux ratios.

preprint2012arXiv

The Hidden Fortress: Structure and substructure of the complex strong lensing cluster SDSS J1029+2623

We present Hubble Space Telescope (HST) ACS and WFC3 observations of SDSS J1029+2623, a three-image quasar lens system produced by a foreground cluster at z=0.584. Our strong lensing analysis reveals 6 additional multiply imaged galaxies. We confirm the complex nature of the mass distribution of the lensing cluster, with a bimodal distribution which deviates from the Chandra X-ray SB distribution. The Einstein radius is estimated to be θ_E=15.2" \pm 0.5" for the quasar redshift of z=2.197. We derive a radial mass distribution from the combination of strong lensing, HST/ACS weak lensing, and Subaru/Suprime-cam weak lensing analysis results, finding a best-fit virial mass of M_vir=(1.5+0.40-0.35) \times 10^14 h^-1 M_sun and a concentration parameter of c_vir=25.7+14.1-7.5. The lensing mass estimate at the outer radius is smaller than the X-ray mass estimate by a factor of ~2. We ascribe this large mass discrepancy to shock heating of the intracluster gas during a merger, which is also suggested by the complex mass and gas distributions and the high value of the concentration parameter. In the HST image, we also identify a probable galaxy, GX, in the vicinity of the faintest quasar image C. The inclusion of GX explains the anomalous flux ratios between the quasar images. The morphology of the highly elongated quasar host galaxy is also well reproduced. The best-fit model suggests large total magnifications of 30 for the quasar and 35 for the quasar host galaxy, and has an AB time delay consistent with the measured value. (Abridged)

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: The Quasar Luminosity Function from Data Release Nine

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $Λ$CDM cosmology with H0 = 70 km/s/Mpc, $Ω_{\rm M}$ = 0.3, and $Ω_Λ$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

preprint2012arXiv

The Sloan Digital Sky Survey Quasar Lens Search. V. Final Catalog from the Seventh Data Release

We present the final statistical sample of lensed quasars from the Sloan Digital Sky Survey (SDSS) Quasar Lens Search (SQLS). The well-defined statistical lens sample consists of 26 lensed quasars brighter than i=19.1 and in the redshift range of 0.6<z<2.2 selected from 50,836 spectroscopically confirmed quasars in the SDSS Data Release 7 (DR7), where we restrict the image separation range to 1"<θ<20" and the i-band magnitude differences in two image lenses to be smaller than 1.25 mag. The SDSS DR7 quasar catalog also contains 36 additional lenses identified with various techniques. In addition to these lensed quasars, we have identified 81 pairs of quasars from follow-up spectroscopy, 26 of which are physically associated binary quasars. The statistical lens sample covers a wide range of image separations, redshifts, and magnitudes, and therefore is suitable for systematic studies of cosmological parameters and surveys of the structure and evolution of galaxies and quasars.

preprint2012arXiv

The Sloan Digital Sky Survey Quasar Lens Search. VI. Constraints on Dark Energy and the Evolution of Massive Galaxies

We present a statistical analysis of the final lens sample from the Sloan Digital Sky Survey Quasar Lens Search (SQLS). The number distribution of a complete subsample of 19 lensed quasars selected from 50,836 source quasars is compared with theoretical expectations, with particular attention to the selection function. Assuming that the velocity function of galaxies does not evolve with redshift, the SQLS sample constrains the cosmological constant to Ω_Λ=0.79^{+0.06}_{-0.07}(stat.)^{+0.06}_{-0.06}(syst.) for a flat universe. The dark energy equation of state is found to be consistent with w=-1 when the SQLS is combined with constraints from baryon acoustic oscillation (BAO) measurements or results from the Wilkinson Microwave Anisotropy Probe (WMAP). We also obtain simultaneous constraints on cosmological parameters and redshift evolution of the galaxy velocity function, finding no evidence for redshift evolution at z<1 in any combinations of constraints. For instance, number density evolution quantified as ν_n=d\lnϕ_*/d\ln(1+z) and the velocity dispersion evolution ν_σ=d\lnσ_*/d\ln(1+z) are constrained to ν_n=1.06^{+1.36}_{-1.39}(stat.)^{+0.33}_{-0.64}(syst.) and ν_σ=-0.05^{+0.19}_{-0.16}(stat.)^{+0.03}_{-0.03}(syst.) respectively when the SQLS result is combined with BAO and WMAP for flat models with a cosmological constant. We find that a significant amount of dark energy is preferred even after fully marginalizing over the galaxy evolution parameters. Thus the statistics of lensed quasars robustly confirm the accelerated cosmic expansion.

preprint2011arXiv

A Catalog of Quasar Properties from SDSS DR7

We present a compilation of properties of the 105,783 quasars in the SDSS Data Release 7 (DR7) quasar catalog. In this value-added product, we compile continuum and emission line measurements around the Halpha, Hbeta, MgII and CIV regions, as well as other quantities such as radio properties, broad absorption line quasar (BALQSO) flags, and disk emitters. We also compile virial black hole mass estimates based on various calibrations. For the fiducial virial mass estimates we use the Vestergaard & Peterson (VP06) calibrations for Hbeta and CIV, and our own calibration for MgII which matches the VP06 Hbeta masses on average. We describe the construction of this catalog, and discuss its limitations. The catalog and its future updates will be made publicly available online.

preprint2011arXiv

A Multiwavelength Study of Binary Quasars and Their Environments

We present Chandra X-ray imaging and spectroscopy for 14 quasars in spatially resolved pairs, part of a complete sample of binary quasars with small transverse separations drawn from Sloan Digital Sky Survey (DR6) photometry. We find no significant difference in X-ray properties when compared with large control samples of isolated quasars. We present infrared photometry from our observations with SWIRC at the MMT, and from the WISE Preliminary Data Release, and fit simple spectral energy distributions to all 14 QSOs. We find preliminary evidence that substantial contributions from star formation are required, but possibly no more so than for isolated X-ray-detected QSOs. Sensitive searches of the X-ray images for extended emission, and the optical images for optical galaxy excess show that these binary QSOs are not preferentially found in rich cluster environments. While larger binary QSO samples with richer far-IR and sub-millimeter multiwavelength data might better reveal signatures of merging and triggering, optical color-selection of QSO pairs may be biased against such signatures. X-ray and/or variability selection of QSO pairs, while challenging, should be attempted. We present in our Appendix a primer on X-ray flux and luminosity calculations.

preprint2011arXiv

A Population of X-ray Weak Quasars: PHL 1811 Analogs at High Redshift

We report the results from Chandra and XMM-Newton observations of a sample of 10 type 1 quasars selected to have unusual UV emission-line properties (weak and blueshifted high-ionization lines; strong UV Fe emission) similar to those of PHL 1811, a confirmed intrinsically X-ray weak quasar. These quasars were identified by the Sloan Digital Sky Survey at high redshift (z~2.2); eight are radio quiet while two are radio intermediate. All of the radio-quiet PHL 1811 analogs are notably X-ray weak by a mean factor of ~13. These sources lack broad absorption lines and have blue UV/optical continua, suggesting they are intrinsically X-ray weak. However, their average X-ray spectrum appears to be harder than those of typical quasars, which may indicate the presence of heavy intrinsic X-ray absorption. Our radio-quiet PHL 1811 analogs support a connection between an X-ray weak spectral energy distribution and PHL 1811-like UV emission lines; this connection provides an economical way to identify X-ray weak type 1 quasars. The fraction of radio-quiet PHL 1811 analogs in the radio-quiet quasar population is estimated to be < 1.2%. We have investigated correlations between relative X-ray brightness and UV emission-line properties for a sample combining radio-quiet PHL 1811 analogs, PHL 1811, and typical type 1 quasars. These correlation analyses suggest that PHL 1811 analogs may have extreme wind-dominated broad emission-line regions. Observationally, radio-quiet PHL 1811 analogs appear to be a subset (~30%) of radio-quiet weak-line quasars. The existence of a subset of quasars in which high-ionization "shielding gas" covers most of the BELR, but little more than the BELR, could potentially unify the PHL 1811 analogs and WLQs. The two radio-intermediate PHL 1811 analogs are X-ray bright. One of them appears to have jet-dominated X-ray emission, while the nature of the other remains unclear.

preprint2011arXiv

Analyzing the Flux Anomalies of the Large-Separation Lensed Quasar SDSS J1029+2623

Using a high resolution radio image, we successfully resolve the two fold image components B and C of the quasar lens system SDSS J1029+2623. The flux anomalies associated with these two components in the optical regime persist, albeit less strongly, in our radio observations, suggesting that the cluster must be modeled by something more than a single central potential. We argue that placing substructure close to one of the components can account for a flux anomaly with negligible changes in the component positions. Our best fit model has a substructure mass of ~10^8 solar masses up to the mass-sheet degeneracy, located roughly 0.1 arcsecs West and 0.1 arcsecs North of component B. We demonstrate that a positional offset between the centers of the source components can explain the differences between the optical and radio flux ratios.

preprint2011arXiv

CIV Emission and the Ultraviolet through X-ray Spectral Energy Distribution of Radio-Quiet Quasars

In the restframe UV, two of the parameters that best characterize the range of emission-line properties in quasar broad emission-line regions are the equivalent width and the blueshift of the CIV line relative to the quasar rest frame. We explore the connection between these emission-line properties and the UV through X-ray spectral energy distribution (SED) for radio-quiet (RQ) quasars. Our sample consists of a heterogeneous compilation of 406 quasars from the Sloan Digital Sky Survey and Palomar-Green survey that have well-measured CIV emission-line and X-ray properties (including 164 objects with measured Gamma). We find that RQ quasars with both strong CIV emission and small CIV blueshifts can be classified as "hard-spectrum" sources that are (relatively) strong in the X-ray as compared to the UV. On the other hand, RQ quasars with both weak CIV emission and large CIV blueshifts are instead "soft-spectrum" sources that are (relatively) weak in the X-ray as compared to the UV. This work helps to further bridge optical/soft X-ray "Eigenvector 1" relationships to the UV and hard X-ray. Based on these findings, we argue that future work should consider systematic errors in bolometric corrections (and thus accretion rates) that are derived from a single mean SED. Detailed analysis of the CIV emission line may allow for SED-dependent corrections to these quantities.

preprint2011arXiv

Discovery of a Radio-Selected z ~ 6 Quasar

We present the discovery of only the second radio-selected, z ~ 6 quasar. We identified SDSS J222843.54+011032.2 (z=5.95) by matching the optical detections of the deep Sloan Digital Sky Survey (SDSS) Stripe 82 with their radio counterparts in the Stripe82 VLA Survey. We also matched the Canadian-France-Hawaiian Telescope Legacy Survey Wide (CFHTLS Wide) with the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey but have yet to find any z ~ 6 quasars in this survey area. The discovered quasar is optically-faint, z = 22.3 and M_{1450} ~ -24.5, but radio-bright, with a flux density of f$_{1.4GHz, peak}$ = 0.31mJy and a radio-loudness of R ~ 1100 (where R = f_{5GHz}/f_{2500}). The i-z color of the discovered quasar places it outside the color selection criteria for existing optical surveys. We conclude by discussing the need for deeper wide-area radio surveys in the context of high-redshift quasars.

preprint2011arXiv

Dusty Tori of Luminous Type 1 Quasars at z \sim 2

We present Spitzer infrared spectra and ultra-violet to mid-infrared spectral energy distributions (SEDs) of 25 luminous type 1 quasars at z \sim 2. In general, the spectra show a bump peaking around 3 μm, and the 10 μm silicate emission feature. The 3 μm emission is identified with hot dust emission at its sublimation temperature. We explore two approaches to modeling the SED: (i) using the Clumpy model SED from Nenkova et al. (2008a), and (ii) the Clumpy model SED, and an additional blackbody component to represent the 3 μm emission. In the first case, a parameter search of \sim 1.25 million Clumpy models shows: (i) if we ignore the UV-to-near-IR SED, models fit the 2-8 μm region well, but not the 10 μm feature; (ii) if we include the UV-to-near-IR SED in the fit, models do not fit the 2-8 μm region. The observed 10 μm features are broader and shallower than those in the best-fit models in the first approach. In the second case, the shape of the 10 μm feature is better reproduced by the Clumpy models. The additional blackbody contribution in the 2-8 μm range allows Clumpy models dominated by cooler temperatures (T < 800K) to better fit the 8-12μm SED. A centrally concentrated distribution of a small number of torus clouds is required in the first case, while in the second case the clouds are more spread out radially. The temperature of the blackbody component is ~ 1200 K as expected for graphite grains.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

preprint2011arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: Quasar Target Selection for Data Release Nine

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg^2, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of a sample of ~150,000 z>2.2 quasars. Along with measuring the angular diameter distance at z\approx2.5, BOSS will provide the first direct measurement of the expansion rate of the Universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars over 2.2 < z < 3.5, where their colors overlap those of stars. During the first year of the BOSS survey, quasar target selection methods were developed and tested to meet the requirement of delivering at least 15 quasars deg^-2 in this redshift range, out of 40 targets deg^-2. To achieve these surface densities, the magnitude limit of the quasar targets was set at g <= 22.0 or r<=21.85. While detection of the BAO signature in the Ly-alpha absorption in quasar spectra does not require a uniform target selection, many other astrophysical studies do. We therefore defined a uniformly-selected subsample of 20 targets deg^-2, for which the selection efficiency is just over 50%. This "CORE" subsample will be fixed for Years Two through Five of the survey. In this paper we describe the evolution and implementation of the BOSS quasar target selection algorithms during the first two years of BOSS operations. We analyze the spectra obtained during the first year. 11,263 new z>2.2 quasars were spectroscopically confirmed by BOSS. Our current algorithms select an average of 15 z > 2.2 quasars deg^-2 from 40 targets deg^-2 using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS quasar target selection. [Abridged]

preprint2011arXiv

The Ultraviolet-to-Mid-Infrared Spectral Energy Distribution of Weak Emission Line Quasars

We present Spitzer Space Telescope photometry of 18 Sloan Digital Sky Survey (SDSS) quasars at 2.7 <= z <= 5.9 which have weak or undetectable high-ionization emission lines in their rest-frame ultraviolet (UV) spectra (hereafter weak-lined quasars, or WLQs). The Spitzer data are combined with SDSS spectra and ground-based, near-infrared (IR) photometry of these sources to produce a large inventory of spectral energy distributions (SEDs) of WLQs across the rest-frame ~0.1-5 mum spectral band. The SEDs of our sources are inconsistent with those of BL Lacertae objects which are dominated by synchrotron emission due to a jet aligned close to our line-of-sight, but are consistent with the SED of ordinary quasars with similar luminosities and redshifts that exhibit a near-to-mid-IR 'bump', characteristic of hot dust emission. This indicates that broad emission lines in WLQs are intrinsically weak, rather than suffering continuum dilution from a jet, and that such sources cannot be selected efficiently from traditional photometric surveys.

preprint2011arXiv

Unification of Luminous Type 1 Quasars through CIV Emission

Using a sample of 30,000 quasars from SDSS-DR7, we explore the range of properties exhibited by high-ionization, broad emission lines, such as CIV 1549. Specifically we investigate the anti-correlation between L_UV and emission line EQW (the Baldwin Effect) and the "blueshifting" of high-ionization emission lines. The blueshift of the CIV emission line is nearly ubiquitous, with a mean shift of 810 km/s for radio-quiet (RQ) quasars and 360 km/s for radio-loud (RL) quasars, and the Baldwin Effect is present in both RQ and RL samples. Composite spectra are constructed as a function of CIV emission line properties in attempt to reveal empirical relationships between different line species and the SED. Within a two-component disk+wind model of the broad emission line region (BELR), where the wind filters the continuum seen by the disk component, we find that RL quasars are consistent with being dominated by the disk component, while BALQSOs are consistent with being dominated by the wind component. Some RQ objects have emission line features similar to RL quasars; they may simply have insufficient black hole (BH) spin to form radio jets. Our results suggest that there could be significant systematic errors in the determination of L_bol and BH mass that make it difficult to place these findings in a more physical context. However, it is possible to classify quasars in a paradigm where the diversity of BELR parameters are due to differences in an accretion disk wind between quasars (and over time); these differences are underlain primarily by the SED, which ultimately must be tied to BH mass and accretion rate.

preprint2010arXiv

A strong redshift dependence of the broad absorption line quasar fraction

We describe the application of non-negative matrix factorisation to generate compact reconstructions of quasar spectra from the Sloan Digital Sky Survey (SDSS), with particular reference to broad absorption line quasars (BALQSOs). BAL properties are measured for SiIV lambda1400, CIV lambda1550, AlIII lambda1860 and MgII lambda2800, resulting in a catalogue of 3547 BALQSOs. Two corrections, based on extensive testing of synthetic BALQSO spectra, are applied in order to estimate the intrinsic fraction of CIV BALQSOs. First, the probability of an observed BALQSO spectrum being identified as such by our algorithm is calculated as a function of redshift, signal-to-noise ratio and BAL properties. Second, the different completenesses of the SDSS target selection algorithm for BALQSOs and non-BAL quasars are quantified. Accounting for these selection effects the intrinsic CIV BALQSO fraction is 41+/-5 per cent. Our analysis of the selection effects allows us to measure the dependence of the intrinsic CIV BALQSO fraction on luminosity and redshift. We find a factor of 3.5+/-0.4 decrease in the intrinsic fraction from the highest redshifts, z~4.0, down to z~2.0. The redshift dependence implies that an orientation effect alone is not sufficient to explain the presence of BAL troughs in some but not all quasar spectra. Our results are consistent with the intrinsic BALQSO fraction having no strong luminosity dependence, although with 3-sigma limits on the rate of change of the intrinsic fraction with luminosity of -6.9 and 7.0 per cent dex^-1 we are unable to rule out such a dependence.

preprint2010arXiv

Dust-free quasars in the early Universe

The most distant quasars known, at redshifts z=6, generally have properties indistinguishable from those of lower-redshift quasars in the rest-frame ultraviolet/optical and X-ray bands. This puzzling result suggests that these distant quasars are evolved objects even though the Universe was only seven per cent of its current age at these redshifts. Recently one z=6 quasar was shown not to have any detectable emission from hot dust, but it was unclear whether that indicated different hot-dust properties at high redshift or if it is simply an outlier. Here we report the discovery of a second quasar without hot-dust emission in a sample of 21 z=6 quasars. Such apparently hot-dust-free quasars have no counterparts at low redshift. Moreover, we demonstrate that the hot-dust abundance in the 21 quasars builds up in tandem with the growth of the central black hole, whereas at low redshift it is almost independent of the black hole mass. Thus z=6 quasars are indeed at an early evolutionary stage, with rapid mass accretion and dust formation. The two hot-dust-free quasars are likely to be first-generation quasars born in dust-free environments and are too young to have formed a detectable amount of hot dust around them.

preprint2010arXiv

Fold Lens Flux Anomalies: A Geometric Approach

We develop a new approach for studying flux anomalies in quadruply-imaged fold lens systems. We show that in the absence of substructure, microlensing, or differential absorption, the expected flux ratios of a fold pair can be tightly constrained using only geometric arguments. We apply this technique to 11 known quadruple lens systems in the radio and infrared, and compare our estimates to the Monte Carlo based results of Keeton, Gaudi, and Petters. We show that a robust estimate for a flux ratio from a smoothly varying potential can be found, and at long wavelengths those lenses deviating from from this ratio almost certainly contain significant substructure.

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Catalog V. Seventh Data Release

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Lens Search. IV. Statistical Lens Sample from the Fifth Data Release

We present the second report of our systematic search for strongly lensed quasars from the data of the Sloan Digital Sky Survey (SDSS). From extensive follow-up observations of 136 candidate objects, we find 36 lenses in the full sample of 77,429 spectroscopically confirmed quasars in the SDSS Data Release 5. We then define a complete sample of 19 lenses, including 11 from our previous search in the SDSS Data Release 3, from the sample of 36,287 quasars with i<19.1 in the redshift range 0.6<z<2.2, where we require the lenses to have image separations of 1"<θ<20" and i-band magnitude differences between the two images smaller than 1.25 mag. Among the 19 lensed quasars, 3 have quadruple-image configurations, while the remaining 16 show double images. This lens sample constrains the cosmological constant to be Ω_Λ=0.84^{+0.06}_{-0.08}(stat.)^{+0.09}_{-0.07}(syst.) assuming a flat universe, which is in good agreement with other cosmological observations. We also report the discoveries of 7 binary quasars with separations ranging from 1.1" to 16.6", which are identified in the course of our lens survey. This study concludes the construction of our statistical lens sample in the full SDSS-I data set.

preprint2010arXiv

Weak-Line Quasars at High Redshift: Extremely High Accretion Rates or Anemic Broad-Line Regions?

We present Gemini-North K-band spectra of two representative members of the class of high-redshift quasars with exceptionally weak rest-frame ultraviolet emission lines (WLQs), SDSS J114153.34+021924.3 at z=3.55 and SDSS J123743.08+630144.9 at z=3.49. In both sources we detect an unusually weak broad H_beta line and we place tight upper limits on the strengths of their [O III] lines. Virial, H_beta-based black-hole mass determinations indicate normalized accretion rates of L/L_Edd=0.4 for these sources, which is well within the range observed for typical quasars with similar luminosities and redshifts. We also present high-quality XMM-Newton imaging spectroscopy of SDSS J114153.34+021924.3 and find a hard-X-ray photon index of Gamma=1.91^{+0.24}_{-0.22} which supports the virial L/L_Edd determination in this source. Our results suggest that the weakness of the broad-emission lines in WLQs is not a consequence of an extreme continuum-emission source but instead due to abnormal broad-emission line region properties.

preprint2009arXiv

A Survey of z~6 Quasars in the SDSS Deep Stripe. II. Discovery of Six Quasars at z_{AB}>21

We present the discovery of six new quasars at z~6 selected from the Sloan Digital Sky Survey (SDSS) southern survey, a deep imaging survey obtained by repeatedly scanning a stripe along the celestial equator. The six quasars are about two magnitudes fainter than the luminous z~6 quasars found in the SDSS main survey and one magnitude fainter than the quasars reported in Paper I (Jiang et al. 2008). Four of them comprise a complete flux-limited sample at 21<z_AB<21.8 over an effective area of 195 deg^2. The other two quasars are fainter than z_AB=22 and are not part of the complete sample. The quasar luminosity function at z~6 is well described as a single power law Φ(L_{1450}) \propto L_{1450}^β over the luminosity range -28<M_{1450}<-25. The best-fitting slope βvaries from -2.6 to -3.1, depending on the quasar samples used, with a statistical error of 0.3-0.4. About 40% of the quasars discovered in the SDSS southern survey have very narrow Lya emission lines, which may indicate small black hole masses and high Eddington luminosity ratios, and therefore short black hole growth time scales for these faint quasars at early epochs.

preprint2009arXiv

Astrometric Redshifts for Quasars

The wavelength dependence of atmospheric refraction causes differential chromatic refraction (DCR), whereby objects imaged at different optical/UV wavelengths are observed at slightly different positions in the plane of the detector. Strong spectral features induce changes in the effective wavelengths of broad-band filters that are capable of producing significant positional offsets with respect to standard DCR corrections. We examine such offsets for broad-emission-line (type 1) quasars from the Sloan Digital Sky Survey (SDSS) spanning 0<z<5 and an airmass range of 1.0 to 1.8. These offsets are in good agreement with those predicted by convolving a composite quasar spectrum with the SDSS bandpasses as a function of redshift and airmass. This astrometric information can be used to break degeneracies in photometric redshifts of quasars (or other emission-line sources) and, for extreme cases, may be suitable for determining "astrometric redshifts". On the SDSS's southern equatorial stripe, where it is possible to average many multi-epoch measurements, more than 60% of quasars have emission-line-induced astrometric offsets larger than the SDSS's relative astrometric errors of 25-35 mas. Folding these astrometric offsets into photometric redshift estimates yields an improvement of 9% within Delta z+/-0.1. Future multi-epoch synoptic surveys such as LSST and Pan-STARRS could benefit from intentionally making ~10 observations at relatively high airmass (AM~1.4) in order to improve their photometric redshifts for quasars.

preprint2009arXiv

Binary Quasars at High Redshift I: 24 New Quasar Pairs at z ~ 3-4

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

preprint2009arXiv

Clustering of Low-Redshift (z <= 2.2) Quasars from the Sloan Digital Sky Survey

We present measurements of the quasar two-point correlation function, ξ_{Q}, over the redshift range z=0.3-2.2 based upon data from the SDSS. Using a homogeneous sample of 30,239 quasars with spectroscopic redshifts from the DR5 Quasar Catalogue, our study represents the largest sample used for this type of investigation to date. With this redshift range and an areal coverage of approx 4,000 deg^2, we sample over 25 h^-3 Gpc^3 (comoving) assuming the current LCDM cosmology. Over this redshift range, we find that the redshift-space correlation function, xi(s), is adequately fit by a single power-law, with s_{0}=5.95+/-0.45 h^-1 Mpc and γ_{s}=1.16+0.11-0.16 when fit over s=1-25 h^-1 Mpc. Using the projected correlation function we calculate the real-space correlation length, r_{0}=5.45+0.35-0.45 h^-1 Mpc and γ=1.90+0.04-0.03, over scales of rp=1-130 h^-1 Mpc. Dividing the sample into redshift slices, we find very little, if any, evidence for the evolution of quasar clustering, with the redshift-space correlation length staying roughly constant at s_{0} ~ 6-7 h^-1 Mpc at z<2.2 (and only increasing at redshifts greater than this). Comparing our clustering measurements to those reported for X-ray selected AGN at z=0.5-1, we find reasonable agreement in some cases but significantly lower correlation lengths in others. We find that the linear bias evolves from b~1.4 at z=0.5 to b~3 at z=2.2, with b(z=1.27)=2.06+/-0.03 for the full sample. We compare our data to analytical models and infer that quasars inhabit dark matter haloes of constant mass M ~2 x 10^12 h^-1 M_Sol from redshifts z~2.5 (the peak of quasar activity) to z~0. [ABRIDGED]

preprint2009arXiv

Eight-Dimensional Mid-Infrared/Optical Bayesian Quasar Selection

We explore the multidimensional, multiwavelength selection of quasars from mid-IR (MIR) plus optical data, specifically from Spitzer-IRAC and the Sloan Digital Sky Survey (SDSS). We apply modern statistical techniques to combined Spitzer MIR and SDSS optical data, allowing up to 8-D color selection of quasars. Using a Bayesian selection method, we catalog 5546 quasar candidates to an 8.0 um depth of 56 uJy over an area of ~24 sq. deg; ~70% of these candidates are not identified by applying the same Bayesian algorithm to 4-color SDSS optical data alone. Our selection recovers 97.7% of known type 1 quasars in this area and greatly improves the effectiveness of identifying 3.5<z<5 quasars. Even using only the two shortest wavelength IRAC bandpasses, it is possible to use our Bayesian techniques to select quasars with 97% completeness and as little as 10% contamination. This sample has a photometric redshift accuracy of 93.6% (Delta Z +/-0.3), remaining roughly constant when the two reddest MIR bands are excluded. While our methods are designed to find type 1 (unobscured) quasars, as many as 1200 of the objects are type 2 (obscured) quasar candidates. Coupling deep optical imaging data with deep mid-IR data could enable selection of quasars in significant numbers past the peak of the quasar luminosity function (QLF) to at least z~4. Such a sample would constrain the shape of the QLF and enable quasar clustering studies over the largest range of redshift and luminosity to date, yielding significant gains in our understanding of quasars and the evolution of galaxies.

preprint2009arXiv

High-Redshift SDSS Quasars with Weak Emission Lines

We identify a sample of 74 high-redshift quasars (z>3) with weak emission lines from the Fifth Data Release of the Sloan Digital Sky Survey and present infrared, optical, and radio observations of a subsample of four objects at z>4. These weak emission-line quasars (WLQs) constitute a prominent tail of the Lya+NV equivalent width distribution, and we compare them to quasars with more typical emission-line properties and to low-redshift active galactic nuclei with weak/absent emission lines, namely BL Lac objects. We find that WLQs exhibit hot (T~1000 K) thermal dust emission and have rest-frame 0.1-5 micron spectral energy distributions that are quite similar to those of normal quasars. The variability, polarization, and radio properties of WLQs are also different from those of BL Lacs, making continuum boosting by a relativistic jet an unlikely physical interpretation. The most probable scenario for WLQs involves broad-line region properties that are physically distinct from those of normal quasars.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2009arXiv

The 2dF-SDSS LRG and QSO Survey: the QSO luminosity function at 0.4<z<2.6

We present the QSO luminosity function of the completed 2dF-SDSS LRG and QSO (2SLAQ) survey, based on QSOs photometrically selected from Sloan Digital Sky Survey imaging data and then observed spectroscopically using the 2dF instrument on the Anglo-Australian Telescope. We analyse 10637 QSOs in the redshift range 0.4<z<2.6 to a g-band flux limit of 21.85 (extinction corrected) and an absolute continuum magnitude of Mg(z=2)<-21.5. This sample covers an area of 191.9 deg^2. The binned QSO luminosity function agrees with that of the brighter SDSS main QSO sample, but extends ~2.5 mags fainter, clearly showing the flattening of the luminosity function towards faint absolute magnitudes. 2SLAQ finds an excess of QSOs compared to the 2dF QSO Redshift Survey at g>20.0, as found previously by Richards et al. (2005). The luminosity function is consistent with other previous, much smaller, samples produced to the depth of 2SLAQ. By combining the 2SLAQ and SDSS QSO samples we produce a QSO luminosity function with an unprecedented combination of precision and dynamic range. With this we are able to accurately constrain both the bright and faint ends of the QSO LF. While the overall trends seen in the evolution of the QSO LF appear similar to pure luminosity evolution, the data show very significant departures from such a model. Most notably we see clear evidence that the number density of faint QSOs peaks at lower redshift than bright QSOs: QSOs with Mg>-23 have space densities which peak at z<1, while QSOs at Mg<-26 peak at z>2. By fitting simple luminosity function models in narrow Mg intervals we find that this downsizing is significant at the 99.98 per cent level (abridged).

preprint2009arXiv

X-ray Insights into the Nature of Weak Emission-Line Quasars at High Redshift

We present Chandra observations of nine high-redshift quasars (z=2.7-5.9) discovered by the Sloan Digital Sky Survey with weak or undetectable high-ionization emission lines in their UV spectra (WLQs). Adding archival X-ray observations of six additional sources of this class has enabled us to place the strongest constraints yet on the X-ray properties of this remarkable class of AGNs. Although our data cannot rule out the possibility that the emission lines are overwhelmed by a relativistically boosted continuum, as manifested by BL Lac objects, we find that WLQs are considerably weaker in the X-ray and radio bands than the majority of BL Lacs found at much lower redshifts. If WLQs are high-redshift BL Lacs, then it is difficult to explain the lack of a large parent population of X-ray and radio bright weak-lined sources at high redshift. We also consider the possibility that WLQs are quasars with extreme properties, and in particular that the emission lines are suppressed by high accretion rates. Using joint spectral fitting of the X-ray spectra of 11 WLQs we find that the mean photon index in the hard X-ray band is consistent with those observed in typical radio-quiet AGNs with no hint of an unusually steep hard-X-ray spectrum. This result poses a challenge to the hypothesis that WLQs have extremely high accretion rates, and we discuss additional observations required to test this idea.

preprint2008arXiv

Efficient Photometric Selection of Quasars from the Sloan Digital Sky Survey: II. ~1,000,000 Quasars from Data Release Six

We present a catalog of 1,172,157 quasar candidates selected from the photometric imaging data of the Sloan Digital Sky Survey (SDSS). The objects are all point sources to a limiting magnitude of i=21.3 from 8417 sq. deg. of imaging from SDSS Data Release 6 (DR6). This sample extends our previous catalog by using the latest SDSS public release data and probing both UV-excess and high-redshift quasars. While the addition of high-redshift candidates reduces the overall efficiency (quasars:quasar candidates) of the catalog to ~80%, it is expected to contain no fewer than 850,000 bona fide quasars -- ~8 times the number of our previous sample, and ~10 times the size of the largest spectroscopic quasar catalog. Cross-matching between our photometric catalog and spectroscopic quasar catalogs from both the SDSS and 2dF Surveys, yields 88,879 spectroscopically confirmed quasars. For judicious selection of the most robust UV-excess sources (~500,000 objects in all), the efficiency is nearly 97% -- more than sufficient for detailed statistical analyses. The catalog's completeness to type 1 (broad-line) quasars is expected to be no worse than 70%, with most missing objects occurring at z<0.7 and 2.5<z<3.0. In addition to classification information, we provide photometric redshift estimates (typically good to Delta z +/- 0.3 [2 sigma]) and cross-matching with radio, X-ray, and proper motion catalogs. Finally, we consider the catalog's utility for determining the optical luminosity function of quasars and are able to confirm the flattening of the bright-end slope of the quasar luminosity function at z~4 as compared to z~2.

preprint2008arXiv

Quasar Clustering from SDSS DR5: Dependences on Physical Properties

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2σ$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

preprint2006arXiv

A Catalog of Broad Absorption Line Quasars from the Sloan Digital Sky Survey Third Data Release

We present a total of 4784 unique broad absorption line quasars from the Sloan Digital Sky Survey Third Data Release. An automated algorithm was used to match a continuum to each quasar and to identify regions of flux at least 10% below the continuum over a velocity range of at least 1000 km/s in the CIV and MgII absorption regions. The model continuum was selected as the best-fit match from a set of template quasar spectra binned in luminosity, emission line width, and redshift, with the power-law spectral index and amount of dust reddening as additional free parameters. We characterize our sample through the traditional ``balnicity'' index and a revised absorption index, as well as through parameters such as the width, outflow velocity, fractional depth and number of troughs. From a sample of 16883 quasars at 1.7<z<4.38, we identify 4386 (26.0%) quasars with broad CIV absorption, of which 1756 (10.4%) satisfy traditional selection criteria. From a sample of 34973 quasars at 0.5<z<2.15, we identify 457 (1.31%) quasars with broad MgII absorption, 191 (0.55%) of which satisfy traditional selection criteria. We also provide a supplementary list of 39 visually identified z>4.38 quasars with broad CIV absorption. We find that broad absorption line quasars may have broader emission lines on average than other quasars.

preprint2006arXiv

QSO Absorption Lines from QSOs

We present the results of a search for metal absorption lines in the spectra of background QSOs whose sightlines pass close to foreground QSOs. We detect MgII(2796,2803) absorption in Sloan Digital Sky Survey (SDSS) spectra of four z>1.5 QSOs whose lines of sight pass within 26-98 kpc of lower redshift (z~0.5-1.5) QSOs. The 100% [4/4 pairs] detection of MgII in the background QSOs is clearly at odds with the incidence of associated (z_abs ~ z_em) systems -- absorbers which exist towards only a few percent of QSOs. Although the quality of our foreground QSO spectra is not as high as the SDSS data, absorption seen towards one of the background QSOs clearly does not show up at the same strength in the spectrum of the corresponding foreground QSO. This implies that the absorbing gas is distributed inhomogeneously around the QSO, presumably as a direct consequence of the anisotropic emission from the central AGN. We discuss possible origins for the MgII lines, including: absorption by gas from the foreground QSO host galaxy; companion galaxies fuelling the QSO through gravitational interactions; and tidal debris left by galaxy mergers or interactions which initiated the QSO activity. No single explanation is entirely satisfactory, and we may well be seeing a mixture of phenomena.

preprint2006arXiv

Quasars Probing Quasars I: Optically Thick Absorbers Near Luminous Quasars

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We search 149 moderate resolution background quasar spectra, from Gemini, Keck, the MMT, and the SDSS to survey Lyman Limit Systems (LLSs) and Damped Ly-alpha systems (DLAs) in the vicinity of 1.8 < z < 4.0 luminous foreground quasars. A sample of 27 new quasar-absorber pairs is uncovered with column densities, 17.2 < log (N_HI/cm^2) < 20.9, and transverse (proper) distances of 22 kpc/h < R < 1.7 Mpc/h, from the foreground quasars. If they emit isotropically, the implied ionizing photon fluxes are a factor of ~ 5-8000 times larger than the ambient extragalactic UV background over this range of distances. The observed probability of intercepting an absorber is very high for small separations: six out of eight projected sightlines with transverse separations R < 150 kpc/h have an absorber coincident with the foreground quasar, of which four have log N_HI > 10^19. The covering factor of log N_HI > 10^19 absorbers is thus ~ 50 % (4/8) on these small scales, whereas < 2% would have been expected at random. There are many cosmological applications of these new sightlines: they provide laboratories for studying fluorescent Ly-alpha recombination radiation from LLSs, constrain the environments, emission geometry, and radiative histories of quasars, and shed light on the physical nature of LLSs and DLAs.

preprint2005arXiv

A Snapshot Survey for Gravitational Lenses Among z>=4.0 Quasars: II. Constraints on the 4.0<z<4.5 Quasar Population

We report on i-band snapshot observations of 157 Sloan Digital Sky Survey (SDSS) quasars at 4<z<5.4 using the Advanced Camera for Surveys on the Hubble Space Telescope (HST) to search for evidence of gravitational lensing of these sources. None of the quasars appear to be strongly lensed and multiply imaged at the angular resolution (~0.1") and sensitivity of HST. The non-detection of strong lensing in these systems constrains the z=4-5 luminosity function to an intrinsic slope of beta>-3.8 (3 sigma), assuming a break in the quasar luminosity function at M_1450^star=-24.5. This constraint is considerably stronger than the limit of beta>-4.63 obtained from the absence of lensing in four z>5.7 quasars. Such constraints are important for our understanding of the true space density of high-redshift quasars and the ionization state of the early universe.

preprint2005arXiv

A Survey of z>5.7 Quasars in the Sloan Digital Sky Survey IV: Discovery of Seven Additional Quasars

We present the discovery of seven quasars at z>5.7, selected from ~2000 deg^2 of multicolor imaging data of the Sloan Digital Sky Survey (SDSS). The new quasars have redshifts z from 5.79 to 6.13. Five are selected as part of a complete flux-limited sample in the SDSS Northern Galactic Cap; two have larger photometric errors and are not part of the complete sample. One of the new quasars, SDSS J1335+3533 (z=5.93), exhibits no emission lines; the 3-sigma limit on the rest-frame equivalent width of Ly alpha+NV line is 5 A. It is the highest redshift lineless quasar known, and could be a gravitational lensed galaxy, a BL Lac object or a new type of quasar. Two new z>6 quasars, SDSS 1250+3130 (z=6.13) and SDSS J1137+3549 (z=6.01), show deep Gunn-Peterson absorption gaps in Ly alpha. These gaps are narrower the complete Gunn-Peterson absorption troughs observed among quasars at z>6.2 and do not have complete Ly beta absorption.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: I. Sample Selection

We have compiled a large sample of low-redshift active galactic nuclei (AGN) identified via their emission line characteristics from the spectroscopic data of the Sloan Digital Sky Survey. Since emission lines are often contaminated by stellar absorption lines, we developed an objective and efficient method of subtracting the stellar continuum from every galaxy spectrum before making emission line measurements. The distribution of the measured H$α$ Full Width at Half Maxima values of emission line galaxies is strongly bimodal, with two populations separated at about 1,200km s$^{-1}$. This feature provides a natural separation between narrow-line and broad-line AGN. The narrow-line AGN are identified using standard emission line ratio diagnostic diagrams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 galaxy spectra selected over 1151 square degrees. This sample is used in a companion paper to determine the emission-line luminosity function of AGN.

preprint2005arXiv

Active Galactic Nuclei in the Sloan Digital Sky Survey: II. Emission-Line Luminosity Function

The emission line luminosity function of active galactic nuclei (AGN) is measured from about 3000 AGN included in the main galaxy sample of the Sloan Digital Sky Survey within a redshift range of $0<z<0.15$. The $\Ha$ and [OIII]$λ5007$ luminosity functions for Seyferts cover luminosity range of $10^{5-9}$$L_\odot$ in H$α$ and the shapes are well fit by broken power laws, without a turnover at fainter nuclear luminosities. Assuming a universal conversion from emission line strength to continuum luminosity, the inferred B band magnitude luminosity function is comparable both to the AGN luminosity function of previous studies and to the low redshift quasar luminosity function derived from the 2dF redshift survey. The inferred AGN number density is approximately 1/5 of all galaxies and about $6\times 10^{-3}$ of the total light of galaxies in the $r$-band comes from the nuclear activity. The numbers of Seyfert 1s and Seyfert 2s are comparable at low luminosity, while at high luminosity, Seyfert 1s outnumber Seyfert 2s by a factor of 2-4. In making the luminosity function measurements, we assumed that the nuclear luminosity is independent of the host galaxy luminosity, an assumption we test {\it a posteriori}, and show to be consistent with the data. Given the relationship between black hole mass and host galaxy bulge luminosity, the lack of correlation between nuclear and host luminosity suggests that the main variable that determines the AGN luminosity is the Eddington ratio, not the black hole mass. This appears to be different from luminous quasars, which are most likely to be shining near the Eddington limit.

preprint2005arXiv

Binary Quasars in the Sloan Digital Sky Survey: Evidence for Excess Clustering on Small Scales

We present a sample of 218 new quasar pairs with proper transverse separations R_prop < 1 Mpc/h over the redshift range 0.5 < z < 3.0, discovered from an extensive follow up campaign to find companions around the Sloan Digital Sky Survey and 2dF Quasar Redshift Survey quasars. This sample includes 26 new binary quasars with separations R_prop < 50 kpc/h (theta < 10 arcseconds), more than doubling the number of such systems known. We define a statistical sample of binaries selected with homogeneous criteria and compute its selection function, taking into account sources of incompleteness. The first measurement of the quasar correlation function on scales 10 kpc/h < R_prop < 400 kpc/h is presented. For R_prop < 40 kpc/h, we detect an order of magnitude excess clustering over the expectation from the large scale R_prop > 3 Mpc/h quasar correlation function, extrapolated down as a power law to the separations probed by our binaries. The excess grows to ~ 30 at R_prop ~ 10 kpc/h, and provides compelling evidence that the quasar autocorrelation function gets progressively steeper on sub-Mpc scales. This small scale excess can likely be attributed to dissipative interaction events which trigger quasar activity in rich environments. Recent small scale measurements of galaxy clustering and quasar-galaxy clustering are reviewed and discussed in relation to our measurement of small scale quasar clustering.

preprint2005arXiv

Discovery of a Fifth Image of the Large Separation Gravitationally Lensed Quasar SDSS J1004+4112

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. A faint point source located 0.2'' from the center of the brightest galaxy in the lensing cluster is detected in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The location and brightness of the point source are consistent with lens model predictions for a lensed image. We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the structure of the lensing critical curves and thereby the mass distribution of the lensing cluster.

preprint2005arXiv

SDSS J0806+2006 and SDSS J1353+1138: Two New Gravitationally Lensed Quasars from the Sloan Digital Sky Survey

We report the discoveries of two, two-image gravitationally lensed quasars selected from the Sloan Digital Sky Survey: SDSS J0806+2006 at z_s=1.540 and SDSS J1353+1138 at z_s=1.629 with image separations of 1.40" and 1.41" respectively. Spectroscopic and optical/near-infrared imaging follow-up observations show that the quasar images have identical redshifts and possess extended objects between the images that are likely to be lens galaxies at z_l~0.6 in SDSS J0806+2006 and z_l~0.3 in SDSS J1353+1138. The field of SDSS J0806+2006 contains several nearby galaxies that may significantly perturb the system, and SDSS J1353+1138 has an extra component near its Einstein ring that is probably a foreground star. Simple mass models with reasonable parameters reproduce the quasar positions and fluxes of both systems.

preprint2005arXiv

SDSSJ102111.02+491330.4: A Newly Discovered Gravitationally Lensed Quasar

We report follow-up observations of two gravitational lens candidates identified in the Sloan Digital Sky Survey (SDSS) dataset. We have confirmed that SDSS J102111.02+491330.4 is a previously unknown gravitationally lensed quasar. This lens system exhibits two images of a $z = 1.72$ quasar, with an image separation of $1{\farcs}14 \pm 0.04$. Optical and near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Observations of SDSS J112012.12+671116.0 indicate that it is more likely a binary quasar than a gravitational lens. This system has two quasars at a redshift of $z = 1.49$, with an angular separation of $1{\farcs}49 \pm 0.02$. However, the two quasars have markedly different SEDs and no lens galaxy is apparent in optical and near-IR images of this system. We also present a list of 31 SDSS lens candidates which follow-up observations have confirmed are \textit{not} gravitational lenses.

preprint2005arXiv

The UV, Optical, and IR Properties of SDSS Sources Detected by GALEX

We discuss the UV, optical, and IR properties of the SDSS sources detected by GALEX as part of its All-sky Imaging Survey Early Release Observations. Virtually all of the GALEX sources in the overlap region are detected by SDSS. GALEX sources represent ~2.5% of all SDSS sources within these fields and about half are optically unresolved. Most unresolved GALEX/SDSS sources are bright blue turn-off thick disk stars and are typically detected only in the GALEX near-UV band. The remaining unresolved sources include low-redshift quasars, white dwarfs, and white dwarf/M dwarf pairs, and these dominate the optically unresolved sources detected in both GALEX bands. Almost all the resolved SDSS sources detected by GALEX are fainter than the SDSS 'main' spectroscopic limit. These sources have colors consistent with those of blue (spiral) galaxies (u-r<2.2), and most are detected in both GALEX bands. Measurements of their UV colors allow much more accurate and robust estimates of star-formation history than are possible using only SDSS data. Indeed, galaxies with the most recent (<20 Myr) star formation can be robustly selected from the GALEX data by requiring that they be brighter in the far-UV than in the near-UV band. However, older starburst galaxies have UV colors similar to AGN, and thus cannot be selected unambiguously on the basis of GALEX fluxes alone. With the aid of 2MASS data, we construct and discuss median 10 band UV-optical-IR spectral energy distributions for turn-off stars, hot white dwarfs, low-redshift quasars, and spiral and elliptical galaxies. We point out the high degree of correlation between the UV color and the contribution of the UV flux to the UV-optical-IR flux of galaxies detected by GALEX.

preprint2004arXiv

Microlensing of the Broad Emission Line Region in the Quadruple Lens SDSS J1004+4112

We present seven epochs of spectroscopy on the quadruply imaged quasar SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. The spectra reveal differences in the emission lines between the lensed images. Specifically, component A showed a strong enhancement in the blue wings of several high-ionization lines relative to component B, which lasted at least 28 days (observed frame) then faded. Since the predicted time delay between A and B is <30 days, our time coverage suggests that the event was not intrinsic to the quasar. We attribute these variations to microlensing of part of the broad emission line region of the quasar, apparently resolving structure in the source plane on a scale of ~10^{16} cm at z=1.734. In addition, we observed smaller differences in the emission line profiles between components A and B that persisted throughout the time span, which may also be due to microlensing or millilensing. Further spectroscopic monitoring of this system holds considerable promise for resolving the structure of the broad emission line region in quasars.

preprint2004arXiv

SDSS J1335+0118: A New Two-Image Gravitational Lens

We report the discovery of the two-image gravitationally lensed quasar SDSS J1335+0118. The object was selected as a lens candidate from the Sloan Digital Sky Survey. The imaging and spectroscopic follow-up observations confirm that the system exhibits two gravitationally lensed images of a quasar at z=1.57. The image separation is 1.56''. We also detect an extended component between the two quasar images, likely the lensing galaxy. Preliminary mass modeling predicts the differential time delay Δt ~ 30 h^{-1} day assuming the redshift of the lens galaxy is 0.5.

preprint2003arXiv

A Snapshot Survey for Gravitational Lenses Among z>=4.0 Quasars: I. The z>5.7 Sample

Over the last few years, the Sloan Digital Sky Survey (SDSS) has discovered several hundred quasars with redshift between 4.0 and 6.4. Including the effects of magnification bias, one expects a priori that an appreciable fraction of these objects are gravitationally lensed. We have used the Advanced Camera for Surveys on the Hubble Space Telescope to carry out a snapshot imaging survey of high-redshift SDSS quasars to search for gravitationally split lenses. This paper, the first in a series reporting the results of the survey, describes snapshot observations of four quasars at z = 5.74, 5.82, 5.99 and 6.30, respectively. We find that none of these objects has a lensed companion within 5 magnitudes with a separation larger than 0.3 arcseconds; within 2.5 magnitudes, we can rule out companions within 0.1 arcseconds. Based on the non-detection of strong lensing in these four systems, we constrain the z~6 luminosity function to a slope of beta>-4.63 (3 sigma), assuming a break in the quasar luminosity function at M_{1450}^*=-24.0. We discuss the implications of this constraint on the ionizing background due to quasars in the early universe. Given that these quasars are not highly magnified, estimates of the masses of their central engines by the Eddington argument must be taken seriously, possibly challenging models of black hole formation.

preprint2003arXiv

SDSS J092455.87+021924.9: an Interesting Gravitationally Lensed Quasar from the Sloan Digital Sky Survey

We report the discovery of a new gravitationally lensed quasar from the Sloan Digital Sky Survey, SDSS J092455.87+021924.9 (SDSS J0924+0219). This object was selected from among known SDSS quasars by an algorithm that was designed to select another known SDSS lensed quasar (SDSS 1226-0006A,B). Five separate components, three of which are unresolved, are identified in photometric follow-up observations obtained with the Magellan Consortium's 6.5m Walter Baade telescope at Las Campanas Observatory. Two of the unresolved components (designated A and B) are confirmed to be quasars with z=1.524; the velocity difference is less than 100 km sec^{-1} according to spectra taken with the W. M. Keck Observatory's Keck II telescope on Mauna Kea. A third stellar component, designated C, has the colors of a quasar with redshift similar to components A and B. The maximum separation of the point sources is 1.78". The other two sources, designated G and D, are resolved. Component G appears to be the best candidate for the lensing galaxy. Although component D is near the expected position of the fourth lensed component in a four image lens system, its properties are not consistent with being the image of a quasar at z~1.5. Nevertheless, the identical redshifts of components A and B and the presence of component C strongly suggest that this object is a gravitational lens. Our observations support the idea that a foreground object reddens the fourth lensed component and that another unmodeled effect (such as micro- or milli-lensing) demagnificates it, but we cannot rule out the possibility that SDSS0924+0219 is an example of the relatively rare class of ``three component'' lens systems.

preprint2003arXiv

SDSS J115517.35+634622.0: A Newly Discovered Gravitationally Lensed Quasar

We report the discovery of SDSSJ115517.35+634622.0, a previously unknown gravitationally lensed quasar. The lens system exhibits two images of a $z = 2.89$ quasar, with an image separation of $1{\farcs}832 \pm 0.007$ . Near-IR imaging of the system reveals the presence of the lensing galaxy between the two quasar images. Based on absorption features seen in the Sloan Digital Sky Survey (SDSS) spectrum, we determine a lens galaxy redshift of $z = 0.1756$. The lens is rather unusual in that one of the quasar images is only $0{\farcs}22\pm0{\farcs}07$ ($\sim 0.1 R_{\rm eff}$) from the center of the lens galaxy and photometric modeling indicates that this image is significantly brighter than predicted by a SIS model. This system was discovered in the course of an ongoing search for strongly lensed quasars in the dataset from the SDSS.