Simulating the Diffuse Neutrino Emission from the Milky Way with GALPROP

We use the GALPROP cosmic ray (CR) propagation framework to model the diffuse neutrino and gamma-ray emissions from the Galaxy. A collection of realistic bounding models are developed and predictions of the resulting neutrino and gamma-ray signals are compared to the IceCube and LHAASO data up to PeV energies. We find that all the GALPROP models are consistent with the neutrino data within uncertainties. They are also consistent with expectations of neutrino emissions derived from LHAASO data when accounting for possible gamma-ray point source contamination. The new models present state-of-the-art predictions for the VHE neutrino emissions from the Galaxy that may be used for future neutrino searches.