We present reverberation mapping results from the first year of combined spectroscopic and photometric observations of the Sloan Digital Sky Survey Reverberation Mapping Project. We successfully recover reverberation time delays between the $g+i$-band emission and the broad H$\beta$ emission line for a total of 44 quasars, and for the broad H$\alpha$ emission line in 18 quasars. Time delays are computed using the JAVELIN and CREAM software and the traditional interpolated cross-correlation function (ICCF): Using well defined criteria, we report measurements of 32 H$\beta$ and 13 H$\alpha$ lags with JAVELIN, 42 H$\beta$ and 17 H$\alpha$ lags with CREAM, and 16 H$\beta$ and 8 H$\alpha$ lags with the ICCF. Lag values are generally consistent among the three methods, though we typically measure smaller uncertainties with JAVELIN and CREAM than with the ICCF, given the more physically motivated light curve interpolation and more robust statistical modeling of the former two methods. The median redshift of our H$\beta$-detected sample of quasars is 0.53, significantly higher than that of the previous reverberation-mapping sample. We find that in most objects, the time delay of the H$\alpha$ emission is consistent with or slightly longer than that of H$\beta$. We measure black hole masses using our measured time delays and line widths for these quasars. These black hole mass measurements are mostly consistent with expectations based on the local M-sigma relationship, and are also consistent with single-epoch black hole mass measurements. This work increases the current sample size of reverberation-mapped active galaxies by about two-thirds and represents the first large sample of reverberation mapping observations beyond the local universe (z < 0.3).

Quasar accretion-disk outflows might play an important role in galaxy evolution, but they are notoriously difficult to study due to line saturation and blending problems in the Ly-alpha forest. We circumvent these problems by constructing median composite spectra of diverse broad absorption lines (BALs) and `mini-BALs' in SDSS-III BOSS quasars at redshifts 2.3 < z < 3.5. Sorting by CIV 1549,1551 absorption-line strength with AlIII 1855,1863 as an additional indicator of low ionisations (LoBALs) we find the following: (1) Deeper and broader BALs are accompanied by weaker HeII 1640 emission lines, consistent with softer ionising spectra producing more effective radiative acceleration. (2) PV 1118,1128 absorption is present with resolved ~1:1 depth ratios in all composites from mini-BALs to strong BALs indicating that line saturation, large total column densities log N_H(cm^-2) > 22.7, and large ionisation parameters log U > -0.5 are common. (3) Different observed depths in saturated lines identify inhomogeneous partial covering on spatial scales <0.006 pc, where weak/low-abundance transitions like PV form in small high-column density clumps while stronger/broader lines like CIV form in larger volumes. (4) Excited-state SiIV* 1073 and CIII* 1176 lines in BAL outflows indicate typical densities n_e > 3 x 10^5 cm^-3 and maximum radial distances R < 23pc from the quasars. (5) For reasonable actual distances, the median BAL outflow has minimum kinetic energy L_K/L > 0.005(R/1.2pc) sufficient (by some estimates) for feedback to galaxy evolution. (6) LoBAL quasars have the largest median outflow column densities, highest velocities, and weakest HeII 1640 emission in our study; they appear to be at one extreme in a distribution of quasar properties where softer ionising spectra drive more powerful outflows.

Spectropolarimetry is a powerful technique that has provided critical support for the geometric unification model of local active galactic nuclei. In this paper, we present optical (rest-frame UV) Keck spectropolarimetry of five luminous obscured (Type 2) and extremely red quasars (ERQs) at z~2.5. Three objects reach polarization fractions of >10% in the continuum. We propose a model in which dust scattering is the dominant scattering and polarization mechanism in our targets, though electron scattering cannot be completely excluded. Emission lines are polarized at a lower level than is the continuum. This suggests that the emission-line region exists on similar spatial scales as the scattering region. In three objects we detect an intriguing 90 degree swing in the polarization position angle as a function of line-of-sight velocity in the emission lines of Ly-alpha, CIV and NV. We interpret this phenomenon in the framework of a geometric model with an equatorial dusty scattering region in which the material is outflowing at several thousand km/sec. Emission lines may also be scattered by dust or resonantly. This model explains several salient features of observations by scattering on scales of a few tens of pc. Our observations provide a tantalizing view of the inner region geometry and kinematics of high-redshift obscured and extremely red quasars. Our data and modeling lend strong support for toroidal obscuration and powerful outflows on the scales of the UV emission-line region, in addition to the larger scale outflows inferred previously from the optical emission-line kinematics.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Quasars may have played a key role in limiting the stellar mass of massive galaxies. Identifying those quasars in the process of removing star formation fuel from their hosts is an exciting ongoing challenge in extragalactic astronomy. In this paper we present X-ray observations of eleven extremely red quasars (ERQs) with $L_{\rm bol}\sim 10^{47}$ erg s$^{-1}$ at $z=1.5-3.2$ with evidence for high-velocity ($v > 1000$ km s$^{-1}$) [OIII]$\lambda$5007\AA\ outflows. X-rays allow us to directly probe circumnuclear obscuration and to measure the instantaneous accretion luminosity. We detect ten out of eleven extremely red quasars available in targeted and archival data. Using a combination of X-ray spectral fitting and hardness ratios, we find that all of the ERQs show signs of absorption in the X-rays with inferred column densities of $N_{\rm H}\approx 10^{23}$ cm$^{-2}$, including four Compton-thick candidates ($N_{\rm H} > 10^{24}$ cm$^{-2}$). We stack the X-ray emission of the seven weakly detected sources, measuring an average column density of $N_{\rm H}\sim 8\times 10^{23}$ cm$^{-2}$. The absorption-corrected (intrinsic) $2-10$ keV X-ray luminosity of the stack is $2.7\times 10^{45}$ erg s$^{-1}$, consistent with X-ray luminosities of type 1 quasars of the same infrared luminosity. Thus, we find that ERQs are a highly obscured, borderline Compton-thick population, and based on optical and infrared data we suggest that these objects are partially hidden by their own equatorial outflows. However, unlike some quasars with known outflows, ERQs do not appear to be intrinsically underluminous in X-rays for their bolometric luminosity. Our observations indicate that low X-rays are not necessary to enable some types of radiatively driven winds.

We present the clustering measurements of quasars in configuration space based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach with a Gaussian Streaming (GS) model for the redshift space distortions of the correlation function and demonstrate its applicability for dark matter halos hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$, we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426 \pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{-13}(r_{s}^{\rm fid}/r_s){\rm km.s}^{-1}.{\rm Mpc}^{-1}$, and the angular diameter distance $D_{A}(z_{\rm eff})=1850^{+90}_{-115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where $r_{s}$ is the sound horizon at the end of the baryon drag epoch and $r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors include both systematic and statistical contributions. The results on the evolution of distances are consistent with the predictions of flat $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmology with Planck parameters, and the measurement of $f\sigma_{8}$ extends the validity of General Relativity (GR) to higher redshifts($z>1$) This paper is released with companion papers using the same sample. The results on the cosmological parameters of the studies are found to be in very good agreement, providing clear evidence of the complementarity and of the robustness of the first full-shape clustering measurements with the eBOSS DR14 quasar sample.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

In these proceedings we describe the WEAVE-QSO survey, which will observe around 400,000 high redshift quasars starting in 2018. This survey is part of a broader WEAVE survey to be conducted at the 4.2m William Herschel Telescope. We will focus on chiefly on the science goals, but will also briefly summarise the target selection methods anticipated and the expected survey plan. Understanding the apparent acceleration in the expansion of the Universe is one of the key scientific challenges of our time. Many experiments have been proposed to study this expansion, using a variety of techniques. Here we describe a survey that can measure this acceleration and therefore help elucidate the nature of dark energy: a survey of the Lyman-alpha forest (and quasar absorption in general) in spectra towards z>2 quasars (QSOs). Further constraints on neutrino masses and warm dark matter are also anticipated. The same data will also shed light on galaxy formation via study of the properties of inflowing/outflowing gas associated with nearby galaxies and in a cosmic web context. Gas properties are sensitive to density, temperature, UV radiation, metallicity and abundance pattern, and so constraint galaxy formation in a variety of ways. WEAVE-QSO will study absorbers with a dynamic range spanning more than 8 orders of magnitude in column density, their thermal broadening, and a host of elements and ionization species. A core principal of the WEAVE-QSO survey is the targeting of QSOs with near 100% efficiency principally through use of the J-PAS (r < 23.2) and Gaia (r < 20) data.

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.

We describe the execution and data reduction of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high-redshift universe with VLT/XSHOOTER" (hereafter `XQ-100'). XQ-100 has produced and made publicly available an homogeneous and high-quality sample of echelle spectra of 100 QSOs at redshifts z~3.5-4.5 observed with full spectral coverage from 315 to 2500 nm at a resolving power ranging from R~4000 to 7000, depending on wavelength. The median signal-to-noise ratios are 33, 25 and 43, as measured at rest-frame wavelengths 1700, 3000 and 3600 Angstrom, respectively. This paper provides future users of XQ-100 data with the basic statistics of the survey, along with details of target selection, data acquisition and data reduction. The paper accompanies the public release of all data products, including 100 reduced spectra. XQ-100 is the largest spectroscopic survey to date of high-redshift QSOs with simultaneous rest-frame UV/optical coverage, and as such enables a wide range of extragalactic research, from cosmology and galaxy evolution to AGN astrophysics.

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.

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We measure the two-point clustering of spectroscopically confirmed quasars from the final sample of the Baryon Oscillation Spectroscopic Survey (BOSS) on comoving scales of 4 < s < 22 Mpc/h. The sample covers 6950 deg^2 (~ 19 (Gpc/h)^3) and, over the redshift range 2.2 < z < 2.8, contains 55,826 homogeneously selected quasars, which is twice as many as in any similar work. We deduce b_Q = 3.54 +/- 0.10 ; the most precise measurement of quasar bias to date at these redshifts. This corresponds to a host halo mass of ~ 2 x 10^12 ~ M_sun/h with an implied quasar duty cycle of ~1 percent. The real-space projected correlation function is well-fit by a power law of index -2 and correlation length r0 = (8.12 +/- 0.22), Mpc/h over scales of 4 < rp < 25 ~ Mpc/h. To better study the evolution of quasar clustering at moderate redshift, we extend the redshift range of our study to z ~ 3.4 and measure the bias and correlation length of three subsamples over 2.2 < z < 3.4. We find no significant evolution of r0 or bias over this range, implying that the host halo mass of quasars decreases somewhat with increasing redshift. We find quasar clustering remains similar over a decade in luminosity, contradicting a scenario in which quasar luminosity is monotonically related to halo mass at z ~ 2.5. Our results are broadly consistent with previous BOSS measurements, but they yield more precise constraints based upon a larger and more uniform data set.

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

In merger-driven models of massive galaxy evolution, the luminous quasar phase is expected to be accompanied by vigorous star formation in quasar host galaxies. In this paper, we use high column density Damped Lyman Alpha (DLA) systems along quasar sight lines as natural coronagraphs to directly study the far-UV (FUV) radiation from the host galaxies of luminous background quasars. We have stacked the spectra of $\sim$2,000 DLA systems (N_HI>10^{20.6} cm^{-2}) with a median absorption redshift < z > = 2.6 selected from quasars observed in the SDSS-III Baryon Oscillation Spectroscopic Survey. We detect residual flux in the dark troughs of the composite DLA spectra. The level of this residual flux significantly exceeds systematic errors in the SDSS fiber sky subtraction; furthermore, the residual flux is strongly correlated with the continuum luminosity of the background quasar, while uncorrelated with DLA column density or metallicity. We conclude that the flux could be associated with the average FUV radiation from the background quasar host galaxies (with medium redshift < z > = 3.1) that is not blocked by the intervening DLA. Assuming all of the detected flux originates from quasar hosts, for the highest quasar luminosity bin (< L >= 2.5x 10^{13} L_sun), the host galaxy has a FUV intensity of 1.5 +/- 0.2 x 10^{40} erg s^{-1} A^{-1}; this corresponds to an unobscured UV star formation rate of 9 M_sun/yr.

We consider how the profile and multi-year variability properties of a large sample of C IV Broad Absorption Line (BAL) troughs change when BALs from Si IV and/or Al III are present at corresponding velocities, indicating that the line-of-sight intercepts at least some lower ionization gas. We derive a number of observational results for C IV BALs separated according to the presence or absence of accompanying lower ionization transitions, including measurements of composite profile shapes, equivalent width (EW), characteristic velocities, composite variation profiles, and EW variability. We also measure the correlations between EW and fractional-EW variability for C IV, Si IV, and Al III. Our measurements reveal the basic correlated changes between ionization level, kinematics, and column density expected in accretion-disk wind models; e.g., lines-of-sight including lower ionization material generally show deeper and broader C IV troughs that have smaller minimum velocities and that are less variable. Many C IV BALs with no accompanying Si IV or Al III BALs may have only mild or no saturation.

We present composite spectra constructed from a sample of 242,150 Lyman-alpha (Lya) forest absorbers at redshifts 2.4<z<3.1 identified in quasar spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) as part of Data Release 9 of the Sloan Digital Sky Survey III. We select forest absorbers by their flux in bins 138 km/s wide (approximately the size of the BOSS resolution element). We split these absorbers into five samples spanning the range of flux -0.05 < F<0.45. Tests on a smaller sample of high-resolution spectra show that our three strongest absorption bins would probe circumgalactic regions (projected separation < 300 proper kpc and |Delta v| < 300km/s) in about 60% of cases for very high signal-to-noise ratio. Within this subset, weakening Lya absorption is associated with decreasing purity of circumgalactic selection once BOSS noise is included. Our weaker two Lya absorption samples are dominated by the intergalactic medium. We present composite spectra of these samples and a catalogue of measured absorption features from HI and 13 metal ionization species, all of which we make available to the community. We compare measurements of seven Lyman series transitions in our composite spectra to single line models and obtain further constraints from their associated excess Lyman limit opacity. This analysis provides results consistent with column densities over the range 14.4 <~ Log (N_HI) <~ 16.45. We compare our measurements of metal absorption to a variety of simple single-line, single-phase models for a preliminary interpretation. Our results imply clumping on scales down to ~30 pc and near-solar metallicities in the circumgalactic samples, while high-ionization metal absorption consistent with typical IGM densities and metallicities is visible in all samples.

We present a study of ~100 high redshift (z~2-4) extremely strong damped Lyman-alpha systems (ESDLA, with N(HI)>0.5x10^22 cm^-2) detected in quasar spectra from the Baryon Oscillation Spectroscopic Survey SDSS-III DR11. We study the neutral hydrogen, metal, and dust content of this elusive population of absorbers and confirm our previous finding that the high column density end of the N(HI) frequency distribution has a relatively shallow slope with power-law index -3.6, similar to what is seen from 21-cm maps in nearby galaxies. The stacked absorption spectrum indicates a typical metallicity ~1/20th solar, similar to the mean metallicity of the overall DLA population. The relatively small velocity extent of the low-ionisation lines suggests that ESDLAs do not arise from large-scale flows of neutral gas. The high column densities involved are in turn more similar to what is seen in DLAs associated with gamma-ray burst afterglows (GRB-DLAs), which are known to occur close to star forming regions. This indicates that ESDLAs arise from lines of sight passing at very small impact parameters from the host galaxy, as observed in nearby galaxies. This is also supported by simple theoretical considerations and recent high-z hydrodynamical simulations. We strongly substantiate this picture by the first statistical detection of Lya emission with <L>~(0.6+/-0.2)x10^42 erg/s in the core of ESDLAs (corresponding to about 0.1 L* at z~2-3), obtained through stacking the fibre spectra (of radius 1" corresponding to ~8 kpc at z~2.5). [truncated]

We present a detailed investigation of the variability of 428 C IV and 235 Si IV Broad Absorption Line (BAL) troughs identified in multi-epoch observations of 291 quasars by the Sloan Digital Sky Survey-I/II/III. These observations primarily sample rest-frame timescales of 1-3 yr over which significant rearrangement of the BAL wind is expected. We derive a number of observational results on, e.g., the frequency of BAL variability, the velocity range over which BAL variability occurs, the primary observed form of BAL-trough variability, the dependence of BAL variability upon timescale, the frequency of BAL strengthening vs. weakening, correlations between BAL variability and BAL-trough profiles, relations between C IV and Si IV BAL variability, coordinated multi-trough variability, and BAL variations as a function of quasar properties. We assess implications of these observational results for quasar winds. Our results support models where most BAL absorption is formed within an order-of-magnitude of the wind-launching radius, although a significant minority of BAL troughs may arise on larger scales. We estimate an average lifetime for a BAL trough along our line-of-sight of a few thousand years. BAL disappearance and emergence events appear to be extremes of general BAL variability, rather than being qualitatively distinct phenomena. We derive the parameters of a random-walk model for BAL EW variability, finding that this model can acceptably describe some key aspects of EW variability. The coordinated trough variability of BAL quasars with multiple troughs suggests that changes in "shielding gas" may play a significant role in driving general BAL variability.

We report the discovery in the Sloan Digital Sky Survey and the SDSS-III Baryon Oscillation Spectroscopic Survey of seventeen broad absorption line (BAL) quasars with high-ionization troughs that include absorption redshifted relative to the quasar rest frame. The redshifted troughs extend to velocities up to v=12,000 km/s and the trough widths exceed 3000 km/s in all but one case. Approximately 1 in 1000 BAL quasars with blueshifted C IV absorption also has redshifted C IV absorption; objects with C IV absorption present only at redshifted velocities are roughly four times rarer. In more than half of our objects, redshifted absorption is seen in C II or Al III as well as C IV, making low-ionization absorption at least ten times more common among BAL quasars with redshifted troughs than among standard BAL quasars. However, the C IV absorption equivalent widths in our objects are on average smaller than those of standard BAL quasars with low-ionization absorption. We consider several possible ways of generating redshifted absorption. The two most likely possibilities may be at work simultaneously, in the same objects or in different ones. Rotationally dominated outflows seen against a quasar's extended continuum source can produce redshifted and blueshifted absorption, but variability consistent with this scenario is seen in only one of the four objects with multiple spectra. The infall of relatively dense and low-ionization gas to radii as small as 400 Schwarzschild radii can in principle explain the observed range of trough profiles, but current models do not easily explain the origin and survival of such gas. Whatever the origin(s) of the absorbing gas in these objects, it must be located at small radii to explain its large redshifted velocities, and thus offers a novel probe of the inner regions of quasars.

We present a measurement of the quasar luminosity function in the range 0.68<z<4 down to extinction corrected magnitude g_dered=22.5, using a simple and well understood target selection technique based on the time-variability of quasars. The completeness of our sample was derived directly from a control sample of quasars, without requiring complex simulations of quasar light-curves or colors. A total of 1877 quasar spectra were obtained from dedicated programs on the Sloan telescope (as part of the SDSS-III/BOSS survey) and on the Multiple Mirror Telescope. They allowed us to derive the quasar luminosity function. It agrees well with previously published results in the redshift range 0.68<z<2.6. Our deeper data allow us to extend the measurement to z=4. We measured quasar densities to g_dered<22.5, obtaining 30 QSO per deg^2 at z<1, 99 QSO per deg^2 for 1<z<2.15, and 47 QSO per deg^2 at z>2.15. Using pure luminosity evolution models, we fitted our LF measurements, and predicted quasar number counts as a function of redshift and observed magnitude. These predictions are useful inputs for future cosmology surveys such as those relying on the observation of quasars to measure baryon acoustic oscillations.

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

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.

(abridged) We describe the automated spectral classification, redshift determination, and parameter measurement pipeline in use for the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III (SDSS-III) as of Data Release 9, encompassing 831,000 moderate-resolution optical spectra. We give a review of the algorithms employed, and describe the changes to the pipeline that have been implemented for BOSS relative to previous SDSS-I/II versions, including new sets of stellar, galaxy, and quasar redshift templates. For the color-selected CMASS sample of massive galaxies at redshift 0.4 <~ z <~ 0.8 targeted by BOSS for the purposes of large-scale cosmological measurements, the pipeline achieves an automated classification success rate of 98.7% and confirms 95.4% of unique CMASS targets as galaxies (with the balance being mostly M stars). Based on visual inspections of a subset of BOSS galaxies, we find that ~0.2% of confidently reported CMASS sample classifications and redshifts are incorrect, and ~0.4% of all CMASS spectra are objects unclassified by the current algorithm which are potentially recoverable. The BOSS pipeline confirms that ~51.5% of the quasar targets have quasar spectra, with the balance mainly consisting of stars. Statistical (as opposed to systematic) redshift errors propagated from photon noise are typically a few tens of km/s for both galaxies and quasars, with a significant tail to a few hundreds of km/s for quasars. We test the accuracy of these statistical redshift error estimates using repeat observations, finding them underestimated by a factor of 1.19 to 1.34 for galaxies, and by a factor of 2 for quasars. We assess the impact of sky-subtraction quality, S/N, and other factors on galaxy redshift success. Finally, we document known issues, and describe directions of ongoing development.

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

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)

The imprint of Baryonic Acoustic Oscillations (BAO) on the matter power spectrum can be constrained using the neutral hydrogen density in the intergalactic medium as a tracer of the matter density. One of the goals of the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS-III) is to derive the Hubble expansion rate and the angular scale from the BAO signal in the IGM. To this aim, the Lyman-alpha forest of 10^5 quasars will be observed in the redshift range 2.2<z<3.5 and over 10,000 deg^2. We simulated the BOSS QSO survey to estimate the statistical accuracy on the BAO scale determination provided by such a large scale survey. In particular, we discuss the effect of the poorly constrained estimate of the unabsorbed intrinsic quasar spectrum. The volume of current N-body simulations being too small for such studies, we resorted to Gaussian random field (GRF) simulations. We validated the use of GRFs by comparing the output of GRF simulations with that of the Horizon N-body simulation with the same initial conditions. Realistic mock samples of QSO Lyman-\alpha forest were generated; their power spectrum was computed and fitted to obtain the BAO scale. The rms of the results for 100 different simulations provides an estimate of the statistical error expected from the BOSS survey. We confirm the results from Fisher matrix estimate. In the absence of error on the unabsorbed quasar spectrum, the BOSS quasar survey should measure the BAO scale with an error of the order of 2.3%, or the transverse and radial BAO scales separately with errors of the order of 6.8% and 3.9%, respectively. The significance of the BAO detection is assessed by an average \Delta\chi^2=17 but for individual realizations \Delta\chi^2 ranges from 2 t o 35. The error on the unabsorbed quasar spectrum increases the error on the BAO scale by 10 to 20% and results in a sub percent bias.

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]

The SDSS-III BOSS Quasar survey will attempt to observe z>2.15 quasars at a density of at least 15 per square degree to yield the first measurement of the Baryon Acoustic Oscillations in the Ly-alpha forest. To help reaching this goal, we have developed a method to identify quasars based on their variability in the u g r i z optical bands. The method has been applied to the selection of quasar targets in the SDSS region known as Stripe 82 (the Southern equatorial stripe), where numerous photometric observations are available over a 10-year baseline. This area was observed by BOSS during September and October 2010. Only 8% of the objects selected via variability are not quasars, while 90% of the previously identified high-redshift quasar population is recovered. The method allows for a significant increase in the z>2.15 quasar density over previous strategies based on optical (ugriz) colors, achieving a density of 24.0 deg^{-2} on average down to g~22 over the 220 deg^2 area of Stripe 82. We applied this method to simulated data from the Palomar Transient Factory and from Pan-STARRS, and showed that even with data that have sparser time sampling than what is available in Stripe 82, including variability in future quasar selection strategies would lead to increased target selection efficiency in the z>2.15 redshift range. We also found that Broad Absorption Line quasars are preferentially present in a variability than in a color selection.

Baryonic Acoustic Oscillations (BAO) and their effects on the matter power spectrum can be studied by using the Lyman-alpha absorption signature of the matter density field along quasar (QSO) lines of sight. A measurement sufficiently accurate to provide useful cosmological constraints requires the observation of ~100000 quasars in the redshift range 2.2<z<3.5 over ~8000 deg2. Such a survey is planned by the Baryon Oscillation Spectroscopic Survey (BOSS) project of the Sloan Digital Sky Survey (SDSS-III).In practice, one needs a stellar rejection of more than two orders of magnitude with a selection efficiency for quasars better than 50% up to magnitudes as large as g ~ 22. To obtain an appropriate target list and estimate quasar redshifts, we have developed an Artificial Neural Networks (NN) with a multilayer perceptron architecture. The input variables are photometric measurements, i.e. the object magnitudes and their errors in the five bands (ugriz) of the SDSS photometry. For target selection, we achieve a non-quasar point-like object rejection of 99.6% and 98.5% for a quasar efficiency of, respectively, 50% and 85%. The photometric redshift precision is of the order of 0.1 over the region relevant for BAO studies.