Theoretical modeling of the redshift-space power spectrum of galaxies is crucially important to correctly extract cosmological information from galaxy redshift surveys. The task is complicated by the nonlinear biasing and redshift space distortion (RSD) effects, which change with halo mass, and by the wide distribution of halo masses and their occupations by galaxies. One of the main modeling challenges is the existence of satellite galaxies that have both radial distribution inside the halos and large virial velocities inside halos, a phenomenon known as the Finger-of-God (FoG) effect. We present a model for the redshift-space power spectrum of galaxies in which we decompose a given galaxy sample into central and satellite galaxies and relate different contributions to the power spectrum to 1-halo and 2-halo terms in a halo model. Our primary goal is to ensure that any parameters that we introduce have physically meaningful values, and are not just fitting parameters. For the lowest order 2-halo terms we use the previously developed RSD modeling of halos in the context of distribution function and perturbation theory approach. This term needs to be multiplied by the effect of radial distances and velocities of satellites inside the halo. To this one needs to add the 1-halo terms, which are nonperturbative. We show that the real space 1-halo terms can be modeled as almost constant, with the finite extent of the satellites inside the halo inducing a small k2R2 term over the range of scales of interest, where R is related to the size of the halo given by its halo mass. We adopt a similar model for FoG in redshift space, ensuring that FoG velocity dispersion is related to the halo mass. For FoG k2 type expansions do not work over the range of scales of interest and FoG resummation must be used instead. We test several simple damping functions to model the velocity dispersion FoG effect. Applying the formalism to mock galaxies modeled after the "CMASS" sample of the BOSS survey, we find that our predictions for the redshift-space power spectra are accurate up to k≃0.4 h Mpc-1 within 1% if the halo power spectrum is measured using N-body simulations and within 3% if it is modeled using perturbation theory.

The Dark Energy Spectroscopic Instrument (DESI), a multiplexed fiber-fed spectrograph, is a Stage-IV ground-based dark energy experiment aiming to measure redshifts for 29 million Emission-Line Galaxies (ELG), 4 million Luminous Red Galaxies (LRG), and 2 million Quasi-Stellar Objects (QSO). The survey design includes a pattern of tiling on the sky, the locations of the fiber positioners in the focal plane of the telescope, and an observation strategy determined by a fiber assignment algorithm that optimizes the allocation of fibers to targets. This strategy allows a given region to be covered on average five times for a five-year survey, with a typical variation of about 1.5 about the mean, which imprints a spatially-dependent pattern on the galaxy clustering. We investigate the systematic effects of the fiber assignment coverage on the anisotropic galaxy clustering of ELGs and show that, in the absence of any corrections, it leads to discrepancies of order ten percent on large scales for the power spectrum multipoles. We introduce a method where objects in a random catalog are assigned a coverage, and the mean density is separately computed for each coverage factor. We show that this method reduces, but does not eliminate the effect. We next investigate the angular dependence of the contaminated signal, arguing that it is mostly localized to purely transverse modes. We demonstrate that the cleanest way to remove the contaminating signal is to perform an analysis of the anisotropic power spectrum P(k,μ) and remove the lowest μ bin, leaving μ > 0 modes accurate at the few-percent level. Here, μ is the cosine of the angle between the line-of-sight and the direction of K→. We also investigate two alternative definitions of the random catalog and show that they are comparable but less effective than the coverage randoms method.

We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface, we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, two- and three-point correlation functions, a friends-of-friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enables nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high-performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at http://nbodykit.readthedocs.io, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys.

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

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147 000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts 0.8 < z < 2.2 and measure their spherically averaged clustering in both configuration and Fourier space. Our observational data set and the 1400 simulated realizations of the data set allow us to detect a preference for BAO that is greater than 2.8σ. We determine the spherically averaged BAO distance to z = 1.52 to 3.8 per cent precision: DV (z = 1.52) = 3843 ± 147 (rd/rd,fid) Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flatΛCDMbest-fitting cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find ΩΛ > 0 at 6.6s significance when testing a ΛCDM model with free curvature.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg2 and volume of 18.7 Gpc3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51 and 0.61. We measure the angular diameter distance DM and Hubble parameterHfrom the baryon acoustic oscillation (BAO) method, in combinationwith a cosmic microwave background prior on the sound horizon scale, after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DMH from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by fσ8(z), from redshift-space distortions (RSD). We combine individual measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method; in particular, the AP measurement from sub-BAO scales sharpens constraints from post-reconstruction BAOs by breaking degeneracy between DM and H. Combined with Planck 2016 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature ωK = 0.0003 ± 0.0026 and a dark energy equation-of-state parameter ω =-1.01 ± 0.06, in strong affirmation of the spatially flat cold dark matter (CDM) model with a cosmological constant (ΛCDM). Our RSD measurements of fσ8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3 ± 1.0 kms-1 Mpc-1 even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8 ± 1.2 kms-1 Mpc-1. Assuming flat ΛCDM, we find ωm = 0.310 ± 0.005 and H0 = 67.6 ± 0.5 kms-1 Mpc-1, and we find a 95 per cent upper limit of 0.16 eV c-2 on the neutrino mass sum.