Galactic diffuse gamma rays meet the PeV frontier

Abstract: The Tibet AS$\gamma$ and LHAASO collaborations recently reported the observation of a $\gamma$-ray diffuse emission with energy up to the PeV from the Galactic plane. We discuss the relevance of non-uniform cosmic-ray transport scenarios and the implications of these results for cosmic-ray physics. We use the {\tt DRAGON} and {\tt HERMES} codes to build high-resolution maps and spectral distributions of that emission for several representative models under the condition that they reproduce a wide set of local cosmic-ray data up to 100 PeV. We show that the energy spectra measured by Tibet AS$\gamma$, LHAASO, ARGO-YBJ and Fermi-LAT in several regions of interest in the sky can all be consistently described in terms of the emission arising by the Galactic cosmic-ray "sea". We also show that all our models are compatible with IceTop $\gamma$-ray upper limits. Our results favor transport models characterized by spatial-dependent diffusion although some degeneracy remains between the choice of the transport scenario and that of the cosmic-ray spectral shape above 10 TeV. We discuss the role of forthcoming measurements in resolving that ambiguity.