Exact joint likelihood of pseudo-Cℓ estimates from correlated Gaussian cosmological fields

We present the exact joint likelihood of pseudo-C_ℓ power spectrum estimates measured from an arbitrary number of Gaussian cosmological fields. Our method is applicable to both spin-0 fields and spin-2 fields, including a mixture of the two, and is relevant to cosmic microwave background (CMB), weak lensing, and galaxy clustering analyses. We show that Gaussian cosmological fields are mixed by a mask in such a way that retains their Gaussianity and derive exact expressions for the covariance of the cut-sky spherical harmonic coefficients, the pseudo-a_ℓms, without making any assumptions about the mask geometry. We then show that each auto or cross-pseudo-C_ℓ estimator can be written as a quadratic form, and apply the known joint distribution of quadratic forms to obtain the exact joint likelihood of a set of pseudo-C_ℓ estimates in the presence of an arbitrary mask. We show that the same formalism can be applied to obtain the exact joint likelihood of quadratic maximum likelihood power spectrum estimates. Considering the polarization of the CMB as an example, we show using simulations that our likelihood recovers the full, exact multivariate distribution of EE, BB, and EB pseudo-C_ℓ power spectra. Our method provides a route to robust cosmological constraints from future CMB and large-scale structure surveys in an era of ever-increasing statistical precision.