Search

Mesoscale convection covers an intermediate scale range between small-scale turbulence and the global organization of the convection flow. It is often characterized by an order of the convection patterns despite very high Rayleigh numbers and strong turbulent fluctuations. In this review, we discuss several aspects of mesoscale convection, which have been investigated by three-dimensional direct numerical simulations. The numerical studies are performed in a characteristic configuration of a plane layer that is heated from below and cooled from above or subject to constant heat flux at the top and bottom boundaries. We discuss the role of the thermal and mechanical boundary conditions for structure formation and study the impact of the domain shape as well as the Prandtl number. With respect to the latter, we focus on low values that arise in astrophysical convection and are partly not anymore accessible in laboratory experiments with liquid metals. Beside these experiments in the Boussinesq approximation, we report studies of non-Boussinesq mesoscale convection. This is done by investigating effects of compressibility and temperature dependence of material properties. The kinetic energy dissipation rate turns out to remain a central quantity for the turbulent mixing in compressible convection. Their different components, statistics, relation to the turbulent viscosity, and the multifractal properties are discussed., Comment: 26 pages, 16 figures (some in reduced quality)

We measure the growth of cosmic density fluctuations on large scales and across the redshift range $0.3

The phonon thermal transport properties of twisted bi-layer graphene are investigated using lattice dynamics and the Boltzmann transport equation. The thermal conductivities of 13.2 and 21.8 twisted configurations are 56% and 36% lower than the untwisted configuration which has a room temperature thermal conductivity of 2260 W/m-K. The reason for this lower thermal conductivity is unraveled from phonon mode-level analysis made possible through the untwisting of layers. Due to a large commensurate unitcell of twisted configuration, the Brillouin zone is folded but this folding has no bearing on the phonon scattering phase space. The major impact of twisting is felt by flexural phonons with out-of-plane vibrations via the change in strength of flexural interatomic interactions and since these flexural phonons carry majority of the heat in bilayer graphene (65% at room temperature), the thermal conductivity is sensitive to layer twisting. Our study suggests that twisting will affect thermal transport only for those materials that have a major contribution from flexural phonon modes.

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a Halo Occupation Distribution (HOD) model upon the \textsc{FastPM} simulations, such that the resulting galaxy clustering reproduces high resolution $N$-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements -- monopole and quadrupole -- to a precision required by the DESI Year 1 Emission Line Galaxy sample down to non-linear scales, i.e. $k<0.5\,h\mathrm{Mpc}$ or $s>10\,\mathrm{Mpc}/h$. Furthermore, we compute covariance matrices based on the resulting \textsc{FastPM} galaxy clustering -- monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate [e.g. power spectrum] on the covariance matrix of the Fourier counterpart [e.g. correlation function]. We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between two to ten per cent Consequently, using the current pipeline we can precisely reproduce the clustering of $N$-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies., Comment: 21 pages, 18 figures, submitted to MNRAS, Zenodo link https://doi.org/10.5281/zenodo.8185822

We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter $\fnl$. Our sample comprises over 12 million targets, covering 14,000 square degrees of the sky, with redshifts in the range $0.2< z < 1.35$. We identify Galactic extinction, survey depth, and astronomical seeing as the primary sources of systematic error, and employ linear regression and artificial neural networks to alleviate non-cosmological excess clustering on large scales. Our methods are tested against simulations with and without $\fnl$ and systematics, showing superior performance of the neural network treatment. The neural network with a set of nine imaging property maps passes our systematic null test criteria, and is chosen as the fiducial treatment. Assuming the universality relation, we find $\fnl = 34^{+24(+50)}_{-44(-73)}$ at 68\%(95\%) confidence. We apply a series of robustness tests (e.g., cuts on imaging, declination, or scales used) that show consistency in the obtained constraints. We study how the regression method biases the measured angular power-spectrum and degrades the $\fnl$ constraining power. The use of the nine maps more than doubles the uncertainty compared to using only the three primary maps in the regression. Our results thus motivate the development of more efficient methods that avoid over-correction, protect large-scale clustering information, and preserve constraining power. Additionally, our results encourage further studies of $\fnl$ with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics and lessen the degradation in the $\fnl$ uncertainty., Comment: 21 pages, 17 figures, 7 tables (Appendix excluded). Published in MNRAS

The One-Percent survey of the Dark Energy Spectroscopic Instrument collected ~ 270k emission line galaxies (ELGs) at 0.8 < z < 1.6. The high completeness of the sample allowed the clustering to be measured down to scales never probed before, 0.04 Mpc/h in rp for the projected 2-point correlation function (2PCF) and 0.17 Mpc/h in galaxy pair separation s for the 2PCF monopole and quadrupole. The most striking feature of the measurements is a strong signal at the smallest scales, below 0.2 Mpc/h in rp and 1 Mpc/h in s. We analyze these data in the halo occupation distribution framework. We consider different distributions for central galaxies, a standard power law for satellites with no condition on the presence of a central galaxy and explore several extensions of these models. For all considered models, the mean halo mass of the sample is found to be log10 ~ 11.9. We obtain a satellite mean occupation function which agrees with physically motivated ELG models only if we introduce central-satellite conformity, meaning that the satellite occupation is conditioned by the presence of central galaxies of the same type. To achieve in addition a good modeling of the clustering between 0.1 and 1 Mpc/h in rp, we allow for ELG positioning outside of the halo virial radius and find 0.5% of ELGs residing in the outskirts of halos. Furthermore, the satellite velocity dispersion inside halos is found to be ~ 30% larger than that of the halo dark matter particles. These are the main findings of our work. We investigate assembly bias as a function of halo concentration, local density or local density anisotropies and observe no significant change in our results. We split the data sample in two redshift bins and report no significant evolution with redshift. Lastly, changing the cosmology in the modeling impacts only slightly our results.

We present the probabilistic stellar mass function (pSMF) of galaxies in the DESI Bright Galaxy Survey (BGS), observed during the One-Percent Survey. The One-Percent Survey was one of DESI's survey validation programs conducted from April to May 2021, before the start of the main survey. It used the same target selection and similar observing strategy as the main survey and successfully observed the spectra and redshifts of 143,017 galaxies in the $r < 19.5$ magnitude-limited BGS Bright sample and 95,499 galaxies in the fainter surface brightness and color selected BGS Faint sample over $z < 0.6$. We derive pSMFs from posteriors of stellar mass, $M_*$, inferred from DESI photometry and spectroscopy using the Hahn et al. (2022a; arXiv:2202.01809) PRObabilistic Value-Added BGS (PROVABGS) Bayesian SED modeling framework. We use a hierarchical population inference framework that statistically and rigorously propagates the $M_*$ uncertainties. Furthermore, we include correction weights that account for the selection effects and incompleteness of the BGS observations. We present the redshift evolution of the pSMF in BGS as well as the pSMFs of star-forming and quiescent galaxies classified using average specific star formation rates from PROVABGS. Overall, the pSMFs show good agreement with previous stellar mass function measurements in the literature. Our pSMFs showcase the potential and statistical power of BGS, which in its main survey will observe >100$\times$ more galaxies. Moreover, we present the statistical framework for subsequent population statistics measurements using BGS, which will characterize the global galaxy population and scaling relations at low redshifts with unprecedented precision., Comment: 25 pages, 12 figures; data used to generate figures is available at https://doi.org/10.5281/zenodo.8018936; submitted to ApJ

We present the first comprehensive Halo Occupation Distribution (HOD) analysis of the DESI One-Percent survey Luminous Red Galaxy (LRG) and Quasi-Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 Mpc/h in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of fsat = 11+-1%, a mean halo mass of log10 Mh = 13.40+0.02-0.02, and a linear bias of blin = 1.93+0.06-0.04. For LRGs in 0.6 < z < 0.8, we find fsat = 14+-1%, log10 Mh = 13.24+0.02-0.02, and blin = 2.08+0.03-0.03. For QSOs, we infer fsat = 3+8-2%, log10 Mh = 12.65+0.09-0.04, and blin = 2.63+0.37-0.26 in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high-fidelity galaxy mocks. We also study the redshift-evolution of the DESI LRG sample from z = 0.4 up to z = 1.1, revealing significant and interesting trends in mean halo mass, linear bias, and satellite fraction., Comment: Submitted to MNRAS, comments welcome

As a key driver of the secular evolution of disc galaxies, bar formation is potentially linked to the surrounding tidal field. We systematically investigate the dependence of bars on both the small (${<}2\,\mathrm{Mpc}/h$) and large-scale (${>}5\,\mathrm{Mpc}/h$) tidal fields using galaxies observed between $0.01{<}z{<}0.11$ by the Sloan Digital Sky Survey (SDSS). We characterise bar strength using the ellipticity of the isophote that corresponds to the bar, $e_{\mathrm{bar}}$, derived from the galaxy image after subtracting the 2D disc component. We demonstrate the efficacy of our bar detection method by performing an extensive comparison with the visual identifications from SDSS and the DESI Legacy Surveys. Using the Yang et al. SDSS group catalogue, we confirm the results from a recent study that the average $e_{\mathrm{bar}}$ of galaxies within interacting clusters is higher than that within isolated ones at $0.01{<}z{<}0.06$, but this small-scale tidal enhancement of bars disappears after we increase the cluster sample by a factor of five to $z{=}0.11$. On large scales, we explore the dependence of $e_{\mathrm{bar}}$ on $\alpha_{5}$, the tidal anisotropy of the density field defined over $5\,\mathrm{Mpc}/h$. We do not detect any such dependence for $98\%$ of the galaxies with $\alpha_{5}{<}10$. Intriguingly, among the $2\%$ with $\alpha_{5}{\ge}10$, we detect some hint of a boost in bar strength in the underdense regions and a suppression in the overdense regions. Combining our results on both scales, we conclude that there is little evidence for the tidal dependence of bar formation in the local Universe, except for the extremely anisotropic environments., Comment: 16 pages, 14 figures. Comments are welcome

Renormalization enables a systematic scale-by-scale analysis of multiscale systems. In this paper, we employ \textit{renormalization group} (RG) to the shell model of turbulence and show that the RG equation is satisfied by $ |u_n|^2 =K_\mathrm{Ko} \epsilon^{2/3} k_n^{-2/3}$ and $ \nu_n = \nu_* \sqrt{K_\mathrm{Ko}} \epsilon^{1/3} k_n^{-4/3}$, where $k_n, u_n $ are the wavenumber and velocity of shell $ n $; $\nu_*, K_\mathrm{Ko}$ are RG and Kolmogorov's constants; and $ \epsilon $ is the energy dissipation rate. We find that $\nu_* \approx 0.5$ and $K_\mathrm{Ko} \approx 1.7$, consistent with earlier RG works on Navier-Stokes equation. We verify the theoretical predictions using numerical simulations., Comment: To appear in Phys. Rev. E

We constrain models of the galaxy distribution in the cosmic web using data from the Galaxy and Mass Assembly (GAMA) survey. We model the redshift-space behaviour of the 2-point correlation function (2pcf) and the recently proposed Voronoi volume function (VVF) -- which includes information beyond 2-point statistics. We extend the standard halo model using extra satellite degrees of freedom and two assembly bias parameters, $\alpha_{\rm cen}$ and $\alpha_{\rm sat}$, which respectively correlate the occupation numbers of central and satellite galaxies with their host halo's tidal environment. We measure $\alpha_{\rm sat}=1.44^{+0.25}_{-0.43}$ and $\alpha_{\rm cen}=-0.79^{+0.29}_{-0.11}$ using a combination of 2pcf and VVF measurements, representing a detection of assembly bias at the 3.3$\sigma$ (2.4$\sigma$) significance level for satellite (central) galaxies. This result remains robust to possible anisotropies in the halo-centric distribution of satellites as well as technicalities of estimating the data covariance. We show that the growth rate ($f\sigma_8$) deduced using models with assembly bias is about 7\% (i.e. $1.5\sigma$) lower than if assembly bias is ignored. When projected onto the $\Omega_m$-$\sigma_8$ plane, the model constraints without assembly bias overlap with Planck expectations, while allowing assembly bias introduces significant tension with Planck, preferring either a lower $\Omega_m$ or a lower $\sigma_8$. Finally, we find that the all-galaxy weak lensing signal is unaffected by assembly bias, but the central and satellite sub-populations individually show significantly different signals in the presence of assembly bias. Our results illustrate the importance of accurately modelling galaxy formation for cosmological inference from future surveys., Comment: 18 pages, 14 figures, 3 tables, minor changes in text and one additional table matches the accepted version

We present the first detection of the baryon acoustic oscillations (BAO) signal obtained using unblinded data collected during the initial two months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). From a selected sample of 261,291 Luminous Red Galaxies spanning the redshift interval 0.4 < z < 1.1 and covering 1651 square degrees with a 57.9% completeness level, we report a ~5 sigma level BAO detection and the measurement of the BAO location at a precision of 1.7%. Using a Bright Galaxy Sample of 109,523 galaxies in the redshift range 0.1 < z < 0.5, over 3677 square degrees with a 50.0% completeness, we also detect the BAO feature at ~3 sigma significance with a 2.6% precision. These first BAO measurements represent an important milestone, acting as a quality control on the optimal performance of the complex robotically-actuated, fiber-fed DESI spectrograph, as well as an early validation of the DESI spectroscopic pipeline and data management system. Based on these first promising results, we forecast that DESI is on target to achieve a high-significance BAO detection at sub-percent precision with the completed 5-year survey data, meeting the top-level science requirements on BAO measurements. This exquisite level of precision will set new standards in cosmology and confirm DESI as the most competitive BAO experiment for the remainder of this decade., Comment: 17 pages, 10 figures, 4 tables. MNRAS, 525, 5406. Author accepted manuscript matching the published version

We present a novel approach to the construction of mock galaxy catalogues for large-scale structure analysis based on the distribution of dark matter halos obtained with effective bias models at the field level. We aim to produce mock galaxy catalogues capable of generating accurate covariance matrices for a number of cosmological probes that are expected to be measured in current and forthcoming galaxy redshift surveys (e.g. two- and three-point statistics). We use the bias assignment method (BAM) to model the statistics of halo distribution through a learning algorithm using a few detailed $N$-body simulations, and approximated gravity solvers based on Lagrangian perturbation theory. Using specific models of halo occupation distributions, we generate galaxy mocks with the expected number density and central-satellite fraction of emission-line galaxies, which are a key target of the DESI experiment. BAM generates mock catalogues with per cent accuracy in a number of summary statistics, such as the abundance, the two- and three-point statistics of halo distributions, both in real and redshift space. In particular, the mock galaxy catalogues display $\sim 3\%-10\%$ accuracy in the multipoles of the power spectrum up to scales of $k\sim 0.4\,h^{-1}{\rm Mpc}$. We show that covariance matrices of two- and three-point statistics obtained with BAM display a similar structure to the reference simulation. BAM offers an efficient way to produce mock halo catalogues with accurate two- and three-point statistics, and is able to generate a variety of multi-tracer catalogues with precise covariance matrices of several cosmological probes. We discuss future developments of the algorithm towards mock production in DESI and other galaxy-redshift surveys. (Abridged), Comment: Accepted for publication at A&A

We describe the spectroscopic data processing pipeline of the Dark Energy Spectroscopic Instrument (DESI), which is conducting a redshift survey of about 40 million galaxies and quasars using a purpose-built instrument on the 4-m Mayall Telescope at Kitt Peak National Observatory. The main goal of DESI is to measure with unprecedented precision the expansion history of the Universe with the Baryon Acoustic Oscillation technique and the growth rate of structure with Redshift Space Distortions. Ten spectrographs with three cameras each disperse the light from 5000 fibers onto 30 CCDs, covering the near UV to near infrared (3600 to 9800 Angstrom) with a spectral resolution ranging from 2000 to 5000. The DESI data pipeline generates wavelength- and flux-calibrated spectra of all the targets, along with spectroscopic classifications and redshift measurements. Fully processed data from each night are typically available to the DESI collaboration the following morning. We give details about the pipeline's algorithms, and provide performance results on the stability of the optics, the quality of the sky background subtraction, and the precision and accuracy of the instrumental calibration. This pipeline has been used to process the DESI Survey Validation data set, and has exceeded the project's requirements for redshift performance, with high efficiency and a purity greater than 99 percent for all target classes., Comment: AJ, revised version, 55 pages, 55 figures, 4 tables

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a 5-year survey that aims to measure the redshifts of tens of millions of galaxies and quasars, including 8 million luminous red galaxies (LRGs) in the redshift range of $0.4

The Dark Energy Spectroscopic Instrument (DESI) will precisely constrain cosmic expansion and the growth of structure by collecting $\sim$40 million extra-galactic redshifts across $\sim$80\% of cosmic history and one third of the sky. The Emission Line Galaxy (ELG) sample, which will comprise about one-third of all DESI tracers, will be used to probe the Universe over the $0.6 < z < 1.6$ range, which includes the $1.1

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

We have studied the galaxy-group cross-correlations in redshift space for the Galaxy And Mass Assembly (GAMA) Survey. We use a set of mock GAMA galaxy and group catalogues to develop and test a novel 'halo streaming' model for redshift-space distortions. This treats 2-halo correlations via the streaming model, plus an empirical 1-halo term derived from the mocks, allowing accurate modelling into the nonlinear regime. In order to probe the robustness of the growth rate inferred from redshift-space distortions, we divide galaxies by colour, and divide groups according to their total stellar mass, calibrated to total mass via gravitational lensing. We fit our model to correlation data, to obtain estimates of the perturbation growth rate, $f\sigma_8$, validating parameter errors via the dispersion between different mock realizations. In both mocks and real data, we demonstrate that the results are closely consistent between different subsets of the group and galaxy populations, considering the use of correlation data down to some minimum projected radius, $r_{\rm min}$. For the mock data, we can use the halo streaming model to below $r_{\rm min} = 5h^{-1}$ Mpc, finding that all subsets yield growth rates within about 3% of each other, and consistent with the true value. For the actual GAMA data, the results are limited by cosmic variance: $f\sigma_8=0.29\pm 0.10$ at an effective redshift of 0.20; but there is every reason to expect that this method will yield precise constraints from larger datasets of the same type, such as the DESI bright galaxy survey., Comment: 17 pages, 12 figures. Matched to the accepted version

The Dark Energy Spectroscopic Instrument (DESI) has embarked on an ambitious five-year survey to explore the nature of dark energy with spectroscopy of 40 million galaxies and quasars. DESI will determine precise redshifts and employ the Baryon Acoustic Oscillation method to measure distances from the nearby universe to z > 3.5, as well as measure the growth of structure and probe potential modifications to general relativity. In this paper we describe the significant instrumentation we developed for the DESI survey. The new instrumentation includes a wide-field, 3.2-deg diameter prime-focus corrector that focuses the light onto 5020 robotic fiber positioners on the 0.812 m diameter, aspheric focal surface. The positioners and their fibers are divided among ten wedge-shaped petals. Each petal is connected to one of ten spectrographs via a contiguous, high-efficiency, nearly 50 m fiber cable bundle. The ten spectrographs each use a pair of dichroics to split the light into three channels that together record the light from 360 - 980 nm with a resolution of 2000 to 5000. We describe the science requirements, technical requirements on the instrumentation, and management of the project. DESI was installed at the 4-m Mayall telescope at Kitt Peak, and we also describe the facility upgrades to prepare for DESI and the installation and functional verification process. DESI has achieved all of its performance goals, and the DESI survey began in May 2021. Some performance highlights include RMS positioner accuracy better than 0.1", SNR per \sqrt{\AA} > 0.5 for a z > 2 quasar with flux 0.28e-17 erg/s/cm^2/A at 380 nm in 4000s, and median SNR = 7 of the [OII] doublet at 8e-17 erg/s/cm^2 in a 1000s exposure for emission line galaxies at z = 1.4 - 1.6. We conclude with highlights from the on-sky validation and commissioning of the instrument, key successes, and lessons learned. (abridged), Comment: 78 pages, 32 figures, submitted to AJ

Cosmological $N$-body simulations provide numerical predictions of the structure of the Universe against which to compare data from ongoing and future surveys, but the growing volume of the Universe mapped by surveys requires correspondingly lower statistical uncertainties in simulations, usually achieved by increasing simulation sizes at the expense of computational power. It was recently proposed to reduce simulation variance without incurring additional computational costs by adopting fixed-amplitude initial conditions. This method has been demonstrated not to introduce bias in various statistics, including the two-point statistics of galaxy samples typically used for extracting cosmological parameters from galaxy redshift survey data, but requires us to revisit current methods for estimating covariance matrices of clustering statistics for simulations. In this work, we find that it is not trivial to construct covariance matrices analytically for fixed-amplitude simulations, but we demonstrate that EZmock (Effective Zel'dovich approximation mock catalogue), the most efficient method for constructing mock catalogues with accurate two- and three-point statistics, provides reasonable covariance matrix estimates for such simulations. We further examine how the variance suppression obtained by amplitude-fixing depends on three-point clustering, small-scale clustering, and galaxy bias, and propose intuitive explanations for the effects we observe based on the EZmock bias model., Comment: 9 pages, 7 figures; in v2 we incorporate minor modifications from peer review and copyediting

Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogs to validate data analysis pipelines and identify potential systematics. We compare three $N$-body simulation codes, ABACUS, GADGET, and SWIFT, to investigate the regimes in which their results agree. We run $N$-body simulations at three different mass resolutions, $6.25\times10^{8}$, $2.11\times10^{9}$, and $5.00\times10^{9}~h^{-1}$M$_{\odot}$, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the DESI statistical error for $s > 20\, h^{-1}$Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of $2.1\times10^{9}~h^{-1}$M$_{\odot}$ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than $20 \, h^{-1}$Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1% for $k \leq 10~h$Mpc$^{-1}$. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-$N$-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy-halo relation will be presented in a follow up study., Comment: 18 pages, 17 figures. Submitted to MNRAS

This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys., Comment: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

The imprints of large-scale structures on the Cosmic Microwave Background can be studied via the CMB lensing and Integrated Sachs-Wolfe (ISW) signals. In particular, the stacked ISW signal around supervoids has been claimed in several works to be anomalously high. In this study, we find cluster and void superstructures using four tomographic redshift bins with $0

We investigate the 3D spin alignment of galaxies with respect to the large-scale filaments using the MaNGA survey. The cosmic web is reconstructed from the Sloan Digital Sky Survey using Disperse and the 3D spins of MaNGA galaxies are estimated using the thin disk approximation with integral field spectroscopy kinematics. Late-type spiral galaxies are found to have their spins parallel to the closest filament's axis. The alignment signal is found to be dominated by low-mass spirals. Spins of S0-type galaxies tend to be oriented preferentially in perpendicular direction with respect to the filament's axis. This orthogonal orientation is found to be dominated by S0s that show a notable misalignment between their kinematic components of stellar and ionised gas velocity fields and/or by low mass S0s with lower rotation support compared to their high mass counterparts. Qualitatively similar results are obtained when splitting galaxies based on the degree of ordered stellar rotation, such that galaxies with high spin magnitude have their spin aligned, and those with low spin magnitude in perpendicular direction to the filaments. In the context of conditional tidal torque theory, these findings suggest that galaxies' spins retain memory of their larger-scale environment. In agreement with measurements from hydrodynamical cosmological simulations, the measured signal at low redshift is weak, yet statistically significant. The dependence of the spin-filament orientation of galaxies on their stellar mass, morphology and kinematics highlights the importance of sample selection to detect the signal., Comment: 9 pages, 5 figures, 1 table. Accepted for publication in MNRAS

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

We use data from the DESI Legacy Survey imaging to probe the galaxy density field in tomographic slices covering the redshift range $0

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 test general relativity (GR) at the effective redshift $\bar{z} \sim 1.5$ by estimating the statistic $E_G$, a probe of gravity, on cosmological scales $19 - 190\,h^{-1}{\rm Mpc}$. This is the highest-redshift and largest-scale estimation of $E_G$ so far. We use the quasar sample with redshifts $0.8 < z < 2.2$ from Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16) as the large-scale structure (LSS) tracer, for which the angular power spectrum $C_\ell^{qq}$ and the redshift-space distortion (RSD) parameter $\beta$ are estimated. By cross correlating with the $\textit{Planck}$ 2018 cosmic microwave background (CMB) lensing map, we detect the angular cross-power spectrum $C_\ell^{\kappa q}$ signal at $12\,\sigma$ significance. Both jackknife resampling and simulations are used to estimate the covariance matrix (CM) of $E_G$ at $5$ bins covering different scales, with the later preferred for its better constraints on the covariances. We find $E_G$ estimates agree with the GR prediction at $1\,\sigma$ level over all these scales. With the CM estimated with $300$ simulations, we report a best-fit scale-averaged estimate of $E_G(\bar{z})=0.30\pm 0.05$, which is in line with the GR prediction $E_G^{\rm GR}(\bar{z})=0.33$ with $\textit{Planck}$ 2018 CMB+BAO matter density fraction $\Omega_{\rm m}=0.31$. The statistical errors of $E_G$ with future LSS surveys at similar redshifts will be reduced by an order of magnitude, which makes it possible to constrain modified gravity models., Comment: 16 pages, 14 figures; references added, matches version accepted by MNRAS

We present a void clustering analysis in configuration-space using the completed Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 samples. These samples consist of Luminous Red Galaxies (LRG) combined with the high redshift tail of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) DR12 CMASS galaxies (called as LRG+CMASS sample), Emission Line Galaxies (ELG) and quasars (QSO). We build void catalogues from the three eBOSS DR16 samples using a ZOBOV-based algorithm, providing 2,814 voids, 1,801 voids and 4,347 voids in the LRG+CMASS, ELG and QSO samples, respectively, spanning the redshift range $0.6

We perform a joint BAO and RSD analysis using the eBOSS DR16 LRG and ELG samples in the redshift range of $z\in[0.6,1.1]$, and detect a RSD signal from the cross power spectrum at a $\sim4\sigma$ confidence level, i.e., $f\sigma_8=0.317\pm0.080$ at $z_{\rm eff}=0.77$. Based on the chained power spectrum, which is a new development in this work to mitigate the angular systematics, we measurement the BAO distances and growth rate simultaneously at two effective redshifts, namely, $D_{\rm M}/r_{\rm d} \ (z=0.70)=17.96\pm0.51, \ D_{\rm H}/r_{\rm d} \ (z=0.70)=21.22\pm1.20, \ f\sigma_8 \ (z=0.70) =0.43\pm0.05$, and $D_{\rm M}/r_{\rm d} \ (z=0.845)=18.90\pm0.78, \ D_{\rm H}/r_{\rm d} \ (z=0.845)=20.91\pm2.86, \ f\sigma_8 \ (z=0.845) =0.30\pm0.08$. Combined with BAO measurements including those from the eBOSS DR16 QSO and Lyman-$\alpha$ sample, our measurement has raised the significance level of a nonzero $\Omega_{\rm \Lambda}$ to $\sim11\sigma$., Comment: 21 pages, 13 figures; version accepted to MNRAS. The data product of this work is publicly available at https://github.com/icosmology/eBOSS_DR16_LRGxELG and at https://www.sdss.org/science/final-bao-and-rsd-measurements/

We present the anisotropic clustering of emission line galaxies (ELGs) from the Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16). Our sample is composed of 173,736 ELGs covering an area of 1170 deg$^2$ over the redshift range $0.6 \leq z \leq 1.1$. We use the Convolution Lagrangian Perturbation Theory in addition to the Gaussian Streaming Redshift-Space Distortions to model the Legendre multipoles of the anisotropic correlation function. We show that the eBOSS ELG correlation function measurement is affected by the contribution of a radial integral constraint that needs to be modelled to avoid biased results. To mitigate the effect from unknown angular systematics, we adopt a modified correlation function estimator that cancels out the angular modes from the clustering. At the effective redshift, $z_{\rm eff}=0.85$, including statistical and systematical uncertainties, we measure the linear growth rate of structure $f\sigma_8(z_{\rm eff}) = 0.35\pm0.10$, the Hubble distance $D_H(z_{\rm eff})/r_{\rm drag} = 19.1^{+1.9}_{-2.1}$ and the comoving angular diameter distance $D_M(z_{\rm eff})/r_{\rm drag} = 19.9\pm1.0$. These results are in agreement with the Fourier space analysis, leading to consensus values of: $f\sigma_8(z_{\rm eff}) = 0.315\pm0.095$, $D_H(z_{\rm eff})/r_{\rm drag} = 19.6^{+2.2}_{-2.1}$ and $D_M(z_{\rm eff})/r_{\rm drag} = 19.5\pm1.0$, consistent with $\Lambda$CDM model predictions with Planck parameters., 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/

Cosmological growth can be measured in the redshift space clustering of galaxies targeted by spectroscopic surveys. Accurate prediction of clustering of galaxies will require understanding galaxy physics which is a very hard and highly non-linear problem. Approximate models of redshift space distortion (RSD) take a perturbative approach to solve the evolution of dark matter and galaxies in the universe. In this paper we focus on eBOSS emission line galaxies (ELGs) which live in intermediate mass haloes. We create a series of mock catalogues using haloes from the Multidark and {\sc Outer Rim} dark matter only N-body simulations. Our mock catalogues include various effects inspired by baryonic physics such as assembly bias and the characteristics of satellite galaxies kinematics, dynamics and statistics deviating from dark matter particles. We analyse these mocks using the TNS RSD model in Fourier space and the CLPT in configuration space. We conclude that these two RSD models provide an unbiased measurement of redshift space distortion within the statistical error of our mocks. We obtain the conservative theoretical systematic uncertainty of $3.3\%$, $1.8\%$ and $1.5\%$ in $f\sigma_8$, $\alpha_{\parallel}$ and $\alpha_{\bot}$ respectively for the TNS and CLPT models. We note that the estimated theoretical systematic error is an order of magnitude smaller than the statistical error of the eBOSS ELG sample and hence are negligible for the purpose of the current eBOSS ELG analysis., Comment: 21 pages, 7 figures and 9 tables, A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found at https://www.sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found at https://www.sdss.org/science/cosmology-results-from-eboss/ . Final published version

We produce 1000 realizations of synthetic clustering catalogues for each type of the tracers used for the baryon acoustic oscillation and redshift space distortion analysis of the Sloan Digital Sky Surveys-IV extended Baryon Oscillation Spectroscopic Survey final data release (eBOSS DR16), covering the redshift range from 0.6 to 2.2, to provide reliable estimates of covariance matrices and test the robustness of the analysis pipeline with respect to observational systematics. By extending the Zel'dovich approximation density field with an effective tracer bias model calibrated with the clustering measurements from the observational data, we accurately reproduce the two- and three-point clustering statistics of the eBOSS DR16 tracers, including their cross-correlations in redshift space with very low computational costs. In addition, we include the gravitational evolution of structures and sample selection biases at different redshifts, as well as various photometric and spectroscopic systematic effects. The agreements on the auto-clustering statistics between the data and mocks are generally within 1 $\sigma$ variances inferred from the mocks, for scales down to a few $h^{-1}\,{\rm Mpc}$ in configuration space, and up to $0.3\,h\,{\rm Mpc}^{-1}$ in Fourier space. For the cross correlations between different tracers, the same level of consistency presents in configuration space, while there are only discrepancies in Fourier space for scales above $0.15\,h\,{\rm Mpc}^{-1}$. The accurate reproduction of the data clustering statistics permits reliable covariances for multi-tracer analysis., Comment: 25 pages, 22 figures; updated to match the version accepted for publication at MNRAS

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 the cosmological analysis of the configuration-space anisotropic clustering in the completed Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 galaxy sample. This sample consists of luminous red galaxies (LRGs) spanning the redshift range $0.6 < z < 1$, at an effective redshift of $z_{\rm eff}=0.698$. It combines 174 816 eBOSS LRGs and 202 642 BOSS CMASS galaxies. We extract and model the baryon acoustic oscillations (BAO) and redshift-space distortions (RSD) features from the galaxy two-point correlation function to infer geometrical and dynamical cosmological constraints. The adopted methodology is extensively tested on a set of realistic simulations. The correlations between the inferred parameters from the BAO and full-shape correlation function analyses are estimated. This allows us to derive joint constraints on the three cosmological parameter combinations: $D_M(z)/r_d$, $D_H(z)/r_d$ and $f\sigma_8(z)$, where $D_M$ is the comoving angular diameter distance, $D_H$ is Hubble distance, $r_d$ is the comoving BAO scale, $f$ is the linear growth rate of structure, and $\sigma_8$ is the amplitude of linear matter perturbations. After combining the results with those from the parallel power spectrum analysis of Gil-Marin et al. 2020, we obtain the constraints: $D_M/r_d = 17.65 \pm 0.30$, $D_H/r_d = 19.77 \pm 0.47$, $f\sigma_8 = 0.473 \pm 0.044$. These measurements are consistent with a flat $\Lambda$CDM model with standard gravity., Comment: 29 pages, 21 figures. Matches accepted version after minor modifications. 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 perform a multi-tracer analysis using the complete Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 luminous red galaxy (LRG) and the DR16 emission line galaxy (ELG) samples in the configuration space, and successfully detect a cross correlation between the two samples, and find the growth rate to be $f\sigma_8=0.342 \pm 0.085$ ($\sim25$ per cent accuracy) from the cross sample alone. We perform a joint measurement of the baryonic acoustic oscillation (BAO) and redshift space distortion (RSD) parameters at a single effective redshift of $z_{\rm eff}= 0.77$, using the auto- and cross-correlation functions of the LRG and ELG samples, and find that the comoving angular diameter distance $D_M(z_{\rm eff})/r_d = 18.85\pm 0.38$, the Hubble distance $D_H(z_{\rm eff})/r_d = 19.64 \pm 0.57$, and $f\sigma_8(z_{\rm eff}) = 0.432 \pm 0.038$, which is consistent with a $\Lambda$CDM model at $68\%$ CL. Compared to the single-tracer analysis on the LRG sample, the Figure of Merit (FoM) of $\alpha_{\perp}, \alpha_{||}$ and $f\sigma_8$ is improved by a factor of $1.11$ in our multi-tracer analysis, and in particular, the statistical uncertainty of $f\sigma_8$ is reduced by $11.6 \%$., Comment: 15 pages, 13 figures, 2 tables; MNRAS accepted; The BAO and RSD measurements and the covariance matrix are made available at https://github.com/ytcosmo/MultiTracerBAORSD/

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/