Multiwavelength investigation of the candidate Galactic PeVatron MGRO J1908+06

The candidate PeVatron MGRO J1908+06, which shows a hard spectrum beyond 100 TeV, is one of the most peculiar -ray sources in the Galactic plane. Its complex morphology and some possible counterparts spatially related with the VHE emission region, preclude to distinguish between a hadronic or leptonic nature of the -ray emission. In this paper we illustrate a new multiwavelength analysis of MGRO J1908+06, with the aim to shed light on its nature and the origin of its ultra high-energy emission. We performed an analysis of the 12CO and 13CO molecular line emission demonstrating the presence of dense molecular clouds spatially correlated with the source region. We also analyzed 12-years of Fermi-LAT data between 10 GeV and 1 TeV finding a counterpart with a hard spectrum (Γ ∼ 1.6). Our reanalysis of XMM-Newton data allowed us to put a more stringent constraint on the X-ray flux from this source. We demonstrate that a single accelerator cannot explain the whole set of multiwavelength data, regardless of whether it accelerates protons or electrons, but a 2-zone model is needed to explain the emission from MGRO J1908+06. The VHE emission seems most likely the superposition of a TeV PWN powered by PSR J1907+0602, in the southern part, and of the interaction between the supernova remnant G40.5-0.5 and the molecular clouds towards the northern region