On the Origin of a Possible Hard Very High Energy Spectrum from M87 Discovered by LHAASO

Recent LHAASO observations hint at potential spectral hardening around 20 TeV in M87’s very high energy emission, suggesting a possible new radiation component. In this work, we construct averaged multiwavelength spectral energy distributions by combining data from Chandra and Swift-UVOT/XRT covering the same period as the LHAASO detection to investigate the origin of this feature. We test several radiation mechanisms, including the pp interaction, proton synchrotron emission, photomeson process, and two-zone leptonic model. We find that only the pion decay gamma rays in pp interactions can interpret this feature in the framework of the one-zone model. With analytical analysis, we prove that proton synchrotron emission cannot generate a hard spectrum above 0.17 TeV. The photomeson model requires an emission zone compressed near the Schwarzschild radius of the central supermassive black hole, incompatible with broadband optical-GeV spectral constraints. In addition, the two-zone leptonic model also emerges as a viable alternative.