Extracting precise cosmology from weak lensing surveys requires modelling the nonlinear matter power spectrum which is suppressed at small scales due to baryonic feedback processes. However hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing WL and Atacama Cosmology Telescope DR5 kinematic SunyaevZeldovich kSZ to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model baryonification. First using WL only we compare the inlineformulatexmath idTM0001 notationLaTeXS8texmathinlineformula constraints using baryonification to a simulationcalibrated halo model a simulationbased emulator model and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of inlineformulatexmath idTM0002 notationLaTeXS8texmathinlineformula and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters with the joint WL  kSZ analysis constraining inlineformulatexmath idTM0003 notationLaTeXS80.8230.0190.020texmathinlineformula. We measure the suppression of the nonlinear matter power spectrum using WL  kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from Xray observations and simulation predictions. We conclude that the WL  kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.
['Astrophysics - Cosmology and Nongalactic Astrophysics', 'Astrophysics - Astrophysics of Galaxies']
Weak lensing combined with the kinetic SunyaevZeldovich effect a study of baryonic feedback
{"Weak lensing combined with the kinetic Sunyaev Zel'dovich effect: A study of baryonic feedback": "L. Bigwood, 1 , 2 ★ A. Amon, 3 , 1 , 2 † A. Schneider, 4 J. Salcido, 5 I. G. McCarthy, 5 C. Preston, 1 D. Sanchez, 6 D. Sijacki, 1 , 2 E. Schaan, 7 , 8 S. Ferraro, 9 , 10 N. Battaglia, 11 A. Chen, 12 , 13 S. Dodelson, 14 , 15 A. Roodman, 16 , 17 A. Pieres, 18 , 19 A. Fert'e, 17 A. Alarcon, 20 , 21 A. Drlica-Wagner, 22 , 23 , 24 A. Choi, 25 A. Navarro-Alsina, 26 A. Campos, 14 , 15 A. J. Ross, 27 A. Carnero Rosell, 28 , 18 B. Yin, 14 B. Yanny, 23 C. S'anchez, 29 C. Chang, 22 , 24 C. Davis, 16 C. Doux, 29 , 30 D. Gruen, 31 E. S. Rykoff, 16 , 17 E. M. Huff, 32 E. Sheldon, 33 F. Tarsitano, 34 F. Andrade-Oliveira, 12 G. M. Bernstein, 29 G. Giannini, 35 , 24 H. T. Diehl, 23 H. Huang, 36 , 37 I. Harrison, 38 I. Sevilla-Noarbe, 39 I. Tutusaus, 40 J. Elvin-Poole, 41 J. McCullough, 16 J. Zuntz, 42 J. Blazek, 43 J. DeRose, 44 J. Cordero, 45 J. Prat, 22 , 46 J. Myles, 3 K. Eckert, 29 K. Bechtol, 47 K. Herner, 23 L. F. Secco, 24 M. Gatti, 29 M. Raveri, 48 M. Carrasco Kind, 49 , 50 M. R. Becker, 20 M. A. Troxel, 51 M. Jarvis, 29 N. MacCrann, 52 O. Friedrich, 53 O. Alves, 12 P.-F. Leget, 16 R. Chen, 51 R. P. Rollins, 45 R. H. Wechsler, 54 , 16 , 17 R. A. Gruendl, 49 , 50 R. Cawthon, 55 S. Allam, 23 S. L. Bridle, 45 S. Pandey, 29 S. Everett, 32 T. Shin, 56 W. G. Hartley, 57 X. Fang, 58 , 36 Y. Zhang, 59 M. Aguena, 18 J. Annis, 23 D. Bacon, 60 E. Bertin, 61 , 62 S. Bocquet, 31 D. Brooks, 63 J. Carretero, 35 F. J. Castander, 64 , 21 L. N. da Costa, 18 M. E. S. Pereira, 65 J. De Vicente, 39 S. Desai, 66 P. Doel, 63 I. Ferrero, 67 B. Flaugher, 23 J. Frieman, 23 , 24 J. Garc'ıa-Bellido, 68 E. Gaztanaga, 64 , 60 , 21 G. Gutierrez, 23 S. R. Hinton, 69 D. L. Hollowood, 70 K. Honscheid, 27 , 71 D. Huterer, 12 D. J. James, 72 K. Kuehn, 73 , 74 O. Lahav, 63 S. Lee, 32 J. L. Marshall, 75 J. Mena-Fern'andez, 76 R. Miquel, 77 , 35 J. Muir, 78 M. Paterno, 23 A. A. Plazas Malag'on, 16 , 17 A. Porredon, 79 A. K. Romer, 80 S. Samuroff, 43 E. Sanchez, 39 D. Sanchez Cid, 39 M. Smith, 81 M. Soares-Santos, 82 , 12 E. Suchyta, 83 M. E. C. Swanson, 49 G. Tarle, 12 C. To, 27 N. Weaverdyck, 58 , 44 J. Weller, 84 , 85 P. Wiseman, 81 and M. Yamamoto 51  \nThe authors' affiliations are shown at the end of this paper.  \n10 April 2024", 'ABSTRACT': "Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev-Zel'dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, 'baryonification'. First, using WL only, we compare the 𝑆 8 constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of 𝑆 8 and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters and shifts 𝑆 8 to 𝑆 8 = 0 . 823 + 0 . 019 -0 . 020 , a higher value than attained using the WL-only analysis. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.  \nKeywords: cosmology: observations - gravitational lensing - large-scale structure of Universe", '1 INTRODUCTION': "The standard cosmological model, Λ CDM, has been very successful when tested against observations of the cosmic microwave background (CMB; Planck Collaboration 2020a; Bennett et al. 2003), the lensing of the CMB at intermediate redshifts (Planck Collaboration 2020b; Madhavacheril et al. 2023) and low-redshift observations of the expansion history as probed by baryon acoustic oscillations (Eisenstein et al. 2005). Measurements of weak galaxy lensing provide a strong test of Λ CDM at relatively small scales in the low-redshift Universe. In order to extract unbiased cosmological constraints from weak lensing, accurate modelling of the non-linear matter distribution at 𝑘 > 0 . 1 ℎ Mpc -1 is crucial. This requires understanding both the non-linear dark matter evolution due to gravity to percent-level accuracy, as well as the impact of the baryons. In particular, a number of physical processes associated with baryons redistribute gas and impact the non-linear matter power spectrum by up to ∼ 30% (see e.g. Chisari et al. 2019; van Daalen et al. 2020). This 'baryonic feedback' encapsulates the energetic effect of Active Galactic Nuclei (AGN) heating gas and ejecting it to the outskirts of groups and clusters, as well as the likely sub-dominant effects of stellar winds, supernovae feedback and gas cooling (van Daalen et al. 2011). At present, astrophysical model uncertainties, such as those due to baryonic feedback, have been shown to limit the precision of weak lensing surveys (Amon et al. 2022; DES & KiDS Collaborations 2023). Therefore, to extract maximal information from weak-lensing data demands improved modelling of baryonic effects and their impact on the matter distribution. Moreover, it has been proposed that the ' 𝑆 8 tension' - the finding that weak lensing constraints on the clustering amplitude parameter, 𝑆 8 = 𝜎 8 ( Ω m / 0 . 3 ) 0 . 5 , are lower than predictions from the CMB could be explained by a suppression of the non-linear matter power spectrum(Amon & Efstathiou 2022; Preston et al. 2023). This could either be caused by a more extreme baryonic feedback effects than hydrodynamical simulations predict or extensions to the standard modelofcosmology.Inordertoisolateadeparturefromthestandard cosmological model, baryonic effects must be better understood.  \nPowerful AGN feedback is believed to have the ability to eject baryons beyond the virial radius of galaxies, redistributing the gas to the outskirts of galaxy groups and clusters (McCarthy et al. 2011; Springel et al. 2018; Dubois et al. 2016; Henden et al. 2018). Indeed, studies of hydrodynamical simulations have demonstrated that AGN feedback alters the total matter distribution relative to dark matter-only simulations, and that it causes a suppression of the power spectrum at scales 0 . 1 ≲ 𝑘 ≲ 10 ℎ Mpc -1 , whereas increased star formation can enhance power on the smallest scales (see Chisari et al. 2019, for a review). These simulations reproduce many of the observed properties of galaxies, including optical properties, galaxy group/cluster profiles, scaling relations and Sunyaev-Zel'dovich counts. However, despite these successes, the scale, amplitude and redshift dependence of the larger-scale power suppression remain largely uncertain, with significant variation between simulations. These differences are direct outcomes of the 'sub-grid' modelling of astrophysical processes, which take place below the resolution scale of the simulation.  \nSpecifically, sub-grid models are required to follow the formation, growth and feedback of black holes, as well as gas cooling, metal enrichment, star formation and associated stellar feedback. The AGN feedback may operate in either kinetic or thermal modes (generally associated with the radio and quasar modes respectively), or alternate between the two depending on the black hole accretion rate in a 'two-mode' feedback scenario (see Sijacki et al. 2007). The  \nfeedback model normally has a number of associated ill-constrained parameters encoding the feedback strength, such as the efficiency of thermal/kinetic coupling, the black hole accretion rate, and, in some models, the minimum heating temperature of gas cells before a feedback event occurs. While physical arguments can be used to narrow the plausible range of some of these parameters, this is normally not sufficiently constraining for precision cosmology purposes. Thus, it is often the case that the parameters are calibrated against key observables.  \nX-ray measurements of the hot gas fractions of groups and clusters within the virial radius are most widely used to benchmark the simulations, along with galaxy stellar mass function (GSMF), star formation history and galaxy sizes (McCarthy et al. 2017; Henden et al. 2018; Schaye et al. 2023; Kugel et al. 2023; Nelson et al. 2019). Even with identical subgrid physics, the box size and resolution of the simulation can also have a non-negligible impact on the matter distribution (van Daalen et al. 2011; Pandey et al. 2023). Indeed, there is a large parameter space of feedback prescriptions, modelling choices and simulation properties that result in significant variation in the suppression of the matter power spectrum (Salcido et al. 2023; Schaye et al. 2023; Hern'andez-Aguayo et al. 2023; Dav'e et al. 2019; Henden et al. 2018).  \nRecent weak lensing analyses have devised various approaches to mitigate the impact of baryonic feedback on cosmological constraints. The DES Y3 cosmic shear analysis opted to discard measurements from the analysis on angular scales that are impacted by baryonic effects from the analysis (Krause et al. 2021; Amon et al. 2023; Secco, Samuroff et al. 2022). Alternatively, baryon feedback has been modelled using a halo model approach (Asgari et al. 2021; Li et al. 2023), and using a halo model that is calibrated to a hydrodynamical simulations (Mead et al. 2021; DES & KiDS Collaborations 2023). More recently, Salcido et al. (2023) have developed an emulator trained using a suite of hydrodynamical simulations with varied feedback efficiencies.  \nInstead of relying on the hydrodynamical simulations, the baryonification model is another approach which shifts particle outputs in gravity-only simulations to attain modified halo profiles, modelling the rearrangement of baryon material caused by feedback effects (Schneider & Teyssier 2015). This approach has been used in Schneider et al. (2022), Chen et al. (2023) and Aric'o et al. (2023). Other approaches include a principal component analysis (Huang et al. 2019).  \nAn alternative approach to using models that are informed by hydrodynamical simulations (and therefore indirectly benchmarked against X-ray data) is to jointly analyse weak lensing data with observations of the gas content in and around galaxy groups and clusters. This has been done to improve cosmological constraints by reducing the model-space of the nuisance parameters (Troster et al. 2020) and to constrain the suppression of the matter power spectrum (Schneider et al. 2022).  \nA highly complementary observable to X-ray is the kinetic Sunyaev-Zel'dovich (kSZ) effect, caused by the Thomson scattering of the CMB photons by free electrons moving with bulk motion in groups and clusters of galaxies relative to the CMB rest-frame (Sunyaev & Zeldovich 1980, 1972). This causes a shift in the CMB temperature while preserving the frequency dependence. If the bulk line-of-sight velocity is known, the kSZ effect directly measures the free electron number density, independent of temperature (Hand et al. 2012; Planck Collaboration 2016; Soergel et al. 2016; Schaan et al. 2021). The kSZ effect is well-suited to probe low density and low temperature environments like the outskirts of galaxies  \nand clusters, whereas X-ray measurements are more sensitive to the inner regions (Amodeo et al. 2021).  \nThe goals of this work are threefold:  \n- I. Wetest the performance of four baryon feedback mitigation strategies for analysing mock and DES Y3 weak lensing data: the DES Y3 scale cut approach (Krause et al. 2021), a halo model approach (Mead et al. 2021) calibrated to a hydrodynamical simulation, an emulator built using a suite of hydrodynamical simulations (Salcido et al. 2023) and the baryonification model (Schneider & Teyssier 2015; Schneider et al. 2019).\n- II. As an alternative to simulation-driven models, we use the most flexible model, baryonification, to jointly analyse the lensing data with the ACT kSZ measurements for improved constraints on the baryonic feedback parameters and therefore the cosmological parameters.\n- III. We use these data to constrain astrophysical observables for the first time, providing a new avenue to benchmark the hydrodynamical simulations.  \nThe paper is structured as follows. Section 2 describes the DES Y3 cosmic shear and ACT DR5 kSZ datasets used in this analysis. Section 3 outlines the modelling of the cosmic WL and kSZ measurements, including the four baryon models.  \nIn Section 4 we summarise the findings of a mock analysis. We compare the cosmological constraints when analysing the shear data with different baryon mitigating models and model complexities in Section 5. In Section 6 we present our constraints on both cosmological and baryonic parameters obtained from a joint WL + kSZ analysis. Finally, in Section 7 we consider our constraints on the observables that simulations benchmark against. We summarise key findings and conclude in Section 8.", '2.1 Dark Energy Survey cosmic shear': "The Dark Energy Survey (DES) has completed six years of photometric observations in the grizY bands, using the 4-metre Blanco telescope located at the Cerro Tololo Inter-American Observatory. The survey spans ∼ 5000deg 2 in the Southern Hemisphere.  \nFor this analysis, we use data taken during the survey's first three years of operation (DES Y3), between 2013 and 2016 (SevillaNoarbe et al. 2021). The DES Y3 footprint covers 4143 deg 2 with a number density of 5 . 59 galaxies arcmin -2 to a depth of 𝑖 ∼ 23.5. The /m.pc/e.pc/t.pc/a.pc/c.pc/a.pc/l.pc/i.pc/b.pc/r.pc/a.pc/t.pc/i.pc/o.pc/n.pc weak lensing catalog has over 100 million galaxies that have passed a raft of validation tests Gatti, Sheldon et al. 2022. The sample has been divided into four redshift bins and the calibrated redshift distributions and associated uncertainty are defined in Myles, Alarcon et al. (2021). Remaining biases in the shape measurement and redshift distributions, primarily due to blending, are calibrated using image simulations, and the associated corrections for each redshift bin are reported in MacCrann et al. (2022). The DES Y3 cosmic shear tomographic two-point correlation functions, 𝜉 ± , are measured in twenty angular logarithmic bins spanning 2.5 to 250.0 arcmin (Amon et al. 2022; Secco, Samuroff et al. 2022).", "2.2 Atacama Cosmology Telescope Kinetic Sunyaev Zel'dovich": "The Atacama Cosmology Telescope (ACT) is a 6-metre millimeter waveband telescope, observing thhe cosmic microwave background  \nFigure 1. The suppression of the matter power spectrum due to baryonic effects predicted by hydrodynamical simulations. We show OWLS-AGN (van Daalen et al. 2011; Schaye et al. 2010, dark green solid line), which informed the scale cuts baryonic feedback mitigation approach. We also plot the matter power spectrum suppression measured in the 400 simulations of the ANTILLESsuite, which calibrated the SP(k) emulator (Salcido et al. 2023, light green region). The range in the suppression predicted by the BAHAMAS suite spanning Θ AGN = 7 . 6 -8 . 0 from which HM20 was calibrated (McCarthy et al. 2017, dark blue region). We further plot FLAMINGO (Schaye et al. 2023, black solid line); BAHAMAS (McCarthy et al. 2017, black dashdotted line); SIMBA (Dav'e et al. 2019, black dotted line), MillenniumTNG (Pakmor et al. 2023, black dashed line) and FABLE (doubledot-dashed line, Henden et al. 2018, Bigwood et al. in prep.). Finally, we show in red the three baryonic feedback scenarios we test in our mock analysis (Section 4.2): BAHAMAS 𝑇 heat = 8 . 0 (dotted red line), cosmo-OWLS 𝑇 heat = 8 . 5 (dashed red line) and 𝐴 mod = 0 . 82 (solid red line) (Preston et al. 2023).  \n<!-- image -->  \n(CMB). Since first light in 2007, it has had three generations of receivers, the most recent of which is the polarisation-sensitive Advanced ACTPol (AdvACT), which extended the frequency coverage to five bands spanning 28 to 230 GHz. The fifth data release (herinafter DR5) coadds maps collected from 2008 to 2018 covering approximately 18000 deg 2 and utilises data from all three generations of receiver (Naess et al. 2020).  \nThis work uses the kinetic Sunyaev Zel'dovich (kSZ) measurements presented in Schaan et al. (2021). These are stacked measurements of the ACT DR5 and Planck CMBtemperature maps at 98 GHz (hereinafter called f90 for consistency with Schaan et al. (2021)) and 150 GHz (f150) with the reconstructed velocities of the spectroscopic BOSS CMASS galaxy catalog. The galaxy sample spans the redshift range 0 . 4 < 𝑧 < 0 . 7 with a median redshift of 𝑧 = 0 . 55. It corresponds to a selection of relatively massive galaxies with a mean stellar mass of log 10 ( 𝑀 star / M ⊙ ) ≈ 11 . 3 and an assumed mean halo mass of ∼ 10 13 M ⊙ (Schaan et al. 2021; Amodeo et al. 2021), though the latter is quite uncertain. Given its importance to the modelling of the stacked kSZ sample, we will discuss the mean halo mass of this sample later in the work.", '2.3 Hydrodynamical simulations': "Throughout the paper, we compare our constraints on the matter power spectrum to predictions from a range of hydrodynamical simulations: the FLAMINGO ( 1 Gpc ) 3 box with baryonic  \nparticle masses of 10 9 𝑀 ⊙ (Schaye et al. 2023, solid line); BAHAMAS ( 400 Mpc / h ) 3 box with Θ AGN = 7 . 8 and 1024 3 dark matter and baryonic particles (McCarthy et al. 2017, dash-dotted line); SIMBA ( 100 Mpc / h ) 3 box with 1024 3 gas elements (Dav'e et al. 2019, dotted line); MillenniumTNG 740 Mpc 3 with baryonic mass resolution of 3 . 1 × 10 7 𝑀 ⊙ (Pakmor et al. 2023, dashed line) and FABLE ( 100 Mpc / h ) 3 box with 1280 3 dark matter particles and 1280 3 baryonic resultion elements (doubledotted-dashed line, Henden et al. 2018, Bigwood et al. in prep.). . These simulations not only span a range of box sizes and resolutions, but also feedback implementations, hydrodynamical schemes and calibration strategies. Fig. 1 shows the suppression of the non-linear matter power spectrum due to baryonic effects, 𝑃 ( 𝑘 )/ 𝑃 DMonly , as predicted by each simulation, demonstrating the spread in the predicted amplitude and scale-dependence. Several additional hydro simulations are also used to calibrate the baryonic feedback models tested throughout this work. We therefore also plot the prediction from the OWLS-AGN ( 100 Mpc / h ) 3 box with 512 3 dark matter and baryonic particles respectively (van Daalen et al. 2011; Schaye et al. 2010, dark green solid line), the span in the suppression predicted by BAHAMAS when Θ AGN is varied within the range 7.6-8.0 (see Section 3.2.2), and range modification measured by the ANTILLES suite of 400 simulations (Salcido et al. 2023, light green region), each with box size ( 100 Mpc / h ) 3 and 256 3 dark matter and baryonic particles respectively. ANTILLES spans significantly more extreme feedback scenarios than the other simulations considered, such that baryonic effects impact the power spectrum with greater amplitude and at smaller 𝑘 scales.", '3.1 Cosmic shear signal': "Theshear two-point correlation functions, 𝜉 ± ( 𝜃 ) , for a given angular separation, 𝜃 , computed for redshift bins 𝑖, 𝑗 , can be related to the 3D non-linear matter power spectrum. First, it can be expressed as a decomposition of the 2D convergence power spectrum 𝐶 𝑖, 𝑗 𝜅 ( ℓ ) at an angular wave number, ℓ , as  \n𝜉 ± ( 𝜃 ) = ∑︁ ℓ 2 ℓ + 1 4 𝜋 𝐺 ± ℓ ( cos 𝜃 ) [ 𝐶 𝑖 𝑗 𝜅, EE ( ℓ ) ± 𝐶 𝑖 𝑗 𝜅, BB ( ℓ )] . (1)  \nWenote that weak gravitational lensing does not produce B-modes. However, we show the more general expression here, as a B-mode contribution from intrinsic alignments is possible. The functions 𝐺 ± ℓ are computed from the Legendre polynomials following Stebbins (1996).  \nUnder the Limber and flat-sky approximations (Limber 1953; LoVerde & Afshordi 2008), we can relate 𝐶 𝑖, 𝑗 𝜅 ( ℓ ) to the 3D nonlinear matter power spectrum, 𝑃 , via  \n𝐶 𝑖 𝑗 𝜅 ( ℓ ) = ∫ 𝜒 H 0 𝑑𝜒 𝑊 𝑖 ( 𝜒 ) 𝑊 𝑗 ( 𝜒 ) 𝜒 2 𝑃 GLYPH<18> 𝑘 = ℓ + 0 . 5 𝜒 ( 𝑧 ) , 𝑧 GLYPH<19> . (2)  \nHere, we assume a spatially flat Universe, 𝜒 is the comoving angular diameter distance and 𝜒 H is the horizon distance. 𝑊 𝑖 ( 𝜒 ) are the lensing efficiency kernels, given by:  \n𝑊 𝑖 ( 𝜒 ) = 3 𝐻 2 0 Ω m 2 𝑐 2 𝜒 𝑎 ( 𝜒 ) ∫ 𝜒 H 𝜒 𝑑𝜒 ' 𝑛 𝑖 ( 𝜒 ' ) 𝜒 ' -𝜒 𝜒 ' , (3)  \nwhere 𝑎 ( 𝜒 ) is the scale factor at comoving distance 𝜒 , and 𝑛 𝑖 ( 𝜒 ) 𝑑𝜒  \nis the effective number density of galaxies in 𝑑𝜒 , normalised to unity.  \nIn this analysis, the linear matter power spectrum is calculated using CAMB (Lewis et al. 2000), and the non-linear correction is determined using HMC/o.pc/d.pc/e.pc2020 Mead et al. (2021). We refer the reader to DES & KiDS Collaborations (2023) for a detailed comparison of cosmic shear analysis choices, which we use as a guide to formulate the baseline model in this work. We assume three neutrino species with two massless states and one massive state with a mass fixed at the minimum mass allowed by oscillation experiments, 𝑚 𝜈 = 0 . 06eV (Patrignani 2016).  \nThe intrinsic alignment (IA) of galaxies with their local environment also contributes to the shear correlation function and must be modelled. We chose to do this using the non-linear linearalignment model (NLA), which describes the linear tidal alignment of galaxies with the density field (Hirata & Seljak 2004), with a non-linear correction to the linear matter power spectrum (Bridle & King 2007). This approach requires two additional free nuisance parameters: 𝐴 IA , modulating the amplitude of the intrinsic alignment model (see equation [3-5] in Bridle & King 2007 for the NLA intrinsic alignment power spectra, 𝐶 GI and 𝐶 II ) and a redshift-dependence parameter, using a power-law with [( 1 + 𝑧 )/( 1 + 0 . 62 )] 𝜂 IA .  \nWe model the uncertainty in mean redshift and the shear calibration for each redshift bin 𝑖 as the free parameters Δ 𝑧 𝑖 and 𝑚 𝑖 , respectively, following Amon et al. (2022) and Secco, Samuroff et al. (2022), and preserve values of the uncertainty determined by Myles, Alarcon et al. (2021) and MacCrann et al. (2022). In the cases where the small-scale measurements are analysed, we refer the reader to Chen et al. (2023) for validation that higher-order cosmic shear modelling corrections remain subdominant.", '3.2 Modelling baryonic feedback for cosmic shear': 'Strategies have been devised to account for baryonic effects when analysing weak lensing measurements in order to extract unbiased cosmological information. In this section we outline the four approaches we investigate in this work: (1) restrict angular scales, (2) a halo model approach, (3) a hydrodynamical simulation-based emulator and (4) an analytical N-body simulation model. In Table 1 we summarise the free parameters of each baryon feedback model and prior choices.', '3.2.1 Restricting angular scales': "The DES Y3 cosmic shear analysis mitigates the impact of baryonic effects by discarding the measurements at small angular scales (Krause et al. 2021; Amon et al. 2023; Secco, Samuroff et al. 2022) and analysing the data assuming a dark matter-only model. In this work, we adopt their ' Λ CDM Optimised' scale cuts, which were designed to minimise the bias on Λ CDM cosmological parameters due to unmodelled baryonic effects to be less than 0.3 𝜎 2D in the Ω m -𝑆 8 parameter space, for a joint lensing and clustering analysis. For a Λ CDM analysis of cosmic shear alone, this corresponds to up to a 0.14 𝜎 2D potential bias. Note that to determine the angular scales to be used, the OWLS-AGN hydrodynamical simulation (van Daalen et al. 2011; Schaye et al. 2010) was chosen as a representative feedback scenario (shown as the green line, Fig 1). Synthetic cosmic shear data vectors were contaminated according to  \n𝑃 𝑏 ( 𝑘, 𝑧 ) = 𝑃 hydro ( 𝑘, 𝑧 ) 𝑃 DMO ( 𝑘, 𝑧 ) 𝑃 ( 𝑘, 𝑧 ) , (4)  \nTable 1. Qualitative descriptions of the parameters associated with the HM20, SP(k) and BCEmu methods, which each model the impact of baryonic feedback on the non-linear power spectrum. We also show the priors utilised on each parameter in this analysis, and the conservative prior alternative if applicable. U [ ] brackets indicate flat uniform priors within the range shown. N( ) brakets indicate Gaussian priors described by their mean and 1 𝜎 width. The BCEmu1 model has log 10 𝑀 c as the only free parameter, BCEmu3 varies log 10 𝑀 c , 𝜃 ej and 𝜂 cga , and all seven parameters are free in BCEmu7. In the case of using the reduced model complexity of BCEmu1 or BCEmu3 models, we also show the values that the BCEmu parameters are fixed to. Note that the 𝑀 200 mass parameter is only included in the BCEmu joint WL + kSZ analysis and the choice of the fixed value is adopted from Schaan et al. (2021).  \nwhere 𝑃 hydro ( 𝑘, 𝑧 ) and 𝑃 DMO ( 𝑘, 𝑧 ) are the full hydrodynamical and dark matter-only matter power spectra from the OWLS-AGN suite and 𝑃 ( 𝑘, 𝑧 ) is the non-linear matter power spectrum.  \nA benefit of this approach is that it is agnostic to the exact shape and physics of the matter power spectrum suppression, once the feedback in the Universe is lower in amplitude and scale extent than the simulation chosen (in this case, OWLS-AGN). However, this approach misses the opportunity to extract high-signal-to-noise information about the underlying cosmological model and the astrophysical effects.", '3.2.2 Halo-model approach: HM20': "HM/c.pc/o.pc/d.pc/e.pc2020, hereinafter HM20, models the non-linear power spectrum and includes a free parameter to modulate the amount of baryonic feedback, Θ AGN = log 10 ( 𝑇 AGN / 𝐾 ) (Mead et al. 2021). This parameter scales the halo concentration and the stellar and gas content, leading to a modification in the overall amplitude and shape of the 'one-halo' term in the halo model. The model is calibrated to fit the power spectrum 'response' (the matter-matter power spectrum divided by the same measurement in an equivalent dark matter-only box) of the BAHAMAS hydrodynamical simulations (McCarthy et al. 2017; van Daalen et al. 2020) in the range Θ AGN = [ 7 . 6 -8 . 0 ] (blue shaded region, Fig. 1). We note that 𝑇 AGN is related to Δ 𝑇 heat , which is the BAHAMAS subgrid heating parameter, where an AGN  \nfeedback event will only occur after the black hole has stored sufficient energy to heat a fixed number of gas particles by Δ 𝑇 heat .  \nWe define two prior ranges for this case. The first spans the range of Θ AGN values that bracket the BAHAMAS Θ AGN = 7 . 6 -8 . 0 simulations that the model was calibrated against. The 'wide prior' chosen here to be Θ AGN = 7 . 3 -9 . 0 extends beyond the calibration range to span more extreme scenarios and allow for a dark matteronly case.  \nThis approach has been shown to be accurate at the level of < 2 . 5%to 𝑘 < 10 ℎ Mpc -1 (Meadetal.2021)whenfittingsimulated power spectra at 𝑧 < 2 with a range of cosmologies thus allowing all measured angular scales of the DES Y3 dataset to be utilised (in this case, to 2 . 5 arcmin). A downside of this model is that it relies on the accuracy of the specific feedback implementation and predicted power suppression of a particular simulation and may not be flexible enough to capture the true scenario.", '3.2.3 Hydrodynamical simulation emulator: SP(k)': 'SP(k) is a flexible empirical model trained on the ANTILLES suite of 400 cosmological hydrodynamical simulations (plotted as the green lines in Fig. 1, Salcido et al. 2023). The model predicts the power spectrum suppression given the baryon fraction-halo mass relation of galaxy groups and clusters as input, building upon previous insights from van Daalen et al. (2020). The ANTILLES', 'L. Bigwood, A. Amon et al.': "suite span a range of feedback scenarios, allowing the emulator to achieve a ∼ 2% level accuracy to describe baryonic effects at scales of up to 𝑘 ≲ 10 ℎ Mpc -1 and redshifts up to 𝑧 = 3.  \nIn particular, SP(k) casts the suppression in terms of the baryon fraction at the optimal mass , ˆ 𝑀 𝑘 , defined as the halo mass that maximizes the strength of the correlation between the suppression of the total matter power spectrum and the total baryon fraction of haloes of different mass. It uses an exponential plateau function to model for the fractional impact of baryons as,  \n𝑃 hydro ( 𝑘 )/ 𝑃 DM ( 𝑘 ) = 𝜆 ( 𝑘, 𝑧 )-[ 𝜆 ( 𝑘, 𝑧 ) -𝜇 ( 𝑘, 𝑧 )] exp GLYPH<2> -𝜈 ( 𝑘, 𝑧 ) ˜ 𝑓 b GLYPH<3> , (5)  \nwhere ˜ 𝑓 b is the baryon fraction at the optimal halo mass normalised by the universal baryon fraction, i.e.:  \n˜ 𝑓 b = 𝑓 b ( ˆ 𝑀 𝑘, SO ( 𝑘, 𝑧 ))/( Ω b / Ω m ) , (6)  \nand ˆ 𝑀 𝑘, SO , 𝜆 ( 𝑘, 𝑧 ) , 𝜇 ( 𝑘, 𝑧 ) and 𝜈 ( 𝑘, 𝑧 ) are functions with best-fit parameters given in Salcido et al. (2023).  \nFor the mass range that can be relatively well probed in current X-ray and Sunyaev-Zel'dovich effect observations (10 13 ≲ 𝑀 200 [ M ⊙ ] ≲ 10 15 ), the total baryon fraction of haloes can be roughly approximated by a power-law with constant slope (e.g. Mulroy et al. 2019; Akino et al. 2022). Salcido et al. (2023) find that a modified version of the functional form 1 presented in Akino et al. (2022) provides a reasonable agreement with simulations up to redshift 𝑧 = 1,  \n𝑓 b /( Ω b / Ω m ) = GLYPH<18> 𝑒 𝛼 100 GLYPH<19> GLYPH<18> 𝑀 500 𝑐 10 14 M ⊙ GLYPH<19> 𝛽 -1 GLYPH<18> 𝐸 ( 𝑧 ) 𝐸 ( 0 . 3 ) GLYPH<19> 𝛾 , (7)  \nwhere 𝛼 sets the power-law normalisation, 𝛽 sets power-law slope, 𝛾 provides the redshift dependence, 𝐸 ( 𝑧 ) is the dimensionless Hubble parameter and 𝑓 b is the baryon fraction measured within 𝑅 500 . We use this function to facilitate the marginalization over the uncertainties in the observed baryon fractions, introducing 𝛼 , 𝛽 and 𝛾 as free parameters in our WL analysis.  \nSP(k) effectively depends only on a single physicallymeaningful parameter, i.e. the baryon fraction. The benefit of this approach is that observational constraints on the baryon fraction could be used to inform the priors used in cosmological analysis. For our study we use two different sets of priors to marginalize over the uncertainties in the observed baryon fractions: wide (conservative) priors consistent with the range of feedback models probed by the ANTILLES simulations used to calibrate the SP(k) model (Salcido et al. 2023), and 'observational' priors that encompass current observational constraints on the baryon fraction - halo mass relation from Akino et al. (2022). Table 1 reports the priors on 𝛼 , 𝛽 and 𝛾 for the two choices.  \nAs with the scale cuts and HM20 model, we caveat that the SP(k) method is limited by the specific feedback implementation used in the hydrodynamical simulation it was trained on, in this case the ANTILLES suite. While ANTILLES is currently the largest and widest suite of hydro simulations in terms of feedback variations, it still may not cover all possible baryonic responses (see e.g. Appendix B4).", '3.2.4 Analytical N-body simulation model: Baryonification': "Baryonification is a method for including the effects of baryonic feedback in dark matter-only 𝑁 -body simulations based on perturbative shifts of particles that mimic the effects of feedback at  \ncosmological scales (Schneider & Teyssier 2015; Schneider et al. 2019; Aric'o et al. 2020). The particles are shifted in order to obtain modified halo profiles that include the presence of gas and stars which are shaped by feedback effects. We provide a summary of the method including some important aspects of the parametrisation and refer the reader to Schneider et al. (2019) and Giri & Schneider (2021) for a more detailed explanation.  \nIn practice, the baryonification (bfc) method relies on a modification of profiles via spherically symmetrical particle shifts around the centres of 𝑁 -body haloes with NFW-like profiles, 𝜌 nfw , following  \n𝜌 nfw ( 𝑟 ) → 𝜌 bfc ( 𝑟 ) = 𝜌 clm ( 𝑟 ) + 𝜌 gas ( 𝑟 ) + 𝜌 cga ( 𝑟 ) . (8)  \nThe final baryonified profiles, 𝜌 bfc , consist of three components: the collisionless matter, gas and central galaxy. The collisionless matter ( 𝜌 clm ) profile is dominated by dark matter but also contains satellite galaxies and halo stars. Its shape is modified with respect to the original NFW shape via adiabatic contraction and expansion (Teyssier et al. 2011). The central galaxy profile, 𝜌 cga , is parametrised as a power law with an exponential cutoff. This component affects only the innermost part of the halo, rather than cosmological scales.  \nAt cosmological scales, baryonic effects are primarily caused by feedback-induced changes of the gas distribution around haloes, described by the gas profile, 𝜌 gas . Motivated by X-ray observations (Eckert et al. 2016; Schneider et al. 2019), these effects are parametrised with five model parameters as  \n𝜌 gas ( 𝑟 ) ∝ Ω b / Ω m -𝑓 star ( 𝑀 ) h 1 + GLYPH<16> 𝑟 𝑟 core GLYPH<17> i 𝛽 ( 𝑀 ) h 1 + GLYPH<16> 𝑟 𝑟 ej GLYPH<17> 𝛾 i 𝛿 -𝛽 ( 𝑀 ) 𝛾 , (9)  \nwhich consists of a cored power-law profile with a truncation at the ejection radius, 𝑟 ej = 𝜃 ej 𝑅 200 , where 𝜃 ej is a dimensionless free model parameter 2 . The shape of the truncation beyond the ejection radius is controlled by the 𝛾 and 𝛿 parameters, where the former defines the abruptness and the latter defines the slope of the function. The core of the profile is fixed to 𝑟 core = 𝜃 core 𝑅 200 with 𝜃 core = 0 . 1. The slope of the power law 𝛽 ( 𝑀 ) is a function that varies with halo mass 3 and is parameterised as  \n𝛽 ( 𝑀 ) = 3 ( 𝑀 / 𝑀 c ) 𝜇 1 + ( 𝑀 / 𝑀 c ) 𝜇 . (10)  \nThe halo mass-dependence of 𝛽 accounts for the fact that AGN feedback is more efficient in removing gas around galaxy groups while large clusters tend to keep most of their gas inside the virial radius. The free model parameters 𝑀 c and 𝜇 thereby define the scale and the abruptness of the transition when the slope of the profile goes from 3 to 0 for decreasing halo masses.  \nThe total fraction of stars, 𝑓 star , and the fraction of stars that belong to the central galaxy, 𝑓 cga , indirectly affect the available gas that can be pushed out by feedback processes. They are parametrised as  \n𝑓 𝑖 ( 𝑀 ) = 0 . 055 GLYPH<18> 𝑀 𝑀 s GLYPH<19> 𝜂 𝑖 , (11)  \nwith 𝑖 = { star , cga } and where 𝑀 s = 2 . 5 × 10 11 M ⊙ /h. The power  \nlaw is constructed to match the Moster relation (Moster et al. 2013), which is reasonably well-understood for galaxy groups and clusters. Note that following Giri & Schneider (2021), we redefine the parameters as 𝜂 ≡ 𝜂 star and 𝛿 𝜂 ≡ 𝜂 cga -𝜂 , star . These two additional bfc parameters, together with the five gas parameters, are summarised in Table 1.  \nAnefficient cosmological analysis will marginalise over a minimumnumberofbaryonic feedback parameters. How many parameters are required is an open research question that ultimately depends on the unknown baryonic feedback realised in nature. For now, we adhere to the requirement that a given parametrisation needs to be able to fit the matter power spectrum suppression predicted by a range of hydrodynamical simulations. Following Giri & Schneider (2021), in this work we consider the models BCEmu7, BCEmu3, and BCEmu1, referring to the number of model parameters varied in the analysis. While in BCEmu7 all parameters introduced above are kept free, BCEmu3 only allows log 10 𝑀 c , 𝜃 ej and 𝜂 𝛿 to vary and BCEmu1 only uses log 10 𝑀 c as a free parameter. The fixed parameters in the BCEmu3 and BCEmu1 models are listed in Table 1 and have been selected so that they provide the best fit to a variety of hydrodynamical simulations (see Giri & Schneider 2021, for more information about the method). Both BCEmu7 and BCEmu3 are able to reproduce the baryonic suppression or the power spectrum predicted by a variety of hydrodynamical simulations to better than one percent. The BCEmu1 model, on the other hand, shows deviations of order five per cent, hinting towards the possibility that one parameter is generally insufficient to describe the variety of existing results from simulations (Giri & Schneider 2021).  \nA key feature of the baryonification model is that it is based on physically motivated profiles around halo centres. The model is not restricted to the power spectrum but can also be used to obtain the three-dimensional baryonified density field and therefore many corresponding summary statistics. In the following we will take advantage of the connection between power spectrum and halo profiles to obtain simultaneous predictions for the cosmic shear and the stacked kinetic Sunyaev-Zel'dovich signal. Compared to other approaches, the baryonification model BCEmu7 is not calibrated to specific hydrodynamical simulations. It rather depends on the parameters of empirical density profiles that are broadly motivated by observations. The model therefore provides a independent check with very different modelling choices and systematics compared to the subgrid modelling in hydrodynamical simulations. We note that despite the model's valuable flexibility, it can result in potentially non-physical scenarios. Other limitations include the fact that the model parameters are currently assumed to have no redshift dependence, as well as the fact that the gas profiles do not separate a hot and cold gas component.", "3.3 Kinetic Sunyaev Zel'dovich signal": "ThekSZmeasurementscanbeusedtoconstrainthegasdensity. This effect arises from the bulk motion of the ionised gas in and around galaxies, galaxy groups and clusters, which imparts a Doppler shift on CMB photons. It preserves the blackbody frequency spectrum of the CMB and instead the thermodynamic temperature fluctuates as  \nΔ 𝑇 kSZ 𝑇 CMB = 𝜎 T 𝑐 ∫ los 𝑒 -𝜏 𝑛 e 𝑣 p 𝑑𝑙 , (12)  \nwhere 𝑇 CMB is the present-day temperature of the CMB, 𝜎 T is the Thomson cross-section, 𝑐 is the speed of light, 𝑛 𝑒 is the freeelectron physical number density, 𝑣 p is the peculiar velocity and  \n𝜏 is the optical depth due to Thomson scattering along the lineof-sight distance, 𝑑𝑙 . Following Schaan et al. (2021) and Amodeo et al. (2021), for the redshift range of the kSZ measurements used in this work, the mean optical depth is observed to be below the percent level (Planck Collaboration 2016) and the CMASS galaxy groups are optically thin, therefore we can assume that the integral 𝑒 -𝜏 ≈ 1. Furthermore, as the measurements are stacked, the velocity field is independently estimated from the large-scale distribution of galaxies via a reconstruction method, thereby eliminating the dependence on the velocity. The peculiar velocity, 𝑣 p , projected along the line-of-sight can be written as the root mean square (RMS) of the peculiar velocity, 𝑣 r , so that the resulting shift in the CMB temperature can be approximated as  \nΔ 𝑇 kSZ 𝑇 CMB = 𝜏 gal ( 𝜃 ) 𝑣 r 𝑐 , (13)  \nwhere 𝜏 gal refers to the contribution of the optical depth to Thomson scattering of the galaxy group considered. For the median redshift of the CMASS sample in the linear approximation, 𝑧 = 0 . 55, the RMS of the peculiar velocities projected along the line of sight is 𝑣 𝑟 = 1 . 06 × 10 -3 𝑐 (Schaan et al. 2021). The uncertainty on the velocity reconstruction is estimated to be less than a few percent, which we ignore given the statistical precision (Schaan et al. 2021; Ried Guachalla et al. 2023; Hadzhiyska et al. 2023).  \nTo model the ACT kSZ measurements it is necessary to convolve equation 13 with the beam profiles utilised at the f90 and f150 frequencies. We follow Schaan et al. (2021) and approximate the beam using a Gaussian with FWHM = 2 . 1 arcmin for the former band and FWHM = 1 . 3 arcmin for the latter. Furthermore, to minimise noise due to degree-scale CMB fluctuations, compensated aperture photometry filters were also used on the observations. We therefore apply the same filter function that was used in the analysis of the data. The smoothing function is + 1 between 𝜃 < 𝜃 d , -1 between 𝜃 d < 𝜃 < √ 2 𝜃 d and 0 otherwise, where 𝜃 d is the aperture radius centred around each galaxy (Schaan et al. 2021, equation 9).  \nTo calculate 𝜏 gal measured within a disk of radius 𝜃 centred on the group or cluster, we assume spherical symmetry and integrate the electron number density, 𝑛 e , over the line-of-sight to a cut-off radius of 𝑅 200 as  \n𝜏 gal ( 𝜃 ) = 2 𝜎 𝑇 ∫ 𝑅 200 0 𝑛 e ( √︃ 𝑙 2 + 𝑑 A ( 𝑧 ) 2 𝜃 2 ) 𝑑𝑙 , (14)  \nwhere 𝑑 A is the angular diameter distance. We assume a fully ionised mediumwithprimordialabundancestodescribetheelectron density in terms of the gas density as  \n𝑛 e ( 𝑟 ) = ( 𝑋 H + 1 ) 2 𝜌 gas ( 𝑟 ) 𝑚 amu , (15)  \nwith 𝑋 H = 0 . 76 being the hydrogen mass fraction and 𝑚 amu the atomic mass unit.  \nIn order to calibrate the model for the kSZ profile, the mean halo mass of the CMASS galaxy sample is needed. Note that because the integrated kSZ signal scales with the gas mass and this quantity approximately tracks the halo mass, it is important that the theoretical predictions are for a sample with the same mean halo mass as the CMASS sample. Given the significant scatter in the literature, we choose to include an additional model parameter in the analysis, 𝑀 h , 200 , corresponding to the mean 𝑀 200 of the CMASS sample, with a prior range provided in Table 1. The justification for this prior choice, and an investigation of its impact are given in Appendix B3.  \nFigure 2. Top: The impact of varying parameters of the baryonification model, BCEmu, on the ratio of the total matter power spectrum compared to a dark-matter only power spectrum, 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) . Each panel varies one baryonification parameter at a time within the prior bounds (reported in Table 2), whilst keeping the remaining six parameters at their fiducial fixed value, corresponding to the fit to a range of hydrodynamical simulations (Giri & Schneider 2021). For reference, various predictions for the suppression of the matter power spectrum from simulations are over-plotted in black: FLAMINGO (Schaye et al. 2023, solid line); BAHAMAS (McCarthy et al. 2017, dash-dotted line); SIMBA (Dav'e et al. 2019, dotted line), MillenniumTNG (Pakmor et al. 2023, dashed line) and FABLE (doubledotted-dashed line, Henden et al. 2018, Bigwood et al. in prep.). Bottom: The stacked kSZ temperature profile at 98GHz as a function of angular radius, 𝜃 , centred on the group or cluster (bottom) for the same baryonification parameters as above. The ACT CMASS measurements at 98GHz are shown as the black data points in the bottom panels and the model profiles are convolved with the f90 beam profile for comparison.  \n<!-- image -->", '3.4 Modelling kSZ with baryonification': "The baryonification method provides a model for the gas density (equation 9) and describes the gas content as measured by kSZ (equations 13-15). We select this model to jointly analyse the kSZ measurements with the lensing data given that it is agnostic to any choice of hydrodynamical simulation.  \nIn Fig. 2, we aim to build a better intuition on the BCEmu model by exploring the impact of each parameter on the suppression of the matter power spectrum, 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) (top row) and the kSZ radial temperature profile 𝑇 kSZ ( 𝜃 ) (bottom row). At a fixed cosmology 4 , we test the dependence on each of the seven BCEmu parameters from left to right within the defined prior range, whilst fixing the remaining six model parameters at their fiducial fixed value (see Table 1).  \nThe 𝑀 c parameter controls the mass or proportion of groups and clusters that have had feedback-induced gas removal and has been previously identified as the most important in the model (Schneider & Teyssier 2015; Giri & Schneider 2021). Indeed, we see that within its prior range, this parameter modulates the amplitude, slope and extent of the suppression of the matter power spectrum at scales 𝑘 ≳ 0 . 1 and the amplitude and slope of the kSZ profile. For clusters of mass greater than 𝑀 c , the slope of the gas profile (equation 10) tends to 𝛽 = 3 and the gas profile approaches a truncated NFW profile. However, for groups of mass smaller than 𝑀 c , the slope of the power law decreases and the gas profile more closely resembles one that has experienced AGN feedback and had gas ejected from the halo. Hence, a larger value of 𝑀 c results in a greater proportion of groups and clusters having gas profiles that mimic the effect of baryonic feedback, leading to a greater suppression of the matter power spectrum and simultaneously, a decrease in the amount of gas within an aperture centred on the galaxy, which gives a smaller kSZ signal.  \nIn the reduced complexity BCEmu1 model, the remaining pa-  \nrameters are kept fixed, although it is clear that their value choice can have a significant impact on the predictions for the matter power spectrum and the kSZ signal. We observe that even at fixed 𝑀 c , a larger value of 𝜃 ej increases the radius that gas is 'ejected to' in the gas profile and effectively causes matter to be redistributed on smaller 𝑘 -scales. This corresponds to a decrease in the amount of gas as measured by the amplitude of the kSZ signal. Similarly, we find that increasing the exponent in the gas profile parameterisation 𝛾 results in a greater suppression of the power spectrum and lower kSZ amplitude, although the impact is on smaller scales ( 𝑘 ≳ 1) than the effect of log 10 𝑀 c and 𝜃 ej . The remaining gas parameters, 𝛿 and 𝜇 , have the reverse effect, i.e., decreasing them results in an enhanced power spectrum suppression and a lower kSZ signal.  \nDecreasing the values of stellar parameters, 𝜂 and 𝜂 𝛿 , makes a less steep stellar-halo mass relation, so that a larger amount of stars condense out of the gas and form galaxies. This causes a boost in the matter power spectrum at large 𝑘 , as we observe in Fig. 2. The kSZ radial temperature profile is only dependent on the gas density profile of groups and clusters (equation 15) and so not directly impacted by the stellar parameters 𝜂 and 𝜂 𝛿 . However, 𝜂 can indirectly have a small impact on the kSZ signal as it alters the number of stars in groups and clusters, and therefore the the reservoir of available gas.  \nWe note two limitations in our modelling of the kSZ signal. While we vary the halo mass of the kSZ sample in our analysis, we assume that the mass dependence of the model is valid beyond the range probed by the kSZ measurement and at varying redshift. The kSZ measurements span a mass range 𝑀 200 ≈ [ 0 . 5 -7 ] × 10 13 𝑀 ⊙ , which is similar to the halo mass range that cosmic shear is most sensitive to ( 𝑀 200 ≈ 10 13 . 5 𝑀 ⊙ ). Nevertheless, cosmic shear still has contributions from higher and lower mass halos (e.g. Salcido et al. 2023).  \nFuture kSZ measurements that span multiple mass and redshift bins will provide a better understanding of the dependence of the kSZ signal on galaxy properties and the suitability of this assumption. Another issue that we do not address in the present study is the  \nTable 2. Summary of cosmological, observational and astrophysical parameters and priors used in the analysis. In the case of flat priors, the prior is bound to the range indicated in the 'value' column, while Gaussian priors are described by their mean and 1 𝜎 width.  \nrole of centrals vs. satellites in the observed stacked kSZ profiles. Our theoretical predictions correspond to central galaxies that are assumed to be perfectly centred within their host haloes and we may reasonably expect some degree of bias (with respect to theoretical predictions) to be introduced by the inclusion of satellites. Without forward modelling the BOSS CMASS selection function, which is beyond the scope of this work, it is difficult to predict the magnitude or sign of this effect. On the one hand, satellites will obviously be mis-centred with respect to their host haloes and one may expect this to lead to a lower kSZ signal. On the other hand, a stellar mass-based selection implies that the satellites will typically be in hosts that are more massive than a host which has a central of similar stellar mass. This will tend to boost the kSZ signal. For the present study we neglect these uncertainties, leaving their careful consideration for a future study.", '4 MODEL PIPELINE AND VALIDATION': 'In this section, we briefly describe the set up of the cosmological inference pipeline (Section 4.1). We validate the robustness of this pipeline using each of the four baryon models with a mock analysis, described in Appendix A. Here, we briefly motivate the choices made in the construction of the mock data (Section 4.2) and summarise the findings (Section 4.3).', '4.1 Inference pipeline': 'To analyse the cosmic shear data we build upon the public DES Y3 cosmological inference pipeline. We utilise one parameter for redshift and shear calibration per tomographic bin, with prior ranges set to those used in DES Y3. Without feedback models, we have 15 parameters that we marginalise over in the analysis. Table 2 summarises the cosmological, observational and astrophysical priors used in this work.  \nParameters are estimated via nesting sampling using M/u.pc/l.pc/t.pc/i.pc/n.pc/e.pc/s.pc/t.pc 5 (Feroz et al. 2009) within the COSMOSIS framework (Zuntz et al. 2015), with the sampler settings listed in Appendix D. We note however that DES & KiDS Collaborations (2023) demonstrates that M/u.pc/l.pc/t.pc/i.pc/n.pc/e.pc/s.pc/t.pc can underestimate the 68% confidence levels for 𝑆 8 , at the level of ∼ 10 -20%, while Polychord (Handley et al. 2015) is more accurate. In agreement with this finding, Appendix D reports that sampling with Multinest instead of Polychord in a WL-only analysis with BCEmu7 results in a 18% smaller 68% confidence level for 𝑆 8 , and a 9% smaller 68% confidence level for 𝑆 8 in a joint WL and kSZ analysis with BCEmu7. We follow Amon et al. (2022) and Secco, Samuroff et al. (2022) when reporting the parameter constraints and quote the mean of the 1D marginal distribution, along with the 68% confidence limit, which defines the area around the peak of the posterior within which 68% of the probability lies 6 .', '4.2 Mock data': 'In order to assess the robustness of our lensing inference pipeline, non-linear power spectrum model and baryon feedback models, we perform analyses using synthetic data and test the ability to recover unbiased cosmology. The mock data was created using the best fit cosmological parameters obtained from the DES Y3 joint lensing and clustering analysis (DES Collaboration 2022) 7 .  \nWe create two dark matter-only mock data vectors, using different models for the non-linear matter power spectrum. The first uses HM20, the same model used to analyse the data throughout this work. This mock is important for testing that the analysis pipeline can accurately recover cosmological parameters before considering feedback effects. The second mock uses the EuclidEmulator2 (EE2), which has been shown to be accurate to 1% for 𝑘 ≤ 10 ℎ Mpc -1 (Euclid Collaboration: Knabenhans et al. 2019; Knabenhans et al. 2021; Adamek et al. 2023). This mock is used to ensure that HM20 provides a sufficiently accurate description for the non-linear matter power spectrum when compared to EE2, for the full angular scale range of the DES Y3 cosmic shear data.  \nTo test the four baryon mitigation strategies, we chose three baryon feedback scenarios to build the mock data and their predictions for the suppression of the matter power spectrum are shown as the red lines in Fig. 1. These choices include extreme scenarios as they are designed to test the flexibility and limits of the baryon models, rather than an attempt to chose the most accurate prescription. First, we consider the upper limit of the BAHAMAS hydrodynamic simulation suite, with 𝑇 heat = 8 . 0 (McCarthy et al. 2017). Next, we consider a mock universe with a more extreme power spectrum prescribed by cosmo-OWLS 𝑇 heat = 8 . 5 (Le Brun et al. 2014), although this simulation does not replicate the local gas fractions in groups and clusters. As it is possible that hydrodynamical simulations do not capture the complexities of feedback, we want to test the ability of the baryon models to accurately capture a scenario that modulates the matter power spectrum with a different shape than that in typical simulations. We consider a mock with a suppression described by the 𝐴 mod parameter (Preston et al. 2023; Amon &  \nL. Bigwood, A. Amon et al.  \nEfstathiou 2022), where 𝐴 mod ≈ 0 . 82 is that required to reconcile the DES Y3 cosmic shear with Planck Λ CDMcosmology.', '4.3 Mock results': "The results of these mock tests are detailed in Appendix A. Here, we summarise the findings, which are shown in Fig. A2:  \n- · When analysing a EE2-generated dark-matter only mock with the HM20darkmatter-only model, 𝑆 8 is over-estimated by ∼ 0.4 𝜎 . (This is reduced to ∼ 0.2 𝜎 when analysing restricted angular scale measurements.) Furthermore, we find that Ω m is underestimated by ∼ 0.9 𝜎 . However, we identify projection effects in this parameter of ∼ 0.5 𝜎 by analysing an HM20 mock with an HM20 model. In this work, we focus on the 𝑆 8 parameter, and we note that further testing is needed to assess the reliability of the Ω m constraints.\n- · The DES Y3 ' Λ CDM-optimised' scale cuts underestimate 𝑆 8 in all three mock baryonic feedback scenarios. This is as expected, as these scale cuts were defined with the OWLS-AGN scenario, which predicts less power suppression than BAHAMAS 8.0, cosmoOWLS 8.5 and 𝐴 mod . While the scale cuts, by design, remove the sensitivity of the analysis to the impact of baryon feedback, this method's success relies on the true feedback scenario to be less extreme than the simulation used to define the cuts.\n- · HM20, as used with their fiducial BAHAMAS-based prior, underestimates 𝑆 8 by ∼ 0.7 𝜎 for the BAHAMAS 8.0 mock, and by more than 1 𝜎 when analyzing a mock with a more extreme feedback scenario. Using a wide prior alleviates this, with the model recovering the true cosmology within ∼ 0.5 𝜎 for all mock scenarios, although with a cost of almost a factor of two in the precision of the 𝑆 8 constraint.\n- · The SP(k) model, used with both a wide and X-ray prior, can recover the input cosmology to within ∼ 0.2 𝜎 for a BAHAMAS 8.0 mock. With the more extreme cosmo-OWLS 8.5 and 𝐴 mod mocks, this model underestimates the value of 𝑆 8 by up to ∼ 0.5 𝜎 . This is as expected, as both cosmo-OWLS 8.5 and 𝐴 mod are outside the expected baryon fractions as compared to observations by Akino et al. (2022).\n- · When allowing only one parameter to vary in BCEmu emulator, BCEmu1, we recover 𝑆 8 to within ∼ 0.2 𝜎 for the BAHAMAS 8.0 mock, and ∼ 0.5 𝜎 for cosmo-OWLS 8.5. When we use the more flexible BCEmu7, all mock scenarios recover the true cosmology within ∼ 0.5 𝜎 , and the error bar on 𝑆 8 is up to 1.5 times wider.  \nOverall, we find that when more restrictive modelling choices are used, we tend to underestimate 𝑆 8 . The bias is worsened when restrictive choices are used to analyse mocks with the more extreme baryonic feedback scenarios, i.e. cosmo-OWLS 8.5 and 𝐴 mod . We note that despite the greater accuracy of using more conservative priors or marginalising over a greater number of baryonic feedback parameters, it is at the expense of the precision. We find that the uncertainty on the constraint of 𝑆 8 can degrade by up to a factor of two when switching to more flexible modelling choices.", '5 RESULTS: ASSESSING MODELS FOR BARYONIC EFFECTS': 'The results of the WL-only DES Year 3 analysis are divided into three sections. In Section 5.1, we present the headline cosmological constraints using the four baryon feedback strategies outlined in Section 3.2. In Section 5.2, we show the constraints on the suppression of the power spectrum. Finally, in 5.3, we explore the dependence of  \nFigure 3. Summaryof the 1D marginalised constraints on 𝑆 8 from analysing the DES Y3 cosmic shear data with different baryonic feedback mitigation strategies. The mean of the 𝑆 8 marginalized posterior is indicated by the symbol and 68% confidence levels are shown as horizontal bars. The primary result from the joint analysis of WL + kSZ using BCEmu7 is represented as the purple shaded region. We compare to the Planck TTTEEE result presented in Efstathiou & Gratton (2021).  \n<!-- image -->  \nour results on the model complexity and prior choices within each strategy.', '5.1 Cosmological parameter constraints': "The 1D marginalised constraints obtained for 𝑆 8 are summarised in Fig. 3 for analyses using all model variants. Here, we compare those from the DES Y3 ' Λ CDM optimised' scale cuts, HM20: BAHAMAS, SP(k): Xray and BCEmu7. For these approaches, the mean marginal values of 𝑆 8 are found with 68% confidence levels to be  \nBCEmu7 : 𝑆 8 = 0 . 818 + 0 . 017 -0 . 024 Spk : Xray : 𝑆 8 = 0 . 806 + 0 . 015 -0 . 013 HM20 : BAHAMAS : 𝑆 8 = 0 . 811 + 0 . 013 -0 . 014 Scale cuts : 𝑆 8 = 0 . 805 + 0 . 018 -0 . 018 . (16)  \nFor reference, we show the result when all angular scales are analysed without any model for baryonic effects and the Planck TTTEEE 8 Λ CDMresult (orange, Efstathiou & Gratton 2021). The 2D marginalised posteriors for 𝑆 8 , Ω m and 𝜎 8 using DES Y3 weak lensing data are plotted in the left panel of Fig. 4, also showing no baryon mitigation (dotted grey), DES Y3 ' Λ CDMoptimised' scale cuts (dashed green), HM20: BAHAMAS (navy), SP(k): Xray (solid green), BCEmu7 (pink) and Planck TTTEEE (orange). Table 3 lists the mean constraints on 𝑆 8 and Ω m and quantifies the goodness of fit to the data of each modelling variant by quoting the minimum and reduced 𝜒 2 , 𝜒 2 min , 𝜒 2 red . We find that each baryon feedback analysis variant demonstrates a suitable fit to the measurements.  \nWhen analysing all angular scales without modelling baryonic feedback, we attain a low value of 𝑆 8 = 0 . 794 + 0 . 013 0 . 011 . This is up  \nFigure 4. Left: The marginalised posteriors for Ω m , 𝜎 8 and 𝑆 8 using the DES Y3 cosmic shear data and different baryon feedback models. We compare the DES Y3 optimised scale cut approach (green dashed) to scenarios where all angular scales of the DES Y3 lensing data are modelled: HM20, using their recommended BAHAMAS-based prior (HM20: BAHAMAS, navy), the seven-parameter baryonification model (BCEmu7, pink) and the SP(k) emulator, using their 'X-ray observational prior' (SPk: Xray, solid green). For reference, we show the case where all scales of the DES data vector are used but baryons are not modelled (gray dotted) and the Planck TTTEEE likelihood (orange, Efstathiou & Gratton 2021). The inner and outer contours show the 68% and 95% confidence levels, respectively. Right: The corresponding constraints on the suppression of the total matter power spectrum compared to a dark matter-only scenario, 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) , at 𝑧 = 0 using the DES Y3 cosmic shear and again, HM20 (HM20: BAHAMAS, navy), BCEmu (BCEmu7, pink) and SP(k) (SPk: Xray, green). The solid lines show the mean suppression and the shaded regions indicate the 68% confidence levels. For reference, various predictions from hydrodynamical simulations are over-plotted in black: FLAMINGO (Schaye et al. 2023, solid line); BAHAMAS (McCarthy et al. 2017, dash-dotted line); SIMBA (Dav'e et al. 2019, dotted line), MillenniumTNG (Pakmor et al. 2023, dashed line) and FABLE (doubledot-dashed line, Henden et al. 2018, Bigwood et al. in prep.).  \n<!-- image -->  \nTable 3. The constraints on 𝑆 8 , Ω m attained for each method of baryonic feedback mitigation. We report the mean value of each parameter with errors given by the 68% confidence levels. We also demonstrate the quality of the fit by reporting 𝜒 2 min , the minimum value of 𝜒 2 , for each analysis variant. We report 𝜒 2 red = 𝜒 2 min / 𝑁 dof , with 𝑁 dof being the number of degrees of freedom 𝑁 dof = 𝑁 dp -𝑁 param , where 𝑁 dp and 𝑁 param are the number of data points and model parameters utilised in the analysis. With the exception of the scale cuts method, angular scales down to 2.5 arcmin of the DES Y3 data were used.  \nto 1 𝜎 lower than constraints attained with modelled baryonic effects, highlighting the importance of mitigating feedback to avoid biased cosmology. When accounting for baryonic effects, we find, in agreement with the mock analysis for the models tested, that the  \nmeasured 𝑆 8 is consistent at the level of 0.6 𝜎 ( ∼ 2%) 9 . However, the errorbar on the 𝑆 8 constraint varies by a factor of 1.5. In more detail,  \nwe see that HM20: BAHAMAS and SP(k):Xray give the tightest constraints on the 𝑆 8 parameter, which are in excellent agreement with each other. This is expected, as these models are calibrated on hydrodynamical simulations informed by X-ray constraints. The BCEmu7analysis gives the highest value of 𝑆 8 . As in the case of the mock analysis, this supports our findings that restrictive modelling choices for baryonic feedback leads to lower value of 𝑆 8 when compared to more flexible models. The flexibility of BCEmu7 comes at a cost, as the uncertainty on 𝑆 8 is a factor of 1.5 larger than that in the HM20: BAHAMAS analysis. This is to be expected given the degeneracy between the extremity of feedback and 𝑆 8 ; greater flexibility in the modelling of baryons inevitably results in a larger errorbar on the 𝑆 8 constraint. The analysis using the DES Y3 scale cuts does not suffer as substantial a loss in constraining power as the BCEmu7 case, but results in the lowest value for 𝑆 8 10 .", '5.2 Power suppression constraints': "Baryon feedback processes modify the gravitational evolution of the cosmic density field and suppress the matter power spectrum compared to a dark matter-only scenario on non-linear scales, as seen in hydrodynamic simulations. This effect has been previously observed by analyses of weak lensing data using variations of the baryonification model (Schneider et al. 2022; Chen et al. 2023; Aric'o et al. 2023). In this section, we constrain the amplitude and scale dependence of the suppression of the matter power spectrum due to baryonic effects, 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) , using the DES Y3 cosmic shear. For the first time, we show the model-dependence of the constraints by considering the model complexity of the baryonification model and the comparison to the SP(k) and HM20 models. For each baryonic feedback model, we record the power spectrum suppression at each step in the chain.  \nIn the right panel of Fig. 4 we plot the mean suppression and the and 68% confidence levels inferred from analyses with our three baseline models: BCEmu7, HM20: BAHAMAS ( Θ AGN = 7 . 6 -8 . 0) and SP(k):Xray. We find that the suppression inferred by the three models are consistent within the 68% confidence limits up to 𝑘 ≈ 3 ℎ /Mpc. However, we find that BCEmu7 allows more extreme suppression of the power spectrum at all non-linear scales. There are substantial differences in the size of the uncertainties, correlated with the flexibility of the model. BCEmu7, the most flexible model, has the largest uncertainty. SP(k):Xray provides the tightest constraints on the power spectrum suppression and constrains a less extreme feedback scenario in terms of the amplitude and the scale extent of the suppression.  \nWe compare our constraints to five hydrodynamical simulations: FLAMINGO, (Schaye et al. 2023, solid line); BAHAMAS (McCarthy et al. 2017, dash-dotted line); SIMBA (Dav'e et al. 2019, dotted line), MillenniumTNG (Pakmor et al. 2023, dashed line) and  \n( 𝜎 2 𝑆 8 ) 2 ] 1 / 2 , where Δ 𝑆 8 is the difference between the respective mean values and 𝜎 1 𝑆 8 , 𝜎 2 𝑆 8 are the 1 -𝜎 errors on 𝑆 8 .  \n10 Note that the DES Y3 ' Λ CDMoptimised' cosmic shear analysis (Amon et al. 2023; Secco, Samuroff et al. 2022) obtains a lower value of 𝑆 8 than that obtained here ( 𝑆 8 = 0 . 772 + 0 . 018 -0 . 017 ). Based on the study of the impact of analysis choices in (DES & KiDS Collaborations 2023), we attribute the difference in our results primarily to the use of HM20 to model the dark matter non-linear matter power spectrum, which was shown to be more accurate than Halofit, as well as the intrinsic alignment model that we chose for this analysis, and the choice to fix the neutrino mass in the analysis. In Appendix C we investigate the impact of these choices further.  \nFABLE (doubledot-dashed line, Henden et al. 2018, Bigwood et al. in prep.). HM20 predicts a feedback strength that encompasses BAHAMAS 𝑇 heat = 7 . 8, which is unsurprising, given that the model is calibrated to span these simulations, but it is notable that the mean constraint is more extreme on all scales. Interestingly, all three models find the mean suppression to be more extreme than FLAMINGO on scales 𝑘 ∼ 0.2-4 ℎ -1 Mpc, at the level of 1 . 1 𝜎 (BCEmu7), 3 . 0 𝜎 (SP(k):Xray) and 1 . 4 𝜎 (HM20:BAHAMAS) at 𝑘 = 2 ℎ -1 Mpc.  \nThe power suppression constraints are broadly consistent with those from previous weak lensing analyses (Schneider et al. 2022) and slightly more extreme than the constraints of Chen et al. (2023); Aric'o et al. (2023); Garc'ıa-Garc'ıa et al. (2024); Terasawa et al. (2024). Here, we make note of some details. Owing to the enhanced statistical power of the DES Y3 data over that of the Kilo-Degree Survey, we find substantially improved constraints from the WLonly analysis compared to Schneider et al. (2022), even though we include an additional parameter for the intrinsic alignment model and eight additional nuisance parameters to account for the uncertainty in the shear and redshift calibration. Furthermore, we use a different baryonification model to that of Chen et al. (2023) and Aric'o et al. (2023) and we do not impose X-ray priors on any of the baryon parameters. One way in which this model is different from BCEmu is that only particles within the virial radius, 𝑅 200 , are displaced (see Grandis et al. 2024, for a more detailed discussion of the model comparison). One implication of the model differences is that it is not straightforward to compare constraints on the 𝑀 c parameter, though we note that in our BCEmu scenario, the posterior on the 𝑀 c parameter is not limited by the upper value of the prior.", '5.3 Impact of model complexity': "Themockanalysisrevealed that marginalising over a greater number of baryonic nuisance parameters, or utilising wider priors on these parameters, generally improved the accuracy of the cosmological constraints. However, we saw that this was at the expense of inflated errors on the cosmological parameters, which is clearly sub-optimal for an effective cosmological analysis. In this section, we explore the impact of altering the complexity and prior choices of each baryon feedback model on the measured cosmological and baryonic constraints when analysing the DES Y3 𝜉 ± measurements.  \nThe upper left panel of Fig. 5 shows the impact of the BAHAMAS-informed prior on the HM20 feedback parameter Θ AGN on the marginalised 𝑆 8 posterior. In the light blue constraint, we extend the prior range for Θ AGN outside of the calibration range to encompass more extreme feedback scenarios, as well as a darkmatter only scenario. The parameters are degenerate and opening up the Θ AGN prior leads to long tails that extend to higher values of 𝑆 8 , such that the mean constraint is > 0.5 𝜎 higher. This illustrates how high values of 𝑆 8 are disfavoured by the restricted prior on the baryonic feedback model, and suggests that weak lensing data may favour a higher value of Θ AGN than the BAHAMAS simulations span. (Although we note that the HM20 model was only calibrated within the BAHAMAS range, so the mapping between the power spectrum suppression and the baryon fraction outside this range is uncertain.) Similarly, for the case shown in the right-hand panel of Fig 5, when the X-ray informed prior is lifted, the data constrains higher values of 𝑆 8 by ∼ 0.5 𝜎 (the posteriors on the SP(k) parameters are shown in Appendix B1). These shifts are consistent with those we determine in the mock analysis (Section 4.2, Appendix A). In Section 7, we discuss the implications of these results on our understanding of the gas models and observations.  \nWe test the impact of limiting the baryonification model com-  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFigure 5. Top: The degeneracy of the marginalised 𝑆 8 posterior with varying the baryon model complexity or prior choices on baryon model parameters. For each panel the inner and outer contours show the 68% and 95% confidence levels, respectively. Left: HM20 𝑆 8 and Θ AGN attained when using a prior bracketing the BAHAMAS simulations (dark blue, 𝜃 AGN = 7 . 6 -8 . 0) and a less-informative prior choice, (light blue, 𝜃 AGN = 7 . 3 -9 . 0). Centre: BCEMu 𝑆 8 and log 10 ( 𝑀 c ) attained for the full seven parameter (pink, BCEmu7), three parameter (purple, BCemu3) and one parameter (blue, BCEmu1) models. Right: SP(k)𝑆 8 and 𝛽 attained for the X-ray informed model (dark green) and a less-informative prior (light green). Bottom: The corresponding constraints on the suppression of the total matter power spectrum compared to a dark matter-only scenario, 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) , at 𝑧 = 0 for each of the models. In general, when allowing for a more flexible model, the constraints indicate more extreme suppression of power, although with degraded constraining power. For each panel the solid lines show the mean suppression predicted, and the shaded regions the 68% confidence levels. For reference, various predictions for the suppression of the matter power spectrum from simulations are over-plotted in black: FLAMINGO (Schaye et al. 2023, solid line); BAHAMAS (McCarthy et al. 2017, dash-dotted line); SIMBA (Dav'e et al. 2019, dotted line), MillenniumTNG (Pakmor et al. 2023, dashed line) and FABLE (doubledot-dashed line, Henden et al. 2018, Bigwood et al. in prep.).  \n<!-- image -->  \nplexity to the one (BCEmu1) and three-parameter (BCEmu3) case, compared to the fiducial BCEmu7. The central panel of Fig. 5 shows the marginalised posteriors on Ω m , 𝑆 8 and the baryonification parameter log 10 𝑀 c . The posteriors on all of the baryonic feedback parameters are shown in Appendix B2. BCEmu7 and BCEmu3 produce comparable constraints on cosmological and feedback parameters, generating a shift in 𝑆 8 of only ∼ 0.2 𝜎 . We also do not see any significant improvement of the precision on these constraints when marginalising over four less baryonic parameters. However, switching to BCEmu1 from BCEmu7 results in a substantially lower values of 𝑆 8 by > 0.5 𝜎 , consistent with the mock analysis. This reduction in the flexibility of the model forces log 10 𝑀 c to a larger value, since the full extremity of feedback has to be captured by only the one parameter. We attribute this shift in the cosmology and feedback parameters to the values at which the remaining six 'underthe-hood' baryonification parameters are fixed to in the BCEmu1 model, quoted in Table 1 11 . Given that the analysis using BCEmu1  \ngives a larger value of log 10 𝑀 c than that attained using BCEmu7, and that for example, the posterior on 𝜃 ej is toward higher values than the fixed 3.5, this implies that the values for the simulationinformed fixed parameters of BCEmu1 represent a less-extreme feedback scenario than those constrained by an analysis of DES cosmic shear.  \nThe lower panels of Fig. 5 illustrate how the corresponding constraints on the power spectrum suppression are sensitive to the restrictiveness of the choices within each model. Analogously to the upper panels, the left and right-most panels test switching to the more conservative prior choices of the HM20 and SP(k) models, i.e. Θ AGN = 7 . 3 -9 . 0 and a prior spanning the feedback landscape of the ANTILLES suite, respectively. The central panel shows the impact of restricting the model complexity of BCEmu, by allowing only one (BCEmu1) out of the full seven (BCEmu7) baryonification parameters to vary in the analysis. As before, we compare the data constraints to predictions from hydrodynamical simulations.  \nsimulations, then summing the likelihoods to find the best-fit parameters to all of the simulations (Giri & Schneider 2021).", '14 L. Bigwood, A. Amon et al.': 'For each model, when we opt for the more flexible modelling choices, the mean constraint on the power spectrum suppression tends to more extreme scenarios at all non-linear 𝑘 with respect to their more restrictive counterpart. In particular, each of the flexible models is most consistent with more extreme feedback scenarios (e.g., SIMBA). This could suggest that the restrictive baryonic modelling choices do not have the flexibility to capture the full extremity of feedback that the data prefers and therefore the higher values of 𝑆 8 . It is clear, however, that the weak lensing data cannot place strong constraints on feedback on its own. The use of the more conservative model choices comes at the expense of reduced constraining power. For example, SP(k) displays over a factor of two increase in the size of the 68% confidence level at all scales when switching from the restrictive to conservative modelling choice. A complementary avenue to simulation-based models that still maximises cosmological constraints is to use the flexible model framework and jointly analyse the lensing with observations of the gas.', '6 RESULTS: JOINT WEAK LENSING + KINETIC SZ': 'In the previous section we demonstrated that weak lensing 𝑆 8 constraints are degenerate with the amount of baryonic feedback. A complementary approach to simulation-informed baryon models is to use a flexible model to jointly analyse cosmic shear with probes of the gas distribution in order to better constrain the model parameters. In this section, we report the results of a joint analysis of the DES Y3 cosmic shear and ACT kSZ measurements, using the BCEmu7 baryon mitigation model, described in Section 3.3.', '6.1 Cosmological parameter constraints': "The left panel of Fig. 6 compares the marginalised posteriors on 𝑆 8 and Ω m attained from the WL-only analysis using BCEmu7 (pink) to those obtained from a joint analysis with kSZ (purple). The Planck TTTEEE Λ CDM posteriors are shown for reference (Efstathiou & Gratton 2021). The mean marginal value of 𝑆 8 for the joint analysis is found with 68% credible levels to be  \nLensing + kSZ 𝑆 8 = 0 . 823 + 0 . 019 -0 . 020 , (17)  \nwhich corresponds to a ∼ 0.2 𝜎 shift towards higher values with respect to the result of the WL-only BCEmu7 analysis of 𝑆 8 = 0 . 818 + 0 . 017 -0 . 024 .  \nAmon & Efstathiou (2022) and Preston et al. (2023) have proposed that the 𝑆 8 tension could be resolved if the non-linear matter power spectrum is suppressed more strongly than is currently assumed in weak lensing analyses, either due to unmodelled baryonic feedback effects or non-standard dark matter. For baryonic feedback to be the source, its effects on the matter power spectrum would be more extreme than is currently predicted by the hydrodynamical simulations. The impact of the kSZ on the 𝑆 8 constraint to reduce the uncertainty and shift the value toward Planck is minor given the low signal-to-noise of the kSZ measurements. However, it motivates us to investigate the impact of the kSZ on the 𝑃 ( 𝑘 ) constraints and discuss the possibility of a more extreme suppression.  \nWefindthat incorporating a joint analysis with kSZ results in a significant improvement in the constraint on log 10 𝑀 c , reducing the uncertainty by a factor of ∼ 3 with respect to the WL-only analysis. The joint kSZ and WL data prefer larger values of 𝜃 ej and 𝛾 , and lower values of 𝜇 (see Fig. B2 and Table B2), suggesting that gas is ejected to larger radii, redistributing matter on larger scales (Fig. 2).  \nThis supports the idea of a more extreme feedback scenario, resulting in higher values of 𝑆 8 . The improvement on the WL constraint on 𝑆 8 with the inclusion of kSZ is ∼ 10%. Although this is modest, it is clear that the parameter space of the baryonification model is better constrained, even in this case of a kSZ measurement with signal-to-noise of ∼ 7.  \nWe compare the best-fit models to the measured DES Y3 cosmic shear two-point correlation functions, 𝜉 ± in Fig. 7. The dark green line indicates the best fit for our reanalysis of DES Y3 using their scale cuts, which are indicated by the greeb shaded region. The purple line shows that for the BCEmu joint analysis of all angular scales of the lensing measurement and the kSZ. The lower panels highlight the fractional residuals between the measurements and the model, ( 𝜉 ± -𝜉 model ± )/ 𝜉 model ± , following the same colour scheme. While the fits are indistinguishable at large scales, at small scales, particularly for 𝜉 -, the predictions differ and the best-fit line for the joint analysis has a lower amplitude. Both of these model choices provide a good fit to the data, although their 𝑆 8 values differ by ∼ 1 𝜎 and their non-linear matter spectrum predictions differ substantially. This highlights the degeneracy between a low𝑆 8 cosmology and a higher𝑆 8 cosmology with baryonic effects modelled on non-linear scales. Fig. 8 displays the joint WL + kSZ BCEmu7 constraints on the stacked kSZ radial temperature profile at 98GHz and 150GHz (purple). We verify that like the DES data vector, the WL + kSZ model provides a good fit to the data. The 𝜒 2 red values reported in Table 3 further show that the best-fit models attained from the shear and WL + kSZ analyses are an equally good fit to the datasets.  \nIt is important to note that the constraints on 𝑆 8 that we obtain are dependent on the choice of NLA as the IA model. Appendix C tests the impact of using the Tidal Alignment and Tidal Torquing ('TATT'; Blazek et al. 2019) superspace IA model. In a WL-only analysis with BCEmu7 and using TATT, we find a value of 𝑆 8 ∼ 0 . 5 𝜎 lower than that we obtain in our fiducial analysis with NLA. This shift is found to be consistent across the various baryon models that we consider, as well as with previous findings (e.g. Secco, Samuroff et al. 2022; DES & KiDS Collaborations 2023). Therefore, when considering the results of this work in the context of the 𝑆 8 tension, it is important to bear in mind the existing uncertainty in IA modelling and the shifts in cosmological parameters that can occur as a result.", '6.2 Power suppression constraints': "We investigate how the constraint on the suppression of the matter power spectrum changes with the addition of the kSZ data. The right panel of Fig. 6 shows the mean constraint and 68% confidence level for the BCEmu7 cosmic shear analysis (pink), as shown previously, compared to the joint WL and kSZ analysis (purple). The joint analysis results in a mean suppression that is more extreme at all displayed 𝑘 -scales. For example, at 𝑘 = 2 ℎ / Mpc, the 1 𝜎 bounds on the power suppression from the WL-only analysis range from 5-20%, and with the inclusion of the kSZ, the suppression ranges from 10-25% . Consistent with Schneider et al. (2022), we find that the inclusion of kSZ data in the analysis improves the constraint from weak lensing only to favour more extreme scenarios. Compared to the hydrodynamical simulations, the mean suppression from the joint analysis is more extreme than all of MillenniumTNG, BAHAMAS, FLAMINGO, FABLE and SIMBA at 𝑘 = 2 ℎ / Mpc by more than ∼ 2 . 5 𝜎 , ∼ 1 . 4 𝜎 , ∼ 1 . 9 𝜎 , ∼ 2 . 4 𝜎 and ∼ 0 . 9 𝜎 respectively. Of the simulations, SIMBA shows the best agreement with the WL and kSZ constraint. We note, however, that for BAHAMAS and FLAMINGO we compare to only their fiducial feedback vari-  \nFigure 6. Left: The marginalised posteriors for Ω m and 𝑆 8 attained modelling baryonic feedback with the BCEmu7 model using the DES Y3 cosmic shear dataset only (pink), or a combined analysis of DES Y3 cosmic shear and ACT DR5 kSZ measurements (purple). The inner and outer contours show the 68% and 95% confidence levels, respectively. We compare to the CMB Λ CDM constraint measured by Planck TTTEEE likelihood (Efstathiou & Gratton 2021). Right: The constraints on the suppression of the total matter power spectrum compared to a dark matter-only scenario, 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) , at 𝑧 = 0 when modelling baryonic feedback with the BCEmu7 model using the DES Y3 cosmic shear dataset only (pink), or a combined analysis of DES Y3 cosmic shear and ACT DR5 kSZ measurements (purple). The solid lines show the mean suppression and the shaded regions indicate the 68% confidence levels. For reference, various predictions for the suppression of the matter power spectrum from simulations are over-plotted in black: FLAMINGO (Schaye et al. 2023, solid line); BAHAMAS (McCarthy et al. 2017, dash-dotted line); SIMBA (Dav'e et al. 2019, dotted line), MillenniumTNG (Pakmor et al. 2023, dashed line) and FABLE (doubledot-dashed line, Henden et al. 2018, Bigwood et al. in prep.).  \n<!-- image -->  \neme feedback variations which are in better agreement with our measurements.  \nFinally, we consider our findings in the context of Amon & Efstathiou (2022); Preston et al. (2023), which proposed that the 𝑆 8 tension could be resolved if a more extreme baryon feedback scenario than that predicted by the state-of-the-art hydrodynamical situations existed. These works analysed the KiDS and DES WL data assuming the Planck cosmology on linear scales, and modulating the non-linear power spectrum suppression via the 𝐴 mod parameter. Preston et al. (2023) therefore quantified the small-scale suppression required to resolve the suppression, which corresponds to 𝐴 mod ≈ 0 . 82. In this scenario, the matter power suppression is surpressed enough to reconcile the difference in 𝑆 8 between DES Y3 cosmic shear and the Planck Λ CDM model. The joint WL + kSZ constraints on the suppression of the matter spectrum that we obtain are consistent with the 𝐴 mod = 0 . 82 prediction.", '7 DISCUSSION: CONSISTENCY OF X-RAY AND KINETIC SZ DATA': "OurWL+kSZapproachiscomplementarytoexistingeffortstoward the goal of a complete model of baryonic feedback that is consistent with a wide range of observables. X-ray measurements of cluster gas mass fractions are the primary observable used to calibrate or benchmark many hydrodynamical simulations (e.g. McCarthy et al. 2017; Henden et al. 2018; Kugel et al. 2023). Observations of the tSZ power spectrum and tSZ flux-halo mass relation have also been used to assess the simulations' realism (McCarthy et al.", '2023; Henden et al. 2018; Schaye et al. 2023; Pakmor et al. 2023).': 'In general, while tSZ and X-ray measurements are more sensitive to the inner regions of galaxy groups and clusters, the kSZ effect is well-suited to probe the outskirts of halos, through its sensitivity to low density and low temperature environments (Schaan et al. 2021). Furthermore, X-ray and tSZ observations are typically derived from massive cluster halos ( 𝑀 500 > 10 14 ), while our kSZ measurements represent halos of mass 𝑀 200 ∼ 10 13 , closer to the mass range that WLis most sensitive to. Since the kSZ effect probes the gas in halos of a different mass regime and on different scales to that which is currently used to calibrate feedback effects in simulations, it may allow new insights to be gained.  \nOur weak lensing constraints on the matter power spectrum suppression point to a feedback scenario that is more extreme than most simulations predict. This observation holds in all three flexible model scenarios tested (Fig. 5). The addition of the kSZ data pushes the mean constraint towards an even more disruptive feedback scenario: at 𝑘 = 2 ℎ /Mpc, the fiducial FLAMINGO simulation is disfavoured at ∼ 2 𝜎 and MillenniumTNG at ∼ 2.5 𝜎 . Beyond the comparison to the simulations, our findings point to a more disruptive feedback scenario than inferred from the predicted 𝑃 ( 𝑘 ) suppression using X-ray gas and stellar fraction observations (Grandis et al. 2024), as well as tSZ-mass relation of clusters (To et al. 2024), which could point to interesting differences between the X-ray view of baryon feedback, compared to that from weak lensing and kSZ. We note that a strong feedback scenario was also found using an analysis of the cross-correlation of diffuse X-ray and weak lensing (Ferreira et al. 2023), and hints of a stronger scenario with the cross-correlation of tSZ with galaxies (Pandey et al. 2023).  \nFig. 9 shows the constraints for the WL-only analysis in pink and the WL + kSZ in purple, with the mean relation indicated by a solid line and 68% confidence level as a shaded region. Measurements of halo masses are generally reliant on the assumption of hydrostatic equilibrium and hence may suffer a bias due to nonthermal support in the halo (see, e.g. Rasia et al. 2006). However, estimates have also been obtained from weak lensing data which do not require this assumption (Akino et al. 2022; Mulroy et al.  \n<!-- image -->  \n-  \n-  \n-  \nFigure 7. The DES Y3 cosmic shear two-point correlation function measurements, 𝜉 ± , as a function of angular separation, 𝜃 , (black data points, upper panels) from Amon et al. (2022) and Secco, Samuroff et al. (2022) for each pair of redshift bins, indicated by the panel label. The measurements are scaled by 𝜃 for visual aid. The error bars are calculated as the square root of the diagonal of the analytic covariance matrix. We show the best fit Λ CDMtheoretical predictions to the large angular scale measurements (indicated by the shaded regions) without a baryon model (green line) and the best fit obtained from jointly analysing all angular scales of the DES shear measurements with ACT kSZ data using the BCEmu7 model (purple line). The residuals between the measurement and best-fit model predictions 𝛿𝜉 ±/ 𝜉 ± are shown in the lower panels.', '7.1 Mass dependence of the total baryon fraction': 'As discussed in Section 3.2.3, the baryon fraction, 𝑓 b , measured in halos of mass 𝑀 500 ≈ 10 14 𝑀 ⊙ can be related directly to the matter power spectrum suppression in a manner which is robust to a number of baryonic feedback prescriptions (van Daalen et al. 2020). In this section, we scrutinise the weak lensing constraints on the halo mass dependence of the total baryon fraction in halos with respect to the cosmic baryon fraction, 𝑓 b /( Ω b / Ω m ) . For a given set of parameters, we can use the BCEmu model to compute the baryon fraction through the summation of the integrated gas and stellar profiles (see equation 9 and equation 11) out to 𝑅 500 . For both the WL-only and WL + kSZ BCEmu7 analyses, we compute the 𝑓 b /( Ω b / Ω m ) -𝑀 500 relation for the cosmology and baryonic feedback parameters sampled at each step of the chain.  \n2019; Hoekstra et al. 2015; Eckert et al. 2016). As a result we compare to the 1 𝜎 region of the X-ray HSC-XXL constraint (Akino et al. 2022) in hatched black. The X-ray data has a dependence on cosmology through 𝐸 ( 𝑧 ) , hence we scale the data to the mean cosmology constrained by the WL + kSZ analysis. We also show the baryon fractions measured in BAHAMAS (McCarthy et al. 2017), FLAMINGO (Schaye et al. 2023), MillenniumTNG (Pakmor et al. 2023) and FABLE (Henden et al. 2018, Bigwood et al. in prep.).  \nOverall, the BCEmu7 WL-only and WL + kSZ constraints are in good agreement. The addition of the kSZ reduces the uncertainty by a factor of ∼ 1 . 5 at 𝑀 500 ≈ 10 13 𝑀 ⊙ and prefers lower baryon fractions for all masses. For large groups with 𝑀 500 ⪆ 5 × 10 13 𝑀 ⊙ , we find good consistency between the 68% confidence levels on 𝑓 b /( Ω b / Ω m ) as predicted by BCEmu7 WL-only and WL + kSZ analyses compared to the X-ray data. However, for lower mass groups with masses 𝑀 500 ≈ 10 13 𝑀 ⊙ , the WL + kSZ data prefers a total baryon fraction that is lower than the X-ray data by ∼ 1 . 4 𝜎 . With the exception of BAHAMAS which is in best agreement, all of the simulations predict higher values of 𝑓 b /( Ω b / Ω m ) for all masses compared to the mean of the WL + kSZ constraint. At 𝑀 500 ∼ 4 × 10 13 𝑀 ⊙ , FLAMINGO predicts ∼ 1 . 5 𝜎 higher values, MilleniumTNG ∼ 3 . 1 𝜎 and FABLE ∼ 1 . 4 𝜎 .  \nIn Appendix B4, we compare the relationship between the  \nFigure 8. Stacked ACT DR5 kSZ temperature profile measurements (black points) as a function of angular radius, 𝜃 , at 98GHz and 150GHz (Schaan et al. 2021) and the best fit from the joint analysis with the DES Y3 cosmic shear using the BCEmu7 model (purple). The model profiles are convolved with the f90 and f150 beam profile for comparison to the 98 GHz and 150 GHz data, respectively.  \n<!-- image -->  \nbaryon fraction and power suppression that our WL + kSZ constrains to that proposed by van Daalen et al. (2020), which is a good fit to many of the hydrodynamical simulations. We note that the BCEmu model does not impose a prior on this relationship and therefore provides a route to place constraints using data. Given that our constraints differ from the simulation-based relationship, we can conclude that either there are unaccounted for systematics in the WLandkSZdata,theBCEmumodelallowsnon-physicalscenarios, or simulations do not currently capture the full possible range of feedback effects, thereby overestimating the relationship (see e.g. Debackere et al. 2020). Future work is needed to further understand this relationship using data.  \nThe SP(k) baryon model also provides a direct mapping from the matter power spectrum suppression to the baryon fraction. In Appendix B1 we show the mean and 68% confidence levels for the baryon fraction attained from the WL-only SP(k) analyses. Similar to the BCEmu case, there is good agreement of the lensing analysis with SP(k):wide and X-ray data at high masses 𝑀 500 ∼ 1 × 10 14 𝑀 ⊙ . In agreement with the findings with BCEmu, we find that lensing prefers slightly lower baryon fractions for halos 𝑀 500 ∼ 10 13 𝑀 ⊙ , with the SP(k):wide analysis lying lower by ∼ 1 . 3 𝜎 . This highlights that the lensing data alone, when analysed with flexible modelling choices, prefers a lower baryon content in halos to that measured by X-ray observations and that predicted by simulations.', '7.2 Mass dependence of the gas mass fraction': "The majority of baryonic mass in galaxy groups and clusters exists as diffuse gas and measurements of the fraction of the halo mass in gas, 𝑀 gas / 𝑀 500 , are also sensitive to the matter power spectrum suppression (Schneider & Teyssier 2015; van Daalen et al. 2020;  \nAric'o et al. 2023). The hot intracluster medium is X-ray luminous, hence studies of the X-ray emissivity of groups and clusters allow measurements of the gas mass to be derived.  \nIn this section, we compare X-ray derived measurements of the 𝑀 gas / 𝑀 500 -𝑀 500 relation to those constrained by the BCEmu7 WL-only and joint WL + kSZ analyses in this work. As with the baryon fractions, we calculate the 𝑀 gas / 𝑀 500 -𝑀 500 relation at each step in the chain by integrating equation 9 to 𝑟 500 . The right panel of Fig. 9 shows the mean and 68% confidence levels on the gas fraction-halo mass relation for both the WL-only (pink) and WL + kSZchains(purple). We compare to the 1 𝜎 region of the X-ray HSCXXL constraint (Akino et al. 2022) in hatched black (also scaled to the mean cosmology obtained in the WL + kSZ analysis) and the baryon fractions measured in BAHAMAS (McCarthy et al. 2017), FLAMINGO (Schaye et al. 2023), MillenniumTNG (Pakmor et al. 2023) and FABLE (Henden et al. 2018, Bigwood et al. in prep.).  \nIn consistency with the baryon fraction result, we find that the BCEmu7WL-onlyandWL+kSZconstraintsareingoodagreement with each other and the X-ray data for massive groups of 𝑀 500 ⪆ 5 × 10 13 𝑀 ⊙ . For groups of lower masses we find that the joint WL + kSZ analysis significantly improves the constraining power on the gas fractions, with the 68% confidence levels on 𝑀 gas / 𝑀 500 at 𝑀 500 ≈ 10 13 𝑀 ⊙ shrinking by a factor of ∼ 2 between the WL-only and WL + kSZ analyses. At 𝑀 500 ≈ 10 13 𝑀 ⊙ , the WL + kSZ data constrains a total baryon fraction that is lower than the X-ray data by ∼ 1.6 𝜎 . The fiducial BAHAMAS simulation is in best agreement with the WL + kSZ constraint, with the remaining simulations predicting higher values of 𝑀 gas , 500 / 𝑀 500 at all halo masses. At 𝑀 500 ∼ 4 × 10 13 𝑀 ⊙ , the fiducial FLAMINGO simulation predicts ∼ 2.3 𝜎 higher values, MilleniumTNG ∼ 7.2 𝜎 and FABLE ∼ 2.1 𝜎 .  \nThe baryon and gas fraction constraints by the WL + kSZ data imply that lower mass groups are expelling a greater amount of baryonic matter due to feedback than predicted by X-ray measurements. When combined with the findings of Section 6.2 that WL + kSZ predicts a greater mean matter power spectrum than all of the simulations, we build a consistent picture: that X-ray observations constrain a weaker feedback scenario than that preferred by WL + kSZ, subject to the uncertainties we have discussed.", '8 SUMMARY AND CONCLUSIONS': 'Weak lensing measurements at small angular scales are statistically powerful. Not only do they offer a window to test Λ CDM in the non-linear regime, but also to constrain baryonic feedback effects and benchmark hydrodynamical simulations. Extracting accurate cosmological information relies on accurate modelling of physical processes associated with baryons which can re-distribute matter on small scales. Otherwise, unmodelled baryonic effects can bias the cosmological constraints from lensing analyses (Amon & Efstathiou 2022; Preston et al. 2023) or severely limit their precision (Amon et al. 2022; DES & KiDS Collaborations 2023). The aim of this work has been two-fold. Firstly, we test four baryon model approaches, comparing their ability to constrain cosmological parameters and the suppression of the matter power spectrum. Secondly, we perform a joint analysis of DES cosmic shear with ACT kSZ measurements. We demonstrate that a combined weak lensing and kSZ analysis provides an exciting opportunity to not only improve constraints on cosmological parameters, but also on astrophysical effects. The main results of this study are:  \n- · We perform the first mock baryon model comparison for cos-', '18 L. Bigwood, A. Amon et al.': "Figure 9. Left: The total baryon fraction, 𝑓 b / Ω b / Ω m , as a function of halo mass, 𝑀 500 , where Ω b / Ω m is the cosmic baryon fraction and 𝑓 b is the fraction of mass in baryons to the total halo mass in groups and clusters. We plot the constraint attained from analysing the DES Y3 cosmic shear data with BCEmu7 (pink) in addition to the constraint attained from the joint WL + kSZ analysis (purple). Solid lines show the mean baryon fraction halo mass relation, and the shaded regions enclose the 1 𝜎 constraints. Right: The fraction of mass in gas to the total halo mass in groups and clusters, 𝑀 gas , 500 / 𝑀 500 , as a function of halo mass, 𝑀 500 . We plot the constraint attained from analysing the DES Y3 cosmic shear data with BCEmu7 (pink) in addition to the constraint attained from the joint WL + kSZ analysis (purple). Solid lines show the mean gas fraction halo mass relation, and the shaded regions enclose the 1 𝜎 constraints. In both panels we plot the X-ray constraints from HSC-XXL 1 𝜎 (Akino et al. 2022) as the black hatched region, scaled to the mean cosmology obtained from the WL + kSZ analysis. For reference, we also plot the measurements from BAHAMAS (McCarthy et al. 2017, dash-dotted line); FLAMINGO (Schaye et al. 2023, solid line); MillenniumTNG (Pakmor et al. 2023, dashed line) and FABLE (doubledot-dashed line, Henden et al. 2018, Bigwood et al. in prep.).  \n<!-- image -->  \n/circledot  \n<!-- image -->  \n/circledot  \nmic shear to assess the robustness of the models to recover the underlying cosmology in a DES Y3-like analysis. We consider three different mock 'feedback scenarios' to test four baryon strategies; a halo model approach, a simulation-based emulator, an analytical N-body simulation model and the approach of discarding small angular scales. In general, using restrictive modelling choices which do not capture the input matter power spectrum suppression can underestimate the recovery of 𝑆 8 . Given the spread in the true suppression of the matter power spectrum as predicted by simulations and the lack of observational constraints on this quantity, model flexibility is crucial to ensure accurate cosmology, but it comes at a cost: the uncertainty on 𝑆 8 can degrade by up to a factor of ∼ 2 with different model choices.  \n- · We analyse the DES Y3 cosmic shear data with the three baryon feedback models, with the scale cut approach, and with no feedback model. Wefindthat each baryon model provides a good fit to the DES Y3 data, but the measured value of 𝑆 8 varies by ∼ 0.5-2 𝜎 and the errorbar by a factor of ∼ 1.5.\n- · We output the posterior for the suppression of the matter power spectrum using the three models, with and without their informative priors. For all three of the models, without their informative priors, the mean suppression of the matter power spectrum constrained using DES Y3 WL is more extreme than the prediction from the hydrodynamical simulations considered .\n- · We jointly analyse the DES Y3 cosmic shear with kSZ measurements from ACT DR5 and CMASS, in order to constrain the BCEmu7 model parameters. We find a slightly higher value of 𝑆 8 = 0 . 823 + 0 . 019 -0 . 020 , compared to the value attained by the WL-only analysis. If instead we analyse the DES Y3 cosmic shear jointly with X-ray baryon fraction constraints using SP(k), we find a lower value of 𝑆 8 = 0 . 806 + 0 . 015 -0 . 013 .\n- · The kSZ significantly improves the constraint on the suppression of the matter power spectrum from WL. The joint WL+kSZ predicts a more extreme suppression of the matter power spectrum  \nthan the WL scenario , with a mean constraint predicting a greater suppression than the fiducial BAHAMAS, MillenniumTNG, fiducial FLAMINGO and FABLE simulations. At 𝑘 ≳ 1 ℎ / Mpc, only SIMBA falls within the 1 𝜎 bound.  \n- · We constrain the baryon fraction-halo mass and the gas fraction-halo mass relations using WL + kSZ data . Both the baryon fraction and gas fraction is consistent with that from Akino et al. (2022) X-ray data within the 68% confidence level for groups 𝑀 500 ≈ 10 14 𝑀 ⊙ . However for lower mass groups 𝑀 500 ≈ 10 13 𝑀 ⊙ the baryon fraction lies ∼ 1.4 𝜎 lower than the X-ray data and the gas fraction lies ∼ 1.6 𝜎 lower.\n- · Our constraints on the matter power spectrum suppression, baryon fractions and gas fractions all point towards a tension between the feedback of groups and clusters predicted by weak lensing + kSZ and X-rays, or X-ray calibrated models.  \nThe next generation of shear surveys, such as Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST; LSST Dark Energy Science Collaboration 2012), Euclid (Laureijs et al. 2011) and the Nancy Grace Roman Space Telescope (Spergel et al. 2015) will deliver unprecedented statistical power to test the cosmological model in the non-linear regime ( 𝑘 > 0 . 1 ℎ Mpc -1 ) using weak lensing. This work highlights the importance of determining an accurate model for baryonic effects, which is flexible enough to not bias cosmological constraints. It will be crucial to test models of feedback with thorough mock analyses. In order for the model uncertainty in weak lensing analyses to not limit the statistical power of the survey, it is crucial to either incorporate external probes of the gas content to constrain the additional baryon parameters, or to reduce the uncertainty in simulation-based priors. The latter requires a consistent picture for baryonic feedback effects, supported by a range of observations.  \nWe demonstrate that joint analyses of gas measurements with weak lensing data not only improve cosmological constraints, but provide valuable constraints on astrophysical feedback models and  \nbenchmark hydrodynamical simulations. We find a consistent picture that could imply that the WL and kSZ data is in tension with the X-ray measurements, and as a result, the predictions from simulations calibrated to X-ray data. This is not particularly surprising as X-ray measurements are generally sensitive to the hot gas content in the inner regions of clusters, compared to the outer regions and lower mass halos that kSZ measurements probe.  \nLooking ahead, kSZ measurements as a function of mass and redshift will provide a handle for improved baryonic feedback models. Spectroscopic galaxy surveys, such as the Dark Energy Spectroscopic Instrument (DESI; DESI Collaboration 2016) and the Prime Focus Spectrograph (PSF; Sugai et al. 2012), will greatly increase the sample size of galaxy catalogues, in combination with the state-of-the-art CMB observations, for example, from Simons Observatory (Ade et al. 2019).", 'ACKNOWLEDGEMENTS': "Leah Bigwood is supported by a Science and Technology Facilities Council studentship. Alexandra Amon was supported by a Kavli Fellowship for part of the duration of this work. Jaime Salcido and Ian McCarthy were supported by an European Research Council (ERC) Consolidator Grant under the European Union's Horizon 2020 research and innovation programme (grant agreement No 769130).  \nAuthor contributions are as follows. L. Bigwood and A. Amon performed the analysis and wrote the bulk of the manuscript. The additional core authors of A. Schneider, J. Salcido and I. McCarthy contributed BCEmu and SP(k) code and continuous insight. C. Preston and D. Sanchez assisted in the early stages, including the mock analysis. D. Sijacki provided useful discussions and detailed comments on the manuscript. E. Schaan, S. Ferraro and N. Battaglia provided the ACT kSZ dataset. The paper has gone through internal review by the DES collaboration, with A. Chen and S. Dodelson serving as excellent internal reviewers. The remaining authors have madecontributions to the DES instrumentation, data collection, data processing and calibration, in addition to other code and validation tests.  \nThis project has used public archival data from the Dark Energy Survey. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and AstroParticle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundac¸˜ao Carlos Chagas Filho de Amparo 'a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cient'ıfico e Tecnol'ogico and the Minist'erio da Ciˆencia, Tecnologia e Inovac¸˜ao, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey.  \nThe Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energ'eticas, Medioambientales y Tecnol'ogicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenossische Technische Hochschule (ETH) Zurich, Fermi National Accelerator Laboratory, the University of  \nIllinois at Urbana-Champaign, the Institut de Ci'encies de l'Espai (IEEC/CSIC), the Institut de F'ısica d'Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universitat Munchen and the associated Excellence Cluster Universe, the University of Michigan, NSF NOIRLab, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, Texas A&M University, and the OzDES Membership Consortium.  \nBased in part on observations at NSF Cerro Tololo InterAmerican Observatory at NSF NOIRLab (NOIRLab Prop. ID 2012B-0001; PI: J. Frieman), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.  \nThe DES data management system is supported by the National Science Foundation under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MICINN under grants ESP201789838, PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCAprogramoftheGeneralitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciˆencia e Tecnologia (INCT) do e-Universo (CNPq grant 465376/2014-2).  \nThis manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics.", 'DATA AVAILABILITY': 'The data underlying this article may be shared on reasonable request to the corresponding authors.', 'REFERENCES': 'Adamek J., et al. 2023, J. Cosmology Astropart. Phys., 2023, 035 Ade P., et al. 2019, J. Cosmology Astropart. Phys., 2019, 056 Akino D., et al. 2022, PASJ, 74, 175 Amodeo S., et al. 2021, Phys. Rev. D, 103, 063514 Amon A., Efstathiou G., 2022, MNRAS, 516, 5355 Amon A., et al. 2022, Phys. Rev. D, 105, 023514 Amon A., et al. 2023, MNRAS, 518, 477 Aric`o G., et al. 2020, MNRAS, 495, 4800 Aric`o G., et al. 2023, A&A, 678, A109 Asgari M., et al. 2021, A&A, 645, A104 Behroozi P., et al. 2019, MNRAS, 488, 3143 Bennett C. L., et al., 2003, Astrophys. J. Suppl., 148, 1 Blazek J. A., et al. 2019, Phys. Rev. D, 100, 103506 Bridle S., King L., 2007, New Journal of Physics, 9, 444 Chen A., et al. 2023, MNRAS, 518, 5340 Chisari N. E., et al. 2019, The Open Journal of Astrophysics, 2, 4 DES KiDS Collaborations 2023, The Open Journal of Astrophysics, 6, 53 DES Collaboration 2022, Phys. Rev. D, 105, 023520 DESI Collaboration 2016, arXiv e-prints, p. arXiv:1611.00036 Dav´e R., et al. 2019, MNRAS, 486, 2827 Debackere S. N. B., Schaye J., Hoekstra H., 2020, MNRAS, 492, 2285 Dubois Y., et al. 2016, MNRAS, 463, 3948', 'APPENDIX A: MOCK ANALYSIS': 'In this section, we provide more detail on the mock analysis presented in Sec. 4.2. We create the synthetic cosmic shear data using the best-fit cosmological parameters attained from the DES Y3 joint lensing and clustering analysis (DES Collaboration 2022) 12 . As described in Sec. 4.2, the mocks were created with HM20 as the dark matter-only non-linear correction to the matter power spectrum, unless otherwise stated. We create data vectors with three variants for the impact on the non-linear power spectrum of baryonic feedback; we use the power suppression predicted by BAHAMAS 8.0, cosmo-OWLS 8.5 and 𝐴 mod = 0 . 82. In this section we present a more detailed summary of the dark matter and baryonic feedback model validation tests ran using synthetic data vectors. We refer to Table A1 which presents the /m.pc/u.pc/l.pc/t.pc/i.pc/n.pc/e.pc/s.pc/t.pc mean values of 𝑆 8 attained for the full mock suite, reporting the relative shifts from the input cosmology. The results are also summarised in in Fig. A2, where we additionally plot the mean and best-fit Ω m and 𝜎 8 for each mock test.  \nWe validate the choice to model the dark matter-only power spectrum with HM20 using dark matter-only mocks created with HM20 and EE2. Here we summarise the findings of mock tests 1-4.  \n- · This test analyses a mock with the same model choices used to create it and is useful for identifying projection effects. We find these to be present for Ω m and 𝜎 8 , which are under- and overestimated, respectively, by ∼ 0.5 𝜎 .\n- · When modelling dark-matter with HM20 in a HM20 generated mock, allowing the neutrino mass to vary with the prior " 𝑚 𝜈 : [0.06,6] eV improves the recovery of Ω m and 𝜎 8 . It however decreases the accuracy in the recovery of 𝑆 8 , resulting in an underestimation of 0 . 5 𝜎 .\n- · Modelling the non-linear dark matter-only power spectrum with HM20inaEE2generated mock over estimates 𝑆 8 by ∼ 0 . 4 𝜎 , which is reduced marginally to ∼ 0 . 2 𝜎 when scale cuts are applied. We note that DES & KiDS Collaborations (2023) find a smaller bias in the mean of ∼ 0 . 1 𝜎 when scale cuts are used. In this mock analysis there are two set-up differences that could explain this: we use a  \n12 Ω m = 0 . 3380, Ω b = 0 . 0450, 10 -9 𝐴 s = 1 . 8418, ℎ = 0 . 615, 𝑛 s = 0 . 949, 𝑆 8 = 0 . 7805, 𝜎 8 = 0 . 7353.  \nTable A1. Summary of the mock tests we performed to validate the modelling of the non-linear matter power spectrum. Mock tests 1-3 validate the non-linear correction to the matter power spectrum due to dark matter only. The remaining tests validate the baryon feedback mitigation strategies we test in this work; scale cuts (4-7), HM20 (8-14), SP(k) (15-20) and BCEmu (21-29). We test analysing the mocks with various prior and analysis choices for each of the four baryon model approaches. If not otherwise specified, the mocks were created with HM20 as the dark matter-only non-linear correction to the matter power spectrum, and three \'baryonic feedback scenarios\'; BAHAMAS, cosmo-OWLS and 𝐴 mod = 0 . 82. The \'mock\' column therefore states the model/baryonic feedback scenario used to create the mock and the \'model\' column labels the model/analysis choices used to analyse the mock\'. For example, a mock of \'EuclidEmu\' and model of HM20-DM\' denotes a mock created with Euclid Emulator for the non-linear correction to the matter power spectrum, and analysed with HM20. We report the 𝑆 8 constraints and 68% confidence level using the mean-marginal approach 𝑆 8 , with Δ 𝑆 8 quantifying the offset from the true 𝑆 8 = 0 . 7805, as a fraction of the 1 𝜎 error.  \nhigher value of Ω m to create the mock and we marginalise over 𝐴 s , instead of 𝑆 8 .  \nNext we validate each of our four baryon models (scale cuts, HM20 with free Θ AGN , BCEmu and SP(k)) and their respective prior and analysis variants. We test the modelling choices on three \'baryonic feedback scenarios\' of increasing \'extremity\' in terms of their impact on the matter power spectrum suppression: BAHAMAS 8.0, cosmo-OWLS 8.5 and 𝐴 mod = 0 . 82. Here we summarise the findings of the mock tests 5-29 and Fig. A1. Note that in the figure, mocks are labelled in the format \'model used to create the mock: model/analysis choices used to analyse the mock\'. For example, \'BAHAMAS: HM20 BAHAMAS\' denotes a mock created  \nwith BAHAMAS-like baryon feedback scenario and analysed with HM20 with the BAHAMAS Θ AGN = 7 . 6 -8 . 0 prior.  \n- · Optimised scale cuts underestimates 𝑆 8 by 0 . 5 -1 . 2 𝜎 in the three feedback scenarios we test, with the tension worsening with more extreme feedback. As noted in Sec. 4.2, this is as expected given our choice of feedback scenarios used to build the synthetic data. The OWLS-AGN scenario, which was used to decide which angular scales of the measurement would be discarded in the DES Y3 analysis predicts a less extreme feedback scenario than the three scenarios used here.\n- · The left panel of Fig. A1 compares prior choices of the HM20 model with Θ AGN as a free parameter. We find that HM20 with the BAHAMAS calibrated prior choice of Θ AGN = 7 . 6 -8 . 0 un-  \n<!-- image -->  \n<!-- image -->  \nFigure A1. Marginalised posteriors for 𝑆 8 and the baryon model parameter using each mock for analyses using the HM20 (upper), BCEmu (middle) and SP(k) (lower) models. The mocks are labelled in the format \'model used to create the mock: model/analysis choices used to analyse the mock\'. The inner and outer contours show the 68% and 95% confidence levels, respectively and the dashed line indicates the input cosmology.  \n<!-- image -->  \nTable B1. The constraints on the SP(k) parameters attained by the WL-only SP(k) analyses with both the wide and X-ray informed prior choices. We report the mean and 68% confidence levels.  \nstimates 𝑆 8 by ∼ 0 . 7 -1 . 2 𝜎 , with the tension worsening with a more extreme feedback scenario. Using the wide prior choice of Θ AGN = 7 . 3 -9 . 0 overestimates 𝑆 8 in BAHAMAS 8.0 and cosmo-OWLS 8.5 by ∼ 0 . 5 𝜎 , but still underestimates 𝑆 8 in an 𝐴 mod = 0 . 82-like feedback scenario by ∼ 0 . 2 𝜎 .  \n- · The right panel of Fig. A1 shows that SP(k) with the X-ray informed prior demonstrates the best recovery of 𝑆 8 in a BAHAMAS 8.0-like feedback scenario of all the baryonic feedback modelling choices. In more extreme feedback scenarios, it however underestimates 𝑆 8 by ∼ 0 . 3 -1 . 6 𝜎 . Using the wide prior choice improves the recovery of 𝑆 8 to be within ∼ 0 . 4 𝜎 for the three feedback scenarios.\n- · The central panel of Fig. A1 demonstrates that BCEmu7 overestimates 𝑆 8 in both BAHAMAS 8.0 and cosmo-OWLS 8.5 baryonic feedback scenarios by ∼ 0 . 4 𝜎 , but underestimates 𝑆 8 in an 𝐴 mod = 0 . 82-like feedback scenario by ∼ 0 . 2 𝜎 . BCEmu1 underestimates 𝑆 8 by 0 . 3 -1 . 2 𝜎 , again with the tension worsening when testing more extreme feedback scenarios.\n- · In the mock tests which over-estimate 𝑆 8 , freeing the neutrino mass with the prior " 𝑚 𝜈 : [0.06,6] eV appears to improve recovery of the input cosmology. The shifts we see in 𝑆 8 are however consistent with the ∼ 0 . 2 𝜎 shift to lower values we find when testing the dark matter-only modelling of HM20 with and without free neutrinos, i.e. tests 1 and 2.  \nUltimately we find that choosing the restrictive analysis variant of a model tends to result in an under-estimation of 𝑆 8 . This motivates the use of more flexible modelling choices in cosmic shear analyses, since these tend to result in improved recovery of 𝑆 8 , however this is at the consequence of a loss in precision.', 'B1 SP(k)': 'The mean and 68% confidence levels on the SP(k) parameters constrained by the WL-only SP(k) analyses with the wide and X-ray informed prior choices are reported in Table B1. The corresponding marginalised posteriors for these analyses are shown in the upper panel of Fig. B1. Since SP(k) directly maps between the baryon fraction in groups and clusters and the matter power spectrum suppression, we also show the constraints on the 𝑓 b / Ω b / Ω m -𝑀 500 relation predicted when analysing the DES cosmic shear with the two SP(k) analysis variants in the lower panel of Fig. B1.', 'B2 Baryonification': "The mean and 68% confidence levels on the baryonification parameters constrained by the WL-only BCEmu1, BCEmu3 and BCEmu7 analyses, in addition to the WL + kSZ BCEmu7 joint analysis (for which the average halo mass of the kSZ sample 𝑀 h , 200 is also reported) are given in Table B2. The corresponding marginalised posteriors for these analyses are shown in Fig. B2.  \n9  \nFigure A2. Summary of the 1D marginalised constraints on 𝑆 8 , Ω m and 𝜎 8 attained by our mock analysis with respect to the input cosmology. We plots the mean and 68% confidence levels as listed in Table A1, with crosses showing the /m.pc/u.pc/l.pc/t.pc/i.pc/n.pc/e.pc/s.pc/t.pc best-fit result. The mocks were created with DES-Y3 covariance, with the input cosmology indicated by the vertical yellow line. The top panel shows the results of validating the non-linear correction to the matter power spectrum due to dark matter only. The mocks are labelled in the format 'model used to create the mock: model/analysis choices used to analyse the mock'. For example, 'EuclidEmu: HM20-DM' denotes a mock created with Euclid Emulator for the non-linear correction to the matter power spectrum, and analysed with HM20. The remaining panels validate the four methods of mitigating for baryonic feedback we test in this work; optimized scale cuts, HM20, BCEmu and SP(k). If not otherwise specified, the mocks were created with HM20 as the dark matter-only non-linear correction to the matter power spectrum, and three 'baryonic feedback scenarios'; BAHAMAS, cosmo-OWLS and 𝐴 mod = 0 . 82. We test analysing the mocks with various prior and analysis choices for each of the four baryon model approaches. Each mock was analysed with HM20 as the dark matter-only non-linear correction to the matter power spectrum.  \n<!-- image -->  \nFigure B1. SP(k) marginalised posteriors when analysing the DES Y3 cosmic shear data with the wide prior on SP(k) (light green) and X-ray informed prior (dark green). We show the 68% and 95% confidence levels of the 𝑆 8 and the three model parameters (upper panel) and the mean total baryon fraction, 𝑓 b / Ω b / Ω m , and 1 𝜎 uncertainty as a function of halo mass, 𝑀 500 (lower). The X-ray prior is derived from HSC-XXL 1 𝜎 constraints (black hatched Akino et al. 2022), in this plot scaled to the mean cosmology obtained from the 𝑆𝑃𝑘 : X-ray analysis.  \n<!-- image -->  \n<!-- image -->  \n/circledot", 'B3 Prior choice for the halo mass of the kSZ sample': "Given the significant scatter in the literature, we choose to include an additional model parameter in the analysis, 𝑀 h , 200 , corresponding to the mean 𝑀 200 of the CMASS sample, with a prior range provided in Table 1. For example, the stacked stellar mass-halo mass relation of Sonnenfeld et al. (2019) of CMASS galaxies derived from HSC galaxy-galaxy lensing measurements, when combined with the stellar mass distribution of CMASS, implies a mean 𝑀 200 ≈ 0 . 5 × 10 13 M ⊙ , whereas the abundance matching methods of Kravtsov et al. (2018), Moster et al. (2018), and Behroozi et al. (2019) imply mean 𝑀 200 values of 2 . 7 × 10 13 M ⊙ , 4 . 4 × 10 13 M ⊙ , and  \nTable B2. The constraints on the baryonification parameters attained by the WL-only BCEmu1, BCEmu3 and BCEmu7 analyses, as well as the BCEmu7 WL + kSZ BCEmu7 analysis. For the WL + kSZ analysis we also report the constraint on the average halo mass of the kSZ sample (10 13 𝑀 ⊙ ). We report the mean and 68% confidence levels.  \n6 . 6 × 10 13 M ⊙ , respectively. A HOD-based analysis of the clustering of BOSS CMASS galaxies by White et al. (2011) found a mean 𝑀 200 of 3 . 6 × 10 13 M ⊙ . Given this large study-to-study variance, we choose a flat prior range of [ 0 . 5 -7 ] × 10 13 M ⊙ and marginalise over this parameter.  \nThe combined WL + kSZ analysis constrains 𝑀 h , 200 = 4 . 098 + 1 . 548 -1 . 639 × 10 13 M ⊙ , indicating we are prior constrained at the 2 𝜎 level. For comparison to the remainder of this work which generally quotes halo masses in 𝑀 500 , this corresponds to 𝑀 ℎ, 500 = 3 . 01 × 10 13 M ⊙ , assuming a NFW profile and a concentration-mass relation from Dutton & Macci'o (2014).  \nIn this appendix we explore the impact on our cosmology and baryon model parameter constraints of choosing a wider prior 𝑀 h , 200 : [ 0 . 8 , 30 ] × 10 13 M ⊙ . We also consider a fixed 𝑀 h , 200 analysis using the mean mass of the CMASS sample as determined by Schaan et al. (2021) and used in (Amodeo et al. 2021): 𝑀 h , 200 = 3 × 10 13 M ⊙ . Fig. B3 shows the marginalised posteriors for Ω m , 𝑆 8 , log 10 ( 𝑀 c ) , 𝜃 ej and 𝑀 h , 200 in a WL + kSZ analysis with the different prior choices on 𝑀 h , 200 . The halo mass, 𝑀 h , 200 , is correlated with log 10 ( 𝑀 c ) and in the case of a wider prior, both 𝑀 h , 200 and log 10 ( 𝑀 c ) prefer higher values. However, 𝑆 8 is relatively stable to the halo mass prior, shifting by 0.3 𝜎 .", 'B4 The relationship between baryon fraction and matter power spectrum suppression': "It has been demonstrated that in hydrodynamical simulations, the mean baryon fraction in halos of 𝑀 ∼ 10 14 M ⊙ can be predictive for the suppression of the matter power spectrum due to baryonic feedback effects, robust to a number of feedback prescriptions (van Daalen et al. 2020). This relationship can be described by an empirical fitting function relating the mean baryon fraction, 𝑓 b / Ω b / Ω m , measured within 𝑅 500 for halos of mass 𝑀 500 = 10 14 𝑀 ⊙ to the suppression of the matter power spectrum 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) at 𝑘 = 1 ℎ Mpc -1 (equation 5 in their work). The best-fit relation, shown as the solid grey line in Fig. B4, was fit to the cosmo-OWLS (Le Brun et al. 2014) and BAHAMAS (McCarthy et al. 2017) simulations and is accurate to 1% for the simulations they test, shown as the grey shaded region in Fig. B4. In this section, we discuss where constraints of our analysis lie in the 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) -𝑓 b / Ω b / Ω m plane with respect to the van Daalen et al. (2020) relation. In Fig. B4 we plot the result for the WL-only and WL + kSZ analyses with BCEmu7 baryon modelling, as well as the WL-only results with SP(k)'s wide and X-ray informed prior choices. For comparison, we also show the simulations discussed throughout this work: FLAMINGO (Schaye et al. 2023); BAHAMAS (McCarthy et al.  \nFigure B2. The marginalised posteriors for 𝑆 8 , the seven baryonification parameters, and the average halo mass og the kSZ sample 𝑀 h , 200 attained by the WL-only BCEmu1 (red solid line), BCEmu3 (light purple dotted line) and BCEmu7 (pink dashed line) analyses, as well as the WL + kSZ BCEmu7 (dark purple line) analysis. The inner and outer contours show the 68% and 95% confidence levels respectively. The grey dashed lines show the values the BCEmu parameters are fixed to in the case of using BCEmu1 and BCEmu3, and the mean halo mass of the CMASS sample reported by Schaan et al. (2021).  \n<!-- image -->", "2017); SIMBA (Dav'e et al. 2019), MillenniumTNG (Pakmor et al. 2023) and FABLE (Henden et al. 2018, Bigwood et al. in prep.).": 'This van Daalen et al. (2020) relationship is at the core of the the SP(k) model (Salcido et al. 2023). Given that the ANTILLES simulations follow this relationship, the mean baryon fraction is used to calibrate their emulator. Naturally, both of our WL-only analyses with SP(k) fall on the van Daalen et al. (2020) relation to within 1%. However, the BCEmu model does not enforce the van Daalen et al. (2020) relationship between the baryon fraction and matter power suppression, and allows for greater flexibility in the  \nimpact of feedback on the matter distribution, including scenarios which, according to the simulations, could be deemed unphysical. Both the WL-only and WL + kSZ analyses with BCEmu give a mean constraint below the van Daalen et al. (2020) relation. That is, they allow for a more extreme matter power spectrum suppression for their predicted mean baryon fraction predicted. In the future, it will be important to test this relationship with a range observations, including those probe lower halo masses (see also Pandey et al. 2023).  \n<!-- image -->  \nFigure B3. Upper: The stacked kSZ temperature profile at 98GHz as a function of angular radius, 𝜃 , centred on the group or cluster (bottom) when varying the mean halo mass of CMASS galaxy sample modelled, 𝑀 h , 200 , within the limits of the fiducial prior choice 𝑀 h , 200 : [ 5 × 10 12 , 7 × 10 13 ] 𝑀 ⊙ . The ACT CMASS measurements at 98GHz are shown as the black data points in the bottom panels and the model profiles are convolved with the f90 beam profile for comparison. Lower: The marginalised posteriors for 𝑆 8 , log 10 ( 𝑀 c ) and 𝑀 h , 200 in a WL + kSZ analysis with different prior choices on 𝑀 h , 200 . We show the prior used in the fiducial analysis 𝑀 h , 200 : [ 5 × 10 12 , 7 × 10 13 ] 𝑀 ⊙ (purple), as well as a wide prior 𝑀 h , 200 : [ 8 × 10 12 , 3 × 10 14 ] 𝑀 ⊙ (light blue) and fixed 𝑀 h , 200 = 3 × 10 13 ] 𝑀 ⊙ (dark blue, also shown as the dashed grey line). The inner and outer contours show the 68% and 95% confidence levels respectively.  \n<!-- image -->', 'APPENDIX C: INTRINSIC ALIGNMENT (IA) MODEL': "In this appendix we consider the impact of the IA model choice on our results. Throughout this work, we use the NLA IA model, which is a sub-space of the TATT model. A WL-only analysis using BCEmu7 with TATT constrains 𝑆 8 = 0 . 802 + 0 . 028 -0 . 024 , lying ∼ 0 . 5 𝜎  \nFigure B4. The relation between the matter power spectrum suppression 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) at 𝑘 = 1 ℎ Mpc -1 and the mean baryon fraction, 𝑓 b / Ω b / Ω m , measured within 𝑅 500 for halos of mass 𝑀 500 = 10 14 𝑀 ⊙ . We plot the empirical best-fit relation of van Daalen et al. (2020) as the grey solid line, with the 1% accuracy on the relation's ability to predict 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) shown as the shaded grey region. We plot the constraints obtained in our analysis in the 𝑃 ( 𝑘 )/ 𝑃 DMonly ( 𝑘 ) -𝑓 b / Ω b / Ω m plane. We show the WL-only (pink starred) and WL + kSZ (pruple starred) analyses with BCEmu7, and the WL-only results with 𝑆𝑃𝑘 , showing both the wide (light green starred) and X-ray (dark green starred) prior choices. We also plot the simulations FLAMINGO (Schaye et al. 2023); BAHAMAS (McCarthy et al. 2017); SIMBA (Dav'e et al. 2019), MillenniumTNG (Pakmor et al. 2023) and FABLE (Henden et al. 2018, Bigwood et al. in prep.) as the black datapoints.  \n<!-- image -->  \nlower than that constrained using NLA. This shift in 𝑆 8 between the two IA models is consistent when using BCEmu1 instead of BCEmu7, and similar to that found in previous work when using HM20 DES & KiDS Collaborations (2023) and scale cuts (Secco et al. 2022; Amon et al. 2023). As in the literature, both IA model choices give comparable 𝜒 2 red values (see Tab.3). In future work, it is important to determine the more accurate IA model. Here, we use this test as validation that our baryon model choice and the results of our baryon model comparison are independent of the choice of IA model.", 'APPENDIX D: SAMPLING ALGORITHM CHOICE': 'In this appendix we compare the cosmological parameter estimates attained using the multinest and polychord samplers. Throughout the work, we use the multinest settings 𝑛 live = 500, efficiency = 0.3, tolerance = 0.1, constant efficiency = False, max. iterations = 50,000 for computing efficiency. For polychord weuse 𝑛 live = 500, tolerance = 0.01, 𝑛 repeats = 60, fast fraction = 0 . 1. Considering the BCEmu7 WL-only analysis we find that sampling with polychord estimates a mean value of 𝑆 8 = 0 . 817 + 0 . 025 0 . 025 , consistent with that attained by multinest of 𝑆 8 = 0 . 818 + 0 . 017 0 . 024 . However, in agreement with the findings of DES & KiDS Collaborations (2023), we find that the 68% confidence level for 𝑆 8 attained using multinest is 18%smaller than that estimated with polychord . Similarly we find a WL + kSZ analysis with polychord returns 𝑆 8 = 0 . 821 + 0 . 020 0 . 023 , consistent with the 𝑆 8 = 0 . 823 + 0 . 019 0 . 020 using multinest but with a 9% larger confidence region.', 'AFFILIATIONS': "- 1 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK\n- 2 Kavli Institute for Cosmology Cambridge, Madingley Road, Cambridge, CB3 OHA\n- 3 Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA\n- 4 Institute for Computational Science, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland\n- 5 Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK\n- 6 Centre for Energy, Environmental and Technological Research, Av. Complutense, 40, Moncloa - Aravaca, 28040 Madrid, Spain\n- 7 Kavli Institute for Particle Astrophysics and Cosmology, 382 Via Pueblo Mall Stanford, CA 94305-4060, USA\n- 8 SLAC National Accelerator Laboratory 2575 Sand Hill Road Menlo Park, California 94025, USA\n- 9 Physics Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA\n- 10 Berkeley Center for Cosmological Physics, Department of Physics, University of California, Berkeley, CA 94720, USA\n- 11 Department of Astronomy, Cornell University, Ithaca, NY 14853, USA\n- 12 Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA\n- 13 Kavli Institute for the Physics and Mathematics of the Universe (WPI), UTIAS, The University of Tokyo, Kashiwa, Chiba 2778583, Japan\n- 14 Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15312, USA\n- 15 NSF AI Planning Institute for Physics of the Future, Carnegie Mellon University, Pittsburgh, PA 15213, USA\n- 16 Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA\n- 17 SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA\n- 18 Laborat'orio Interinstitucional de e-Astronomia - LIneA, Rua Gal. Jos'e Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil\n- 19 Observat'orio Nacional, Rua Gal. Jos'e Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil\n- 20 Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA\n- 21 Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s/n, 08193 Barcelona, Spain\n- 22 Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA\n- 23 Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA\n- 24 Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA\n- 25 NASA Goddard Space Flight Center, 8800 Greenbelt Rd, Greenbelt, MD 20771, USA\n- 26 Instituto de F'ısica Gleb Wataghin, Universidade Estadual de Campinas, 13083-859, Campinas, SP, Brazil\n- 27 Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA\n- 28 Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain\n- 29 Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA\n- 30 Universit'e Grenoble Alpes, CNRS, LPSC-IN2P3, 38000  \nGrenoble, France  \n- 31 University Observatory, Faculty of Physics, LudwigMaximilians-Universitat, Scheinerstr. 1, 81679 Munich, Germany 32 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA\n- 33 Brookhaven National Laboratory, Bldg 510, Upton, NY 11973, USA\n- 34 Department of Physics, ETH Zurich, Wolfgang-Pauli-Strasse 16, CH-8093 Zurich, Switzerland\n- 35 Institut de F'ısica d'Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona) Spain\n- 36 Department of Astronomy/Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA\n- 37 Department of Physics, University of Arizona, Tucson, AZ 85721, USA\n- 38 School of Physics and Astronomy, Cardiff University, CF24 3AA, UK\n- 39 Centro de Investigaciones Energ'eticas, Medioambientales y Tecnol'ogicas (CIEMAT), Madrid, Spain\n- 40 Institut de Recherche en Astrophysique et Plan'etologie (IRAP), Universit'e de Toulouse, CNRS, UPS, CNES, 14 Av. Edouard Belin, 31400 Toulouse, France\n- 41 Department of Physics and Astronomy, University of Waterloo, 200 University Ave W, Waterloo, ON N2L 3G1, Canada\n- 42 Institute for Astronomy, University of Edinburgh, Edinburgh EH9 3HJ, UK\n- 43 Department of Physics, Northeastern University, Boston, MA 02115, USA\n- 44 Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA\n- 45 Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK\n- 46 Nordita, KTH Royal Institute of Technology and Stockholm University, Hannes Alfv'ens vag 12, SE-10691 Stockholm, Sweden 47 Physics Department, 2320 Chamberlin Hall, University of WI\n- Wisconsin-Madison, 1150 University Avenue Madison, 53706-1390\n- 48 Department of Physics, University of Genova and INFN, Via Dodecaneso 33, 16146, Genova, Italy\n- 49 Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA\n- 50 Department of Astronomy, University of Illinois at UrbanaChampaign, 1002 W. Green Street, Urbana, IL 61801, USA\n- 51 Department of Physics, Duke University Durham, NC 27708, USA\n- 52 Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 0WA, UK\n- 53 Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK\n- 54 Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA\n- 55 Physics Department, William Jewell College, Liberty, MO, 64068\n- 56 Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794, USA\n- 57 Department of Astronomy, University of Geneva, ch. d' ' Ecogia 16, CH-1290 Versoix, Switzerland\n- 58 Department of Astronomy, University of California, Berkeley, 501 Campbell Hall, Berkeley, CA 94720, USA", '28 L. Bigwood, A. Amon et al.': "- 59 Cerro Tololo Inter-American Observatory, NSFs National Optical-Infrared Astronomy Research Laboratory, Casilla 603, La Serena, Chile\n- 60 Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth, PO1 3FX, UK\n- 61 CNRS, UMR 7095, Institut d'Astrophysique de Paris, F-75014, Paris, France\n- 62 Sorbonne Universit'es, UPMC Univ Paris 06, UMR 7095, Institut d'Astrophysique de Paris, F-75014, Paris, France\n- 63 Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK\n- 64 Institut d'Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain\n- 65 Hamburger Sternwarte, Universitat Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany\n- 66 Department of Physics, IIT Hyderabad, Kandi, Telangana 502285, India\n- 67 Institute of Theoretical Astrophysics, University of Oslo. P.O. Box 1029 Blindern, NO-0315 Oslo, Norway\n- 68 Instituto de Fisica Teorica UAM/CSIC, Universidad Autonoma de Madrid, 28049 Madrid, Spain\n- 69 School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia\n- 70 Santa Cruz Institute for Particle Physics, Santa Cruz, CA 95064, USA\n- 71 Department of Physics, The Ohio State University, Columbus, OH 43210, USA\n- 72 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA\n- 73 Australian Astronomical Optics, Macquarie University, North Ryde, NSW 2113, Australia\n- 74 Lowell Observatory, 1400 Mars Hill Rd, Flagstaff, AZ 86001, USA\n- 75 George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA\n- 76 LPSC Grenoble - 53, Avenue des Martyrs 38026 Grenoble, France\n- 77 Instituci'o Catalana de Recerca i Estudis Avanc¸ats, E-08010 Barcelona, Spain\n- 78 Perimeter Institute for Theoretical Physics, 31 Caroline St. North, Waterloo, ON N2L 2Y5, Canada\n- 79 Ruhr University Bochum, Faculty of Physics and Astronomy, Astronomical Institute, German Centre for Cosmological Lensing, 44780 Bochum, Germany\n- 80 Department of Physics and Astronomy, Pevensey Building, University of Sussex, Brighton, BN1 9QH, UK\n- 81 School of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK 82\n- 83 Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831\n- 84 Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany\n- 85 Universitats-Sternwarte, Fakultat fur Physik, LudwigMaximilians Universitat Munchen, Scheinerstr. 1, 81679 Munchen, Germany"}