Cosmological Predictions for Minor Axis Stellar Density Profiles in the Inner Regions of Milky Way–mass Galaxies

ΛCDM cosmology predicts the hierarchical formation of galaxies, which build up mass by merger events and accreting smaller systems. The stellar halo of the Milky Way (MW) has proven to be useful a tool for tracing this accretion history. However, most of this work has focused on the outer halo where dynamical times are large and the dynamical properties of accreted systems are preserved. In this work, we investigate the inner galaxy regime, where dynamical times are relatively small and systems are generally completely phase mixed. Using the FIRE-2 and Auriga cosmological zoom-in simulation suites of MW-mass galaxies, we find the stellar density profiles along the minor axis (perpendicular to the galactic disk) within the Navarro–Frenk–White scale radii (R ≈ 15 kpc) are best described as an exponential disk with scale height < 0.3 kpc and a power-law component with slope α ≈ −4. The stellar density amplitude and slope for the power-law component are not significantly correlated with metrics of the galaxy's accretion history. Instead, we find the stellar profiles strongly correlate with the dark matter profile. Across simulation suites, the galaxies studied in this work have a stellar-to-dark-matter mass ratio that decreases as 1/r² along the minor axis.