The Legacy Surveys, a combination of three ground-based imaging surveys, have mapped 16,000 deg$^2$ in three optical bands ($g$, $r$, and $z$) to a depth 1--$2$~mag deeper than the Sloan Digital Sky Survey (SDSS). Our work addresses one of the major challenges of wide-field imaging surveys conducted at ground-based observatories: the varying depth that results from varying observing conditions at Earth-bound sites. To mitigate these effects, two of the Legacy Surveys (the Dark Energy Camera Legacy Survey, or DECaLS; and the Mayall $z$-band Legacy Survey, or MzLS) employed a unique strategy to dynamically adjust the exposure times as rapidly as possible in response to the changing observing conditions. We present the tiling and observing strategies used by these surveys. We demonstrate that the tiling and dynamic observing strategies jointly result in a more uniform-depth survey that has higher efficiency for a given total observing time compared with the traditional approach of using fixed exposure times.

The Lyman-$\alpha$ (hereafter Ly-$\alpha$) forest is a probe of large-scale matter density fluctuations at high redshift, $z > 2.1$. It consists of HI absorption spectra along individual lines-of-sight. If the line-of-sight density is large enough, 3D maps of HI absorption can be inferred by tomographic reconstruction. In this article, we investigate the Ly-$\alpha$ forest available in the Stripe 82 field ($220\,\mathrm{deg^{2}}$), based on the quasar spectra from SDSS Data Release DR16. The density of observed quasar spectra is $37\,\mathrm{deg^{-2}}$ with a mean pixel signal-to-noise ratio of two per angstrom. This study provides an intermediate case between the average SDSS density and that of the much denser but smaller CLAMATO survey. We derive a 3D map of large-scale matter fluctuations from these data, using a Wiener filter technique. The total volume of the map is $0.94\,\mathrm{h^{-3} Gpc^{3}}$. Its resolution is $13\,\mathrm{h^{-1} Mpc}$, which is related to the mean transverse distance between nearest lines-of-sight. From this map, we provide a catalog of voids and protocluster candidates in the cosmic web. The map-making and void catalog are compared to simulated eBOSS Stripe 82 observations. A stack over quasar positions provides a visualization of the Ly-$\alpha$ quasar cross-correlation. This tomographic reconstruction constitutes the largest-volume high-redshift 3D map of matter fluctuations.

Understanding the links between the activity of supermassive black holes (SMBH) at the centres of galaxies and their host dark matter haloes is a key question in modern astrophysics. The final data release of the SDSS-IV eBOSS provides the largest contemporary spectroscopic sample of galaxies and QSOs. Using this sample and covering the redshift interval $z=0.7-1.1$, we have measured the clustering properties of the eBOSS QSOs, Emission Line Galaxies (ELGs) and Luminous Red Galaxies (LRGs). We have also measured the fraction of QSOs as a function of the overdensity defined by the galaxy population. Using these measurements, we investigate how QSOs populate and sample the galaxy population, and how the host dark-matter haloes of QSOs sample the underlying halo distribution. We find that the probability of a galaxy hosting a QSO is independent of the host dark matter halo mass of the galaxy. We also find that about 60\% of eBOSS QSOs are hosted by LRGs and about 20-40\% of QSOs are hosted by satellite galaxies. We find a slight preference for QSOs to populate satellite galaxies over central galaxies. This is connected to the host halo mass distribution of different types of galaxies. Based on our analysis, QSOs should be hosted by a very broad distribution of haloes, and their occurrence should be modulated only by the efficiency of galaxy formation processes.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We search in the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys for new strong lensing systems by using deep residual neural networks, building on previous work presented in Huang et al. These surveys together cover approximately one third of the sky visible from the northern hemisphere, reaching a z-band magnitude of ~22.5. We compile a training sample that consists of known lensing systems as well as non-lenses in the Legacy Surveys and the Dark Energy Survey. After applying our trained neural networks to the survey data, we visually inspect and rank images with probabilities above a threshold. Here we present 1014 new strong lens candidates.

High-redshift quasars typically have their redshift determined from rest-frame ultraviolet (UV) emission lines. However, these lines, and more specifically the prominent C IV $\lambda 1549$ emission line, are typically blueshifted yielding highly uncertain redshift estimates compared to redshifts determined from rest-frame optical emission lines. We present near-infrared spectroscopy of 18 luminous quasars at $2.15 < z < 3.70$ that allows us to obtain reliable systemic redshifts for these sources. Together with near-infrared spectroscopy of an archival sample of 44 quasars with comparable luminosities and redshifts, we provide prescriptions for correcting UV-based redshifts. Our prescriptions reduce velocity offsets with respect to the systemic redshifts by $\sim140$ km s$^{-1}$ and reduce the uncertainty on the UV-based redshift by $\sim25\%$ with respect to the best method currently used for determining such values. We also find that the redshifts determined from the Sloan Digital Sky Survey Pipeline for our sources suffer from significant uncertainties, which cannot be easily mitigated. We discuss the potential of our prescriptions to improve UV-based redshift corrections given a much larger sample of high redshift quasars with near-infrared spectra.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

We combine the most precise small scale ($< 100\, \rm h^{-1}kpc$) quasar clustering constraintsto date with recent measurements at large scales ($> 1\, \rm h^{-1}Mpc$) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) to better constrain the satellite fraction of quasars at $z\sim 1.5$ in the halo occupation formalism. We build our Halo Occupation Distribution (HOD) framework based on commonly used analytic forms for the one and two-halo terms with two free parameters: the minimum halo mass that hosts a central quasar and the fraction of satellite quasars that are within one halo. Inspired by recent studies that propose a steeper density profile for the dark matter haloes that host quasars, we explore HOD models at kiloparsec scales and best-fit parameters for models with $10\times$ higher concentration parameter. We find that an HOD model with a satellite fraction of $f_{\rm sat} = 0.071_{-0.004}^{+0.009}$ and minimum mass of $\rm M_{m} = 2.31_{-0.38}^{+0.41} \times 10^{12}\, \, \rm h^{-1} M_{\odot}$ for the host dark matter haloes best describes quasar clustering (on all scales) at $z \sim 1.5$. Our results are marginally inconsistent with earlier work that studied brighter quasars, hinting at a luminosity-dependence to the one-halo term.

We examine multiple AGN systems (triples and quadruples, in particular) in the \texttt{MassiveBlackII} simulation over a redshift range of $0.06\lesssim z \lesssim 4$. We identify AGN systems (with bolometric luminosity $L_{\mathrm{bol}}>10^{42}~\mathrm{ergs/sec}$) at different scales~(defined by the maximum distance between member AGNs) to determine the AGN multiplicity functions. This is defined as the volume/ surface density of AGN systems per unit \textit{richness} $R$, the number of AGNs in a system. We find that gravitationally bound multiple AGN systems tend to populate scales of $\lesssim0.7~\mathrm{cMpc}/h$; this corresponds to angular separations of $\lesssim100~\mathrm{arcsec}$ and a line of sight velocity difference $\lesssim200~\mathrm{km/sec}$. The simulation contains $\sim 10$ and $\sim100$ triples/quadruples per $\mathrm{deg}^2$ up to depths of DESI ($g\lesssim24$) and LSST ($g\lesssim26$) imaging respectively; at least $20\%$ of these should be detectable in spectroscopic surveys. The simulated quasar ($L_{\mathrm{bol}}>10^{44}~\mathrm{ergs/sec}$) triples and quadruples predominantly exist at $1.5\lesssim z \lesssim 3$. Their members have black hole masses $10^{6.5}\lesssim M_{bh}\lesssim 10^{9}~M_{\odot}/h$ and live in separate (one central and multiple satellite) galaxies with stellar masses $10^{10}\lesssim M_{*}\lesssim 10^{12}~M_{\odot}/h$. They live in the most massive haloes (for e.g. $\sim 10^{13}~M_{\odot}/h$ at $z=2.5$; $\sim 10^{14}~M_{\odot}/h$ at $z=1$) in the simulation. Their detections provide an exciting prospect for understanding massive black hole growth and their merger rates in galaxies in the era of multi-messenger astronomy.

We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The LRG sample contains 2.7 million objects over the redshift range of $0.4 < z < 0.9$ over 5655 sq. degrees. We have developed new photometric redshift (photo-$z$) estimates using the Legacy Survey DECam and WISE photometry, with $\sigma_{\mathrm{NMAD}} = 0.02$ precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a 5-parameter HOD model fit using the MultiDark halo catalog, we find that there is little evolution in HOD parameters except at the highest-redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of MCMC fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogs, and present a new technique for finding best-fit parameters in this situation. Accompanying this paper we have released the PRLS (Photometric Redshifts for the Legacy Surveys) catalog of photo-$z$'s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 dataset. This catalog provides accurate photometric redshifts for objects with $z < 21$ over more than 16,000 square degrees of sky.

Clustering measurements of obscured and unobscured quasars show that obscured quasars reside in more massive dark matter halos than their unobscured counterparts. These results are inconsistent with simple unified (torus) scenarios, but might be explained by models in which the distribution of obscuring material depends on Eddington ratio or galaxy stellar mass. We test these possibilities by constructing simple physical models to compare to observed AGN populations. We find that previously observed relationships between obscuration and Eddington ratio or stellar mass are not sufficient reproduce the observed quasar clustering results ($\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.94 ^{+ 0.10}_{- 0.11}$ and $\langle \log M_{\text{halo}}/M_{\odot} \rangle = 12.49 ^{+ 0.08}_{- 0.08}$ for obscured and unobscured populations, respectively) while maintaining the observed fraction of obscured quasars (30-65$\%$). This work suggests that evolutionary models, in which obscuration evolves on the typical timescale for black hole growth, are necessary to understand the observed clustering of mid-IR selected quasars.

Stars that escape globular clusters form tidal tails that are predominantly shaped by the global distribution of mass in the Galaxy, but also preserve a historical record of small-scale perturbations. Using deep $grz$ photometry from DECaLS, we present highly probable members of the tidal tails associated with the disrupting globular cluster Palomar 5. These data yield the cleanest view of a stellar stream beyond $\sim20\,\rm kpc$ and reveal: (1) a wide, low surface-brightness extension of the leading tail; (2) significant density variations along the stream; and (3) sharp changes in the direction of both the leading and the trailing tail. In the fiducial Milky Way model, a rotating bar perturbs the Palomar 5 tails and can produce streams with similar width and density profiles to those observed. However, the deviations of the stream track in this simple model do not match those observed in the Palomar 5 trailing tail, indicating the need for an additional source of perturbation. These discoveries open up the possibility of measuring the population of perturbers in the Milky Way, including dark-matter subhalos, with an ensemble of stellar streams and deep photometry alone.

Merger simulations predict that tidally induced gas inflows can trigger kpc-scale dual active galactic nuclei (dAGN) in heavily obscured environments. Previously with the Very Large Array, we have confirmed four dAGN with redshifts between $0.04 < z < 0.22$ and projected separations between 4.3 and 9.2 kpc in the SDSS Stripe 82 field. Here, we present $Chandra$ X-ray observations that spatially resolve these dAGN and compare their multi-wavelength properties to those of single AGN from the literature. We detect X-ray emission from six of the individual merger components and obtain upper limits for the remaining two. Combined with previous radio and optical observations, we find that our dAGN have properties similar to nearby low-luminosity AGN, and they agree well with the black hole fundamental plane relation. There are three AGN-dominated X-ray sources, whose X-ray hardness-ratio derived column densities show that two are unobscured and one is obscured. The low obscured fraction suggests these dAGN are no more obscured than single AGN, in contrast to the predictions from simulations. These three sources show an apparent X-ray deficit compared to their mid-infrared continuum and optical [OIII] line luminosities, suggesting higher levels of obscuration, in tension with the hardness-ratio derived column densities. Enhanced mid-infrared and [OIII] luminosities from star formation may explain this deficit. There is ambiguity in the level of obscuration for the remaining five components since their hardness ratios may be affected by non-nuclear X-ray emissions, or are undetected altogether. They require further observations to be fully characterized.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

We present the first attempt at measuring the baryonic acoustic oscillations (BAO) in the large scale cross-correlation between the magnesium-II doublet (MgII) flux transmission field and the position of quasar and galaxy tracers. The MgII flux transmission continuous field at $0.3 < z < 1.6$ is measured from 500,589 quasar spectra obtained in the Baryonic Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS). The position of 246,697 quasar tracers and 1,346,776 galaxy tracers are extracted from the Sloan Digital Sky Survey (SDSS) I, II, BOSS and eBOSS catalogs. In addition to measuring the cosmological BAO scale and the biased matter density correlation, this study allows tests and improvements to cosmological Lyman-$\alpha$ analyses. A feature consistent with that of the BAO is detected at a significance of $\Delta \chi^{2} = 7.25$. The measured MgII linear transmission bias parameters are $b_{\mathrm{MgII(2796)}} (z = 0.59) = (-6.82 \pm 0.54) \, \times 10^{-4}$ and $b_{\mathrm{MgII(2804)}} (z = 0.59) = (-5.55 \pm 0.46) \, \times 10^{-4}$, and the MgI bias is $b_{\mathrm{MgI(2853)}} (z = 0.59) = (-1.48 \pm 0.24) \, \times 10^{-4}$. Their redshift evolution is characterized by the power-law index: $\gamma_{\mathrm{Mg}} = 3.36 \pm 0.46$. These measurements open a new window towards using BAO from flux transmission at $z < 1$ in the final eBOSS sample and in the upcoming sample from the Dark Energy Spectroscopic Instrument.

We perform a semi-automated search for strong gravitational lensing systems in the 9,000 deg$^2$ Dark Energy Camera Legacy Survey (DECaLS), part of the DESI Legacy Imaging Surveys (Dey et al.). The combination of the depth and breadth of these surveys are unparalleled at this time, making them particularly suitable for discovering new strong gravitational lensing systems. We adopt the deep residual neural network architecture (He et al.) developed by Lanusse et al. for the purpose of finding strong lenses in photometric surveys. We compile a training set that consists of known lensing systems in the Legacy Surveys and DES as well as non-lenses in the footprint of DECaLS. In this paper we show the results of applying our trained neural network to the cutout images centered on galaxies typed as ellipticals (Lang et al.) in DECaLS. The images that receive the highest scores (probabilities) are visually inspected and ranked. Here we present 335 candidate strong lensing systems, identified for the first time.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

We examine recent high-precision measurements of small-scale quasar clustering (at $z\sim0.5-2$ on scales of $\sim25~\mathrm{kpc/h}$) from the SDSS in the context of the MassiveBlackII (MBII) cosmological hydrodynamic simulation and conditional luminosity function (CLF) modeling. At these high luminosities ($g < 20.85$ quasars), the MBII simulation volume ($100~\mathrm{cMpc}/h$ comoving boxsize) has only 3 quasar pairs at distances of $1-4$ Mpc. The black-hole masses for the pairs range between $M_{bh}\sim1-3\times 10^{9}~M_{\odot}/h$ and the quasar hosts are haloes of $M_h\sim1-3\times10^{14}~M_{\odot}/h$. Such pairs show signs of recent major mergers in the MBII simulation. By modeling the central and satellite AGN CLFs as log-normal and Schechter distributions respectively (as seen in MBII AGNs), we arrive at CLF models which fit the simulation predictions and observed luminosity function and the small-scale clustering measured for the SDSS sample. The small-scale clustering of our mock quasars is well-explained by central-satellite quasar pairs that reside in $M_h>10^{14}~M_{\odot}/h$ dark matter haloes. For these pairs, satellite quasar luminosity is similar to that of central quasars. Our CLF models imply a relatively steep increase in the maximum satellite luminosity, $L^*_{\mathrm{sat}}$, in haloes of $M_h>10^{14}~M_{\odot}/h$ with associated larger values of $L^*_{\mathrm{sat}}$ at higher redshift. This leads to increase in the satellite fraction that manifests itself in an enhanced clustering signal at $\lesssim$ 1 Mpc/h. For the ongoing eBOSS-CORE sample, we predict $\sim 200-500$ quasar pairs at $z\sim1.5$ (with $M_h \gtrsim10^{13}~M_{\odot}/h$ and $M_{bh} \gtrsim10^{8}~M_{\odot}/h$) at $\sim25~\mathrm{kpc}$ scales. Such a sample would be $\gtrsim10$ times larger than current pair samples.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

We measure the average deflection of cosmic microwave background photons by quasars at $\langle z \rangle =1.7$. Our sample is selected from the Sloan Digital Sky Survey to cover the redshift range $0.9\leq z\leq2.2$ with absolute i-band magnitudes of $M_i\leq-24$ (K-corrected to z=2). A stack of nearly 200,000 targets reveals an 8$\sigma$ detection of Planck's estimate of the lensing convergence towards the quasars. We fit the signal with a model comprising a Navarro-Frenk-White density profile and a 2-halo term accounting for correlated large scale structure, which dominates the observed signal. The best-fitting model is described by an average halo mass $\log_{10}(M_{\rm h}/h^{-1}M_\odot)=12.6\pm0.2$ and linear bias $b=2.7\pm0.3$ at $z=1.7$, in excellent agreement with clustering studies. We also report of a hint, at a 90% confidence level, of a correlation between the convergence amplitude and luminosity, indicating that quasars brighter than $M_i\lesssim -26$ reside in halos of typical mass ${M_{\rm h}\approx 10^{13}\,h^{-1}M_\odot}$, scaling roughly as ${M_{\rm h}\propto L_{\rm opt}^{3/4}}$ at ${M_i\lesssim-24}$, in good agreement with physically-motivated quasar demography models. Although we acknowledge this luminosity dependence is a marginal result, the observed $M_{\rm h}$-$L_{\rm opt}$ relationship could be interpreted as a reflection of the cutoff in the distribution of black hole accretion rates towards high Eddington ratios: the weak trend of $M_{\rm h}$ with $L_{\rm opt}$ observed at low luminosity becomes stronger for the most powerful quasars, which tend to be accreting close to the Eddington limit.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We present new multi-broadband transit photometry of HD 189733b obtained at the Wyoming Infrared Observatory. Using an ensemble of five-band Sloan filter observations across multiple transits we have created an "ultra-low" resolution transmission spectrum to discern the nature of the exoplanet atmosphere. The observations were taken over three transit events and total 108 u', 120 g', 120 r', 110 i', and 116 z' images with an average exposure cadence of seven minutes for an entire series. The analysis was performed with a Markov-Chain Monte-Carlo method assisted by a Gaussian processes regression model. We find the apparent planet radius to increase from 0.154 +0.000920-0.00096 R* at z'-band to 0.157 +0.00074-0.00078 R* at u'-band. Whether this apparent radius implies an enhanced Rayleigh scattering or clear or grey planet atmosphere is highly dependent on stellar spot modeling assumptions, but our results are consistent with the literature for HD 189733b. This set of observations demonstrates the ability of our 2.3-m ground-based observatory to measure atmospheres of large exoplanets.

A primary aim of the Nuclear Spectroscopic Telescope Array (NuSTAR) mission is to find and characterize heavily obscured Active Galactic Nuclei (AGNs). Based on mid-infrared photometry from the Wide-Field Infrared Survey Explorer (WISE) and optical photometry from the Sloan Digital Sky Surveys, we have selected a large population of luminous obscured AGN (i.e., obscured quasars). Here we report NuSTAR observations of four WISE-selected heavily obscured quasars for which we have optical spectroscopy from the Southern African Large Telescope and W. M. Keck Observatory. Optical diagnostics confirm that all four targets are AGNs. With NuSTAR hard X-ray observations, three of the four objects are undetected, while the fourth has a marginal detection. We confirm that these objects have observed hard X-ray (10-40 keV) luminosities at or below ~ 10^43 erg s^-1. We compare X-ray and IR luminosities to obtain estimates of the hydrogen column densities (N_H) based on the suppression of the hard X-ray emission. We estimate N_H of these quasars to be at or larger than 10^25 cm^-2, confirming that WISE and optical selection can identify very heavily obscured quasars that may be missed in X-ray surveys, and do not contribute significantly to the cosmic X-ray background. From the optical Balmer decrements, we found that our three extreme obscured targets lie in highly reddened host environments. This galactic extinction is not adequate to explain the more obscured AGN, but it may imply a different scale of obscuration in the galaxy.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8<z<2.2$. We report a joint measurement of $f\sigma_8$ and two-dimensional BAO parameters $D_{\rm A}$ and $H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of dark energy is supported by observations at a $7.4\sigma$ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

We report on an X-ray and optical/UV study of eight Broad Absorption Line (BAL) to non-BAL transforming quasars at $z\,\approx\,$1.7-2.2 over 0.29-4.95 rest-frame years with at least three spectroscopic epochs for each quasar from the SDSS, BOSS, $Gemini$, and ARC 3.5-m telescopes. New $Chandra$ observations obtained for these objects show their values of $\alpha_{\rm ox}$ and $\Delta{\alpha}_{\rm ox}$, as well as their spectral energy distributions, are consistent with those of non-BAL quasars. Moreover, our targets have X-ray spectral shapes that are, on average, consistent with weakened absorption with an effective power-law photon index of $\Gamma_{\rm eff}\,=\,1.69^{+0.25}_{-0.25}$. The newer $Gemini$ and ARC 3.5-m spectra reveal that the BAL troughs have remained absent since the BOSS observations where the BAL disappearance was discovered. The X-ray and optical/UV results in tandem are consistent with at least the X-ray absorbing material moving out of the line-of-sight, leaving an X-ray unabsorbed non-BAL quasar. The UV absorber might have become more highly ionized (in a shielding-gas scenario) or also moved out of the line-of-sight (in a wind-clumping scenario).

This paper presents the second data release (DR2) of the Beijing-Arizona Sky Survey (BASS). BASS is an imaging survey of about 5400 deg$^2$ in $g$ and $r$ bands using the 2.3 m Bok telescope. DR2 includes the observations as of July 2017 obtained by BASS and Mayall $z$-band Legacy Survey (MzLS). This is our first time to include the MzLS data covering the same area as BASS. BASS and MzLS have respectively completed about 72% and 76% of their observations. The two surveys will be served for the spectroscopic targeting of the upcoming Dark Energy Spectroscopic Instrument. Both BASS and MzLS data are reduced by the same pipeline. We have updated the basic data reduction and photometric methods in DR2. In particular, source detections are performed on stacked images, and photometric measurements are co-added from single-epoch images based on these sources. The median 5$\sigma$ depths with corrections of the Galactic extinction are 24.05, 23.61, and 23.10 mag for $g$, $r$, and $z$ bands, respectively. The DR2 data products include stacked images, co-added catalogs, and single-epoch images and catalogs. The BASS website (http://batc.bao.ac.cn/BASS/) provides detailed information and links to download the data.

We explore the cosmological implications of anisotropic clustering measurements of the quasar sample from Data Release 14 of the Sloan Digital Sky Survey IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in configuration space. The $\sim 147,000$ quasar sample observed by eBOSS offers a direct tracer of the density field and bridges the gap of previous BAO measurements between redshift $0.8<z<2.2$. By analysing the two-point correlation function characterized by clustering wedges $\xi_{\rm w_i}(s)$ and multipoles $\xi_{\ell}(s)$, we measure the angular diameter distance, Hubble parameter and cosmic structure growth rate. We define a systematic error budget for our measurements based on the analysis of $N$-body simulations and mock catalogues. Based on the DR14 large scale structure quasar sample at the effective redshift $z_{\rm eff}=1.52$, we find the growth rate of cosmic structure $f\sigma_8(z_{\rm eff})=0.396\pm 0.079$, and the geometric parameters $D_{\rm V}(z)/r_{\rm d}=26.47\pm 1.23$, and $F_{\rm AP}(z)=2.53\pm 0.22$, where the uncertainties include both statistical and systematic errors. These values are in excellent agreement with the best-fitting standard ${\rm \Lambda CDM}$ model to the latest cosmic microwave background data from Planck.

We present a measurement of the two-point autocorrelation function of photometrically-selected, high-$z$ quasars over $\sim$ 100 deg$^2$ on the Sloan Digitial Sky Survey Stripe 82 field. Selection is performed using three machine-learning algorithms, trained on known high-$z$ quasar colors, in a six-dimensional, optical/mid-infrared color space. Optical data from the Sloan Digitial Sky Survey is combined with overlapping deep mid-infrared data from the \emph{Spitzer} IRAC Equatorial Survey and the \emph{Spitzer}-HETDEX Exploratory Large-area survey. The selected quasar sample consists of 1378 objects and contains both spectroscopically-confirmed quasars and photometrically-selected quasar candidates. These objects span a redshift range of $2.9 \leq z \leq 5.1$ and are generally fainter than $i=20.2$; a regime which has lacked sufficient number density to perform autocorrelation function measurements of photometrically-classified quasars. We compute the angular correlation function of these data, marginally detecting quasar clustering. We fit a single power-law with an index of $\delta = 1.39 \pm 0.618$ and amplitude of $\theta_0 = 0.71 \pm 0.546$ arcmin. A dark-matter model is fit to the angular correlation function to estimate the linear bias. At the average redshift of our survey ($\langle z \rangle = 3.38$) the bias is $b = 6.78 \pm 1.79$. Using this bias, we calculate a characteristic dark-matter halo mass of 1.70--9.83$\times 10^{12}h^{-1} M_{\odot}$. Our bias estimate suggests that quasar feedback intermittently shuts down the accretion of gas onto the central super-massive black hole at early times. If confirmed, these results hint at a level of luminosity dependence in the clustering of quasars at high-$z$.

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.

The rare close projection of two quasars on the sky provides the opportunity to study the host galaxy environment of a foreground quasar in absorption against the continuum emission of a background quasar. For over a decade the "Quasars probing quasars" series has utilized this technique to further the understanding of galaxy formation and evolution in the presence of a quasar at z>2, resolving scales as small as a galactic disc and from bound gas in the circumgalactic medium to the diffuse environs of intergalactic space. Presented here, is the public release of the quasar pair spectral database utilized in these studies. In addition to projected pairs at z>2, the database also includes quasar pair members at z<2, gravitational lens candidates and quasars closely separated in redshift that are useful for small-scale clustering studies. In total the database catalogs 5627 distinct objects, with 4083 lying within 5' of at least one other source. A spectral library contains 3582 optical and near-infrared spectra for 3028 of the cataloged sources. As well as reporting on 54 newly discovered quasar pairs, we outline the key contributions made by this series over the last ten years, summarize the imaging and spectroscopic data used for target selection, discuss the target selection methodologies, describe the database content and explore some avenues for future work. Full documentation for the spectral database, including download instructions, are supplied at the following URL http://specdb.readthedocs.io/en/latest/

We model the projected angular two-point correlation function (2PCF) of obscured and unobscured quasars selected using the Wide-field Infrared Survey Explorer (WISE), at a median redshift of $z \sim 1$ using a five-parameter Halo Occupation Distribution (HOD) parameterization, derived from a cosmological hydrodynamic simulation by Chatterjee et al. The HOD parameterization was previously used to model the 2PCF of optically selected quasars and X-ray bright active galactic nuclei (AGN) at $z \sim 1$. The current work shows that a single HOD parameterization can be used to model the population of different kinds of AGN in dark matter halos suggesting the universality of the relationship between AGN and their host dark matter halos. Our results show that the median halo mass of central quasar hosts increases from optically selected ($4.1^{+0.3}_{-0.4} \times 10^{12} \; h^{-1} \; {M_{sun}}$) and infra-red (IR) bright unobscured populations ($6.3^{+6.2}_{-2.3} \times 10^{12} \; h^{-1} \; {M_{sun}}$) to obscured quasars ($10.0^{+2.6}_{-3.7} \times 10^{12} \; h^{-1} \; {M_{sun}}$), signifying an increase in the degree of clustering. The projected satellite fractions also increase from optically bright to obscured quasars and tend to disfavor a simple `orientation only' theory of active galactic nuclei unification. Our results also show that future measurements of the small-scale clustering of obscured quasars can constrain current theories of galaxy evolution where quasars evolve from an IR- bright obscured phase to the optically bright unobscured phase.

We characterise the distribution of quasars within dark matter halos using a direct measurement technique for the first time at redshifts as high as $z \sim 1$. Using the Planck Sunyaev-Zeldovich (SZ) catalogue for galaxy groups and the Sloan Digital Sky Survey (SDSS) DR12 quasar dataset, we assign host clusters/groups to the quasars and make a measurement of the mean number of quasars within dark matter halos as a function of halo mass. We find that a simple power-law fit of $\log\left <N\right> = (2.11 \pm 0.01) \log (M) -(32.77 \pm 0.11)$ can be used to model the quasar fraction in dark matter halos. This suggests that the quasar fraction increases monotonically as a function of halo mass even to redshifts as high as $z\sim 1$.

We present a measurement of the anisotropic and isotropic Baryon Acoustic Oscillations (BAO) from the extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample with optimal redshift weights. Applying the redshift weights improves the constraint on the BAO dilation parameter $\alpha(z_{\rm eff})$ by 17\%. We reconstruct the evolution history of the BAO distance indicators in the redshift range of $0.8<z<2.2$. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

We present the Data Release 14 Quasar catalog (DR14Q) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). This catalog includes all SDSS-IV/eBOSS objects that were spectroscopically targeted as quasar candidates and that are confirmed as quasars via a new automated procedure combined with a partial visual inspection of spectra, have luminosities $M_{\rm i} \left[ z=2 \right] < -20.5$ (in a $\Lambda$CDM cosmology with $H_0 = 70 \ {\rm km \ s^{-1} \ Mpc ^{-1}}$, $\Omega_{\rm M} = 0.3$, and $\Omega_{\rm \Lambda} = 0.7$), and either display at least one emission line with a full width at half maximum (FWHM) larger than $500 \ {\rm km \ s^{-1}}$ or, if not, have interesting/complex absorption features. The catalog also includes previously spectroscopically-confirmed quasars from SDSS-I, II and III. The catalog contains 526,356 quasars 144,046 are new discoveries since the beginning of SDSS-IV) detected over 9,376 deg$^2$ (2,044 deg$^2$ having new spectroscopic data available) with robust identification and redshift measured by a combination of principal component eigenspectra. The catalog is estimated to have about 0.5% contamination. The catalog identifies 21,877 broad absorption line quasars and lists their characteristics. For each object, the catalog presents SDSS five-band CCD-based photometry with typical accuracy of 0.03 mag. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

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 analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample (DR14Q). We measure the redshift space distortions using the power spectrum monopole, quadrupole and hexadecapole inferred from 148,659 quasars between redshifts 0.8 and 2.2 covering a total sky footprint of 2112.9 deg$^2$. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f\sigma_8$, and the Alcock-Paczynski dilation scales which allow constraints to be placed on the angular diameter distance $D_A(z)$ and the Hubble $H(z)$ parameter. At the effective redshift of $z_{\rm eff}=1.52$, $f\sigma_8(z_{\rm eff})=0.420\pm0.076$, $H(z_{\rm eff})=[162\pm 12]\, (r_s^{\rm fid}/r_s)\,{\rm km\, s}^{-1}{\rm Mpc}^{-1}$, and $D_A(z_{\rm eff})=[1.85\pm 0.11]\times10^3\,(r_s/r_s^{\rm fid})\,{\rm Mpc}$, where $r_s$ is the comoving sound horizon at the baryon drag epoch and the superscript `fid' stands for its fiducial value. The errors take into account the full error budget, including systematics and statistical contributions. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model inferred from Planck measurements. Finally, we compare our measurements with other eBOSS companion papers and find excellent agreement, demonstrating the consistency and complementarity of the different methods used for analysing the data.

As astronomers increasingly exploit the information available in the time domain, spectroscopic variability in particular opens broad new channels of investigation. Here we describe the selection algorithms for all targets intended for repeat spectroscopy in the Time Domain Spectroscopic Survey (TDSS), part of the extended Baryon Oscillation Spectroscopic Survey within the Sloan Digital Sky Survey-IV. Also discussed are the scientific rationale and technical constraints leading to these target selections. The TDSS includes a large "Repeat Quasar Spectroscopy" (RQS) program delivering ~13,000 repeat spectra of confirmed SDSS quasars, and several smaller "Few-Epoch Spectroscopy" (FES) programs targeting specific classes of quasars as well as stars. The RQS program aims to provide a large and diverse quasar data set for studying variations in quasar spectra on timescales of years, a comparison sample for the FES quasar programs, and opportunity for discovering rare, serendipitous events. The FES programs cover a wide variety of phenomena in both quasars and stars. Quasar FES programs target broad absorption line quasars, high signal-to-noise ratio normal broad line quasars, quasars with double-peaked or very asymmetric broad emission line profiles, binary supermassive black hole candidates, and the most photometrically variable quasars. Strongly variable stars are also targeted for repeat spectroscopy, encompassing many types of eclipsing binary systems, and classical pulsators like RR Lyrae. Other stellar FES programs allow spectroscopic variability studies of active ultracool dwarf stars, dwarf carbon stars, and white dwarf/M dwarf spectroscopic binaries. We present example TDSS spectra and describe anticipated sample sizes and results.

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.

We present observed mid-infrared and optical colors and composite spectral energy distributions (SEDs) of type 1 (broad-line) and 2 (narrow-line) quasars selected from Sloan Digital Sky Survey (SDSS) spectroscopy. A significant fraction of powerful quasars are obscured by dust, and are difficult to detect in optical photometric or spectroscopic surveys. However these may be more easily identified on the basis of mid-infrared (MIR) colors and SEDs. Using samples of SDSS type 1 type 2 matched in redshift and [OIII] luminosity, we produce composite rest-frame 0.2-15 micron SEDs based on SDSS, UKIDSS, and Wide-Field Infrared Survey Explorer (WISE) photometry and perform model fits using simple galaxy and quasar SED templates. The SEDs of type 1 and 2 quasars are remarkably similar, with the differences explained primarily by the extinction of the quasar component in the type 2 systems. For both types of quasar, the flux of the AGN relative to the host galaxy increases with AGN luminosity (L_[OIII]) and redder observed MIR color, but we find only weak dependencies of the composite SEDs on mechanical jet power as determined through radio luminosity. We conclude that luminous quasars can be effectively selected using simple MIR color criteria similar to those identified previously (W1-W2 > 0.7 [Vega]), although these criteria miss many heavily obscured objects. Obscured quasars can be further identified based on optical-IR colors (for example, (u-W3 [AB]) > 1.4(W1-W2 [Vega])+3.2). These results illustrate the power of large statistical studies of obscured quasars selected on the basis of mid-IR and optical photometry.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Quasar redshifts require the best possible precision and accuracy for a number of applications, such as setting the velocity scale for outflows as well as measuring small-scale quasar-quasar clustering. The most reliable redshift standard in luminous quasars is arguably the narrow [OIII] $\lambda\lambda$4959,5007 emission line doublet in the rest-frame optical. We use previously published [OIII] redshifts obtained using near-infrared spectra in a sample of 45 high-redshift (z > 2.2) quasars to evaluate redshift measurement techniques based on rest-frame ultraviolet spectra. At redshifts above z = 2.2 the MgII $\lambda$2798 emission line is not available in observed-frame optical spectra, and the most prominent unblended and unabsorbed spectral feature available is usually CIV $\lambda$1549. Peak and centroid measurements of the CIV profile are often blueshifted relative to the rest-frame of the quasar, which can significantly bias redshift determinations. We show that redshift determinations for these high-redshift quasars are significantly correlated with the emission-line properties of CIV (i.e., the equivalent width, or EW, and the full width at half maximum, or FWHM) as well as the luminosity, which we take from the Sloan Digital Sky Survey Data Release 7. We demonstrate that empirical corrections based on multiple regression analyses yield significant improvements in both the precision and accuracy of the redshifts of the most distant quasars and are required to establish consistency with redshifts determined in more local quasars.

In this study, we probe the transition to cosmic homogeneity in the Large Scale Structure (LSS) of the Universe using the CMASS galaxy sample of BOSS spectroscopic survey which covers the largest effective volume to date, $3\ h^{-3}\ \mathrm{Gpc}^3$ at $0.43 \leq z \leq 0.7$. We study the scaled counts-in-spheres, $\mathcal{N}(<r)$, and the fractal correlation dimension, $\mathcal{D}_2(r)$, to assess the homogeneity scale of the universe using a $Landy\ \&\ Szalay$ inspired estimator. Defining the scale of transition to homogeneity as the scale at which $\mathcal{D}_2(r)$ reaches 3 within $1\%$, i.e. $\mathcal{D}_2(r)>2.97$ for $r>\mathcal{R}_H$, we find $\mathcal{R}_H = (63.3\pm0.7) \ h^{-1}\ \mathrm{Mpc}$, in agreement at the percentage level with the predictions of the $\Lambda$CDM model $\mathcal{R}_H=62.0\ h^{-1}\ \mathrm{Mpc}$. Thanks to the large cosmic depth of the survey, we investigate the redshift evolution of the transition to homogeneity scale and find agreement with the $\Lambda$CDM prediction. Finally, we find that $\mathcal{D}_2$ is compatible with $3$ at scales larger than $300\ h^{-1}\ $Mpc in all redshift bins. These results consolidate the Cosmological Principle and represent a precise consistency test of the $\Lambda CDM$ model.

Recent work has found evidence for a difference in the bias and dark matter halo masses of WISE-selected obscured and unobscured quasars, implying a distinction between these populations beyond random line-of-sight effects. However, the significance of this difference in the most up-to-date measurements is relatively weak, at $\sim$2$\sigma$ for individual measurements but bolstered by agreement from different techniques, including angular clustering and cross-correlations with cosmic microwave background (CMB) lensing maps. Here, we expand the footprint of previous work, aiming to improve the precision of both methods. In this larger area we correct for position dependent selection effects, in particular fluctuations of the WISE-selected quasar density as a function of Galactic latitude. We also measure the cross-correlation of the obscured and unobscured samples and confirm that they are well-matched in redshift, both centred at $z=1$. Combined with very similar detection fractions and magnitude distributions in the long-wavelength WISE bands, this redshift match strongly supports the fact that IR selection identifies obscured and unobscured quasars of similar bolometric luminosity. Finally, we perform cross-correlations with confirmed spectroscopic quasars, again confirming the results from other methods --- obscured quasars reside in haloes a factor of 3 times more massive than unobscured quasars. This difference is significant at the $\sim$5$\sigma$ level when the measurements are combined, strong support for the idea that obscuration in at least some quasars is tied to the larger environment, and may have an evolutionary component.

We study the first year of the eBOSS quasar sample in the redshift range $0.9<z<2.2$ which includes 68,772 homogeneously selected quasars. We show that the main source of systematics in the evaluation of the correlation function arises from inhomogeneities in the quasar target selection, particularly related to the extinction and depth of the imaging data used for targeting. We propose a weighting scheme that mitigates these systematics. We measure the quasar correlation function and provide the most accurate measurement to date of the quasar bias in this redshift range, $b_Q = 2.45 \pm 0.05$ at $\bar z=1.55$, together with its evolution with redshift. We use this information to determine the minimum mass of the halo hosting the quasars and the characteristic halo mass, which we find to be both independent of redshift within statistical error. Using a recently-measured quasar-luminosity-function we also determine the quasar duty cycle. The size of this first year sample is insufficient to detect any luminosity dependence to quasar clustering and this issue should be further studied with the final $\sim$500,000 eBOSS quasar sample.

SPIDERS (SPectroscopic IDentification of eROSITA Sources) is an SDSS-IV survey running in parallel to the eBOSS cosmology project. SPIDERS will obtain optical spectroscopy for large numbers of X-ray-selected AGN and galaxy cluster members detected in wide area eROSITA, XMM-Newton and ROSAT surveys. We describe the methods used to choose spectroscopic targets for two sub-programmes of SPIDERS: X-ray selected AGN candidates detected in the ROSAT All Sky and the XMM-Newton Slew surveys. We have exploited a Bayesian cross-matching algorithm, guided by priors based on mid-IR colour-magnitude information from the WISE survey, to select the most probable optical counterpart to each X-ray detection. We empirically demonstrate the high fidelity of our counterpart selection method using a reference sample of bright well-localised X-ray sources collated from XMM-Newton, Chandra and Swift-XRT serendipitous catalogues, and also by examining blank-sky locations. We describe the down-selection steps which resulted in the final set of SPIDERS-AGN targets put forward for spectroscopy within the eBOSS/TDSS/SPIDERS survey, and present catalogues of these targets. We also present catalogues of ~12000 ROSAT and ~1500 XMM-Newton Slew survey sources which have existing optical spectroscopy from SDSS-DR12, including the results of our visual inspections. On completion of the SPIDERS program, we expect to have collected homogeneous spectroscopic redshift information over a footprint of ~7500 deg$^2$ for >85 percent of the ROSAT and XMM-Newton Slew survey sources having optical counterparts in the magnitude range 17<r<22.5, producing a large and highly complete sample of bright X-ray-selected AGN suitable for statistical studies of AGN evolution and clustering.

We present the first discoveries from a survey of $z\gtrsim6$ quasars using imaging data from the DECam Legacy Survey (DECaLS) in the optical, the UKIRT Deep Infrared Sky Survey (UKIDSS) and a preliminary version of the UKIRT Hemisphere Survey (UHS) in the near-IR, and ALLWISE in the mid-IR. DECaLS will image 9000 deg$^2$ of sky down to $z_{\rm AB}\sim23.0$, and UKIDSS and UHS, which will map the northern sky at $0<DEC<+60^{\circ}$, reaching $J_{\rm VEGA}\sim19.6$ (5-$\sigma$). The combination of these datasets allows us to discover quasars at redshift $z\gtrsim7$ and to conduct a complete census of the faint quasar population at $z\gtrsim6$. In this paper, we report on the selection method of our search, and on the initial discoveries of two new, faint $z\gtrsim6$ quasars and one new $z=6.63$ quasar in our pilot spectroscopic observations. The two new $z\sim6$ quasars are at $z=6.07$ and $z=6.17$ with absolute magnitudes at rest-frame wavelength 1450 \AA\ being $M_{1450}=-25.83$ and $M_{1450}=-25.76$, respectively. These discoveries suggest that we can find quasars close to or fainter than the break magnitude of the Quasar Luminosity Function (QLF) at $z\gtrsim6$. The new $z=6.63$ quasar has an absolute magnitude of $M_{1450}=-25.95$. This demonstrates the potential of using the combined DECaLS and UKIDSS/UHS datasets to find $z\gtrsim7$ quasars. Extrapolating from previous QLF measurements, we predict that these combined datasets will yield $\sim200$ $z\sim6$ quasars to $z_{\rm AB} < 21.5$, $\sim1{,}000$ $z\sim6$ quasars to $z_{\rm AB}<23$, and $\sim 30$ quasars at $z>6.5$ to $J_{\rm VEGA}<19.5$.

We present the most precise estimate to date of the clustering of quasars on very small scales, based on a sample of 47 binary quasars with magnitudes of $g<20.85$ and proper transverse separations of $\sim 25\,h^{-1}$\,kpc. Our sample of binary quasars, which is about 6 times larger than any previous spectroscopically confirmed sample on these scales, is targeted using a Kernel Density Estimation technique (KDE) applied to Sloan Digital Sky Survey (SDSS) imaging over most of the SDSS area. Our sample is "complete" in that all of the KDE target pairs with $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc in our area of interest have been spectroscopically confirmed from a combination of previous surveys and our own long-slit observational campaign. We catalogue 230 candidate quasar pairs with angular separations of $<8\arcsec$, from which our binary quasars were identified. We determine the projected correlation function of quasars ($\bar W_{\rm p}$) in four bins of proper transverse scale over the range $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc. The implied small-scale quasar clustering amplitude from the projected correlation function, integrated across our entire redshift range, is $A=24.1\pm3.6$ at $\sim 26.6 ~h^{-1}$\,kpc. Our sample is the first spectroscopically confirmed sample of quasar pairs that is sufficiently large to study how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ kpc. We find that empirical descriptions of how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ Mpc also adequately describe the evolution of quasar clustering at $\sim 25 ~h^{-1}$ kpc.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We report on the diversity in quasar spectra from the Baryon Oscillation Spectroscopic Survey. After filtering the spectra to mitigate selection effects and Malmquist bias associated with a nearly flux-limited sample, we create high signal-to-noise ratio composite spectra from 58,656 quasars (2.1 \le z \le 3.5), binned by luminosity, spectral index, and redshift. With these composite spectra, we confirm the traditional Baldwin effect (BE, i.e., the anticorrelation of C IV equivalent width (EW) and luminosity) that follows the relation W_\lambda \propto L^{\beta_w} with slope \beta_w = -0.35 \pm 0.004, -0.35 \pm 0.005, and -0.41 \pm 0.005 for z = 2.25, 2.46, and 2.84, respectively. In addition to the redshift evolution in the slope of the BE, we find redshift evolution in average quasar spectral features at fixed luminosity. The spectroscopic signature of the redshift evolution is correlated at 98% with the signature of varying luminosity, indicating that they arise from the same physical mechanism. At a fixed luminosity, the average C IV FWHM decreases with increasing redshift and is anti-correlated with C IV EW. The spectroscopic signature associated with C IV FWHM suggests that the trends in luminosity and redshift are likely caused by a superposition of effects that are related to black hole mass and Eddington ratio. The redshift evolution is the consequence of a changing balance between these two quantities as quasars evolve toward a population with lower typical accretion rates at a given black hole mass.

The BOSS quasar sample is used to study cosmic homogeneity with a 3D survey in the redshift range $2.2<z<2.8$. We measure the count-in-sphere, $N(<\! r)$, i.e. the average number of objects around a given object, and its logarithmic derivative, the fractal correlation dimension, $D_2(r)$. For a homogeneous distribution $N(<\! r) \propto r^3$ and $D_2(r)=3$. Due to the uncertainty on tracer density evolution, 3D surveys can only probe homogeneity up to a redshift dependence, i.e. they probe so-called "spatial isotropy". Our data demonstrate spatial isotropy of the quasar distribution in the redshift range $2.2<z<2.8$ in a model-independent way, independent of any FLRW fiducial cosmology, resulting in $3-\langle D_2 \rangle < 1.7 \times 10^{-3}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc for the quasar distribution. If we assume that quasars do not have a bias much less than unity, this implies spatial isotropy of the matter distribution on large scales. Then, combining with the Copernican principle, we finally get homogeneity of the matter distribution on large scales. Alternatively, using a flat $\Lambda$CDM fiducial cosmology with CMB-derived parameters, and measuring the quasar bias relative to this $\Lambda$CDM model, our data provide a consistency check of the model, in terms of how homogeneous the Universe is on different scales. $D_2(r)$ is found to be compatible with our $\Lambda$CDM model on the whole $10<r<1200 \, h^{-1}$Mpc range. For the matter distribution we obtain $3-\langle D_2 \rangle < 5 \times 10^{-5}$ (2 $\sigma$) over the range $250<r<1200 \, h^{-1}$Mpc, consistent with homogeneity on large scales.

We present two wide-field catalogs of photometrically-selected emission line galaxies (ELGs) at z=0.8 covering about 2800 deg^2 over the south galactic cap. The catalogs were obtained using a Fisher discriminant technique described in a companion paper. The two catalogs differ by the imaging used to define the Fisher discriminant: the first catalog includes imaging from the Sloan Digital Sky Survey and the Wide-Field Infrared Survey Explorer, the second also includes information from the South Galactic Cap U-band Sky Survey (SCUSS). Containing respectively 560,045 and 615,601 objects, they represent the largest ELG catalogs available today and were designed for the ELG programme of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We study potential sources of systematic variation in the angular distribution of the selected ELGs due to fluctuations of the observational parameters. We model the influence of the observational parameters using a multivariate regression and implement a weighting scheme that allows effective removal of all of the systematic errors induced by the observational parameters. We show that fluctuations in the imaging zero-points of the photometric bands have minor impact on the angular distribution of objects in our catalogs. We compute the angular clustering of both catalogs and show that our weighting procedure effectively removes spurious clustering on large scales. We fit a model to the small scale angular clustering, showing that the selections have similar biases of 1.35/D_a(z) and 1.28/D_a(z). Both catalogs are publicly available.

We describe the redmonster automated redshift measurement and spectral classification software designed for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). We describe the algorithms, the template standard and requirements, and the newly developed galaxy templates to be used on eBOSS spectra. We present results from testing on early data from eBOSS, where we have found a 90.5% automated redshift and spectral classification success rate for the luminous red galaxy sample (redshifts 0.6$\lesssim z \lesssim$1.0). The redmonster performance meets the eBOSS cosmology requirements for redshift classification and catastrophic failures, and represents a significant improvement over the previous pipeline. We describe the empirical processes used to determine the optimum number of additive polynomial terms in our models and an acceptable $\Delta\chi_r^2$ threshold for declaring statistical confidence. Statistical errors on redshift measurement due to photon shot noise are assessed, and we find typical values of a few tens of km s$^{-1}$. An investigation of redshift differences in repeat observations scaled by error estimates yields a distribution with a Gaussian mean and standard deviation of $\mu\sim$0.01 and $\sigma\sim$0.65, respectively, suggesting the reported statistical redshift uncertainties are over-estimated by $\sim$54%. We assess the effects of object magnitude, signal-to-noise ratio, fiber number, and fiber head location on the pipeline's redshift success rate. Finally, we describe directions of ongoing development.

Recent measurements of the dark matter halo masses of infrared-selected obscured quasars are in tension --- some indicate that obscured quasars have higher halo mass compared to their unobscured counterparts, while others find no difference. The former result is inconsistent with the simplest models of quasar unification that rely solely on viewing angle, while the latter may support such models. Here, using empirical relationships between dark matter halo and supermassive black hole masses, we provide a simple evolutionary picture that naturally explains these findings and is motivated by more sophisticated merger-driven quasar fueling models. The model tracks the growth rate of haloes, with the black hole growing in spurts of quasar activity in order to "catch-up" with the $M_{\textrm{BH}}$ - $M_{\textrm{star}}$ - $M_{\textrm{halo}}$ relationship. The first part of the quasar phase is obscured and is followed by an unobscured phase. Depending on the luminosity limit of the sample, driven by observational selection effects, a difference in halo masses may or may not be significant. For high luminosity samples, the difference can be large (a few to 10 times higher masses in obscured quasars), while for lower luminosity samples the halo mass difference is very small, much smaller than current observational constraints. Such a simple model provides a qualitative explanation for the higher mass haloes of obscured quasars, as well as rough quantitative agreement with seemingly disparate results.

We present the Data Release 12 Quasar catalog (DR12Q) from the Baryon Oscillation Spectroscopic Survey (BOSS) of the SDSS-III. This catalog includes all SDSS-III/BOSS objects that were spectroscopically targeted as quasar candidates during the full survey and that are confirmed as quasars via visual inspection of the spectra, have luminosities Mi[z=2]<-20.5 (in a LCDM cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega _L=0.7), and either display at least one emission line with a full width at half maximum (FWHM)larger than 500 km/s or, if not, have interesting/complex absorption features. The catalog also includes previously known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 297,301 quasars detected over 9,376 square degrees with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with z>2.15 is about an order of magnitude greater than the number of z>2.15 quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 29,580 broad absorption line quasars and lists their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry together with some information on the optical morphology and the selection criteria. When available, the catalog also provides information on the optical variability of quasars using SDSS and PTF multi-epoch photometry. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra, covering the wavelength region 3,600-10,500A at a spectral resolution in the range 1,300<R<2,500, can be retrieved from the SDSS Catalog Archive Server.

As star-forming clouds collapse, the gas within them fragments to ever-smaller masses. Naively one might expect this process to continue down to the smallest mass that is able to radiate away its binding energy on a dynamical timescale, the opacity limit for fragmentation, at $\sim 0.01$ $M_\odot$. However, the observed peak of the initial mass function (IMF) lies a factor of $20-30$ higher in mass, suggesting that some other mechanism halts fragmentation before the opacity limit is reached. In this paper we analyse radiation-magnetohydrodynamic simulations of star cluster formation in typical Milky Way environments in order to determine what physical process limits fragmentation in them. We examine the regions in the vicinity of stars that form in the simulations to determine the amounts of mass that are prevented from fragmenting by thermal and magnetic pressure. We show that, on small scales, thermal pressure enhanced by stellar radiation heating is the dominant mechanism limiting the ability of the gas to further fragment. In the brown dwarf mass regime, $\sim 0.01$ $M_\odot$, the typical object that forms in the simulations is surrounded by gas whose mass is several times its own that is unable to escape or fragment, and instead is likely to accrete. This mechanism explains why $\sim 0.01$ $M_\odot$ objects are rare: unless an outside agent intervenes (e.g., a shock strips away the gas around them), they will grow by accreting the warmed gas around them. In contrast, by the time stars grow to masses of $\sim 0.2$ $M_\odot$, the mass of heated gas is only tens of percent of the central star mass, too small to alter its final mass by a large factor. This naturally explains why the IMF peak is at $\sim 0.2$ $M_\odot$.

The Baryon Oscillation Spectroscopic Survey (BOSS) has collected more than 150,000 $2.1 \leq z \leq 3.5$ quasar spectra since 2009. Using this unprecedented sample, we create a composite spectrum in the rest-frame of 102,150 quasar spectra from 800 \AA\ to 3300 \AA\ at a signal-to-noise ratio close to 1000 per pixel ($\Delta v$ of 69 km~s$^{-1}$). Included in this analysis is a correction to account for flux calibration residuals in the BOSS spectrophotometry. We determine the spectral index as a function of redshift of the full sample, warp the composite spectrum to match the median spectral index, and compare the resulting spectrum to SDSS photometry used in target selection. The quasar composite matches the color of the quasar population to within 0.02 magnitudes in $g-r$, 0.03 magnitudes in $r-i$, and 0.01 magnitudes in $i-z$ over the redshift range $2.2<z<2.6$. The composite spectrum deviates from the imaging photometry by 0.05 magnitudes around $z = 2.7$, likely due to differences in target selection as the quasar colors become similar to the stellar locus at this redshift. Finally, we characterize the line features in the high signal-to-noise composite and identify nine faint lines not found in the previous composite spectrum from SDSS.

Recent studies have found that obscured quasars cluster more strongly and are thus hosted by dark matter haloes of larger mass than their unobscured counterparts. These results pose a challenge for the simplest unification models, in which obscured objects are intrinsically the same as unobscured sources but seen through a dusty line of sight. There is general consensus that a structure like a "dusty torus" exists, meaning that this intrinsic similarity is likely the case for at least some subset of obscured quasars. However, the larger host halo masses of obscured quasars implies that there is a second obscured population that has an even higher clustering amplitude and typical halo mass. Here, we use simple assumptions about the host halo mass distributions of quasars, along with analytical methods and cosmological $N$-body simulations to isolate the signal from this population. We provide values for the bias and halo mass as a function of the fraction of the "non-torus obscured" population. Adopting a reasonable value for this fraction of $\sim$25% implies a non-torus obscured quasar bias that is much higher than the observed obscured quasar bias, because a large fraction of the obscured population shares the same clustering strength as the unobscured objects. For this non-torus obscured population, we derive a bias of $\sim$3, and typical halo masses of $\sim3\times10^{13}$ M$_{\odot}/h$ at $z=1$. These massive haloes are likely the descendants of high-mass unobscured quasars at high redshift, and will evolve into members of galaxy groups at $z=0$.

Upgrades at the 2.3 meter Wyoming Infrared Observatory telescope have provided the capability for fully-remote operations by a single operator from the University of Wyoming campus. A line-of-sight 300 Megabit/s 11 GHz radio link provides high-speed internet for data transfer and remote operations that include several real-time video feeds. Uninterruptable power is ensured by a 10 kVA battery supply for critical systems and a 55 kW autostart diesel generator capable of running the entire observatory for up to a week. Construction of a new four-element prime-focus corrector with fused-silica elements allows imaging over a 40' field-of-view with a new 4096x4096 UV-sensitive prime-focus camera and filter wheel. A new telescope control system facilitates the remote operations model and provides 20'' rms pointing over the usable sky. Taken together, these improvements pave the way for a new generation of sky surveys supporting space-based missions and flexible-cadence observations advancing emerging astrophysical priorities such as planet detection, quasar variability, and long-term time-domain campaigns.

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.

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 present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

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.

We report the discovery of a quasar pair at z=5 separated by 21 arcsec. Both objects were identified as quasar candidates using simple color selection techniques applied to photometric catalogs from the CFHT Legacy Survey (CFHTLS). Spectra obtained with the MMT present no discernible offset in redshift between the two objects; on the other hand, there are clear differences in the emission line profiles and in the multiwavelength spectral energy distributions that strongly disfavor the hypothesis that they are gravitationally lensed images of a single quasar. Both quasars are surprisingly bright given their proximity (a projected separation of ~135 kpc), with i=19.4 and i=21.4. Previous measurements of the luminosity function demonstrate that luminous quasars are extremely rare at z=5; the existence of this pair suggests that quasars have strong small-scale clustering at high redshift. Assuming a real-space correlation function of the form $\xi(r) \propto (r/r_0)^{-2}$, this discovery implies a correlation length $r_0 > 20 h^{-1}$ Mpc, consistent with a rapid strengthening of quasar clustering at high redshift as seen in previous observations and predicted by theoretical models where feedback effects are inefficient at shutting down black hole growth at high redshift.

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.

The SDSS-IV/eBOSS has an extensive quasar program that combines several selection methods. Among these, the photometric variability technique provides highly uniform samples, unaffected by the redshift bias of traditional optical-color selections, when $z= 2.7 - 3.5$ quasars cross the stellar locus or when host galaxy light affects quasar colors at $z < 0.9$. Here, we present the variability selection of quasars in eBOSS, focusing on a specific program that led to a sample of 13,876 quasars to $g_{\rm dered}=22.5$ over a 94.5 deg$^2$ region in Stripe 82, an areal density 1.5 times higher than over the rest of the eBOSS footprint. We use these variability-selected data to provide a new measurement of the quasar luminosity function (QLF) in the redshift range $0.68<z<4.0$. Our sample is denser, reaches deeper than those used in previous studies of the QLF, and is among the largest ones. At the faint end, our QLF extends to $M_g(z\!=\!2)=-21.80$ at low redshift and to $M_g(z\!=\!2)=-26.20$ at $z\sim 4$. We fit the QLF using two independent double-power-law models with ten free parameters each. The first model is a pure luminosity-function evolution (PLE) with bright-end and faint-end slopes allowed to be different on either side of $z=2.2$. The other is a simple PLE at $z<2.2$, combined with a model that comprises both luminosity and density evolution (LEDE) at $z>2.2$. Both models are constrained to be continuous at $z=2.2$. They present a flattening of the bright-end slope at large redshift. The LEDE model indicates a reduction of the break density with increasing redshift, but the evolution of the break magnitude depends on the parameterization. The models are in excellent accord, predicting quasar counts that agree within 0.3\% (resp., 1.1\%) to $g<22.5$ (resp., $g<23$). The models are also in good agreement over the entire redshift range with models from previous studies.

As part of the Sloan Digital Sky Survey IV the extended Baryon Oscillation Spectroscopic Survey (eBOSS) will improve measurements of the cosmological distance scale by applying the Baryon Acoustic Oscillation (BAO) method to quasar samples. eBOSS will adopt two approaches to target quasars over 7500 sq. deg. First, a "CORE" quasar sample will combine optical selection in ugriz using a likelihood-based routine called XDQSOz, with a mid-IR-optical color-cut. eBOSS CORE selection (to g < 22 OR r < 22) should return ~ 70 quasars per sq. deg. at redshifts 0.9 < z < 2.2 and ~7 z > 2.1 quasars per sq. deg. Second, a selection based on variability in multi-epoch imaging from the Palomar Transient Factory should recover an additional ~3-4 z > 2.1 quasars per sq. deg. to g < 22.5. A linear model of how imaging systematics affect target density recovers the angular distribution of eBOSS CORE quasars over 96.7% (76.7%) of the SDSS North (South) Galactic Cap area. The eBOSS CORE quasar sample should thus be sufficiently dense and homogeneous over 0.9 < z < 2.2 to yield the first few-percent-level BAO constraint near z~1.5. eBOSS quasars at z > 2.1 will be used to improve BAO measurements in the Lyman-alpha Forest. Beyond its key cosmological goals, eBOSS should be the next-generation quasar survey, comprising > 500,000 new quasars and > 500,000 uniformly selected spectroscopically confirmed 0.9 < z < 2.2 quasars. At the conclusion of eBOSS, the SDSS will have provided unique spectra of over 800,000 quasars.

Representing simultaneous black hole accretion during a merger, binary active galactic nuclei (AGNs) could provide valuable observational constraints to models of galaxy mergers and AGN triggering. High-resolution radio interferometer imaging offers a promising method to identify a large and uniform sample of binary AGNs, because it probes a generic feature of nuclear activity and is free from dust obscuration. Our previous search yielded 52 strong candidates of kpc-scale binaries over the 92 deg$^2$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 area with 2"-resolution Very Large Array (VLA) images. Here we present 0.3"-resolution VLA 6 GHz observations for six candidates that have complete optical spectroscopy. The new data confirm the binary nature of four candidates and identify the other two as line-of-sight projections of radio structures from single AGNs. The four binary AGNs at $z \sim 0.1$ reside in major mergers with projected separations of 4.2-12 kpc. Optical spectral modeling shows that their hosts have stellar masses between $10.3 < log(M_\star/M_\odot) < 11.5$ and velocity dispersions between $120 < \sigma_\star < 320$ km/s. The radio emission is compact ($<$0.4") and show steep spectrum ($-1.8 < \alpha < -0.5$) at 6 GHz. The host galaxy properties and the Eddington-scaled accretion rates broadly correlate with the excitation state, similar to the general radio-AGN population at low redshifts. Our estimated binary AGN fraction indicates that simultaneous accretion occurs $>23^{+15}_{-8}$% of the time when a kpc-scale galaxy pair is detectable as a radio-AGN. The high duty cycle of the binary phase strongly suggests that major mergers can trigger and synchronize black hole accretion.

Using the most recent releases of WISE and Planck data, we perform updated measurements of the bias and typical dark matter halo mass of infrared-selected obscured and unobscured quasars, using the angular autocorrelation function and cosmic microwave background (CMB) lensing cross-correlations. Since our recent work of this kind, the WISE Allwise catalogue was released with improved photometry, and the Planck mission was completed and released improved products. These new data provide a more reliable measurement of the quasar bias and provide an opportunity to explore the role of changing survey pipelines in results downstream. We present a comparison of IR color-selected quasars, split into obscured and unobscured populations based on optical-IR colors, selected from two versions of the WISE data. Which combination of data is used impacts the final results, particularly for obscured quasars, both because of mitigation of some systematics and because the newer catalogue provides a slightly different sample. We show that Allwise data is superior in several ways, though there may be some systematic trends with Moon contamination that were not present in the previous catalogue. We opt currently for the most conservative sample that meet our selection criteria in both the previous and new WISE catalogues. We measure a higher bias and halo mass for obscured quasars ($b_{\textrm{obsc}} \sim 2.1$, $b_{\textrm{unob}} \sim 1.8$) --- at odds with simple orientation models --- but at a reduced significance ($\sim$1.5$\sigma$) as compared to our work with previous survey data.

Particle methods are a ubiquitous tool for solving the Vlasov-Poisson equation in comoving coordinates, which is used to model the gravitational evolution of dark matter in an expanding universe. However, these methods are known to produce poor results on idealized test problems, particularly at late times, after the particle trajectories have crossed. To investigate this, we have performed a series of one- and two-dimensional "Zel'dovich Pancake" calculations using the popular Particle-in-Cell (PIC) method. We find that PIC can indeed converge on these problems provided the following modifications are made. The first modification is to regularize the singular initial distribution function by introducing a small but finite artificial velocity dispersion. This process is analogous to artificial viscosity in compressible gas dynamics, and, as with artificial viscosity, the amount of regularization can be tailored so that its effect outside of a well-defined region - in this case, the high-density caustics - is small. The second modification is the introduction of a particle remapping procedure that periodically re-expresses the dark matter distribution function using a new set of particles. We describe a remapping algorithm that is third-order accurate and adaptive in phase space. This procedure prevents the accumulation of numerical errors in integrating the particle trajectories from growing large enough to significantly degrade the solution. Once both of these changes are made, PIC converges at second order on the Zel'dovich Pancake problem, even at late times, after many caustics have formed. Furthermore, the resulting scheme does not suffer from the unphysical, small-scale "clumping" phenomenon known to occur on the Pancake problem when the perturbation wave vector is not aligned with one of the Cartesian coordinate axes.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

The intermittent dissipation of interstellar turbulence is an important energy source in the diffuse ISM. Though on average smaller than the heating rates due to cosmic rays and the photoelectric effect on dust grains, the turbulent cascade can channel large amounts of energy into a relatively small fraction of the gas that consequently undergoes significant heating and chemical enrichment. In particular, this mechanism has been proposed as a solution to the long-standing problem of the high abundance of CH+ along diffuse molecular sight lines, which steady-state, low temperature models under-produce by over an order of magnitude. While much work has been done on the structure and chemistry of these small-scale dissipation zones, comparatively little attention has been paid to relating these zones to the properties of the large-scale turbulence. In this paper, we attempt to bridge this gap by estimating the temperature and CH+ column density along diffuse molecular sight-lines by post-processing 3-dimensional MHD turbulence simulations. Assuming reasonable values for the cloud density (30 / cm^3), size (20 pc), and velocity dispersion (2.3 km / s), we find that our computed abundances compare well with CH+ column density observations, as well as with observations of emission lines from rotationally excited H2 molecules.

We describe the algorithm used to select the Luminous Red Galaxy (LRG) sample for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV) using photometric data from both the SDSS and the Wide-Field Infrared Survey Explorer (WISE). LRG targets are required to meet a set of color selection criteria and have z-band and i-band MODEL magnitudes z < 19.95 and 19.9 < i < 21.8, respectively. Our algorithm selects roughly 50 LRG targets per square degree, the great majority of which lie in the redshift range 0.6 < z < 1.0 (median redshift 0.71). We demonstrate that our methods are highly effective at eliminating stellar contamination and lower-redshift galaxies. We perform a number of tests using spectroscopic data from SDSS-III/BOSS to determine the redshift reliability of our target selection and its ability to meet the science requirements of eBOSS. The SDSS spectra are of high enough signal-to-noise ratio that at least 89% of the target sample yields secure redshift measurements. We also present tests of the uniformity and homogeneity of the sample, demonstrating that it should be clean enough for studies of the large-scale structure of the universe at higher redshifts than SDSS-III/BOSS LRGs reached.

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

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 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.

Extreme deconvolution (XD) of broad-band photometric data can both separate stars from quasars and generate probability density functions for quasar redshifts, while incorporating flux uncertainties and missing data. Mid-infrared photometric colors are now widely used to identify hot dust intrinsic to quasars, and the release of all-sky WISE data has led to a dramatic increase in the number of IR-selected quasars. Using forced-photometry on public WISE data at the locations of SDSS point sources, we incorporate this all-sky data into the training of the XDQSOz models originally developed to select quasars from optical photometry. The combination of WISE and SDSS information is far more powerful than SDSS alone, particularly at $z>2$. The use of SDSS$+$WISE photometry is comparable to the use of SDSS$+$ultraviolet$+$near-IR data. We release a new public catalogue of 5,537,436 (total; 3,874,639 weighted by probability) potential quasars with probability $P_{\textrm{QSO}} > 0.2$. The catalogue includes redshift probabilities for all objects. We also release an updated version of the publicly available set of codes to calculate quasar and redshift probabilities for various combinations of data. Finally, we demonstrate that this method of selecting quasars using WISE data is both more complete and efficient than simple WISE color-cuts, especially at high redshift. Our fits verify that above $z \sim 3$ WISE colors become bluer than the standard cuts applied to select quasars. Currently, the analysis is limited to quasars with optical counterparts, and thus cannot be used to find highly obscured quasars that WISE color-cuts identify in significant numbers.

Many astronomers now participate in large international collaborations, and it is important to examine whether these structures foster a scientific climate that is inclusive and diverse. The Committee on the Participation of Women in the Sloan Digital Sky Survey (CPWS) was formed to evaluate the demographics and gender climate within SDSS-IV, one of the largest and most geographically distributed astronomical collaborations. In April 2014, the CPWS administered a demographic survey to establish a baseline for the incipient SDSS-IV. We received responses from 250 participants (46% of the active membership). Half of the survey respondents were located in the US or Canada and 30% were based in Europe. Eleven percent of survey respondents considered themselves to be an ethnic minority at their current institution. Twenty-five percent of the SDSS-IV collaboration members are women, a fraction that is consistent with the US astronomical community, but substantially higher than the fraction of women in the IAU (16%). Approximately equal fractions of men and women report holding positions of leadership. When binned by academic age and career level, men and women also assume leadership roles at approximately equal rates, in a way that increases steadily for both genders with increasing seniority. In this sense, SDSS-IV has been successful in recruiting leaders that are representative of the collaboration. Yet, more progress needs to be made towards achieving gender balance and increasing diversity in the field of astronomy, and there is still room for improvement in the membership and leadership of SDSS-IV. For example, at the highest level of SDSS-IV leadership, women disproportionately assume roles related to education and public outreach. The CPWS plans to use these initial data to establish a baseline for tracking demographics over time as we work to assess and improve the climate of SDSS-IV.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

Using data from the All Wavelength Extended Groth Strip International Survey (AEGIS) we statistically detect the extended X-ray emission in the interstellar medium (ISM)/intra-cluster medium (ICM) in both active and normal galaxies at 0.3 <= z <= 1.3. For both active galactic nuclei (AGN) host galaxy and normal galaxy samples that are matched in restframe color, luminosity, and redshift distribution, we tentatively detect excess X-ray emission at scales of 1--10 arcsec at a few sigma significance in the surface brightness profiles. The exact significance of this detection is sensitive to the true characterization of Chandra's point spread function. The observed excess in the surface brightness profiles is suggestive of lower extended emission in AGN hosts compared to normal galaxies. This is qualitatively similar to theoretical predictions of the X-ray surface brightness profile from AGN feedback models, where feedback from AGN is likely to evacuate the gas from the center of the galaxy/cluster. We propose that AGN that are intrinsically under-luminous in X-rays, but have equivalent bolometric luminosities to our sources will be the ideal sample to study more robustly the effect of AGN feedback on diffuse ISM/ICM gas.

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.

Using data from the DEEP2 galaxy redshift survey and the All Wavelength Extended Groth Strip International Survey we obtain stacked X-ray maps of galaxies at 0.7 < z < 1.0 as a function of stellar mass. We compute the total X-ray counts of these galaxies and show that in the soft band (0.5--2,kev) there exists a significant correlation between galaxy X-ray counts and stellar mass at these redshifts. The best-fit relation between X-ray counts and stellar mass can be characterized by a power law with a slope of 0.58 +/- 0.1. We do not find any correlation between stellar mass and X-ray luminosities in the hard (2--7,kev) and ultra-hard (4--7,kev) bands. The derived hardness ratios of our galaxies suggest that the X-ray emission is degenerate between two spectral models, namely point-like power-law emission and extended plasma emission in the interstellar medium. This is similar to what has been observed in low redshift galaxies. Using a simple spectral model where half of the emission comes from power-law sources and the other half from the extended hot halo we derive the X-ray luminosities of our galaxies. The soft X-ray luminosities of our galaxies lie in the range 10^39-8x10^40, ergs/s. Dividing our galaxy sample by the criteria U-B > 1, we find no evidence that our results for X-ray scaling relations depend on optical color.

(Abridged) We detect the large-scale structure of Lya emission in the Universe at redshifts z=2-3.5 by measuring the cross-correlation of Lya surface brightness with quasars in SDSS/BOSS. We use a million spectra targeting Luminous Red Galaxies at z<0.8, after subtracting a best fit model galaxy spectrum from each one, as an estimate of the high-redshift Lya surface brightness. The quasar-Lya emission cross-correlation we detect has a shape consistent with a LambdaCDM model with Omega_M =0.30^+0.10-0.07. The predicted amplitude of this cross-correlation is proportional to the product of the mean Lya surface brightness, <mu_alpha>, the amplitude of mass fluctuations, and the quasar and Lya emission bias factors. Using known values, we infer <mu_alpha>(b_alpha/3) = (3.9 +/- 0.9) x 10^-21 erg/s cm^-2 A^-1 arcsec^-2, where b_alpha is the Lya emission bias factor. If the dominant sources of Lya emission are star forming galaxies, we infer rho_SFR = (0.28 +/- 0.07) (3/b_alpha) /yr/Mpc^3 at z=2-3.5. For b_alpha=3, this value is a factor of 21-35 above previous estimates from individually detected Lya emitters, although consistent with the total rho_SFR derived from dust-corrected, continuum UV surveys. 97% of the Lya emission in the Universe at these redshifts is therefore undetected in previous surveys of Lya emitters. Our measurement is much greater than seen from stacking analyses of faint halos surrounding previously detected Lya emitters, but we speculate that it arises from similar Lya halos surrounding all luminous star-forming galaxies. We also detect redshift space anisotropy of the quasar-Lya emission cross-correlation, finding evidence at the 3.0 sigma level that it is radially elongated, consistent with distortions caused by radiative-transfer effects (Zheng et al. (2011)). Our measurements represent the first application of the intensity mapping technique to optical observations.

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky, and delivered over a trillion pixels of imaging data. We present the large-scale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that have been classified from this imaging, representing the highest density of quasars ever studied for clustering measurements. This data set spans ~11,000 square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large volume and medium density of tracers should facilitate high-precision cosmological constraints. We measure the angular clustering of photometrically classified quasars using an optimal quadratic estimator in four redshift slices with an accuracy of ~25% over a bin width of l ~10 - 15 on scales corresponding to matter-radiation equality and larger (l ~ 2 - 30). Observational systematics can strongly bias clustering measurements on large scales, which can mimic cosmologically relevant signals such as deviations from Gaussianity in the spectrum of primordial perturbations. We account for systematics by employing a new method recently proposed by Agarwal et al. (2014) to the clustering of photometrically classified quasars. We carefully apply our methodology to mitigate known observational systematics and further remove angular bins that are contaminated by unknown systematics. Combining quasar data with the photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et al. (2012), and marginalizing over all bias and shot noise-like parameters, we obtain a constraint on local primordial non-Gaussianity of fNL = -113+/-154 (1\sigma error). [Abridged]

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

We report a detection of the baryon acoustic oscillation (BAO) feature in the flux-correlation function of the Ly{\alpha} forest of high-redshift quasars with a statistical significance of five standard deviations. The study uses 137,562 quasars in the redshift range $2.1\le z \le 3.5$ from the Data Release 11 (DR11) of the Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-III. This sample contains three times the number of quasars used in previous studies. The measured position of the BAO peak determines the angular distance, $D_A(z=2.34)$ and expansion rate, $H(z=2.34)$, both on a scale set by the sound horizon at the drag epoch, $r_d$. We find $D_A/r_d=11.28\pm0.65(1\sigma)^{+2.8}_{-1.2}(2\sigma)$ and $D_H/r_d=9.18\pm0.28(1\sigma)\pm0.6(2\sigma)$ where $D_H=c/H$. The optimal combination, $\sim D_H^{0.7}D_A^{0.3}/r_d$ is determined with a precision of $\sim2\%$. For the value $r_d=147.4~{\rm Mpc}$, consistent with the CMB power spectrum measured by Planck, we find $D_A(z=2.34)=1662\pm96(1\sigma)~{\rm Mpc}$ and $H(z=2.34)=222\pm7(1\sigma)~{\rm km\,s^{-1}Mpc^{-1}}$. Tests with mock catalogs and variations of our analysis procedure have revealed no systematic uncertainties comparable to our statistical errors. Our results agree with the previously reported BAO measurement at the same redshift using the quasar-Ly{\alpha} forest cross-correlation. The auto-correlation and cross-correlation approaches are complementary because of the quite different impact of redshift-space distortion on the two measurements. The combined constraints from the two correlation functions imply values of $D_A/r_d$ and $D_H/r_d$ that are, respectively, 7% low and 7% high compared to the predictions of a flat $\Lambda$CDM cosmological model with the best-fit Planck parameters. With our estimated statistical errors, the significance of this discrepancy is $\approx 2.5\sigma$.

We use spectroscopy of close pairs of quasars to study diffuse gas in the circumgalactic medium (CGM) surrounding a sample of 40 Damped Lya systems (DLAs). The primary sightline in each quasar pair probes an intervening DLA in the redshift range 1.6 < z_DLA < 3.6, such that the second quasar sightline then probes Lya, CII, SiII, and CIV absorption in the CGM transverse to the DLA to projected distances $R_{\perp} < 300$ kpc. Analysis of the Lya profiles in these CGM sightlines constrains the covering fraction (f_C) of optically thick HI (having column density N_HI > 10^17.2 cm^-2) to be greater than ~30% within $R_{\perp} < 200$ kpc of DLAs. Strong SiII 1526 absorption with equivalent width W_1526 > 0.2 Ang occurs with an incidence f_C (W_1526 > 0.2 Ang) = 20(+12/-8)% within $R_{\perp}<100$ kpc, indicating that low-ionization metal absorption associated with DLAs probes material at a physical distance R_3D < 30 kpc. However, we find that strong CIV 1548 absorption is ubiquitous in these environments (f_C (W_1548 > 0.2 Ang) = 57(+12/-13)% within $R_{\perp} < 100$ kpc), and in addition exhibits a high degree of kinematic coherence on scales up to ~175 kpc. We infer that this high-ionization material arises predominantly in large, quiescent structures extending beyond the scale of the DLA host dark matter halos rather than in ongoing galactic winds. The Lya equivalent width in the DLA-CGM is anticorrelated with $R_{\perp}$ at >98% confidence, suggesting that DLAs arise close to the centers of their host halos rather than on their outskirts. Finally, the average Lya, CII and CIV equivalent widths are consistent with those measured around z~2 Lyman Break Galaxies. Assuming that DLAs trace a galaxy population with lower masses and luminosities, this finding implies that the absorption strength of cool circumgalactic material has a weak dependence on dark matter halo mass for M_h < 10^12 M_sun.

Galaxy mergers play an important role in the growth of galaxies and their supermassive black holes. Simulations suggest that tidal interactions could enhance black hole accretion, which can be tested by the fraction of binary active galactic nuclei (AGNs) among galaxy mergers. But determining the fraction requires a statistical sample of binaries. We have identified kpc-scale binary AGNs directly from high-resolution radio imaging. Inside the 92 square deg covered by the high-resolution Very Large Array survey of the Sloan Digital Sky Survey (SDSS) Stripe 82 field, we identified 22 grade A and 30 grade B candidates of binary radio AGNs with angular separations less than 5" (10 kpc at z = 0.1). Eight of the candidates have optical spectra for both components from the SDSS spectroscopic surveys and our Keck program. Two grade B candidates are projected pairs, but the remaining six candidates are all compelling cases of binary AGNs based on either emission line ratios or the excess in radio power compared to the H-alpha-traced star formation rate. Only two of the six binaries were previously discovered by an optical spectroscopic search. Based on these results, we estimate that ~60% of our binary candidates would be confirmed once we obtain complete spectroscopic information. We conclude that wide-area high-resolution radio surveys offer an efficient method to identify large samples of binary AGNs. These radio-selected binary AGNs complement binaries identified at other wavelengths and are useful for understanding the triggering mechanisms of black hole accretion.

We cross-correlate a cosmic microwave background (CMB) lensing map with the projected space densities of quasars to measure the bias and halo masses of a quasar sample split into obscured and unobscured populations, the first application of this method to distinct quasar subclasses. Several recent studies of the angular clustering of obscured quasars have shown that these objects likely reside in higher-mass halos compared to their unobscured counterparts. This has important implications for models of the structure and geometry of quasars, their role in growing supermassive black holes, and mutual quasar/host galaxy evolution. However, the magnitude and significance of this difference has varied from study to study. Using data from \planck, \wise, and SDSS, we follow up on these results using the independent method of CMB lensing cross-correlations. The region and sample are identical to that used for recent angular clustering measurements, allowing for a direct comparison of the CMB-lensing and angular clustering methods. At $z \sim 1$, we find that the bias of obscured quasars is $b_q = 2.57 \pm 0.24$, while that of unobscured quasars is $b_q = 1.89 \pm 0.19$. This corresponds to halo masses of $\log (M_h / M_{\odot} h^{-1}) = 13.24_{-0.15}^{+0.14}$ (obscured) and $\log (M_h / M_{\odot} h^{-1}) = 12.71_{-0.13}^{+0.15}$ (unobscured). These results agree well with with those from angular clustering (well within $1\sigma$), and confirm that obscured quasars reside in host halos $\sim$3 times as massive as halos hosting unobscured quasars. This implies that quasars spend a significant portion of their lifetime in an obscured state, possibly more than one half of the entire active phase.

Blazars are classically divided into the BL Lac (BLL) and Flat-Spectrum Radio Quasar (FSRQ) subclasses, corresponding to radiatively inefficient and efficient accretion regimes, respectively, largely based on the equivalent width (EW) of their optical broad emission lines (BEL). However, EW-based classification criteria are not physically motivated, and a few blazars have previously 'transitioned' from one subclass to the other. We present the first systematic search for these transition blazars in a sample of 602 unique pairs of repeat spectra of 354 blazars in SDSS, finding six clear cases. These transition blazars have bolometric Eddington ratios of ~0.3 and low-frequency synchrotron peaks, and are thus FSRQ-like. We show that the strong EW variability (up to an unprecedented factor of >60) is due to swamping of the BELs from variability in jet continuum emission, which is stronger in amplitude and shorter in timescale than typical blazars. Although these transition blazars appear to switch between FSRQ and BLL according to the phenomenologically-based EW scheme, we show that they are most likely rare cases of FSRQs with radiatively efficient accretion flows and especially strongly-beamed jets. These results have implications for the decrease of the apparent BLL population at high-redshifts, and may add credence to claims of a negative BLL redshift evolution.

We present the results of an optical spectroscopic survey of a sample of 40 candidate obscured quasars identified on the basis of their mid-infrared emission detected by the Wide-Field Infrared Survey Explorer (WISE). Optical spectra for this survey were obtained using the Robert Stobie Spectrograph (RSS) on the Southern African Large Telescope (SALT). Our sample was selected with WISE colors characteristic of AGNs, as well as red optical to mid-IR colors indicating that the optical/UV AGN continuum is obscured by dust. We obtain secure redshifts for the majority of the objects that comprise our sample (35/40), and find that sources that are bright in the WISE W4 (22$\mu$m) band are typically at moderate redshift (<z> = 0.35$) while sources fainter in W4 are at higher redshifts (<z> = 0.73$). The majority of the sources have narrow emission lines, with optical colors and emission line ratios of our WISE-selected sources that are consistent with the locus of AGN on the rest-frame $g-z$ color vs. [NeIII]$\lambda$3869 / [OII]$\lambda\lambda$3726+3729 line ratio diagnostic diagram. We also use empirical AGN and galaxy templates to model the spectral energy distributions (SEDs) for the objects in our sample, and find that while there is significant variation in the observed SEDs for these objects, the majority require a strong AGN component. Finally, we use the results from our analysis of the optical spectra and the SEDs to compare our selection criteria to alternate criteria presented in the literature. These results verify the efficacy of selecting luminous obscured AGNs based on their WISE colors.

Studies of radio-loud (RL) broad absorption line (BAL) quasars indicate that popular orientation-based BAL models fail to account for all observations. Are these results extendable to radio-quiet (RQ) BAL quasars? Comparisons of RL and RQ BAL quasars show that many of their properties are quite similar. Here we extend these analyses to the rest-frame ultraviolet (UV) spectral properties, using a sample of 73 RL and 473 RQ BAL quasars selected from the Sloan Digital Sky Survey (SDSS). Each RQ quasar is individually matched to a RL quasar in both redshift (over the range $1.5 < z < 3.5$) and continuum luminosity. We compare several continuum, emission line, and absorption line properties, as well as physical properties derived from these measurements. Most properties in the samples are statistically identical, though we find slight differences in the velocity structure of the BALs that cause apparent differences in CIV emission line properties. Differences in the velocities may indicate an interaction between the radio jets and the absorbing material. We also find that UV FeII emission is marginally stronger in RL BAL quasars. All of these differences are subtle, so in general we conclude that RL and RQ BAL QSOs are not fundamentally different objects, except in their radio properties. They are therefore likely to be driven by similar physical phenomena, suggesting that results from samples of RL BAL quasars can be extended to their RQ counterparts.

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 the first cross-correlation measurement between Sloan Digital Sky Survey (SDSS) Type 1 quasars and the cosmic infrared background (CIB) measured by Herschel. The distribution of the quasars at 0.15<z<3.5 covers the redshift range where we expect most of the CIB to originate. We detect the sub-mm emission of the quasars, which dominates on small scales, as well as correlated emission from dusty star-forming galaxies (DSFGs) dominant on larger scales. The mean sub-mm flux densities of the DR7 quasars (median redshift <z>=1.4) is $11.1^{+1.6}_{-1.4}$, $7.1^{+1.6}_{-1.3}$ and $3.6^{+1.4}_{-1.0}$ mJy at 250, 350 and 500 microns, respectively, while the mean sub-mm flux densities of the DR9 quasars (<z>=2.5) is $5.7^{+0.7}_{-0.6}$, $5.0^{+0.8}_{-0.7}$ and $1.8^{+0.5}_{-0.4}$ mJy. We find that the correlated sub-mm emission includes both the emission from satellite DSFGs in the same halo as the central quasar and the emission from DSFGs in separate halos (correlated with the quasar-hosting halo). The amplitude of the one-halo term is ~10 times smaller than the sub-mm emission of the quasars, implying the the satellites have a lower star-formation rate than the quasars. The satellite fraction for the DR7 quasars is $0.008^{+0.008}_{-0.005}$ and the host halo mass scale for the central and satellite quasars is $10^{12.36\pm0.87}$ M$_{\odot}$ and $10^{13.60\pm0.38}$ M$_{\odot}$, respectively. The satellite fraction of the DR9 quasars is $0.065^{+0.021}_{-0.031}$ and the host halo mass scale for the central and satellite quasars is $10^{12.29\pm0.62}$ M$_{\odot}$ and $10^{12.82\pm0.39}$ M$_{\odot}$, respectively. Thus, the typical halo environment of the SDSS Type 1 quasars is found to be similar to that of DSFGs, which supports the generally accepted view that dusty starburst and quasar activity are evolutionarily linked.

Recent studies of luminous infrared-selected active galactic nuclei (AGN) suggest that the reddest, most obscured objects display a higher angular clustering amplitude, and thus reside in higher-mass dark matter halos. This is a direct contradiction to the prediction of the simplest unification-by-orientation models of AGN and quasars. However, clustering measurements depend strongly on the "mask" that removes low-quality data and describes the sky and selection function. We find that applying a robust, conservative mask to WISE-selected quasars yields a weaker but still significant difference in the bias between obscured and unobscured quasars. These findings are consistent with results from previous Spitzer surveys, and removes any scale dependence of the bias. For obscured quasars with $\langle z \rangle = 0.99$ we measure a bias of $b_q = 2.67 \pm 0.16$, corresponding to a halo mass of $\log (M_h / M_{\odot} h^{-1}) = 13.3 \pm 0.1$, while for unobscured sources with $\langle z \rangle = 1.04$ we find $b_q = 2.04 \pm 0.17$ with a halo mass $\log (M_h / M_{\odot} h^{-1} )= 12.8 \pm 0.1$. This improved measurement indicates that WISE-selected obscured quasars reside in halos only a few times more massive than the halos of their unobscured counterparts, a reduction in the factor of $\sim$10 larger halo mass as has been previously reported using WISE-selected samples. Additionally, an abundance matching analysis yields lifetimes for both obscured and unobscured quasar phases on the order of a few 100 Myr ($\sim$ 1\% of the Hubble time) --- however, the obscured phase lasts roughly twice as long, in tension with many model predictions.

We present and release co-added images of the Sloan Digital Sky Survey (SDSS) Stripe 82. Stripe 82 covers an area of 300 deg^2 on the Celestial Equator, and has been repeatedly scanned 70-90 times in the ugriz bands by the SDSS imaging survey. By making use of all available data in the SDSS archive, our co-added images are optimized for depth. Input single-epoch frames were properly processed and weighted based on seeing, sky transparency, and background noise before co-addition. The resultant products are co-added science images and their associated weight images that record relative weights at individual pixels. The depths of the co-adds, measured as the 5 sigma detection limits of the aperture (3.2 arcsec diameter) magnitudes for point sources, are roughly 23.9, 25.1, 24.6, 24.1, and 22.8 AB magnitudes in the five bands, respectively. They are 1.9-2.2 mag deeper than the best SDSS single-epoch data. The co-added images have good image quality, with an average point-spread function FWHM of ~1 arcsec in the r, i, and z bands. We also release object catalogs that were made with SExtractor. These co-added products have many potential uses for studies of galaxies, quasars, and Galactic structure. We further present and release near-IR J-band images that cover ~90 deg^2 of Stripe 82. These images were obtained using the NEWFIRM camera on the NOAO 4-m Mayall telescope, and have a depth of about 20.0--20.5 Vega magnitudes (also 5 sigma detection limits for point sources).

We measure the large-scale cross-correlation of quasars with the Lyman alpha forest absorption, using over 164,000 quasars from Data Release 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We extend the previous study of roughly 60,000 quasars from Data Release 9 to larger separations, allowing a measurement of the Baryonic Acoustic Oscillation (BAO) scale along the line of sight $c/(H(z=2.36) ~ r_s) = 9.0 \pm 0.3$ and across the line of sight $D_A(z=2.36) / ~ r_s = 10.8 \pm 0.4$, consistent with CMB and other BAO data. Using the best fit value of the sound horizon from Planck data ($r_s=147.49 Mpc$), we can translate these results to a measurement of the Hubble parameter of $H(z=2.36) = 226 \pm 8 km/s / Mpc$ and of the angular diameter distance of $D_A(z=2.36) = 1590 \pm 60 Mpc$. The measured cross-correlation function and an update of the code to fit the BAO scale (baofit) are made publicly available.

We present JHK near-infrared (NIR) spectroscopy of 25 candidate Type II quasars selected from the Sloan Digital Sky Survey, using Triplespec on the Apache Point Observatory 3.5m telescope, FIRE at the Magellan/Baade 6.5m telescope, and GNIRS on Gemini. At redshifts of 2 < z < 3.4, our NIR spectra probe the rest-frame optical region of these targets, which were initially selected to have strong lines of CIV and Ly alpha, with FWHM<2000 km/s from the SDSS pipeline. We use the [OIII]5007 line shape as a model for the narrow line region emission, and find that \halpha\ consistently requires a broad component with FWHMs ranging from 1000 to 7500 km/s. Interestingly, the CIV lines also require broad bases, but with considerably narrower widths of 1000 to 4500 km/s. Estimating the extinction using the Balmer decrement and also the relationship in lower-z quasars between rest equivalent width and luminosity in the [OIII] line, we find typical A_V values of 0-2 mag, which naturally explain the attenuated CIV lines relative to Halpha. We propose that our targets are moderately obscured quasars. We also describe one unusual object with three distinct velocity peaks in its [OIII] spectrum.

Quasar luminosity functions are a fundamental probe of the growth and evolution of supermassive black holes. Measuring the intrinsic luminosity function is difficult in practice, due to a multitude of observational and systematic effects. As sample sizes increase and measurement errors drop, characterizing the systematic effects is becoming more important. It is well known that the continuum emission from the accretion disk of quasars is anisotropic --- in part due to its disk-like structure --- but current luminosity function calculations effectively assume isotropy over the range of unobscured lines of sight. Here, we provide the first steps in characterizing the effect of random quasar orientations and simple models of anisotropy on observed luminosity functions. We find that the effect of orientation is not insignificant and exceeds other potential corrections such as those from gravitational lensing of foreground structures. We argue that current observational constraints may overestimate the intrinsic luminosity function by as much as a factor of ~2 on the bright end. This has implications for models of quasars and their role in the Universe, such as quasars' contribution to cosmological backgrounds.

We examine the spatial extent of the narrow-line regions (NLRs) of a sample of 30 luminous obscured quasars at $0.4 < z < 0.7$ observed with spatially resolved Gemini-N GMOS long-slit spectroscopy. Using the [OIII]$\lambda5007$ emission feature, we estimate the size of the NLR using a cosmology-independent measurement: the radius where the surface brightness falls to 10$^{-15}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. We then explore the effects of atmospheric seeing on NLR size measurements and conclude that direct measurements of the NLR size from observed profiles are too large by 0.1 - 0.2 dex on average, as compared to measurements made to best-fit S\'{e}rsic or Voigt profiles convolved with the seeing. These data, which span a full order of magnitude in IR luminosity ($\log{(L_{8 \mu \mathrm{m}} / \mathrm{erg\, s}^{-1})} = 44.4 - 45.4$) also provide strong evidence that there is a flattening of the relationship between NLR size and AGN luminosity at a seeing-corrected size of $\sim 7$ kpc. The objects in this sample have high luminosities which place them in a previously under-explored portion of the size-luminosity relationship. These results support the existence of a maximal size of the narrow-line region around luminous quasars; beyond this size either there is not enough gas, or the gas is over-ionized and does not produce enough [OIII]$\lambda5007$ emission.

Photometric large scale structure (LSS) surveys probe the largest volumes in the Universe, but are inevitably limited by systematic uncertainties. Imperfect photometric calibration leads to biases in our measurements of the density fields of LSS tracers such as galaxies and quasars, and as a result in cosmological parameter estimation. Earlier studies have proposed using cross-correlations between different redshift slices or cross-correlations between different surveys to reduce the effects of such systematics. In this paper we develop a method to characterize unknown systematics. We demonstrate that while we do not have sufficient information to correct for unknown systematics in the data, we can obtain an estimate of their magnitude. We define a parameter to estimate contamination from unknown systematics using cross-correlations between different redshift slices and propose discarding bins in the angular power spectrum that lie outside a certain contamination tolerance level. We show that this method improves estimates of the bias using simulated data and further apply it to photometric luminous red galaxies in the Sloan Digital Sky Survey as a case study.

We present three Orion simulations of star cluster formation in a 1000 Msun, turbulent molecular cloud clump, including the effects of radiative transfer, protostellar outflows, and magnetic fields. Our simulations all use self-consistent turbulent initial conditions and vary the mean mass-to-flux ratio relative to the critical value over 2, 10, and infinity to gauge the influence of magnetic fields on star cluster formation. We find, in good agreement with previous studies, that magnetic fields of typically observed strengths lower the star formation rate by a factor of 2.4 and reduce the amount of fragmentation by a factor of 2 relative to the zero-field case. We also find that the field increases the characteristic sink particle mass, again by a factor of 2.4. The magnetic field also increases the degree of clustering in our simulations, such that the maximum stellar densities in the strong field case are higher than the others by again a factor of 2. This clustering tends to encourage the formation of multiple systems, which are more common in the rad-MHD runs than the rad-hydro run. The companion frequency in our simulations is consistent with observations of multiplicity in Class I sources, particularly for the strong field case. Finally, we find evidence of primordial mass segregation in our simulations reminiscent of that observed in star clusters like the Orion Nebula Cluster.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present the Data Release 10 Quasar (DR10Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the first 2.5 years of the survey and that are confirmed as quasars via visual inspection of the spectra. The catalog also includes known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 166,583 quasars (74,454 are new discoveries since SDSS-DR9) detected over 6,373 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with $z>2.15$ (117,668) is $\sim$5 times greater than the number of $z>2.15$ quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII, MgII). The catalog identifies 16,461 broad absorption line quasars and gives their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry with typical accuracy of 0.03 mag and information on the optical morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3,600-10,500\AA\ at a spectral resolution in the range 1,300$<$R$<$2,500; the spectra can be retrieved from the SDSS Catalog Archive Server. We also provide a supplemental list of an additional 2,376 quasars that have been identified among the galaxy targets of the SDSS-III/BOSS.

We have developed two independent methods to measure the one-dimensional power spectrum of the transmitted flux in the Lyman-$\alpha$ forest. The first method is based on a Fourier transform, and the second on a maximum likelihood estimator. The two methods are independent and have different systematic uncertainties. The determination of the noise level in the data spectra was subject to a novel treatment, because of its significant impact on the derived power spectrum. We applied the two methods to 13,821 quasar spectra from SDSS-III/BOSS DR9 selected from a larger sample of over 60,000 spectra on the basis of their high quality, large signal-to-noise ratio, and good spectral resolution. The power spectra measured using either approach are in good agreement over all twelve redshift bins from $<z> = 2.2$ to $<z> = 4.4$, and scales from 0.001 $\rm(km/s)^{-1}$ to $0.02 \rm(km/s)^{-1}$. We determine the methodological and instrumental systematic uncertainties of our measurements. We provide a preliminary cosmological interpretation of our measurements using available hydrodynamical simulations. The improvement in precision over previously published results from SDSS is a factor 2--3 for constraints on relevant cosmological parameters. For a $\Lambda$CDM model and using a constraint on $H_0$ that encompasses measurements based on the local distance ladder and on CMB anisotropies, we infer $\sigma_8 =0.83\pm0.03$ and $n_s= 0.97\pm0.02$ based on \ion{H}{i} absorption in the range $2.1<z<3.7$.

Recent work on quasar clustering suggests a degeneracy in the halo occupation distribution constrained from two-point correlation functions. To break this degeneracy, we make the first empirical measurement of the mean occupation function (MOF) of quasars at $z \sim 0.2$ by matching quasar positions with groups and clusters identified in the MaxBCG sample. We fit two models to the MOF, a power law and a 4-parameter model. The number distribution of quasars in host halos is close to Poisson, and the slopes of the MOF obtained from our best-fit models (for the power law case) favor a MOF that monotonically increases with halo mass. The best-fit slopes are $0.53 \pm 0.04$ and $1.03\pm 1.12$ for the power law model and the 4-parameter model, respectively. We measure the radial distribution of quasars within dark matter halos and find it to be adequately described by a power law with a slope $-2.3 \pm 0.4$. We measure the conditional luminosity function (CLF) of quasars and show that there is no evidence that quasar luminosity depends on host halo mass, similar to the inferences drawn from clustering measurements. We also measure the conditional black hole mass function (CMF) of our quasars. Although the results are consistent with no dependence on halo mass, we observe a slight indication of downsizing of the black hole mass function. The lack of halo mass dependence in the CLF and CMF shows that quasars residing in galaxy clusters have characteristic luminosity and black hole mass scales.

We measure the cross-power spectrum of the projected mass density as traced by the convergence of the cosmic microwave background lensing field from the South Pole Telescope (SPT) and a sample of Type 1 and 2 (unobscured and obscured) quasars at z~1 selected with the Wide-field Infrared Survey Explorer, over 2500 deg^2. The cross-power spectrum is detected at ~7-sigma, and we measure a linear bias b=1.67+/-0.24, consistent with clustering analyses. Using an independent lensing map, derived from Planck observations, to measure the cross-spectrum, we find excellent agreement with the SPT analysis. The bias of the combined sample of Type 1 and 2 quasars determined in this work is similar to that previously determined for Type 1 quasars alone; we conclude that that obscured and unobscured quasars must trace the matter field in a similar way. This result has implications for our understanding of quasar unification and evolution schemes.

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 measure the large-scale cross-correlation of quasars with the Lyman \alpha\ forest absorption in redshift space, using ~ 60000 quasar spectra from Data Release 9 (DR9) of the Baryon Oscillation Spectroscopic Survey (BOSS). The cross-correlation is detected over a wide range of scales, up to comoving separations r of 80 Mpc/h. For r > 15 Mpc/h, we show that the cross-correlation is well fitted by the linear theory prediction for the mean overdensity around a quasar host halo in the standard \Lambda CDM model, with the redshift distortions indicative of gravitational evolution detected at high confidence. Using previous determinations of the Lyman \alpha\ forest bias factor obtained from the Lyman \alpha\ autocorrelation, we infer the quasar bias factor to be b_q = 3.64^+0.13_-0.15 at a mean redshift z=2.38, in agreement with previous measurements from the quasar auto-correlation. We also obtain a new estimate of the Lyman \alpha\ forest redshift distortion factor, \beta_F = 1.1 +/- 0.15, slightly larger than but consistent with the previous measurement from the Lyman \alpha\ forest autocorrelation. The simple linear model we use fails at separations r < 15 Mpc/h, and we show that this may reasonably be due to the enhanced ionization due to radiation from the quasars. We also provide the expected correction that the mass overdensity around the quasar implies for measurements of the ionizing radiation background from the line-of-sight proximity effect.

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We use a sample of 650 projected quasar pairs to study the HI Lya absorption transverse to luminous, z~2 quasars at proper separations of 30kpc < R < 1Mpc. In contrast to measurements along the line-of-sight, regions transverse to quasars exhibit enhanced HI Lya absorption and a larger variance than the ambient intergalactic medium, with increasing absorption and variance toward smaller scales. Analysis of composite spectra reveals excess absorption characterized by a Lya equivalent width profile W = 2.3A (R/100kpc)^-0.46. We also observe a high (~60%) covering factor of strong, optically thick HI absorbers (HI column log NHI > 17.3) at separations R<200kpc, which decreases to ~20% at R~1Mpc, but still represents a significant excess over the cosmic average. This excess of optically thick absorption can be described by a quasar-absorber cross-correlation function xi_QA(r) = (r/r_0)^gamma with a large correlation length r_0 = 12.5+2.7-1.4 Mpc/h (comoving) and gamma = 1.68+0.14-0.30. The HI absorption measured around quasars exceeds that of any previously studied population, consistent with quasars being hosted by massive dark matter halos Mhalo~10^12.5 Msun at z~2.5. The environments of these massive halos are highly biased towards producing optically thick gas, and may even dominate the cosmic abundance of Lyman limit systems and hence the intergalactic opacity to ionizing photons at z~2.5. The anisotropic absorption around quasars implies the transverse direction is much less likely to be illuminated by ionizing radiation than the line-of-sight, which we interpret in terms of the same obscuration effects frequently invoked in unified models of active galactic nuclei.

We searched quasar spectra from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) for the rare occurrences where a strong damped Lyman-alpha absorber (DLA) blocks the Broad Line Region emission from the quasar and acts as a natural coronagraph to reveal narrow Ly\alpha\ emission from the host galaxy. We define a statistical sample of 31 DLAs in Data Release 9 (DR9) with log N(HI) > 21.3 cm^-2 located at less than 1500 km s^-1 from the quasar redshift. In 25% (8) of these DLAs, a strong narrow Ly\alpha\ emission line is observed with flux ~25 x 10^-17 erg s^-1 cm^-2 on average. For DLAs without this feature in their troughs, the average 3-\sigma\ upper limit is < 0.8 x 10^-17 erg s^-1 cm^-2. Our statistical sample is nearly 2.5 times larger than the anticipated number of intervening DLAs in DR9 within 1500 km s^-1 of the quasar redshift. We also define a sample of 26 DLAs from DR9 and DR10 with narrow Ly\alpha\ emission detected and no limit on the HI column density to better characterize properties of the host galaxy emission. Analyzing the statistical sample, we do not find substantial differences in the kinematics, metals, or reddening for the two populations with and without emission detected. The highly symmetric narrow Ly\alpha\ emission line profile centered in the HI trough indicates that the emitting region is separate from the absorber. The luminosity of the narrow Ly\alpha\ emission peaks is intermediate between that of Lyman-alpha emitters and radio galaxies, implying that the Ly\alpha\ emission is predominantly due to ionizing radiation from the AGN. Galaxies neighboring the quasar host are likely responsible for the majority (> 75%) of these DLAs, with only a minority (< 25%) arising from HI clouds located in the AGN host galaxy.

At low redshifts, dust-obscured quasars often have strong yet narrow permitted lines in the rest-frame optical and ultraviolet, excited by the central active nucleus, earning the designation Type II quasars. We present a sample of 145 candidate Type II quasars at redshifts between 2 and 4.3, encompassing the epoch at which quasar activity peaked in the universe. These objects, selected from the quasar sample of the Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey III, are characterized by weak continuum in the rest-frame ultraviolet (typical continuum magnitude of i \approx 22) and strong lines of CIV and Ly \alpha, with Full Width at Half Maximum less than 2000 kms-1. The continuum magnitudes correspond to an absolute magnitude of -23 or brighter at redshift 3, too bright to be due exclusively to the host galaxies of these objects. Roughly one third of the objects are detected in the shorter-wavelength bands of the WISE survey; the spectral energy distributions (SEDs) of these objects appear to be intermediate between classic Type I and Type II quasars seen at lower redshift. Five objects are detected at rest frame 6\mu m by Spitzer, implying bolometric luminosities of several times 10^46 erg s-1. We have obtained polarization measurements for two objects; they are roughly 3% polarized. We suggest that these objects are luminous quasars, with modest dust extinction (A_V ~ 0.5 mag), whose ultraviolet continuum also includes a substantial scattering contribution. Alternatively, the line of sight to the central engines of these objects may be partially obscured by optically thick material.

We perform halo occupation distribution (HOD) modeling of the projected two-point correlation function (2PCF) of high-redshift (z~1.2) X-ray-bright active galactic nuclei (AGN) in the XMM-COSMOS field measured by Allevato et al. The HOD parameterization is based on low-luminosity AGN in cosmological simulations. At the median redshift of z~1.2, we derive a median mass of (1.02+0.21/-0.23)x10^{13} Msun/h for halos hosting central AGN and an upper limit of ~10% on the AGN satellite fraction. Our modeling results indicate (at the 2.5-sigma level) that X-ray AGN reside in more massive halos compared to more bolometrically luminous, optically-selected quasars at similar redshift. The modeling also yields constraints on the duty cycle of the X-ray AGN, and we find that at z~1.2 the average duration of the X-ray AGN phase is two orders of magnitude longer than that of the quasar phase. Our inferred mean occupation function of X-ray AGN is similar to recent empirical measurements with a group catalog and suggests that AGN halo occupancy increases with increasing halo mass. We project the XMM-COSMOS 2PCF measurements to forecast the required survey parameters needed in future AGN clustering studies to enable higher precision HOD constraints and determinations of key physical parameters like the satellite fraction and duty cycle. We find that N^{2}/A~5x10^{6} deg^{-2} (with N the number of AGN in a survey area of A deg^{2}) is sufficient to constrain the HOD parameters at the 10% level, which is easily achievable by upcoming and proposed X-ray surveys.

We present spatially resolved long-slit spectroscopy from the Southern African Large Telescope (SALT) to examine the spatial extent of the narrow-line regions (NLRs) of a sample of 8 luminous obscured quasars at 0.10 < z < 0.43. Our results are consistent with an observed shallow slope in the relationship between NLR size and L_[OIII], which has been interpreted to indicate that NLR size is limited by the density and ionization state of the NLR gas rather than the availability of ionizing photons. We also explore how the NLR size scales with a more direct measure of instantaneous AGN power using mid-IR photometry from WISE, which probes warm to hot dust near the central black hole and so, unlike [OIII], does not depend on the properties of the NLR. Using our results as well as samples from the literature, we obtain a power-law relationship between NLR size and L_8micron that is significantly steeper than that observed for NLR size and L_[OIII]. We find that the size of the NLR goes approximately as L^(1/2)_8micron, as expected from the simple scenario of constant-density clouds illuminated by a central ionizing source. We further see tentative evidence for a flattening of the relationship between NLR size and L_8micron at the high luminosity end, and propose that we are seeing a limiting NLR size of 10 - 20 kpc, beyond which the availability of gas to ionize becomes too low. We find that L_[OIII] ~ L_8micron^(1.4), consistent with a picture in which the L_[OIII] is dependent on the volume of the NLR. These results indicate that high-luminosity quasars have a strong effect in ionizing the available gas in a galaxy.

Unification by orientation is a ubiquitous concept in the study of active galactic nuclei. A gold standard of the orientation paradigm is the hypothesis that radio galaxies and radio-loud quasars are intrinsically the same, but are observed over different ranges of viewing angles. Historically, strong support for this model was provided by the projected sizes of radio structure in luminous radio galaxies, which were found to be significantly larger than those of quasars, as predicted due to simple geometric projection. Recently, this test of the simplest prediction of orientation-based models has been revisited with larger samples that cover wider ranges of fundamental properties---and no clear difference in projected sizes of radio structure is found. Cast solely in terms of viewing angle effects, these results provide convincing evidence that unification of these objects solely through orientation fails. However, it is possible that conflicting results regarding the role orientation plays in our view of radio sources simply result from insufficient sampling of their intrinsic size distribution. We test this possibility using Monte-Carlo simulations constrained by real sample sizes and properties. We develop models for the real intrinsic size distribution of radio sources, simulate observations by randomly sampling intrinsic sizes and viewing angles, and analyze how likely each sample is to support or dispute unification by orientation. We find that, while it is possible to reconcile conflicting results purely within a simple, orientation-based framework, it is very unlikely. We analyze the effects that sample size, relative numbers of radio galaxies and quasars, the critical angle that separates the two subclasses, and the shape of the intrinsic size distribution have on this type of test.

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.

The flat spectrum radio quasar (FSRQ) PKS 0208-512 underwent three outbursts at the optical-near-infrared (OIR) wavelengths during 2008-2011. The second OIR outburst did not have a gamma-ray counterpart despite being comparable in brightness and temporal extent to the other two. We model the time variable spectral energy distribution of PKS 0208-512 during those three flaring episodes with leptonic models to investigate the physical mechanism that can produce this anomalous flare. We show that the redder-when-brighter spectral trend in the OIR bands can be explained by the superposition of a fixed thermal component from the accretion disk and a synchrotron component of fixed shape and variable normalization. We estimate the accretion disk luminosity at L_d ~8 X 10^45 erg/s. Using the observed variability timescale in the OIR band t_{var,obs} ~2 d and the X-ray luminosity L_X ~3.5 X 10^45 erg/s, we constrain the location of the emitting region to distance scales that are broadly comparable with the dusty torus. We show that variations in the Compton dominance parameter by a factor of ~4 --- which may result in the anomalous outburst --- can be relatively easily accounted for by moderate variations in the magnetic field strength or the location of the emission region. Since such variations appear to be rare among FSRQs, we propose that most gamma-ray/OIR flares in these objects are produced in jet regions where the magnetic field and external photon fields vary similarly with distance along the jet, e.g., u_B ~u_ext ~r^{-2}.

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.

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-alpha forest flux fluctuations at a redshift z=2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard LCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100x(alpha_iso-1) = -1.6 ^{+2.0+4.3+7.4}_{-2.0-4.1-6.8} (stat.) +/- 1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100x(alpha_par-1) = -1.3 ^{+3.5+7.6 +12.3}_{-3.3-6.7-10.2} (stat.) +/- 2.0 (syst.) and 100x(alpha_perp-1) = -2.2 ^{+7.4+17}_{-7.1-15} +/- 3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ~ -0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-alpha forest sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the \Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range $2.1\le z \le 3.5$ from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift $z=2.3$, we measure the monopole and quadrupole components of the correlation function for separations in the range $20\hMpc<r<200\hMpc$. A peak in the correlation function is seen at a separation equal to $(1.01\pm0.03)$ times the distance expected for the BAO peak within a concordance $\Lambda$CDM cosmology. This first detection of the BAO peak at high redshift, when the universe was strongly matter dominated, results in constraints on the angular diameter distance $\da$ and the expansion rate $H$ at $z=2.3$ that, combined with priors on $H_0$ and the baryon density, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background (CMB) observations, this result implies $H(z=2.3)=(224\pm8){\rm km\,s^{-1}Mpc^{-1}}$, indicating that the time derivative of the cosmological scale parameter $\dot{a}=H(z=2.3)/(1+z)$ is significantly greater than that measured with BAO at $z\sim0.5$. This demonstrates that the expansion was decelerating in the range $0.7<z<2.3$, as expected from the matter domination during this epoch. Combined with measurements of $H_0$, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe.

We present a set of 3-dimensional, radiation-magnetohydrodynamic calculations of the gravitational collapse of massive (300 Msun), star-forming molecular cloud cores. We show that the combined effects of magnetic fields and radiative feedback strongly suppress core fragmentation, leading to the production of single star systems rather than small clusters. We find that the two processes are efficient at suppressing fragmentation in different regimes, with the feedback most effective in the dense, central region and the magnetic field most effective in more diffuse, outer regions. Thus, the combination of the two is much more effective at suppressing fragmentation than either one considered in isolation. Our work suggests that typical massive cores, which have mass-to-flux ratios of about 2 relative to critical, likely form a single star system, but that cores with weaker fields may form a small star cluster. This result helps us understand why the observed relationship between the core mass function and the stellar initial mass function holds even for ~100 Msun cores with many thermal Jeans masses of material. We also demonstrate that a ~40 AU Keplerian disk is able to form in our simulations, despite the braking effect caused by the strong magnetic field.

We describe fitting methods developed to analyze fluctuations in the Lyman-{\alpha} forest and measure the parameters of baryon acoustic oscillations (BAO). We apply our methods to BOSS Data Release 9. Our method is based on models of the three-dimensional correlation function in physical coordinate space, and includes the effects of redshift-space distortions, anisotropic non-linear broadening, and broadband distortions. We allow for independent scale factors along and perpendicular to the line of sight to minimize the dependence on our assumed fiducial cosmology and to obtain separate measurements of the BAO angular and relative velocity scales. Our fitting software and the input files needed to reproduce our main BOSS Data Release 9 results are publicly available.

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 measurement of the two-point cross-correlation function (CCF) of 8,198 Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) quasars and 349,608 DR10 CMASS galaxies from the Baryonic Oscillation Spectroscopic Survey (BOSS) at redshift <z>~0.5 (0.3<z<0.9). The cross-correlation function can be reasonably well fit by a power-law model xi_QG(r)=(r/r_0)^(-gamma) on projected scales of r_p=2-25 Mpc/h with r_0=6.61+-0.25 Mpc/h and gamma=1.69+-0.07. We estimate a quasar linear bias of b_Q=1.38+-0.10 at <z>=0.53 from the CCF measurements. This linear bias corresponds to a characteristic host halo mass of ~4x10^12 M_sun/h, compared to ~10^13 M_sun/h characteristic host halo mass for CMASS galaxies. We divide the quasar sample in luminosity and constrain the luminosity dependence of quasar bias to be db_Q/dlogL=0.20+-0.34 or 0.11+-0.32 (depending on different luminosity divisions) for quasar luminosities -23.5>M_i(z=2)>-25.5, implying a weak luminosity dependence of quasar clustering for the bright end of the quasar population at <z>~0.5. We compare our measurements with theoretical predictions, Halo Occupation Distribution (HOD) models and mock catalogs. These comparisons suggest quasars reside in a broad range of host halos, and the host halo mass distributions significantly overlap with each other for quasars at different luminosities, implying a poor correlation between halo mass and instantaneous quasar luminosity. We also find that the quasar HOD parameterization is largely degenerate such that different HODs can reproduce the CCF equally well, but with different outcomes such as the satellite fraction and host halo mass distribution. These results highlight the limitations and ambiguities in modeling the distribution of quasars with the standard HOD approach and the need for additional information in populating quasars in dark matter halos with HOD. [Abridged]

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.

If broad absorption line (BAL) quasars represent a high covering fraction evolutionary state (even if this is not the sole factor governing the presence of BALs), it is expected that they should show an excess of mid-infrared radiation compared to normal quasars. Some previous studies have suggested that this is not the case. We perform the first analysis of the IR properties of radio-loud BAL quasars, using IR data from WISE and optical (rest-frame ultraviolet) data from SDSS, and compare the BAL quasar sample with a well-matched sample of unabsorbed quasars. We find a statistically significant excess in the mid- to near-infrared luminosities of BAL quasars, particularly at rest-frame wavelengths of 1.5 and 4 microns. Our sample was previously used to show that BALs are observed along many lines of sight towards quasars, but with an overabundance of more edge-on sources, suggesting that orientation factors into the appearance of BALs. The evidence here---of a difference in IR luminosities between BAL quasars and unabsorbed quasars---may be ascribed to evolution. This suggests that a merging of the current BAL paradigms is needed to fully describe the class.

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

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

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.

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

We present the first results from an ongoing survey for Damped Lyman-alpha systems (DLAs) in the spectra of z>2 quasars observed in the course of the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey (SDSS) III. Our full (non-statistical) sample, based on Data Release 9, comprises 12,081 systems with log N(HI)>=20, out of which 6,839 have log N(HI)>=20.3. This is the largest DLA sample ever compiled, superseding that from SDSS-II by a factor of seven. Using a statistical sub-sample and estimating systematics from realistic mock data, we probe the N(HI) distribution at <z> = 2.5. Contrary to what is generally believed, the distribution extends beyond 10^22 cm^-2 with a moderate slope of index\approx-3.5. This result matches surprisingly well the opacity-corrected distribution observed at z = 0. The cosmological mass density of neutral gas in DLAs is found to be Omega_g_DLA~10^-3, evolving only mildly over the past 12 billion years.

The temperature of the gas in molecular clouds is a key determinant of the characteristic mass of star formation. Ambipolar diffusion (AD) is considered one of the most important heating mechanisms in weakly ionized molecular clouds. In this work, we study the AD heating rate using 2-fluid turbulence simulations and compare it with the overall heating rate due to turbulent dissipation. We find that for observed molecular clouds, which typically have Alfven Mach numbers of ~1 (Crutcher 1999) and AD Reynolds numbers of ~20 (McKee et al. 2010), about 70% of the total turbulent dissipation is in the form of AD heating. AD has an important effect on the length scale where energy is dissipated: when AD heating is strong, most of the energy in the cascade is removed by ion-neutral drift, with a comparatively small amount of energy making it down to small scales. We derive a relation for the AD heating rate that describes the results of our simulations to within a factor of two. Turbulent dissipation, including AD heating, is generally less important that cosmic-ray heating in molecular clouds, although there is substantial scatter in both.

(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 present 21 examples of C IV Broad Absorption Line (BAL) trough disappearance in 19 quasars selected from systematic multi-epoch observations of 582 bright BAL quasars (1.9 < z < 4.5) by the Sloan Digital Sky Survey-I/II (SDSS-I/II) and SDSS-III. The observations span 1.1-3.9 yr rest-frame timescales, longer than have been sampled in many previous BAL variability studies. On these timescales, ~2.3% of C IV BAL troughs disappear and ~3.3% of BAL quasars show a disappearing trough. These observed frequencies suggest that many C IV BAL absorbers spend on average at most a century along our line of sight to their quasar. Ten of the 19 BAL quasars showing C IV BAL disappearance have apparently transformed from BAL to non-BAL quasars; these are the first reported examples of such transformations. The BAL troughs that disappear tend to be those with small-to-moderate equivalent widths, relatively shallow depths, and high outflow velocities. Other non-disappearing C IV BALs in those nine objects having multiple troughs tend to weaken when one of them disappears, indicating a connection between the disappearing and non-disappearing troughs, even for velocity separations as large as 10000-15000 km/s. We discuss possible origins of this connection including disk-wind rotation and changes in shielding gas.