Data compression and covariance matrix inspection: cosmic shear

Physical Review D American Physical Society 103:10 (2021) 103535

Authors:

Tassia Ferreira, Tianqing Zhang, Nianyi Chen, Scott Dodelson

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.

The SAMI Galaxy Survey: a statistical approach to an optimal classification of stellar kinematics in galaxy surveys

Monthly Notices of the Royal Astronomical Society Oxford University Press 505:2 (2021) 3078-3106

Authors:

Jesse van de Sande, Sam P Vaughan, Luca Cortese, Nicholas Scott, Joss Bland-Hawthorn, Scott M Croom, Claudia DP Lagos, Sarah Brough, Julia J Bryant, Julien Devriendt, Yohan Dubois, Francesco D'Eugenio, Caroline Foster, Amelia Fraser-McKelvie, Katherine E Harborne, Jon S Lawrence, Sree Oh, Matt S Owers, Adriano Poci, Rhea-Silvia Remus, Samuel N Richards, Felix Schulze, Sarah M Sweet, Mathew R Varidel, Charlotte Welker

Abstract:

Large galaxy samples from multi-object IFS surveys now allow for a statistical analysis of the z~0 galaxy population using resolved kinematics. However, the improvement in number statistics comes at a cost, with multi-object IFS surveys more severely impacted by the effect of seeing and lower signal-to-noise. We present an analysis of ~1800 galaxies from the SAMI Galaxy Survey and investigate the spread and overlap in the kinematic distributions of the spin parameter proxy $\lambda_{Re}$ as a function of stellar mass and ellipticity. For SAMI data, the distributions of galaxies identified as regular and non-regular rotators with $kinemetry$ show considerable overlap in the $\lambda_{Re}$-$\varepsilon_e$ diagram. In contrast, visually classified galaxies (obvious and non-obvious rotators) are better separated in $\lambda_{Re}$ space, with less overlap of both distributions. Then, we use a Bayesian mixture model to analyse the $\lambda_{Re}$-$\log(M_*/M_{\odot})$ distribution. As a function of mass, we investigate whether the data are best fit with a single kinematic distribution or with two. Below $\log(M_*/M_{\odot})$~10.5 a single beta distribution is sufficient to fit the complete $\lambda_{Re}$ distribution, whereas a second beta distribution is required above $\log(M_*/M_{\odot})$~10.5 to account for a population of low-$\lambda_{Re}$ galaxies, presenting the cleanest separation of the two populations. We apply the same analysis to mock-observations from cosmological simulations. The mixture model predicts a bimodal $\lambda_{Re}$ distribution for all simulations, albeit with different positions of the $\lambda_{Re}$ peaks and with different ratios of both populations. Our analysis validates the conclusions from previous, smaller IFS surveys, but also demonstrates the importance of using kinematic selection criteria that are dictated by the quality of the observed or simulated data.

The SAMI Galaxy Survey: Towards an Optimal Classification of Galaxy Stellar Kinematics

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2021)

Authors:

Jesse van de Sande, Sam P Vaughan, Luca Cortese, Nicholas Scott, Joss Bland-Hawthorn, Scott M Croom, Claudia DP Lagos, Sarah Brough, Julia J Bryant, Julien Devriendt, Yohan Dubois, Francesco D'Eugenio, Caroline Foster, Amelia Fraser-McKelvie, Katherine E Harborne, Jon S Lawrence, Sree Oh, Matt S Owers, Adriano Poci, Rhea-Silvia Remus, Samuel N Richards, Felix Schulze, Sarah M Sweet, Mathew R Varidel, Charlotte Welker

Abstract:

Large galaxy samples from multi-object IFS surveys now allow for a statistical analysis of the z~0 galaxy population using resolved kinematics. However, the improvement in number statistics comes at a cost, with multi-object IFS surveys more severely impacted by the effect of seeing and lower signal-to-noise. We present an analysis of ~1800 galaxies from the SAMI Galaxy Survey and investigate the spread and overlap in the kinematic distributions of the spin parameter proxy $\lambda_{Re}$ as a function of stellar mass and ellipticity. For SAMI data, the distributions of galaxies identified as regular and non-regular rotators with $kinemetry$ show considerable overlap in the $\lambda_{Re}$-$\varepsilon_e$ diagram. In contrast, visually classified galaxies (obvious and non-obvious rotators) are better separated in $\lambda_{Re}$ space, with less overlap of both distributions. Then, we use a Bayesian mixture model to analyse the $\lambda_{Re}$-$\log(M_*/M_{\odot})$ distribution. As a function of mass, we investigate whether the data are best fit with a single kinematic distribution or with two. Below $\log(M_*/M_{\odot})$~10.5 a single beta distribution is sufficient to fit the complete $\lambda_{Re}$ distribution, whereas a second beta distribution is required above $\log(M_*/M_{\odot})$~10.5 to account for a population of low-$\lambda_{Re}$ galaxies, presenting the cleanest separation of the two populations. We apply the same analysis to mock-observations from cosmological simulations. The mixture model predicts a bimodal $\lambda_{Re}$ distribution for all simulations, albeit with different positions of the $\lambda_{Re}$ peaks and with different ratios of both populations. Our analysis validates the conclusions from previous, smaller IFS surveys, but also demonstrates the importance of using kinematic selection criteria that are dictated by the quality of the observed or simulated data.

The growth of density perturbations in the last $\sim$10 billion years from tomographic large-scale structure data

(2021)

Authors:

Carlos Garc铆a-Garc铆a, Jaime Ruiz Zapatero, David Alonso, Emilio Bellini, Pedro G Ferreira, Eva-Maria Mueller, Andrina Nicola, Pilar Ruiz-Lapuente

Probing galaxy bias and intergalactic gas pressure with KiDS Galaxies-tSZ-CMB lensing cross-correlations

Astronomy & Astrophysics EDP Sciences (2021)

Authors:

Z Yan, L van Waerbeke, T Tr枚ster, Ah Wright, D Alonso, M Asgari, M Bilicki, T Erben, S Gu, C Heymans, H Hildebrandt, G Hinshaw, N Koukoufilippas, A Kannawadi, K Kuijken, et al

Abstract:

We constrain the redshift dependence of gas pressure bias $\left\langle b_{y} P_{\mathrm{e}}\right\rangle$ (bias-weighted average electron pressure), which characterises the thermodynamics of intergalactic gas, through a combination of cross-correlations between galaxy positions and the thermal Sunyaev-Zeldovich (tSZ) effect, as well as galaxy positions and the gravitational lensing of the cosmic microwave background (CMB). The galaxy sample is from the fourth data release of the Kilo-Degree Survey (KiDS). The tSZ $y$ map and the CMB lensing map are from the {\textit{Planck}} 2015 and 2018 data releases, respectively. The measurements are performed in five redshift bins with $z\lesssim1$. With these measurements, combining galaxy-tSZ and galaxy-CMB lensing cross-correlations allows us to break the degeneracy between galaxy bias and gas pressure bias, and hence constrain them simultaneously. In all redshift bins, the best-fit values of $\bpe$ are at a level of $\sim 0.3\, \mathrm{meV/cm^3}$ and increase slightly with redshift. The galaxy bias is consistent with unity in all the redshift bins. Our results are not sensitive to the non-linear details of the cross-correlation, which are smoothed out by the {\textit{Planck}} beam. Our measurements are in agreement with previous measurements as well as with theoretical predictions. We also show that our conclusions are not changed when CMB lensing is replaced by galaxy lensing, which shows the consistency of the two lensing signals despite their radically different redshift ranges. This study demonstrates the feasibility of using CMB lensing to calibrate the galaxy distribution such that the galaxy distribution can be used as a mass proxy without relying on the precise knowledge of the matter distribution.