Impact of cosmology dependence of baryonic feedback in weak lensing

Robust modeling of non-linear scales is critical for accurate cosmological inference in Stage IV surveys. For weak lensing analyses in particular, a key challenge arises from the incomplete understanding of how non-gravitational processes, such as supernovae and active galactic nuclei — collectively known as baryonic feedback — affect the matter distribution. Several existing methods for modeling baryonic feedback treat it independently from the underlying cosmology, an assumption which has been found to be inaccurate by hydrodynamical simulations. In this work, we examine the impact of this coupling between baryonic feedback and cosmology on parameter inference at LSST Y1 precision. We build mock 3×2pt data vectors using the Magneticum suite of hydrodynamical simulations, which span a wide range of cosmologies while keeping subgrid parameters fixed. We perform simulated likelihood analyses for two baryon mitigation techniques: (i) the Principal Component Analysis (PCA) method which identifies eigenmodes for capturing the effect baryonic feedback on the data vector and (ii) HMCode2020 [1] which analytically models the modification in the matter distribution using a halo model approach. Our results show that the PCA method is more robust than HMCode2020 with biases in Ω_m-S ₈ up to 0.3σ and 0.6σ, respectively, for large deviations from the baseline cosmology. For HMCode2020, the bias correlates with the input cosmology while for PCA we find no such correlation.