Stellar Population Spectra Incorporating Detailed Binary Evolution using POSYDON

Abstract: The accuracy of stellar population properties inferred through spectral energy distribution fitting hinges on the reliability of the underlying spectral models. Binary interactions are fundamental for massive star evolution, and ignoring their spectral contribution can lead to incorrect results. We use the POSYDON binary population synthesis code to generate spectral models of stellar populations that include binaries at solar metallicity. Our framework incorporates a collection of spectral libraries that is designed to address key outcomes of binary stellar evolution like Wolf-Rayet stars, stripped helium stars, and a treatment for stellar mergers. Our models confirm previous results showing that the inclusion of binary interactions has a significant effect on the UV and ionizing regime of the integrated spectrum. In particular we find that Wolf-Rayet stars and other massive stars dominate the production of ionizing radiation at earlier times, but after $\simeq$16 Myr stripped stars produced through mass transfer begin to dominate. Furthermore, we show that the production of ionizing He II photons is especially sensitive to the underlying population of stripped stars. While our results currently focus on high-mass stars ($\ge4~M_{\odot}$) at Solar metallicity, they provide the framework for binary spectral synthesis across a range of metallicities and masses and lay the foundation for calculations of the emergent emission-line spectra in the UV, optical, and IR regimes. We make the spectral models from this work publicly available for use in a format that can be integrated into fitting codes.