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.