Search

We extend the evolution mapping approach, originally proposed by Sanchez (2022) to describe non-linear matter density fluctuations, to statistics of the cosmic velocity field. This framework classifies cosmological parameters into shape parameters, which determine the shape of the linear matter power spectrum, $P_L(k, z)$, and evolution parameters, which control its amplitude at any given redshift. Evolution mapping leverages the fact that density fluctuations in cosmologies with identical shape parameters but different evolution parameters exhibit remarkably similar non-linear evolutions when expressed as a function of the clustering amplitude. We use a suite of N-body simulations sharing identical shape parameters but spanning a wide range of evolution parameters. Using an efficient method for estimating the volume-weighted velocity field based on the Voronoi tesselation of the simulation particles, we study the non-linear evolution of the power spectra of the velocity divergence, $P_{\theta\theta}(k)$, and its cross-power spectrum with the density field, $P_{\delta\theta}(k)$. By analysing snapshots at redshifts where the linear matter perturbations have the same amplitude, we demonstrate that evolution mapping accurately applies to $P_{\theta\theta}(k)$ and $P_{\delta\theta}(k)$. Deviations at small scales can be modelled in terms of differences in the suppression factor, $g(a) = D(a)/a$, akin to those observed for the density field. Evolution mapping simplifies the description of the cosmological dependence of non-linear density and velocity statistics, streamlining the sampling of large cosmological parameter spaces for the analysis of cosmological observables., Comment: 10 pages, 4 figures, including an appendix. Prepared for submission to MNRAS, comments welcome

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

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

We present the Texas Euclid Survey for Lyman-Alpha (TESLA), a spectroscopic survey in the 10 square degree of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyman-alpha emitters (LAEs) correlate with Lyman-alpha emission to understand the escape of Lyman alpha from galaxies at redshifts 2 -- 3.5. We present an analysis of 43 LAEs performed in the NEP field using early data from the TESLA survey. We use Subaru Hyper Suprime-Cam imaging in the grizy-bands, Spitzer/IRAC channels 1 and 2 from the Hawaii 20 square degree (H20) survey and spectra acquired by the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby-Eberly Telescope. We perform spectral energy distribution (SED) fitting to compute the galaxy properties of 43 LAEs, and study correlations between stellar mass, star formation rate (SFR), and dust, to the Lyman-alpha rest-frame equivalent widths (EW). We uncover marginal (1 sigma significance) correlations between stellar mass and Lyman-alpha EW, and star formation rate (SFR) and Lyman-alpha EW, with a Spearman correlation coefficient of -0.$34_{-.14}^{+.17}$ and -0.$37_{-.14}^{+.16}$ respectively. We show that the Lyman-alpha distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of 150 Angstrom. Once complete the TESLA survey will enable the study of ~ thousands of LAEs to explore correlations between galaxy properties and Lyman-alpha EW. The large sample size will allow the construction of a predictive model for the Lyman-alpha EW as a function of SED-derived galaxy properties, which could be used to improve Lyman-alpha based constraints on reionization.

We analyse the full shape of anisotropic clustering measurements from the extended Baryon Oscillation Spectroscopic survey (eBOSS) quasar sample together with the combined galaxy sample from the Baryon Oscillation Spectroscopic Survey (BOSS). We obtain constraints on the cosmological parameters independent of the Hubble parameter $h$ for the extensions of the $\Lambda$CDM models, focusing on cosmologies with free dark energy equation of state parameter $w$. We combine the clustering constraints with those from the latest CMB data from Planck to obtain joint constraints for these cosmologies for $w$ and the additional extension parameters - its time evolution $w_{\rm{a}}$, the physical curvature density $\omega_{K}$ and the neutrino mass sum $\sum m_{\nu}$. Our joint constraints are consistent with flat $\Lambda$CDM cosmological model within 68\% confidence limits. We demonstrate that the Planck data are able to place tight constraints on the clustering amplitude today, $\sigma_{12}$, in cosmologies with varying $w$ and present the first constraints for the clustering amplitude for such cosmologies, which is found to be slightly higher than the $\Lambda$CDM value. Additionally, we show that when we vary $w$ and allow for non-flat cosmologies and the physical curvature density is used, Planck prefers a curved universe at $4\sigma$ significance, which is $\sim2\sigma$ higher than when using the relative curvature density $\Omega_{\rm{K}}$. Finally, when $w$ is varied freely, clustering provides only a modest improvement (of 0.021 eV) on the upper limit of $\sum m_{\nu}$., Comment: 12 pages, 6 figures, submitted to MNRAS

We present the results of a stellar-population analysis of Lyman-alpha emitting galaxies (LAES) in GOODS-N at 1.9 < z < 3.5 spectroscopically identified by the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We provide a method for connecting emission-line detections from the blind spectroscopic survey to imaging counterparts, a crucial tool needed as HETDEX builds a massive database of ~1 million Lyman-alpha detections. Using photometric data spanning as many as 11 filters covering 0.4-4.5 microns from the Hubble and Spitzer Space Telescopes, we study the objects' global properties and explore which properties impact the strength of Lyman-alpha emission. We measure a median stellar mass of 0.8 (^+2.9_-0.5) x 10^9 Msol and conclude that the physical properties of HETDEX spectroscopically-selected LAEs are comparable to LAEs selected by previous deep narrow band studies. We find that stellar mass and star formation rate correlate strongly with the Lyman-alpha equivalent width. We then use a known sample of z>7 LAEs to perform a proto-study of predicting Lyman-alpha emission from galaxies in the Epoch of Reionization, finding agreement at the 1-sigma level between prediction and observation for the majority of strong emitters., Comment: 29 pages, 22 figures, Accepted to the Astrophysical Journal

In this paper we present COMET, a Gaussian process emulator of the galaxy power spectrum multipoles in redshift-space. The model predictions are based on one-loop perturbation theory and we consider two alternative descriptions of redshift-space distortions: one that performs a full expansion of the real- to redshift-space mapping, as in recent effective field theory models, and another that preserves the non-perturbative impact of small-scale velocities by means of an effective damping function. The outputs of COMET can be obtained at arbitrary redshifts (up to $z \sim 3$), for arbitrary fiducial background cosmologies, and for a large parameter space that covers the shape parameters $\omega_c$, $\omega_b$, and $n_s$, as well as the evolution parameters $h$, $A_s$, $\Omega_K$, $w_0$, and $w_a$. This flexibility does not impair COMET's accuracy, since we exploit an exact degeneracy between the evolution parameters that allows us to train the emulator on a significantly reduced parameter space. While the predictions are sped up by at least two orders of magnitude, validation tests reveal an accuracy of $0.1\,\%$ for the monopole and quadrupole ($0.3\,\%$ for the hexadecapole), or alternatively, better than $0.25\,\sigma$ for all three multipoles in comparison to statistical uncertainties expected for the Euclid survey with a tenfold increase in volume. We show that these differences translate into shifts in mean posterior values that are at most of the same size, meaning that COMET can be used with the same confidence as the exact underlying models. COMET is a publicly available Python package that also provides the tree-level bispectrum multipoles in redshift-space and Gaussian covariance matrices., Comment: 19 pages, 11 figures; for the COMET Python package, see https://gitlab.com/aegge/comet-emu; v2: added discussion on binning and discreteness effects, minor improvements to match published version

We report and analyse the presence of foregrounds in the cosmic microwave background (CMB) radiation associated to extended galactic halos. Using the cross correlation of Planck and WMAP maps and the 2MRS galaxy catalogue, we find that the mean temperature radial profiles around nearby galaxies at $cz\le 4500~\rm{km~s^{-1}}$ show a statistically significant systematic decrease of $\sim 15~\mu \rm{K}$ extending up to several galaxy radii. This effect strongly depends on the galaxy morphological type at scales within several tens of times the galaxy size, becoming nearly independent of galaxy morphology at larger scales. The effect is significantly stronger for the more extended galaxies, with galaxy clustering having a large impact on the results. Our findings indicate the presence of statistically relevant foregrounds in the CMB maps that should be considered in detailed cosmological studies. Besides, we argue that these can be used to explore the intergalactic medium surrounding bright late-type galaxies and allow for diverse astrophysical analyses., Comment: 10 pages, 9 figures, submitted to MNRAS

We present the analysis of the full shape of anisotropic clustering measurement from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) quasar sample together with the combined galaxy sample from the Baryon Oscillation Spectroscopic Survey (BOSS), re-analysed using an updated recipe for the non-linear matter power spectrum and the non-local bias parameters. We obtain constraints for flat $\Lambda$CDM cosmologies, focusing on the cosmological parameters that are independent of the Hubble parameter $h$. Our recovered value for the RMS linear perturbation theory variance as measured on the scale of $12\,{\rm Mpc}$ is $\sigma_{12}=0.805\pm 0.049$, while using the traditional reference scale of $8\,h^{-1}{\rm Mpc}$ gives $\sigma_{8}=0.815\pm 0.044$. We quantify the agreement between our measurements and the latest CMB data from Planck using the suspiciousness metric, and find them to be consistent within $0.64 \pm 0.03\sigma$. Combining our clustering constraints with the $3\times2$pt data sample from the Dark Energy Survey (DES) Year 1 release slightly degrades this agreement to the level of $1.54 \pm 0.08\sigma$, while still showing an overall consistency with Planck. We furthermore study the effect of imposing a Planck - like prior on the parameters that define the shape of the linear matter power spectrum, and find significantly tighter constraints on the parameters that control the evolution of density fluctuations. In particular, the combination of low-redshift data sets prefers a value of the physical dark energy density $\omega_{\rm DE}=0.335 \pm 0.011$, which is 1.7$\sigma$ higher than the one preferred by Planck., Comment: 14 pages, 10 figures, matches the version accepted by MNRAS

We present a new approach to describe statistics of the non-linear matter density field that exploits a degeneracy in the impact of different cosmological parameters on the linear dimensionless matter power spectrum, $\Delta^2_{\rm L}(k)$. We classify all cosmological parameters into two groups, shape parameters, which determine the shape of $\Delta^2_{\rm L}(k)$, and evolution parameters, which only affect its amplitude at any given redshift. With this definition, the time evolution of $\Delta^2_{\rm L}(k)$ in models with identical shape parameters but different evolution parameters can be mapped from one to the other by relabelling the redshifts that correspond to the same clustering amplitude, which we characterize by the linear mass fluctuation in spheres of radius $12\,{\rm Mpc}$, $\sigma_{12}(z)$. We use N-body simulations to show that the same evolution mapping relation gives a good description of the non-linear power spectrum, the halo mass function, or the full density field. The deviations from the exact degeneracy are the result of the different structure formation histories experienced by each model to reach the same clustering amplitude and can be accurately described in terms of differences in the suppression factor $g(a) = D(a)/a$. These relations can be used to drastically reduce the number of parameters required to describe the cosmology dependence of the power spectrum. We show how this can help to speed up the inference of parameter constraints from cosmological observations. We also present a new design of an emulator of the non-linear power spectrum whose predictions can be adapted to an arbitrary choice of evolution parameters and redshift., Comment: 13 pages, 11 figures, replaced to match version accepted by MNRAS

This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above the shot-noise level identified with a spherical void finder. Under this mapping, void properties are affected by three effects: (i) a systematic expansion as a consequence of the distortions induced by galaxy dynamics, (ii) the Alcock-Paczynski volume effect, which manifests as an overall expansion or contraction depending on the fiducial cosmology, and (iii) a systematic off-centring along the line of sight as a consequence of the distortions induced by void dynamics. We found that correlations are also affected by an additional source of distortions: the ellipticity of voids. This is the first time that distortions due to the off-centring and ellipticity effects are detected and quantified. With a simplified test, we verified that the Gaussian streaming model is still robust provided all these effects are taken into account, laying the foundations for improvements in current models in order to obtain unbiased cosmological constraints from spectroscopic surveys. Besides this practical importance, this analysis also encodes key information about the structure and dynamics of the Universe at the largest scales. Furthermore, some of the effects constitute cosmological probes by themselves, as is the case of the void ellipticity., Comment: 14 pages, 9 figures, published in MNRAS, accepted after minor comments

The construction of catalogues of a particular type of galaxy can be complicated by interlopers contaminating the sample. In spectroscopic galaxy surveys this can be due to the misclassification of an emission line; for example in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) low redshift [OII] emitters may make up a few percent of the observed Ly${\alpha}$ emitter (LAE) sample. The presence of contaminants affects the measured correlation functions and power spectra. Previous attempts to deal with this using the cross-correlation function have assumed sources at a fixed redshift, or not modelled evolution within the adopted redshift bins. However, in spectroscopic surveys like HETDEX, where the contamination fraction is likely to be redshift dependent, the observed clustering of misclassified sources will appear to evolve strongly due to projection effects, even if their true clustering does not. We present a practical method for accounting for the presence of contaminants with redshift-dependent contamination fractions and projected clustering. We show using mock catalogues that our method, unlike existing approaches, yields unbiased clustering measurements from the upcoming HETDEX survey in scenarios with redshift-dependent contamination fractions within the redshift bins used. We show our method returns auto-correlation functions with systematic biases much smaller than the statistical noise for samples with at least as high as 7 per cent contamination. We also present and test a method for fitting for the redshift-dependent interloper fraction using the LAE-[OII] galaxy cross-correlation function, which gives less biased results than assuming a single interloper fraction for the whole sample., Comment: This is a pre-copy edited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record Farrow et al, MNRAS, 2021 is available online at: https://academic.oup.com/mnras/advance-article-abstract/doi/10.1093/mnras/stab1986/6322852

We investigate the impact of different assumptions in the modeling of one-loop galaxy bias on the recovery of cosmological parameters, as a follow up of the analysis done in the first paper of the series at fixed cosmology. We use three different synthetic galaxy samples whose clustering properties match the ones of the CMASS and LOWZ catalogues of BOSS and the SDSS Main Galaxy Sample. We investigate the relevance of allowing for either short range non-locality or scale-dependent stochasticity by fitting the real-space galaxy auto power spectrum or the combination of galaxy-galaxy and galaxy-matter power spectrum. From a comparison among the goodness-of-fit ($\chi^2$), unbiasedness of cosmological parameters (FoB), and figure-of-merit (FoM), we find that a four-parameter model (linear, quadratic, cubic non-local bias, and constant shot-noise) with fixed quadratic tidal bias provides a robust modelling choice for the auto power spectrum of the three samples, up to $k_{\rm max}=0.3\,h\,\mathrm{Mpc}^{-1}$ and for an effective volume of $6\,h^{-3}\,\mathrm{Gpc}^3$. Instead, a joint analysis of the two observables fails at larger scales, and a model extension with either higher derivatives or scale-dependent shot-noise is necessary to reach a similar $k_{\rm max}$, with the latter providing the most stable results. These findings are obtained with three, either hybrid or perturbative, prescriptions for the matter power spectrum, \texttt{RESPRESSO}, gRPT and EFT. In all cases, the inclusion of scale-dependent shot-noise increases the range of validity of the model in terms of FoB and $\chi^2$. Interestingly, these model extensions with additional free parameters do not necessarily lead to an increase in the maximally achievable FoM for the cosmological parameters $\left(h,\,\Omega_ch^2,\,A_s\right)$, which are generally consistent to those of the simpler model at smaller $k_{\rm max}$., Comment: 24 pages, 8 figures

We present a joint likelihood analysis of the real-space power spectrum and bispectrum measured from a variety of halo and galaxy mock catalogs. A novel aspect of this work is the inclusion of nonlinear triangle configurations for the bispectrum, made possible by a complete next-to-leading order ("one-loop") description of galaxy bias, as is already common practice for the power spectrum. Based on the goodness-of-fit and the unbiasedness of the parameter posteriors, we accomplish a stringent validation of this model compared to the leading order ("tree-level") bispectrum. Using measurement uncertainties that correspond to an effective survey volume of $6\,(\mathrm{Gpc}/h)^3$, we determine that the one-loop corrections roughly double the applicable range of scales, from $\sim 0.17\,h/\mathrm{Mpc}$ (tree-level) to $\sim 0.3\,h/\mathrm{Mpc}$. This converts into a $1.5 - 2$x improvement on constraints of the linear bias parameter at fixed cosmology, and a $1.5 - 2.4$x shrinkage of uncertainties on the amplitude of fluctuations $A_s$, which clearly demonstrates the benefit of extracting information from nonlinear scales despite having to marginalize over a larger number of bias parameters. Besides, our precise measurements of galaxy bias parameters up to fourth order allow for thorough comparisons to coevolution relations, showing excellent agreement for all contributions generated by the nonlocal action of gravity. Using these relations in the likelihood analysis does not compromise the model validity and is crucial for obtaining the quoted improvements on $A_s$. We also analyzed the impact of higher-derivative and scale-dependent stochastic terms, finding that for a subset of our tracers the former can boost the performance of the tree-level model with constraints on $A_s$ that are only slightly degraded compared to the one-loop model., Comment: 33 pages, 13 figures, supporting Mathematica notebook available at https://doi.org/10.5281/zenodo.4529886, Python package for Delaunay binning correction available at https://github.com/aegge/BispTools

Minkowski functionals (MFs) are a set of statistics that characterise the geometry and topology of the cosmic density field and contain complementary information to the standard two-point analyses. We present MEDUSA, an implementation of an accurate method for estimating the MFs of three-dimensional point distributions. These estimates are inferred from triangulated isodensity surfaces that are constructed from the Delaunay tessellation of the input point sample. Contrary to previous methods, MEDUSA can account for periodic boundary conditions, which is crucial for the analysis of N-body simulations. We validate our code against several test samples with known MFs, including Gaussian random fields with a ${\Lambda}$CDM power spectrum, and find excellent agreement with the theory predictions. We use MEDUSA to measure the MFs of synthetic galaxy catalogues constructed from N-body simulations. Our results show clearly non-Gaussian signatures that arise from the non-linear gravitational evolution of the density field. We find that, although redshift-space distortions significantly change our MFs estimates, their impact is considerably reduced if these measurements are expressed as a function of the volume-filling fraction. We also show that the effect of Alcock-Paczynski (AP) distortions on the MFs can be described by scaling them with different powers of the isotropic AP parameter $q$ defined in terms of the volume-averaged distance $D_{\rm V}(z)$. Thus the MFs estimates by MEDUSA are useful probes of non-linearities in the density field, and the expansion and growth of structure histories of the Universe., Comment: 14 pages, 16 figures, submitted to MNRAS

All estimators of the two-point correlation function are based on a random catalogue, a set of points with no intrinsic clustering following the selection function of a survey. High-accuracy estimates require the use of large random catalogues, which imply a high computational cost. We propose to replace the standard random catalogues by glass-like point distributions or glass catalogues, which are characterized by a power spectrum $P(k)\propto k^4$ and exhibit significantly less power than a Poisson distribution with the same number of points on scales larger than the mean inter-particle separation. We show that these distributions can be obtained by iteratively applying the technique of Zeldovich reconstruction commonly used in studies of baryon acoustic oscillations (BAO). We provide a modified version of the widely used Landy-Szalay estimator of the correlation function adapted to the use of glass catalogues and compare its performance with the results obtained using random samples. Our results show that glass-like samples do not add any bias with respect to the results obtained using Poisson distributions. On scales larger than the mean inter-particle separation of the glass catalogues, the modified estimator leads to a significant reduction of the variance of the Legendre multipoles $\xi_\ell(s)$ with respect to the standard Landy-Szalay results with the same number of points. The size of the glass catalogue required to achieve a given accuracy in the correlation function is significantly smaller than when using random samples. Even considering the small additional cost of constructing the glass catalogues, their use could help to drastically reduce the computational cost of configuration-space clustering analysis of future surveys while maintaining high-accuracy requirements., Comment: 9 pages, 7 figures, minor changes to match version accepted by MNRAS

We present constraints on extensions to the flat $\Lambda$CDM cosmological model by varying the spatial curvature $\Omega_K$, the sum of the neutrino masses $\sum m_\nu$, the dark energy equation of state parameter $w$, and the Hu-Sawicki $f(R)$ gravity $f_{R0}$ parameter. With the combined $3\times2$pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing from the overlap between KiDS-1000, BOSS, and the spectroscopic 2-degree Field Lensing Survey (2dFLenS), we find results that are fully consistent with a flat $\Lambda$CDM model with $\Omega_K=0.011^{+0.054}_{-0.057}$, $\sum m_\nu<1.76$ eV (95% CL), and $w=-0.99^{+0.11}_{-0.13}$. The $f_{R0}$ parameter is unconstrained in our fully non-linear $f(R)$ cosmic shear analysis. Considering three different model selection criteria, we find no clear preference for either the fiducial flat $\Lambda$CDM model or any of the considered extensions. Besides extensions to the flat $\Lambda$CDM parameter space, we also explore restrictions to common subsets of the flat $\Lambda$CDM parameter space by fixing the amplitude of the primordial power spectrum to the Planck best-fit value, as well as adding external data from supernovae and lensing of the CMB. Neither the beyond-$\Lambda$CDM models nor the imposed restrictions explored in this analysis are able to resolve the $\sim 3\sigma$ tension in $S_8$ between the $3\times2$pt constraints and Planck, with the exception of $w$CDM, where the $S_8$ tension is resolved. The tension in the $w$CDM case persists, however, when considering the joint $S_8$-$w$ parameter space. The joint flat $\Lambda$CDM CMB lensing and $3\times2$pt analysis is found to yield tight constraints on $\Omega_{\rm m}=0.307^{+0.008}_{-0.013}$, $\sigma_8=0.769^{+0.022}_{-0.010}$, and $S_8=0.779^{+0.013}_{-0.013}$., Comment: 15 pages, 9 figures, 3 tables, accepted in A&A. This paper concludes the KiDS-1000 series of papers: Heymans, Tr\"oster et al. (arXiv:2007.15632), Asgari et al. (arXiv:2007.15633), Hildebrandt et al. (arXiv:2007.15635), Joachimi et al. (arXiv:2007.01844), and Giblin et al. (arXiv:2007.01845). Data products, likelihoods, and posteriors can be found http://kids.strw.leidenuniv.nl/DR4/lensing.php

We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS and the spectroscopic 2-degree Field Lensing Survey (2dFLenS). This combination of large-scale structure probes breaks the degeneracies between cosmological parameters for individual observables, resulting in a constraint on the structure growth parameter $S_8=\sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.766^{+0.020}_{-0.014}$, that has the same overall precision as that reported by the full-sky cosmic microwave background observations from Planck. The recovered $S_8$ amplitude is low, however, by $8.3 \pm 2.6$ % relative to Planck. This result builds from a series of KiDS-1000 analyses where we validate our methodology with variable depth mock galaxy surveys, our lensing calibration with image simulations and null-tests, and our optical-to-near-infrared redshift calibration with multi-band mock catalogues and a spectroscopic-photometric clustering analysis. The systematic uncertainties identified by these analyses are folded through as nuisance parameters in our cosmological analysis. Inspecting the offset between the marginalised posterior distributions, we find that the $S_8$-difference with Planck is driven by a tension in the matter fluctuation amplitude parameter, $\sigma_8$. We quantify the level of agreement between the CMB and our large-scale structure constraints using a series of different metrics, finding differences with a significance ranging between $\sim\! 3\,\sigma$, when considering the offset in $S_{8}$, and $\sim\! 2\,\sigma$, when considering the full multi-dimensional parameter space., Comment: 25 pages, 11 figures, 5 tables, A&A accepted, including a new appendix on Intrinsic Alignments. The KiDS-1000 data products are available for download at http://kids.strw.leidenuniv.nl/DR4/lensing.php. This data release includes open source software, the shear-photo-z catalogue, the cosmic shear and 3x2pt data vectors and covariances, and posteriors in the form of Multinest chains

Voids are promising cosmological probes. Nevertheless, every cosmological test based on voids must necessarily employ methods to identify them in redshift space. Therefore, redshift-space distortions (RSD) and the Alcock-Paczynski effect (AP) have an impact on the void identification process itself generating distortion patterns in observations. Using a spherical void finder, we developed a statistical and theoretical framework to describe physically the connection between the identification in real and redshift space. We found that redshift-space voids above the shot noise level have a unique real-space counterpart spanning the same region of space, they are systematically bigger and their centres are preferentially shifted along the line of sight. The expansion effect is a by-product of RSD induced by tracer dynamics at scales around the void radius, whereas the off-centring effect constitutes a different class of RSD induced at larger scales by the global dynamics of the whole region containing the void. The volume of voids is also altered by the fiducial cosmology assumed to measure distances, this is the AP change of volume. These three systematics have an impact on cosmological statistics. In this work, we focus on the void size function. We developed a theoretical framework to model these effects and tested it with a numerical simulation, recovering the statistical properties of the abundance of voids in real space. This description depends strongly on cosmology. Hence, we lay the foundations for improvements in current models of the abundance of voids in order to obtain unbiased cosmological constraints from redshift surveys., Comment: 16 pages, 12 figures, accepted in MNRAS

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

We 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 test the regime of validity of one-loop galaxy bias for a wide variety of biased tracers. Our most stringent test asks the bias model to simultaneously match the galaxy-galaxy and galaxy-mass spectrum, using the measured nonlinear matter spectrum from the simulations to test one-loop effects from the bias expansion alone. In addition, we investigate the relevance of short-range nonlocality and halo exclusion through higher-derivative and scale-dependent noise terms, as well as the impact of using co-evolution relations to reduce the number of free fitting parameters. From comparing validity and merit of these assumptions we find that a four-parameter model (linear, quadratic, cubic nonlocal bias, and constant shot noise) with fixed quadratic tidal bias provides a robust modeling choice for the auto power spectrum of the less massive halos in our set of samples and their galaxy populations (up to $k_{\mathrm{max}} = 0.35\,h/\mathrm{Mpc}$ for a sample volume of $6\,(\mathrm{Gpc}/h)^3$). For the more biased tracers it is most beneficial to include scale-dependent noise. This is also the preferred option when considering combinations of the auto and cross power spectrum, which might be relevant in joint studies of galaxy clustering and weak lensing. We also test the use of perturbation theory to account for matter loops through gRPT, EFT and the hybrid approach RESPRESSO. While all these have similar performance, we find the latter to be the best in terms of validity and recovered mean posterior values, in accordance with it being based partially on simulations., Comment: 30 pages, 12 figures; v2: includes gRPT expressions and some minor edits, matches published version

It is common to express cosmological measurements in units of $h^{-1}{\rm Mpc}$. Here, we review some of the complications that originate from this practice. A crucial problem caused by these units is related to the normalization of the matter power spectrum, which is commonly characterized in terms of the linear-theory rms mass fluctuation in spheres of radius $8\,h^{-1}{\rm Mpc}$, $\sigma_8$. This parameter does not correctly capture the impact of $h$ on the amplitude of density fluctuations. We show that the use of $\sigma_8$ has caused critical misconceptions for both the so-called $\sigma_8$ tension regarding the consistency between low-redshift probes and cosmic microwave background data, and the way in which growth-rate estimates inferred from redshift-space distortions are commonly expressed. We propose to abandon the use of $h^{-1}{\rm Mpc}$ units in cosmology and to characterize the amplitude of the matter power spectrum in terms of $\sigma_{12}$, defined as the mass fluctuation in spheres of radius $12\,{\rm Mpc}$, whose value is similar to the standard $\sigma_8$ for $h\sim 0.67$., Comment: 6 pages, 4 figures. Minor changes. Matches version published in PRD

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17)., Comment: DR16 release: Monday Dec 9th 2019. This is the alphabetical order SDSS-IV collaboration data release paper. 25 pages, 6 figures, accepted by ApJS on 11th May 2020. Minor changes clarify or improve text and figures relative to v1

We reanalyse the anisotropic galaxy clustering measurement from the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating that using the full shape information provides cosmological constraints that are comparable to other low-redshift probes. We find $\Omega_\mathrm{m} = 0.317^{+0.015}_{-0.019}$, $\sigma_8 = 0.710\pm 0.049$, and $h = 0.704\pm 0.024$ for flat $\Lambda$CDM cosmologies using uninformative priors on $\Omega_\mathrm{c}h^2$, $100\theta_\mathrm{MC}$, $\ln 10^{10} A_{s}$, and $n_{s}$, and a prior on $\Omega_\mathrm{b}h^2$ that is much wider than current constraints. We quantify the agreement between the Planck 2018 constraints from the cosmic microwave background and BOSS, finding the two data sets to be consistent within a flat $\Lambda$CDM cosmology using the Bayes factor as well as the prior-insensitive suspiciousness statistic. Combining two low-redshift probes, we jointly analyse the clustering of BOSS galaxies with weak lensing measurements from the Kilo-Degree Survey (KV450). The combination of BOSS and KV450 improves the measurement by up to 45%, constraining $\sigma_8 = 0.702\pm 0.029$ and $S_8 = \sigma_8\sqrt{\Omega_\mathrm{m}/0.3} = 0.728\pm 0.026$. Over the full 5D parameter space, the odds in favour of single cosmology describing galaxy clustering, lensing, and the cosmic microwave background are $7\pm2$. The suspiciousness statistic signals a $2.1\pm0.3\sigma$ tension between the combined low-redshift probes and measurements from the cosmic microwave background., Comment: 9 pages, 9 figures. Accepted in A&A

The forthcoming generation of galaxy redshift surveys will sample the large-scale structure of the Universe over unprecedented volumes with high-density tracers. This advancement will make robust measurements of three-point clustering statistics possible. In preparation for this improvement, we investigate how several methodological choices can influence inferences based on the bispectrum about galaxy bias and shot noise. We first measure the real-space bispectrum of dark-matter haloes extracted from 298 N-body simulations covering a volume of approximately $1000 h^{-3} \mathrm{Gpc}^3$. We then fit a series of theoretical models based on tree-level perturbation theory to the numerical data. To achieve this, we estimate the covariance matrix of the measurement errors by using 10,000 mock catalogues generated with the Pinocchio code. We study how the model constraints are influenced by the binning strategy for the bispectrum configurations and by the form of the likelihood function. We also use Bayesian model-selection techniques to single out the optimal theoretical description of our data. We find that a three-parameter bias model combined with Poissonian shot noise is necessary to model the halo bispectrum up to scales of $k_\mathrm{max}\lesssim 0.08 h \mathrm{Mpc}^{-1}$, although fitting formulae that relate the bias parameters can be helpful to reduce the freedom of the model without compromising accuracy. Our data clearly disfavour local Eulerian and local Lagrangian bias models and do not require corrections to Poissonian shot noise. We anticipate that model-selection diagnostics will be particularly useful to extend the analysis to smaller scales as, in this case, the number of model parameters will grow significantly large., Comment: 38 pages, 17 figures; major revisions: inclusion of open triangle bins in all results, new priors for the shot-noise parameters, and additional comparison with Gaussian variance

In this paper we demonstrate that the information encoded in \emph{one} single (sufficiently large) $N$-body simulation can be used to reproduce arbitrary numbers of halo catalogues, using approximated realisations of dark matter density fields with different initial conditions. To this end we use as a reference one realisation (from an ensemble of $300$) of the Minerva $N$-body simulations and the recently published Bias Assignment Method to extract the local and non-local bias linking the halo to the dark matter distribution. We use an approximate (and fast) gravity solver to generate $300$ dark matter density fields from the down-sampled initial conditions of the reference simulation and sample each of these fields using the halo-bias and a kernel, both calibrated from the arbitrarily chosen realisation of the reference simulation. We show that the power spectrum, its variance and the three-point statistics are reproduced within $\sim 2\%$ (up to $k\sim1.0\,h\,{\rm Mpc}^{-1}$), $\sim 5-10\%$ and $\sim 10\%$, respectively. Using a model for the real space power spectrum (with three free bias parameters), we show that the covariance matrices obtained from our procedure lead to parameter uncertainties that are compatible within $\sim 10\%$ with respect to those derived from the reference covariance matrix, and motivate approaches that can help to reduce these differences to $\sim 1\%$. Our method has the potential to learn from one simulation with moderate volumes and high-mass resolution and extrapolate the information of the bias and the kernel to larger volumes, making it ideal for the construction of mock catalogues for present and forthcoming observational campaigns such as Euclid or DESI., Comment: Accepted for publication at MNRAS

Leveraging the Baryon Acoustic Oscillations (BAO) feature present in clustering 2-point statistics, we aim to measure cosmological distances independently of the underlying background cosmological model. However this inference is complicated by late-time non-linearities that introduce model and tracer dependencies in the clustering correlation function and power spectrum, which must be properly accounted for. With this in mind, we introduce the "Purely-Geometric-BAO," which provides a rigorous tool to measure cosmological distances without assuming a specific background cosmology. We focus on the 2-point clustering correlation function monopole, and show how to implement such an inference scheme employing two different methodologies: the Linear Point standard ruler (LP) and correlation-function model-fitting (CF-MF). For the first time we demonstrate how, by means of the CF-MF, we can measure very precisely the sound-horizon/isotropic-volume-distance ratio, $r_{d}/D_{V}(\bar{z})$, while correctly propagating all the uncertainties. Using synthetic data, we compare the outcomes of the two methodologies, and find that the LP provides up to $50\%$ more precise measurements than the CF-MF. Finally, we test a procedure widely employed in BAO analyses: fitting the 2-point function while fixing the cosmological and the non-linear-damping parameters at fiducial values. We find that this underestimates the distance errors by nearly a factor of $2$. We thus recommend that this practice be reconsidered, whether for parameter determination or model selection., Comment: 12 pages, 2 tables

We present a new cosmological test using the distribution of galaxies around cosmic voids without assuming a fiducial cosmology. The test is based on a physical model for the void-galaxy cross-correlation function projected along and perpendicular to the line of sight. We treat correlations in terms of void-centric angular distances and redshift differences between void-galaxy pairs, hence it is not necessary to assume a fiducial cosmology. This model reproduces the coupled dynamical (Kaiser effect, RSD) and geometrical (Alcock-Paczynski effect, GD) distortions that affect the correlation measurements. It also takes into account the scale mixing due to the projection ranges in both directions. The model is general, so it can be applied to an arbitrary cylindrical binning scheme, not only in the case of the projected correlations. It primarily depends on two cosmological parameters: $\Omega_m$, the matter fraction of the Universe today (sensitive to GD), and $\beta$, the ratio between the growth rate factor of density perturbations and the tracer bias (sensitive to RSD). In the context of the new generation of galaxy spectroscopic surveys, we calibrated the test using the Millennium XXL simulation for different redshifts. The method successfully recovers the cosmological parameters. We studied the effect of measuring with different projection ranges, finding robust results up to wide ranges. The resulting data covariance matrices are relatively small, which reduces the noise in the Gaussian likelihood analysis and will allow the usage of a smaller number of mock catalogues. The performance evaluated in this work indicates that the developed method is a promising test to be applied on real data., Comment: 13 pages. 7 Figures. Accepted in MNRAS

Measurements of redshift space distortions (RSD) provide a means to test models of gravity on large-scales. We use mock galaxy catalogues constructed from large N-body simulations of standard and modified gravity models to measure galaxy clustering in redshift space. We focus our attention on two of the most representative and popular families of modified gravity models: the Hu \& Sawicki $f(R)$ gravity and the normal branch of the DGP model. The galaxy catalogues are built using a halo occupation distribution (HOD) prescription with the HOD parameters in the modified gravity models tuned to match with the number density and the real-space clustering of {\sc boss-cmass} galaxies. We employ two approaches to model RSD: the first is based on linear perturbation theory and the second models non-linear effects on small-scales by assuming standard gravity and including biasing and RSD effects. We measure the monopole to real-space correlation function ratio, the quadrupole to monopole ratio, clustering wedges and multipoles of the correlation function and use these statistics to find the constraints on the distortion parameter, $\beta$. We find that the linear model fails to reproduce the N-body simulation results and the true value of $\beta$ on scales $s < 40\Mpch$, while the non-linear modelling of RSD recovers the value of $\beta$ on the scales of interest for all models. RSD on large scales ($s\gtrsim20$-$40\Mpch$) have been found to show significant deviations from the prediction of standard gravity in the DGP models. However, the potential to use RSD to constrain $f(R)$ models is less promising, due to the different screening mechanism in this model., Comment: 18 pages, 11 figures, submitted to MNRAS

We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method., Comment: Additional results with respect to v1, new section and new figures added. 25 pages, 1 table 18 figures

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance., Comment: 20 pages, 16 figures, replaced to match accepted MNRAS version. Results on parameter errors changed

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others., Comment: 23 pages, 11 figures. Replaced to match accepted MNRAS version. Included Kullback-Leibler divergence

We present the clustering measurements of quasars in configuration space based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach with a Gaussian Streaming (GS) model for the redshift space distortions of the correlation function and demonstrate its applicability for dark matter halos hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$, we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426 \pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{-13}(r_{s}^{\rm fid}/r_s){\rm km.s}^{-1}.{\rm Mpc}^{-1}$, and the angular diameter distance $D_{A}(z_{\rm eff})=1850^{+90}_{-115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where $r_{s}$ is the sound horizon at the end of the baryon drag epoch and $r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors include both systematic and statistical contributions. The results on the evolution of distances are consistent with the predictions of flat $\Lambda$-Cold Dark Matter ($\Lambda$-CDM) cosmology with Planck parameters, and the measurement of $f\sigma_{8}$ extends the validity of General Relativity (GR) to higher redshifts($z>1$) This paper is released with companion papers using the same sample. The results on the cosmological parameters of the studies are found to be in very good agreement, providing clear evidence of the complementarity and of the robustness of the first full-shape clustering measurements with the eBOSS DR14 quasar sample., Comment: 25 pages, 22 figures

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8

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

We analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample (DR14Q). We measure the redshift space distortions using the power spectrum monopole, quadrupole and hexadecapole inferred from 148,659 quasars between redshifts 0.8 and 2.2 covering a total sky footprint of 2112.9 deg$^2$. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, $f\sigma_8$, and the Alcock-Paczynski dilation scales which allow constraints to be placed on the angular diameter distance $D_A(z)$ and the Hubble $H(z)$ parameter. At the effective redshift of $z_{\rm eff}=1.52$, $f\sigma_8(z_{\rm eff})=0.420\pm0.076$, $H(z_{\rm eff})=[162\pm 12]\, (r_s^{\rm fid}/r_s)\,{\rm km\, s}^{-1}{\rm Mpc}^{-1}$, and $D_A(z_{\rm eff})=[1.85\pm 0.11]\times10^3\,(r_s/r_s^{\rm fid})\,{\rm Mpc}$, where $r_s$ is the comoving sound horizon at the baryon drag epoch and the superscript `fid' stands for its fiducial value. The errors take into account the full error budget, including systematics and statistical contributions. These results are in full agreement with the current $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model inferred from Planck measurements. Finally, we compare our measurements with other eBOSS companion papers and find excellent agreement, demonstrating the consistency and complementarity of the different methods used for analysing the data., Comment: 37 pages, 25 figures, 14 tables. Data and covariance will be made public after the acceptance of the paper. Submitted to MNRAS

We explore the cosmological implications of anisotropic clustering measurements of the quasar sample from Data Release 14 of the Sloan Digital Sky Survey IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in configuration space. The $\sim 147,000$ quasar sample observed by eBOSS offers a direct tracer of the density field and bridges the gap of previous BAO measurements between redshift $0.8

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V., Comment: SDSS-IV collaboration alphabetical author data release paper. DR14 happened on 31st July 2017. 19 pages, 5 figures. Accepted by ApJS on 28th Nov 2017 (this is the "post-print" and "post-proofs" version; minor corrections only from v1, and most of errors found in proofs corrected)

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature., Comment: Accepted by MNRAS; BAO distance likelihood available in source files 'QSOv1.9fEZmock_BAOchi2.dat'; full set of data to be public eventually from SDSS website

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016., Comment: Published in Astronomical Journal

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey., Comment: 27 pages, 3 figures and one table. A supplementary document is included. The BOSS DR12 BAO data used in the work can be downloaded from the SDSS website https://sdss3.org//science/boss_publications.php (check https://data.sdss.org/sas/dr12/boss/papers/clustering/Zhao_etal_2016_DR12COMBINED_tomoBAO_powspec.tar.gz for the tomographic BAO measurement presented in https://arxiv.org/abs/1607.03153)

We use baryon acoustic oscillation and redshift space distortion from the completed Baryon Oscillation Spectroscopic Survey, corresponding to data release 12 of the Sloan Digital Sky Survey, combined sample analysis in combination with cosmic microwave background, supernova and redshift space distortion measurements from additional spectroscopic surveys to test deviations from general relativity. We present constraints on several phenomenological models of modified gravity: First, we parametrise the growth of structure using the growth index $\gamma$, finding $\gamma=0.566\pm0.058$ (68\% C.L.). Second, we modify the relation of the two Newtonian potentials by introducing two additional parameters, $G_M$ and $G_L$. In this approach, $G_M$ refers to modifications of the growth of structure whereas $G_L$ to modification of the lensing potential. We consider a power law to model the redshift dependency of $G_M$ and $G_L$ as well as binning in redshift space, introducing four additional degrees of freedom, $G_M(z<0.5)$, $G_M(z>0.5)$, $G_L(z<0.5)$, $G_L(z>0.5)$. At 68\% C.L. we measure $G_M=0.980 \pm 0.096$ and $G_L=1.082 \pm 0.060$ for a linear model, $G_M= 1.01\pm 0.36$ and $G_L=1.31 \pm 0.19$ for a cubic model as well as $G_M(z<0.5)=1.26 \pm 0.32$, $G_M(z>0.5)=0.986\pm0.022$, $G_L(z<0.5)=1.067\pm0.058$ and $G_L(z>0.5)=1.037\pm0.029$. Thirdly, we investigate general scalar tensor theories of gravity, finding the model to be mostly unconstrained by current data. Assuming a one-parameter $f(R)$ model we can constrain $B_0< 7.7 \times 10^{-5}$ (95\% C.L). For all models we considered we find good agreement with general relativity., Comment: 10 pages, 7 figures, submitted to MNRAS

We investigate the potential sources of theoretical systematics in the anisotropic Baryon Acoustic Oscillation (BAO) distance scale measurements from the clustering of galaxies in configuration space using the final Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). We perform a detailed study of the impact on BAO measurements from choices in the methodology such as fiducial cosmology, clustering estimators, random catalogues, fitting templates, and covariance matrices. The theoretical systematic uncertainties in BAO parameters are found to be 0.002 in the isotropic dilation $\alpha$ and 0.003 in the quadrupolar dilation $\epsilon$. The leading source of systematic uncertainty is related to the reconstruction techniques. Theoretical uncertainties are sub-dominant compared with the statistical uncertainties for BOSS survey, accounting $0.2\sigma_{stat}$ for $\alpha$ and $0.25\sigma_{stat}$ for $\epsilon$ ($\sigma_{\alpha,stat} \sim$0.010 and $\sigma_{\epsilon,stat}\sim$ 0.012 respectively). We also present BAO-only distance scale constraints from the anisotropic analysis of the correlation function. Our constraints on the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$, including both statistical and theoretical systematic uncertainties, are 1.5\% and 2.8\% at $z_{\rm eff}=0.38$, 1.4\% and 2.4\% at $z_{\rm eff}=0.51$, and 1.7\% and 2.6\% at $z_{\rm eff}=0.61$. This paper is part of a set that analyzes the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are cross-checked with other BAO analysis in \citet{Acacia16}. The systematic error budget concerning the methodology on post-reconstruction BAO analysis presented here is used in \citet{Acacia16} to produce the final cosmological constraints from BOSS., Comment: 38 pages, 22 figures, accepted for publication in MNRAS