Abstract:Context. Supernova remnants (SNRs) are likely to be significant sources of cosmic rays up to the knee of the local cosmic-ray (CR) spectrum. They produce gamma rays in the very-high-energy (VHE) ( E > 0.1 TeV) range mainly via two mechanisms: hadronic interactions of accelerated protons with the interstellar medium and leptonic interactions of accelerated electrons with soft photons. Observations with current instruments have lead to the detection of about a dozen SNRs emitting VHE gamma rays and future instruments should significantly increase this number. Yet, the details of particle acceleration at SNRs and of the mechanisms producing VHE gamma-rays at SNRs remain poorly understood. Aims. We aim to study the population of SNRs detected in the TeV range and its properties and confront it to simulated samples in order to address fundamental questions concerning particle acceleration at SNR shocks. Such questions concern the spectrum of accelerated particles, the efficiency of particle acceleration, and the gamma-ray emission being dominated by hadronic or leptonic interactions. Methods. By means of Monte Carlo methods, we simulated the population of SNRs in the gamma-ray domain and confronted our simulations to the catalogue of sources from the High Energy Stereoscopic System (H.E.S.S.) Galactic Plane Survey (HGPS). Results. We systematically explored the parameter space defined in our model, including for example, the slope of accelerated particles α , the electron-to-proton ratio K ep , and the efficiency of particle acceleration ξ . In particular, we found possible sets of parameters for which ≳90% of Monte Carlo realisations are found to be in agreement with the HGPS. These parameters are typically found at 4.2 ≳ α ≳ 4.1, 10 −5 ≲ K ep ≲ 10 −4.5 , and 0.03 ≲ ξ ≲ 0.1. We are able to strongly argue against some regions of the parameter space describing the population of Galactic SNRs in the TeV range, such as α ≲ 4.05, α ≳ 4.35, or K ep ≳ 10 −3 . Conclusions. Our model is so far able to explain the SNR population of the HGPS. Our approach, when confronted with the results of future systematic surveys, such as the Cherenkov Telescope Array Observatory, will help remove degeneracy from the solutions and to better understand particle acceleration at SNR shocks in the Galaxy.

Supernova Remnant G106.3+2.7: a Likely Proton PeVatron

Revealing a Peculiar Supernova Remnant G106.3+2.7 As a Petaelectronvolt Proton Accelerator with X-ray Observations

On the Hard Gamma-Ray Spectrum of the Potential PeVatron Supernova Remnant G106.3+2.7

On the hard γ-ray spectrum of the supernova remnant G106.3+2.7

Supernova remnants as cosmic ray accelerators. SNR IC 443

PeVatron candidate SNR G106.3+2.7 in a low-density cavity: a multiwavelength test

The Role of the Diffusive Protons in the Gamma-Ray Emission of Supernova Remnant Rx J1713.7-3946-A Two-Zone Model

Non-thermal emission from young supernova remnants in dense circumstellar environments

High-energy gamma-ray study of the dynamically young SNR G150.3+4.5

The population of Galactic supernova remnants in the TeV range

Search for synchrotron emission from secondary electrons of proton-proton interaction in Galactic PeVatron candidate HESS J1641$-$463

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Is Supernova Remnant Cassiopeia A a PeVatron?

The Maximum Energy of Shock-Accelerated Cosmic Rays

SNR G54.1+0.3, a PeVatron candidate unveiled by LHAASO

Search for Synchrotron Emission from Secondary Electrons of ProtonвҖ“Proton Interactions in Galactic PeVatron Candidate HESS J1641вҖ“463

Revisiting the Distance, Environment, and Supernova Properties of SNR G57.2+0.8 That Hosts SGR 1935+2154

Supernova Remnants in Gamma Rays

The Γ-Ray Emission Produced by Protons That Escape from Supernova Remnant G349.7+0.2