2022ApJ...930...69S

We investigate whether the considerable diversity in the satellite populations of nearby Milky Way (MW)-mass galaxies is connected with the diversity in their host's merger histories. Analyzing eight nearby galaxies with extensive observations of their satellite populations and stellar halos, we characterize each galaxy's merger history using the metric of its most dominant merger, M_*,Dom, defined as the greater of either its total accreted stellar mass or most massive current satellite. We find an unexpectedly tight relationship between these galaxies' number of M_V < - 9 satellites within 150 kpc (N_Sat) and M_*,Dom. This relationship remains even after accounting for differences in galaxy mass. Using the star formation and orbital histories of satellites around the MW and M81, we demonstrate that both likely evolved along the M_*,Dom-N_Sat relation during their current dominant mergers with the Large Magellanic Cloud and M82, respectively. We investigate the presence of this relation in galaxy formation models, including using the Feedback In Realistic Environments (FIRE) simulations to directly compare to the observations. We find no relation between M_*,Dom and N_Sat in FIRE, and a universally large scatter in N_Sat with M_*,Dom across simulations-in direct contrast with the tightness of the empirical relation. This acute difference in the observed and predicted scaling relation between two fundamental galaxy properties signals that current simulations do not sufficiently reproduce diverse merger histories and their effects on satellite populations. Explaining the emergence of this relation is therefore essential for obtaining a complete understanding of galaxy formation.