We examine the cosmological implications of measurements of the void-galaxy cross-correlation at redshift $z=0.57$ combined with baryon acoustic oscillation (BAO) data at $0.1<z<2.4$. We find direct evidence of the late-time acceleration due to dark energy at $>10\sigma$ significance from these data alone, independent of the cosmic microwave background and supernovae. Using a nucleosynthesis prior on $\Omega_bh^2$, we measure the Hubble constant to be $H_0=72.3\pm1.9\;{\rm km\,s}^{-1}{\rm Mpc}^{-1}$ from BAO+voids at $z<2$, and $H_0=69.0\pm1.2\;{\rm km\,s}^{-1}{\rm Mpc}^{-1}$ when adding Lyman-$\alpha$ BAO at $z=2.34$, both independent of the CMB. Adding voids to CMB, BAO and supernova data greatly improves measurement of the dark energy equation of state, increasing the figure of merit by >40%, but remaining consistent with flat flat $\Lambda$ cold dark matter.

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

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

We present evidence of anisotropic galaxy assembly bias in the Baryon Oscillation Spectroscopic Survey Data Release 12 galaxy sample at a level exceeding $5\sigma$. We use measurements of the line-of-sight velocity dispersion $\sigma_\star$ and stellar mass $M_\star$ to perform a simple split into subsamples of galaxies. We show that the amplitude of the monopole and quadrupole moments of the power spectrum depend differently on $\sigma_\star$ and $M_\star$, allowing us to split the galaxy sample into subsets with matching monopoles but significantly different quadrupoles on all scales. Combining data from the LOWZ and CMASS NGC galaxy samples, we find $>5\sigma$ evidence for anisotropic bias on scales $k<0.15\,h\,\rm{Mpc}^{-1}$. We also examine splits using other observed properties. For galaxy samples split using $M_\star$ and projected size $R_0$, we find no significant evidence of anisotropic bias. Galaxy samples selected using additional properties exhibit strongly varying degrees of anisotropic assembly bias, depending on which combination of properties is used to split into subsets. This may explain why previous searches for this effect using the Fundamental Plane found inconsistent results. We conclude that any selection of a galaxy sample that depends on $\sigma_\star$ can give biased and incorrect Redshift Space Distortion measurements.

Existing cosmological simulation methods lack a high degree of parallelism due to the long-range nature of the gravitational force, which limits the size of simulations that can be run at high resolution. To solve this problem, we propose a new, perfectly parallel approach to simulate cosmic structure formation, based on the spatial COmoving Lagrangian Acceleration (sCOLA) framework. Building upon a hybrid analytical and numerical description of particles' trajectories, our algorithm allows an efficient tiling of a cosmological volume, where the dynamics within each tile is computed independently. As a consequence, the degree of parallelism is equal to the number of tiles. We optimise the accuracy of sCOLA by the use of a buffer region around tiles, and of appropriate Dirichlet boundary conditions around sCOLA boxes. As a result, we show that cosmological simulations at the degree of accuracy required for the analysis of the next generation of surveys can be run in drastically reduced wall-clock times and with very low memory requirements. The perfect scalability of our algorithm unlocks profoundly new possibilities of computing larger and higher-resolution cosmological simulations, taking advantage of a variety of hardware architectures.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

We measure the imprint of primordial baryon acoustic oscillations (BAO) in the correlation function of Ly$\alpha$ absorption in quasar spectra from the Baryon Oscillation Spectroscopic Survey (BOSS) and the extended BOSS (eBOSS) in Data Release 14 (DR14) of the Sloan Digital Sky Survey (SDSS)-IV. In addition to 179,965 spectra with absorption in the Lyman-$\alpha$ (Ly$\alpha$) region, we use, for the first time, Ly$\alpha$ absorption in the Lyman-$\beta$ region of 56,154 spectra. We measure the Hubble distance, $D_H$, and the comoving angular diameter distance, $D_M$, relative to the sound horizon at the drag epoch $r_d$ at an effective redshift $z=2.34$. Using a physical model of the correlation function outside the BAO peak, we find $D_H(2.34)/r_d=8.86\pm 0.29$ and $D_M(2.34)/r_d=37.41\pm 1.86$, within 1$\sigma$ from the flat-$\Lambda$CDM model consistent with CMB anisotropy measurements. With the addition of polynomial "broadband" terms, the results remain within one standard deviation of the CMB-inspired model. Combined with the quasar-Ly$\alpha$ cross-correlation measurement presented in a companion paper Blomqvist19, the BAO measurements at $z=2.35$ are within 1.7$\sigma$ of the predictions of this model.

We present a measurement of the baryon acoustic oscillation (BAO) scale at redshift $z=2.35$ from the three-dimensional correlation of Lyman-$\alpha$ (Ly$\alpha$) forest absorption and quasars. The study uses 266,590 quasars in the redshift range $1.77<z<3.5$ from the Sloan Digital Sky Survey (SDSS) Data Release 14 (DR14). The sample includes the first two years of observations by the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS), providing new quasars and re-observations of BOSS quasars for improved statistical precision. Statistics are further improved by including Ly$\alpha$ absorption occurring in the Ly$\beta$ wavelength band of the spectra. From the measured BAO peak position along and across the line of sight, we determined the Hubble distance $D_{H}$ and the comoving angular diameter distance $D_{M}$ relative to the sound horizon at the drag epoch $r_{d}$: $D_{H}(z=2.35)/r_{d}=9.20\pm 0.36$ and $D_{M}(z=2.35)/r_{d}=36.3\pm 1.8$. These results are consistent at $1.5\sigma$ with the prediction of the best-fit spatially-flat cosmological model with the cosmological constant reported for the Planck (2016) analysis of cosmic microwave background anisotropies. Combined with the Ly$\alpha$ auto-correlation measurement presented in a companion paper, the BAO measurements at $z=2.34$ are within $1.7\sigma$ of the predictions of this model.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

We study the effect of large-scale tidal fields on internal halo properties using a set of N-body simulations. We measure significant cross-correlations between large-scale tidal fields and several non-scalar halo properties: shapes, velocity dispersion, and angular momentum. Selection effects that couple to these non-scalar halo properties can produce anisotropic clustering even in real-space. We investigate the size of this effect and show that it can produce a non-zero quadrupole similar in size to the one generated by linear redshift-space distortions (RSD). Finally, we investigate the clustering properties of halos identified in redshift-space and find enormous deviations from the standard linear RSD model, again caused by anisotropic assembly bias. These effects could contaminate the values of cosmological parameters inferred from the observed redshift-space clustering of galaxies, groups, or 21cm emission from atomic hydrogen, if their selection depends on properties affected by halo assembly bias. We briefly discuss ways in which this effect can be measured in existing and future large-scale structure surveys.

We present a measurement of the anisotropic void-galaxy cross-correlation function in the CMASS galaxy sample of the BOSS DR12 data release. We perform a joint fit to the data for redshift space distortions (RSD) due to galaxy peculiar velocities and anisotropies due to the Alcock-Paczynski (AP) effect, for the first time using a velocity field reconstruction technique to remove the complicating effects of RSD in the void centre positions themselves. Fits to the void-galaxy function give a 1% measurement of the AP parameter combination $D_A(z)H(z)/c = 0.4367\pm 0.0045$ at redshift $z=0.57$, where $D_A$ is the angular diameter distance and $H$ the Hubble parameter, exceeding the precision obtainable from baryon acoustic oscillations (BAO) by a factor of ~3.5 and free of systematic errors. From voids alone we also obtain a 10% measure of the growth rate, $f\sigma_8(z=0.57)=0.501\pm0.051$. The parameter degeneracies are orthogonal to those obtained from galaxy clustering. Combining void information with that from BAO and galaxy RSD in the same CMASS sample, we measure $D_A(0.57)/r_s=9.383\pm 0.077$ (at 0.8% precision), $H(0.57)r_s=(14.05\pm 0.14)\;10^3$ kms$^{-1}$Mpc$^{-1}$ (1%) and $f\sigma_8=0.453\pm0.022$ (4.9%), consistent with cosmic microwave background (CMB) measurements from Planck. These represent a factor \sim2 improvement in precision over previous results through the inclusion of void information. Fitting a flat cosmological constant $\Lambda$CDM model to these results in combination with Planck CMB data, we find up to an 11% reduction in uncertainties on $H_0$ and $\Omega_m$ compared to use of the corresponding BOSS consensus values. Constraints on extended models with non-flat geometry and a dark energy of state that differs from $w=-1$ show an even greater improvement.

Standard analysis pipelines for measurements of Baryon Acoustic Oscillations (BAO) in galaxy surveys make use of a fiducial cosmological model to guide the data compression required to transform from observed redshifts and angles to the measured angular and radial BAO peak positions. In order to remove any dependence on the fiducial cosmology from the results, all models compared to the data should mimic the compression and its dependence on the fiducial model. In practice, approximations are made when testing models: (1) There is assumed to be no residual dependence on the fiducial cosmology after reconstruction, (2) differences in the distance--redshift relationship are assumed to match a linear scaling, and (3) differences in clustering between true and fiducial models are assumed to be removed by the free parameters used to null the non-BAO signal. We test these approximations using the current standard measurement procedure with a set of halo catalogs from the {\sc Aemulus} suite of $N$-body simulations, which span a range of $w\mathrm{CDM}$ cosmological models. We focus on reconstruction of the primordial BAO and locating the BAO. For the range of $w\mathrm{CDM}$ cosmologies covered by the {\sc Aemulus} suite, we find no evidence for systematic errors in the measured BAO shift parameters $\alpha_{\parallel}$ and $\alpha_{\bot}$ to $< 0.1\%$. However, the measured errors $\sigma_{\alpha_{\parallel}}$ and $\sigma_{\alpha_{\bot}}$ show a notable absolute increase by up to $+0.001$ and $+0.002$ respectively in the case that the fiducial cosmology does not match the truth. These effects on the inferred BAO scale will be important given the precision of measurements expected from future surveys including DESI, Euclid, and WFIRST.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

We describe the first public data release of the Dark Energy Survey, DES DR1, consisting of reduced single epoch images, coadded images, coadded source catalogs, and associated products and services assembled over the first three years of DES science operations. DES DR1 is based on optical/near-infrared imaging from 345 distinct nights (August 2013 to February 2016) by the Dark Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread function of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a photometric precision of < 1% in all bands, and an astrometric precision of 151 mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10 is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1 includes nearly 400M distinct astronomical objects detected in ~10,000 coadd tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures. Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects, respectively, following a basic object quality selection. These data are accessible through a range of interfaces, including query web clients, image cutout servers, jupyter notebooks, and an interactive coadd image visualization tool. DES DR1 constitutes the largest photometric data set to date at the achieved depth and photometric precision.

Likelihood fitting to two-point clustering statistics made from galaxy surveys usually assumes a multivariate normal distribution for the measurements, with justification based on the central limit theorem given the large number of overdensity modes. However, this assumption cannot hold on the largest scales where the number of modes is low. Whilst more accurate distributions have previously been developed in idealized cases, we derive a procedure suitable for analysing measured monopole power spectra with window effects, stochastic shot noise, and the dependence of the covariance matrix on the model being fitted all taken into account. A data transformation is proposed to give an approximately Gaussian likelihood, with a variance--correlation decomposition of the covariance matrix to account for its cosmological dependence. By comparing with the modified-$t$ likelihood derived under the usual normality assumption, we find in numerical tests that our new procedure gives more accurate constraints on the local non-Gaussianity parameter $f_\mathrm{NL}$, which is sensitive to the large-scale power. A simple data analysis pipeline is provided for straightforward application of this new approach in preparation for forthcoming large galaxy surveys such as DESI and Euclid.

We present constraints on local primordial non-Gaussianity (PNG), parametrized through $f^{\rm loc}_{\rm NL}$, using the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample. We measure and analyze the anisotropic clustering of the quasars in Fourier space, testing for the scale-dependent bias introduced by primordial non-Gaussianity on large scales. We derive and employ a power spectrum estimator using optimal weights that account for the redshift evolution of the PNG signal. We find constraints of $-51<f^{\rm loc}_{\rm NL}<21$ at 95% confidence level. These are amont the tightest constraints from Large Scale Structure (LSS) data. Our redshift weighting improves the error bar by 15% in comparison to the unweighted case. If quasars have lower response to PNG, the constraint degrades to $-81<f^{\rm loc}_{\rm NL}<26$, with a 40% improvement over the standard approach. We forecast that the full eBOSS dataset could reach $\sigma_{f^{\rm loc}_{\rm NL}}\simeq 5\text{-}8$ using optimal methods and full range of scales.

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

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.

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

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 nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

The Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) will be a survey of bright, low redshift galaxies, which is planned to cover an area of ~14,000 sq deg in 3 passes. Each pass will cover the survey area with ~2000 pointings, each of area ~8 sq deg. The BGS is currently proposed to consist of a bright high priority sample to an r-band magnitude limit r ~ 19.5, with a fainter low priority sample to r ~ 20. The geometry of the DESI fibre positioners in the focal plane of the telescope affects the completeness of the survey, and has a non-trivial impact on clustering measurements. Using a BGS mock catalogue, we show that completeness due to fibre assignment primarily depends on the surface density of galaxies. Completeness is high (>95%) in low density regions, but very low (<10%) in the centre of massive clusters. We apply the pair inverse probability (PIP) weighting correction to clustering measurements from a BGS mock which has been through the fibre assignment algorithm. This method is only unbiased if it is possible to observe every galaxy pair. To facilitate this, we randomly promote a small fraction of the fainter sample to be high priority, and dither the set of tile positions by a small angle. We show that inverse pair weighting combined with angular upweighting provides an unbiased correction to galaxy clustering measurements for the complete 3 pass survey, and also after 1 pass, which is highly incomplete.

We define and characterise a sample of 1.3 million galaxies extracted from the first year of Dark Energy Survey data, optimised to measure Baryon Acoustic Oscillations in the presence of significant redshift uncertainties. The sample is dominated by luminous red galaxies located at redshifts $z \gtrsim 0.6$. We define the exact selection using color and magnitude cuts that balance the need of high number densities and small photometric redshift uncertainties, using the corresponding forecasted BAO distance error as a figure-of-merit in the process. The typical photo-$z$ uncertainty varies from $2.3\%$ to $3.6\%$ (in units of 1+$z$) from $z=0.6$ to $1$, with number densities from $200$ to $130$ galaxies per deg$^2$ in tomographic bins of width $\Delta z = 0.1$. Next we summarise the validation of the photometric redshift estimation. We characterise and mitigate observational systematics including stellar contamination, and show that the clustering on large scales is robust in front of those contaminants. We show that the clustering signal in the auto-correlations and cross-correlations is generally consistent with theoretical models, which serves as an additional test of the redshift distributions.

Redshift space distortions within voids provide a unique method to test for environmental dependence of the growth rate of structures in low density regions, where effects of modified gravity theories might be important. We derive a linear theory model for the redshift space void-galaxy correlation that is valid at all pair separations, including deep within the void, and use this to obtain expressions for the monopole $\xi^s_0$ and quadrupole $\xi^s_2$ contributions. Our derivation highlights terms that have previously been neglected but are important within the void interior. As a result our model differs from previous works and predicts new physical effects, including a change in the sign of the quadrupole term within the void radius. We show how the model can be generalised to include a velocity dispersion. We compare our model predictions to measurements of the correlation function using mock void and galaxy catalogues modelled after the BOSS CMASS galaxy sample using the Big MultiDark $N$-body simulation, and show that the linear model with dispersion provides an excellent fit to the data at all scales, $0\leq s\leq120\;h^{-1}$Mpc. While the RSD model matches simulations, the linear bias approximation does not hold within voids, and care is needed in fitting for the growth rate. We show that fits to the redshift space correlation recover the growth rate $f(z=0.52)$ to a precision of $2.7\%$ using the simulation volume of $(2.5\;h^{-1}\mathrm{Gpc})^3$.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We present angular diameter distance measurements obtained by locating the BAO scale in the distribution of galaxies selected from the first year of Dark Energy Survey data. We consider a sample of over 1.3 million galaxies distributed over a footprint of 1318 deg$^2$ with $0.6 < z_{\rm photo} < 1$ and a typical redshift uncertainty of $0.03(1+z)$. This sample was selected, as fully described in a companion paper, using a color/magnitude selection that optimizes trade-offs between number density and redshift uncertainty. We investigate the BAO signal in the projected clustering using three conventions, the angular separation, the co-moving transverse separation, and spherical harmonics. Further, we compare results obtained from template based and machine learning photometric redshift determinations. We use 1800 simulations that approximate our sample in order to produce covariance matrices and allow us to validate our distance scale measurement methodology. We measure the angular diameter distance, $D_A$, at the effective redshift of our sample divided by the true physical scale of the BAO feature, $r_{\rm d}$. We obtain close to a 4 per cent distance measurement of $D_A(z_{\rm eff}=0.81)/r_{\rm d} = 10.75\pm 0.43 $. These results are consistent with the flat $\Lambda$CDM concordance cosmological model supported by numerous other recent experimental results. All data products are publicly available here: https://des.ncsa.illinois.edu/releases/y1a1/bao

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<z<2.2$. We report a joint measurement of $f\sigma_8$ and two-dimensional BAO parameters $D_{\rm A}$ and $H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of dark energy is supported by observations at a $7.4\sigma$ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

Redshift space distortions (RSD) in the void-galaxy correlation $\xi^s$ provide information on the linear growth rate of structure in low density environments. Accurate modelling of these RSD effects can also allow the use of voids in competitive Alcock-Paczynski measurements. Linear theory models of $\xi^s$ are able to provide extremely good descriptions of simulation data on all scales provided the real space void positions are known. However, by reference to simulation data we demonstrate the failure of the assumptions implicit in current models of $\xi^s$ for voids identified directly in redshift space, as would be simplest using real observational data. To overcome this problem we instead propose using a density-field reconstruction method based on the Zeldovich approximation to recover the real space void positions from redshift space data. We show that this recovers the excellent agreement between theory and data for $\xi^s$. Performing the reconstruction requires an input cosmological model so, to be self-consistent, we have to perform reconstruction for every model to be tested. We apply this method to mock galaxy and void catalogues in the Big MultiDark $N$-body simulation and consistently recover the fiducial growth rate to a precision of $3.4\%$ using the simulation volume of $(2.5\;h^{-1}\mathrm{Gpc})^3$.

The Emission Line Galaxy survey made by the Dark Energy Spectroscopic Instrument (DESI) survey will be created from five passes of the instrument on the sky. On each pass, the constrained mobility of the ends of the fibres in the DESI focal plane means that the angular-distribution of targets that can be observed is limited. Thus, the clustering of samples constructed using a limited number of passes will be strongly affected by missing targets. In two recent papers, we showed how the effect of missing galaxies can be corrected when calculating the correlation function using a weighting scheme for pairs. Using mock galaxy catalogues we now show that this method provides an unbiased estimator of the true correlation function for the DESI survey after any number of passes. We use multiple mocks to determine the expected errors given one to four passes, compared to an idealised survey observing an equivalent number of randomly selected targets. On BAO scales, we find that the error is a factor 2 worse after one pass, but that after three or more passes, the errors are very similar. Thus we find that the fibre assignment strategy enforced by the design of DESI will not affect the cosmological measurements to be made by the survey, and can be removed as a potential risk for this experiment.

These are advanced lecture notes covering recent developments in the methodology used to analyse galaxy surveys. The focus is particularly on direct measurements of the galaxy power spectrum although I also discuss its Fourier transform, the correlation function for comparison. These 2-point statistics, under the assumption that the overdensity field has Gaussian statistics on large-scales, contain the majority of the cosmological signal available from the galaxy distribution. Recent developments in multipole measurements, dealing with systematics, convolving theoretical models with the survey window function, the approximation of covariance matrices, and weighting schemes for measuring evolution with redshift are considered. The focus is on analytic explanation of the issues involved rather than on recent analyses or simulation results.

We develop a practical methodology to remove modes from a galaxy survey power spectrum that are associated with systematic errors. We apply this to the BOSS CMASS sample, to see if it removes the excess power previously observed beyond the best-fit $\Lambda$CDM model on very large scales. We consider several possible sources of data contamination, and check whether they affect the number of targets that can be observed and the power spectrum measurements. We describe a general framework for how such knowledge can be transformed into template fields. Mode subtraction can then be used to remove these systematic contaminants at least as well as applying corrective weighting to the observed galaxies, but benefits from giving an unbiased power. Even after applying templates for all known systematics, we find a large-scale power excess, but this is reduced compared with that observed using the weights provided by the BOSS team. This excess is at much larger scales than the BAO scale and does not affect the main results of BOSS. However, it will be important for the measurement of a scale-dependent bias due to primordial non-Gaussianity. The excess is beyond that allowed by any simple model of non-Gaussianity matching Planck data, and is not matched in other surveys. We show that this power excess can further be reduced but is still present using "phenomenological" templates, designed to consider further potentially unknown sources of systematic contamination. As all discrepant angular modes can be removed using "phenomenological" templates, the potentially remaining contaminant acts radially.

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

We measure the clustering of DES Year 1 galaxies that are intended to be combined with weak lensing samples in order to produce precise cosmological constraints from the joint analysis of large-scale structure and lensing correlations. Two-point correlation functions are measured for a sample of $6.6 \times 10^{5}$ luminous red galaxies selected using the \textsc{redMaGiC} algorithm over an area of $1321$ square degrees, in the redshift range $0.15 < z < 0.9$, split into five tomographic redshift bins. The sample has a mean redshift uncertainty of $\sigma_{z}/(1+z) = 0.017$. We quantify and correct spurious correlations induced by spatially variable survey properties, testing their impact on the clustering measurements and covariance. We demonstrate the sample's robustness by testing for stellar contamination, for potential biases that could arise from the systematic correction, and for the consistency between the two-point auto- and cross-correlation functions. We show that the corrections we apply have a significant impact on the resultant measurement of cosmological parameters, but that the results are robust against arbitrary choices in the correction method. We find the linear galaxy bias in each redshift bin in a fiducial cosmology to be $b(z$=$0.24)=1.40 \pm 0.08$, $b(z$=$0.38)=1.61 \pm 0.05$, $b(z$=$0.53)=1.60 \pm 0.04$ for galaxies with luminosities $L/L_*>$$0.5$, $b(z$=$0.68)=1.93 \pm 0.05$ for $L/L_*>$$1$ and $b(z$=$0.83)=1.99 \pm 0.07$ for $L/L_*$$>1.5$, broadly consistent with expectations for the redshift and luminosity dependence of the bias of red galaxies. We show these measurements to be consistent with the linear bias obtained from tangential shear measurements.

We study the methodology and potential theoretical systematics of measuring Baryon Acoustic Oscillations (BAO) using the angular correlation functions in tomographic bins. We calibrate and optimize the pipeline for the Dark Energy Survey Year 1 dataset using 1800 mocks. We compare the BAO fitting results obtained with three estimators: the Maximum Likelihood Estimator (MLE), Profile Likelihood, and Markov Chain Monte Carlo. The fit results from the MLE are the least biased and their derived 1-$\sigma$ error bar are closest to the Gaussian distribution value after removing the extreme mocks with non-detected BAO signal. We show that incorrect assumptions in constructing the template, such as mismatches from the cosmology of the mocks or the underlying photo-$z$ errors, can lead to BAO angular shifts. We find that MLE is the method that best traces this systematic biases, allowing to recover the true angular distance values. In a real survey analysis, it may happen that the final data sample properties are slightly different from those of the mock catalog. We show that the effect on the mock covariance due to the sample differences can be corrected with the help of the Gaussian covariance matrix or more effectively using the eigenmode expansion of the mock covariance. In the eigenmode expansion, the eigenmodes are provided by some proxy covariance matrix. The eigenmode expansion is significantly less susceptible to statistical fluctuations relative to the direct measurements of the covariance matrix because of the number of free parameters is substantially reduced

Mock catalogues are a crucial tool in the analysis of galaxy surveys data, both for the accurate computation of covariance matrices, and for the optimisation of analysis methodology and validation of data sets. In this paper, we present a set of 1800 galaxy mock catalogues designed to match the Dark Energy Survey Year-1 BAO sample (Crocce et al. 2017) in abundance, observational volume, redshift distribution and uncertainty, and redshift dependent clustering. The simulated samples were built upon HALOGEN (Avila et al. 2015) halo catalogues, based on a $2LPT$ density field with an exponential bias. For each of them, a lightcone is constructed by the superposition of snapshots in the redshift range $0.45<z<1.4$. Uncertainties introduced by so-called photometric redshifts estimators were modelled with a \textit{double-skewed-Gaussian} curve fitted to the data. We also introduce a hybrid HOD-HAM model with two free parameters that are adjusted to achieve a galaxy bias evolution $b(z_{\rm ph})$ that matches the data at the 1-$\sigma$ level in the range $0.6<z_{\rm ph}<1.0$. We further analyse the galaxy mock catalogues and compare their clustering to the data using the angular correlation function $ w(\theta)$, the comoving transverse separation clustering $\xi_{\mu<0.8}(s_{\perp})$ and the angular power spectrum $C_\ell$.

We use data from the first-year (Y1) observations of the DES collaboration to measure the galaxy angular power spectrum (APS), and search for its BAO feature using a template-fitting method. We test our methodology in a sample of 1800 DES Y1-like mock catalogs. The APS is measured with the pseudo-$C_\ell$ method, using pixelized maps constructed from the mock catalogs and the DES mask. The covariance matrix of the $C_\ell$'s in these tests are also obtained from the mock catalogs. We use templates to model the measured spectra and estimate template parameters firstly from the $C_\ell$'s of the mocks using two different methods, a maximum likelihood estimator and a MCMC, finding consistent results with a good reduced $\chi^2$. Robustness tests are performed to estimate the impact of different choices of settings used in our analysis. After these tests on mocks, we apply our method to a galaxy sample constructed from DES Y1 data specifically for LSS studies. This catalog comprises galaxies within an effective area of 1318 deg$^2$ and $0.6<z<1.0$. We fit the observed spectra with our optimized templates, considering models with and without BAO features. We find that the DES Y1 data favors a model with BAO at the $2.6\,\sigma$ C.L. with a best-fit shift parameter of $\alpha=1.023\pm 0.047$. However, the goodness-of-fit is somewhat poor, with $\chi^2/$(dof) = 1.49. We identify a possible cause of this issue and show that using a theoretical covariance matrix obtained from $C_\ell$'s that are better adjusted to data results in an improved value of $\chi^2/$(dof) = 1.36 which is similar to the value obtained with the real-space analysis. Our results correspond to a distance measurement of $D_A(z_{\rm eff} = 0.81)/r_d = 10.65 \pm 0.49$, consistent with the main DES BAO findings. This is a companion paper to the main DES BAO article showing the details of the harmonic space analysis.

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<z<2.2$. By analysing the two-point correlation function characterized by clustering wedges $\xi_{\rm w_i}(s)$ and multipoles $\xi_{\ell}(s)$, we measure the angular diameter distance, Hubble parameter and cosmic structure growth rate. We define a systematic error budget for our measurements based on the analysis of $N$-body simulations and mock catalogues. Based on the DR14 large scale structure quasar sample at the effective redshift $z_{\rm eff}=1.52$, we find the growth rate of cosmic structure $f\sigma_8(z_{\rm eff})=0.396\pm 0.079$, and the geometric parameters $D_{\rm V}(z)/r_{\rm d}=26.47\pm 1.23$, and $F_{\rm AP}(z)=2.53\pm 0.22$, where the uncertainties include both statistical and systematic errors. These values are in excellent agreement with the best-fitting standard ${\rm \Lambda CDM}$ model to the latest cosmic microwave background data from Planck.

We present measurements of the cross-correlation of the triply-ionized carbon (CIV) forest with quasars using Sloan Digital Sky Survey Data Release 14. The study exploits a large sample of new quasars from the first two years of observations by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The CIV forest is a weaker tracer of large-scale structure than the Ly$\alpha$ forest, but benefits from being accessible at redshifts $z<2$ where the quasar number density from eBOSS is high. Our data sample consists of 287,651 CIV forest quasars in the redshift range $1.4<z<3.5$ and 387,315 tracer quasars with $1.2<z<3.5$. We measure large-scale correlations from CIV absorption occuring in three distinct quasar rest-frame wavelength bands of the spectra referred to as the CIV forest, the SiIV forest and the Ly$\alpha$ forest. From the combined fit to the quasar-CIV cross-correlations for the CIV forest and the SiIV forest, the CIV redshift-space distortion parameter is $\beta_{\rm CIV}=0.27_{\ -0.14}^{\ +0.16}$ and its combination with the CIV linear transmission bias parameter is $b_{\rm CIV}(1+\beta_{\rm CIV})=-0.0183_{\ -0.0014}^{\ +0.0013}$ ($1\sigma$ statistical error) at the mean redshift $z=2.00$. Splitting the sample at $z=2.2$ to constrain the bias evolution with redshift yields the power-law exponent $\gamma=0.60\pm0.63$, indicating a significantly weaker redshift-evolution than for the Ly$\alpha$ forest linear transmission bias. We demonstrate that CIV absorption has the potential to be used as a probe of baryon acoustic oscillations (BAO). While the current data set is insufficient for a detection of the BAO peak feature, the final quasar samples for redshifts $1.4<z<2.2$ from eBOSS and the Dark Energy Spectroscopic Instrument (DESI) are expected to provide measurements of the isotropic BAO scale to $\sim7\%$ and $\sim3\%$ precision, respectively, at $z\simeq1.6$.

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.

Low redshift measurements of Baryon Acoustic Oscillations (BAO) test the late time evolution of the Universe and are a vital probe of Dark Energy. Over the past decade both the 6-degree Field Galaxy Survey (6dFGS) and Sloan Digital Sky Survey (SDSS) have provided important distance constraints at $z < 0.3$. In this paper we re-evaluate the cosmological information from the BAO detection in 6dFGS making use of HOD populated COLA mocks for a robust covariance matrix and taking advantage of the now commonly implemented technique of density field reconstruction. For the 6dFGS data, we find consistency with the previous analysis, and obtain an isotropic volume averaged distance measurement of $D_{V}(z_{\mathrm{eff}}=0.097) = 372\pm17(r_{s}/r_{s}^{\mathrm{fid}})\,\mathrm{Mpc}$, which has a non-Gaussian likelihood outside the $1\sigma$ region. We combine our measurement from both the post-reconstruction clustering of 6dFGS and SDSS MGS offering the most robust constraint to date in this redshift regime, $D_{V}(z_{\mathrm{eff}}=0.122)=539\pm17(r_{s}/r^{\mathrm{fid}}_{s})\,\mathrm{Mpc}$. These measurements are consistent with standard $\Lambda\mathrm{CDM}$ and after fixing the standard ruler using a Planck prior on $\Omega_{m}h^{2}$, the joint analysis gives $H_{0}=64.0\pm3.5\,\mathrm{kms}^{-1}\mathrm{Mpc}^{-1}$. In the near future both the Taipan Galaxy Survey and the Dark Energy Spectroscopic Instrument (DESI) will improve this measurement to $1\%$ at low redshift.

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.

We present a measurement of the anisotropic and isotropic Baryon Acoustic Oscillations (BAO) from the extended Baryon Oscillation Spectroscopic Survey Data Release 14 quasar sample with optimal redshift weights. Applying the redshift weights improves the constraint on the BAO dilation parameter $\alpha(z_{\rm eff})$ by 17\%. We reconstruct the evolution history of the BAO distance indicators in the redshift range of $0.8<z<2.2$. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

We present the Data Release 14 Quasar catalog (DR14Q) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). This catalog includes all SDSS-IV/eBOSS objects that were spectroscopically targeted as quasar candidates and that are confirmed as quasars via a new automated procedure combined with a partial visual inspection of spectra, have luminosities $M_{\rm i} \left[ z=2 \right] < -20.5$ (in a $\Lambda$CDM cosmology with $H_0 = 70 \ {\rm km \ s^{-1} \ Mpc ^{-1}}$, $\Omega_{\rm M} = 0.3$, and $\Omega_{\rm \Lambda} = 0.7$), and either display at least one emission line with a full width at half maximum (FWHM) larger than $500 \ {\rm km \ s^{-1}}$ or, if not, have interesting/complex absorption features. The catalog also includes previously spectroscopically-confirmed quasars from SDSS-I, II and III. The catalog contains 526,356 quasars 144,046 are new discoveries since the beginning of SDSS-IV) detected over 9,376 deg$^2$ (2,044 deg$^2$ having new spectroscopic data available) with robust identification and redshift measured by a combination of principal component eigenspectra. The catalog is estimated to have about 0.5% contamination. The catalog identifies 21,877 broad absorption line quasars and lists their characteristics. For each object, the catalog presents SDSS five-band CCD-based photometry with typical accuracy of 0.03 mag. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

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.

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.

We present an anisotropic analysis of Baryon Acoustic Oscillation (BAO) signal from the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample. The sample consists of 147,000 quasars distributed over a redshift range of $0.8 < z < 2.2$. We apply the redshift weights technique to the clustering of quasars in this sample and achieve a 4.6 per cent measurement of the angular distance measurement $D_M$ at $z = 2.2$ and Hubble parameter $H$ at $z=0.8$. We parameterize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. The coefficients of this expansion are used to reconstruct the distance-redshift relation and obtain distance and Hubble parameter measurements at all redshifts within the redshift range of the sample. Reporting the result at two characteristic redshifts, we determine $D_M(z=1) = 3405\pm305 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=1) = 120.7\pm 7.3 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$ and $D_M(z=2) = 5325\pm249 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=2) = 189.9\pm 32.9 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$. These measurements are highly correlated. We assess the outlook of BAO analysis from the final quasar sample by testing the method on a set of mocks that mimic the noise level in the final sample. We demonstrate on these mocks that redshift weighting shrinks the measurement error by over 25 per cent on average. We conclude redshift weighting can bring us closer to the cosmological goal of the final quasar sample.

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.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 sample includes 80,118 Luminous Red Galaxies. By combining these galaxies with the high-redshift tail of the BOSS galaxy sample, we form a sample of LRGs at an effective redshift $z=0.72$, covering an effective volume of 0.9~Gpc$^3$. We introduce new techniques to account for spurious fluctuations caused by targeting and by redshift failures which were validated on a set of mock catalogs. This analysis is sufficient to provide a $2.6$\% measurement of spherically averaged BAO, $D_V(z=0.72) = 2353^{+63}_{-61} (r_d/r_{d,\rm{fid}}) h^{-1}$Mpc, at 2.8$\sigma$ of significance. Together with the recent quasar-based BAO measurement at $z=1.5$, and forthcoming Emission Line Galaxy-based measurements, this measurement demonstrates that eBOSS is fulfilling its remit of extending the range of redshifts covered by such measurements, laying the ground work for forthcoming surveys such as the Dark Energy Spectroscopic Survey and Euclid.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is one of the first of a new generation of galaxy redshift surveys that will cover a large range in redshift with sufficient resolution to measure the baryon acoustic oscillations (BAO) signal. For surveys covering a large redshift range we can no longer ignore cosmological evolution, meaning that either the redshift shells analysed have to be significantly narrower than the survey, or we have to allow for the averaging over evolving quantities. Both of these have the potential to remove signal: analysing small volumes increases the size of the Fourier window function, reducing the large-scale information, while averaging over evolving quantities can, if not performed carefully, remove differential information. It will be important to measure cosmological evolution from these surveys to explore and discriminate between models. We apply a method to optimally extract this differential information to mock catalogues designed to mimic the eBOSS quasar sample. By applying a set of weights to extract redshift space distortion measurements as a function of redshift, we demonstrate an analysis that does not invoke the problems discussed above. We show that our estimator gives unbiased constraints.

We determine an optimized clustering statistic to be used for galaxy samples with significant redshift uncertainty, such as those that rely on photometric redshifts. To do so, we study the baryon acoustic oscillation (BAO) information content as a function of the orientation of galaxy clustering modes with respect to their angle to the line-of-sight (LOS). The clustering along the LOS, as observed in a redshift-space with significant redshift uncertainty, has contributions from clustering modes with a range of orientations with respect to the true LOS. For redshift uncertainty $\sigma_z \geq 0.02(1+z)$ we find that while the BAO information is confined to transverse clustering modes in the true space, it is spread nearly evenly in the observed space. Thus, measuring clustering in terms of the projected separation (regardless of the LOS) is an efficient and nearly lossless compression of the signal for $\sigma_z \geq 0.02(1+z)$. For reduced redshift uncertainty, a more careful consideration is required. We then use more than 1700 realizations (combining two separate sets) of galaxy simulations mimicking the Dark Energy Survey Year 1 sample to validate our analytic results and optimized analysis procedure. We find that using the correlation function binned in projected separation, we can achieve uncertainties that are within 10 per cent of those predicted by Fisher matrix forecasts. We predict that DES Y1 should achieve a 5 per cent distance measurement using our optimized methods. We expect the results presented here to be important for any future BAO measurements made using photometric redshift data.

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.

We present a measurement of galaxy-galaxy lensing around a magnitude-limited ($i_{AB} < 22.5$) sample of galaxies from the Dark Energy Survey Science Verification (DES-SV) data. We split these lenses into three photometric-redshift bins from 0.2 to 0.8, and determine the product of the galaxy bias $b$ and cross-correlation coefficient between the galaxy and dark matter overdensity fields $r$ in each bin, using scales above 4 Mpc/$h$ comoving, where we find the linear bias model to be valid given our current uncertainties. We compare our galaxy bias results from galaxy-galaxy lensing with those obtained from galaxy clustering (Crocce et al. 2016) and CMB lensing (Giannantonio et al. 2016) for the same sample of galaxies, and find our measurements to be in good agreement with those in Crocce et al. (2016), while, in the lowest redshift bin ($z\sim0.3$), they show some tension with the findings in Giannantonio et al. (2016). We measure $b\cdot r$ to be $0.87\pm 0.11$, $1.12 \pm 0.16$ and $1.24\pm 0.23$, respectively for the three redshift bins of width $\Delta z = 0.2$ in the range $0.2<z <0.8$, defined with the photometric-redshift algorithm BPZ. Using a different code to split the lens sample, TPZ, leads to changes in the measured biases at the 10-20\% level, but it does not alter the main conclusion of this work: when comparing with Crocce et al. (2016) we do not find strong evidence for a cross-correlation parameter significantly below one in this galaxy sample, except possibly at the lowest redshift bin ($z\sim 0.3$), where we find $r = 0.71 \pm 0.11$ when using TPZ, and $0.83 \pm 0.12$ with BPZ.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We perform a detailed analysis of the covariance matrix of the spherically averaged galaxy power spectrum and present a new, practical method for estimating this within an arbitrary survey without the need for running mock galaxy simulations that cover the full survey volume. The method uses theoretical arguments to modify the covariance matrix measured from a set of small-volume cubic galaxy simulations, which are computationally cheap to produce compared to larger simulations and match the measured small-scale galaxy clustering more accurately than is possible using theoretical modelling. We include prescriptions to analytically account for the window function of the survey, which convolves the measured covariance matrix in a non-trivial way. We also present a new method to include the effects of supersample covariance and modes outside the small simulation volume which requires no additional simulations and still allows us to scale the covariance matrix. As validation, we compare the covariance matrix estimated using our new method to that from a brute force calculation using 500 simulations originally created for analysis of the Sloan Digital Sky Survey Main Galaxy Sample (SDSS-MGS). We find excellent agreement on all scales of interest for large scale structure analysis, including those dominated by the effects of the survey window, and on scales where theoretical models of the clustering normally break-down, but the new method produces a covariance matrix with significantly better signal-to-noise. Although only formally correct in real-space, we also discuss how our method can be extended to incorporate the effects of Redshift Space Distortions.

In order to be efficient, spectroscopic galaxy redshift surveys do not obtain redshifts for all galaxies in the population targeted. The missing galaxies are often clustered, commonly leading to a lower proportion of successful observations in dense regions. One example is the close-pair issue for SDSS spectroscopic galaxy surveys, which have a deficit of pairs of observed galaxies with angular separation closer than the hardware limit on placing neighbouring fibers. Spatially clustered missing observations will exist in the next generations of surveys. Various schemes have previously been suggested to mitigate these effects, but none works for all situations. We argue that the solution is to link the missing galaxies to those observed with statistically equivalent clustering properties, and that the best way to do this is to rerun the targeting algorithm, varying the angular position of the observations. Provided that every pair has a non-zero probability of being observed in one realisation of the algorithm, then a pair-upweighting scheme linking targets to successful observations, can correct these issues. We present such a scheme, and demonstrate its validity using realisations of an idealised simple survey strategy.

Three dimensional galaxy clustering measurements provide a wealth of cosmological information. However, obtaining spectra of galaxies is expensive, and surveys often only measure redshifts for a subsample of a target galaxy population. Provided that the spectroscopic data is representative, we argue that angular pair upweighting should be used in these situations to improve the 3D clustering measurements. We present a toy model showing mathematically how such a weighting can improve measurements, and provide a practical example of its application using mocks created for the Baryon Oscillation Spectroscopic Survey (BOSS). Our analysis of mocks suggests that, if an angular clustering measurement is available over twice the area covered spectroscopically, weighting gives a $\sim$10-20% reduction of the variance of the monopole correlation function on the BAO scale.

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.

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.

Previous studies using numerical simulations have demonstrated that the shape of the cosmic web can be described by studying the Lagrangian displacement field. We extend these analyses, showing that it is now possible to perform a Lagrangian description of cosmic structure in the nearby Universe based on large-scale structure observations. Building upon recent Bayesian large-scale inference of initial conditions, we present a cosmographic analysis of the dark matter distribution and its evolution, referred to as the dark matter phase-space sheet, in the nearby universe as probed by the Sloan Digital Sky Survey main galaxy sample. We consider its stretchings and foldings using a tetrahedral tessellation of the Lagrangian lattice. The method provides extremely accurate estimates of nearby density and velocity fields, even in regions of low galaxy density. It also measures the number of matter streams, and the deformation and parity reversals of fluid elements, which were previously thought inaccessible using observations. We illustrate the approach by showing the phase-space structure of known objects of the nearby Universe such as the Sloan Great Wall, the Coma cluster and the Bo\"otes void. We dissect cosmic structures into four distinct components (voids, sheets, filaments, and clusters), using the Lagrangian classifiers DIVA, ORIGAMI, and a new scheme which we introduce and call LICH. Because these classifiers use information other than the sheer local density, identified structures explicitly carry physical information about their formation history. Accessing the phase-space structure of dark matter in galaxy surveys opens the way for new confrontations of observational data and theoretical models. We have made our data products publicly available.

We use the final catalogue of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to measure the power spectrum of the galaxy distribution at high redshift, presenting results that extend beyond $z=1$ for the first time. We apply an FFT technique to four independent sub-volumes comprising a total of $51,728$ galaxies at $0.6<z<1.1$ (out of the nearly $90,000$ included in the whole survey). We concentrate here on the shape of the direction-averaged power spectrum in redshift space, explaining the level of modelling of redshift-space anisotropies and the anisotropic survey window function that are needed to deduce this in a robust fashion. We then use covariance matrices derived from a large ensemble of mock datasets in order to fit the spectral data. The results are well matched by a standard $\Lambda$CDM model, with density parameter $\Omega_M h =\smash{0.227^{+0.063}_{-0.050}}$ and baryon fraction $\smash{f_B=\Omega_B/\Omega_M=0.220^{+0.058}_{-0.072}}$. These inferences from the high-$z$ galaxy distribution are consistent with results from local galaxy surveys, and also with the Cosmic Microwave Background. Thus the $\Lambda$CDM model gives a good match to cosmic structure at all redshifts so far accessible to observational study.

We study the first year of the eBOSS quasar sample in the redshift range $0.9<z<2.2$ which includes 68,772 homogeneously selected quasars. We show that the main source of systematics in the evaluation of the correlation function arises from inhomogeneities in the quasar target selection, particularly related to the extinction and depth of the imaging data used for targeting. We propose a weighting scheme that mitigates these systematics. We measure the quasar correlation function and provide the most accurate measurement to date of the quasar bias in this redshift range, $b_Q = 2.45 \pm 0.05$ at $\bar z=1.55$, together with its evolution with redshift. We use this information to determine the minimum mass of the halo hosting the quasars and the characteristic halo mass, which we find to be both independent of redshift within statistical error. Using a recently-measured quasar-luminosity-function we also determine the quasar duty cycle. The size of this first year sample is insufficient to detect any luminosity dependence to quasar clustering and this issue should be further studied with the final $\sim$500,000 eBOSS quasar sample.

We investigate the cosmological implications of studying galaxy clustering using a tomographic approach applied to the final BOSS DR12 galaxy sample, including both auto- and cross-correlation functions between redshift shells. We model the signal of the full shape of the angular correlation function, $\omega(\theta)$, in redshift bins using state-of-the-art modelling of non-linearities, bias and redshift-space distortions. We present results on the redshift evolution of the linear bias of BOSS galaxies, which cannot be obtained with traditional methods for galaxy-clustering analysis. We also obtain constraints on cosmological parameters, combining this tomographic analysis with measurements of the cosmic microwave background (CMB) and type Ia supernova (SNIa). We explore a number of cosmological models, including the standard $\Lambda$CDM model and its most interesting extensions, such as deviations from $w_\rm{DE} = -1$, non-minimal neutrino masses, spatial curvature and deviations from general relativity using the growth-index $\gamma$ parametrisation. These results are, in general, comparable to the most precise present-day constraints on cosmological parameters, and show very good agreement with the standard model. In particular, combining CMB, $\omega(\theta)$ and SNIa, we find a value of $w_\rm{DE}$ consistent with $-1$ to a precision better than 5\% when it is assumed to be constant in time, and better than 6\% when we also allow for a spatially-curved Universe.

We describe the algorithm used to select the Emission Line Galaxy (ELG) sample at $z \sim 0.85$ for the extended Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey IV, using photometric data from the DECam Legacy Survey. Our selection is based on a selection box in the $g-r$ vs. $r-z$ colour-colour space and a cut on the $g$-band magnitude, to favour galaxies in the desired redshift range with strong [OII] emission. It provides a target density of 200 deg$^{-2}$ on the North Galactic Cap (NGC) and of 240 deg$^{-2}$ on the South Galactic Cap (SGC), where we use a larger selection box because of deeper imaging. We demonstrate that this selection passes the eBOSS requirements in terms of homogeneity. About 50,000 ELGs have been observed since the observations have started in 2016, September. These roughly match the expected redshift distribution, though the measured efficiency is slightly lower than expected. The efficiency can be increased by enlarging the redshift range and with incoming pipeline improvement. The cosmological forecast based on these first data predict $\sigma_{D_V}/D_V = 0.023$, in agreement with previous forecasts. Lastly, we present the stellar population properties of the ELG SGC sample. Once observations are completed, this sample will be suited to provide a cosmological analysis at $z \sim 0.85$, and will pave the way for the next decade of massive spectroscopic cosmological surveys, which heavily rely on ELGs. The target catalogue over the SGC will be released along with DR14.

We present a measurement of the linear growth rate of structure, \textit{f} from the Sloan Digital Sky Survey III (SDSS III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) using Convolution Lagrangian Perturbation Theory (CLPT) with Gaussian Streaming Redshift-Space Distortions (GSRSD) to model the two point statistics of BOSS galaxies in DR12. The BOSS-DR12 dataset includes 1,198,006 massive galaxies spread over the redshift range $0.2 < z < 0.75$. These galaxy samples are categorized in three redshift bins. Using CLPT-GSRSD in our analysis of the combined sample of the three redshift bins, we report measurements of $f \sigma_8$ for the three redshift bins. We find $f \sigma_8 = 0.430 \pm 0.054$ at $z_{\rm eff} = 0.38$, $f \sigma_8 = 0.452 \pm 0.057$ at $z_{\rm eff} = 0.51$ and $f \sigma_8 = 0.457 \pm 0.052$ at $z_{\rm eff} = 0.61$. Our results are consistent with the predictions of Planck $\Lambda$CDM-GR. Our constraints on the growth rates of structure in the Universe at different redshifts serve as a useful probe, which can help distinguish between a model of the Universe based on dark energy and models based on modified theories of gravity. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. 2016 to produce the final cosmological constraints from BOSS.

We perform a tomographic baryon acoustic oscillations analysis using the two-point galaxy correlation function measured from the combined sample of BOSS DR12, which covers the redshift range of $0.2<z<0.75$. Splitting the sample into multiple overlapping redshift slices to extract the redshift information of galaxy clustering, we obtain a measurement of $D_A(z)/r_d$ and $H(z)r_d$ at nine effective redshifts with the full covariance matrix calibrated using MultiDark-Patchy mock catalogues. Using the reconstructed galaxy catalogues, we obtain the precision of $1.3\%-2.2\%$ for $D_A(z)/r_d$ and $2.1\%-6.0\%$ for $H(z)r_d$. To quantify the gain from the tomographic information, we compare the constraints on the cosmological parameters using our 9-bin BAO measurements, the consensus 3-bin BAO and RSD measurements at three effective redshifts in \citet{Alam2016}, and the non-tomographic (1-bin) BAO measurement at a single effective redshift. Comparing the 9-bin with 1-bin constraint result, it can improve the dark energy Figure of Merit by a factor of 1.24 for the Chevallier-Polarski-Linder parametrisation for equation of state parameter $w_{\rm DE}$. The errors of $w_0$ and $w_a$ from 9-bin constraints are slightly improved when compared to the 3-bin constraint result.

Galaxy clustering data from current and upcoming large scale structure surveys can provide strong constraints on primordial non-Gaussianity through the scale-dependent halo bias. To fully exploit the information from galaxy surveys, optimal analysis methods need to be developed and applied to the data. Since the halo bias is sensitive to local non-Gaussianity predominately at large scales, the volume of a given survey is crucial. Consequently, for such analyses we do not want to split into redshift bins, which would lead to information loss due to edge effects, but instead analyse the full sample. We present an optimal technique to directly constrain local non-Gaussianity parametrised by $f_\mathrm{NL}^\mathrm{loc}$, from galaxy clustering by applying redshift weights to the galaxies. We derive a set of weights to optimally measure the amplitude of local non-Gaussianity, $f_\mathrm{NL}^\mathrm{loc}$, discuss the redshift weighted power spectrum estimators, outline the implementation procedure and test our weighting scheme against Lognormal catalogs for two different surveys: the quasar sample of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the emission line galaxy sample of the Dark Energy Spectroscopic Instrument (DESI) survey. We find an improvement of 30 percent for eBOSS and 6 percent for DESI compared to the standard Feldman, Kaiser $\&$ Peacock weights, although these predictions are sensitive to the bias model assumed.

We consider the application of relative self-calibration using overlap regions to spectroscopic galaxy surveys that use slit-less spectroscopy. This method is based on that developed for the SDSS by Padmanabhan at al. (2008) in that we consider jointly fitting and marginalising over calibrator brightness, rather than treating these as free parameters. However, we separate the calibration of the detector-to-detector from the full-focal-plane exposure-to-exposure calibration. To demonstrate how the calibration procedure will work, we simulate the procedure for a potential implementation of the spectroscopic component of the wide Euclid survey. We study the change of coverage and the determination of relative multiplicative errors in flux measurements for different dithering configurations. We use the new method to study the case where the flat-field across each exposure or detector is measured precisely and only exposure-to-exposure or detector-to-detector variation in the flux error remains. We consider several base dither patterns and find that they strongly influence the ability to calibrate, using this methodology. To enable self-calibration, it is important that the survey strategy connects different observations with at least a minimum amount of overlap, and we propose an "S"-pattern for dithering that fulfills this requirement. The final survey strategy adopted by Euclid will have to optimise for a number of different science goals and requirements. The large-scale calibration of the spectroscopic galaxy survey is clearly cosmologically crucial, but is not the only one.

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We perform a tomographic baryon acoustic oscillations (BAO) analysis using the monopole, quadrupole and hexadecapole of the redshift-space galaxy power spectrum measured from the pre-reconstructed combined galaxy sample of the completed Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR)12 covering the redshift range of $0.20<z<0.75$. By allowing for overlap between neighbouring redshift slices, we successfully obtained the isotropic and anisotropic BAO distance measurements within nine redshift slices to a precision of $1.5\%-3.4\%$ for $D_V/r_d$, $1.8\% -4.2\%$ for $D_A/r_d$ and $3.7\% - 7.5\%$ for $H \ r_d$, depending on effective redshifts. We provide our BAO measurement of $D_A/r_d$ and $H \ r_d$ with the full covariance matrix, which can be used for cosmological implications. Our measurements are consistent with those presented in \citet{Acacia}, in which the BAO distances are measured at three effective redshifts. We constrain dark energy parameters using our measurements, and find an improvement of the Figure-of-Merit of dark energy in general due to the temporal BAO information resolved. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

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.

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 present the full public data release (PDR-2) of the VIMOS Public Extragalactic Redshift Survey (VIPERS), performed at the ESO VLT. We release redshifts, spectra, CFHTLS magnitudes and ancillary information (as masks and weights) for a complete sample of 86,775 galaxies (plus 4,732 other objects, including stars and serendipitous galaxies); we also include their full photometrically-selected parent catalogue. The sample is magnitude limited to i_AB < 22.5, with an additional colour-colour pre-selection devised as to exclude galaxies at z<0.5. This practically doubles the effective sampling of the VIMOS spectrograph over the range 0.5<z<1.2 (reaching 47% on average), yielding a final median local galaxy density close to 5 10^{-3} h^3 Mpc}^{-3}. The total area spanned by the final data set is ~ 23.5 deg^2, corresponding to 288 VIMOS fields with marginal overlaps, split over two regions within the CFHTLS-Wide W1 and W4 equatorial fields (at R.A. ~2 and ~22 hours, respectively). Spectra were observed at a resolution R=220, covering a wavelength range 5500-9500 Angstrom. Data reduction and redshift measurements were performed through a fully automated pipeline; all redshift determinations were then visually validated and assigned a quality flag. Measurements with a quality flag >= 2 are shown to have a confidence level of 96% or larger and make up 88% of all measured galaxy redshifts (76,552 out of 86,775), constituting the VIPERS prime catalogue for statistical investigations. For this sample the rms redshift error, estimated using repeated measurements of about 3,000 galaxies, is found to be sigma_z = 0.00054(1+z). All data are available at http://vipers.inaf.it and on the ESO Archive.

We present two wide-field catalogs of photometrically-selected emission line galaxies (ELGs) at z=0.8 covering about 2800 deg^2 over the south galactic cap. The catalogs were obtained using a Fisher discriminant technique described in a companion paper. The two catalogs differ by the imaging used to define the Fisher discriminant: the first catalog includes imaging from the Sloan Digital Sky Survey and the Wide-Field Infrared Survey Explorer, the second also includes information from the South Galactic Cap U-band Sky Survey (SCUSS). Containing respectively 560,045 and 615,601 objects, they represent the largest ELG catalogs available today and were designed for the ELG programme of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We study potential sources of systematic variation in the angular distribution of the selected ELGs due to fluctuations of the observational parameters. We model the influence of the observational parameters using a multivariate regression and implement a weighting scheme that allows effective removal of all of the systematic errors induced by the observational parameters. We show that fluctuations in the imaging zero-points of the photometric bands have minor impact on the angular distribution of objects in our catalogs. We compute the angular clustering of both catalogs and show that our weighting procedure effectively removes spurious clustering on large scales. We fit a model to the small scale angular clustering, showing that the selections have similar biases of 1.35/D_a(z) and 1.28/D_a(z). Both catalogs are publicly available.

We extract cosmological information from the anisotropic power spectrum measurements from the recently completed Baryon Oscillation Spectroscopic Survey (BOSS), extending the concept of clustering wedges to Fourier space. Making use of new FFT-based estimators, we measure the power spectrum clustering wedges of the BOSS sample by filtering out the information of Legendre multipoles l > 4. Our modelling of these measurements is based on novel approaches to describe non-linear evolution, bias, and redshift-space distortions, which we test using synthetic catalogues based on large-volume N-body simulations. We are able to include smaller scales than in previous analyses, resulting in tighter cosmological constraints. Using three overlapping redshift bins, we measure the angular diameter distance, the Hubble parameter, and the cosmic growth rate, and explore the cosmological implications of our full shape clustering measurements in combination with CMB and SN Ia data. Assuming a {\Lambda}CDM cosmology, we constrain the matter density to {\Omega}_m = 0.311 -0.010 +0.009 and the Hubble parameter to H_0 = 67.6 -0.6 +0.7 km s^-1 Mpc^-1, at a confidence level (CL) of 68 per cent. We also allow for non-standard dark energy models and modifications of the growth rate, finding good agreement with the {\Lambda}CDM paradigm. For example, we constrain the equation-of-state parameter to w = -1.019 -0.039 +0.048. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. 2016 to produce the final cosmological constraints from BOSS.

The low statistical errors on cosmological parameters promised by future galaxy surveys will only be realised with the development of new, fast, analysis methods that reduce potential systematic problems to low levels. We present an efficient method for measuring the evolution of the growth of structure using Redshift Space Distortions (RSD), that removes the need to make measurements in redshift shells. We provide sets of galaxy-weights that cover a wide range in redshift, but are optimised to provide differential information about cosmological evolution. These are derived to optimally measure the coefficients of a parameterisation of the redshift-dependent matter density, which provides a framework to measure deviations from the concordance $\Lambda$CDM cosmology, allowing for deviations in both geometric and/or growth. We test the robustness of the weights by comparing with alternative schemes and investigate the impact of galaxy bias. We extend the results to measure the combined anisotropic Baryon Acoustic Oscillation (BAO) and RSD signals.

We measure and analyse the bispectrum of the final, Data Release 12, galaxy sample provided by the Baryon Oscillation Spectroscopic Survey, splitting by selection algorithm into LOWZ and CMASS galaxies. The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. Combining the power spectrum, measured relative to the line-of-sight, with the spherically averaged bispectrum, we are able to constrain the product of the growth of structure parameter, $f$, and the amplitude of dark matter density fluctuations, $\sigma_8$, along with the geometric Alcock-Paczynski parameters, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch, $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon, $D_A(z)/r_s(z_d)$. After combining pre-reconstruction RSD analyses of the power spectrum monopole, quadrupole and bispectrum monopole; with post-reconstruction analysis of the BAO power spectrum monopole and quadrupole, we find $f(z_{\rm LOWZ})\sigma_8(z_{\rm LOWZ})=0.427\pm 0.056$, $D_A(z_{\rm LOWZ})/r_s(z_d)=6.60 \pm 0.13$, $H(z_{\rm LOWZ})r_s(z_d)=(11.55\pm 0.38)10^3\,{\rm kms}^{-1}$ for the LOWZ sample, and $f(z_{\rm CMASS})\sigma_8(z_{\rm CMASS})=0.426\pm 0.029$, $D_A(z_{\rm CMASS})/r_s(z_d)=9.39 \pm 0.10$, $H(z_{\rm CMASS})r_s(z_d)=(14.02\pm 0.22)10^3\,{\rm kms}^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 and DR12 measurements. Combining our dataset with {\it Planck15} we perform a null test of General Relativity (GR) through the $\gamma$-parametrisation finding $\gamma=0.733^{+0.068}_{-0.069}$, which is $\sim2.7\sigma$ away from the GR predictions.

We present baryon acoustic oscillation (BAO) scale measurements determined from the clustering of 1.2 million massive galaxies with redshifts 0.2 < z < 0.75 distributed over 9300 square degrees, as quantified by their redshift-space correlation function. In order to facilitate these measurements, we define, describe, and motivate the selection function for galaxies in the final data release (DR12) of the SDSS III Baryon Oscillation Spectroscopic Survey (BOSS). This includes the observational footprint, masks for image quality and Galactic extinction, and weights to account for density relationships intrinsic to the imaging and spectroscopic portions of the survey. We simulate the observed systematic trends in mock galaxy samples and demonstrate that they impart no bias on baryon acoustic oscillation (BAO) scale measurements and have a minor impact on the recovered statistical uncertainty. We measure transverse and radial BAO distance measurements in 0.2 < z < 0.5, 0.5 < z < 0.75, and (overlapping) 0.4 < z < 0.6 redshift bins. In each redshift bin, we obtain a precision that is 2.7 per cent or better on the radial distance and 1.6 per cent or better on the transverse distance. The combination of the redshift bins represents 1.8 per cent precision on the radial distance and 1.1 per cent precision on the transverse distance. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

We develop a model for the redshift-space correlation function, valid for both dark matter particles and halos on scales $>5\,h^{-1}$Mpc. In its simplest formulation, the model requires the knowledge of the first three moments of the line-of-sight pairwise velocity distribution plus two well-defined dimensionless parameters. The model is obtained by extending the Gaussian-Gaussianity prescription for the velocity distribution, developed in a previous paper, to a more general concept allowing for local skewness, which is required to match simulations. We compare the model with the well known Gaussian streaming model and the more recent Edgeworth streaming model. Using N-body simulations as a reference, we show that our model gives a precise description of the redshift-space clustering over a wider range of scales. We do not discuss the theoretical prescription for the evaluation of the velocity moments, leaving this topic to further investigation.

We assess and develop techniques to remove contaminants when calculating the 3D galaxy power spectrum. We separate the process into three separate stages: (i) removing the contaminant signal, (ii) estimating the uncontaminated cosmological power spectrum, (iii) debiasing the resulting estimates. For (i), we show that removing the best-fit contaminant mode subtraction), and setting the contaminated components of the covariance to be infinite (mode deprojection) are mathematically equivalent. For (ii), performing a Quadratic Maximum Likelihood (QML) estimate after mode deprojection gives an optimal unbiased solution, although it requires the manipulation of large $N_{\rm mode}^2$ matrices ($N_{\rm mode}$ being the total number of modes)}, which is unfeasible for recent 3D galaxy surveys. Measuring a binned average of the modes for (ii) as proposed by \citet*[FKP]{Feldman} is faster and simpler, but is sub-optimal and gives rise to a biased solution. We present a method to debias the resulting FKP measurements that does not require any large matrix calculations. We argue that the sub-optimality of the FKP estimator compared with the QML estimator, caused by contaminants is less severe than that commonly ignored due to the survey window.

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).

We present the large-scale 3-point correlation function (3PCF) of the SDSS DR12 CMASS sample of $777,202$ Luminous Red Galaxies, the largest-ever sample used for a 3PCF or bispectrum measurement. We make the first high-significance ($4.5\sigma$) detection of Baryon Acoustic Oscillations (BAO) in the 3PCF. Using these acoustic features in the 3PCF as a standard ruler, we measure the distance to $z=0.57$ to $1.7\%$ precision (statistical plus systematic). We find $D_{\rm V}= 2024\pm29\;{\rm Mpc\;(stat)}\pm20\;{\rm Mpc\;(sys)}$ for our fiducial cosmology (consistent with Planck 2015) and bias model. This measurement extends the use of the BAO technique from the 2-point correlation function (2PCF) and power spectrum to the 3PCF and opens an avenue for deriving additional cosmological distance information from future large-scale structure redshift surveys such as DESI. Our measured distance scale from the 3PCF is fairly independent from that derived from the pre-reconstruction 2PCF and is equivalent to increasing the length of BOSS by roughly 10\%; reconstruction appears to lower the independence of the distance measurements. Fitting a model including tidal tensor bias yields a moderate significance ($2.6\sigma)$ detection of this bias with a value in agreement with the prediction from local Lagrangian biasing.

We search for a galaxy clustering bias due to a modulation of galaxy number with the baryon-dark matter relative velocity resulting from recombination-era physics. We find no detected signal and place the constraint $b_v < 0.01$ on the relative velocity bias for the CMASS galaxies. This bias is an important potential systematic of Baryon Acoustic Oscillation (BAO) method measurements of the cosmic distance scale using the 2-point clustering. Our limit on the relative velocity bias indicates a systematic shift of no more than $0.3\%$ rms in the distance scale inferred from the BAO feature in the BOSS 2-point clustering, well below the $1\%$ statistical error of this measurement. This constraint is the most stringent currently available and has important implications for the ability of upcoming large-scale structure surveys such as DESI to self-protect against the relative velocity as a possible systematic.

We develop a new methodology called double-probe analysis with the aim of minimizing informative priors in the estimation of cosmological parameters. We extract the dark-energy-model-independent cosmological constraints from the joint data sets of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy sample and Planck cosmic microwave background (CMB) measurement. We measure the mean values and covariance matrix of $\{R$, $l_a$, $\Omega_b h^2$, $n_s$, $log(A_s)$, $\Omega_k$, $H(z)$, $D_A(z)$, $f(z)\sigma_8(z)\}$, which give an efficient summary of Planck data and 2-point statistics from BOSS galaxy sample, where $R=\sqrt{\Omega_m H_0^2}\,r(z_*)$, and $l_a=\pi r(z_*)/r_s(z_*)$, $z_*$ is the redshift at the last scattering surface, and $r(z_*)$ and $r_s(z_*)$ denote our comoving distance to $z_*$ and sound horizon at $z_*$ respectively. The advantage of this method is that we do not need to put informative priors on the cosmological parameters that galaxy clustering is not able to constrain well, i.e. $\Omega_b h^2$ and $n_s$. Using our double-probe results, we obtain $\Omega_m=0.304\pm0.009$, $H_0=68.2\pm0.7$, and $\sigma_8=0.806\pm0.014$ assuming $\Lambda$CDM; and $\Omega_k=0.002\pm0.003$ and $w=-1.00\pm0.07$ assuming o$w$CDM. The results show no tension with the flat $\Lambda$CDM cosmological paradigm. By comparing with the full-likelihood analyses with fixed dark energy models, we demonstrate that the double-probe method provides robust cosmological parameter constraints which can be conveniently used to study dark energy models. We extend our study to measure the sum of neutrino mass and obtain $\Sigma m_\nu<0.10/0.22$ (68\%/95\%) assuming $\Lambda$CDM and $\Sigma m_\nu<0.26/0.52$ (68\%/95\%) assuming $w$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

We analyse the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS and LOWZ galaxy sample to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalised growth rate $f(z)\sigma_8(z)$, and the physical matter density $\Omega_mh^2$. We adopt wide and flat priors on all model parameters in order to ensure the results are those of a `single-probe' galaxy clustering analysis. We also marginalise over three nuisance terms that account for potential observational systematics affecting the measured monopole. However, such Monte Carlo Markov Chain analysis is computationally expensive for advanced theoretical models, thus we develop a new methodology to speed up our analysis. We obtain $\{D_A(z)r_{s,fid}/r_s$Mpc, $H(z)r_s/r_{s,fid}$kms$^{-1}$Mpc$^{-1}$, $f(z)\sigma_8(z)$, $\Omega_m h^2\}$ = $\{956\pm28$ , $75.0\pm4.0$ , $0.397 \pm 0.073$, $0.143\pm0.017\}$ at $z=0.32$ and $\{1421\pm23$, $96.7\pm2.7$ , $0.497 \pm 0.058$, $0.137\pm0.015\}$ at $z=0.59$ where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$Mpc for the fiducial cosmology in this study. In addition, we divide the galaxy sample into four redshift bins to increase the sensitivity of redshift evolution. However, we do not find improvements in terms of constraining dark energy model parameters. Combining our measurements with Planck data, we obtain $\Omega_m=0.306\pm0.009$, $H_0=67.9\pm0.7$kms$^{-1}$Mpc$^{-1}$, and $\sigma_8=0.815\pm0.009$ assuming $\Lambda$CDM; $\Omega_k=0.000\pm0.003$ assuming oCDM; $w=-1.01\pm0.06$ assuming $w$CDM; and $w_0=-0.95\pm0.22$ and $w_a=-0.22\pm0.63$ assuming $w_0w_a$CDM. Our results show no tension with the flat $\Lambda$CDM cosmological paradigm. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS.

We investigate the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) sample, which consists of $1\,198\,006$ galaxies in the redshift range $0.2 < z < 0.75$ and a sky coverage of $10\,252\,$deg$^2$. We analyse this dataset in Fourier space, using the power spectrum multipoles to measure Redshift-Space Distortions (RSD) simultaneously with the Alcock-Paczynski (AP) effect and the Baryon Acoustic Oscillation (BAO) scale. We include the power spectrum monopole, quadrupole and hexadecapole in our analysis and compare our measurements with a perturbation theory based model, while properly accounting for the survey window function. To evaluate the reliability of our analysis pipeline we participate in a mock challenge, which resulted in systematic uncertainties significantly smaller than the statistical uncertainties. While the high-redshift constraint on $f\sigma_8$ at $z_{\rm eff}=0.61$ indicates a small ($\sim 1.4\sigma$) deviation from the prediction of the Planck $\Lambda$CDM model, the low-redshift constraint is in good agreement with Planck $\Lambda$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in~\citet{Alam2016} to produce the final cosmological constraints from BOSS.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum.

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in Fourier-space, using the power spectrum monopole and quadrupole. The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z < 0.75$. We divide this dataset into three (overlapping) redshift bins with the effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as $\sim 1000$ MultiDark-Patchy mock catalogues, which mimic the BOSS-DR12 target selection. We apply density field reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can separate the line-of-sight and angular modes, which allows us to constrain the angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately. We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and $2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high ($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and $0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring the Alcock-Paczynski effect. The detection significance of the BAO signal is of the order of $8\sigma$ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within $\Lambda$CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in~\citet{Alam2016} to produce the final cosmological constraints from BOSS.

We explore the cosmological implications of anisotropic clustering measurements in configuration space of the final galaxy samples from Data Release 12 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We implement a new detailed modelling of the effects of non-linearities, galaxy bias and redshift-space distortions that can be used to extract unbiased cosmological information from our measurements for scales $s \gtrsim 20\,h^{-1}{\rm Mpc}$. We combined the galaxy clustering information from BOSS with the latest cosmic microwave background (CMB) observations and Type Ia supernovae samples and found no significant evidence for a deviation from the $\Lambda$CDM cosmological model. In particular, these data sets can constrain the dark energy equation of state parameter to $w_{\rm DE}=-0.996\pm0.042$ when assumed time-independent, the curvature of the Universe to $\Omega_{k}=-0.0007\pm 0.0030$ and the sum of the neutrino masses to $\sum m_{\nu} < 0.25\,{\rm eV}$ at 95 per cent CL. We explore the constraints on the growth rate of cosmic structures assuming $f(z)=\Omega_{\rm m}(z)^\gamma$ and obtain $\gamma = 0.609\pm 0.079$, in good agreement with the predictions of general relativity of $\gamma=0.55$. We compress the information of our clustering measurements into constraints on the parameter combinations $D_{\rm V}(z)/r_{\rm d}$, $F_{\rm AP}(z)$ and $f\sigma_8(z)$ at the effective redshifts of $z=0.38$, $0.51$ and $0.61$ with their respective covariance matrices and find good agreement with the predictions for these parameters obtained from the best-fitting $\Lambda$CDM model to the CMB data from the Planck satellite. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

The cosmological information contained in anisotropic galaxy clustering measurements can often be compressed into a small number of parameters whose posterior distribution is well described by a Gaussian. We present a general methodology to combine these estimates into a single set of consensus constraints that encode the total information of the individual measurements, taking into account the full covariance between the different methods. We illustrate this technique by applying it to combine the results obtained from different clustering analyses, including measurements of the signature of baryon acoustic oscillations (BAO) and redshift-space distortions (RSD), based on a set of mock catalogues of the final SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). Our results show that the region of the parameter space allowed by the consensus constraints is smaller than that of the individual methods, highlighting the importance of performing multiple analyses on galaxy surveys even when the measurements are highly correlated. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The methodology presented here is used in Alam et al. (2016) to produce the final cosmological constraints from BOSS.

We use the first release of the VImos Public Extragalactic Redshift Survey of galaxies (VIPERS) of ~50,000 objects to measure the biasing relation between galaxies and mass in the redshift range z=[0.5,1.1]. We estimate the 1-point distribution function [PDF] of VIPERS galaxies from counts in cells and, assuming a model for the mass PDF, we infer their mean bias relation. The reconstruction of the bias relation from PDFs is performed through a novel method that accounts for Poisson noise, redshift distortions, inhomogeneous sky coverage and other selection effects. With this procedure we constrain galaxy bias and its deviations from linearity down to scales as small as 4 Mpc/h and out to z=1.1. We detect small (~3%) but significant deviations from linear bias. The mean biasing function is close to linear in regions above the mean density. The mean slope of the biasing relation is a proxy to the linear bias parameter. It increases both with luminosity, in agreement with results of previous analyses, and with redshift. However, we detect a strong bias evolution only for z>0.9 in agreement with some, but not all, previous studies. We also detected a significant increase of the bias with the scale, from 4 to 8 Mpc/h, now seen for the first time out to z=1. The amplitude of nonlinearity depends on redshift, luminosity and on scales but no clear trend is detected. Thanks to the large cosmic volume probed by VIPERS we find that the mismatch between the previous estimates of bias at z~1 from zCOSMOS and VVDS-Deep galaxy samples is fully accounted for by cosmic variance. The results of our work confirm the importance of going beyond the over-simplistic linear bias hypothesis showing that non-linearities can be accurately measured through the applications of the appropriate statistical tools to existing datasets like VIPERS.

The Sloan Digital Sky Survey IV extended Baryonic Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) will observe 195,000 emission-line galaxies (ELGs) to measure the Baryonic Acoustic Oscillation standard ruler (BAO) at redshift 0.9. To test different ELG selection algorithms, 9,000 spectra were observed with the SDSS spectrograph as a pilot survey based on data from several imaging surveys. First, using visual inspection and redshift quality flags, we show that the automated spectroscopic redshifts assigned by the pipeline meet the quality requirements for a reliable BAO measurement. We also show the correlations between sky emission, signal-to-noise ratio in the emission lines, and redshift error. Then we provide a detailed description of each target selection algorithm we tested and compare them with the requirements of the eBOSS experiment. As a result, we provide reliable redshift distributions for the different target selection schemes we tested. Finally, we determine an target selection algorithms that is best suited to be applied on DECam photometry because they fulfill the eBOSS survey efficiency requirements.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

[abridged] We present an anisotropic analysis of the baryonic acoustic oscillation (BAO) scale in the twelfth and final data release of the Baryonic Oscillation Spectroscopic Survey (BOSS). We independently analyse the LOWZ and CMASS galaxy samples: the LOWZ sample contains contains 361 762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$; the CMASS sample consists of 777 202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. We extract the BAO peak position from the monopole power spectrum moment, $\alpha_0$, and from the $\mu^2$ moment, $\alpha_2$, where $\mu$ is the cosine of the angle to the line-of-sight. The $\mu^2$-moment provides equivalent information to that available in the quadrupole but is simpler to analyse. After applying a reconstruction algorithm to reduce the BAO suppression by bulk motions, we measure the BAO peak position in the monopole and $\mu^2$-moment, which are related to radial and angular shifts in scale. We report $H(z_{\rm LOWZ})r_s(z_d)=(11.60\pm0.60)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.66\pm0.16$ with a cross-correlation coefficient of $r_{HD_A}=0.41$, for the LOWZ sample; and $H(z_{\rm CMASS})r_s(z_d)=(14.56\pm0.37)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.42\pm0.13$ with a cross-correlation coefficient of $r_{HD_A}=0.47$, for the CMASS sample. We combine these results with the measurements of the BAO peak position in the monopole and quadrupole correlation function of the same dataset \citep[][companion paper]{Cuestaetal2015} and report the consensus values: $H(z_{\rm LOWZ})r_s(z_d)=(11.63\pm0.69)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm LOWZ})/r_s(z_d)=6.67\pm0.15$ with $r_{HD_A}=0.35$ for the LOWZ sample; $H(z_{\rm CMASS})r_s(z_d)=(14.67\pm0.42)\cdot10^3 {\rm km}s^{-1}$ and $D_A(z_{\rm CMASS})/r_s(z_d)=9.47\pm0.12$ with $r_{HD_A}=0.52$ for the CMASS sample.

We measure and analyse the clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) relative to the line-of-sight (LOS), for LOWZ and CMASS galaxy samples drawn from the final Data Release 12 (DR12). The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm lowz}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm cmass}=0.57$. From the power spectrum monopole and quadrupole moments around the LOS, we measure the growth of structure parameter $f$ times the amplitude of dark matter density fluctuations $\sigma_8$ by modeling the Redshift-Space Distortion signal. When the geometrical Alcock-Paczynski effect is also constrained from the same data, we find joint constraints on $f\sigma_8$, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon $D_A(z)/r_s(z_d)$. We find $f(z_{\rm lowz})\sigma_8(z_{\rm lowz})=0.394\pm0.062$, $D_A(z_{\rm lowz})/r_s(z_d)=6.35\pm0.19$, $H(z_{\rm lowz})r_s(z_d)=(11.41\pm 0.56)\,{10^3\rm km}s^{-1}$ for the LOWZ sample, and $f(z_{\rm cmass})\sigma_8(z_{\rm cmass})=0.444\pm0.038$, $D_A(z_{\rm cmass})/r_s(z_d)=9.42\pm0.15$, $H(z_{\rm cmass})r_s(z_d)=(13.92 \pm 0.44)\, {10^3\rm km}s^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 measurements. Assuming the Hubble parameter and angular distance parameter are fixed at fiducial $\Lambda$CDM values, we find $f(z_{\rm lowz})\sigma_8(z_{\rm lowz})=0.485\pm0.044$ and $f(z_{\rm cmass})\sigma_8(z_{\rm cmass})=0.436\pm0.022$ for the LOWZ and CMASS samples, respectively.

Sound waves from the primordial fluctuations of the Universe imprinted in the large-scale structure, called baryon acoustic oscillations (BAOs), can be used as standard rulers to measure the scale of the Universe. These oscillations have already been detected in the distribution of galaxies. Here we propose to measure BAOs from the troughs (minima) of the density field. Based on two sets of accurate mock halo catalogues with and without BAOs in the seed initial conditions, we demonstrate that the BAO signal cannot be obtained from the clustering of classical disjoint voids, but is clearly detected from overlapping voids. The latter represent an estimate of all troughs of the density field. We compute them from the empty circumsphere centers constrained by tetrahedra of galaxies using Delaunay triangulation. Our theoretical models based on an unprecedented large set of detailed simulated void catalogues are remarkably well confirmed by observational data. We use the largest recently publicly available sample of luminous red galaxies from SDSS-III BOSS DR11 to unveil for the first time a >3$\sigma$ BAO detection from voids in observations. Since voids are nearly isotropically expanding regions, their centers represent the most quiet places in the Universe, keeping in mind the cosmos origin and providing a new promising window in the analysis of the cosmological large-scale structure from galaxy surveys.

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

Accurate statistical measurement with large imaging surveys has traditionally required throwing away a sizable fraction of the data. This is because most measurements have have relied on selecting nearly complete samples, where variations in the composition of the galaxy population with seeing, depth, or other survey characteristics are small. We introduce a new measurement method that aims to minimize this wastage, allowing precision measurement for any class of stars or galaxies detectable in an imaging survey. We have implemented our proposal in Balrog, a software package which embeds fake objects in real imaging in order to accurately characterize measurement biases. We demonstrate this technique with an angular clustering measurement using Dark Energy Survey (DES) data. We first show that recovery of our injected galaxies depends on a wide variety of survey characteristics in the same way as the real data. We then construct a flux-limited sample of the faintest galaxies in DES, chosen specifically for their sensitivity to depth and seeing variations. Using the synthetic galaxies as randoms in the standard Landy-Szalay correlation function estimator suppresses the effects of variable survey selection by at least two orders of magnitude. With this correction, our measured angular clustering is found to be in excellent agreement with that of a matched sample drawn from much deeper, higher-resolution space-based Cosmological Evolution Survey (COSMOS) imaging; over angular scales of $0.004^{\circ} < \theta < 0.2^{\circ}$, we find a best-fit scaling amplitude between the DES and COSMOS measurements of $1.00 \pm 0.09$. We expect this methodology to be broadly useful for extending the statistical reach of measurements in a wide variety of coming imaging surveys.

We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range $0.2 < z < 1.2$, a cross-correlation signal is detected at $6 \sigma$ and $4\sigma$ with SPT and Planck respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant ($>$$2 \sigma$) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the $z<1.2$ universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck cosmology, a $1.7\sigma$ deviation.

We reproduce the galaxy clustering catalogue from the SDSS-III Baryon Oscillation Spectroscopic Survey Final Data Release (BOSS DR11&DR12) with high fidelity on all relevant scales in order to allow a robust analysis of baryon acoustic oscillations and redshift space distortions. We have generated (6,000) 12,288 MultiDark PATCHY BOSS (DR11) DR12 light-cones corresponding to an effective volume of $\sim192,000\,[h^{-1}\,{\rm Gpc}]^3$ (the largest ever simulated volume), including cosmic evolution in the redshift range from 0.15 to 0.75. The mocks have been calibrated using a reference galaxy catalogue based on the halo abundance matching modelling of the BOSS DR11&DR12 galaxy clustering data and on the data themselves. The production follows three steps. First, we apply the PATCHY code to generate a dark matter field and an object distribution including nonlinear stochastic galaxy bias. Secondly, we run the halo/stellar distribution reconstruction HADRON code to assign masses to the various objects. This step uses the mass distribution as a function of local density and non-local indicators (i.e., tidal field tensor eigenvalues and relative halo exclusion separation for massive objects) from the reference simulation applied to the corresponding PATCHY dark matter and galaxy distribution. Finally, we apply the SUGAR code to build the light cones. The resulting MultiDark PATCHY mock light cones reproduce the number density, selection function, survey geometry, and in general within 1 $\sigma$, for arbitrary stellar mass bins, the power spectrum up to $k=0.3\,h\,{\rm Mpc}^{-1}$, the two-point correlation functions down to a few Mpc scales, and the three-point statistics of the BOSS DR11&DR12 galaxy samples.

We present distance scale measurements from the baryon acoustic oscillation signal in the CMASS and LOWZ samples from the Data Release 12 of the Baryon Oscillation Spectroscopic Survey (BOSS). The total volume probed is 14.5 Gpc$^3$, a 10 per cent increment from Data Release 11. From an analysis of the spherically averaged correlation function, we infer a distance to $z=0.57$ of $D_V(z)r^{\rm fid}_{\rm d}/r_ {\rm d}=2028\pm21$ Mpc and a distance to $z=0.32$ of $D_V(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1264\pm22$ Mpc assuming a cosmology in which $r^{\rm fid}_{\rm d}=147.10$ Mpc. From the anisotropic analysis, we find an angular diameter distance to $z=0.57$ of $D_{\rm A}(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1401\pm21$ Mpc and a distance to $z=0.32$ of $981\pm20$ Mpc, a 1.5 per cent and 2.0 per cent measurement respectively. The Hubble parameter at $z=0.57$ is $H(z)r_{\rm d}/r^{\rm fid}_{\rm d}=100.3\pm3.7$ km s$^{-1}$ Mpc$^{-1}$ and its value at $z=0.32$ is $79.2\pm5.6$ km s$^{-1}$ Mpc$^{-1}$, a 3.7 per cent and 7.1 per cent measurement respectively. These cosmic distance scale constraints are in excellent agreement with a $\Lambda$CDM model with cosmological parameters released by the recent Planck 2015 results.

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

We study the clustering of galaxies detected at $i<22.5$ in the Science Verification observations of the Dark Energy Survey (DES). Two-point correlation functions are measured using $2.3\times 10^6$ galaxies over a contiguous 116 deg$^2$ region in five bins of photometric redshift width $\Delta z = 0.2$ in the range $0.2 < z < 1.2.$ The impact of photometric redshift errors are assessed by comparing results using a template-based photo-$z$ algorithm (BPZ) to a machine-learning algorithm (TPZ). A companion paper (Leistedt et al 2015) presents maps of several observational variables (e.g. seeing, sky brightness) which could modulate the galaxy density. Here we characterize and mitigate systematic errors on the measured clustering which arise from these observational variables, in addition to others such as Galactic dust and stellar contamination. After correcting for systematic effects we measure galaxy bias over a broad range of linear scales relative to mass clustering predicted from the Planck $\Lambda$CDM model, finding agreement with CFHTLS measurements with $\chi^2$ of 4.0 (8.7) with 5 degrees of freedom for the TPZ (BPZ) redshifts. We test a "linear bias" model, in which the galaxy clustering is a fixed multiple of the predicted non-linear dark-matter clustering. The precision of the data allow us to determine that the linear bias model describes the observed galaxy clustering to $2.5\%$ accuracy down to scales at least $4$ to $10$ times smaller than those on which linear theory is expected to be sufficient.

We report a measurement of the large-scale 3-point correlation function of galaxies using the largest dataset for this purpose to date, 777, 202 Luminous Red Galaxies in the Sloan Digital Sky Survey Baryon Acoustic Oscillation Spectroscopic Survey (SDSS BOSS) DR12 CMASS sample. This work exploits the novel algorithm of Slepian & Eisenstein (2015b) to compute the multipole moments of the 3PCF in $\mathcal{O}(N^2)$ time, with $N$ the number of galaxies. Leading-order perturbation theory models the data well in a compressed basis where one triangle side is integrated out. We also present an accurate and computationally efficient means of estimating the covariance matrix. With these techniques the redshift-space linear and non-linear bias are measured, with 2.6% precision on the former if $\sigma_8$ is fixed. The data also indicates a $2.8\sigma$ preference for the BAO, confirming the presence of BAO in the 3-point function.

We present the results of the first test plates of the extended Baryon Oscillation Spectroscopic Survey. This paper focuses on the emission line galaxies (ELG) population targetted from the Dark Energy Survey (DES) photometry. We analyse the success rate, efficiency, redshift distribution, and clustering properties of the targets. From the 9000 spectroscopic redshifts targetted, 4600 have been selected from the DES photometry. The total success rate for redshifts between 0.6 and 1.2 is 71\% and 68\% respectively for a bright and faint, on average more distant, samples including redshifts measured from a single strong emission line. We find a mean redshift of 0.8 and 0.87, with 15 and 13\% of unknown redshifts respectively for the bright and faint samples. In the redshift range 0.6<z<1.2, for the most secure spectroscopic redshifts, the mean redshift for the bright and faint sample is 0.85 and 0.9 respectively. Star contamination is lower than 2\%. We measure a galaxy bias averaged on scales of 1 and 10~Mpc/h of 1.72 \pm 0.1 for the bright sample and of 1.78 \pm 0.12 for the faint sample. The error on the galaxy bias have been obtained propagating the errors in the correlation function to the fitted parameters. This redshift evolution for the galaxy bias is in agreement with theoretical expectations for a galaxy population with MB-5\log h < -21.0. We note that biasing is derived from the galaxy clustering relative to a model for the mass fluctuations. We investigate the quality of the DES photometric redshifts and find that the outlier fraction can be reduced using a comparison between template fitting and neural network, or using a random forest algorithm.

As part of the Sloan Digital Sky Survey IV the extended Baryon Oscillation Spectroscopic Survey (eBOSS) will improve measurements of the cosmological distance scale by applying the Baryon Acoustic Oscillation (BAO) method to quasar samples. eBOSS will adopt two approaches to target quasars over 7500 sq. deg. First, a "CORE" quasar sample will combine optical selection in ugriz using a likelihood-based routine called XDQSOz, with a mid-IR-optical color-cut. eBOSS CORE selection (to g < 22 OR r < 22) should return ~ 70 quasars per sq. deg. at redshifts 0.9 < z < 2.2 and ~7 z > 2.1 quasars per sq. deg. Second, a selection based on variability in multi-epoch imaging from the Palomar Transient Factory should recover an additional ~3-4 z > 2.1 quasars per sq. deg. to g < 22.5. A linear model of how imaging systematics affect target density recovers the angular distribution of eBOSS CORE quasars over 96.7% (76.7%) of the SDSS North (South) Galactic Cap area. The eBOSS CORE quasar sample should thus be sufficiently dense and homogeneous over 0.9 < z < 2.2 to yield the first few-percent-level BAO constraint near z~1.5. eBOSS quasars at z > 2.1 will be used to improve BAO measurements in the Lyman-alpha Forest. Beyond its key cosmological goals, eBOSS should be the next-generation quasar survey, comprising > 500,000 new quasars and > 500,000 uniformly selected spectroscopically confirmed 0.9 < z < 2.2 quasars. At the conclusion of eBOSS, the SDSS will have provided unique spectra of over 800,000 quasars.