DELVE Milky Way Satellite Galaxy Census I: Satellite Population and Survey Selection Function in DES, DELVE, and Pan-STARRS

The properties of Milky Way satellite galaxies have important implications for galaxy formation, reionization, and the fundamental physics of dark matter. However, the population of Milky Way satellites includes the faintest known galaxies, and current observations are incomplete. To understand the impact of observational selection effects on the known satellite population, we perform rigorous, quantitative estimates of the Milky Way satellite galaxy detection efficiency in three wide-field survey datasets: the Dark Energy Survey Year 6, the DECam Local Volume Exploration Data Release 3, and the Pan-STARRS1 Data Release 1. Together, these surveys cover $\sim$13,600 deg$^2$ to $g \sim 24.0$ and $\sim$27,700 deg$^2$ to $g \sim 22.5$, spanning $\sim$91% of the high-Galactic-latitude sky ($|b| \geq 15^\circ$). We apply multiple detection algorithms over the combined footprint and recover 49 known satellites above a strict census detection threshold. To characterize the sensitivity of our census, we run our detection algorithms on a large set of simulated galaxies injected into the survey data, which allows us to develop models that predict the detectability of satellites as a function of their properties. We then fit an empirical model to our data and infer the luminosity function, radial distribution, and size-luminosity relation of Milky Way satellite galaxies. Our empirical model predicts a total of $265^{+79}_{-47}$ satellite galaxies with $-20 \leq M_V \leq 0$, half-light radii of $15 \leq r_{1/2} (\rm pc) \leq 3000$, and galactocentric distances of $10 \leq D_{\rm GC} (\rm kpc) \leq 300$. We also identify a mild anisotropy in the angular distribution of the observed galaxies, at a significance of $\sim$$2\sigma$, which can be attributed to the clustering of satellites associated with the LMC.