2024ApJ...970...70W

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.