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.
ABSTRACT
In this paper, we predict the covariance matrices of both the power spectrum and the bispectrum, including full non-Gaussian contributions, redshift space distortions, linear bias effects, ...and shot-noise corrections, using perturbation theory (PT). To quantify the redshift-space distortion effect, we focus mainly on the monopole and quadrupole components of both the power and bispectra. We, for the first time, compute the 5- and 6-point spectra to predict the cross-covariance between the power and bispectra, and the autocovariance of the bispectrum in redshift space. We test the validity of our calculations by comparing them with the covariance matrices measured from the MultiDark-Patchy mock catalogues that are designed to reproduce the galaxy clustering measured from the Baryon Oscillation Spectroscopic Survey Data Release 12. We argue that the simple, leading-order PT works because the shot-noise corrections for the Patchy mocks are more dominant than other higher order terms we ignore. In the meantime, we confirm some discrepancies in the comparison, especially of the cross-covariance. We discuss potential sources of such discrepancies. We also show that our PT model reproduces well the cumulative signal-to-noise ratio of the power spectrum and the bispectrum as a function of maximum wavenumber, implying that our PT model captures successfully essential contributions to the covariance matrices.
Low-redshift measurements of baryon acoustic oscillations (BAO) measure the late-time evolution of the Universe and are a vital probe of dark energy. Over the past decade both the 6-degree Field ...Galaxy Survey (6dFGS) and Sloan Digital Sky Survey Main Galaxy Sample (SDSS MGS) have provided important distance constraints at Formula: see text < 0.3. In this paper we re-evaluate the cosmological information from the BAO detection in 6dFGS making use of halo occupation distribution (HOD)-populated COmoving Lagrangian Acceleration (COLA) mocks for an accurate covariance matrix and take advantage of the now commonly implemented technique of density field reconstruction. For the 6dFGS data, we find consistency with the previous analysis, and obtain an isotropic volume-averaged distance measurement of Formula: see text but with a highly non-Gaussian likelihood. We combine our measurement from both the post-reconstruction clustering of 6dFGS and SDSS MGS offering an updated constraint in this redshift regime, Formula: see text. This measurement tightens the constraint in comparison to the result from SDSS MGS alone, especially at the 2σ and higher significance levels. These measurements are consistent with standard Λcold dark matter (ΛCDM) and after fixing the standard ruler using a Planck prior on Ω
, the joint analysis gives Formula: see text. This result is consistent with other BAO and Cosmic microwave background (CMB) studies but is in >2σ tension with supernova distance ladder measurements. In the near future both the Taipan Galaxy Survey and the Dark Energy Spectroscopic Instrument (DESI) will improve this measurement to Formula: see text at low redshift.
Modeling the reconstructed BAO in Fourier space Seo, Hee-Jong; Beutler, Florian; Ross, Ashley J ...
Monthly notices of the Royal Astronomical Society,
08/2016, Letnik:
460, Številka:
3
Journal Article
Recenzirano
Odprti dostop
The density field reconstruction technique, which partially reverses the non-linear degradation of the Baryon acoustic oscillation (BAO) feature in the galaxy redshift surveys, has been successful in ...substantially improving the cosmology constraints from recent surveys such as Baryon Oscillation Spectroscopic Survey (BOSS). We estimate the efficiency of the method as a function of various reconstruction details. To directly quantify the BAO information in non-linear density fields before and after reconstruction, we calculate the cross-correlations (i.e. propagators) of the pre(post)-reconstructed density field with the initial linear field using a mock sample that mimics the clustering of the BOSS galaxies. The results directly provide the BAO damping as a function of wavenumber that can be implemented into the Fisher matrix analysis. We focus on investigating the dependence of the propagator on a choice of smoothing filters and on two major different conventions of the redshift-space density field reconstruction that have been used in literature. By estimating the BAO signal to noise for each case, we predict constraints on the angular diameter distance and Hubble parameter using the Fisher matrix analysis. We thus determine an optimal Gaussian smoothing filter scale for the signal-to-noise level of the BOSS CMASS. We also present appropriate BAO fitting models for different reconstruction methods based on the first- and second-order Lagrangian perturbation theory in Fourier space. Using the mock data, we show that the modified BAO fitting model can substantially improve the accuracy of the BAO position in the best fits as well as the goodness of the fits.