An empirical approach to model selection: weak lensing and intrinsic alignments

• Published in: Monthly notices of the Royal Astronomical Society Vol. 525; no. 2; pp. 1885 - 1901
• Main Authors: Campos, A, Samuroff, S, Mandelbaum, R
• Format: Journal Article
• Language: English
• Published: United States Oxford University Press 21.08.2023
• Subjects: ASTRONOMY AND ASTROPHYSICS; cosmological parameters; cosmology: observations; gravitational lensing: weak; methods: statistical; PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; methods: statistical; cosmology: observations; gravitational lensing: weak; cosmological parameters
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1093/mnras/stad2213

ABSTRACT In cosmology, we routinely choose between models to describe our data, and can incur biases due to insufficient models or lose constraining power with overly complex models. In this paper, we propose an empirical approach to model selection that explicitly balances parameter bias against model complexity. Our method uses synthetic data to calibrate the relation between bias and the χ2 difference between models. This allows us to interpret χ2 values obtained from real data (even if catalogues are blinded) and choose a model accordingly. We apply our method to the problem of intrinsic alignments – one of the most significant weak lensing systematics, and a major contributor to the error budget in modern lensing surveys. Specifically, we consider the example of the Dark Energy Survey Year 3 (DES Y3), and compare the commonly used non-linear alignment (NLA) and tidal alignment and tidal torque (TATT) models. The models are calibrated against bias in the Ωm–S8 plane. Once noise is accounted for, we find that it is possible to set a threshold Δχ2 that guarantees an analysis using NLA is unbiased at some specified level Nσ and confidence level. By contrast, we find that theoretically defined thresholds (based on, e.g. p-values for χ2) tend to be overly optimistic, and do not reliably rule out cosmological biases up to ∼1–2σ. Considering the real DES Y3 cosmic shear results, based on the reported difference in χ2 from NLA and TATT analyses, we find a roughly $30{{\ \rm per\ cent}}$ chance that were NLA to be the fiducial model, the results would be biased (in the Ωm–S8 plane) by more than 0.3σ. More broadly, the method we propose here is simple and general, and requires a relatively low level of resources. We foresee applications to future analyses as a model selection tool in many contexts.