Analytic Model for Covariance Matrices of the 2-, 3-, and 4-Point Correlation Functions in the Gaussian Random Field Approximation
Abstract
Analyses of the galaxy N-Point Correlation Functions (NPCFs) have a large number of degrees of freedom, meaning one cannot directly estimate an invertible covariance matrix purely from mock catalogs, as has been the standard approach for the 2PCF and power spectrum. Instead, analyses use templates based on assuming a Gaussian Random Field density with the true, Boltzmann-solver-computed power spectrum. The resulting covariance matrices are sparse but have notable internal structure. To understand this structure better, we seek a fully analytic, closed-form covariance matrix template, using a power law power spectrum $P(k) \propto 1/k$. We obtain a simple closed-form solution for the covariance of the 2PCF, as well as closed-form solutions for the fundamental building blocks of the covariance matrices for the 3PCF, 4PCF, and beyond. We use our results to present a clearer picture of the covariance matrices' structure and sparsity, corresponding to triangular and non-triangular regions. This will be useful in guiding future NPCF analyses with spectroscopic surveys such as DESI, Euclid, Roman, and SPHEREx.