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

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey.

We present catalogues of cosmic voids identified in the distribution of Luminous Red Galaxies (LRGs) and Quasi Stellar Objects (QSOs) in the fourteenth data release (DR14) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We perform a multivariate analysis to assess the level of contamination in these catalogues by spurious Poisson underdensities. We find that the LRG void catalogue is largely free from contamination but that the QSO catalogue may be heavily contaminated. We analyse the multipoles of the void-galaxy cross-correlation function in these catalogues to obtain constraints on the growth rate of structure around voids. We find a value of $\beta(z=0.703)=0.58^{+0.33}_{-0.28}$ for the LRG voids and $\beta(z=1.53)=0.15^{+0.13}_{-0.12}$ for the QSO voids.

We use the statistics of the VIPERS survey to investigate the relation between the surface mean stellar mass density Sigma=Mstar/(2*pi*Re^2) of massive passive galaxies (MPGs, Mstar>10^11 Msun) and their environment in the redshift range 0.5<z<0.8. Passive galaxies were selected on the basis of their NUVrK colors (~900 objects), and the environment was defined as the galaxy density contrast, delta, using the fifth nearest-neighbor approach. The analysis of Sigma vs. delta was carried out in two stellar mass bins. In galaxies with Mstar<2*10^11 Msun, no correlation between Sigma and delta is observed. This implies that the accretion of satellite galaxies, which is more frequent in denser environments and efficient in reducing the galaxy Sigma, is not relevant in the formation and evolution of these systems. Conversely, in galaxies with Mstar>2*10^11 Msun, we find an excess of MPGs with low Sigma and a deficit of high-Sigma MPGs in the densest regions wrt other environments. We interpret this result as due to the migration of some high-Sigma MPGs (<1% of the total population of MPGs) into low-Sigma MPGs, probably through mergers or cannibalism of small satellites. In summary, our results imply that the accretion of satellite galaxies has a marginal role in the mass-assembly history of most MPGs. We have previously found that the number density of VIPERS massive star-forming galaxies (MSFGs) declines rapidily from z=0.8 to z=0.5, which mirrors the rapid increase in the number density of MPGs. This indicates that the MSFGs at z>0.8 migrate to the MPG population. Here, we investigate the Sigma-delta relation of MSFGs at z>0.8 and find that it is consistent within 1 sigma with that of low-Sigma MPGs at z<0.8. Thus, the results of this and our previous paper show that MSFGs at z>0.8 are consistent in terms of number and environment with being the progenitors of low-Sigma MPGs at z<0.8.

Do void statistics contain information beyond the tracer 2-point correlation function? Yes! As we vary the sum of the neutrino masses, we find void statistics contain information absent when using just tracer 2-point statistics. Massive neutrinos uniquely affect cosmic voids. We explore their impact on void clustering using both the DEMNUni and MassiveNuS simulations. For voids, neutrino effects depend on the observed void tracers. As the neutrino mass increases, the number of small voids traced by cold dark matter particles increases and the number of large voids decreases. Surprisingly, when massive, highly biased, halos are used as tracers, we find the opposite effect. The scale at which voids cluster, as well as the void correlation, is similarly sensitive to the sum of neutrino masses and the tracers. This scale dependent trend is not due to simulation volume or halo density. The interplay of these signatures in the void abundance and clustering leaves a distinct fingerprint that could be detected with observations and potentially help break degeneracies between different cosmological parameters. This paper paves the way to exploit cosmic voids in future surveys to constrain the mass of neutrinos.

In this paper, we study the normalised characteristic scale of transition to cosmic homogeneity, $\mathcal{R}_H/d_V$, as a cosmological probe. We use a compilation of the SDSS galaxy samples, comprising more than $10^6$ galaxies in the redshift range $0.17 \leq z \leq 2.2$ within the largest comoving volume to date, $\sim 8 h^{-3}Gpc^3$. We show that these samples can be described by a single bias model as a function of redshift. By combining our measurements with prior Cosmic Microwave Background and Lensing information from the Planck satellite, we constrain the total matter density ratio of the universe $\Omega_m = 0.340 \pm 0.029$ and the Dark Energy density ratio $\Omega_{\Lambda} = 0.668 \pm 0.023$. Our results are compatible with a flat $\Lambda$CDM model. These results show the complementarity of the normalised homogeneity scale with other cosmological probes and open new roads to cosmometry.

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

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

The aim of this study is to distinguish genuine cosmic voids, found in a galaxy catalog by the void finder ZOBOV-VIDE, from under-dense regions in a Poisson distribution of objects. For this purpose, we perform two multivariate analyses using the following physical void characteristics: volume, redshift, density contrast, minimum density, contrast significance and number of member galaxies of the void. The multivariate analyses are trained on a catalog of voids obtained from a random Poisson distribution of points, used as background, and a catalog of voids identified in a mock galaxy catalog, used as signal. The classifications are then applied to voids extracted from the Data Release 12 sample of Luminous Red Galaxies in the redshift range 0.45 $\leq$ z $\leq$ 0.7 from the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (SDSS BOSS DR12 CMASS). Our results show that the resulting void catalog is nearly free of contamination by Poisson noise. We also study the effect of tracer sparsity and bias on the classification efficiencies.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

In this paper, we study the characteristic scale of transition to cosmic homogeneity of the universe, $\mathcal{R}_H$, as a standard ruler, to constrain cosmological parameters on mock galaxy catalogues. We use mock galaxy catalogues that simulate the CMASS galaxy sample of the BOSS survey in the redshift range $0.43 \leq z \leq 0.7$. In each redshift bin we obtain the homogeneity scale, defined as the scale at which the universe becomes homogeneous to $1\%$, i.e. $D_2(\mathcal{R}_H) = 2.97$. With a simple Fisher analysis, we find that the performance of measuring the cosmological parameters with either the position of the BAO peak or the homogeneity scale is comparable. We show that $\mathcal{R}_H$ has a dependence on the galaxy bias. If the accuracy and precision of this bias is achieved to $1\%$, as expected for future surveys, then $\mathcal{R}_H$ is a competitive standard ruler.

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.

Using an unconventional single line diagnostic that unambiguously identifies AGNs in composite galaxies we report statistical differences in the properties (stellar age, [OII] luminosity, colour) between active and inactive galaxies at 0.62<z<1.2 extracted from the VIMOS Public Extragalactic Redshift Survey (VIPERS). The nuclear activity is probed by the high-ionization [NeV] emission line and along with their parent samples, the galaxies are properly selected according to their stellar mass, redshift, and colour distributions. We report younger underlying stellar ages and higher [OII] luminosities of active galaxies in the green valley and in the blue cloud compared to control samples. We observe higher fractions of green galaxies hosting AGN activity at progressively bluer (r-K) colours. Depending on the location of the host galaxy in the NUVrK colour diagram we find higher AGN fractions in massive blue galaxies and in the least massive red galaxies, in agreement with the picture that black holes vary their properties when hosted in either star-forming or passive galaxies. Exactly where the fast quenching processes are expected to play a role, we identify a novel class of active galaxies in the blue cloud with signatures typical for a suddenly suppression of their star formation activity after a burst happening in the recent past. Their optical spectra resemble those of post-starburst galaxies, that would never be identified in a spectroscopic search using classical selection techniques. Broadly, these active galaxies selected on the [NeV] line are not commonly represented in shallow X-ray, mid-IR, or classical line diagnostics. If we consider that our results are limited by the shallow observational limits and rapid AGN variability, the impact of AGN feedback on galaxy formation and evolution may represent an important channel of fast-transiting galaxies moving to the red sequence.

[Abridged] We use the final data of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to investigate the effect of environment on the evolution of galaxies between $z=0.5$ and $z=0.9$. We characterise local environment in terms of the density contrast smoothed over a cylindrical kernel, the scale of which is defined by the distance to the $5^{th}$ nearest neighbour. We find that more massive galaxies tend to reside in higher-density environments over the full redshift range explored. Defining star-forming and passive galaxies through their (NUV$-r$) vs ($r-K$) colours, we then quantify the fraction of star-forming over passive galaxies, $f_{\rm ap}$, as a function of environment at fixed stellar mass. $f_{\rm ap}$ is higher in low-density regions for galaxies with masses ranging from $\log(\mathcal{M}/\mathcal{M}_\odot)=10.38$ (the lowest value explored) to at least $\log(\mathcal{M}/\mathcal{M}_\odot)\sim11.3$, although with decreasing significance going from lower to higher masses. This is the first time that environmental effects on high-mass galaxies are clearly detected at redshifts as high as $z\sim0.9$. We compared these results to VIPERS-like galaxy mock catalogues based on the galaxy formation model of De Lucia & Blaizot. The model correctly reproduces $f_{\rm ap}$ in low-density environments, but underpredicts it at high densities. The discrepancy is particularly strong for the lowest-mass bins. We find that this discrepancy is driven by an excess of low-mass passive satellite galaxies in the model. Looking at the accretion history of these model galaxies, i.e. the times when they become satellites, a better (yet not perfect) agreement with observations can be obtained in high density regions by assuming either that a not-negligible fraction of satellites is destroyed, or that their quenching time-scale is longer than $\sim 2$ Gyr.

We present the first quantitative detection of large-scale filamentary structure at $z \simeq 0.7$ in the large cosmological volume probed by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We use simulations to show the capability of VIPERS to recover robust topological features in the galaxy distribution, in particular the filamentary network. We then investigate how galaxies with different stellar masses and stellar activities are distributed around the filaments and find a significant segregation, with the most massive or quiescent galaxies being closer to the filament axis than less massive or active galaxies. The signal persists even after down-weighting the contribution of peak regions. Our results suggest that massive and quiescent galaxies assemble their stellar mass through successive mergers during their migration along filaments towards the nodes of the cosmic web. On the other hand, low-mass star-forming galaxies prefer the outer edge of filaments, a vorticity rich region dominated by smooth accretion, as predicted by the recent spin alignment theory. This emphasizes the role of large scale cosmic flows in shaping galaxy properties.

We identified voids in the completed VIMOS Public Extragalactic Redshift Survey (VIPERS), using an algorithm based on searching for empty spheres. We measured the cross-correlation between the centres of voids and the complete galaxy catalogue. The cross-correlation function exhibits a clear anisotropy in both VIPERS fields (W1 and W4), which is characteristic of linear redshift space distortions. By measuring the projected cross-correlation and then deprojecting it we are able to estimate the undistorted cross-correlation function. We propose that given a sufficiently well measured cross-correlation function one should be able to measure the linear growth rate of structure by applying a simple linear Gaussian streaming model for the redshift space distortions (RSD). Our study of voids in 306 mock galaxy catalogues mimicking the VIPERS fields would suggest that VIPERS is capable of measuring $\beta$ with an error of around $25\%$. Applying our method to the VIPERS data, we find a value for the redshift space distortion parameter, $\beta = 0.423^{+0.104}_{-0.108}$, which given the bias of the galaxy population we use gives a linear growth rate of $f\sigma_8 = 0.296^{+0.075}_{-0.078}$ at $z = 0.727$. These results are consistent with values observed in parallel VIPERS analysis using standard techniques.

We use the final data from the VIPERS redshift survey to extract an unparalleled sample of more than 2000 massive M > 10^11 M_sun passive galaxies (MPGs) at redshift 0.5 < z < 1.0, based on their NUVrK colours. This enables us to investigate how the population of these objects was built up over cosmic time. We find that the evolution of the number density depends on the galaxy mean surface stellar mass density, Sigma. In particular, dense (Sigma > 2000 M_sun pc^-2) MPGs show a constant comoving number density over this redshift range, whilst this increases by a factor ~ 4 for the least dense objects, defined as having Sigma < 1000 M_sun pc^-2. We estimate stellar ages for the MPG population both fitting the Spectral Energy Distribution (SED) and through the D4000_n index, obtaining results in good agreement. Our findings are consistent with passive ageing of the stellar content of dense MPGs. We show that at any redshift the less dense MPGs are younger than dense ones and that their stellar populations evolve at a slower rate than predicted by passive evolution. This points to a scenario in which the overall population of MPGs was built up over the cosmic time by continuous addition of less dense galaxies: on top of an initial population of dense objects that passively evolves, new, larger, and younger MPGs continuously join the population at later epochs. Finally, we demonstrate that the observed increase in the number density of MPGs is totally accounted for by the observed decrease in the number density of correspondingly massive star forming galaxies (i.e. all the non-passive M > 10^11 M_sun objects). Such systems observed at z ~ 1 in VIPERS, therefore, represent the most plausible progenitors of the subsequent emerging class of larger MPGs.

It is well known that the probability distribution function (PDF) of galaxy density contrast is approximately lognormal; whether the PDF of mass fluctuations derived from weak lensing convergence (kappa_WL) is lognormal is less well established. We derive PDFs of the galaxy and projected matter density distributions via the Counts in Cells (CiC) method. We use maps of galaxies and weak lensing convergence produced from the Dark Energy Survey (DES) Science Verification data over 139 deg^2. We test whether the underlying density contrast is well described by a lognormal distribution for the galaxies, the convergence and their joint PDF. We confirm that the galaxy density contrast distribution is well modeled by a lognormal PDF convolved with Poisson noise at angular scales from 10-40 arcmin (corresponding to physical scales of 3-10 Mpc). We note that as kappa_WL is a weighted sum of the mass fluctuations along the line of sight, its PDF is expected to be only approximately lognormal. We find that the kappa_WL distribution is well modeled by a lognormal PDF convolved with Gaussian shape noise at scales between 10 and 20 arcmin, with a best-fit chi^2/DOF of 1.11 compared to 1.84 for a Gaussian model, corresponding to p-values 0.35 and 0.07 respectively, at a scale of 10 arcmin. Above 20 arcmin a simple Gaussian model is sufficient. The joint PDF is also reasonably fitted by a bivariate lognormal. As a consistency check we compare the variances derived from the lognormal modelling with those directly measured via CiC. Our methods are validated against maps from the MICE Grand Challenge N-body simulation.

Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectrum, the effective power spectrum deviates more significantly from the $\Lambda$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150\% while the deviation of the conventional matter power spectrum is less than 50\%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $\Lambda$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies. In addition, based on our findings, we also discuss several possible methods of making use of the differences between the conventional and effective dark matter power spectra in $f(R)$ gravity to discriminate the theory from the $\Lambda$CDM model.

Cosmic structures leave an imprint on the microwave background radiation through the integrated Sachs-Wolfe effect. We construct a template map of the linear signal using the SDSS-III Baryon Acoustic Oscillation Survey at redshift 0.43 < z < 0.65. We verify the imprint of this map on the Planck CMB temperature map at the 97% confidence level and show consistency with the density-temperature cross-correlation measurement. Using this ISW reconstruction as a template we investigate the presence of ISW sources and further examine the properties of the Granett-Neyrinck-Szapudi supervoid and supercluster catalogue. We characterise the three-dimensional density profiles of these structures for the first time and demonstrate that they are significant structures. Model fits demonstrate that the supervoids are elongated along the line-of-sight and we suggest that this special orientation may be picked out by the void-finding algorithm in photometric redshift space. We measure the mean temperature profiles in Planck maps from public void and cluster catalogues. In an attempt to maximise the stacked ISW signal we construct a new catalogue of super-structures based upon local peaks and troughs of the gravitational potential. However, we do not find a significant correlation between these structures and the CMB temperature.

We introduce the idea of {\it effective} dark matter halo catalog in $f(R)$ gravity, which is built using the {\it effective} density field. Using a suite of high resolution N-body simulations, we find that the dynamical properties of halos, such as the distribution of density, velocity dispersion, specific angular momentum and spin, in the effective catalog of $f(R)$ gravity closely mimic those in the $\Lambda$CDM model. Thus, when using effective halos, an $f(R)$ model can be viewed as a $\Lambda$CDM model. This effective catalog therefore provides a convenient way for studying the baryonic physics, the galaxy halo occupation distribution and even semi-analytical galaxy formation in $f(R)$ cosmologies.

We reexamine the screening mechanism in $f(R)$ gravity using N-body simulations. By explicitly examining the relation between the extra scalar field $\delta f_R$ and the gravitational potential $\phi$ in the perturbed Universe, we find that the relation between these two fields plays an important role in understanding the screening mechanism. We show that the screening mechanism in $f(R)$ gravity depends mainly on the depth of the potential well, and find a useful condition for identifying unscreened halos in simulations. We also discuss the potential application of our results to real galaxy surveys.

The characterisation of cosmic voids gives unique information about the large-scale distribution of galaxies, their evolution and the cosmological model. We identify and characterise cosmic voids in the VIMOS Public Extragalactic Redshift Survey (VIPERS) at redshift 0.55 < z < 0.9. A new void search method is developed based upon the identification of empty spheres that fit between galaxies. The method can be used to characterise the cosmic voids despite the presence of complex survey boundaries and internal gaps. We investigate the impact of systematic observational effects and validate the method against mock catalogues. We measure the void size distribution and the void-galaxy correlation function. We construct a catalogue of voids in VIPERS. The distribution of voids is found to agree well with the distribution of voids found in mock catalogues. The void-galaxy correlation function shows indications of outflow velocity from the voids.

The degeneracy between geometric (Alcock-Paczynski) and dynamic (redshift space) distortions in the pattern of the galaxy distribution has been a long standing problem in the study of the large scale structure of the universe. We examine the possibility of lifting this degeneracy and constraining cosmological parameters by using the Baryon Acoustic Oscillation (BAO) scale as a feature of known physical size, the sound horizon r_s ~ 150 Mpc. We callibrate this scale with the equivalent feature in the Cosmic Microwave Background (CMB). First, we construct a toy model of a power spectrum which includes the BAO as well as geometric and dynamic distortions. By adding a prior onto the sound horizon of ~1% we show, using a Fisher matrix analysis, that error ellipses for line of sight and tangential distortion parameters shrink by a factor of two for a 20(Gpc/h)^3 `DESpec/BigBOSS'-like galaxy survey including shot noise. This improvement is even more marked in smaller surveys. We also carry out a Monte Carlo Nested Sampling analysis on our parameter space. We find that Monte Carlo and Fisher methods can agree reasonably well for surveys with large volume but differ greatly for small volume surveys.

We present equations for the gravitational lensing flexion expected for an elliptical lens mass distribution. These can be reduced to one-dimensional finite integrals, thus saving significant computing time over a full two-dimensional calculation. We estimate constraints on galaxy halo ellipticities for a range of potential future surveys, finding that the constraints from the two different types of flexion are comparable and are up to two orders of magnitude tighter than those from shear. Flexion therefore appears to be a very promising potential tool for constraining the shapes of galaxy haloes from future surveys.