In this paper we study the universality of the results of SQUEzE, a software package to classify quasar spectra and estimate their redshifts. The code is presented in \cite{Perez-Rafols+2019}. We test the results against changes on signal-to-noise, spectral resolution, wavelength coverage, and quasar brightness. We find that SQUEzE levels of performance (quantified with purity and completeness) are stable to spectra that have a noise dispersion 4 times that of our standard test sample, BOSS. We also find that the performance remains unchanged if pixels of width 25\AA are considered, and decreases by $\sim2\%$ for pixels of width 100\AA. We see no effect when analyzing subsets of different quasar brightness, and we establish that the blue part (up to 7000\AA) of the spectra is sufficient for the classification. Finally, we compare our suite of tests with samples of spectra expected from WEAVE-QSO and DESI, and narrow-band imaging from J-PAS. We conclude that SQUEzE will perform similarly when confronted with the demands of these future surveys as when applied to current BOSS (and eBOSS) data. {\it Keywords: cosmology: observations - quasar: emission lines - quasar: absorption lines}

We present SQUEzE, a software package to classify quasar spectra and estimate their redshifts. SQUEzE is a random forest classifier operating on the parameters of candidate emission peaks identified in the spectra. We test the performance of the algorithm using visually inspected data from BOSS as a truth table. Only 4\% of the sample ($\sim$6,800 quasars and $\sim$11,520 contaminants) is needed for converged training in recommended choices of the confidence threshold ($0.2<p_{\rm min}<0.7$). For an operational mode which balances purity and completeness ($p_{\rm min}=0.28$) we recover a purity of $96.81\pm0.39\%$ ($99.30\pm0.14\%$ for quasars with $z \geq 2.1$) and a completeness of $96.83\pm0.30\%$ ($98.42\pm0.15\%$ for quasars with $z \geq 2.1$). SQUEzE can be used to obtain a $\approx$100\% pure sample of $z \geq 2.1$. quasars (with $\approx$96\% completeness) by using a confidence threshold of $p_{\rm min}=0.7$. The estimated redshift error is $1,500{\rm \thinspace km/s}$ and we recommend that SQUEzE be used in conjunction with an additional step of redshift tuning to achieve maximum precision. We find that SQUEzE achieves the necessary performance to replace visual inspection in BOSS-like spectroscopic surveys of quasars with subsequent publications in this series exploring expectations for future surveys and alternative methods. Keywords: cosmology: observations - quasar: emission lines - quasar: absorption lines

The statistical power of Lyman-${\alpha}$ forest Baryon Acoustic Oscillation (BAO) measurements is set to increase significantly in the coming years as new instruments such as the Dark Energy Spectroscopic Instrument deliver progressively more constraining data. Generating mock datasets for such measurements will be important for validating analysis pipelines and evaluating the effects of systematics. With such studies in mind, we present LyaCoLoRe: a package for producing synthetic Lyman-${\alpha}$ forest survey datasets for BAO analyses. LyaCoLoRe transforms initial Gaussian random field skewers into skewers of transmitted flux fraction via a number of fast approximations. In this work we explain the methods of producing mock datasets used in LyaCoLoRe, and then measure correlation functions on a suite of realisations of such data. We demonstrate that we are able to recover the correct BAO signal, as well as large-scale bias parameters similar to literature values. Finally, we briefly describe methods to add further astrophysical effects to our skewers - high column density systems and metal absorbers - which act as potential complications for BAO analyses.

We measure the imprint of primordial baryon acoustic oscillations (BAO) in the correlation function of Ly$\alpha$ absorption in quasar spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS) in Data Release 14 (DR14) of the Sloan Digital Sky Survey (SDSS)-IV. In addition to 179,965 spectra with absorption in the Lyman-$\alpha$ (Ly$\alpha$) region, we use, for the first time, Ly$\alpha$ absorption in the Lyman-$\beta$ region of 56,154 spectra. We measure the Hubble distance, $D_H$, and the comoving angular diameter distance, $D_M$, relative to the sound horizon at the drag epoch $r_d$ at an effective redshift $z=2.34$. Using a physical model of the correlation function outside the BAO peak, we find $D_H(2.34)/r_d=8.86\pm 0.29$ and $D_M(2.34)/r_d=37.41\pm 1.86$, within 1$\sigma$ from the flat-$\Lambda$CDM model consistent with CMB anisotropy measurements. With the addition of polynomial "broadband" terms, the results remain within one standard deviation of the CMB-inspired model. Combined with the quasar-Ly$\alpha$ cross-correlation measurement presented in a companion paper Blomqvist19, the BAO measurements at $z=2.35$ are within 1.7$\sigma$ of the predictions of this model.

We present a measurement of the baryon acoustic oscillation (BAO) scale at redshift $z=2.35$ from the three-dimensional correlation of Lyman-$\alpha$ (Ly$\alpha$) forest absorption and quasars. The study uses 266,590 quasars in the redshift range $1.77<z<3.5$ from the Sloan Digital Sky Survey (SDSS) Data Release 14 (DR14). The sample includes the first two years of observations by the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS), providing new quasars and re-observations of BOSS quasars for improved statistical precision. Statistics are further improved by including Ly$\alpha$ absorption occurring in the Ly$\beta$ wavelength band of the spectra. From the measured BAO peak position along and across the line of sight, we determined the Hubble distance $D_{H}$ and the comoving angular diameter distance $D_{M}$ relative to the sound horizon at the drag epoch $r_{d}$: $D_{H}(z=2.35)/r_{d}=9.20\pm 0.36$ and $D_{M}(z=2.35)/r_{d}=36.3\pm 1.8$. These results are consistent at $1.5\sigma$ with the prediction of the best-fit spatially-flat cosmological model with the cosmological constant reported for the Planck (2016) analysis of cosmic microwave background anisotropies. Combined with the Ly$\alpha$ auto-correlation measurement presented in a companion paper, the BAO measurements at $z=2.34$ are within $1.7\sigma$ of the predictions of this model.

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 measure the cosmological bias factor of DLAs from their cross-correlation with the Ly$\alpha$ forest absorption, as a function of the DLA metal strength, defined from an average of equivalent widths of the strongest detectable low-ionization metal lines. A clear increase of the bias factor with metal strength is detected, as expected from a relation of metallicity and velocity dispersion with host halo mass. The relation is stronger after the metal strength is corrected for the HI column density, to make it more related to metallicity instead of metal column density. After correcting for the effects of measurement errors of the metal strength parameter, we find that the bias factor of DLAs with the weakest metal lines is close to unity, consistent with an origin in dwarf galaxies with host halo masses $\sim10^{10} M_{\odot}$, whereas the most metal rich DLAs have a bias factor as large as $b_{\rm DLA} \sim 3$, indicative of massive galaxies or galaxy groups in host halos with masses $\sim 10^{12} M_{\odot}$. Our result confirms the physical origin of the relation of bias factors measured from cross-correlation studies to the host halos of the absorbers.

We present the first paper of the series Origin of Metals around Galaxies (OMG) aimed to study the origin of the metals observed in the circumgalactic and intergalactic media. In this work we extract and build the catalogs of metal absorbers that will be used in future analyses, and make our results publicly available to the community. We design a fully automatic algorithm to search for absorption metal-line doublets of the species CIV, NV, SiIV and MgII in high-resolution ($R\gtrsim30\,000$) quasar spectra without human intervention, and apply it to the high-resolution and signal-to-noise ratio spectra of 690 quasars, observed with the UVES and HIRES instruments. We obtain $5\,656$ CIV doublets, $7\,919$ doublets of MgII, $2\,258$ of SiIV, and 239 of NV, constituting the largest high-resolution metal-doublet samples to date, and estimate the dependence of their completeness and purity on various doublet parameters such as equivalent width and redshift, using real and artificial quasar spectra. The catalogs include doublets with rest-frame line equivalent widths down to a few ${\rm m\AA}$, all detected at a significance above 3$\sigma$, and covering the redshifts between $1<z \lesssim 5$, properties that make them useful for a wide range of chemical evolution studies.

We present measurements of the cross-correlation of the triply-ionized carbon (CIV) forest with quasars using Sloan Digital Sky Survey Data Release 14. The study exploits a large sample of new quasars from the first two years of observations by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The CIV forest is a weaker tracer of large-scale structure than the Ly$\alpha$ forest, but benefits from being accessible at redshifts $z<2$ where the quasar number density from eBOSS is high. Our data sample consists of 287,651 CIV forest quasars in the redshift range $1.4<z<3.5$ and 387,315 tracer quasars with $1.2<z<3.5$. We measure large-scale correlations from CIV absorption occuring in three distinct quasar rest-frame wavelength bands of the spectra referred to as the CIV forest, the SiIV forest and the Ly$\alpha$ forest. From the combined fit to the quasar-CIV cross-correlations for the CIV forest and the SiIV forest, the CIV redshift-space distortion parameter is $\beta_{\rm CIV}=0.27_{\ -0.14}^{\ +0.16}$ and its combination with the CIV linear transmission bias parameter is $b_{\rm CIV}(1+\beta_{\rm CIV})=-0.0183_{\ -0.0014}^{\ +0.0013}$ ($1\sigma$ statistical error) at the mean redshift $z=2.00$. Splitting the sample at $z=2.2$ to constrain the bias evolution with redshift yields the power-law exponent $\gamma=0.60\pm0.63$, indicating a significantly weaker redshift-evolution than for the Ly$\alpha$ forest linear transmission bias. We demonstrate that CIV absorption has the potential to be used as a probe of baryon acoustic oscillations (BAO). While the current data set is insufficient for a detection of the BAO peak feature, the final quasar samples for redshifts $1.4<z<2.2$ from eBOSS and the Dark Energy Spectroscopic Instrument (DESI) are expected to provide measurements of the isotropic BAO scale to $\sim7\%$ and $\sim3\%$ precision, respectively, at $z\simeq1.6$.

We present a measurement of baryon acoustic oscillations (BAO) in the cross-correlation of quasars with the Ly$\alpha$-forest flux-transmission at a mean redshift $z=2.40$. The measurement uses the complete SDSS-III data sample: 168,889 forests and 234,367 quasars from the SDSS Data Release DR12. In addition to the statistical improvement on our previous study using DR11, we have implemented numerous improvements at the analysis level allowing a more accurate measurement of this cross-correlation. We also developed the first simulations of the cross-correlation allowing us to test different aspects of our data analysis and to search for potential systematic errors in the determination of the BAO peak position. We measure the two ratios $D_{H}(z=2.40)/r_{d} = 9.01 \pm 0.36$ and $D_{M}(z=2.40)/r_{d} = 35.7 \pm 1.7$, where the errors include marginalization over the non-linear velocity of quasars and the metal - quasar cross-correlation contribution, among other effects. These results are within $1.8\sigma$ of the prediction of the flat-$\Lambda$CDM model describing the observed CMB anisotropies. We combine this study with the Ly$\alpha$-forest auto-correlation function [2017A&A...603A..12B], yielding $D_{H}(z=2.40)/r_{d} = 8.94 \pm 0.22$ and $D_{M}(z=2.40)/r_{d} = 36.6 \pm 1.2$, within $2.3\sigma$ of the same flat-$\Lambda$CDM model.

We present a measurement of the DLA mean bias from the cross-correlation of DLA and the Ly$\alpha$ forest, updating earlier results of Font-Ribera et al. 2012 with the final BOSS Data Release and an improved method to address continuum fitting corrections. Our cross-correlation is well fitted by linear theory with the standard $\Lambda CDM$ model, with a DLA bias of $b_{\rm DLA} = 1.99\pm 0.11$; a more conservative analysis, which removes DLA in the Ly$\beta$ forest and uses only the cross-correlation at $r> 10{\rm h^{-1}\,Mpc}$, yields $b_{\rm DLA} = 2.00\pm 0.19$. This assumes the cosmological model from \cite{Planck2015} and the Ly$\alpha$ forest bias factors of Bautista et al. 2017, and includes only statistical errors obtained from bootstrap analysis. The main systematic errors arise from possible impurities and selection effects in the DLA catalogue, and from uncertainties in the determination of the Ly$\alpha$ forest bias factors and a correction for effects of high column density absorbers. We find no dependence of the DLA bias on column density or redshift. The measured bias value corresponds to a host halo mass $\sim 4\cdot10^{11} {\rm M_{\odot}}$ if all DLA were hosted in halos of a similar mass. In a realistic model where host halos over a broad mass range have a DLA cross section $\Sigma(M_h) \propto M_h^{\alpha}$ down to $M_h > M_{\rm min} =10^{8.5} {\rm M_{\odot}}$, we find that $\alpha > 1$ is required to have $b_{\rm DLA}> 1.7$, implying a steeper relation or higher value of $M_{\rm min}$ than is generally predicted in numerical simulations of galaxy formation.

We study the mean absorption spectrum of the Damped Lyman alpha population at $z\sim 2.6$ by stacking normalized, rest-frame shifted spectra of $\sim 27\,000$ DLAs from the DR12 of BOSS/SDSS-III. We measure the equivalent widths of 50 individual metal absorption lines in 5 intervals of DLA hydrogen column density, 5 intervals of DLA redshift, and overall mean equivalent widths for an additional 13 absorption features from groups of strongly blended lines. The mean equivalent width of low-ionization lines increases with $N_{\rm HI}$, whereas for high-ionization lines the increase is much weaker. The mean metal line equivalent widths decrease by a factor $\sim 1.1-1.5$ from $z\sim2.1$ to $z \sim 3.5$, with small or no differences between low- and high-ionization species. We develop a theoretical model, inspired by the presence of multiple absorption components observed in high-resolution spectra, to infer mean metal column densities from the equivalent widths of partially saturated metal lines. We apply this model to 14 low-ionization species and to AlIII, SIII, SiIII, CIV, SiIV, NV and OVI. We use an approximate derivation for separating the equivalent width contributions of several lines to blended absorption features, and infer mean equivalent widths and column densities from lines of the additional species NI, ZnII, CII${}^{*}$, FeIII, and SIV. Several of these mean column densities of metal lines in DLAs are obtained for the first time; their values generally agree with measurements of individual DLAs from high-resolution, high signal-to-noise ratio spectra when they are available.

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.

The existence of a cosmic neutrino background can be probed indirectly by CMB experiments, not only by measuring the background density of radiation in the universe, but also by searching for the typical signatures of the fluctuations of free-streaming species in the temperature and polarisation power spectrum. Previous studies have already proposed a rather generic parametrisation of these fluctuations, that could help to discriminate between the signature of ordinary free-streaming neutrinos, or of more exotic dark radiation models. Current data are compatible with standard values of these parameters, which seems to bring further evidence for the existence of a cosmic neutrino background. In this work, we investigate the robustness of this conclusion under various assumptions. We generalise the definition of an effective sound speed and viscosity speed to the case of massive neutrinos or other dark radiation components experiencing a non-relativistic transition. We show that current bounds on these effective parameters do not vary significantly when considering an arbitrary value of the particle mass, or extended cosmological models with a free effective neutrino number, dynamical dark energy or a running of the primordial spectrum tilt. We conclude that it is possible to make a robust statement about the detection of the cosmic neutrino background by CMB experiments.

We present a measurement of the cross-correlation of MgII absorption and massive galaxies, using the DR11 main galaxy sample of the Baryon Oscillation Spectroscopic Survey of SDSS-III (CMASS galaxies), and the DR7 quasar spectra of SDSS-II. The cross-correlation is measured by stacking quasar absorption spectra shifted to the redshift of galaxies that are within a certain impact parameter bin of the quasar, after dividing by a quasar continuum model. This results in an average MgII equivalent width as a function of impact parameter from a galaxy, ranging from 50 kpc to more than 10 Mpc in proper units, which includes all MgII absorbers. We show that special care needs to be taken to use an unbiased quasar continuum estimator, to avoid systematic errors in the measurement of the mean stacked MgII equivalent width. The measured cross-correlation follows the expected shape of the galaxy correlation function, although measurement errors are large. We use the cross-correlation amplitude to derive the bias factor of MgII absorbers, finding bMgII = 2.33 \pm? 0.19, where the error accounts only for the statistical uncertainty in measuring the mean equivalent width. This bias factor is larger than that obtained in previous studies and may be affected by modeling uncertainties that we discuss, but if correct it suggests that MgII absorbers at redshift z \simeq 0:5 are spatially distributed on large scales similarly to the CMASS galaxies in BOSS. Keywords: galaxies: haloes, galaxies: formation, quasars: absorption lines, large-scale structure of universe

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.

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.