Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNear optimal bispectrum estimators for large-scale structure
174
0
0.0
(
0
)
Added by
Marcel M. Schmittfull
Publication date
2014
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Clustering of large-scale structure provides significant cosmological information through the power spectrum of density perturbations. Additional information can be gained from higher-order statistics like the bispectrum, especially to break the degeneracy between the linear halo bias $b_1$ and the amplitude of fluctuations $sigma_8$. We propose new simple, computationally inexpensive bispectrum statistics that are near optimal for the specific applications like bias determination. Corresponding to the Legendre decomposition of nonlinear halo bias and gravitational coupling at second order, these statistics are given by the cross-spectra of the density with three quadratic fields: the squared density, a tidal term, and a shift term. For halos and galaxies the first two have associated nonlinear bias terms $b_2$ and $b_{s^2}$, respectively, while the shift term has none in the absence of velocity bias (valid in the $k rightarrow 0$ limit). Thus the linear bias $b_1$ is best determined by the shift cross-spectrum, while the squared density and tidal cross-spectra mostly tighten constraints on $b_2$ and $b_{s^2}$ once $b_1$ is known. Since the form of the cross-spectra is derived from optimal maximum-likelihood estimation, they contain the full bispectrum information on bias parameters. Perturbative analytical predictions for their expectation values and covariances agree with simulations on large scales, $klesssim 0.09h/mathrm{Mpc}$ at $z=0.55$ with Gaussian $R=20h^{-1}mathrm{Mpc}$ smoothing, for matter-matter-matter, and matter-matter-halo combinations. For halo-halo-halo cross-spectra the model also needs to include corrections to the Poisson stochasticity.
rate research
Read More
We study perturbation theory for large-scale structure in the most general scalar-tensor theories propagating a single scalar degree of freedom, which include Horndeski theories and beyond. We model the parameter space using the effective field theory of dark energy. For Horndeski theories, the gravitational field and fluid equations are invariant under a combination of time-dependent transformations of the coordinates and fields. This symmetry allows one to construct a physical adiabatic mode which fixes the perturbation-theory kernels in the squeezed limit and ensures that the well-known consistency relations for large-scale structure, originally derived in general relativity, hold in modified gravity as well. For theories beyond Horndeski, instead, one generally cannot construct such an adiabatic mode. Because of this, the perturbation-theory kernels are modified in the squeezed limit and the consistency relations for large-scale structure do not hold. We show, however, that the modification of the squeezed limit depends only on the linear theory. We investigate the observational consequences of this violation by computing the matter bispectrum. In the squeezed limit, the largest effect is expected when considering the cross-correlation between different tracers. Moreover, the individual contributions to the 1-loop matter power spectrum do not cancel in the infrared limit of the momentum integral, modifying the power spectrum on non-linear scales.
When analyzing the galaxy bispectrum measured from spectroscopic surveys, it is imperative to account for the effects of non-uniform survey geometry. Conventionally, this is done by convolving the theory model with the the window function; however, the computational expense of this prohibits full exploration of the bispectrum likelihood. In this work, we provide a new class of estimators for the unwindowed bispectrum; a quantity that can be straightforwardly compared to theory. This builds upon the work of Philcox (2021) for the power spectrum, and comprises two parts (both obtained from an Edgeworth expansion): a cubic estimator applied to the data, and a Fisher matrix, which deconvolves the bispectrum components. In the limit of weak non-Gaussianity, the estimator is minimum-variance; furthermore, we give an alternate form based on FKP weights that is close-to-optimal and easy to compute. As a demonstration, we measure the binned bispectrum monopole of a suite of simulations both using conventional estimators and our unwindowed equivalents. Computation times are comparable, except that the unwindowed approach requires a Fisher matrix, computable in an additional $mathcal{O}(100)$ CPU-hours. Our estimator may be straightforwardly extended to measure redshift-space distortions and the components of the bispectrum in arbitrary separable bases. The techniques of this work will allow the bispectrum to straightforwardly included in the cosmological analysis of current and upcoming survey data.
We study the matter bispectrum of the large-scale structure by comparing different perturbative and phenomenological models with measurements from $N$-body simulations obtained with a modal bispectrum estimator. Using shape and amplitude correlators, we directly compare simulated data with theoretical models over the full three-dimensional domain of the bispectrum, for different redshifts and scales. We review and investigate the main perturbative methods in the literature that predict the one-loop bispectrum: standard perturbation theory, effective field theory, resummed Lagrangian and renormalised perturbation theory, calculating the latter also at two loops for some triangle configurations. We find that effective field theory (EFT) succeeds in extending the range of validity furthest into the mildly nonlinear regime, albeit at the price of free extra parameters requiring calibration on simulations. For the more phenomenological halo model, we confirm that despite its validity in the deeply nonlinear regime it has a deficit of power on intermediate scales, which worsens at higher redshifts; this issue is ameliorated, but not solved, by combined halo-perturbative models. We show from simulations that in this transition region there is a strong squeezed bispectrum component that is significantly underestimated in the halo model at earlier redshifts. We thus propose a phenomenological method for alleviating this deficit, which we develop into a simple phenomenological three-shape benchmark model based on the three fundamental shapes we have obtained from studying the halo model. When calibrated on the simulations, this three-shape benchmark model accurately describes the bispectrum on all scales and redshifts considered, providing a prototype bispectrum HALOFIT-like methodology that could be used to describe and test parameter dependencies.
Higher-order clustering statistics, like the galaxy bispectrum, can add complementary cosmological information to what is accessible with two-point statistics, like the power spectrum. While the standard way of measuring the bispectrum involves estimating a bispectrum value in a large number of Fourier triangle bins, the compressed modal bispectrum approximates the bispectrum as a linear combination of basis functions and estimates the expansion coefficients on the chosen basis. In this work, we compare the two estimators by using parallel pipelines to analyze the real-space halo bispectrum measured in a suite of $N$-body simulations corresponding to a total volume of $sim 1{,}000 ,h^{-3},{rm Gpc}^3$, with covariance matrices estimated from 10,000 mock halo catalogs. We find that the modal bispectrum yields constraints that are consistent and competitive with the standard bispectrum analysis: for the halo bias and shot noise parameters within the tree-level halo bispectrum model up to $k_{rm max} approx 0.06 , (0.10) ,h,{rm Mpc}^{-1}$, only 6 (10) modal expansion coefficients are necessary to obtain constraints equivalent to the standard bispectrum estimator using $sim$ 20 to 1,600 triangle bins, depending on the bin width. For this work, we have implemented a modal estimator pipeline using Markov Chain Monte Carlo simulations for the first time, and we discuss in detail how the parameter posteriors and modal expansion are robust to, or sensitive to, several user settings within the modal bispectrum pipeline. The combination of the highly efficient compression that is achieved and the large number of mock catalogs available allows us to quantify how our modal bispectrum constraints depend on the number of mocks that are used to estimate covariance matrices and the functional form of the likelihood.
The jackknife method gives an internal covariance estimate for large-scale structure surveys and allows model-independent errors on cosmological parameters. Using the SDSS-III BOSS CMASS sample, we study how the jackknife size and number of resamplings impact the precision of the covariance estimate on the correlation function multipoles and the error on the inferred baryon acoustic scale. We compare the measurement with the MultiDark Patchy mock galaxy catalogues, and we also validate it against a set of log-normal mocks with the same survey geometry. We build several jackknife configurations that vary in size and number of resamplings. We introduce the Hartlap factor in the covariance estimate that depends on the number of jackknife resamplings. We also find that it is useful to apply the tapering scheme to estimate the precision matrix from a limited number of resamplings. The results from CMASS and mock catalogues show that the error estimate of the baryon acoustic scale does not depend on the jackknife scale. For the shift parameter $alpha$, we find an average error of 1.6%, 2.2% and 1.2%, respectively from CMASS, Patchy and log-normal jackknife covariances. Despite these uncertainties fluctuate significantly due to some structural limitations of the jackknife method, our $alpha$ estimates are in reasonable agreement with published pre-reconstruction analyses. Jackknife methods will provide valuable and complementary covariance estimates for future large-scale structure surveys.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
