English

Renormalization group computation of likelihood functions for cosmological data sets

Cosmology and Nongalactic Astrophysics 2019-03-06 v1

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 N\sim N time, even for problems with no symmetry, compared to N3N^3 for brute force linear algebra (where NN is the number of data points -- to be fair, methods already exist to solve the Gaussian problem in at worst NlogNN \log N 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.

Keywords

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

R2 v1 2026-06-23T04:45:42.765Z