skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Friday, February 6 until 10:00 AM ET on Saturday, February 7 due to maintenance. We apologize for the inconvenience.


Title: Zooming by in the CARPoolGP Lane: New CAMELS-TNG Simulations of Zoomed-in Massive Halos
Abstract Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with nontrivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, in particular for halos in the high-mass end of the mass function. In this work, we develop a novel sampling and reduced variance regression method,CARPoolGP, which leverages built-in correlations between samples in different locations of high dimensional parameter spaces to provide an efficient way to explore parameter space and generate low-variance emulations of summary statistics. We use this method to extend the Cosmology and Astrophysics with machinE Learning Simulations to include a set of 768 zoom-in simulations of halos in the mass range of 1013–1014.5Mh−1that span a 28-dimensional parameter space in the IllustrisTNG model. With these simulations and the CARPoolGP emulation method, we explore parameter trends in the ComptonY–M, black hole mass–halo mass, and metallicity–mass relations, as well as thermodynamic profiles and quenched fractions of satellite galaxies. We use these emulations to provide a physical picture of the complex interplay between supernova and active galactic nuclei feedback. We then use emulations of theY–Mrelation of massive halos to perform Fisher forecasts on astrophysical parameters for future Sunyaev–Zeldovich observations and find a significant improvement in forecasted constraints. We publicly release both the simulation suite and CARPoolGP software package.  more » « less
Award ID(s):
2108678 2108944 2307419
PAR ID:
10524237
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
OIP Science
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
968
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
11
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The properties of warm-hot gas around ∼L*galaxies can be studied with absorption lines from highly ionized metals. We predict Neviiicolumn densities from cosmological zoom-in simulations of halos with masses in ∼1012and ∼1013Mfrom the Feedback in Realistic Environments (FIRE) project. Neviiitraces the volume-filling, virial-temperature gas in ∼1012Mhalos. In ∼1013Mhalos the Neviiigas is clumpier, and biased toward the cooler part of the warm-hot phase. We compare the simulations to observations from the COS Absorption Survey of Baryon Harbors (or CASBaH) and COS Ultraviolet Baryon Survey (or CUBS). We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass–halo-mass relation, especially atM≳ 1010.5M. Median Neviiicolumns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Neviiiprofiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) ∼0.1, indicating that observations are consistent with plausible circumgalactic medium temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or active galactic nuclei feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Neviiicolumns. The circumgalactic Neviiiobservations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties. 
    more » « less
  2. Abstract We present a study on the inference of cosmological and astrophysical parameters using stacked galaxy cluster profiles. Utilizing the CAMELS-zoomGZ simulations, we explore how various cluster properties—such as X-ray surface brightness, gas density, temperature, metallicity, and Compton-y profiles—can be used to predict parameters within the 28-dimensional parameter space of the IllustrisTNG model. Through neural networks, we achieve a high correlation coefficient of 0.97 or above for all cosmological parameters, including Ωm,H0, andσ8, and over 0.90 for the remaining astrophysical parameters, showcasing the effectiveness of these profiles for parameter inference. We investigate the impact of different radial cuts, with bins ranging from 0.1R200cto 0.7R200c, to simulate current observational constraints. Additionally, we perform a noise sensitivity analysis, adding up to 40% Gaussian noise (corresponding to signal-to-noise ratios as low as 2.5), revealing that key parameters such as Ωm,H0, and the initial mass function slope remain robust even under extreme noise conditions. We also compare the performance of full radial profiles against integrated quantities, finding that profiles generally lead to more accurate parameter inferences. Our results demonstrate that stacked galaxy cluster profiles contain crucial information on both astrophysical processes within groups and clusters and the underlying cosmology of the Universe. This underscores their significance for interpreting the complex data expected from next-generation surveys and reveals, for the first time, their potential as a powerful tool for parameter inference. 
    more » « less
  3. Abstract We discover analytic equations that can infer the value of Ωmfrom the positions and velocity moduli of halo and galaxy catalogs. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from GadgetN-body simulations to perform field-level likelihood-free inference, and show that our model can infer Ωmwith ∼6% accuracy from halo catalogs of thousands ofN-body simulations run with six different codes: Abacus, CUBEP3M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict Ωmfrom halo catalogs of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that, by tuning a single free parameter, our equations can also infer the value of Ωmfrom galaxy catalogs of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogs from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and Ωm, one that is not affected by galaxy formation physics down to scales as small as 10h−1kpc. 
    more » « less
  4. Abstract As cosmic microwave background (CMB) photons traverse the universe, anisotropies can be induced via Thomson scattering (proportional to the electron density; optical depth) and inverse Compton scattering (proportional to the electron pressure; thermal Sunyaev–Zel’dovich effect). Measurements of anisotropy in optical depthτand Comptonyparameters are imprinted by the galaxies and galaxy clusters and are thus sensitive to the thermodynamic properties of the circumgalactic medium and intergalactic medium. We use an analytic halo model to predict the power spectrum of the optical depth (ττ), the cross-correlation between the optical depth and the Comptonyparameter (τy), and the cross-correlation between the optical depth and galaxy clustering (τg), and compare this model to cosmological simulations. We constrain the optical depths of halos atz≲ 3 using a technique originally devised to constrain patchy reionization at a higher redshift range. The forecasted signal-to-noise ratio is 2.6, 8.5, and 13, respectively, for a CMB-S4-like experiment and a Vera C. Rubin Observatory–like optical survey. We show that a joint analysis of these probes can constrain the amplitude of the density profiles of halos to 6.5% and the pressure profiles to 13%. These constraints translate to astrophysical parameters, such as the gas mass fraction,fg, which can be constrained to 5.3% uncertainty atz∼ 0. The cross-correlations presented here are complementary to other CMB and galaxy cross-correlations since they do not require spectroscopic galaxy redshifts and are another example of how such correlations are a powerful probe of the astrophysics of galaxy evolution. 
    more » « less
  5. Abstract The shape and orientation of dark matter (DM) halos are sensitive to the microphysics of the DM particles, yet in many mass models, the symmetry axes of the Milky Way’s DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well motivated for the inner DM halo, but not for the outer halo. We use zoomed-in cosmological baryonic simulations from the Latte suite of FIRE-2 Milky Way–mass galaxies to explore the evolution of the DM halo’s orientation with radius and time, with or without a major merger with a Large Magellanic Cloud analog, and when varying the DM model. In three of the four cold DM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20° beyond about 30 galactocentric kpc, reaching a maximum of 30°–90°, depending on the individual halo’s formation history. In identical simulations using a model of self-interacting DM withσ= 1 cm2g−1, the halo remains aligned with the stellar disk out to ∼200–400 kpc. Interactions with massive satellites (M≳ 4 × 1010Mat pericenter;M≳ 3.3 × 1010Mat infall) affect the orientation of the halo significantly, aligning the halo’s major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of self-interacting DM, if the effects of massive satellites can be accounted for. 
    more » « less