Publications

The circumgalactic medium (CGM) in $\gtrsim 10^{12}\ M_\odot$ halos is dominated by a hot phase ($T \gtrsim 10^{6}$ K). While many models exist for the hot gas structure, there is as yet no consensus. We compare cooling flow models, in which the hot CGM flows inward due to radiative cooling, to the CGM of $\sim 10^{12}-10^{13}\ M_\odot$ halos in galaxy formation simulations from the FIRE project at $z\sim0$. The simulations include realistic cosmological evolution and feedback from stars but neglect AGN feedback. At both mass scales, CGM inflows are typically dominated by the hot phase rather than by the ‘precipitation’ of cold gas. Despite being highly idealized, we find that cooling flows describe $\sim 10^{13}\ M_\odot$ halos very well, with median agreement in the density and temperature profiles of $\sim 20\%$ and $\sim 10\%$, respectively. This indicates that stellar feedback has little impact on CGM scales in those halos. For $\sim 10^{12}\ M_\odot$ halos, the thermodynamic profiles are also accurately reproduced in the outer CGM. For some of these lower-mass halos, cooling flows significantly overpredict the hot gas density in the inner CGM. This could be due to multidimensional angular momentum effects not well captured by our 1D cooling flow models and/or to the larger cold gas fractions in these regions. Turbulence, which contributes $\sim 10-40\%$ of the total pressure, must be included to accurately reproduce the temperature profiles. Overall, cooling flows predict entropy profiles in better agreement with the FIRE simulations than other idealized models in the literature.

Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. Active galactic nuclei (AGN) feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of Feedback In Realistic Environments (FIRE)-3 cosmological zoom-in simulations with AGN feedback evolved to $z \sim 0$, examining the impact of AGN feedback on a set of galaxies with halos in the mass range $10^{12}–10^{13}\ M_\odot$. These simulations combine detailed stellar and ISM physics with multichannel AGN feedback including radiative feedback, mechanical outflows, and, in some simulations, cosmic rays (CRs). We find that massive ($> L_\ast$) galaxies in these simulations can match local scaling relations including the stellar mass–halo mass relation and the $M_{\mathrm{BH}}–\sigma$ relation; in the stronger model with CRs, they also match the size–mass relation and the Faber–Jackson relation. Many of the massive galaxies in the simulations with AGN feedback have quenched star formation and elliptical morphologies, in qualitative agreement with observations. In contrast, simulations at the massive end without AGN feedback produce galaxies that are too massive and form stars too rapidly, are order-of-magnitude too compact, and have velocity dispersions well above Faber–Jackson. Despite these successes, the AGN models analyzed do not produce uniformly realistic galaxies when the feedback parameters are held constant: While the stronger model produces the most realistic massive galaxies, it tends to overquench the lower-mass galaxies. This indicates that further refinements of the AGN modeling are needed.

The properties of warm-hot gas around $\sim L_\ast$ galaxies can be studied with absorption lines from highly ionized metals. We predict Ne VIII column densities from cosmological zoom-in simulations of halos with masses in $\sim 10^{12}$ and $\sim 10^{13}\ M_\odot$ from the Feedback in Realistic Environments (FIRE) project. Ne VIII traces the volume-filling, virial-temperature gas in $\sim 10^{12}\ M_\odot$ halos. In $\sim 10^{13}\ M_\odot$ halos the Ne VIII gas 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 at $M_* \gtrsim 10^{10.5}\ M_\odot$. Median Ne VIII columns 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 Ne VIII profiles 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) $\sim 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 Ne VIII columns. The circumgalactic Ne VIII observations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties.

In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000$h^{-1}$Mpc)$^3$ and evolves close to two trillion particles, corresponding to a mass resolution of $m_p\sim 4.6\cdot 10^7 h^{-1}$M$_\odot$. These specifications enable comprehensive investigations of the galaxy-halo connection, capturing halos down to small masses. Further, the large volume resolves scales typical of modern surveys with good statistical coverage of high mass halos. The simulation was carried out on the GPU-accelerated system Summit, one of the fastest supercomputers currently available. We provide specifics about the Farpoint run and present an initial set of results. The high mass resolution facilitates precise measurements of important global statistics, such as the halo concentration-mass relation and the correlation function down to small scales. Selected subsets of the simulation data products are publicly available via the HACC Simulation Data Portal.

This paper introduces Subhalo Mass-loss Analysis using Core Catalogs (SMACC). SMACC adds a mass model to substructure merger trees based on halo “core tracking.” Our approach avoids the need for running expensive subhalo finding algorithms and instead uses subhalo mass-loss modeling to assign masses to halo cores. We present details of the SMACC methodology and demonstrate its excellent performance in describing halo substructure and its evolution. Validation of the approach is carried out using cosmological simulations at significantly different resolutions. We apply SMACC to the 1.24-trillion-particle Last Journey simulation and construct core catalogs with the additional mass information. These catalogs can be readily used as input to semianalytic models or subhalo abundance matching approaches to determine approximate galaxy distributions, as well as for in-depth studies of small-scale structure evolution.

The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have been run in situ to enable a diverse set of science projects, and at the same time, to keep the resulting data amount manageable. Analysis outputs have been generated starting at redshift $z \sim 10$ to allow for construction of synthetic galaxy catalogs using a semi-analytic modeling approach in post-processing. As part of our in situ analysis pipeline we employ a new method for tracking halo sub-structures, introducing the concept of subhalo cores. The production of multi-wavelength synthetic sky maps is facilitated by generating particle lightcones in situ, also beginning at $z \sim 10$. We provide an overview of the simulation set-up and the generated data products; a first set of analysis results is presented. A subset of the data is publicly available.