Showing 1–22 of 22 results for author: Anbajagane, D

Dark Energy Survey Year 3 results: simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps II. Cosmological results

Authors: M. Gatti, G. Campailla, N. Jeffrey, L. Whiteway, A. Porredon, J. Prat, J. Williamson, M. Raveri, B. Jain, V. Ajani, G. Giannini, M. Yamamoto, C. Zhou, J. Blazek, D. Anbajagane, S. Samuroff, T. Kacprzak, A. Alarcon, A. Amon, K. Bechtol, M. Becker, G. Bernstein, A. Campos, C. Chang, R. Chen , et al. (77 additional authors not shown)

Abstract: We present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years (Y3) of the Dark Energy Survey (DES). We implement: 1) second and third moments; 2) wavelet phase harmonics; 3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven $νw$CDM cosm… ▽ More We present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years (Y3) of the Dark Energy Survey (DES). We implement: 1) second and third moments; 2) wavelet phase harmonics; 3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven $νw$CDM cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the Figure of Merit $\textrm{FoM}(S_8, Ω_{\textrm{m}})$ is improved by 70 percent ($Λ$CDM) and 90 percent ($w$CDM). When all the summary statistics are combined, we achieve a 2 percent constraint on the amplitude of fluctuations parameter $S_8 \equiv σ_8 (Ω_{\textrm{m}}/0.3)^{0.5}$, obtaining $S_8 = 0.794 \pm 0.017$ ($Λ$CDM) and $S_8 = 0.817 \pm 0.021$ ($w$CDM). The constraints from different statistics are shown to be internally consistent (with a $p$-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the DES Y3 data, finding good consistency; we also compare with results from external datasets, such as \planck{} constraints from the Cosmic Microwave Background, finding statistical agreement, with discrepancies no greater than $<2.2σ$. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: 24 pages, 13 figures, to be submitted to PRD. Comments welcome!

Deciphering baryonic feedback with galaxy clusters

Authors: Chun-Hao To, Shivam Pandey, Elisabeth Krause, Nihar Dalal, Dhayaa Anbajagane, David H. Weinberg

Abstract: Upcoming cosmic shear analyses will precisely measure the cosmic matter distribution at low redshifts. At these redshifts, the matter distribution is affected by galaxy formation physics, primarily baryonic feedback from star formation and active galactic nuclei. Employing measurements from the Magneticum and IllustrisTNG simulations and a dark matter + baryon (DMB) halo model, this paper demonstr… ▽ More Upcoming cosmic shear analyses will precisely measure the cosmic matter distribution at low redshifts. At these redshifts, the matter distribution is affected by galaxy formation physics, primarily baryonic feedback from star formation and active galactic nuclei. Employing measurements from the Magneticum and IllustrisTNG simulations and a dark matter + baryon (DMB) halo model, this paper demonstrates that Sunyaev-Zel'dovich (SZ) effect observations of galaxy clusters, whose masses have been calibrated using weak gravitational lensing, can constrain the baryonic impact on cosmic shear with statistical and systematic errors subdominant to the measurement errors of DES-Y3 and LSST-Y1. We further dissect the contributions from different scales and halos with different masses to cosmic shear, highlighting the dominant role of SZ clusters at scales critical for cosmic shear analyses. These findings suggest a promising avenue for future joint analyses of Cosmic Microwave Background (CMB) and lensing surveys. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: To be submitted to JCAP. Comments are welcome!

GODMAX: Modeling gas thermodynamics and matter distribution using JAX

Authors: Shivam Pandey, Jaime Salcido, Chun-Hao To, J. Colin Hill, Dhayaa Anbajagane, Eric J. Baxter, Ian G. McCarthy

Abstract: We introduce GODMAX (Gas thermODynamics and Matter distribution using jAX), a novel code designed to calculate correlations between the cosmological matter distribution and various gas thermodynamic quantities. Utilizing the extensive ANTILLES suite of 200 hydrodynamical simulations with a diverse range of baryonic feedback strengths, we jointly fit the 3D profiles of total matter distribution, el… ▽ More We introduce GODMAX (Gas thermODynamics and Matter distribution using jAX), a novel code designed to calculate correlations between the cosmological matter distribution and various gas thermodynamic quantities. Utilizing the extensive ANTILLES suite of 200 hydrodynamical simulations with a diverse range of baryonic feedback strengths, we jointly fit the 3D profiles of total matter distribution, electron density, and electron pressure across various halo masses and redshifts. By accommodating significant variations in gas profiles expected due to baryonic feedback, solving exact hydrostatic equilibrium equation and offering flexible modeling of non-thermal pressure support, GODMAX has the capability to jointly fit all these profiles within the measurement uncertainties. This advancement enables, for the first time, robust joint analyses of multiple cosmic probes, including the kinetic and thermal Sunyaev-Zel'dovich effect, weak lensing, and X-ray observations. Furthermore, the model accurately captures correlations between the total matter power suppression due to baryonic feedback and local average thermodynamic quantities, such as the baryon fraction and integrated tSZ effect, in high-mass halos, aligning with observations from hydrodynamical simulations. Looking ahead, we forecast the expected constraints on cosmological and baryonic parameters from upcoming weak lensing catalogs from the LSST and tSZ maps from the Simons Observatory. This analysis underscores the importance of cross-correlations between weak lensing and tSZ in enhancing parameter constraints by resolving major systematic uncertainties due to baryonic physics. The GODMAX code leverages the JAX library, resulting in a fully differentiable halo model with native GPU compilation support. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 17 pages, 9 figures

Cosmological Constraints from Combining Galaxy Surveys and Gravitational Wave Observatories

Authors: E. L. Gagnon, D. Anbajagane, J. Prat, C. Chang, J. Frieman

Abstract: Spatial variations in survey properties due to selection effects generate substantial systematic errors in large-scale structure measurements in optical galaxy surveys on very large scales. On such scales, the statistical sensitivity of optical surveys is also limited by their finite sky coverage. By contrast, gravitational wave (GW) sources appear to be relatively free of these issues, provided t… ▽ More Spatial variations in survey properties due to selection effects generate substantial systematic errors in large-scale structure measurements in optical galaxy surveys on very large scales. On such scales, the statistical sensitivity of optical surveys is also limited by their finite sky coverage. By contrast, gravitational wave (GW) sources appear to be relatively free of these issues, provided the angular sensitivity of GW experiments can be accurately characterized. We quantify the expected cosmological information gain from combining the forecast LSST 3$\times$2pt analysis (combination of three 2-point correlations of galaxy density and weak lensing shear fields) with the large-scale auto-correlation of GW sources from proposed next-generation GW experiments. We find that in $Λ$CDM and $w$CDM models, there is no significant improvement in cosmological constraints from combining GW with LSST 3$\times$2pt over LSST alone, due to the large shot noise for the former; however, this combination does enable a $\sim6\%$ constraint on the linear galaxy bias of GW sources. More interestingly, the optical-GW data combination provides tight constraints on models with primordial non-Gaussianity (PNG), due to the predicted scale-dependent bias in PNG models on large scales. Assuming that the largest angular scales that LSST will probe are comparable to those in Stage III surveys ($\ell_{\rm min}\sim50$), the inclusion of next-generation GW measurements could improve constraints on the PNG parameter $f_{\rm NL}$ by up to a factor of $\simeq6.6$ compared to LSST alone, yielding $σ(f_{\rm NL})=8.5$. These results assume the expected capability of a network of Einstein Telescope-like GW observatories, with a detection rate of $10^6$ events/year. We investigate the sensitivity of our results to different assumptions about future GW detectors as well as different LSST analysis choices. △ Less

Submitted 26 December, 2023; originally announced December 2023.

Comments: 16 pages, 8 figures, 5 tables

Subhalos in Galaxy Clusters: Coherent Accretion and Internal Orbits

Authors: Chi Han, Kuan Wang, Camille Avestruz, Dhayaa Anbajagane

Abstract: Subhalo dynamics in galaxy cluster host halos govern the observed distribution and properties of cluster member galaxies. We use the IllustrisTNG simulation to investigate the accretion and orbits of subhalos found in cluster-size halos. We find that the median change in the major axis direction of cluster-size host halos is approximately $80$ degrees between $a\sim0.1$ and present-day. We identif… ▽ More Subhalo dynamics in galaxy cluster host halos govern the observed distribution and properties of cluster member galaxies. We use the IllustrisTNG simulation to investigate the accretion and orbits of subhalos found in cluster-size halos. We find that the median change in the major axis direction of cluster-size host halos is approximately $80$ degrees between $a\sim0.1$ and present-day. We identify coherent regions in the angular distribution of subhalo accretion, and $\sim 68\%$ of accreted subhalos enter their host halo through $\sim38\%$ of the surface area at the virial radius. The majority of galaxy clusters in the sample have $\sim2$ such coherent regions, likely corresponding to cluster-feeding filaments. We further measure angular orbits of subhalos with respect to the host major axis and use a clustering algorithm to identify distinct orbit modes with varying oscillation timescales. The orbit modes correlate with subhalo accretion conditions. Subhalos in orbit modes with shorter oscillations tend to have lower peak masses and accretion directions somewhat more aligned with the major axis. One orbit mode, exhibiting the least oscillatory behavior, largely consists of subhalos that accrete near the plane perpendicular to the host halo major axis. Our findings are consistent with expectations from inflow from major filament structures and internal dynamical friction: most subhalos accrete through filaments, and more massive subhalos experience fewer orbits after accretion. Our work offers a unique quantification of subhalo dynamics that can be connected to how the intracluster medium strips and quenches cluster galaxies. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 10 pages, 8 figures, to be submitted to ApJ, comments welcome

Merger Response of Halo Anisotropy Properties

Authors: Kuan Wang, Philip Mansfield, Dhayaa Anbajagane, Camille Avestruz

Abstract: Anisotropy properties -- halo spin, shape, position offset, velocity offset, and orientation -- are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of haloes. These properties reflect the dynamical state of haloes, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy pro… ▽ More Anisotropy properties -- halo spin, shape, position offset, velocity offset, and orientation -- are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of haloes. These properties reflect the dynamical state of haloes, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy properties in response to merger activity using the IllustrisTNG simulations. We find that the response trajectories of the anisotropy properties significantly deviate from secular evolution. These trajectories have the same qualitative features and timescales across a wide range of merger and host properties. We propose explanations for the behaviour of these properties and connect their evolution to the relevant stages of merger dynamics. We measure the relevant dynamical timescales. We also explore the dependence of the strength of the response on time of merger, merger ratio, and mass of the main halo. These results provide insight into the physics of halo mergers and their effects on the statistical behaviour of halo properties. This study paves the way towards a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of haloes. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 12+3 pages, 5+2 figures. Fig. 4 and 5 are the main figures. To be submitted to MNRAS, comments welcome

Baryonic Imprints on DM Halos: the concentration-mass relation and its dependence on halo and galaxy properties

Authors: Mufan Shao, Dhayaa Anbajagane

Abstract: The halo concentration-mass relation has ubiquitous use in modeling the matter field for cosmological and astrophysical analyses, and including the imprints from galaxy formation physics is tantamount to its robust usage. Many analyses, however, probe the matter around halos selected by a given halo/galaxy property -- rather than by halo mass -- and the imprints under each selection choice can be… ▽ More The halo concentration-mass relation has ubiquitous use in modeling the matter field for cosmological and astrophysical analyses, and including the imprints from galaxy formation physics is tantamount to its robust usage. Many analyses, however, probe the matter around halos selected by a given halo/galaxy property -- rather than by halo mass -- and the imprints under each selection choice can be different. We employ the CAMELS simulation suite to quantify the astrophysics and cosmology dependence of the concentration-mass relation, $c_{\rm vir}-M_{\rm vir}$, when selected on five properties: (i) velocity dispersion, (ii) formation time, (iii) halo spin, (iv) stellar mass, and (v) gas mass. We construct simulation-informed nonlinear models for all properties as a function of halo mass, redshift, and six cosmological/astrophysical parameters, with a mass range $M_{\rm vir} \in [10^{11}, 10^{14.5}] M_\odot/h$. There are many mass-dependent imprints in all halo properties, with clear connections across different properties and non-linear couplings between the parameters. Finally, we extract the $c_{\rm vir}-M_{\rm vir}$ relation for subsamples of halos that have scattered above/below the mean property-$M_{\rm vir}$ relation for a chosen property. Selections on gas mass or stellar mass have a significant impact on the astrophysics/cosmology dependence of $c_{\rm vir}$, while those on any of the other three properties have a significant (mild) impact on the cosmology (astrophysics) dependence. We show that ignoring such selection effects can lead to errors of $\approx 25\%$ in baryon imprint modelling of $c_{\rm vir}$. Our nonlinear model for all properties is made publicly available. △ Less

Submitted 23 April, 2024; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: V2: Version accepted to Theoj

Dark Energy Survey Year 3 results: simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps I: validation on simulations

Authors: M. Gatti, N. Jeffrey, L. Whiteway, J. Williamson, B. Jain, V. Ajani, D. Anbajagane, G. Giannini, C. Zhou, A. Porredon, J. Prat, M. Yamamoto, J. Blazek, T. Kacprzak, S. Samuroff, A. Alarcon, A. Amon, K. Bechtol, M. Becker, G. Bernstein, A. Campos, C. Chang, R. Chen, A. Choi, C. Davis , et al. (76 additional authors not shown)

Abstract: Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end simulation-based analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogues of the Dark Energy Survey. We implement: 1) second and third moments; 2) wavelet phase… ▽ More Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end simulation-based analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogues of the Dark Energy Survey. We implement: 1) second and third moments; 2) wavelet phase harmonics (WPH); 3) the scattering transform (ST). Our analysis is fully based on simulations, it spans a space of seven $νw$CDM cosmological parameters, and it forward models the most relevant sources of systematics of the data (masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration). We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a likelihood-free-inference approach. We validate the pipeline extensively, and we find that WPH exhibits the strongest performance when combined with second moments, followed by ST. and then by third moments. The combination of all the different statistics further enhances constraints with respect to second moments, up to 25 per cent, 15 per cent, and 90 per cent for $S_8$, $Ω_{\rm m}$, and the Figure-Of-Merit ${\rm FoM_{S_8,Ω_{\rm m}}}$, respectively. We further find that non-Gaussian statistics improve constraints on $w$ and on the amplitude of intrinsic alignment with respect to second moments constraints. The methodological advances presented here are suitable for application to Stage IV surveys from Euclid, Rubin-LSST, and Roman with additional validation on mock catalogues for each survey. In a companion paper we present an application to DES Year 3 data. △ Less

Submitted 4 November, 2023; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 25 pages, 18 figures. Comments welcome!

Primordial non-Gaussianities with weak lensing: Information on non-linear scales in the Ulagam full-sky simulations

Authors: Dhayaa Anbajagane, Chihway Chang, Hayden Lee, Marco Gatti

Abstract: Primordial non-Gaussianities (PNGs) are signatures in the density field that encode particle physics processes from the inflationary epoch. Such signatures have been extensively studied using the Cosmic Microwave Background, through constraining the amplitudes, $f^{X}_{\rm NL}$, with future improvements expected from large-scale structure surveys; specifically, the galaxy correlation functions. We… ▽ More Primordial non-Gaussianities (PNGs) are signatures in the density field that encode particle physics processes from the inflationary epoch. Such signatures have been extensively studied using the Cosmic Microwave Background, through constraining the amplitudes, $f^{X}_{\rm NL}$, with future improvements expected from large-scale structure surveys; specifically, the galaxy correlation functions. We show that weak lensing fields can be used to achieve competitive and complementary constraints. This is shown via the new Ulagam suite of N-body simulations, a subset of which evolves primordial fields with four types of PNGs. We create full-sky lensing maps and estimate the Fisher information from three summary statistics measured on the maps: the moments, the cumulative distribution function, and the 3-point correlation function. We find that the year 10 sample from the Rubin Observatory Legacy Survey of Space and Time (LSST) can constrain PNGs to $σ(f^{\rm\,eq}_{\rm NL}) \approx 110$, $σ(f^{\rm\,or,lss}_{\rm NL}) \approx 120$, $σ(f^{\rm\,loc}_{\rm NL}) \approx 40$. For the former two, this is better than or comparable to expected galaxy clustering-based constraints from the Dark Energy Spectroscopic Instrument (DESI). The PNG information in lensing fields is on non-linear scales and at low redshifts ($z \lesssim 1.25$), with a clear origin in the evolution history of massive halos. The constraining power degrades by $\sim\!\!60\%$ under scale cuts of $\gtrsim 20{\,\rm Mpc}$, showing there is still significant information on scales mostly insensitive to small-scale systematic effects (e.g. baryons). We publicly release the Ulagam suite to enable more survey-focused analyses. △ Less

Submitted 29 March, 2024; v1 submitted 3 October, 2023; originally announced October 2023.

Comments: [v2]: Version accepted to JCAP

Cosmological shocks around galaxy clusters: A coherent investigation with DES, SPT & ACT

Authors: D. Anbajagane, C. Chang, E. J. Baxter, S. Charney, M. Lokken, M. Aguena, S. Allam, O. Alves, A. Amon, R. An, F. Andrade-Oliveira, D. Bacon, N. Battaglia, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, L. Bleem, S. Bocquet, J. R. Bond, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, R. Chen, A. Choi , et al. (89 additional authors not shown)

Abstract: We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around… ▽ More We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around $10^5$ clusters with mass and redshift ranges $10^{13.7} < M_{\rm 200m}/M_\odot < 10^{15.5}$ and $0.1 < z < 2$, and the total sky coverage of the maps is $\approx 15,000 \,\,{\rm deg}^2$. We find a clear pressure deficit at $R/R_{\rm 200m}\approx 1.1$ in SZ profiles around both ACT and SPT clusters, estimated at $6σ$ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2-5 compared to predictions, and we discuss possible interpretations of this behavior. An oriented stacking of clusters -- where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogs -- shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius. △ Less

Submitted 12 December, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: [v2]: Version accepted to MNRAS

Beyond the 3rd moment: A practical study of using lensing convergence CDFs for cosmology with DES Y3

Authors: D. Anbajagane, C. Chang, A. Banerjee, T. Abel, M. Gatti, V. Ajani, A. Alarcon, A. Amon, E. J. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Chen, A. Choi, C. Davis, J. DeRose, H. T. Diehl, S. Dodelson, C. Doux, A. Drlica-Wagner, K. Eckert, J. Elvin-Poole , et al. (73 additional authors not shown)

Abstract: Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at… ▽ More Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at multiple scales as a summary of the galaxy lensing convergence field. Using a suite of N-body lightcone simulations, we show the CDFs' constraining power is modestly better than that of the 2nd and 3rd moments of the field, as they approximately capture the information from all moments of the field in a concise data vector. We then study the practical aspects of applying the CDFs to observational data, using the first three years of the Dark Energy Survey (DES Y3) data as an example, and compute the impact of different systematics on the CDFs. The contributions from the point spread function are 2-3 orders of magnitude below the cosmological signal, while those from reduced shear approximation contribute $\lesssim 1\%$ to the signal. Source clustering effects and baryon imprints contribute $1-10\%$. Enforcing scale cuts to limit systematics-driven biases in parameter constraints degrades these constraints a noticeable amount, and this degradation is similar for the CDFs and the moments. We also detect correlations between the observed convergence field and the shape noise field at $13σ$. We find that the non-Gaussian correlations in the noise field must be modeled accurately to use the CDFs, or other statistics sensitive to all moments, as a rigorous cosmology tool. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 21 pages, 12 figures

Baryonic Imprints on DM Halos: The concentration-mass relation in the CAMELS simulations

Authors: Mufan Shao, Dhayaa Anbajagane, Chihway Chang

Abstract: The physics of baryons in halos, and their subsequent influence on the total matter phase space, has a rich phenomenology and must be well understood in order to pursue a vast set of questions in both cosmology and astrophysics. We use the CAMELS simulation suite to quantify the impact of four different galaxy formation parameters/processes (as well as two cosmological parameters) on the concentra… ▽ More The physics of baryons in halos, and their subsequent influence on the total matter phase space, has a rich phenomenology and must be well understood in order to pursue a vast set of questions in both cosmology and astrophysics. We use the CAMELS simulation suite to quantify the impact of four different galaxy formation parameters/processes (as well as two cosmological parameters) on the concentration-mass relation, $c_{\rm vir} - M_{\rm vir}$. We construct a simulation-informed nonlinear model for concentration as a function of halo mass, redshift, and 6 cosmological/astrophysical parameters. This is done for two galaxy formation models, IllustrisTNG and SIMBA, using 1000 simulations of each. We extract the imprints of galaxy formation across a wide range in mass $M_{\rm vir} \in [10^{11}, 10^{14.5}] M_{\rm \odot}/h$ and in redshift $z \in [0,6]$ finding many strong mass- and redshift-dependent features. Comparisons between the IllustrisTNG and SIMBA results show the astrophysical model choices cause significant differences in the mass and redshift dependence of these baryon imprints. Finally, we use existing observational measurements of $c_{\rm vir} - M_{\rm vir}$ to provide rough limits on the four astrophysical parameters. Our nonlinear model is made publicly available and can be used to include CAMELS-based baryon imprints in any halo model-based analysis. △ Less

Submitted 17 August, 2023; v1 submitted 12 December, 2022; originally announced December 2022.

Comments: [v1]: 8 figures, 14 pages [v2]: version accepted in MNRAS

A multi-simulation study of relativistic SZ temperature scalings in galaxy clusters and groups

Authors: Elizabeth Lee, Dhayaa Anbajagane, Priyanka Singh, Jens Chluba, Daisuke Nagai, Scott T. Kay, Weiguang Cui, Klaus Dolag, Gustavo Yepes

Abstract: The Sunyaev-Zeldovich (SZ) effect is a powerful tool in modern cosmology. With future observations promising ever improving SZ measurements, the relativistic corrections to the SZ signals from galaxy groups and clusters are increasingly relevant. As such, it is important to understand the differences between three temperature measures: (a) the average relativistic SZ (rSZ) temperature, (b) the mas… ▽ More The Sunyaev-Zeldovich (SZ) effect is a powerful tool in modern cosmology. With future observations promising ever improving SZ measurements, the relativistic corrections to the SZ signals from galaxy groups and clusters are increasingly relevant. As such, it is important to understand the differences between three temperature measures: (a) the average relativistic SZ (rSZ) temperature, (b) the mass-weighted temperature relevant for the thermal SZ (tSZ) effect, and (c) the X-ray spectroscopic temperature. In this work, we compare these cluster temperatures, as predicted by the {\sc Bahamas} \& {\sc Macsis}, {\sc Illustris-TNG}, {\sc Magneticum}, and {\sc The Three Hundred Project} simulations. Despite the wide range of simulation parameters, we find the SZ temperatures are consistent across the simulations. We estimate a $\simeq 10\%$ level correction from rSZ to clusters with $Y\simeq10^{-4}$~Mpc$^{-2}$. Our analysis confirms a systematic offset between the three temperature measures; with the rSZ temperature $\simeq 20\%$ larger than the other measures, and diverging further at higher redshifts. We demonstrate that these measures depart from simple self-similar evolution and explore how they vary with the defined radius of haloes. We investigate how different feedback prescriptions and resolution affect the observed temperatures, and discover the SZ temperatures are rather insensitive to these details. The agreement between simulations indicates an exciting avenue for observational and theoretical exploration, determining the extent of relativistic SZ corrections. We provide multiple simulation-based fits to the scaling relations for use in future SZ modelling. △ Less

Submitted 26 September, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: Accepted by MNRAS

Correlations of Dark Matter, Gas and Stellar Profiles in Dark Matter Halos

Authors: Arya Farahi, Daisuke Nagai, Dhayaa Anbajagane

Abstract: Halos of similar mass and redshift exhibit a large degree of variability in their differential properties, such as dark matter, hot gas, and stellar mass density profiles. This variability is an indicator of diversity in the formation history of these dark matter halos that is reflected in the coupling of scatters about the mean relations. In this work, we show that the strength of this coupling d… ▽ More Halos of similar mass and redshift exhibit a large degree of variability in their differential properties, such as dark matter, hot gas, and stellar mass density profiles. This variability is an indicator of diversity in the formation history of these dark matter halos that is reflected in the coupling of scatters about the mean relations. In this work, we show that the strength of this coupling depends on the scale at which halo profiles are measured. By analyzing the outputs of the IllustrisTNG hydrodynamical cosmological simulations we report the radial- and mass-dependent couplings between the dark matter, hot gas, and stellar mass radial density profiles utilizing the population diversity in dark matter halos. We find that for the same mass halos the scatters in the density of baryons and dark matter are strongly coupled at large scales ($r>R_{200}$); but the coupling between gas and dark matter density profiles fades near the core of halos ($r < 0.3 R_{200}$). We then show that the correlation between halo profile and integrated quantities induces a radius-dependent additive bias in the profile observables of halos when halos are selected on properties other than their mass. We discuss the impact of this effect on cluster abundance and cross-correlations cosmology with multi-wavelength cosmological surveys. △ Less

Submitted 31 May, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: 8 figures, 4 appendices, data in Table 2 (Accepted)

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

Authors: Weiguang Cui, Romeel Dave, Alexander Knebe, Elena Rasia, Meghan Gray, Frazer Pearce, Chris Power, Gustavo Yepes, Dhayaa Anbajagane, Daniel Ceverino, Ana Contreras-Santos, Daniel de Andres, Marco De Petris, Stefano Ettori, Roan Haggar, Qingyang Li, Yang Wang, Xiaohu Yang, Stefano Borgani, Klaus Dolag, Ying Zu, Ulrike Kuchner, Rodrigo Cañas, Antonio Ferragamo, Giulia Gianfagna

Abstract: We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art ga… ▽ More We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments. △ Less

Submitted 31 May, 2022; v1 submitted 28 February, 2022; originally announced February 2022.

Comments: 20 pages, 12 figures, matched to published version in MNRAS

The scatter in the galaxy-halo connection: a machine learning analysis

Authors: Richard Stiskalek, Deaglan J. Bartlett, Harry Desmond, Dhayaa Anbajagane

Abstract: We apply machine learning, a powerful method for uncovering complex correlations in high-dimensional data, to the galaxy-halo connection of cosmological hydrodynamical simulations. The mapping between galaxy and halo variables is stochastic in the absence of perfect information, but conventional machine learning models are deterministic and hence cannot capture its intrinsic scatter. To overcome t… ▽ More We apply machine learning, a powerful method for uncovering complex correlations in high-dimensional data, to the galaxy-halo connection of cosmological hydrodynamical simulations. The mapping between galaxy and halo variables is stochastic in the absence of perfect information, but conventional machine learning models are deterministic and hence cannot capture its intrinsic scatter. To overcome this limitation, we design an ensemble of neural networks with a Gaussian loss function that predict probability distributions, allowing us to model statistical uncertainties in the galaxy-halo connection as well as its best-fit trends. We extract a number of galaxy and halo variables from the Horizon-AGN and IllustrisTNG100-1 simulations and quantify the extent to which knowledge of some subset of one enables prediction of the other. This allows us to identify the key features of the galaxy-halo connection and investigate the origin of its scatter in various projections. We find that while halo properties beyond mass account for up to 50 per cent of the scatter in the halo-to-stellar mass relation, the prediction of stellar half-mass radius or total gas mass is not substantially improved by adding further halo properties. We also use these results to investigate semi-analytic models for galaxy size in the two simulations, finding that assumptions relating galaxy size to halo size or spin are not successful. △ Less

Submitted 8 June, 2022; v1 submitted 28 February, 2022; originally announced February 2022.

Comments: 20 pages, 11 figures. Accepted in MNRAS

KLLR: A scale-dependent, multivariate model class for regression analysis

Authors: Arya Farahi, Dhayaa Anbajagane, August Evrard

Abstract: The underlying physics of astronomical systems governs the relation between their measurable properties. Consequently, quantifying the statistical relationships between system-level observable properties of a population offers insights into the astrophysical drivers of that class of systems. While purely linear models capture behavior over a limited range of system scale, the fact that astrophysic… ▽ More The underlying physics of astronomical systems governs the relation between their measurable properties. Consequently, quantifying the statistical relationships between system-level observable properties of a population offers insights into the astrophysical drivers of that class of systems. While purely linear models capture behavior over a limited range of system scale, the fact that astrophysics is ultimately scale-dependent implies the need for a more flexible approach to describing population statistics over a wide dynamic range. For such applications, we introduce and implement a class of Kernel-Localized Linear Regression (KLLR) models. KLLR is a natural extension to the commonly-used linear models that allows the parameters of the linear model -- normalization, slope, and covariance matrix -- to be scale-dependent. KLLR performs inference in two steps: (1) it estimates the mean relation between a set of independent variables and a dependent variable and; (2) it estimates the conditional covariance of the dependent variables given a set of independent variables. We demonstrate the model's performance in a simulated setting and showcase an application of the proposed model in analyzing the baryonic content of dark matter halos. As a part of this work, we publicly release a Python implementation of the KLLR method. △ Less

Submitted 26 April, 2022; v1 submitted 20 February, 2022; originally announced February 2022.

Comments: The code is publicly available at https://github.com/afarahi/kllr and can be installed through `pip install kllr`

Shocks in the Stacked Sunyaev-Zel'dovich Profiles of Clusters II: Measurements from SPT-SZ + Planck Compton-y Map

Authors: D. Anbajagane, C. Chang, B. Jain, S. Adhikari, E. J. Baxter, B. A. Benson, L. E. Bleem, S. Bocquet, M. S. Calzadilla, J. E. Carlstrom, C. L. Chang, R. Chown, T. M. Crawford, A. T. Crites, W. Cui, T. de Haan, L. Di Mascolo, M. A. Dobbs, W. B. Everett, E. M. George, S. Grandis, N. W. Halverson, G. P. Holder, W. L. Holzapfel, J. D. Hrubes , et al. (21 additional authors not shown)

Abstract: We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg^2 SPT-SZ survey on the locations of clusters identified in that same dataset. The sample contains 516 clusters with mean mass <M200m> = 1e14.9 msol and redshift <z> = 0.55. We analyz… ▽ More We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg^2 SPT-SZ survey on the locations of clusters identified in that same dataset. The sample contains 516 clusters with mean mass <M200m> = 1e14.9 msol and redshift <z> = 0.55. We analyze in parallel a set of zoom-in hydrodynamical simulations from The Three Hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = $1.08 \pm 0.09$, measured at $3.1σ$ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = $4.58 \pm 1.24$ at $2.0σ$ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyze clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius (Rsh_acc) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, Rsh_acc/Rsp > $2.16 \pm 0.59$. △ Less

Submitted 16 May, 2022; v1 submitted 8 November, 2021; originally announced November 2021.

Comments: [v1]: 8 Figures, 16 Pages in Main text. [v2]: Added text to discussion. Version accepted in MNRAS

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Authors: Dhayaa Anbajagane, Han Aung, August E. Evrard, Arya Farahi, Daisuke Nagai, David J. Barnes, Weiguang Cui, Klaus Dolag, Ian G. McCarthy, Elena Rasia, Gustavo Yepes

Abstract: Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the rat… ▽ More Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available. △ Less

Submitted 21 December, 2021; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: [v1]: 8 figures, 14 pages in main text + 5 figures, 3 pages in Appendix. [v2]: Minor text edits. Accepted in MNRAS

Baryonic Imprints on DM Halos: Population Statistics from Dwarf Galaxies to Galaxy Clusters

Authors: Dhayaa Anbajagane, August E. Evrard, Arya Farahi

Abstract: In a purely cold dark matter universe, the initial matter power spectrum and its subsequent gravitational growth contain no special mass- or time-scales, and so neither do the emergent population statistics of internal dark matter (DM) halo properties. Using 1.5 million halos from three IllustrisTNG realizations of a LambdaCDM universe, we show that galaxy formation physics drives non-monotonic fe… ▽ More In a purely cold dark matter universe, the initial matter power spectrum and its subsequent gravitational growth contain no special mass- or time-scales, and so neither do the emergent population statistics of internal dark matter (DM) halo properties. Using 1.5 million halos from three IllustrisTNG realizations of a LambdaCDM universe, we show that galaxy formation physics drives non-monotonic features ("wiggles") into DM property statistics across six decades in halo mass, from dwarf galaxies to galaxy clusters. We characterize these features by extracting the halo mass-dependent statistics of five DM halo properties -- velocity dispersion, NFW concentration, density- and velocity-space shapes, and formation time -- using kernel-localized linear regression (KLLR). Comparing precise estimates of normalizations, slopes, and covariances between realizations with and without galaxy formation, we find systematic deviations across all mass-scales, with maximum deviations of 25% at the Milky-Way mass of 1e12 Msun. The mass-dependence of the wiggles is set by the interplay between different cooling and feedback mechanisms, and we discuss its observational implications. The property covariances depend strongly on halo mass and physics treatment, but the correlations are mostly robust. Using multivariate KLLR and interpretable machine learning, we show the halo concentration and velocity-space shape are principal contributors, at different mass, to the velocity dispersion variance. Statistics of mass accretion rate and DM surface pressure energy are provided in an appendix. We publicly release halo property catalogs and KLLR parameters for the TNG runs at twenty epochs up to z = 12. △ Less

Submitted 20 November, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: [v1]: 10 figures, 16 pages in Main text, and 5 figures, 5 pages in Appendix. [v2]: Minor text edit. Accepted in MNRAS

A Test of the Standard Cosmological Model with Geometry and Growth

Authors: Uendert Andrade, Dhayaa Anbajagane, Rodrigo von Marttens, Dragan Huterer, Jailson Alcaniz

Abstract: We perform a general test of the $Λ{\rm CDM}$ and $w {\rm CDM}$ cosmological models by comparing constraints on the geometry of the expansion history to those on the growth of structure. Specifically, we split the total matter energy density, $Ω_M$, and (for $w {\rm CDM}$) dark energy equation of state, $w$, into two parameters each: one that captures the geometry, and another that captures the gr… ▽ More We perform a general test of the $Λ{\rm CDM}$ and $w {\rm CDM}$ cosmological models by comparing constraints on the geometry of the expansion history to those on the growth of structure. Specifically, we split the total matter energy density, $Ω_M$, and (for $w {\rm CDM}$) dark energy equation of state, $w$, into two parameters each: one that captures the geometry, and another that captures the growth. We constrain our split models using current cosmological data, including type Ia supernovae, baryon acoustic oscillations, redshift space distortions, gravitational lensing, and cosmic microwave background (CMB) anisotropies. We focus on two tasks: (i) constraining deviations from the standard model, captured by the parameters $ΔΩ_M \equiv Ω_M^{\rm grow}-Ω_M^{\rm geom}$ and $Δw \equiv w^{\rm grow}-w^{\rm geom}$, and (ii) investigating whether the $S_8$ tension between the CMB and weak lensing can be translated into a tension between geometry and growth, i.e. $ΔΩ_M \neq 0$, $Δw \neq 0$. In both the split $Λ{\rm CDM}$ and $w {\rm CDM}$ cases, our results from combining all data are consistent with $ΔΩ_M = 0$ and $Δw = 0$. If we omit BAO/RSD data and constrain the split $w {\rm CDM}$ cosmology, we find the data prefers $Δw<0$ at $3.6σ$ significance and $ΔΩ_M>0$ at $4.2σ$ evidence. We also find that for both CMB and weak lensing, $ΔΩ_M$ and $S_8$ are correlated, with CMB showing a slightly stronger correlation. The general broadening of the contours in our extended model does alleviate the $S_8$ tension, but the allowed nonzero values of $ΔΩ_M$ do not encompass the $S_8$ values that would point toward a mismatch between geometry and growth as the origin of the tension. △ Less

Submitted 24 November, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 27 pages, 10 figures. References updated; matches version published in JCAP

Journal ref: JCAP 11 (2021) 014

Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

Authors: Dhayaa Anbajagane, August E. Evrard, Arya Farahi, David J. Barnes, Klaus Dolag, Ian G. McCarthy, Dylan Nelson, Annalisa Pillepich

Abstract: We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generate… ▽ More We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $γ= -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $γ\simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 17 pages, 15 figures