We model the hadronic gamma-ray production in the Galaxy using an updated phenomenological description of the cosmic ray spectra and different Galactic gas templates. We find a persistent discrepancy in both normalization and spectral shape between the predicted gamma-ray flux and observations by LHAASO. We conclude that the local cosmic-ray proton flux measured by LHAASO overshoots the Galactic gamma-ray emission observed in the TeV-PeV sky, calling for a revision of current cosmic ray models.
Castro et al. (Tue,) studied this question.