Renormalization group computation of likelihood functions for cosmological data sets
Abstract
I show how a renormalization group (RG) method can be used to incrementally integrate the information in cosmological large-scale structure data sets (including CMB, galaxy redshift surveys, etc.). I show numerical tests for Gaussian fields, where the method allows arbitrarily close to exact computation of the likelihood function in order time, even for problems with no symmetry, compared to for brute force linear algebra (where is the number of data points -- to be fair, methods already exist to solve the Gaussian problem in at worst time, and this method will not necessarily be faster in practice). The method requires no sampling or other Monte Carlo (random) element. Non-linearity/non-Gaussianity can be accounted for to the extent that terms generated by integrating out small scale modes can be projected onto a sufficient basis, e.g., at least in the sufficiently perturbative regime. The formulas to evaluate are straightforward and require no understanding of quantum field theory, but this paper may also serve as a pedagogical introduction to Wilsonian RG for astronomers.
Cite
@article{arxiv.1810.08454,
title = {Renormalization group computation of likelihood functions for cosmological data sets},
author = {Patrick McDonald},
journal= {arXiv preprint arXiv:1810.08454},
year = {2019}
}
Comments
13 pg, 4 fig