Dr. Tassia Ferreira

Abstract:

Covariance matrices are important tools for obtaining reliable parameter constraints. Advancements in cosmological surveys lead to larger data vectors and, consequently, increasingly complex covariance matrices, whose number of elements grows as the square of the size of the data vector. The most straightforward way of comparing these matrices, in terms of their ability to produce parameter constraints, involves a full cosmological analysis, which can be very computationally expensive. Using the concept and construction of compression schemes, which have become increasingly popular, we propose a fast and reliable way of comparing covariance matrices. The basic idea is to focus only on the portion of the covariance matrix that is relevant for the parameter constraints and quantify, via a fast Monte Carlo simulation, the difference of a second candidate matrix from the baseline one. To test this method, we apply it to two covariance matrices that were used to analyse the cosmic shear measurements for the Dark Energy Survey Year 1. We found that the uncertainties on the parameters change by 2.6 per cent, a figure in agreement with the full cosmological analysis. While our approximate method cannot replace a full analysis, it may be useful during the development and validation of codes that estimate covariance matrices. Our method takes roughly 100 times less CPUh than a full cosmological analysis.

Abstract:

Covariance matrices are among the most difficult pieces of end-to-end cosmological analyses. In principle, for two-point functions, each component involves a four-point function, and the resulting covariance often has hundreds of thousands of elements. We investigate various compression mechanisms capable of vastly reducing the size of the covariance matrix in the context of cosmic shear statistics. This helps identify which of its parts are most crucial to parameter estimation. We start with simple compression methods, by isolating and “removing” 200 modes associated with the lowest eigenvalues, then those with the lowest signal-to-noise ratio, before moving on to more sophisticated schemes like compression at the tomographic level and, finally, with the massively optimized parameter estimation and data compression (MOPED). We find that, while most of these approaches prove useful for a few parameters of interest, like Ωm, the simplest yield a loss of constraining power on the intrinsic alignment (IA) parameters as well as S8. For the case considered—cosmic shear from the first year of data from the Dark Energy Survey—only MOPED was able to replicate the original constraints in the 16-parameter space. Finally, we apply a tolerance test to the elements of the compressed covariance matrix obtained with MOPED and confirm that the IA parameter AIA is the most susceptible to inaccuracies in the covariance matrix.