-
Most active galactic nuclei (AGNs) are radio-quiet, and the origin of their
radio emission is not well-understood. One hypothesis is that this radio
emission is a by-product of quasar-driven winds. In this paper, we present the
radio properties of 108 extremely red quasars (ERQs) at $z=2-4$. ERQs are among
the most luminous quasars ($L_{bol} \sim 10^{47-48}$ erg/s) in the Universe,
with signatures of extreme ($\gg 1000$ km/s) outflows in their
[OIII]$\lambda$5007 \AA\ emission, making them the best subjects to seek the
connection between radio and outflow activity. All ERQs but one are unresolved
in the radio on $\sim 10$ kpc scales, and the median radio luminosity of ERQs
is $\nu L_\nu [{\rm 6\,GHz}] = 10^{41.0}$ erg/s, in the radio-quiet regime, but
one to two orders of magnitude higher than that of other quasar samples. The
radio spectra are steep, with a mean spectral index $\langle \alpha \rangle =
-1.0$. In addition, ERQs neatly follow the extrapolation of the low-redshift
correlation between radio luminosity and the velocity dispersion of
[OIII]-emitting ionized gas. Uncollimated winds, with a power of one per cent
of the bolometric luminosity, can account for all these observations. Such
winds would interact with and shock the gas around the quasar and in the host
galaxy, resulting in acceleration of relativistic particles and the consequent
synchrotron emission observed in the radio. Our observations support the
picture in which ERQs are signposts of extremely powerful episodes of quasar
feedback, and quasar-driven winds as a contributor of the radio emission in the
intermediate regime of radio luminosity $\nu L_\nu = 10^{39}-10^{42}$ erg/s.
-
AGN exhibit rapid, high amplitude stochastic flux variations across the
entire electromagnetic spectrum on timescales ranging from hours to years. The
cause of this variability is poorly understood. We present a Green's
Function-based method for using variability to (1) measure the time-scales on
which flux perturbations evolve and (2) characterize the driving flux
perturbations. We model the observed light curve of an AGN as a linear
differential equation driven by stochastic impulses. We analyze the light curve
of the Kepler AGN Zw 229-15 and find that the observed variability behavior can
be modeled as a damped harmonic oscillator perturbed by a colored noise
process. The model powerspectrum turns over on time-scale $385$~d. On shorter
time-scales, the log-powerspectrum slope varies between $2$ and $4$, explaining
the behavior noted by previous studies. We recover and identify both the
$5.6$~d and $67$~d timescales reported by previous work using the Green's
Function of the C-ARMA equation rather than by directly fitting the
powerspectrum of the light curve. These are the timescales on which flux
perturbations grow, and on which flux perturbations decay back to the
steady-state flux level respectively. We make the software package KALI used to
study light curves using our method available to the community.
-
The CIV broad emission line is visible in optical spectra to redshifts
exceeding z~5. CIV has long been known to exhibit significant displacements to
the blue and these `blueshifts' almost certainly signal the presence of strong
outflows. As a consequence, single-epoch virial black hole (BH) mass estimates
derived from CIV velocity-widths are known to be systematically biased compared
to masses from the hydrogen Balmer lines. Using a large sample of 230
high-luminosity (log $L_{\rm Bol}$ = 45.5-48 erg/s), redshift 1.5<z<4.0 quasars
with both CIV and Balmer line spectra, we have quantified the bias in CIV BH
masses as a function of the CIV blueshift. CIV BH masses are shown to be a
factor of five larger than the corresponding Balmer-line masses at CIV
blueshifts of 3000 km/s and are over-estimated by almost an order of magnitude
at the most extreme blueshifts, >5000 km/s. Using the monotonically increasing
relationship between the CIV blueshift and the mass ratio BH(CIV)/BH(H$\alpha$)
we derive an empirical correction to all CIV BH-masses. The scatter between the
corrected CIV masses and the Balmer masses is 0.24 dex at low CIV blueshifts
(~0 km/s) and just 0.10 dex at high blueshifts (~3000 km/s), compared to 0.40
dex before the correction. The correction depends only on the CIV line
properties - i.e. full-width at half maximum and blueshift - and can therefore
be applied to all quasars where CIV emission line properties have been
measured, enabling the derivation of un-biased virial BH mass estimates for the
majority of high-luminosity, high-redshift, spectroscopically confirmed quasars
in the literature.
-
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.
-
Red quasars are candidate young objects in an early transition stage of
massive galaxy evolution. Our team recently discovered a population of
extremely red quasars (ERQs) in the Baryon Oscillation Spectroscopic Survey
(BOSS) that has a suite of peculiar emission-line properties including large
rest equivalent widths (REWs), unusual "wingless" line profiles, large NV/Lya,
NV/CIV, SiIV/CIV and other flux ratios, and very broad and blueshifted [OIII]
5007. Here we present a new catalog of CIV and NV emission-line data for
216,188 BOSS quasars to characterize the ERQ line properties further. We show
that they depend sharply on UV-to-mid-IR color, secondarily on REW(CIV), and
not at all on luminosity or the Baldwin Effect. We identify a "core" sample of
97 ERQs with nearly uniform peculiar properties selected via i-W3 > 4.6 (AB)
and REW(CIV) > 100 A at redshifts 2.0-3.4. A broader search finds 235 more red
quasars with similar unusual characteristics. The core ERQs have median
luminosity log L (ergs/s) ~ 47.1, sky density 0.010 deg^-2, surprisingly
flat/blue UV spectra given their red UV-to-mid-IR colors, and common outflow
signatures including BALs or BAL-like features and large CIV emission-line
blueshifts. Their SEDs and line properties are inconsistent with normal quasars
behind a dust reddening screen. We argue that the core ERQs are a unique
obscured quasar population with extreme physical conditions related to powerful
outflows across the line-forming regions. Patchy obscuration by small dusty
clouds could produce the observed UV extinctions without substantial UV
reddening.
-
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$\alpha$ 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$\beta$, or
H$\alpha$, 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.
-
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{\sigma} depths of 6.13 {\mu}Jy (21.93 AB magnitude)
and 5.75 {\mu}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 {\mu}m-only
detection catalog containing 6.1 million sources, a 4.5 {\mu}m-only detection
catalog containing 6.5 million sources, and a dual-band detection catalog
containing 5.4 million sources.
-
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.
-
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.
-
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.
-
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.
-
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 ($\Delta 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 $\delta$-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.
-
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.
-
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(\delta 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 ($\gamma = 2.90$ to within $10$ percent) for the re-processed
light curve as compared to the light curve found on MAST ($\gamma = 2.65$ to
within $10$ percent)---both inconsistent with a damped random walk which
requires $\gamma = 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.
-
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 $\gamma$ and turnover timescale $\tau$. 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 $\tau$ 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.
-
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 (${\alpha}_{\lambda}$) 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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.