Search

Over the next decade, the astronomical community will be commissioning multiple wide-field observatories well suited for studying stellar halos in both integrated light and resolved stars. In preparation for this, we use five high-resolution cosmological simulations of Milky Way–like galaxies from the FOGGIE suite to explore the properties and components of stellar halos. These simulations are run with high time (5 Myr) and stellar mass (1000 M _⊙ ) resolution to better model the properties and origins of low-density regions like stellar halos. We find that the FOGGIE stellar halos have masses, metallicity gradients, and surface brightness profiles that are consistent with observations. In agreement with other simulations, the FOGGIE stellar halos receive 30%–40% of their mass from in situ stars. However, this population is more centrally concentrated in the FOGGIE simulations and therefore does not contribute excess light to the halo outskirts. The remaining stars are accreted from ∼10–50 other galaxies, with the majority of the accreted mass originating in two to four galaxies. While the inner halo ( r < 50 kpc) of each FOGGIE galaxy has a large number of contributors, the halo outskirts of three of the five galaxies are primarily made up of stars from only a few contributors. We predict that upcoming wide-field observatories, like the Nancy Grace Roman Space Telescope, will probe stellar halos around Milky Way–like galaxies out to ∼100 kpc in integrated light and will be able to distinguish the debris of dwarf galaxies with extended star formation histories from the underlying halo with resolved color–magnitude diagrams.

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium, with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly resolved L ^* galaxies from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM (Δ x ≲ 1 kpc h ^−1 and M ≃ 200 M _⊙ ). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of ≲5 kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond ∼0.25 R _200 , or ∼50 kpc at z = 0. Within ∼0.25 R _200 , turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the nonequilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies, even when averaging out small-scale variations.