A photometric redshift sample of Luminous Red Galaxies (hereafter LRGs) obtained from The DECam Legacy Survey (DECaLS) is analysed to probe cosmic distances by exploiting the wedge approach of the two-point correlation function. Although the cosmological information is highly contaminated by the uncertainties existing in the photometric redshifts from the galaxy map, an angular diameter distance can be probed at the perpendicular configuration in which the measured correlation function is minimally contaminated. An ensemble of wedged correlation functions selected up to a given threshold based on having the least contamination was studied in the previous work (arXiv:1903.09651v2 [astro-ph.CO]) using simulations, and the extracted cosmological information was unbiased within this threshold. We apply the same methodology for analysing the LRG sample from DECaLS which will provide the optical imaging for targeting two-thirds of the DESI footprint and measure the angular diameter distances at $z=0.69$ and $z=0.87$ to be $D_{A}(0.697)=(1499 \pm 77\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ and $D_{A}(0.874)=(1680 \pm 109\,\mathrm{Mpc})(r_{d}/r_{d,fid})$ with a fractional error of 5.14% and 6.48% respectively. We obtain a value of $H_{0}=67.59\pm5.52$ km/s/Mpc which supports the $H_0$ measured by all other BAO results and is consistent with $\Lambda$CDM model.

We develop a methodology to optimise the measurement of Baryon Acoustic Oscillation (BAO) from a given galaxy sample. In our previous work, we demonstrated that one can measure BAO from tracers in under-dense regions (voids). In this study, we combine the over-dense and under-dense tracers (galaxies & voids) to obtain better constraints on the BAO scale. To this end, we modify the de-wiggled BAO model with an additional parameter to describe both the BAO peak and the underlying exclusion pattern of void 2PCFs. We show that after applying BAO reconstruction to galaxies, the BAO peak scale of both galaxies and voids are unbiased using the modified model. Furthermore, we use a new 2PCF estimator for a multi-tracer analysis with galaxies and voids. In simulations, the joint sample improves by about 10% the constraint for the post-reconstruction BAO peak position compared to the result from galaxies alone, which is equivalent to an enlargement of the survey volume by 20%. Applying this method to the BOSS DR12 data, we have an 18% improvement for the low-z sample (0.2<z<0.5), but a worse constraint for the high-z sample (0.5<z<0.75), which is consistent with statistical fluctuations for the current survey volume. Future larger samples will give more robust improvements due to less statistical fluctuations.

We present a comparison between approximated methods for the construction of mock catalogs based on the halo-bias mapping technique. To this end, we use as reference a high resolution $N$-body simulation of 3840$^3$ dark matter particles on a 400$h^{-1}\rm{Mpc}$ cube box from the Multidark suite. In particular, we explore parametric versus non-parametric bias mapping approaches and compare them at reproducing the halo distribution in terms of the two and three point statistics down to $\sim 10^8\,{\rm M}_{\odot}\,h^{-1}$ halo masses. Our findings demonstrate that the parametric approach remains inaccurate even including complex deterministic and stochastic components. On the contrary, the non-parametric one is indistinguishable from the reference $N$-body calculation in the power-spectrum beyond $k=1\,h\,{\rm Mpc}^{-1}$, and in the bispectrum for typical configurations relevant to baryon acoustic oscillation analysis. We conclude, that approaches which extract the full bias information from $N$-body simulations in a non-parametric fashion are ready for the analysis of the new generation of large scale structure surveys.

We perform a simulation with Galacticus, a semi-analytical galaxy formation model, to predict the number counts of H$\alpha$ and [OIII] emitting galaxies. With a state-of-the-art N-body simulation, UNIT, we first calibrate Galacticus with the current observation of H$\alpha$ luminosity function. The resulting model coupled with a dust attenuation model, can reproduce the current observations, including the H$\alpha$ luminosity function from HiZELS and number density from WISP. We extrapolate the model prediction to higher redshift and the result is found to be consistent with previous investigations. We then use the same galaxy formation model to predict the number counts for [OIII] emitting galaxies. The result provides further validation of our galaxy formation model and dust model. We present number counts of H$\alpha$ and [OIII] emission line galaxies for three different line flux limits: $5\times10^{-17}$erg/s/cm$^{2}$, $1\times10^{-16}$erg/s/cm$^{2}$ (6.5$\sigma$ nominal depth for WFIRST GRS), and $2\times10^{-16}$erg/s/cm$^{2}$ (3.5$\sigma$ depth of Euclid GRS). At redshift $2<z<3$, our model predicts that WFIRST can observe hundreds of [OIII] emission line galaxies per square degree with a line flux limit of $1\times10^{-16}$erg/s/cm$^{2}$. This will provide accurate measurement of large scale structure to probe dark energy over a huge cosmic volume to an unprecedented high redshift. Finally, we compare the flux ratio of H$\alpha$/[OIII] within the redshift range of $0<z<3$. Our results show the known trend of increasing H$\alpha$/[O III] flux ratio with H$\alpha$ flux at low redshift, which becomes a weaker trend at higher redshifts.

We present the measurements and modelling of the small-to-intermediate scale (0.1--25 Mpc/h) projected and three-dimensional (3D) redshift-space two-point correlation functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7. We find a clear dependence of galaxy clustering on luminosity in both projected and redshift spaces, generally being stronger for more luminous samples. The measurements are successfully interpreted within the halo occupation distribution (HOD) framework with central and satellite velocity bias parameters to describe galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. In agreement with previous studies, we find that more luminous galaxies reside in more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD constraints. Moreover, we find that luminous central galaxies are not at rest at the halo centres, with the velocity dispersion about 30% that of the dark matter. Such a relative motion may reflect the consequence of galaxy and halo mergers, and we find that central galaxies in lower mass haloes tend to be more relaxed with respect to their host haloes. The motion of satellite galaxies in luminous samples is consistent with their following that of the dark matter. For faint samples, satellites tends to have slower motion, with velocity dispersion inside haloes about 85% that of the dark matter. We discuss possible applications of the velocity bias constraints on studying galaxy evolution and cosmology. In the appendix, we characterize the distribution of galaxy redshift measurement errors, which is well described by a Gaussian-convolved double exponential distribution.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

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 present the UNIT $N$-body cosmological simulations project, designed to provide precise predictions for nonlinear statistics of the galaxy distribution. We focus on characterizing statistics relevant to emission line and luminous red galaxies in the current and upcoming generation of galaxy surveys. We use a suite of precise particle mesh simulations (FastPM) as well as with full $N$-body calculations with a mass resolution of $\sim 1.2\times10^9\,h^{-1}$M$_{\odot}$ to investigate the recently suggested technique of Angulo & Pontzen 2016 to suppress the variance of cosmological simulations We study redshift space distortions, cosmic voids, higher order statistics from $z=2$ down to $z=0$. We find that both two- and three-point statistics are unbiased. Over the scales of interest for baryon acoustic oscillations and redshift-space distortions, we find that the variance is greatly reduced in the two-point statistics and in the cross correlation between halos and cosmic voids, but is not reduced significantly for the three-point statistics. We demonstrate that the accuracy of the two-point correlation function for a galaxy survey with effective volume of 20 ($h^{-1}$Gpc)$^3$ is improved by about a factor of 40, indicating that two pairs of simulations with a volume of 1 ($h^{-1}$Gpc)$^3$ lead to the equivalent variance of $\sim$150 such simulations. The $N$-body simulations presented here thus provide an effective survey volume of about seven times the effective survey volume of DESI or Euclid. The data from this project, including dark matter fields, halo catalogues, and their clustering statistics, are publicly available at http://www.unitsims.org.

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.

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.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

We perform a tomographic structure growth and expansion rate analysis using the monopole, quadrupole and hexadecapole of the redshift-space galaxy power spectrum derived from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 combined sample, which covers the redshift range of $0.20<z<0.75$. By allowing for overlap between neighbouring redshift slices in order to extract information on the light-cone, we successfully obtain joint BAO and RSD constraints with a precision of $2-3\%$ for $D_A$, $3-10\%$ for $H$ and $9-12\%$ for $f\sigma_8$ with a redshift resolution of $\Delta z\sim0.04$. Our measurement is consistent with that presented in arXiv:1709.05173, where the analysis is performed in configuration space. We apply our measurement to constrain the $f(R)$ gravity model, and find that the 95\% CL upper limit of ${\rm log_{10}}B_0$ can be reduced by 11\% by our tomographic BAO and RSD measurement.

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.

We perform a tomographic analysis of structure growth and expansion rate from the anisotropic galaxy clustering of the combined sample of Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12, which covers the redshift range of $0.2<z<0.75$. In order to extract the redshift information of anisotropic galaxy clustering, we analyse this data set in nine overlapping redshift slices in configuration space and perform the joint constraints on the parameters $(D_V, F_{\mathrm{AP}}, f\sigma_8)$ using the correlation function multipoles. The analysis pipeline is validated using the MultiDark-Patchy mock catalogues. We obtain a measurement precision of $1.5\%-2.9\%$ for $D_V$, $5.2\%-9\%$ for $F_{\mathrm{AP}}$ and $13.3\%-24\%$ for $f \sigma_8$, depending on the effective redshift of the slices. We report a joint measurement of $(D_V, F_{\mathrm{AP}}, f\sigma_8)$ with the full covariance matrix in nine redshift slices. We use our joint BAO and RSD measurement combined with external datasets to constrain the gravitational growth index $\gamma$, and find $\gamma=0.656 \pm 0.057$, which is consistent with the $\Lambda$CDM prediction within 95\% CL.

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.

We further develop the description of redshift space distortions within the Effective Field Theory of Large Scale Structures. First, we generalize the counterterms to include the effect of baryonic physics and primordial non-Gaussianity. Second, we evaluate the IR-resummation of the dark matter power spectrum in redshift space. This requires us to identify a controlled approximation that makes the numerical evaluation straightforward and efficient. Third, we compare the predictions of the theory at one loop with the power spectrum from numerical simulations up to $\ell=6$. We find that the IR-resummation allows us to correctly reproduce the BAO peak. The $k$-reach, or equivalently the precision for a given $k$, depends on additional counterterms that need to be matched to simulations. Since the non-linear scale for the velocity is expected to be longer than the one for the overdensity, we consider a minimal and a non-minimal set of counterterms. The quality of our numerical data makes it hard to firmly establish the performance of the theory at high wavenumbers. Within this limitation, we find that the theory at redshift $z=0.56$ and up to $\ell=2$ matches the data to percent level approximately up to $k \sim 0.13 \, h { \rm Mpc^{-1}}$ or $k \sim 0.18 \, h { \rm Mpc^{-1}}$, depending on the number of counterterms used, with potentially large improvement over former analytical techniques.

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 study on the robustness of the covariance matrix estimation for galaxy clustering measurements depending on the cosmological parameters and galaxy bias. To this end, we have produced 9000 galaxy mock catalogues relying on the effective Zel'dovich approximation implemented in the EZmocks computer code, using different input cosmological models and bias parameters. The reference catalogue has also been produced with this code making our study insensitive to the approximation at least on a relative-qualitative level. Our findings indicate that the covariance matrix is insensitive to the input power spectrum (including $\sigma_8$), as long as the 2- and 3-point galaxy clustering measurements agree with the given data. In fact, the covariance matrix shows a bias at small scales ($r\lesssim40 h^{-1}$Mpc) when the chosen galaxy bias parameters yield a 3-point statistics, which is not compatible with the reference one within the error bars, even though the 2-point statistics agree within 1%. Nevertheless, the error becomes negligible at large scales making the covariance matrix still reliable for data analysis using only measurements in that regime (e.g., measuring baryon acoustic oscillations). High precision in cosmological parameter estimation is expected for covariance matrices extracted from mock galaxy catalogues which take accurately into account both the 2- and the 3- point statistics. This is independent of whether this is achieved by using the right cosmology and galaxy bias (which are not a priori known) or just any combination of both fitting the net observed galaxy clustering.

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

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.

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.

Reliable extraction of cosmological information from clustering measurements of galaxy surveys requires estimation of the error covariance matrices of observables. The accuracy of covariance matrices is limited by our ability to generate sufficiently large number of independent mock catalogs that can describe the physics of galaxy clustering across a wide range of scales. Furthermore, galaxy mock catalogs are required to study systematics in galaxy surveys and to test analysis tools. In this investigation, we present a fast and accurate approach for generation of mock catalogs for the upcoming galaxy surveys. Our method relies on low-resolution approximate gravity solvers to simulate the large scale dark matter field, which we then populate with halos according to a flexible nonlinear and stochastic bias model. In particular, we extend the \textsc{patchy} code with an efficient particle mesh algorithm to simulate the dark matter field (the \textsc{FastPM} code), and with a robust MCMC method relying on the \textsc{emcee} code for constraining the parameters of the bias model. Using the halos in the BigMultiDark high-resolution $N$-body simulation as a reference catalog, we demonstrate that our technique can model the bivariate probability distribution function (counts-in-cells), power spectrum, and bispectrum of halos in the reference catalog. Specifically, we show that the new ingredients permit us to reach percentage accuracy in the power spectrum up to $k\sim 0.4\; \,h\,{\rm Mpc}^{-1}$ (within 5\% up to $k\sim 0.6\; \,h\,{\rm Mpc}^{-1}$) with accurate bispectra improving previous results based on Lagrangian perturbation theory.

We present measurements of the anisotropic galaxy clustering from the Data Release 9 (DR9) CMASS sample of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We analyze the broad-range shape of the monopole and quadrupole correlation functions to obtain constraints, at the effective redshift $z=0.57$ of the sample, on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)\sigma_8(z)$, the physical matter density $\Omega_m h^2$, and the biased amplitude of matter fluctuation b\sigma_8(z). We obtain {$H(0.57)$, $D_A(0.57)$, $f(0.57)\sigma_8(0.57)$, $\Omega_m h^2$, $b\sigma_8(0.57)$} = {$87.6_{-6.8}^{+6.7}$, $1396\pm73$, $0.126_{-0.010}^{+0.008}$, $1.19\pm0.14$, $0.428\pm0.066$} and their covariance matrix as well. The parameters which are not well constrained by our of galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets (i.e., CMB) are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e. CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. We show that the major power on constraining dark energy from the anisotropic galaxy clustering signal, as compared to the angular-averaged one (monopole), arises from including the normalized growth rate $f(z)\sigma_8(z)$. In the case of the wCDM cosmological model our single-probe CMASS constraints, combined with CMB (WMAP9+SPT), yield a value for the dark energy equation of state parameter of $w=-0.90\pm0.11$. Therefore, it is important to include $f(z)\sigma_8(z)$ while investigating the nature of dark energy with current and upcoming large-scale galaxy surveys.

We present a new methodology to generate mock halo or galaxy catalogues, which have accurate clustering properties, nearly indistinguishable from full $N$-body solutions, in terms of the one-point, two-point, and three-point statistics. In particular, the agreement is remarkable, within $1\%$ up to $k=0.55$ $h$Mpc$^{-1}$ and down to $r=10$ $h^{-1}$Mpc, for the power spectrum and two-point correlation function respectively, while the bispectrum agrees in general within $20\%$ for different scales and shapes. Our approach is based on the Zel'dovich approximation, however, effectively including with the simple prescriptions the missing physical ingredients, and stochastic scale-dependent, non-local and nonlinear biasing contributions. The computing time and memory required to produce one mock is similar to that using the log-normal model. With high accuracy and efficiency, the effective Zel'dovich approximation mocks (EZmocks) provide a reliable and practical method to produce massive mock galaxy catalogues for the analysis of large-scale structure measurements.

Cosmic voids found in galaxy surveys are defined based on the galaxy distribution in redshift space. We show that the large scale distribution of voids in redshift space traces the fluctuations in the dark matter density field \delta(k) (in Fourier space with \mu being the line of sight projected k-vector): \delta_v^s(k) = (1 + \beta_v \mu^2) b^s_v \delta(k), with a beta factor that will be in general different than the one describing the distribution of galaxies. Only in case voids could be assumed to be quasi-local transformations of the linear (Gaussian) galaxy redshift space field, one gets equal beta factors \beta_v=\beta_g=f/b_g with f being the growth rate, and b_g, b^s_v being the galaxy and void bias on large scales defined in redshift space. Indeed, in our mock void catalogs we measure void beta factors being in good agreement with the galaxy one. Further work needs to be done to confirm the level of accuracy of the beta factor equality between voids and galaxies, but in general the void beta factor needs to be considered as a free parameter for RSD studies.

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.

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

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.

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

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.

Current and future large redshift surveys, as the Sloan Digital Sky Survey IV extended Baryon Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) or the Dark Energy Spectroscopic Instrument (DESI), will use emission-line galaxies (ELG) to probe cosmological models by mapping the large-scale structure of the Universe in the redshift range $0.6 < z < 1.7$. With current data, we explore the halo-galaxy connection by measuring three clustering properties of $g$-selected ELGs as matter tracers in the redshift range $0.6 < z < 1$: (i) the redshift-space two-point correlation function using spectroscopic redshifts from the BOSS ELG sample and VIPERS; (ii) the angular two-point correlation function on the footprint of the CFHT-LS; (iii) the galaxy-galaxy lensing signal around the ELGs using the CFHTLenS. We interpret these observations by mapping them onto the latest high-resolution MultiDark Planck N-body simulation, using a novel (Sub)Halo-Abundance Matching technique that accounts for the ELG incompleteness. ELGs at $z\sim0.8$ live in halos of $(1\pm 0.5)\times10^{12}\,h^{-1}$M$_{\odot}$ and 22.5$\pm2.5$% of them are satellites belonging to a larger halo. The halo occupation distribution of ELGs indicates that we are sampling the galaxies in which stars form in the most efficient way, according to their stellar-to-halo mass ratio.

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

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

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

With the largest spectroscopic galaxy survey volume drawn from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), we can extract cosmological constraints from the measurements of redshift and geometric distortions at quasi-linear scales (e.g. above 50 $h^{-1}$Mpc). We analyze the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS galaxy sample, at the effective redshift $z=0.59$, to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)\sigma_8(z)$, and the physical matter density $\Omega_mh^2$. We obtain robust measurements by including a polynomial as the model for the systematic errors, and find it works very well against the systematic effects, e.g., ones induced by stars and seeing. We provide accurate measurements $\{D_A(0.59)r_{s,fid}/r_s$ $\rm Mpc$, $H(0.59)r_s/r_{s,fid}$ $km s^{-1} Mpc^{-1}$, $f(0.59)\sigma_8(0.59)$, $\Omega_m h^2\}$ = $\{1427\pm26$, $97.3\pm3.3$, $0.488 \pm 0.060$, $0.135\pm0.016\}$, where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$ Mpc is the sound scale of the fiducial cosmology used in this study. The parameters which are not well constrained by our galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets, e.g., cosmic microwave background (CMB), are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e., CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. The uncertainty on the dark energy equation of state parameter, $w$, from CMB+CMASS is about 8 per cent. The uncertainty on the curvature fraction, $\Omega_k$, is 0.3 per cent. We do not find deviation from flat $\Lambda$CDM.

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

We have developed an analytical method based on forward-modeling techniques to characterize the high-mass end of the red sequence (RS) galaxy population at redshift $z\sim0.55$, from the DR10 BOSS CMASS spectroscopic sample, which comprises $\sim600,000$ galaxies. The method, which follows an unbinned maximum likelihood approach, allows the deconvolution of the intrinsic CMASS colour-colour-magnitude distributions from photometric errors and selection effects. This procedure requires modeling the covariance matrix for the i-band magnitude, g-r colour and r-i colour using Stripe 82 multi-epoch data. Our results indicate that the error-deconvolved intrinsic RS distribution is consistent, within the photometric uncertainties, with a single point ($<0.05~{\rm{mag}}$) in the colour-colour plane at fixed magnitude, for a narrow redshift slice. We have computed the high-mass end ($^{0.55}M_i \lesssim -22$) of the $^{0.55}i$-band RS Luminosity Function (RS LF) in several redshift slices within the redshift range $0.52<z<0.63$. In this narrow redshift range, the evolution of the RS LF is consistent, within the uncertainties in the modeling, with a passively-evolving model with $\Phi_* = (7.248 \pm 0.204) \times10^{-4}$ Mpc$^{-3}$ mag$^{-1}$, fading at a rate of $1.5\pm0.4$ mag per unit redshift. We report RS completeness as a function of magnitude and redshift in the CMASS sample, which will facilitate a variety of galaxy-evolution and clustering studies using BOSS. Our forward-modeling method lays the foundations for future studies using other dark-energy surveys like eBOSS or DESI, which are affected by the same type of photometric blurring/selection effects.

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $\Lambda$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$\sigma$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

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 model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $\chi^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies provides a standard ruler to measure the accelerated expansion of the Universe. To extract all available information about dark energy, it is necessary to measure a standard ruler in the local, z<0.2, universe where dark energy dominates most the energy density of the Universe. Though the volume available in the local universe is limited, it is just big enough to measure accurately the long 100 Mpc/h wave-mode of the BAO. Using cosmological N-body simulations and approximate methods based on Lagrangian perturbation theory, we construct a suite of a thousand light-cones to evaluate the precision at which one can measure the BAO standard ruler in the local universe. We find that using the most massive galaxies on the full sky (34,000 sq. deg.), i.e. a K(2MASS)<14 magnitude-limited sample, one can measure the BAO scale up to a precision of 4\% and 1.2\% using reconstruction). We also find that such a survey would help to detect the dynamics of dark energy.Therefore, we propose a 3-year long observational project, named the Low Redshift survey at Calar Alto (LoRCA), to observe spectroscopically about 200,000 galaxies in the northern sky to contribute to the construction of aforementioned galaxy sample. The suite of light-cones is made available to the public.

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.

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 investigate the necessary methodology to optimally measure the baryon acoustic oscillation (BAO) signal, from voids based on galaxy redshift catalogues. To this end, we study the dependency of the BAO signal on the population of voids classified by their sizes. We find for the first time the characteristic features of the correlation function of voids including the first robust detection of BAOs in mock galaxy catalogues. These show an anti-correlation around the scale corresponding to the smallest size of voids in the sample (the void exclusion effect), and dips at both sides of the BAO peak, which can be used to determine the significance of the BAO signal without any priori model. Furthermore, our analysis demonstrates that there is a scale dependent bias for different populations of voids depending on the radius, with the peculiar property that the void population with the largest BAO significance corresponds to tracers with approximately zero bias on the largest scales. We further investigate the methodology on an additional set of 1,000 realistic mock galaxy catalogues reproducing the SDSS-III/BOSS CMASS DR11 data, to control the impact of sky mask and radial selection function. Our solution is based on generating voids from randoms including the same survey geometry and completeness, and a post-processing cleaning procedure in the holes and at the boundaries of the survey. The methodology and optimal selection of void populations validated in this work have been used to perform the first BAO detection from voids in observations, presented in a companion paper.

We present a novel parameter-free cosmological void finder (\textsc{dive}, Delaunay TrIangulation Void findEr) based on Delaunay Triangulation (DT), which efficiently computes the empty spheres constrained by a discrete set of tracers. We define the spheres as DT voids, and describe their properties, including an universal density profile together with an intrinsic scatter. We apply this technique on 100 halo catalogues with volumes of 2.5\,$h^{-1}$Gpc side each, with a bias and number density similar to the BOSS CMASS Luminous Red Galaxies, performed with the \textsc{patchy} code. Our results show that there are two main species of DT voids, which can be characterised by the radius: they have different responses to halo redshift space distortions, to number density of tracers, and reside in different dark matter environments. Based on dynamical arguments using the tidal field tensor, we demonstrate that large DT voids are hosted in expanding regions, whereas the haloes used to construct them reside in collapsing ones. Our approach is therefore able to efficiently determine the troughs of the density field from galaxy surveys, and can be used to study their clustering. We further study the power spectra of DT voids, and find that the bias of the two populations are different, demonstrating that the small DT voids are essentially tracers of groups of haloes.

The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

We present a Bayesian reconstruction method which maps a galaxy distribution from redshift-space to real-space inferring the distances of the individual galaxies. The method is based on sampling density fields assuming a lognormal prior with a likelihood given by the negative binomial distribution function modelling stochastic bias. We assume a deterministic bias given by a power law relating the dark matter density field to the expected halo or galaxy field. Coherent redshift-space distortions are corrected in a Gibbs-sampling procedure by moving the galaxies from redshift-space to real-space according to the peculiar motions derived from the recovered density field using linear theory with the option to include tidal field corrections from second order Lagrangian perturbation theory. The virialised distortions are corrected by sampling candidate real-space positions (being in the neighbourhood of the observations along the line of sight), which are compatible with the bulk flow corrected redshift-space position adding a random dispersion term in high density collapsed regions. The latter are defined according to the eigenvalues of the Hessian. This approach presents an alternative method to estimate the distances to galaxies using the three dimensional spatial information, and assuming isotropy. Hence the number of applications is very broad. In this work we show the potential of this method to constrain the growth rate up to $k$ ~ 0.3 $h$ Mpc$^{-1}$. Furthermore it could be useful to correct for photo-metric redshift errors, and to obtain improved BAO reconstructions.

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

We explore the bluer star-forming population of the Sloan Digital Sky Survey (SDSS) III/BOSS CMASS DR11 galaxies at $z>0.55$ to quantify their differences, in terms of redshift-space distortions and large-scale bias, with respect to the luminous red galaxy sample. We perform a qualitative analysis to understand the significance of these differences and whether we can model and reproduce them in mock catalogs. Specifically, we measure galaxy clustering in CMASS on small and intermediate scales ($r\lesssim 50\,h^{-1}$Mpc) by computing the two-point correlation function $-$ both projected and redshift-space $-$ of these galaxies, and a new statistic, $\Sigma(\pi)$, able to provide robust information about redshift-space distortions and large-scale bias. We interpret our clustering measurements by adopting a Halo Occupation Distribution (HOD) scheme that maps them onto high-resolution N-body cosmological simulations to produce suitable mock galaxy catalogs. The traditional HOD prescription can be applied to the red and the blue samples, independently, but this approach is unphysical since it allows the same mock galaxies to be either red or blue. To overcome this failure, we modify the standard formulation and infer the red and the blue mock catalogs directly from the full one, so that they are complementary and non-overlapping. This separation is performed by matching the observed CMASS red and blue galaxy fractions and produces reliable and accurate models.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We study the relation between halo mass and its environment from a probabilistic perspective. We find that halo mass depends not only on local dark matter density, but also on non-local quantities such as the cosmic web environment and the halo-exclusion effect. Given these accurate relations, we have developed the HADRON-code (Halo mAss Distribution ReconstructiON), a technique which permits us to assign halo masses to a distribution of haloes in three-dimensional space. This can be applied to the fast production of mock galaxy catalogues, by assigning halo masses, and reproducing accurately the bias for different mass cuts. The resulting clustering of the halo populations agree well with that drawn from the BigMultiDark $N$-body simulation: the power spectra are within 1-$\sigma$ up to scales of $k=0.2\,h\,{\rm Mpc}^{-1}$, when using augmented Lagrangian perturbation theory based mock catalogues. Only the most massive haloes show a larger deviation. For these, we find evidence of the halo-exclusion effect. A clear improvement is achieved when assigning the highest masses to haloes with a minimum distance separation. We also compute the 2- and 3-point correlation functions, and find an excellent agreement with $N$-body results. Our work represents a quantitative application of the cosmic web classification. It can have further interesting applications in the multi-tracer analysis of the large-scale structure for future galaxy surveys.

Observations of redshift-space distortions in spectroscopic galaxy surveys offer an attractive method for observing the build-up of cosmological structure, which depends both on the expansion rate of the Universe and our theory of gravity. In preparation for analysis of redshift-space distortions from the Baryon Oscillation Spectroscopic Survey (BOSS) final data release we compare a number of analytic and phenomenological `streaming' models, specified in configuration space, to mock catalogs derived in different ways from several N-body simulations. The galaxies in each mock catalog have properties similar to those of the higher redshift galaxies measured by BOSS but differ in the details of how small-scale velocities and halo occupancy are determined. We find that all of the analytic models fit the simulations over a limited range of scales while failing at small scales. We discuss which models are most robust and on which scales they return reliable estimates of the rate of growth of structure: we find that models based on some form of resummation can fit our N-body data for BOSS-like galaxies above $30\,h^{-1}$Mpc well enough to return unbiased parameter estimates.

Future dark energy experiments will require better and more accurate theoretical predictions for the baryonic acoustic oscillations (BAO) signature in the spectrum of cosmological perturbations. Here, we use large N-body simulations of the \LambdaCDM Planck cosmology to study any possible systematic shifts and damping in BAO due to the impact of nonlinear gravitational growth of structure, scale dependent and non-local bias, and redshift-space distortions. The effect of cosmic variance is largely reduced by dividing the tracer power spectrum by that from a BAO-free simulation starting with the same phases. This permits us to study with unprecedented accuracy (better than 0.02% for dark matter and 0.07% for low-bias halos) small shifts of the pristine BAO wavenumbers towards larger k, and non-linear damping of BAO wiggles in the power spectrum of dark matter and halo populations in the redshift range z=0-1. For dark matter, we provide an accurate parametrization of the evolution of \alpha as a function of the linear growth factor D(z). For halo samples, with bias ranging from 1.2 to 2.8, we measure a typical BAO shift of ~0.25%, observed in real-space, which does not show an appreciable evolution with redshift within the uncertainties. Moreover, we report a constant shift as a function of halo bias. We find a different evolution of the damping of the acoustic feature in all halo samples as compared to dark matter with haloes suffering less damping, and also find some weak dependence on bias. A larger BAO shift and damping is measured in redshift-space which can be well explained by linear theory due to redshift-space distortions. A clear modulation in phase with the acoustic scale is observed in the scale-dependent halo bias due to the presence of the baryonic acoustic oscillations.

We investigate the cosmological implications of the latest growth of structure measurement from the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 with particular focus on the sum of the neutrino masses, $\sum m_{\nu}$. We examine the robustness of the cosmological constraints from the Baryon Acoustic Oscillation (BAO) scale, the Alcock-Paczynski effect and redshift-space distortions ($D_V/r_s$, $F_{\rm AP}$, $f\sigma_8$) of \citet{Beutler:2013b}, when introducing a neutrino mass in the power spectrum template. We then discuss how the neutrino mass relaxes discrepancies between the Cosmic Microwave Background (CMB) and other low-redshift measurements within $\Lambda$CDM. Combining our cosmological constraints with WMAP9 yields $\sum m_{\nu} = 0.36\pm0.14\,$eV ($68\%$ c.l.), which represents a $2.6\sigma$ preference for non-zero neutrino mass. The significance can be increased to $3.3\sigma$ when including weak lensing results and other BAO constraints, yielding $\sum m_{\nu} = 0.35\pm0.10\,$eV ($68\%$ c.l.). However, combining CMASS with Planck data reduces the preference for neutrino mass to $\sim 2\sigma$. When removing the CMB lensing effect in the Planck temperature power spectrum (by marginalising over $A_{\rm L}$), we see shifts of $\sim 1\sigma$ in $\sigma_8$ and $\Omega_m$, which have a significant effect on the neutrino mass constraints. In case of CMASS plus Planck without the $A_{\rm L}$-lensing signal, we find a preference for a neutrino mass of $\sum m_{\nu} = 0.34\pm 0.14\,$eV ($68\%$ c.l.), in excellent agreement with the WMAP9+CMASS value. The constraint can be tightened to $3.4\sigma$ yielding $\sum m_{\nu} = 0.36\pm 0.10\,$eV ($68\%$ c.l.) when weak lensing data and other BAO constraints are included.

We present a method to produce mock galaxy catalogues with efficient perturbation theory schemes, which match the number density, power spectra and bispectra in real and in redshift space from N-body simulations. The essential contribution of this work is the way in which we constrain the bias parameters in the PATCHY-code. In addition of aiming at reproducing the two-point statistics, we seek the set of bias parameters, which constrain the univariate halo probability distribution function (PDF) encoding higher-order correlation functions. We demonstrate that halo catalogues based on the same underlying dark matter field with a fix halo number density, and accurately matching the power spectrum (within 2%), can lead to very different bispectra depending on the adopted halo bias model. A model ignoring the shape of the halo PDF can lead to deviations up to factors of 2. The catalogues obtained additionally constraining the shape of the halo PDF can significantly lower the discrepancy in the three-point statistics, yielding closely unbiased bispectra both in real and in redshift space; which are in general compatible with those corresponding to an N-body simulation within 10% (deviating at most up to 20%). Our calculations show that the constant linear bias of ~2 for Luminous Red Galaxy (LRG) like galaxies seen in the power spectrum, mainly comes from sampling halos in high density peaks, choosing a high density threshold rather than from a factor multiplying the dark matter density field. Our method contributes towards an efficient modelling of the halo/galaxy distribution required to estimate uncertainties in the clustering measurements from galaxy redshift surveys. We have also demonstrated that it represents a powerful tool to test various bias models.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present improved methodology for including covariance matrices in the error budget of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering measurements, revisiting Data Release 9 (DR9) analyses, and describing a method that is used in DR10/11 analyses presented in companion papers. The precise analysis method adopted is becoming increasingly important, due to the precision that BOSS can now reach: even using as many as 600 mock catalogues to estimate covariance of 2-point clustering measurements can still lead to an increase in the errors of ~20%, depending on how the cosmological parameters of interest are measured. In this paper we extend previous work on this contribution to the error budget, deriving formulae for errors measured by integrating over the likelihood, and to the distribution of recovered best-fit parameters fitting the simulations also used to estimate the covariance matrix. Both are situations that previous analyses of BOSS have considered. We apply the formulae derived to Baryon Acoustic Oscillation (BAO) and Redshift-Space Distortion (RSD) measurements from BOSS in our companion papers. To further aid these analyses, we consider the optimum number of bins to use for 2-point measurements using the monopole power spectrum or correlation function for BAO, and the monopole and quadrupole moments of the correlation function for anisotropic-BAO and RSD measurements.

We present a method to measure the Hubble parameter $H(z)$ and the angular diameter distance $D_A(z)$ simultaneously from the two-dimensional matter power spectrum from galaxy surveys with broad sky coverage. We validate this method by applying it to the LasDamas mock galaxy catalogs. Then we apply this method to Sloan Digital Sky Survey (SDSS) Data Release 7 and obtain measurements of $\Omega_mh^2=0.1268 \pm 0.0085$, $H(z=0.35)=81.3\pm 3.8$km/s/Mpc, $D_A(z=0.35) = 1037\pm44$Mpc, without assuming a dark energy model or a flat universe. We also find that the derived parameters $H(0.35)r_s(z_d)/c=0.0431 \pm 0.0018$ and $D_A(0.35)/r_s(z_d)=6.48 \pm 0.25$. These are in excellent agreement with similar measurements from the two-dimensional correlation function of the same data.

We present a simple and efficient phenomenological model for the two-dimensional two-point galaxy correlation function that works well over a wide range of scales, from large scales down to scales as small as 25Mpc/h. Our model incorporates nonlinear effects, a scale-dependent galaxy bias on small scales, and allows the redshift-space distortions to be scale and direction dependent. We validate our model using LasDamas mock catalogs, and apply it to the Sloan Digital Sky Survey (SDSS) DR7 Luminous Red Galaxies (LRGs). Using only the monopole and quadrupole of the correlation function measured from the SDSS DR7 LRGs, we obtain improved measurements H(z)r_s(z_d)/c=0.0433\pm 0.0042, D_A(z)/r_s(z_d)=6.59\pm 0.46, and f(z)sigma_8(z)=0.429\pm 0.089 at z=0.35, using the scale range of 25<s<120Mpc/h. We expect our results and model to be useful in tightening dark energy and gravity constraints from the full analysis of current and future galaxy clustering data.

Galaxy clustering data can be used to measure the cosmic expansion history H(z), the angular-diameter distance D_A(z), and the linear redshift-space distortion parameter beta(z). Here we present a method for using effective multipoles of the galaxy two-point correlation function (\xi_0(s), \xi_2(s), \xi}_4(s), and \xi_6(s), with s denoting the comoving separation) to measure H(z), D_A(z)$, and beta(z), and validate it using LasDamas mock galaxy catalogs. Our definition of effective multipoles explicitly incorporates the discreteness of measurements, and treats the measured correlation function and its theoretical model on the same footing. We find that for the mock data, \xi_0+\xi_2+\xi_4 captures nearly all the information, and gives significantly stronger constraints on H(z), D_A(z), and beta(z), compared to using only \xi_0+\xi_2. We apply our method to the sample of luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) without assuming a dark energy model or a flat Universe. We find that \xi}_4(s) deviates on scales of s<60Mpc/h from the measurement from mock data (in contrast to \xi_0(s), \xi_2(s), and \xi_6(s)), thus we only use \xi_0+\xi_2 for our fiducial constraints. We obtain {H(0.35), D_A(0.35), Omega_mh^2, beta(z)} = {79.6_{-8.7}^{+8.3} km/s/Mpc, 1057_{-87}^{+88}Mpc, 0.103\pm0.015, 0.44\pm0.15} using \xi_0+\xi_2. We find that H(0.35)r_s(z_d)/c and D_A(0.35)/r_s(z_d) (where r_s(z_d) is the sound horizon at the drag epoch) are more tightly constrained: {H(0.35)r_s(z_d)/c, D_A(0.35)/r_s(z_d)} = {0.0437_{-0.0043}^{+0.0041}, 6.48_{-0.43}^{+0.44}\} using \xi_0+\xi_2.

We present measurements of the angular diameter distance to and Hubble parameter at z=0.57 from the measurement of the baryon acoustic peak in the correlation of galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey. Our analysis is based on a sample from Data Release 9 of 264,283 galaxies over 3275 square degrees in the redshift range 0.43<z<0.70. We use two different methods to provide robust measurement of the acoustic peak position across and along the line of sight in order to measure the cosmological distance scale. We find D_A(0.57) = 1408 +/- 45 Mpc and H(0.57) = 92.9 +/- 7.8 km/s/Mpc for our fiducial value of the sound horizon. These results from the anisotropic fitting are fully consistent with the analysis of the spherically averaged acoustic peak position presented in Anderson et al, 2012. Our distance measurements are a close match to the predictions of the standard cosmological model featuring a cosmological constant and zero spatial curvature.

We analyze the 2D correlation function of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of massive galaxies of the ninth data release to measure cosmic expansion H and the angular diameter distance D_A at a mean redshift of <z> = 0.57. We apply, for the first time, a new correlation function technique called clustering wedges. Using a physically motivated model, the anisotropic baryonic acoustic feature in the galaxy sample is detected at a significance level of 4.7 sigma compared to a featureless model. The baryonic acoustic feature is used to obtain model independent constraints cz/H/r_s = 12.28 +- 0.82 (6.7 per-cent accuracy) and D_A/r_s = 9.05 +- 0.27 (3.0 per-cent) with a correlation coefficient of -0.5, where r_s is the sound horizon scale at the end of the baryonic drag era. We conduct thorough tests on the data and 600 simulated realizations, finding robustness of the results regardless of the details of the analysis method. Combining with r_s constraints from the Cosmic Microwave Background we obtain H(0.57) = 90.8 +- 6.2 kms-1Mpc-1 and D_A(0.57) = 1386 +- 45 Mpc. We use simulations to forecast results of the final BOSS CMASS data set. We apply the reconstruction technique on the simulations demonstrating that the sharpening of the anisotropic baryonic acoustic feature should improve the detection as well as tighten constraints of H and D_A by 30 per-cent on average.

We explore the cosmological implications of the clustering wedges, xi_perp(s) and xi_para(s), of the CMASS Data Release 9 (DR9) sample of the Baryon Oscillation Spectroscopic Survey (BOSS). These clustering wedges are defined by averaging the full two-dimensional correlation function, xi(mu,s), over the ranges 0<mu<0.5 and 0.5<mu<1, respectively. These measurements allow us to constrain the parameter combinations D_A(z)/r_s(z_d)=9.03 +- 0.21 and cz/(r_s(z_d)H(z)) = 12.14 +- 0.43 at the mean redsfhit of the sample, z=0.57. We combine the information from the clustering wedges with recent measurements of CMB, BAO and type Ia supernovae to obtain constraints on the cosmological parameters of the standard LCDM model and a number of potential extensions. The information encoded in the clustering wedges is most useful when the dark energy equation of state is allowed to deviate from its standard LCDM value. The combination of all datasets shows no evidence of a deviation from a constant dark energy equation of state, in which case we find w_DE = -1.013 +- 0.064, in complete agreement with a cosmological constant. We explore potential deviations from general relativity by constraining the growth rate f(z)=d ln D(a)/ d ln a, in which case the combination of the CMASS clustering wedges with CMB data implies f(z=0.57)=0.719 +- 0.094, in accordance with the predictions of GR. Our results clearly illustrate the additional constraining power of anisotropic clustering measurements with respect to that of angle-averaged quantities.

Galaxy redshift surveys are becoming increasingly important as a dark energy probe. We improve the forecasting of dark energy constraints from galaxy redshift surveys by using the "dewiggled" galaxy power spectrum, P_{dw}(k), in the Fisher matrix calculations. Since P_{dw}(k) is a good fit to real galaxy clustering data over most of the scale range of interest, our approach is more realistic compared to previous work in forecasting dark energy constraints from galaxy redshift surveys. We find that our new approach gives results in excellent agreement when compared to the results from the actual data analysis of the clustering of the Sloan Digital Sky Survey DR7 luminous red galaxies. We provide forecasts of the dark energy constraints from a plausible Stage IV galaxy redshift survey.

We present a method for measuring the Hubble parameter, H(z), and angular diameter distance, D_A(z), from the two-dimensional two-point correlation function, and validate it using LasDamas mock galaxy catalogs. Applying our method to the sample of luminous red galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we measure H(z=0.35)=82.1_{-4.9}^{+4.8} km s^{-1}Mpc^{-1}, D_A(z=0.35)=1048_{-58}^{+60} Mpc without assuming a dark energy model or a flat Universe. We find that the derived measurements of H(0.35)r_s(z_d)/c and D_A(0.35)/r_s(z_d) (where r_s(z_d) is the sound horizon at the drag epoch) are nearly uncorrelated, have tighter constraints and are more robust with respect to possible systematic effects. Our galaxy clustering measurements of {H(0.35)r_s(z_d)/c, D_A(0.35)/r_s(z_d)}={0.0434\pm 0.0018,6.60\pm 0.26} (with the correlation coefficient r = 0.0604) can be used to combine with cosmic microwave background and any other cosmological data sets to constrain dark energy. Our results represent the first measurements of H(z) and D_A(z) (or H(z)r_s(z_d)/c and D_A(0.35)/r_s(z_d)) from galaxy clustering data. Our work has significant implications for future surveys in establishing the feasibility of measuring both H(z) and D_A(z) from galaxy clustering data.

We present a comprehensive analysis of a bright, long duration (T90 ~ 257 s) GRB 110205A at redshift z= 2.22. The optical prompt emission was detected by Swift/UVOT, ROTSE-IIIb and BOOTES telescopes when the GRB was still radiating in the gamma-ray band. Nearly 200 s of observations were obtained simultaneously from optical, X-ray to gamma-ray, which makes it one of the exceptional cases to study the broadband spectral energy distribution across 6 orders of magnitude in energy during the prompt emission phase. By fitting the time resolved prompt spectra, we clearly identify, for the first time, an interesting two-break energy spectrum, roughly consistent with the standard GRB synchrotron emission model in the fast cooling regime. Although the prompt optical emission is brighter than the extrapolation of the best fit X/gamma-ray spectra, it traces the gamma-ray light curve shape, suggesting a relation to the prompt high energy emission. The synchrotron + SSC scenario is disfavored by the data, but the models invoking a pair of internal shocks or having two emission regions can interpret the data well. Shortly after prompt emission (~ 1100 s), a bright (R = 14.0) optical emission hump with very steep rise (alpha ~ 5.5) was observed which we interpret as the emission from the reverse shock. It is the first time that the rising phase of a reverse shock component has been closely observed. The full optical and X-ray afterglow lightcurves can be interpreted within the standard reverse shock (RS) + forward shock (FS) model. In general, the high quality prompt emission and afterglow data allow us to apply the standard fireball shock model to extract valuable information about the GRB including the radiation mechanism, radius of prompt emission R, initial Lorentz factor of the outflow, the composition of the ejecta, as well as the collimation angle and the total energy budget.

We measure the effective distance to z=0.35, D_V(0.35), from the overall shape of the spherically-averaged two-point correlation function of the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) luminous red galaxy (LRG) sample. We find D_V(0.35)=1428_{-73}^{+74} without assuming a dark energy model or a flat Universe. We find that the derived measurement of r_s(z_d)/D_V(0.35)=0.1143 \pm 0.0030 (the ratio of the sound horizon at the drag epoch to the effective distance to z=0.35) is more tightly constrained and more robust with respect to possible systematic effects. It is also nearly uncorrelated with \Omega_m h^2. Combining our results with the cosmic microwave background and supernova data, we obtain \Omega_k=-0.0032^{+0.0074}_{-0.0072} and w=-1.010^{+0.046}_{-0.045} (assuming a constant dark energy equation of state). By scaling the spherically-averaged correlation function, we find the Hubble parameter H(0.35)=83^{+13}_{-15} km s^{-1}Mpc^{-1} and the angular diameter distance D_A(0.35)=1089^{+93}_{-87} Mpc. We use LasDamas SDSS mock catalogs to compute the covariance matrix of the correlation function, and investigate the use of lognormal catalogs as an alternative. We find that the input correlation function can be accurately recovered from lognormal catalogs, although they give larger errors on all scales (especially on small scales) compared to the mock catalogs derived from cosmological N-body simulations.

We examine how dark energy constraints from current observational data depend on the analysis methods used: the analysis of Type Ia supernovae (SNe Ia), and that of galaxy clustering data. We generalize the flux-averaging analysis method of SNe Ia to allow correlated errors of SNe Ia, in order to reduce the systematic bias due to weak lensing of SNe Ia. We find that flux-averaging leads to larger errors on dark energy and cosmological parameters if only SN Ia data are used. When SN Ia data (the latest compilation by the SNLS team) are combined with WMAP 7 year results (in terms of our Gaussian fits to the probability distributions of the CMB shift parameters), the latest Hubble constant (H_0) measurement using the Hubble Space Telescope (HST), and gamma ray burst (GRB) data, flux-averaging of SNe Ia increases the concordance with other data, and leads to significantly tighter constraints on the dark energy density at z=1, and the cosmic curvature \Omega_k. The galaxy clustering measurements of H(z=0.35)r_s(z_d) and r_s(z_d)/D_A(z=0.35) (where H(z) is the Hubble parameter, D_A(z) is the angular diameter distance, and r_s(z_d) is the sound horizon at the drag epoch) by Chuang & Wang (2011) are consistent with SN Ia data, given the same pirors (CMB+H_0+GRB), and lead to significantly improved dark energy constraints when combined. Current data are fully consistent with a cosmological constant and a flat universe.

It is the common consensus that the expansion of a universe always slows down if the gravity provided by the energy sources therein is attractive and accordingly one needs to invoke dark energy as a source of anti-gravity for understanding the cosmic acceleration. To examine this point we find counter-examples for a spherically symmetric dust fluid described by the Lemaitre-Tolman-Bondi solution without singularity. Thus, the validity of this naive consensus is indeed doubtful and the effects of inhomogeneities should be restudied. These counter-intuitive examples open a new perspective on the understanding of the evolution of our universe.