Search

Dark Energy Spectroscopic Instrument (DESI) uses more than 2.4 million Emission Line Galaxies (ELGs) for 3D large-scale structure (LSS) analyses in its Data Release 1 (DR1). Such large statistics enable thorough research on systematic uncertainties. In this study, we focus on spectroscopic systematics of ELGs. The redshift success rate ($f_{\rm goodz}$) is the relative fraction of secure redshifts among all measurements. It depends on observing conditions, thus introduces non-cosmological variations to the LSS. We, therefore, develop the redshift failure weight ($w_{\rm zfail}$) and a per-fibre correction ($\eta_{\rm zfail}$) to mitigate these dependences. They have minor influences on the galaxy clustering. For ELGs with a secure redshift, there are two subtypes of systematics: 1) catastrophics (large) that only occur in a few samples; 2) redshift uncertainty (small) that exists for all samples. The catastrophics represent 0.26\% of the total DR1 ELGs, composed of the confusion between O\,\textsc{ii} and sky residuals, double objects, total catastrophics and others. We simulate the realistic 0.26\% catastrophics of DR1 ELGs, the hypothetical 1\% catastrophics, and the truncation of the contaminated $1.31

We present a Neural Network based emulator for the galaxy redshift-space power spectrum that enables several orders of magnitude acceleration in the galaxy clustering parameter inference, while preserving 3$\sigma$ accuracy better than 0.5\% up to $k_{\mathrm{max}}$=0.25$h^{-1}Mpc$ within $\Lambda$CDM and around 0.5\% $w_0$-$w_a$CDM. Our surrogate model only emulates the galaxy bias-invariant terms of 1-loop perturbation theory predictions, these terms are then combined analytically with galaxy bias terms, counter-terms and stochastic terms in order to obtain the non-linear redshift space galaxy power spectrum. This allows us to avoid any galaxy bias prescription in the training of the emulator, which makes it more flexible. Moreover, we include the redshift $z \in [0,1.4]$ in the training which further avoids the need for re-training the emulator. We showcase the performance of the emulator in recovering the cosmological parameters of $\Lambda$CDM by analysing the suite of 25 AbacusSummit simulations that mimic the DESI Luminous Red Galaxies at $z=0.5$ and $z=0.8$, together as the Emission Line Galaxies at $z=0.8$. We obtain similar performance in all cases, demonstrating the reliability of the emulator for any galaxy sample at any redshift in $0 < z < 1.4$

Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained challenging. We present a simulation-based model that is able to capture the cosmological dependence of the full shape of the density-split clustering (DSC) statistics down to intra-halo scales. Our models are based on neural-network emulators that are trained on high-fidelity mock galaxy catalogues within an extended-$\Lambda$CDM framework, incorporating the effects of redshift-space, Alcock-Paczynski distortions and models of the halo-galaxy connection. Our models reach sub-percent level accuracy down to $1\,h^{-1}{\rm Mpc}$ and are robust against different choices of galaxy-halo connection modelling. When combined with the galaxy 2PCF, DSC can tighten the constraints on $\omega_{\rm cdm}$, $\sigma_8$, and $n_s$ by factors of 2.9, 1.9, and 2.1, respectively, compared to a 2PCF-only analysis. DSC additionally puts strong constraints on environment-based assembly bias parameters. Our code is made publicly available on Github., Comment: Submitted to MNRAS. Source code to generate the figures available in the captions. Updated to add missing references

We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1\,h^{-1}{\rm Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-$\Lambda$CDM framework, including redshift-space and Alcock-Paczynski distortions. Our base-$\Lambda$CDM analysis finds $\omega_{\rm cdm} = 0.1201\pm 0.0022$, $\sigma_8 = 0.792\pm 0.034$, and $n_s = 0.970\pm 0.018$, corresponding to $f\sigma_8 = 0.462\pm 0.020$ at $z \approx 0.525$, which is in agreement with Planck 2018 predictions and various clustering studies in the literature. We test single-parameter extensions to base-$\Lambda$CDM, varying the running of the spectral index, the dark energy equation of state, and the density of massless relic neutrinos, finding no compelling evidence for deviations from the base model. We model the galaxy-halo connection using a halo occupation distribution framework, finding signatures of environment-based assembly bias in the data. We validate our pipeline against mock catalogues that match the clustering and selection properties of CMASS, showing that we can recover unbiased cosmological constraints even with a volume 84 times larger than the one used in this study., Comment: Submitted to MNRAS. Source code to generate the figures available in the captions. Updated to add missing references and fix legend of Fig. 6

We present an approach for accurate estimation of the covariance of 2-point correlation functions that requires fewer mocks than the standard mock-based covariance. This can be achieved by dividing a set of mocks into jackknife regions and fitting the correction term first introduced in Mohammad & Percival (2022), such that the mean of the jackknife covariances corresponds to the one from the mocks. This extends the model beyond the shot-noise limited regime, allowing it to be used for denser samples of galaxies. We test the performance of our fitted jackknife approach, both in terms of accuracy and precision, using lognormal mocks with varying densities and approximate EZmocks mimicking the DESI LRG and ELG samples in the redshift range of z = [0.8, 1.2]. We find that the Mohammad-Percival correction produces a bias in the 2-point correlation function covariance matrix that grows with number density and that our fitted jackknife approach does not. We also study the effect of the covariance on the uncertainty of cosmological parameters by performing a full-shape analysis. We find that our fitted jackknife approach based on 25 mocks is able to recover unbiased and as precise cosmological parameters as the ones obtained from a covariance matrix based on 1000 or 1500 mocks, while the Mohammad-Percival correction produces uncertainties that are twice as large. The number of mocks required to obtain an accurate estimation of the covariance for 2-point correlation function is therefore reduced by a factor of 40-60., Comment: 13 pages, 14 figures, accepted to MNRAS

In this series of papers we present an emulator-based halo model for the non-linear clustering of galaxies in modified gravity cosmologies. In the first paper, we present emulators for the following halo properties: the halo mass function, concentration-mass relation and halo-matter cross-correlation function. The emulators are trained on data extracted from the \textsc{FORGE} and \textsc{BRIDGE} suites of $N$-body simulations, respectively for two modified gravity (MG) theories: $f(R)$ gravity and the DGP model, varying three standard cosmological parameters $\Omega_{\mathrm{m0}}, H_0, \sigma_8$, and one MG parameter, either $\bar{f}_{R0}$ or $r_{\mathrm{c}}$. Our halo property emulators achieve an accuracy of $\lesssim 1\%$ on independent test data sets. We demonstrate that the emulators can be combined with a galaxy-halo connection prescription to accurately predict the galaxy-galaxy and galaxy-matter correlation functions using the halo model framework., Comment: 18 pages, 15 figures

The dependence of galaxy clustering on local density provides an effective method for extracting non-Gaussian information from galaxy surveys. The two-point correlation function (2PCF) provides a complete statistical description of a Gaussian density field. However, the late-time density field becomes non-Gaussian due to non-linear gravitational evolution and higher-order summary statistics are required to capture all of its cosmological information. Using a Fisher formalism based on halo catalogues from the Quijote simulations, we explore the possibility of retrieving this information using the density-split clustering (DS) method, which combines clustering statistics from regions of different environmental density. We show that DS provides more precise constraints on the parameters of the $\nu \Lambda$CDM model compared to the 2PCF, and we provide suggestions for where the extra information may come from. DS improves the constraints on the sum of neutrino masses by a factor of $7$ and by factors of 4, 3, 3, 6, and 5 for $\Omega_{\rm m}$, $\Omega_{\rm b}$, $h$, $n_s$, and $\sigma_8$, respectively. We compare DS statistics when the local density environment is estimated from the real or redshift-space positions of haloes. The inclusion of DS autocorrelation functions, in addition to the cross-correlation functions between DS environments and haloes, recovers most of the information that is lost when using the redshift-space halo positions to estimate the environment. We discuss the possibility of constructing simulation-based methods to model DS clustering statistics in different scenarios., Comment: Author's accepted manuscript

Over the next five years, the Dark Energy Spectroscopic Instrument (DESI) will use 10 spectrographs with 5000 fibers on the 4m Mayall Telescope at Kitt Peak National Observatory to conduct the first Stage-IV dark energy galaxy survey. At $z < 0.6$, the DESI Bright Galaxy Survey (BGS) will produce the most detailed map of the Universe during the dark energy dominated epoch with redshifts of >10 million galaxies over 14,000 deg$^2$. In this work, we present and validate the final BGS target selection and survey design. From the Legacy Surveys, BGS will target a $r < 19.5$ magnitude-limited sample (BGS Bright); a fainter $19.5 < r < 20.175$ sample, color-selected to have high redshift efficiency (BGS Faint); and a smaller low-z quasar sample. BGS will observe these targets using exposure times, scaled to achieve uniform completeness, and visit each point on the footprint three times. We use observations from the Survey Validation programs conducted prior to the main survey along with realistic simulations to show that BGS can complete its strategy and make optimal use of `bright' time. We demonstrate that BGS targets have stellar contamination <1% and that their densities do not depend strongly on imaging properties. We also confirm that BGS Bright will achieve >80% fiber assignment efficiency. Finally, we show that BGS Bright and Faint will achieve >95% redshift success rates with no significant dependence on observing conditions. BGS meets the requirements for an extensive range of scientific applications. BGS will yield the most precise Baryon Acoustic Oscillations and Redshift-Space Distortions measurements at $z < 0.4$. It also presents opportunities to exploit new methods that require highly complete and dense galaxy samples (e.g. N-point statistics, multi-tracers). BGS further provides a powerful tool to study galaxy populations and the relations between galaxies and dark matter., Comment: AJ, submitted, 34 pages, 22 figures, one of a suite of 8 papers detailing targeting for DESI

In this series of papers, we present a simulation-based model for the non-linear clustering of galaxies based on separate modelling of clustering in real space and velocity statistics. In the first paper, we present an emulator for the real-space correlation function of galaxies, whereas the emulator of the real-to-redshift space mapping based on velocity statistics is presented in the second paper. Here, we show that a neural network emulator for real-space galaxy clustering trained on data extracted from the Dark Quest suite of N-body simulations achieves sub-per cent accuracies on scales $1 < r < 30 $ $h^{-1} \,\mathrm{Mpc}$, and better than $3\%$ on scales $r < 1$ $h^{-1}\mathrm{Mpc}$ in predicting the clustering of dark-matter haloes with number density $10^{-3.5}$ $(h^{-1}\mathrm{Mpc})^{-3}$, close to that of SDSS LOWZ-like galaxies. The halo emulator can be combined with a galaxy-halo connection model to predict the galaxy correlation function through the halo model. We demonstrate that we accurately recover the cosmological and galaxy-halo connection parameters when galaxy clustering depends only on the mass of the galaxies' host halos. Furthermore, the constraining power in $\sigma_8$ increases by about a factor of $2$ when including scales smaller than $5$ $h^{-1} \,\mathrm{Mpc}$. However, when mass is not the only property responsible for galaxy clustering, as observed in hydrodynamical or semi-analytic models of galaxy formation, our emulator gives biased constraints on $\sigma_8$. This bias disappears when small scales ($r < 10$ $h^{-1}\mathrm{Mpc}$) are excluded from the analysis. This shows that a vanilla halo model could introduce biases into the analysis of future datasets.

The use of realistic mock galaxy catalogues is essential in the preparation of large galaxy surveys, in order to test and validate theoretical models and to assess systematics. We present an updated version of the mock catalogue constructed from the Millennium-XXL simulation, which uses a halo occupation distribution (HOD) method to assign galaxies r-band magnitudes and g-r colours. We have made several modifications to the mock to improve the agreement with measurements from the SDSS and GAMA surveys. We find that cubic interpolation, which was used to build the original halo lightcone, produces extreme velocities between snapshots. Using linear interpolation improves the correlation function quadrupole measurements on small scales. We also update the g-r colour distributions so that the observed colours better agree with measurements from GAMA data, particularly for faint galaxies. As an example of the science that can be done with the mock, we investigate how the luminosity function depends on environment and colour, and find good agreement with measurements from the GAMA survey. This full-sky mock catalogue is designed for the ongoing Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS), and is complete to a magnitude limit r=20.2., Comment: 15 pages, 11 figures. Accepted for publication in MNRAS

Realistic lightcone mocks are important in the clustering analyses of large galaxy surveys. For simulations where only the snapshots are available, it is common to create approximate lightcones by joining together the snapshots in spherical shells. We assess the two-point clustering measurements of central galaxies in approximate lightcones built from the Millennium-XXL simulation, which are constructed using different numbers of snapshots. The monopole and quadrupole of the real-space correlation function is strongly boosted on small scales below 1 Mpc/h, due to some galaxies being duplicated at the boundaries between snapshots in the lightcone. When more snapshots are used, the total number of duplicated galaxies is approximately constant, but they are pushed to smaller separations. The effect of this in redshift space is small, as long as the snapshots are cut into shells in real space. Randomly removing duplicated galaxies is able to reduce the excess clustering signal. Including satellite galaxies will reduce the impact of the duplicates, since many small-scale pairs come from satellites in the same halo. Galaxies that are missing from the lightcone at the boundaries can be added to the lightcone by having a small overlap between each shell. This effect will impact analyses that use very small-scale clustering measurements, and when using mocks to test the impact of fibre collisions., Comment: 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

We investigate the large-scale clustering of the final spectroscopic sample of quasars from the recently completed extended Baryon Oscillation Spectroscopic Survey (eBOSS). The sample contains $343708$ objects in the redshift range $0.8

We characterise the selection cuts and clustering properties of a magnitude-limited sample of bright galaxies that is part of the Bright Galaxy Survey (BGS) of the Dark Energy Spectroscopic Instrument (DESI) using the ninth data release of the Legacy Imaging Surveys (DR9). We describe changes in the DR9 selection compared to the DR8 one as explored in Ruiz-Macias et al. (2021). We also compare the DR9 selection in three distinct regions: BASS/MzLS in the north Galactic Cap (NGC), DECaLS in the NGC, and DECaLS in the south Galactic Cap (SGC). We investigate the systematics associated with the selection and assess its completeness by matching the BGS targets with the Galaxy and Mass Assembly (GAMA) survey. We measure the angular clustering for the overall bright sample (r $\leq$ 19.5) and as function of apparent magnitude and colour. This enables to determine the clustering strength and slope by fitting a power-law model that can be used to generate accurate mock catalogues for this tracer. We use a counts-in-cells technique to explore higher-order statistics and cross-correlations with external spectroscopic data sets in order to check the evolution of the clustering with redshift and the redshift distribution of the BGS targets using clustering-redshifts. While this work validates the properties of the BGS bright targets, the final target selection pipeline and clustering properties of the entire DESI BGS will be fully characterised and validated with the spectroscopic data of Survey Validation., Comment: 16 pages, 17 figures

Shortly after its discovery, General Relativity (GR) was applied to predict the behavior of our Universe on the largest scales, and later became the foundation of modern cosmology. Its validity has been verified on a range of scales and environments from the Solar system to merging black holes. However, experimental confirmations of GR on cosmological scales have so far lacked the accuracy one would hope for -- its applications on those scales being largely based on extrapolation and its validity sometimes questioned in the shadow of the unexpected cosmic acceleration. Future astronomical instruments surveying the distribution and evolution of galaxies over substantial portions of the observable Universe, such as the Dark Energy Spectroscopic Instrument (DESI), will be able to measure the fingerprints of gravity and their statistical power will allow strong constraints on alternatives to GR. In this paper, based on a set of $N$-body simulations and mock galaxy catalogs, we study the predictions of a number of traditional and novel estimators beyond linear redshift distortions in two well-studied modified gravity models, chameleon $f(R)$ gravity and a braneworld model, and the potential of testing these deviations from GR using DESI. These estimators employ a wide array of statistical properties of the galaxy and the underlying dark matter field, including two-point and higher-order statistics, environmental dependence, redshift space distortions and weak lensing. We find that they hold promising power for testing GR to unprecedented precision. The major future challenge is to make realistic, simulation-based mock galaxy catalogs for both GR and alternative models to fully exploit the statistic power of the DESI survey and to better understand the impact of key systematic effects. Using these, we identify future simulation and analysis needs for gravity tests using DESI., Comment: 69 pages, 36 figures, 6 tables; accepted for publication in JCAP. The modified gravity HOD galaxy catalogues used in this paper can be downloaded from http://icc.dur.ac.uk/data/, by searching for the title of this paper from the dataset table

The Dark Energy Spectroscopic Instrument (DESI) will execute a nearly magnitude-limited survey of low redshift galaxies ($0.05 \leq z \leq 0.4$, median $z \approx 0.2$). Clustering analyses of this Bright Galaxy Survey (BGS) will yield the most precise measurements to date of baryon acoustic oscillations and redshift-space distortions at low redshift. DESI BGS will comprise two target classes: (i) BRIGHT ($r<19.5$~mag), and (ii) FAINT ($19.5

We search for the Baryon Acoustic Oscillations in the projected cross-correlation function binned into transverse comoving radius between the SDSS-IV DR16 eBOSS quasars and a dense photometric sample of galaxies selected from the DESI Legacy Imaging Surveys. We estimate the density of the photometric sample of galaxies in this redshift range to be about 2900 deg$^{-2}$, which is deeper than the official DESI ELG selection, and the density of the spectroscopic sample is about 20 deg$^{-2}$. In order to mitigate the systematics related to the use of different imaging surveys close to the detection limit, we use a neural network approach that accounts for complex dependencies between the imaging attributes and the observed galaxy density. We find that we are limited by the depth of the imaging surveys which affects the density and purity of the photometric sample and its overlap in redshift with the quasar sample, which thus affects the performance of the method. When cross-correlating the photometric galaxies with quasars in $0.6 \leq z \leq 1.2$, the cross-correlation function can provide better constraints on the comoving angular distance, $D_{\rm M}$ (6\% precision) compared to the constraint on the spherically-averaged distance $D_{\rm V}$ (9\% precision) obtained from the auto-correlation. Although not yet competitive, this technique will benefit from the arrival of deeper photometric data from upcoming surveys which will enable it to go beyond the current limitations we have identified in this work., Comment: 21 pages, 25 figures

We present the steps taken to produce a reliable and complete input galaxy catalogue for the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) using the photometric Legacy Survey DR8 DECam. We analyse some of the main issues faced in the selection of targets for the DESI BGS, such as star-galaxy separation, contamination by fragmented stars and bright galaxies. Our pipeline utilizes a new way to select BGS galaxies using Gaia photometry and we implement geometrical and photometric masks that reduce the number of spurious objects. The resulting catalogue is cross-matched with the Galaxy And Mass Assembly (GAMA) survey to assess the completeness of the galaxy catalogue and the performance of the target selection. We also validate the clustering of the sources in our BGS catalogue by comparing with mock catalogues and the Sloan Digital Sky Survey (SDSS) data. Finally, the robustness of the BGS selection criteria is assessed by quantifying the dependence of the target galaxy density on imaging and other properties. The largest systematic correlation we find is a 7 per cent suppression of the target density in regions of high stellar density., Comment: 23 pages, 18 Figures

We measure the anisotropic clustering of the quasar sample from Data Release 16 (DR16) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS). A sample of $343,708$ spectroscopically confirmed quasars between redshift $0.8

We measure the clustering of quasars of the final data release (DR16) of eBOSS. The sample contains $343\,708$ quasars between redshifts $0.8\leq z\leq2.2$ over $4699\,\mathrm{deg}^2$. We calculate the Legendre multipoles (0,2,4) of the anisotropic power spectrum and perform a BAO and a Full-Shape (FS) analysis at the effective redshift $z{\rm eff}=1.480$. The errors include systematic errors that amount to 1/3 of the statistical error. The systematic errors comprise a modelling part studied using a blind N-Body mock challenge and observational effects studied with approximate mocks to account for various types of redshift smearing and fibre collisions. For the BAO analysis, we measure the transverse comoving distance $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.60\pm{0.90}$ and the Hubble distance $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=13.34\pm{0.60}$. This agrees with the configuration space analysis, and the consensus yields: $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.69\pm{0.80}$ and $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=13.26\pm{0.55}$. In the FS analysis, we fit the power spectrum using a model based on Regularised Perturbation Theory, which includes Redshift Space Distortions and the Alcock-Paczynski effect. The results are $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.68\pm{0.90}$ and $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=13.52\pm{0.51}$ and we constrain the linear growth rate of structure $f(z_{\rm eff})\sigma_8(z_{\rm eff})=0.476\pm{0.047}$. Our results agree with the configuration space analysis. The consensus analysis of the eBOSS quasar sample yields: $D_{\rm M}(z_{\rm eff})/r_{\rm drag}=30.21\pm{0.79}$, $D_{\rm H}(z_{\rm eff})/r_{\rm drag}=3.23\pm{0.47}$ and $f(z_{\rm eff})\sigma_8(z_{\rm eff})=0.462\pm{0.045}$ and is consistent with a flat $\Lambda {\rm CDM}$ cosmological model using Planck results., Comment: 22 pages, 16 figures, 12 tables A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found here: https://sdss.org/science/cosmology-results-from-eboss/

We present the final Sloan Digital Sky Survey IV (SDSS-IV) quasar catalog from Data Release 16 of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). This catalog comprises the largest selection of spectroscopically confirmed quasars to date. The full catalog includes two sub-catalogs: a "superset" of all SDSS-IV/eBOSS objects targeted as quasars containing 1,440,615 observations and a quasar-only catalog containing 750,414 quasars, including 225,082 new quasars appearing in an SDSS data release for the first time, as well as known quasars from SDSS-I/II/III. We present automated identification and redshift information for these quasars alongside data from visual inspections for 320,161 spectra. The quasar-only catalog is estimated to be 99.8% complete with 0.3% to 1.3% contamination. Automated and visual inspection redshifts are supplemented by redshifts derived via principal component analysis and emission lines. We include emission line redshifts for H$\alpha$, H$\beta$, Mg II, C III], C IV, and Ly$\alpha$. Identification and key characteristics generated by automated algorithms are presented for 99,856 Broad Absorption Line quasars and 35,686 Damped Lyman Alpha quasars. In addition to SDSS photometric data, we also present multi-wavelength data for quasars from GALEX, UKIDSS, WISE, FIRST, ROSAT/2RXS, XMM-Newton, and Gaia. Calibrated digital optical spectra for these quasars can be obtained from the SDSS Science Archive Server., Comment: 32 pages, 13 figures, 6 tables. Accepted to ApJS. Catalog files are available at https://data.sdss.org/sas/dr16/eboss/qso/DR16Q/ . A summary of all SDSS BAO and RSD measurements with legacy figures can be found at https://www.sdss.org/science/final-bao-and-rsd-measurements/ while full cosmological interpretation of these can be found at https://www.sdss.org/science/cosmology-results-from-eboss/

We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Ly$\alpha$ forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, $r_d$, from eight different samples and six measurements of the growth rate parameter, $f\sigma_8$, from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, $\Lambda$CDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization the BAO data provide nearly an order of magnitude improvement on curvature constraints. The RSD measurements indicate a growth rate that is consistent with predictions from Planck primary data and with General Relativity. When combining the results of SDSS BAO and RSD with external data, all multiple-parameter extensions remain consistent with a $\Lambda$CDM model. Regardless of cosmological model, the precision on $\Omega_\Lambda$, $H_0$, and $\sigma_8$, remains at roughly 1\%, showing changes of less than 0.6\% in the central values between models. The inverse distance ladder measurement under a o$w_0w_a$CDM yields $H_0= 68.20 \pm 0.81 \, \rm km\, s^{-1} Mpc^{-1}$, remaining in tension with several direct determination methods. (abridged), Comment: A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://www.sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found here: https://www.sdss.org/science/cosmology-results-from-eboss/

We present large-scale structure catalogs from the completed extended Baryon Oscillation Spectroscopic Survey (eBOSS). Derived from Sloan Digital Sky Survey (SDSS) -IV Data Release 16 (DR16), these catalogs provide the data samples, corrected for observational systematics, and random positions sampling the survey selection function. Combined, they allow large-scale clustering measurements suitable for testing cosmological models. We describe the methods used to create these catalogs for the eBOSS DR16 Luminous Red Galaxy (LRG) and Quasar samples. The quasar catalog contains 343,708 redshifts with $0.8 < z < 2.2$ over 4,808\,deg$^2$. We combine 174,816 eBOSS LRG redshifts over 4,242\,deg$^2$ in the redshift interval $0.6 < z < 1.0$ with SDSS-III BOSS LRGs in the same redshift range to produce a combined sample of 377,458 galaxy redshifts distributed over 9,493\,deg$^2$. Improved algorithms for estimating redshifts allow that 98 per cent of LRG observations result in a successful redshift, with less than one per cent catastrophic failures ($\Delta z > 1000$ ${\rm km~s}^{-1}$). For quasars, these rates are 95 and 2 per cent (with $\Delta z > 3000$ ${\rm km~s}^{-1}$). We apply corrections for trends between the number densities of our samples and the properties of the imaging and spectroscopic data. For example, the quasar catalog obtains a $\chi^2$/DoF$= 776/10$ for a null test against imaging depth before corrections and a $\chi^2$/DoF$=6/8$ after. The catalogs, combined with careful consideration of the details of their construction found here-in, allow companion papers to present cosmological results with negligible impact from observational systematic uncertainties., Comment: Matches version accepted by MNRAS, very minor changes. A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://www.sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found here: https://sdss.org/science/cosmology-results-from-eboss/

The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock-Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a series of HOD mock catalogues generated using the OuterRim N-body simulation. We test three models on a series of non-blind mocks, in the OuterRim cosmology, and blind mocks, which have been rescaled to new cosmologies, and investigate the effects of redshift smearing and catastrophic redshifts. We find that for the non-blind mocks, the models are able to recover $f\sigma_8$ to within 3% and $\alpha_\parallel$ and $\alpha_\bot$ to within 1%. The scatter in the measurements is larger for the blind mocks, due to the assumption of an incorrect fiducial cosmology. From this mock challenge, we find that all three models perform well, with similar systematic errors on $f\sigma_8$, $\alpha_\parallel$ and $\alpha_\bot$ at the level of $\sigma_{f\sigma_8}=0.013$, $\sigma_{\alpha_\parallel}=0.012$ and $\sigma_{\alpha_\bot}=0.008$. The systematic error on the combined consensus is $\sigma_{f\sigma_8}=0.011$, $\sigma_{\alpha_\parallel}=0.008$ and $\sigma_{\alpha_\bot}=0.005$, which is used in the final DR16 analysis. For BAO fits in configuration and Fourier space, we take conservative systematic errors of $\sigma_{\alpha_\parallel}=0.010$ and $\sigma_{\alpha_\bot}=0.007$., Comment: 25 pages, 17 figures, 8 tables, updated to match the version accepted for publication in MNRAS. A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://sdss.org/science/final-bao-and-rsd-measurements/. The full cosmological interpretation of these measurements can be found here: https://sdss.org/science/cosmology-results-from-eboss/

To understand the nature of the accelerated expansion of the Universe, we need to combine constraints on the expansion rate and growth of structure. The growth rate is usually extracted from three dimensional galaxy maps by exploiting the effects of peculiar motions on galaxy clustering. However, theoretical models of the probability distribution function (PDF) of galaxy pairwise peculiar velocities are not accurate enough on small scales to reduce the error on theoretical predictions to the level required to match the precision expected for measurements from future surveys. Here, we improve the modelling of the pairwise velocity distribution by using the Skew-T PDF, which has nonzero skewness and kurtosis. Our model accurately reproduces the redshift-space multipoles (monopole, quadrupole and hexadecapole) predicted by N-body simulations, above scales of about $10\,h^{-1}{\rm Mpc}$. We illustrate how a Taylor expansion of the streaming model can reveal the contributions of the different moments to the clustering multipoles, which are independent of the shape of the velocity PDF. The Taylor expansion explains why the Gaussian streaming model works well in predicting the first two redshift-space multipoles, although the velocity PDF is non-Gaussian even on large scales. Indeed, any PDF with the correct first two moments would produce precise results for the monopole down to scales of about $10\,h^{-1}{\rm Mpc}$, and for the quadrupole down to about $30\,h^{-1}{\rm Mpc}$. An accurate model for the hexadecapole needs to include higher-order moments., Comment: Accepted for publication in MNRAS

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)., Comment: DR16 release: Monday Dec 9th 2019. This is the alphabetical order SDSS-IV collaboration data release paper. 25 pages, 6 figures, accepted by ApJS on 11th May 2020. Minor changes clarify or improve text and figures relative to v1

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

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., Comment: 25 pages, 22 figures

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

We measure the growth rate and its evolution using the anisotropic clustering of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample, which includes $148\,659$ quasars covering the wide redshift range of $0.8 < z < 2.2$ and a sky area of $2112.90$ $\rm deg^2$. To optimise measurements we deploy a redshift-dependent weighting scheme, which allows us to avoid binning, and perform the data analysis consistently including the redshift evolution across the sample. We perform the analysis in Fourier space, and use the redshift evolving power spectrum multipoles to measure the redshift space distortion parameter $f\sigma_8$ and parameters controlling the anisotropic projection of the cosmological perturbations. We measure $f \sigma_8(z=1.52)=0.43 \pm 0.05 $ and $df\sigma_8/dz (z=1.52)= - 0.16 \pm 0.08$, consistent with the expectation for a $\Lambda$CDM cosmology as constrained by the Planck experiment., Comment: 11 pages

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., Comment: 37 pages, 25 figures, 14 tables. Data and covariance will be made public after the acceptance of the paper. Submitted to MNRAS

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

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., Comment: accepted by A&A

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., Comment: SDSS-IV collaboration alphabetical author data release paper. DR14 happened on 31st July 2017. 19 pages, 5 figures. Accepted by ApJS on 28th Nov 2017 (this is the "post-print" and "post-proofs" version; minor corrections only from v1, and most of errors found in proofs corrected)

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., Comment: Accepted by MNRAS; BAO distance likelihood available in source files 'QSOv1.9fEZmock_BAOchi2.dat'; full set of data to be public eventually from SDSS website

We study the first year of the eBOSS quasar sample in the redshift range $0.9

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., Comment: Published in Astronomical Journal

In current and future surveys, quasars play a key role. The new data will extend our knowledge of the Universe as it will be used to better constrain the cosmological model at redshift $z>1$ via baryon acoustic oscillation and redshift space distortion measurements. Here, we present the first clustering study of quasars observed by the extended Baryon Oscillation Spectroscopic Survey. We measure the clustering of $\sim 70,000$ quasars located in the redshift range $0.9