We present the GEOMAX algorithm and its Python implementation for a two-step compression of bispectrum measurements. The first step groups bispectra by the geometric properties of their arguments; the second step then maximises the Fisher information with respect to a chosen set of model parameters in each group. The algorithm only requires the derivatives of the data vector with respect to the parameters and a small number of mock data, producing an effective, non-linear compression. By applying GEOMAX to bispectrum monopole measurements from BOSS DR12 CMASS redshift-space galaxy clustering data, we reduce the $68\%$ credible intervals for the inferred parameters $\left(b_1,b_2,f,\sigma_8\right)$ by $\left(50.4\%,56.1\%,33.2\%,38.3\%\right)$ with respect to standard MCMC on the full data vector. We run the analysis and comparison between compression methods over one hundred galaxy mocks to test the statistical significance of the improvements. On average GEOMAX performs $\sim15\%$ better than geometrical or maximal linear compression alone and is consistent with being lossless. Given its flexibility, the GEOMAX approach has the potential to optimally exploit three-point statistics of various cosmological probes like weak lensing or line-intensity maps from current and future cosmological data-sets such as DESI, Euclid, PFS and SKA.

Detecting primordial non-Gaussianity on mildly non-linear scales requires precise modelling of late-time structure formation. Accurately predicting the impact of non-linear gravitational collapse, non-linear tracer bias and baryonic physics on the variance and higher order moments of the cosmic density field is challenging, as they strongly depend on the tails of the probability distribution function (PDF) of density fluctuations. A way around this problem is to directly analyse the bulk of the PDF instead. For this purpose we devise a new method to predict the impact of general non-Gaussian initial conditions on the late-time density PDF. With this formalism we show that - even when marginalizing over potential ignorance of the amplitude and slope of the non-linear power spectrum - an analysis of the PDF at mildly non-linear densities can measure the amplitude of different primordial bispectrum shapes to an accuracy of $\Delta f_{\mathrm{NL}}^{\mathrm{loc}} = \pm 3.1\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{equi}} = \pm 10.0\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{ortho}} = \pm 17.0\ $. This assumes a joint analysis of the PDF on smoothing scales of $15$Mpc/$h$ and $30$Mpc/$h$ in a survey volume of $V=100(\mathrm{Gpc}/h)^3$ at $z=1$, analysing only densities of $\delta(15\mathrm{Mpc}/h) \in [-0.4, 0.5]$ ($\approx 87\%$ of probability) and $\delta(30\mathrm{Mpc}/h) \in [-0.3, 0.4]$ ($\approx 95\%$ of probability). Note that a formalism closely related to ours was already successfully applied to observational data \citep{Gruen2018, Friedrich2018}, demonstrating that the methodology developed here can indeed be carried over to real data analysis.

We investigate various phenomenological schemes for the rapid generation of 3D mock galaxy catalogues with a given power spectrum and bispectrum. We apply the fast bispectrum estimator \MODALLSS{} to these mock galaxy catalogues and compare to $N$-body simulation data analysed with the halo-finder \texttt{ROCKSTAR} (our benchmark data). We propose an assembly bias model for populating parent halos with subhalos by using a joint lognormal-Gaussian probability distribution for the subhalo occupation number and the halo concentration. This prescription enabled us to recover the benchmark power spectrum from $N$-body simulations to within 1\% and the bispectrum to within 4\% across the entire range of scales of the simulation. A small further boost adding an extra galaxy to all parent halos above the mass threshold $M>2\times10^{14}\,h^{-1} M_\odot$ obtained a better than 1\% fit to both power spectrum and bispectrum in the range $K/3<1.1\,h\,\text{Mpc}^{-1}$, where $K=k_1+k_2+k_3$. This statistical model should be applicable to fast dark matter codes, allowing rapid generation of mock catalogues which simultaneously reproduce the halo power spectrum and bispectrum obtained from $N$-body simulations. We also investigate alternative schemes using the Halo Occupation Distribution (HOD) which depend only on halo mass, but these yield results deficient in both the power spectrum (2\%) and the bispectrum (>4\%) at $k,K/3 \approx 0.2\,h\,\text{Mpc}^{-1}$, with poor scaling for the latter. Efforts to match the power spectrum by modifying the standard four-parameter HOD model result in overboosting the bispectrum (with a 10\% excess). We also characterise the effect of changing the halo profile on the power spectrum and bispectrum.

We investigate the impact of prior models on the upper bound of the sum of neutrino masses, $\sum m_{\nu}$. We use data from Large Scale Structure of galaxies, Cosmic Microwave Background, Type Ia SuperNovae, and Big Bang Nucleosynthesis. We probe physically motivated neutrino mass models (respecting oscillation experiment constraints) and compare them to constraints using standard cosmological approximations. The former give a consistent upper bound of $\sum m_{\nu} \lesssim 0.26$ eV ($95\%$ CI) and yields a strong competitive upper bound for the lightest neutrino mass species, $m_0^{\nu} < 0.086$ eV ($95\%$ CI). By contrast one of the approximations, which is somewhat inconsistent with oscillation experiments, yields an upper bound of $\sum m_{\nu} \lesssim 0.15$ eV ($95\%$ CI), which differs substantially from the former upper bound. We, therefore, argue that cosmological neutrino mass and hierarchy determination should be pursued using physically motivated models since approximations might lead to incorrect and nonphysical upper bounds.

We constrain cosmological parameters by analysing the angular power spectra of the Baryon Oscillation Spectroscopic Survey DR12 galaxies, a spectroscopic follow-up of around 1.3 million SDSS galaxies over 9,376 deg$^2$ with an effective volume of $\sim 6.5$ (Gpc $h^{-1}$)$^3$ in the redshift range $0.15 \leq z < 0.80$. We split this sample into 13 tomographic bins ($\Delta z = 0.05$); angular power spectra were calculated using a Pseudo-$C_{\ell}$ estimator, and covariance matrices were estimated using log-normal simulated maps. Cosmological constraints obtained from these data were combined with constraints from Planck CMB experiment as well as the JLA supernovae compilation. Considering a $w$CDM cosmological model measured on scales up to $k_{max} = 0.07h$ Mpc$^{-1}$, we constrain a constant dark energy equation-of-state with a $\sim 4\%$ error at the 1-$\sigma$ level: $w_0 = -0.993^{+0.046}_{-0.043}$, together with $\Omega_m = 0.330\pm 0.012$, $\Omega_b = 0.0505 \pm 0.002$, $S_8 \equiv \sigma_8 \sqrt{\Omega_m/0.3} = 0.863 \pm 0.016$, and $h = 0.661 \pm 0.012$. For the same combination of datasets, but now considering a $\Lambda$CDM model with massive neutrinos and the same scale cut, we find: $\Omega_m = 0.328 \pm 0.009$, $\Omega_b = 0.05017^{+0.0009}_{-0.0008}$, $S_8 = 0.862 \pm 0.017$, and $h = 0.663^{+0.006}_{-0.007}$ and a 95\% credible interval (CI) upper limit of $\sum m_{\nu} < 0.14$ eV for a normal hierarchy. These results are competitive if not better than standard analyses with the same dataset, and demonstrate this should be a method of choice for future surveys, opening the door for their full exploitation in cross-correlations probes.

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

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

We present a novel method to compress galaxy clustering three-point statistics and apply it to redshift space galaxy bispectrum monopole measurements from BOSS DR12 CMASS data considering a $k$-space range of $0.03-0.12\,h/\mathrm{Mpc}$. The method consists in binning together bispectra evaluated at sets of wave-numbers forming closed triangles with similar geometrical properties: the area, the cosine of the largest angle and the ratio between the cosines of the remaining two angles. This enables us to increase the number of bispectrum measurements for example by a factor of $23$ over the standard binning (from 116 to 2734 triangles used), which is otherwise limited by the number of mock catalogues available to estimate the covariance matrix needed to derive parameter constraints. The $68\%$ credible intervals for the inferred parameters $\left(b_1,b_2,f,\sigma_8\right)$ are thus reduced by $\left(-39\%,-49\%,-29\%,-22\%\right)$, respectively. We find very good agreement with the posteriors recently obtained by alternative maximal compression methods. This new method does not require the a-priori computation of the data-vector covariance matrix and has the potential to be directly applicable to other three-point statistics (e.g. glaxy clustering, weak gravitational lensing, 21 cm emission line) measured from future surveys such as DESI, Euclid, PFS and SKA.

We apply two compression methods to the galaxy power spectrum monopole/quadrupole and bispectrum monopole measurements from the BOSS DR12 CMASS sample. Both methods reduce the dimension of the original data-vector to the number of cosmological parameters considered, using the Karhunen-Lo\`eve algorithm with an analytic covariance model. In the first case, we infer the posterior through MCMC sampling from the likelihood of the compressed data-vector (MC-KL). The second, faster option, works by first Gaussianising and then orthogonalising the parameter space before the compression; in this option (G-PCA) we only need to run a low-resolution preliminary MCMC sample for the Gaussianization to compute our posterior. Both compression methods accurately reproduce the posterior distributions obtained by standard MCMC sampling on the CMASS dataset for a $k$-space range of $0.03-0.12\,h/\mathrm{Mpc}$. The compression enables us to increase the number of bispectrum measurements by a factor of $\sim 23$ over the standard binning (from 116 to 2734 triangles used), which is otherwise limited by the number of mock catalogues available. This reduces the $68\%$ credible intervals for the parameters $\left(b_1,b_2,f,\sigma_8\right)$ by $\left(-24.8\%,-52.8\%,-26.4\%,-21\%\right)$, respectively. The best-fit values we obtain are $(b_1=2.31\pm0.17,b_2=0.77\pm0.19,$ $f(z_{\mathrm{CMASS}})=0.67\pm0.06,\sigma_8(z_{\mathrm{CMASS}})=0.51\pm0.03)$. Using these methods for future redshift surveys like DESI, Euclid and PFS will drastically reduce the number of simulations needed to compute accurate covariance matrices and will facilitate tighter constraints on cosmological parameters.

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

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

We explore two methods of compressing the redshift space galaxy power spectrum and bispectrum with respect to a chosen set of cosmological parameters. Both methods involve reducing the dimension of the original data-vector ( e.g. 1000 elements ) to the number of cosmological parameters considered ( e.g. seven ) using the Karhunen-Lo\`eve algorithm. In the first case, we run MCMC sampling on the compressed data-vector in order to recover the one-dimensional (1D) and two-dimensional (2D) posterior distributions. The second option, approximately 2000 times faster, works by orthogonalising the parameter space through diagonalisation of the Fisher information matrix before the compression, obtaining the posterior distributions without the need of MCMC sampling. Using these methods for future spectroscopic redshift surveys like DESI, EUCLID and PFS would drastically reduce the number of simulations needed to compute accurate covariance matrices with minimal loss of constraining power. We consider a redshift bin of a DESI-like experiment. Using the power spectrum combined with the bispectrum as a data-vector, both compression methods on average recover the 68% credible regions to within 0.7% and 2% of those resulting from standard MCMC sampling respectively. These confidence intervals are also smaller than the ones obtained using only the power spectrum by (81%, 80%, 82%) respectively for the bias parameter b_1, the growth rate f and the scalar amplitude parameter A_s.

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

We present measurements of the 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 general parallelized and easy-to-use code to perform numerical simulations of structure formation using the COLA (COmoving Lagrangian Acceleration) method for cosmological models that exhibit scale-dependent growth at the level of first and second order Lagrangian perturbation theory. For modified gravity theories we also include screening using a fast approximate method that covers all the main examples of screening mechanisms in the literature. We test the code by comparing it to full simulations of two popular modified gravity models, namely $f(R)$ gravity and nDGP, and find good agreement in the modified gravity boost-factors relative to $\Lambda$CDM even when using a fairly small number of COLA time steps.

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

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

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

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.

We present a new method to measure the redshift-dependent galaxy bias by combining information from the galaxy density field and the weak lensing field. This method is based on Amara et al. (2012), where they use the galaxy density field to construct a bias-weighted convergence field kg. The main difference between Amara et al. (2012) and our new implementation is that here we present another way to measure galaxy bias using tomography instead of bias parameterizations. The correlation between kg and the true lensing field k allows us to measure galaxy bias using different zero-lag correlations, such as <kgk>/<kk> or <kgkg>/<kgk>. Our method measures the linear bias factor on linear scales under the assumption of no stochasticity between galaxies and matter. We use the MICE simulation to measure the linear galaxy bias for a flux-limited sample (i < 22.5) in tomographic redshift bins using this method. This paper is the first that studies the accuracy and systematic uncertainties associated with the implementation of the method, and the regime where it is consistent with the linear galaxy bias defined by projected 2-point correlation functions (2PCF). We find that our method is consistent with linear bias at the percent level for scales larger than 30 arcmin, while nonlinearities appear at smaller scales. This measurement is a good complement to other measurements of bias, since it does not depend strongly on sigma8 as the 2PCF measurements. We apply this method to the Dark Energy Survey Science Verification data in a follow-up paper.

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

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)

Robust measurements based on current large-scale structure surveys require precise knowledge of statistical and systematic errors. This can be obtained from large numbers of realistic mock galaxy catalogues that mimic the observed distribution of galaxies within the survey volume. To this end we present a fast, distributed-memory, planar-parallel code, L-PICOLA, which can be used to generate and evolve a set of initial conditions into a dark matter field much faster than a full non-linear N-Body simulation. Additionally, L-PICOLA has the ability to include primordial non-Gaussianity in the simulation and simulate the past lightcone at run-time, with optional replication of the simulation volume. Through comparisons to fully non-linear N-Body simulations we find that our code can reproduce the $z=0$ power spectrum and reduced bispectrum of dark matter to within 2% and 5% respectively on all scales of interest to measurements of Baryon Acoustic Oscillations and Redshift Space Distortions, but 3 orders of magnitude faster. The accuracy, speed and scalability of this code, alongside the additional features we have implemented, make it extremely useful for both current and next generation large-scale structure surveys. L-PICOLA is publicly available at https://cullanhowlett.github.io/l-picola

We examine the cosmological implications of the measurements of the linear growth rate of cosmological structure obtained in a companion paper from the power spectrum and bispectrum monopoles of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey Data, Release 11, CMASS galaxies. This measurement was of $f^{0.43}\sigma_8$, where $\sigma_8$ is the amplitude of dark matter density fluctuations, and $f$ is the linear growth rate, at the effective redshift of the survey, $z_{\rm eff}=0.57$. In conjunction with Cosmic Microwave Background (CMB) data, interesting constraints can be placed on models with non-standard neutrino properties and models where gravity deviates from general relativity on cosmological scales. In particular, the sum of the masses of the three species of the neutrinos is constrained to $m_\nu<0.49\,{\rm eV}$ (at 95\% confidence level) when the $f^{0.43}\sigma_8$ measurement is combined with state-of-the-art CMB measurements. Allowing the effective number of neutrinos to vary as a free parameter does not significantly change these results. When we combine the measurement of $f^{0.43}\sigma_8$ with the complementary measurement of $f\sigma_8$ from the monopole and quadrupole of the two-point correlation function we are able to obtain an independent measurements of $f$ and $\sigma_8$. We obtain $f=0.63\pm0.16$ and $\sigma_8=0.710\pm0.086$ (68\% confidence level). This is the first time when these parameters have been able to be measured independently using the redshift-space power spectrum and bispectrum measurements from galaxy clustering data only.

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.

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 analyse the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 sample, which consists of $690 827$ galaxies in the redshift range $0.43 < z < 0.70$ and has a sky coverage of $8 498$ deg$^2$ corresponding to an effective volume of $\sim6\,\rm{Gpc}^3$. We fit the Fourier space statistics, the power spectrum and bispectrum monopoles to measure the linear and quadratic bias parameters, $b_1$ and $b_2$, for a non-linear non-local bias model, the growth of structure parameter $f$ and the amplitude of dark matter density fluctuations parametrised by $\sigma_8$. We obtain $b_1(z_{\rm eff})^{1.40}\sigma_8(z_{\rm eff})=1.672\pm 0.060$ and $b_2^{0.30}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.579\pm0.082$ at the effective redshift of the survey, $z_{\rm eff}=0.57$. The main cosmological result is the constraint on the combination $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.582\pm0.084$, which is complementary to $f\sigma_8$ constraints obtained from 2-point redshift space distortion analyses. A less conservative analysis yields $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.584\pm0.051$. We ensure that our result is robust by performing detailed systematic tests using a large suite of survey galaxy mock catalogs and N-body simulations. The constraints on $f^{0.43}\sigma_8$ are useful for setting additional constrains on neutrino mass, gravity, curvature as well as the number of neutrino species from galaxy surveys analyses (as presented in a companion paper).

Anisotropic measurements of the Baryon Acoustic Oscillation (BAO) feature within a galaxy survey enable joint inference about the Hubble parameter $H(z)$ and angular diameter distance $D_A(z)$. These measurements are typically obtained from moments of the measured 2-point clustering statistics, with respect to the cosine of the angle to the line of sight $\mu$. The position of the BAO features in each moment depends on a combination of $D_A(z)$ and $H(z)$, and measuring the positions in two or more moments breaks this parameter degeneracy. We derive analytic formulae for the parameter combinations measured from moments given by Legendre polynomials, power laws and top-hat Wedges in $\mu$, showing explicitly what is being measured by each in real-space for both the correlation function and power spectrum, and in redshift-space for the power spectrum. The large volume covered by modern galaxy samples means that the correlation function can be well approximated as having no correlations at different $\mu$ on the BAO scale, and that the errors on this scale are approximately independent of $\mu$. Using these approximations, we derive the information content of various moments. We show that measurements made using either the monopole and quadrupole, or the monopole and $\mu^2$ power-law moment, are optimal for anisotropic BAO measurements, in that they contain all of the available information using two moments, the minimal number required to measure both $H(z)$ and $D_A(z)$. We test our predictions using 600 mock galaxy samples, matched to the SDSS-III Baryon Oscillation Spectroscopic Survey CMASS sample, finding a good match to our analytic predictions. Our results should enable the optimal extraction of information from future galaxy surveys such as eBOSS, DESI and Euclid.

In this work we study the relevance of the cosmic web and substructures on the matter and lensing power spectra measured from halo mock catalogues extracted from the N-body simulations. Since N-body simulations are computationally expensive, it is common to use faster methods that approximate the dark matter field as a set of halos. In this approximation, we replace mass concentrations in N-body simulations by a spherically symmetric Navarro-Frenk-White halo density profile. We also consider the full mass field as the sum of two distinct fields: dark matter halos ($M>9\times 10^{12}~M_{\odot}$/h) and particles not included into halos. Mock halos reproduce well the matter power spectrum, but underestimate the lensing power spectrum on large and small scales. For sources at $z_{\rm s}=1$ the lensing power spectrum is underestimated by up to 40% at $\ell\approx 10^4$ with respect to the simulated halos. The large scale effect can be alleviated by combining the mock catalogue with the dark matter distribution outside the halos. In addition, to evaluate the contribution of substructures we have smeared out the intra-halo substructures in a N-body simulation while keeping the halo density profiles unchanged. For the matter power spectrum the effect of this smoothing is only of the order of 5%, but for lensing substructures are much more important: for $\ell\approx 10^4$ the internal structures contribute 30% of the total spectrum. These findings have important implications in the way mock catalogues have to be created, suggesting that some approximate methods currently used for galaxy surveys will be inadequate for future weak lensing surveys.

Reconstructing an estimate of linear Baryon Acoustic Oscillations (BAO) from an evolved galaxy field has become a standard technique in recent analyses. By partially removing non-linear damping caused by bulk motions, the real-space BAO peak in the correlation function is sharpened, and oscillations in the power spectrum are visible to smaller scales. In turn these lead to stronger measurements of the BAO scale. Future surveys are being designed assuming that this improvement has been applied, and this technique is therefore of critical importance for future BAO measurements. A number of reconstruction techniques are available, but the most widely used is a simple algorithm that decorrelates large-scale and small-scale modes approximately removing the bulk-flow displacements by moving the overdensity field (Eisenstein et al. 2007; Padmanabhan, White & Cohn 2009). We consider the practical implementation of this algorithm, looking at the efficiency of reconstruction as a function of the assumptions made for the bulk-flow scale, the shot noise level in a random catalogue used to quantify the mask and the method used to estimate the bulk-flow shifts. We also examine the efficiency of reconstruction against external factors including galaxy density, volume and edge effects, and consider their impact for future surveys. Throughout we make use of the mocks catalogues created for the Baryon Oscillation Spectroscopic Survey (BOSS) Date Release 11 samples covering 0.43 < z < 0.7 (CMASS) and 0.15 < z < 0.43 (LOWZ), to empirically test these changes.

We measure Redshift-Space Distortions (RSD) in the two-point correlation function of a sample of $63,163$ spectroscopically identified galaxies with $z < 0.2$, an epoch where there are currently only limited measurements, from the Sloan Digital Sky Survey (SDSS) Data Release 7 Main Galaxy Sample. Our sample, which we denote MGS, covers 6,813 deg$^2$ with an effective redshift $z_{eff}=0.15$ and is described in our companion paper (Paper I), which concentrates on BAO measurements. In order to validate the fitting methods used in both papers, and derive errors, we create and analyse 1000 mock catalogues using a new algorithm called PICOLA to generate accurate dark matter fields. Haloes are then selected using a friends-of-friends algorithm, and populated with galaxies using a Halo-Occupation Distribution fitted to the data. Using errors derived from these mocks, we fit a model to the monopole and quadrupole moments of the MGS correlation function. If we assume no Alcock-Paczynski (AP) effect (valid at $z=0.15$ for any smooth model of the expansion history), we measure the amplitude of the velocity field, $f\sigma_{8}$, at $z=0.15$ to be $0.49_{-0.14}^{+0.15}$. We also measure $f\sigma_{8}$ including the AP effect. This latter measurement can be freely combined with recent Cosmic Microwave Background results to constrain the growth index of fluctuations, $\gamma$. Assuming a background $\Lambda$CDM cosmology and combining with current Baryon Acoustic Oscillation data we find $\gamma = 0.64 \pm 0.09$, which is consistent with the prediction of General Relativity ($\gamma \approx 0.55$), though with a slight preference for higher $\gamma$ and hence models with weaker gravitational interactions.

We create a sample of spectroscopically identified galaxies with $z < 0.2$ from the Sloan Digital Sky Survey (SDSS) Data Release 7, covering 6813 deg$^2$. Galaxies are chosen to sample the highest mass haloes, with an effective bias of 1.5, allowing us to construct 1000 mock galaxy catalogs (described in Paper II), which we use to estimate statistical errors and test our methods. We use an estimate of the gravitational potential to "reconstruct" the linear density fluctuations, enhancing the Baryon Acoustic Oscillation (BAO) signal in the measured correlation function and power spectrum. Fitting to these measurements, we determine $D_{V}(z_{\rm eff}=0.15) = (664\pm25)(r_d/r_{d,{\rm fid}})$ Mpc; this is a better than 4 per cent distance measurement. This "fills the gap" in BAO distance ladder between previously measured local and higher redshift measurements, and affords significant improvement in constraining the properties of dark energy. Combining our measurement with other BAO measurements from BOSS and 6dFGS galaxy samples provides a 15 per cent improvement in the determination of the equation of state of dark energy and the value of the Hubble parameter at $z=0$ ($H_0$). Our measurement is fully consistent with the Planck results and the $\Lambda$CDM concordance cosmology, but increases the tension between Planck$+$BAO $H_0$ determinations and direct $H_0$ measurements.

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.

We use the observed anisotropic clustering of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 CMASS sample to measure the linear growth rate of structure, the Hubble expansion rate and the comoving distance scale. Our sample covers 8498 ${\rm deg}^2$ and encloses an effective volume of 6.0 ${\rm Gpc}^3$ at an effective redshift of $\bar{z} = 0.57$. We find $f\sigma_8 = 0.441 \pm 0.044$, $H = 93.1 \pm 3.0\ {\mathrm{km}\ \mathrm{s}^{-1} \mathrm{Mpc}^{-1}}$ and $D_{\rm A} = 1380 \pm 23\ {\rm Mpc}$ when fitting the growth and expansion rate simultaneously. When we fix the background expansion to the one predicted by spatially-flat $\Lambda$CDM model in agreement with recent Planck results, we find $f\sigma_8 = 0.447 \pm 0.028$ (6 per cent accuracy). While our measurements are generally consistent with the predictions of $\Lambda$CDM and General Relativity, they mildly favor models in which the strength of gravitational interactions is weaker than what is predicted by General Relativity. Combining our measurements with recent cosmic microwave background data results in tight constraints on basic cosmological parameters and deviations from the standard cosmological model. Separately varying these parameters, we find $w = -0.983 \pm 0.075$ (8 per cent accuracy) and $\gamma = 0.69 \pm 0.11$ (16 per cent accuracy) for the effective equation of state of dark energy and the growth rate index, respectively. Both constraints are in good agreement with the standard model values of $w=-1$ and $\gamma = 0.554$.

We use cosmological N-body simulations to investigate whether measurements of the moments of large-scale structure can yield constraints on primordial non-Gaussianity. We measure the variance, skewness, and kurtosis of the evolved density field from simulations with Gaussian and three different non-Gaussian initial conditions: a local model with f_NL=100, an equilateral model with f_NL=-400, and an orthogonal model with f_NL=-400. We show that the moments of the dark matter density field differ significantly between Gaussian and non-Gaussian models. We also make the measurements on mock galaxy catalogs that contain galaxies with clustering properties similar to those of luminous red galaxies (LRGs). We find that, in the case of skewness and kurtosis, galaxy bias reduces the detectability of non-Gaussianity, though we can still clearly discriminate between different models in our simulation volume. However, in the case of the variance, galaxy bias greatly amplifies the detectability of non-Gaussianity. In all cases we find that redshift distortions do not significantly affect the detectability. When we restrict our measurements to volumes equivalent to the Sloan Digital Sky Survey II (SDSS-II) or Baryon Oscillation Spectroscopic Survey (BOSS) samples, the probability of detecting a departure from the Gaussian model is high by using measurements of the variance, but very low by using only skewness and kurtosis measurements. We find that skewness and kurtosis measurements are never likely to yield useful constraints on primordial non-Gaussianity, but future surveys should be large enough to place meaningful constraints using measurements of the galaxy variance.(Abridged)

We measure the sum of the neutrino particle masses using the three-dimensional galaxy power spectrum of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) CMASS galaxy sample. Combined with the cosmic microwave background (CMB), supernova (SN) and additional baryonic acoustic oscillation (BAO) data, we find upper 95 percent confidence limits of the neutrino mass $\Sigma m_{\nu}<0.340$ eV within a flat $\Lambda$CDM background, and $\Sigma m_{\nu}<0.821$ eV, assuming a more general background cosmological model. The number of neutrino species is measured to be $N_{\rm eff}=4.308\pm0.794$ and $N_{\rm eff}=4.032^{+0.870}_{-0.894}$ for these two cases respectively. We study and quantify the effect of several factors on the neutrino measurements, including the galaxy power spectrum bias model, the effect of redshift-space distortion, the cutoff scale of the power spectrum, and the choice of additional data. The impact of neutrinos with unknown masses on other cosmological parameter measurements is investigated. The fractional matter density and the Hubble parameter are measured to be $\Omega_M=0.2796\pm0.0097$, $H_0=69.72^{+0.90}_{-0.91}$ km/s/Mpc (flat $\Lambda$CDM) and $\Omega_M=0.2798^{+0.0132}_{-0.0136}$, $H_0=73.78^{+3.16}_{-3.17}$ km/s/Mpc (more general background model). Based on a Chevallier-Polarski-Linder (CPL) parametrisation of the equation-of-state $w$ of dark energy, we find that $w=-1$ is consistent with observations, even allowing for neutrinos. Similarly, the curvature \Omega_K and the running of the spectral index $\alpha_s$ are both consistent with zero. The tensor-to-scaler ratio is constrained down to $r<0.198$ (95 percent CL, flat $\Lambda$ CDM) and $r<0.440$ (95 percent CL, more general background model).

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 one thousand mock galaxy catalogues for the analysis of the Low Redshift Sample (LOWZ, effective redshift z ~ 10.32) of the Baryon Oscillation Spectroscopic Survey Data Releases 10 and 11. These mocks have been created following the PTHalos method of Manera13 et al. (2013) revised to include new developments. The main improvement is the introduction of a redshift dependence in the Halo Occupation Distribution in order to account for the change of the galaxy number density with redshift. These mock catalogues are used in the analyses of the LOWZ galaxy clustering by the BOSS collaboration.

We present the distance measurement to z = 0.32 using the 11th data release of the Sloan Digital Sky Survey-III Baryon Acoustic Oscillation Survey (BOSS). We use 313,780 galaxies of the low-redshift (LOWZ) sample over 7,341 square-degrees to compute $D_V = (1264 \pm 25)(r_d/r_{d,fid})$ - a sub 2% measurement - using the baryon acoustic feature measured in the galaxy two-point correlation function and power-spectrum. We compare our results to those obtained in DR10. We study observational systematics in the LOWZ sample and quantify potential effects due to photometric offsets between the northern and southern Galactic caps. We find the sample to be robust to all systematic effects found to impact on the targeting of higher-redshift BOSS galaxies, and that the observed north-south tensions can be explained by either limitations in photometric calibration or by sample variance, and have no impact on our final result. Our measurement, combined with the baryonic acoustic scale at z = 0.57, is used in Anderson et al. (2013a) to constrain cosmological parameters.

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.

Extraction of the Baryon Acoustic Oscillations (BAO) to percent level accuracy is challenging and demands an understanding of many potential systematic to an accuracy well below 1 per cent, in order ensure that they do not combine significantly when compared to statistical error of the BAO measurement. Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) SDSS Data Release Eleven (DR11) reaches a distance measurement with $\sim 1\%$ statistical error and this prompts an extensive search for all possible sub-percent level systematic errors which could be safely ignored previously. In this paper, we analyze the potential systematics in BAO fitting methodology using mocks and data from BOSS DR10 and DR11. We demonstrate the robustness of the fiducial multipole fitting methodology to be at $0.1\%-0.2\%$ level with a wide range of tests in mock galaxy catalogs pre- and post-reconstruction. We also find the DR10 and DR11 data from BOSS to be robust against changes in methodology at similar level. This systematic error budget is incorporated into the the error budget of Baryon Oscillation Spectroscopic Survey (BOSS) DR10 and DR11 BAO measurements. Of the wide range of changes we have investigated, we find that when fitting pre-reconstructed data or mocks, the following changes have the largest effect on the best fit values of distance measurements both parallel and perpendicular to the line of sight: (a) Changes in non-linear correlation function template; (b) Changes in fitting range of the correlation function; (c) Changes to the non-linear damping model parameters. The priors applied do not matter in the estimates of the fitted errors as long as we restrict ourselves to physically meaningful fitting regions.[abridged]

We explore the cosmological implications of the angle-averaged correlation function, xi(s), and the clustering wedges, xi_perp(s) and xi_para(s), of the LOWZ and CMASS galaxy samples from Data Release 10 and 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. Our results show no significant evidence for a deviation from the standard LCDM model. The combination of the information from our clustering measurements with recent data from the cosmic microwave background is sufficient to constrain the curvature of the Universe to Omega_k = 0.0010 +- 0.0029, the total neutrino mass to Sum m_nu < 0.23 eV (95% confidence level), the effective number of relativistic species to N_eff=3.31 +- 0.27, and the dark energy equation of state to w_DE = -1.051 +- 0.076. These limits are further improved by adding information from type Ia supernovae and baryon acoustic oscillations from other samples. In particular, this data set combination is completely consistent with a time-independent dark energy equation of state, in which case we find w_DE=-1.024 +- 0.052. We explore the constraints on the growth-rate of cosmic structures assuming f(z)=Omega_m(z)^gamma and obtain gamma=0.69 +- 0.15, in agreement with the predictions from general relativity of gamma=0.55.

We study the clustering of galaxies, as a function of their colour, from Data Release Ten (DR10) of the SDSS-III Baryon Oscillation Spectroscopic Survey. We select 122,967 galaxies with 0.43 < z < 0.7 into a "Blue" sample and 131,969 into a "Red" sample based on k+e corrected (to z=0.55) r-i colours and i band magnitudes. The samples are chosen to each contain more than 100,000 galaxies, have similar redshift distributions, and maximize the difference in clustering amplitude. The Red sample has a 40% larger bias than the Blue (b_Red/b_Blue = 1.39+-0.04), implying the Red galaxies occupy dark matter halos with an average mass that is 0.5 log Mo greater. Spherically averaged measurements of the correlation function, \xi 0, and the power spectrum are used to locate the position of the baryon acoustic oscillation (BAO) feature of both samples. Using \xi 0, we obtain distance scales, relative to our reference LCDM cosmology, of 1.010+-0.027 for the Red sample and 1.005+-0.031 for the Blue. After applying reconstruction, these measurements improve to 1.013+/-0.020 for the Red sample and 1.008+-0.026 for the Blue. For each sample, measurements of \xi 0 and the second multipole moment, \xi 2, of the anisotropic correlation function are used to determine the rate of structure growth, parameterized by f\sigma 8. We find f\sigma 8,Red = 0.511+-0.083, f\sigma 8,Blue = 0.509+/-0.085, and f\sigma 8,Cross = 0.423+-0.061 (from the cross-correlation between the Red and Blue samples). We use the covariance between the bias and growth measurements obtained from each sample and their cross-correlation to produce an optimally-combined measurement of f\sigma 8,comb = 0.443+-0.055. In no instance do we detect significant differences in distance scale or structure growth measurements obtained from the Blue and Red samples.

Measuring the two-point correlation function of the galaxies in the Universe gives access to the underlying dark matter distribution, which is related to cosmological parameters and to the physics of the primordial Universe. The estimation of the correlation function for current galaxy surveys makes use of the Landy-Szalay estimator, which is supposed to reach minimal variance. This is only true, however, for a vanishing correlation function. We study the Landy-Szalay estimator when these conditions are not fulfilled and propose a new estimator that provides the smallest variance for a given survey geometry. Our estimator is a linear combination of ratios between paircounts of data and/or random catalogues (DD, RR and DR). The optimal combination for a given geometry is determined by using lognormal mock catalogues. The resulting estimator is biased in a model-dependent way, but we propose a simple iterative procedure for obtaining an unbiased model- independent estimator.Our method can be easily applied to any dataset and requires few extra mock catalogues compared to the standard Landy-Szalay analysis. Using various sets of simulated data (lognormal, second-order LPT and N-Body), we obtain a 20-25% gain on the error bars on the two-point correlation function for the SDSS geometry and $\Lambda$CDM correlation function. When applied to SDSS data (DR7 and DR9), we achieve a similar gain on the correlation functions, which translates into a 10-15% improvement over the estimation of the densities of matter $\Omega_m$ and dark energy $\Omega_\Lambda$ in an open $\Lambda$CDM model. The constraints derived from DR7 data with our estimator are similar to those obtained with the DR9 data and the Landy-Szalay estimator, which covers a volume twice as large and has a density that is three times higher.

We present results on the clustering of 282,068 galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) sample of massive galaxies with redshifts 0.4<z<0.7 which is part of the Sloan Digital Sky Survey III project. Our results cover a large range of scales from ~0.5 to ~90 Mpc/h. We compare these estimates with the expectations of the flat LCDM cosmological model with parameters compatible with WMAP7 data. We use the MultiDark cosmological simulation together with a simple halo abundance matching technique, to estimate galaxy correlation functions, power spectra, abundance of subhaloes and galaxy biases. We find that the LCDM model gives a reasonable description to the observed correlation functions at z~0.5, which is a remarkably good agreement considering that the model, once matched to the observed abundance of BOSS galaxies, does not have any free parameters. However, we find a deviation (>~10%) in the correlation functions for scales less than ~1 Mpc/h and ~10-40 Mpc/h. A more realistic abundance matching model and better statistics from upcoming observations are needed to clarify the situation. We also estimate that about 12% of the "galaxies" in the abundance-matched sample are satellites inhabiting central haloes with mass M>~1e14 M_sun/h. Using the MultiDark simulation we also study the real space halo bias b(r) of the matched catalogue finding that b=2.00+/-0.07 at large scales, consistent with the one obtained using the measured BOSS projected correlation function. Furthermore, the linear large-scale bias depends on the number density n of the abundance-matched sample as b=-0.048-(0.594+/-0.02)*log(n/(h/Mpc)^3). Extrapolating these results to BAO scales we measure a scale-dependent damping of the acoustic signal produced by non-linear evolution that leads to ~2-4% dips at ~3 sigma level for wavenumbers k>~0.1 h/Mpc in the linear large-scale bias.

We measure the luminosity and color dependence and the redshift evolution of galaxy clustering in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Ninth Data Release. We focus on the projected two-point correlation function (2PCF) of subsets of its CMASS sample, which includes about 260,000 galaxies over ~3,300 sq. deg in the redshift range 0.43<z<0.7. To minimize the selection effect on galaxy clustering, we construct well-defined luminosity and color subsamples by carefully accounting for the CMASS galaxy selection cuts. The 2PCF of the whole CMASS sample, if approximated by a power-law, has a correlation length of r_0=7.93\pm0.06Mpc/h and an index of \gamma=1.85\pm0.01. Clear dependences on galaxy luminosity and color are found for the projected 2PCF in all redshift bins, with more luminous and redder galaxies generally exhibiting stronger clustering and steeper 2PCF. The color dependence is also clearly seen for galaxies within the red sequence, consistent with the behavior of SDSS-II main sample galaxies at lower redshifts. At a given luminosity (k+e corrected), no significant evolution of the projected 2PCFs with redshift is detected for red sequence galaxies. We also construct galaxy samples of fixed number density at different redshifts, using redshift-dependent magnitude thresholds. The clustering of these galaxies in the CMASS redshift range is found to be consistent with that predicted by passive evolution. Our measurements of the luminosity and color dependence and redshift evolution of galaxy clustering will allow for detailed modeling of the relation between galaxies and dark matter halos and new constraints on galaxy formation and evolution.

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

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 use the joint measurement of geometry and growth from anisotropic galaxy clustering in the Baryon Oscillation Spectroscopic Survey Data Release 9 CMASS sample reported by Reid et al. to constrain dark energy properties and possible deviations from the General Relativity. Assuming GR and taking a prior on the linear matter power spectrum at high redshift from the cosmic microwave background (CMB), anisotropic clustering of the CMASS DR9 galaxies alone constrains $\Omega_{\rm m} = 0.308 \pm 0.022$ and $100\Omega_{\rm k} = 5.9 \pm 4.8$ for $w = -1$, or $w = -0.91 \pm 0.12$ for $\Omega_k = 0$. When combined with the full CMB likelihood, the addition of the anisotropic clustering measurements to the spherically-averaged BAO location increases the constraining power on dark energy by a factor of 4 in a flat CDM cosmology with constant dark energy equation of state $w$ (giving $w = -0.87 \pm 0.05$). This impressive gain depends on our measurement of both the growth of structure and Alcock-Paczynski effect, and is not realised when marginalising over the amplitude of redshift space distortions. Combining with both the CMB and Supernovae Type Ia (SNeIa), we find $\Omega_{\rm m} = 0.281 \pm 0.014$ and $1000\Omega_{\rm k}=-9.2\pm5.0$ for $w = -1$, or $w_0 = -1.13 \pm 0.12$ and $w_{\rm a}=0.65 \pm 0.36$ assuming $\Omega_k = 0$. Finally, when a $\Lambda$CDM background expansion is assumed, the combination of our estimate of the growth rate with previous growth measurements provides tight constraints on the parameters describing possible deviations from GR giving $\gamma = 0.64 \pm 0.05$. For one parameter extensions of the flat $\Lambda$CDM model, we find a $\sim 2\sigma$ preference either for $w > -1$ or slower growth than in GR. However, the data is fully consistent with the concordance model, and the evidence for these additional parameters is weaker than $2\sigma$.

We analyze the density field of 264,283 galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) and included in the SDSS data release nine (DR9). In total, the SDSS DR9 BOSS data includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of more than 3,000 deg^2. We measure the power spectrum of these galaxies with redshifts 0.43 < z < 0.7 in order to constrain the amount of local non-Gaussianity, f_NL,local, in the primordial density field, paying particular attention to the impact of systematic uncertainties. The BOSS galaxy density field is systematically affected by the local stellar density and this influences the ability to accurately measure f_NL,local. In the absence of any correction, we find (erroneously) that the probability that f_NL,local is greater than zero, P(f_NL,local >0), is 99.5%. After quantifying and correcting for the systematic bias and including the added uncertainty, we find -45 < f_NL,local < 195 at 95% confidence, and P(f_NL,local >0) = 91.0%. A more conservative approach assumes that we have only learned the k-dependence of the systematic bias and allows any amplitude for the systematic correction; we find that the systematic effect is not fully degenerate with that of f_NL,local, and we determine that -82 < f_NL,local < 178 (at 95% confidence) and P(f_NL,local >0) = 68%. This analysis demonstrates the importance of accounting for the impact of Galactic foregrounds on f_NL,local measurements. We outline the methods that account for these systematic biases and uncertainties. We expect our methods to yield robust constraints on f_NL,local for both our own and future large-scale-structure investigations.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

We obtain constraints on the variation of the fundamental constants from the full shape of the redshift-space correlation function of a sample of luminous galaxies drawn from the Data Release 9 of the Baryonic Oscillations Spectroscopic Survey. We combine this information with data from recent CMB, BAO and H_0 measurements. We focus on possible variations of the fine structure constant \alpha and the electron mass m_e in the early universe, and study the degeneracies between these constants and other cosmological parameters, such as the dark energy equation of state parameter w_DE, the massive neutrinos fraction f_\nu, the effective number of relativistic species N_eff, and the primordial helium abundance Y_He. When only one of the fundamental constants is varied, our final bounds are \alpha / \alpha_0 = 0.9957_{-0.0042}^{+0.0041} and m_e /(m_e)_0 = 1.006_{-0.013}^{+0.014}. For their joint variation, our results are \alpha / \alpha_0 = 0.9901_{-0.0054}^{+0.0055} and m_e /(m_e)_0 = 1.028 +/- 0.019. Although when m_e is allowed to vary our constraints on w_DE are consistent with a cosmological constant, when \alpha is treated as a free parameter we find w_DE = -1.20 +/- 0.13; more than 1 \sigma away from its standard value. When f_\nu and \alpha are allowed to vary simultaneously, we find f_\nu < 0.043 (95% CL), implying a limit of \sum m_\nu < 0.46 eV (95% CL), while for m_e variation, we obtain f_nu < 0.086 (95% CL), which implies \sum m_\nu < 1.1 eV (95% CL). When N_eff or Y_He are considered as free parameters, their simultaneous variation with \alpha provides constraints close to their standard values (when the H_0 prior is not included in the analysis), while when m_e is allowed to vary, their preferred values are significantly higher. In all cases, our results are consistent with no variations of \alpha or m_e at the 1 or 2 \sigma level.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We explore the benefits of using a passively evolving population of galaxies to measure the evolution of the rate of structure growth between z=0.25 and z=0.65 by combining data from the SDSS-I/II and SDSS-III surveys. The large-scale linear bias of a population of dynamically passive galaxies, which we select from both surveys, is easily modeled. Knowing the bias evolution breaks degeneracies inherent to other methodologies, and decreases the uncertainty in measurements of the rate of structure growth and the normalization of the galaxy power-spectrum by up to a factor of two. If we translate our measurements into a constraint on sigma_8(z=0) assuming a concordance cosmological model and General Relativity (GR), we find that using a bias model improves our uncertainty by a factor of nearly 1.5. Our results are consistent with a flat Lambda Cold Dark Matter model and with GR.

We obtain constraints on cosmological parameters from the spherically averaged redshift-space correlation function of the CMASS Data Release 9 (DR9) sample of the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this information with additional data from recent CMB, SN and BAO measurements. Our results show no significant evidence of deviations from the standard flat-Lambda CDM model, whose basic parameters can be specified by Omega_m = 0.285 +- 0.009, 100 Omega_b = 4.59 +- 0.09, n_s = 0.96 +- 0.009, H_0 = 69.4 +- 0.8 km/s/Mpc and sigma_8 = 0.80 +- 0.02. The CMB+CMASS combination sets tight constraints on the curvature of the Universe, with Omega_k = -0.0043 +- 0.0049, and the tensor-to-scalar amplitude ratio, for which we find r < 0.16 at the 95 per cent confidence level (CL). These data show a clear signature of a deviation from scale-invariance also in the presence of tensor modes, with n_s <1 at the 99.7 per cent CL. We derive constraints on the fraction of massive neutrinos of f_nu < 0.049 (95 per cent CL), implying a limit of sum m_nu < 0.51 eV. We find no signature of a deviation from a cosmological constant from the combination of all datasets, with a constraint of w_DE = -1.033 +- 0.073 when this parameter is assumed time-independent, and no evidence of a departure from this value when it is allowed to evolve as w_DE(a) = w_0 + w_a (1 - a). The achieved accuracy on our cosmological constraints is a clear demonstration of the constraining power of current cosmological observations.

We analyze the density field of galaxies observed by the Sloan Digital Sky Survey (SDSS)-III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Nine (DR9). DR9 includes spectroscopic redshifts for over 400,000 galaxies spread over a footprint of 3,275 deg^2. We identify, characterize, and mitigate the impact of sources of systematic uncertainty on large-scale clustering measurements, both for angular moments of the redshift-space correlation function and the spherically averaged power spectrum, P(k), in order to ensure that robust cosmological constraints will be obtained from these data. A correlation between the projected density of stars and the higher redshift (0.43 < z < 0.7) galaxy sample (the `CMASS' sample) due to imaging systematics imparts a systematic error that is larger than the statistical error of the clustering measurements at scales s > 120h^-1Mpc or k < 0.01hMpc^-1. We find that these errors can be ameliorated by weighting galaxies based on their surface brightness and the local stellar density. We use mock galaxy catalogs that simulate the CMASS selection function to determine that randomly selecting galaxy redshifts in order to simulate the radial selection function of a random sample imparts the least systematic error on correlation function measurements and that this systematic error is negligible for the spherically averaged correlation function. The methods we recommend for the calculation of clustering measurements using the CMASS sample are adopted in companion papers that locate the position of the baryon acoustic oscillation feature (Anderson et al. 2012), constrain cosmological models using the full shape of the correlation function (Sanchez et al. 2012), and measure the rate of structure growth (Reid et al. 2012). (abridged)

We present a comprehensive study of 250,000 galaxies targeted by the Baryon Oscillation Spectroscopic Survey (BOSS) up to z ~ 0.7 with the specific goal of identifying and characterising a population of galaxies that has evolved without significant merging. We compute a likelihood that each BOSS galaxy is a progenitor of the Luminous Red Galaxies (LRGs) sample, targeted by SDSS-I/II up z ~ 0.5, by using the fossil record of LRGs and their inferred star-formation histories, metallicity histories and dust content. We determine merger rates, luminosity growth rates and the evolution of the large-scale clustering between the two surveys, and we investigate the effect of using different stellar population synthesis models in our conclusions. We demonstrate that our sample is slowly evolving (of the order of 2 +/- 1.5% per Gyr by merging). Our conclusions refer to the bright and massive end of the galaxy population, with Mi0.55 < -22, and M* > 1E11.2 Msolar, corresponding roughly to 95% and 40% of the LRGs and BOSS galaxy populations, respectively. Our analysis further shows that any possible excess of flux in BOSS galaxies, when compared to LRGs, from potentially unresolved targets at z ~ 0.55 must be less than 1% in the r0.55-band (approximately equivalent to the g-band in the rest-frame of galaxies at z=0.55). When weighting the BOSS galaxies based on the predicted properties of the LRGs, and restricting the analysis to the reddest BOSS galaxies, we find an evolution of the large-scale clustering that is consistent with dynamical passive evolution, assuming a standard cosmology. We conclude that our likelihoods give a weighted sample that is as clean and as close to passive evolution (in dynamical terms, i.e. no or negligible merging) as possible, and that is optimal for cosmological studies.

We present a fast method of producing mock galaxy catalogues that can be used to compute covariance matrices of large-scale clustering measurements and test the methods of analysis. Our method populates a 2nd-order Lagrangian Perturbation Theory (2LPT) matter field, where we calibrate masses of dark matter halos by detailed comparisons with N-body simulations. We demonstrate the clustering of halos is recovered at ~10 per cent accuracy. We populate halos with mock galaxies using a Halo Occupation Distribution (HOD) prescription, which has been calibrated to reproduce the clustering measurements on scales between 30 and 80 Mpc/h. We compare the sample covariance matrix from our mocks with analytic estimates, and discuss differences. We have used this method to make catalogues corresponding to Data Release 9 of the Baryon Oscillation Spectroscopic Survey (BOSS),producing 600 mock catalogues of the "CMASS" galaxy sample. These mocks enabled detailed tests of methods and errors that formed an integral part of companion analyses of these galaxy data.

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy 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 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

We analyze the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264,283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3,275 square degrees. Both peculiar velocities and errors in the assumed redshift-distance relation ("Alcock-Paczynski effect") generate correlations between clustering amplitude and orientation with respect to the line-of-sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broadband shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z=0.57, the Hubble expansion rate at z=0.57 (92.4 +/- 4.5 km/s/Mpc), and the growth rate of structure at that same redshift (d sigma8/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a LCDM expansion history, our growth constraint tightens to d sigma8/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of DA, H, and growth all separately require dark energy at z > 0.57, and when combined imply \Omega_{\Lambda} = 0.74 +/- 0.016, independent of the Universe's evolution at z<0.57. In our companion paper (Samushia et al. prep), we explore further cosmological implications of these observations.

We study the scale-dependence of halo bias in generic (non-local) primordial non-Gaussian (PNG) initial conditions of the type motivated by inflation, parametrized by an arbitrary quadratic kernel. We first show how to generate non-local PNG initial conditions with minimal overhead compared to local PNG models for a general class of primordial bispectra that can be written as linear combinations of separable templates. We run cosmological simulations for the local, and non-local equilateral and orthogonal models and present results on the scale-dependence of halo bias. We also derive a general formula for the Fourier-space bias using the peak-background split (PBS) in the context of the excursion set approach to halos and discuss the difference and similarities with the known corresponding result from local bias models. Our PBS bias formula generalizes previous results in the literature to include non-Markovian effects and non-universality of the mass function and are in better agreement with measurements in numerical simulations than previous results for a variety of halo masses, redshifts and halo definitions. We also derive for the first time quadratic bias results for arbitrary non-local PNG, and show that non-linear bias loops give small corrections at large-scales. The resulting well-behaved perturbation theory paves the way to constrain non-local PNG from measurements of the power spectrum and bispectrum in galaxy redshift surveys.

We explore the biasing in the clustering statistics of halos as compared to dark matter (DM) in simulations. We look at the second and third order statistics at large scales of the (intermediate) MICEL1536 simulation and also measure directly the local bias relation h = f({\delta}) between DM fluctuations, {\delta}, smoothed over a top-hat radius Rs at a point in the simulation and its corresponding tracer h (i.e. halos) at the same point. This local relation can be Taylor expanded to define a linear (b1) and non-linear (b2) bias parameters. The values of b1 and b2 in the simulation vary with Rs approaching a constant value around Rs > 30 - 60 Mpc/h. We use the local relation to predict the clustering of the tracer in terms of the one of DM. This prediction works very well (about percent level) for the halo 2-point correlation {\xi}(r_12) for r_12 > 15 Mpc/h, but only when we use the biasing values that we found at very large smoothing radius Rs > 30 - 60 Mpc/h. We find no effect from stochastic or next to leading order terms in the f({\delta}) expansion. But we do find some discrepancies in the 3-point function that needs further understanding. We also look at the clustering of the smoothed moments, the variance and skewness which are volume average correlations and therefore include clustering from smaller scales. In this case, we find that both next to leading order and discreetness corrections (to the local model) are needed at the 10 - 20% level. Shot-noise can be corrected with a term {\sigma}e^2/n where {\sigma}e^2 < 1, i.e., always smaller than the Poisson correction. We also compare these results with the peak-background split predictions from the measured halo mass function. We find 5-10% systematic (and similar statistical) errors in the mass estimation when we use the halo model biasing predictions to calibrate the mass.

The peak-background split argument is commonly used to relate the abundance of dark matter halos to their spatial clustering. Testing this argument requires an accurate determination of the halo mass function. We present a Maximum Likelihood method for fitting parametric functional forms to halo abundances which differs from previous work because it does not require binned counts. Our conclusions do not depend on whether we use our method or more conventional ones. In addition, halo abundances depend on how halos are defined. Our conclusions do not depend on the choice of link length associated with the friends-of-friends halo-finder, nor do they change if we identify halos using a spherical overdensity algorithm instead. The large scale halo bias measured from the matter-halo cross spectrum b_x and the halo autocorrelation function b_xi (on scales k~0.03h/Mpc and r ~50 Mpc/h) can differ by as much as 5% for halos that are significantly more massive than the characteristic mass M*. At these large masses, the peak background split estimate of the linear bias factor b1 is 3-5% smaller than b_xi, which is 5% smaller than b_x. We discuss the origin of these discrepancies: deterministic nonlinear local bias, with parameters determined by the peak-background split argument, is unable to account for the discrepancies we see. A simple linear but nonlocal bias model, motivated by peaks theory, may also be difficult to reconcile with our measurements. More work on such nonlocal bias models may be needed to understand the nature of halo bias at this level of precision.

Galaxy surveys provide a large-scale view of the universe that typically has a limited line-of-sight or redshift resolution. The lack of radial accuracy in these surveys can be modelled by picturing the universe as a set of concentric radial shells of finite width around the observer, i.e, an onion-like structure. We present a new N-body simulation with 2048^3 particles developed at the Marenostrum supercomputer with the GADGET-2 code. Using the lightcone output we build a set of angular maps that mimic this onion-like representation of the universe. The onion maps are a highly compressed version of the raw data (i.e. a factor >1000 smaller size for arcminute resolution maps) and they provide a new and powerful tool to exploit large scale structure observations. We introduce two basic applications of these maps that are especially useful for constraning dark energy properties: the baryon acoustic oscillations (BAO) in the galaxy power spectrum and all-sky maps of the weak lensing distortion. Using the weak lensing maps, we measure the convergence power spectra and compare it to halo fit predictions. As a further application, we compute the variance and higher-order moments of the maps. We show that sampling variance on scales of few degrees is quite large, resulting in a significant (25% at 10 arminute scales) bias in the variance. We caution that current lensing surveys such as the COSMOS HST should take into account this bias and extra sampling error in their clustering analyses and inferred cosmological parameter constraints. Finally, we test the importance of projection effects in the weak lensing mass reconstruction. On the mean, the mass calibration works well but it exhibits a large non-Gaussian scatter what could induce a large bias in the recovered mass function.

Constraining cosmological parameters from measurements of the Integrated Sachs-Wolfe effect requires developing robust and accurate methods for computing statistical errors in the cross-correlation between maps. This paper presents a detailed comparison of such error estimation applied to the case of cross-correlation of Cosmic Microwave Background (CMB) and large-scale structure data. We compare theoretical models for error estimation with montecarlo simulations where both the galaxy and the CMB maps vary around a fiducial auto-correlation and cross-correlation model which agrees well with the current concordance LCDM cosmology. Our analysis compares estimators both in harmonic and configuration (or real) space, quantifies the accuracy of the error analysis and discuss the impact of partial sky survey area and the choice of input fiducial model on dark-energy constraints. We show that purely analytic approaches yield accurate errors even in surveys that cover only 10% of the sky and that parameter constraints strongly depend on the fiducial model employed. Alternatively, we discuss the advantages and limitations of error estimators that can be directly applied to data. In particular, we show that errors and covariances from the Jack-Knife method agree well with the theoretical approaches and simulations. We also introduce a novel method in real space that is computationally efficient and can be applied to real data and realistic survey geometries. Finally, we present a number of new findings and prescriptions that can be useful for analysis of real data and forecasts, and present a critical summary of the analyses done to date.

We cross-correlate the third-year WMAP data with galaxy samples extracted from the SDSS DR4 (SDSS4) covering 13% of the sky, increasing by a factor of 3.7 the volume sampled in previous analyses. The new measurements confirm a positive cross-correlation with higher significance (total signal-to-noise of about 4.7). The correlation as a function of angular scale is well fitted by the integrated Sachs-Wolfe (ISW) effect for LCDM flat FRW models with a cosmological constant. The combined analysis of different samples gives Omega_L=0.80-0.85$ (68% Confidence Level, CL) or $0.77-0.86$ (95% CL). We find similar best fit values for Omega_L for different galaxy samples with median redshifts of z ~0.3 and z ~0.5, indicating that the data scale with redshift as predicted by the LCDM cosmology (with equation of state parameter w=-1). This agreement is not trivial, but can not yet be used to break the degeneracy constraints in the w versus Omega_L plane using only the ISW data.

Cluster number counts can be used to test dark energy models. We investigate dark energy candidates which are coupled to dark matter. We analyze the cluster number counts dependence on the amount of dark matter coupled to dark energy. Further more, we study how dark energy inhomogeneities affect cluster abundances. It is shown that increasing the coupling reduces significantly the cluster number counts, and that dark energy inhomogeneities increases cluster abundances. Wiggles in cluster number counts are shown to be a specific signature of coupled dark energy models. Future observations will possibly detect such oscillations and discriminate among the different dark energy models.

Recent studies by a number of independent collaborations, have correlated the CMB temperatures measured by the WMAP satellite with different galaxy surveys that trace the matter distribution with light from the whole range of the electromagnetic spectrum: radio, far-infrared, optical and X-ray surveys. The new data systematically finds positive correlations, indicating a rapid slow down in the growth of structure in the universe. Individual cross-correlation measurements are of low significance, but we show that combining data at different redshifts introduces important new constraints. Contrary to what happens at low redshifts, for a fixed $\Omm$, the higher the dark energy contend, $\Ol$, the lower the ISW cross-correlation amplitude. At 68% confidence level, the data finds new independent evidence of dark energy: $\Ol =0.42-1.22$ . It also confirms, to higher significance, the presence of a large dark matter component: $\Omm =0.18-0.34$, exceeding the density of baryonic matter, but far from the critical value. Combining these new constraints with the prior of a flat universe, or the prior of an accelerating universe provides strong new evidence for a dark cosmos. Combination with supernova data yields $\Ol = 0.71 \pm 0.13$, $\Omm = 0.29 \pm 0.04$. If we also assume a flat universe, we find $\Ol = 0.70 \pm 0.05$ and $w = -1.02 \pm 0.17$ for a constant dark energy equation of state.

The growth of large scale structure is studied in a universe containing both cold dark matter (CDM) and generalized Chaplygin gas (GCg). GCg is assumed to contribute only to the background evolution of the universe while the CDM component collapses and forms structures. We present some new analytical as well as numerical results for linear and non-linear growth in such model. The model passes the standard cosmological distance test without the need of a cosmological constant (LCDM). But we find that the scenario is severely constrained by current observations of large scale structure. Any small deviations of the GCg parameters away from the standard Lambda dominated cosmology (LCDM) produces substantial suppression for the growth of structures.

We study the growth of large scale structure in two recently proposed non-standard cosmological models: the brane induced gravity model of Dvali, Gabadadze and Porrati (DGP) and the Cardassian models of Freese and Lewis. A general formalism for calculating the growth of fluctuations in models with a non-standard Friedman equation and a normal continuity equation of energy density is discussed. Both linear and non-linear growth are studied, together with their observational signatures on higher order statistics and abundance of collapsed objects. In general, models which show similar cosmic acceleration at z ~ 1, can produce quite different normalization for large scale density fluctuations, ie sigma_8, cluster abundance or higher order statistics, such as the normalized skewness S_3, which is independent of the linear normalization. For example, for a flat universe with Omega_M ~ 0.22, DGP and standard Cardassian cosmologies predict about 2 and 3 times more clusters respectively than the standard Lambda model at z=1.5. When normalized to CMB fluctuations the sigma_8 amplitude turns out to be lower by a few tens of a percent. We also find that, for a limited red-shift range, the linear growth rate can be faster in some models (eg modified polytropic Cardassian with q>1) than in the Einstein-deSitter universe. The value of the skewness S_3 is found to have up to ~10 percent variations (up or down) from model to model.