Galactic cosmic rays reach energies of at least a few Peta-electronvolts (1 PeV =1015 electron volts)1 . This implies our Galaxy contains PeV accelerators (PeVatrons), but all proposed models of Galactic cosmic-ray accelerators encounter non-trivial difficulties at exactly these energies 2 . Tens of Galactic accelerators capable of accelerating particle to tens of TeV (1 TeV =10 12 electron volts) energies were inferred from recent gamma-ray observations 3 . None of the currently known accelerators, however, not even the handful of shell-type supernova remnants commonly believed to supply most Galactic cosmic rays, have shown the characteristic tracers of PeV particles: power-law spectra of gamma rays extending without a cutoff or a spectral break to tens of TeV 4 . Here we report deep gamma-ray observations with arcminute angular resolution of the Galactic Centre regions, which show the expected tracer of the presence of PeV particles within the central 10 parsec of the Galaxy. We argue that the supermassive black hole Sagittarius A* is linked to this PeVatron. Sagittarius A* went through active phases in the past, as demonstrated by X-ray outbursts 5 and an outflow from the Galactic Center 6 . Although its current rate of particle acceleration is not sufficient to provide a substantial contribution to Galactic cosmic rays, Sagittarius A* could have plausibly been more active over the last 10 6−7 years, and therefore should be considered as a viable alternative to supernova remnants as a source of PeV Galactic cosmic rays.The large photon statistics accumulated over the last 10 years of observations with the High Energy Stereoscopic System (H.E.S.S.), together with improvements in the methods of data analysis, allow for a deep study of the properties of the diffuse very-high-energy (VHE; more than 100 GeV) emission of the central molecular zone. This region surrounding the Galactic Centre contains predominantly molecular gas and extends (in projection) out to r∼250 pc at positive galactic longitudes and r∼150 pc at negative longitudes. The map of the central molecular zone as seen in VHE γ-rays (Fig. 1) shows a strong (although not linear; see below) correlation between the brightness distribution of VHE γ-rays and the locations of massive gas-rich complexes. This points towards a hadronic origin of the diffuse emission 7 , where the γ-rays result from the interactions of relativistic protons with the ambient gas. The second important mechanism of production of VHE γ-rays 3 is the inverse Compton scattering of electrons. However, the severe radiative losses suffered by multi-TeV electrons in the Galactic Centre region prevent them from propagating over scales comparable to the size of the central molecular zone, thus disfavouring a leptonic origin of the γ-rays (see discussion in Methods and Extended Data Figures 1 and 2). The location and the particle injection rate history of the cosmic-ray accelerator(s), responsible for the relativistic protons, determine the spatial distribution of these cosmic rays which...

A deep survey of the Large Magellanic Cloud at ∼ 0.1−100 TeV photon energies with the Cherenkov Telescope Array is planned. We assess the detection prospects based on a model for the emission of the galaxy, comprising the four known TeV emitters, mock populations of sources, and interstellar emission on galactic scales. We also assess the detectability of 30 Doradus and SN 1987A, and the constraints that can be derived on the nature of dark matter. The survey will allow for fine spectral studies of N 157B, N 132D, LMC P3, and 30 Doradus C, and half a dozen other sources should be revealed, mainly pulsar-powered objects. The remnant from SN 1987A could be detected if it produces cosmic-ray nuclei with a flat power-law spectrum at high energies, or with a steeper index 2.3 − 2.4 pending a flux increase by a factor > 3 − 4 over ∼ 2015 − 2035. Large-scale interstellar emission remains mostly out of reach of the survey if its > 10 GeV spectrum has a soft photon index ∼ 2.7, but degree-scale 0.1 − 10 TeV pion-decay emission could be detected if the cosmic-ray spectrum hardens above >100 GeV. The 30 Doradus star-forming region is detectable if acceleration efficiency is on the order of 1 − 10% of the mechanical luminosity and diffusion is suppressed by two orders of magnitude within < 100 pc. Finally, the survey could probe the canonical velocity-averaged cross section for self-annihilation of weakly interacting massive particles for cuspy Navarro-Frenk-White profiles.

The nine-year H.E.S.S. Galactic Plane Survey (HGPS) has yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also present several upper limits for regions around pulsars without a detected TeV wind nebula. Our data exhibit a correlation of TeV surface brightness with pulsar spin-down powerĖ. This seems to be caused both by an increase of extension with decreasingĖ, and hence with time, compatible with a power law RPWN(Ė) ∼Ė −0.65±0.20 , and by a mild decrease of TeV gamma-ray luminosity with decreasingĖ, compatible with L1−10 TeV ∼Ė 0.59±0.21 . We also find that the offsets of pulsars with respect to the wind nebula centre with ages around 10 kyr are frequently larger than can be plausibly explained by pulsar proper motion and could be due to an asymmetric environment. In the present data, it seems that a large pulsar offset is correlated with a high apparent TeV efficiency L1−10 TeV /Ė. In addition to 14 HGPS sources considered firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that are likely TeV pulsar wind nebula candidates. Using a model that subsumes the present common understanding of the very high-energy radiative evolution of pulsar wind nebulae, we find that the trends and variations of the TeV observables and limits can be reproduced to a good level, drawing a consistent picture of present-day TeV data and theory.

Gamma-ray bursts (GRBs), which are bright flashes of gamma rays from extragalactic sources followed by fading afterglow emission, are associated with stellar core collapse events. We report the detection of very-high-energy (VHE) gamma rays from the afterglow of GRB 190829A, between 4 and 56 hours after the trigger, using the High Energy Stereoscopic System (H.E.S.S.). The low luminosity and redshift of GRB 190829A reduce both internal and external absorption, allowing determination of its intrinsic energy spectrum. Between energies of 0.18 and 3.3 tera–electron volts, this spectrum is described by a power law with photon index of 2.07 ± 0.09, similar to the x-ray spectrum. The x-ray and VHE gamma-ray light curves also show similar decay profiles. These similar characteristics in the x-ray and gamma-ray bands challenge GRB afterglow emission scenarios.

The diffuse very high−energy (VHE, > 100 GeV) γ-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow the dense matter distribution in the Central Molecular Zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic Centre (GC), where the flux rapidly decreases. This was initially interpreted as the result of a burst−like injection of energetic particles 10 4 years ago, but a recent more sensitive H.E.S.S. analysis revealed that the cosmic−ray (CR) density profile drops with the distance to the centre, making data compatible with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis to obtain for the first time a detailed characterisation of the correlation with matter and to search for additional features and individual γ-ray sources in the inner 200 pc. Taking advantage of 250 hours of H.E.S.S. data and improved analysis techniques we perform a detailed morphology study of the diffuse VHE emission observed from the GC ridge and reconstruct its total spectrum. To test the various contributions to the total γ-ray emission, we use an iterative 2D maximum likelihood approach that allows us to build a phenomenological model of the emission by summing a number of different spatial components. We show that the emission correlated with dense matter covers the full CMZ and that its flux is about half the total diffuse emission flux. We also detect some emission at higher latitude likely produced by hadronic collisions of CRs in less dense regions of the GC interstellar medium. We detect an additional emission component centred on the GC and extending over about 15 pc that is consistent with the existence of a strong CR density gradient and confirms the presence of a CR accelerator at the very centre of our Galaxy. We show that the spectrum of the full ridge diffuse emission is compatible with the one previously derived from the central regions, suggesting that a single population of particles fills the entire CMZ. Finally, we report the discovery of a VHE γ-ray source near the GC radio arc and argue that it is produced by the pulsar wind nebula candidate G0.13−0.11.