Beecroft Institute for Particle Astrophysics and Cosmology

Abstract:

Galactic scale tests have proven to be powerful tools in constraining fundamental physics in previously under-explored regions of parameter space. The astrophysical regime which they probe is inherently complicated, and the inference methods used to make these constraints should be robust to baryonic effects. Previous analyses have assumed simple empirical models for astrophysical noise without detailed calibration or justification. We outline a framework for assessing the reliability of such methods by constructing and testing more advanced baryonic models using cosmological hydrodynamical simulations. As a case study, we use the Horizon-AGN simulation to investigate warping of stellar disks and offsets between gas and stars within galaxies, which are powerful probes of screened fifth forces. We show that the degree of `U'-shaped warping of galaxies is well modelled by Gaussian random noise, but that the magnitude of the gas-star offset is correlated with the virial radius of the host halo. By incorporating this correlation we confirm recent results ruling out astrophysically relevant Hu-Sawicki $f(R)$ gravity, and identify a $\sim 30\%$ systematic uncertainty due to baryonic physics. Such an analysis must be performed case-by-case for future galactic tests of fundamental physics.

Abstract:

<jats:p>We measured the cross-correlation between galaxy weak lensing data from the Kilo Degree Survey (KiDS-1000, DR4) and cosmic microwave background (CMB) lensing data from the Atacama Cosmology Telescope (ACT, DR4) and the <jats:italic>Planck</jats:italic> Legacy survey. We used two samples of source galaxies, selected with photometric redshifts, (0.1 &lt; <jats:italic>z</jats:italic><jats:sub>B</jats:sub> &lt; 1.2) and (1.2 &lt; <jats:italic>z</jats:italic><jats:sub>B</jats:sub> &lt; 2), which produce a combined detection significance of the CMB lensing and weak galaxy lensing cross-spectrum of 7.7<jats:italic>σ</jats:italic>. With the lower redshift galaxy sample, for which the cross-correlation was detected at a significance of 5.3<jats:italic>σ</jats:italic>, we present joint cosmological constraints on the matter density parameter, Ω<jats:sub>m</jats:sub>, and the matter fluctuation amplitude parameter, <jats:italic>σ</jats:italic><jats:sub>8</jats:sub>, marginalising over three nuisance parameters that model our uncertainty in the redshift and shear calibration as well as the intrinsic alignment of galaxies. We find our measurement to be consistent with the best-fitting flat ΛCDM cosmological models from both <jats:italic>Planck</jats:italic> and KiDS-1000. We demonstrate the capacity of CMB weak lensing cross-correlations to set constraints on either the redshift or shear calibration by analysing a previously unused high-redshift KiDS galaxy sample (1.2 &lt; <jats:italic>z</jats:italic><jats:sub>B</jats:sub> &lt; 2), with the cross-correlation detected at a significance of 7<jats:italic>σ</jats:italic>. This analysis provides an independent assessment for the accuracy of redshift measurements in a regime that is challenging to calibrate directly owing to known incompleteness in spectroscopic surveys.</jats:p>

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.

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.