Improved renormalization group computation of likelihood functions for cosmological data sets

Abstract

Evaluation of likelihood functions for cosmological large scale structure data sets (including CMB, galaxy redshift surveys, etc.) naturally involves marginalization, i.e., integration, over an unknown underlying random signal field. Recently, I showed how a renormalization group method can be used to carry out this integration efficiently by first integrating out the smallest scale structure, i.e., localized structure on the scale of differences between nearby data cells, then combining adjacent cells in a coarse graining step, then repeating this process over and over until all scales have been integrated. Here I extend the formulation in several ways in order to reduce the prefactor on the method's linear scaling with data set size. The key improvement is showing how to integrate out the difference between specific adjacent cells before summing them in the coarse graining step, compared to the original formulation in which small-scale fluctuations were integrated more generally. I suggest some other improvements in details of the scheme, including showing how to perform the integration around a maximum likelihood estimate for the underlying random field. In the end, an accurate likelihood computation for a million-cell Gaussian test data set runs in two minutes on my laptop, with room for further optimization and straightforward parallelization.

10 pg, 2 fig

References in corpus (23)

• Clustering of dark matter tracers: generalizing bias for the coming era of precision LSS
• Detection of Gravitational Lensing in the Cosmic Microwave Background
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in Fourier-space
• The Scale Dependence of Halo and Galaxy Bias: Effects in Real Space
• Understanding higher-order nonlocal halo bias at large scales by combining the power spectrum with the bispectrum
• Primordial non-Gaussianity: large-scale structure signature in the perturbative bias model
• Nonlinear Perturbation Theory Integrated with Nonlocal Bias, Redshift-space Distortions, and Primordial Non-Gaussianity
• Dark matter clustering: a simple renormalization group approach
• Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
• Perturbation Theory Reloaded II: Non-linear Bias, Baryon Acoustic Oscillations and Millennium Simulation In Real Space
• Large-N expansions applied to gravitational clustering
• Baryonic Acoustic Oscillations via the Renormalization Group
• Scale-Dependent Galaxy Bias from Massive Particles with Spin during Inflation
• GridSPT: Grid-based calculation for perturbation theory of large-scale structure
• Solving the Vlasov equation in two spatial dimensions with the Schrödinger method
• Information theory for fields
• How to estimate the 3D power spectrum of the Lyman-$\alpha$ forest
• The Generation of Vorticity in Cosmological N-body Simulations
• Renormalized Newtonian Cosmic Evolution with Primordial Non-Gaussianity
• Large-scale structure perturbation theory without losing stream crossing
• Halo Sampling, Local Bias and Loop Corrections
• Lensing smoothing of BAO wiggles
• Renormalization group computation of likelihood functions for cosmological data sets

Cited by in corpus (2)

• Testing Quadratic Maximum Likelihood estimators for forthcoming Stage-IV weak lensing surveys
• AKRA 2.0: Accurate Kappa Reconstruction Algorithm for masked shear catalog