Yanbo Wang, Zhenglong Wang, Rui Zhang, Yi Zhang

The Galactic diffuse gamma-ray emission is conventionally modeled as the product of cosmic-ray interactions with the interstellar medium. However, the cumulative contribution of stellar atmospheres acting as hadronic interaction targets remains an unexplored multi-messenger background. In this work, we present the first systematic evaluation of this stellar diffuse emission by coupling MESA stellar evolution profiles and magnetic-field-modulated cosmic-ray transport with a 3D Galactic population synthesis framework. We find that the cumulative stellar contribution to the Galactic diffuse gamma-ray flux is negligible at 1 TeV, and the associated diffuse neutrino flux ($\sim 10^{-16}\;\mathrm{TeV\;cm^{-2}\;s^{-1}\;sr^{-1}}$) remains orders of magnitude below current IceCube limits. Nevertheless, at ultra-high energies ($>10\;\mathrm{TeV}$), this emission establishes an irreducible local background that overtakes the strongly attenuated extragalactic isotropic gamma-ray background. Our results demonstrate that the Galactic stellar ensemble is a strictly sub-dominant background, indicating that stellar subtraction templates are not required for identifying Galactic PeVatrons or constraining dark matter annihilation.