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

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

This document summarizes the science cases related to cosmology studies with the MaunaKea Spectroscopic Explorer (MSE), a highly-multiplexed (4332 fibers), wide FOV (1.5 sq deg), large aperture (11.25 m in diameter), optical/NIR (360nm to 1300nm) facility. The MSE High-z Cosmology Survey is designed to probe a large volume of the Universe with a galaxy density sufficient to measure the extremely-large-scale density fluctuations required to explore primordial non-Gaussianity and therefore inflation. We expect a measurement of the local parameter $f_{NL}$ to a precision $\sigma(f_{NL}) = 1.8$. Combining the MSE High-z Cosmology Survey data with data from a next generation CMB stage 4 experiment and existing DESI data will provide the first $5\sigma$ confirmation of the neutrino mass hierarchy from astronomical observations. In addition, the Baryonic Acoustic Oscillations (BAO) observed within the sample will provide measurements of the distance-redshift relationship in six different redshift bins between $z=1.6$ and 4.0, each with an accuracy of $\sim0.6\%$. The simultaneous measurements of Redshift Space Distortions (RSD) constrain the amplitude of the fluctuations, at a level ranging from $1.9\%$ to $3.6\%$. The proposed survey covers 10,000 ${\rm deg}^2$, measuring redshifts for three classes of target objects: Emission Line Galaxies (ELGs) with $1.6<z<2.4$, Lyman Break Galaxies (LBGs) with $2.4<z<4.0$, and quasars $2.1<z<3.5$. The ELGs and LBGs will be used as direct tracers of the underlying density field, while the Lyman-$\alpha$ forests in the quasar spectra will be utilized to probe structure. Exposures of duration 1,800sec will guarantee a redshift determination efficiency of $90\%$ for ELGS and at least $50\%$ for LBGs. The survey will represent 100 nights per year for a 5-year MSE program. Finally, three ideas for additional projects of cosmological interest are proposed.

The VIPERS galaxy survey has measured the clustering of $0.5<z<1.2$ galaxies, enabling a number of measurements of galaxy properties and cosmological redshift-space distortions (RSD). Because the measurements were made using one-pass of the VIMOS instrument on the Very Large Telescope (VLT), the galaxies observed only represent approximately 47\% of the parent target sample, with a distribution imprinted with the pattern of the VIMOS slitmask. Correcting for the effect on clustering has previously been achieved using an approximate approach developed using mock catalogues. Pairwise inverse probability (PIP) weighting has recently been proposed by Bianchi & Percival to correct for missing galaxies, and we apply it to mock VIPERS catalogues to show that it accurately corrects the clustering for the VIMOS effects, matching the clustering measured from the observed sample to that of the parent. We then apply PIP-weighting to the VIPERS data, and fit the resulting monopole and quadrupole moments of the galaxy two-point correlation function with respect to the line-of-sight, making measurements of RSD. The results are close to previous measurements, showing that the previous approximate methods used by the VIPERS team are sufficient given the errors obtained on the RSD parameter.

Galaxy redshift surveys are one of the pillars of the current standard cosmological model and remain a key tool in the experimental effort to understand the origin of cosmic acceleration. To this end, the next generation of surveys aim at achieving sub-percent precision in the measurement of the equation of state of dark energy $w(z)$ and the growth rate of structure $f(z)$. This however requires comparable control over systematic errors, stressing the need for improved modelling methods. In this contribution we review at the introductory level some highlights of the work done in this direction by the {\it Darklight} project. Supported by an ERC Advanced Grant, {\it Darklight} developed novel techniques for clustering analysis, which were tested through numerical simulations before being finally applied to galaxy data as in particular those of the recently completed VIPERS redshift survey. We focus in particular on: (a) advances on estimating the growth rate of structure from redshift-space distortions; (b) parameter estimation through global Bayesian reconstruction of the density field from survey data; (c) impact of massive neutrinos on large-scale structure measurements. Overall, {\it Darklight} has contributed to paving the way for forthcoming high-precision experiments, such as {\it Euclid}, the next ESA cosmological mission.

We use the VIPERS final data release to investigate the performance of colour-selected populations of galaxies as tracers of linear large-scale motions. We empirically select volume-limited samples of blue and red galaxies as to minimise the systematic error on the estimate of the growth rate $f\sigma_8$ from the anisotropy of the two-point correlation function. To this end, rather than rigidly splitting the sample into two colour classes we define the red/blue fractional contribution of each object through a weight based on the $(U-V)$ colour distribution. Using mock surveys that are designed to reproduce the observed properties of VIPERS galaxies, we find the systematic error in recovering the fiducial value of $f\sigma_8$ to be minimized when using a volume-limited sample of luminous blue galaxies. We model non-linear corrections via the Scoccimarro extension of the Kaiser model, finding systematic errors on $f\sigma_8$ of below $1-2\%$, using scales as small as 5 $h^{-1}\mathrm{Mpc}$. We interpret this result as indicating that selection of luminous blue galaxies maximises the fraction that are central objects in their dark matter haloes; this in turn minimises the contribution to the measured $\xi(r_p,\pi)$ from the 1-halo term, which is dominated by non-linear motions. The gain is inferior if one uses the full magnitude-limited sample of blue objects, consistent with the presence of a significant fraction of blue, fainter satellites dominated by non-streaming, orbital velocities. We measure a value of $f\sigma_8=0.45 \pm 0.11$ over the single redshift range $0.6\le z\le 1.0$, corresponding to an effective redshift for the blue galaxies $\left<z\right>=0.85$. Including in the likelihood the potential extra information contained in the blue-red galaxy cross-correlation function does not lead to an appreciable improvement in the error bars, while it increases the systematic error.

We present measurements of the growth rate of cosmological structure from the modelling of the anisotropic galaxy clustering measured in the final data release of the VIPERS survey. The analysis is carried out in configuration space and based on measurements of the first two even multipole moments of the anisotropic galaxy auto-correlation function, in two redshift bins spanning the range $0.5 < z < 1.2$. We provide robust and cosmology-independent corrections for the VIPERS angular selection function, allowing recovery of the underlying clustering amplitude at the percent level down to the Mpc scale. We discuss several improvements on the non-linear modelling of redshift-space distortions (RSD) and perform detailed tests of a variety of approaches against a set of realistic VIPERS-like mock realisations. This includes using novel fitting functions to describe the velocity divergence and density power spectra $P_{\theta\theta}$ and $P_{\delta\theta}$ that appear in RSD models. These tests show that we are able to measure the growth rate with negligible bias down to separations of $5h^{-1}Mpc$. Interestingly, the application to real data shows a weaker sensitivity to the details of non-linear RSD corrections compared to mock results. We obtain consistent values for the growth rate times the matter power spectrum normalisation parameter of $f\sigma_8=0.55\pm 0.12$ and $0.40\pm0.11$ at effective redshifts of $z = 0.6$ and $z=0.86$ respectively. These results are in agreement with standard cosmology predictions assuming Einstein gravity in a $\Lambda \rm{CDM}$ background.

We present the latest version of Pinocchio, a code that generates catalogues of DM haloes in an approximate but fast way with respect to an N-body simulation. This code version extends the computation of particle and halo displacements up to 3rd-order Lagrangian Perturbation Theory (LPT), in contrast with previous versions that used Zeldovich approximation (ZA). We run Pinocchio on the same initial configuration of a reference N-body simulation, so that the comparison extends to the object-by-object level. We consider haloes at redshifts 0 and 1, using different LPT orders either for halo construction - where displacements are needed to decide particle accretion onto a halo or halo merging - or to compute halo final positions. We compare the clustering properties of Pinocchio haloes with those from the simulation by computing the power spectrum and 2-point correlation function (2PCF) in real and redshift space (monopole and quadrupole), the bispectrum and the phase difference of halo distributions. We find that 2LPT and 3LPT give noticeable improvement. 3LPT provides the best agreement with N-body when it is used to displace haloes, while 2LPT gives better results for constructing haloes. At the highest orders, linear bias is typically recovered at a few per cent level. In Fourier space and using 3LPT for halo displacements, the halo power spectrum is recovered to within 10 per cent up to $k_{max}\sim0.5\ h/$Mpc. The results presented in this paper have interesting implications for the generation of large ensemble of mock surveys aimed at accurately compute covariance matrices for clustering statistics.

We investigate the use of the cross-correlation between galaxies and galaxy groups to measure redshift-space distortions (RSD) and thus probe the growth rate of cosmological structure. This is compared to the classical approach based on using galaxy auto-correlation. We make use of realistic simulated galaxy catalogues that have been constructed by populating simulated dark matter haloes with galaxies through halo occupation prescriptions. We adapt the classical RSD dispersion model to the case of the group-galaxy cross-correlation function and estimate the RSD parameter $\beta$ by fitting both the full anisotropic correlation function $\xi(r_p,\pi)$ and its multipole moments. In addition, we define a modified version of the latter statistics by truncating the multipole moments to exclude strongly non-linear distortions at small transverse scales. We fit these three observable quantities in our set of simulated galaxy catalogues and estimate statistical and systematic errors on $\beta$ for the case of galaxy-galaxy, group-group, and group-galaxy correlation functions. When ignoring off-diagonal elements of the covariance matrix in the fitting, the truncated multipole moments of the group-galaxy cross-correlation function provide the most accurate estimate, with systematic errors below 3% when fitting transverse scales larger than $10Mpc/h$. Including the full data covariance enlarges statistical errors but keep unchanged the level of systematic error. Although statistical errors are generally larger for groups, the use of group-galaxy cross-correlation can potentially allow the reduction of systematics while using simple linear or dispersion models.