The nearby radio galaxy Centaurus A belongs to a class of Active Galaxies that are very luminous at radio wavelengths. The majority of these galaxies show collimated relativistic outflows known as jets, that extend over hundreds of thousands of parsecs for the most powerful sources. Accretion of matter onto the central super-massive black hole is believed to fuel these jets and power their emission, with the radio emission being related to the synchrotron radiation of relativistic electrons in magnetic fields. The origin of the extended X-ray emission seen in the kiloparsec-scale jets from these sources is still a matter of debate, although Cen A's X-ray emission has been suggested to originate in electron synchrotron processes. The other possible explanation is Inverse Compton (IC) scattering with CMB soft photons. Synchrotron radiation needs ultra-relativistic electrons ($\sim50$ TeV), and given their short cooling times, requires some continuous re-acceleration mechanism to be active. IC scattering, on the other hand, does not require very energetic electrons, but requires jets that stay highly relativistic on large scales ($\geq$1 Mpc) and that remain well-aligned with the line of sight. Some recent evidence disfavours inverse Compton-CMB models, although other evidence seems to be compatible with them. In principle, the detection of extended gamma-ray emission, directly probing the presence of ultra-relativistic electrons, could distinguish between these options, but instruments have hitherto been unable to resolve the relevant structures. At GeV energies there is also an unusual spectral hardening in Cen A, whose explanation is unclear. Here we report observations of Cen A at TeV energies that resolve its large-scale jet. We interpret the data as evidence for the acceleration of ultra-relativistic electrons in the jet, and favour the synchrotron explanation for the X-rays.

The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution, out-of-equilibrium ionization state and gas chemistry during the collapse of a prestellar core. To compute the grain-grain collisional rate, we consider models for both random and systematic, size-dependent, velocities. We include grain growth through grain-grain coagulation and ice accretion, but ignore grain fragmentation. Starting with a MRN (Mathis et al. 1977) size distribution, we find that coagulation in grain-grain collisions generated by hydrodynamical turbulence is not efficient at removing the smallest grains, and as a consequence does not affect much the evolution of the Hall and ambipolar diffusion MHD resistivities which still severly drop during the collapse like in models without coagulation. The inclusion of systematic velocities, possibly induced by the presence of ambipolar diffusion, increases the coagulation rate between small and large grains, removing small grains earlier in the collapse and therefore limiting the drop in the Hall and ambipolar diffusion resistivities. At intermediate densities ($n_{\rm H} \sim 10^8\,$cm$^{-3}$), the Hall and ambipolar diffusion resistivities are found to be higher by 1 to 2 orders of magnitude in models with coagulation than in models where coagulation is ignored, and also higher than in a toy model without coagulation where all grains smaller than $0.1\,{\mu}$m would have been removed in the parent cloud before the collapse. When grain drift velocities induced by ambipolar diffusion are included, dust coagulation ... (abridged)

The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $\gamma$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which lasted more than 100 days, were somewhat unexpected. Current theoretical models attempt to explain this temporal behavior as either originating from a relativistic off-axis jet or a kilonova-like outflow. In either scenario, there is some ambiguity regarding how much energy is transported in the non-thermal electrons versus the magnetic field of the emission region. Combining the VLA (radio) and Chandra (X-ray) measurements with observations in the GeV-TeV domain can help break this ambiguity, almost independently of the assumed origin of the emission. Here we report for the first time on deep H.E.S.S. observations of GW170817 / GRB 170817A between 124 and 272 days after the BNS merger with the full H.E.S.S. array of telescopes, as well as on an updated analysis of the prompt (<5 days) observations with the upgraded H.E.S.S. phase-I telescopes. We discuss implications of the H.E.S.S. measurement for the magnetic field in the context of different source scenarios.

We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $\gamma$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE $\gamma$-ray regime, we deduce a 95% confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98, and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results the redshift of KUV 00311-1938 is constrained to $0.51 \leq z < 0.98$ and for PKS 1440-389 to $0.14 \lessapprox z < 0.53$.

Some core-collapse supernovae are likely to be efficient cosmic-ray accelerators up to the PeV range, and therefore, to potentially play an important role in the overall Galactic cosmic-ray population. The TeV gamma-ray domain can be used to study particle acceleration in the multi-TeV and PeV range. This motivates the study of the detectability of such supernovae by current and future gamma-ray facilities. The gamma-ray emission of core-collapse supernovae strongly depends on the level of the two-photon annihilation process: high-energy gamma-ray photons emitted at the expanding shock wave following the supernova explosion can interact with soft photons from the supernova photosphere through the pair production channel, thereby strongly suppressing the flux of gamma rays leaving the system. In the case of SN 1993J, whose photospheric and shock-related parameters are well measured, we calculate the temporal evolution of the expected gamma-ray attenuation by accounting for the temporal and geometrical effects. We find the attenuation to be of about $10$ orders of magnitude in the first few days after the SN explosion. The probability of detection of a supernova similar to SN 1993J with the Cherenkov Telescope Array is highest if observations are performed either earlier than 1 day, or later than 10 days after the explosion, when the gamma-ray attenuation decreases to about $2$ orders of magnitude.

Non-thermal particles and high-energy radiation can play a role in the dynamical processes in star-forming regions and provide an important piece of the multiwavelength observational picture of their structure and components. Powerful stellar winds and supernovae in compact clusters of massive stars and OB associations are known to be favourable sites of high-energy particle acceleration and sources of non-thermal radiation and neutrinos. Namely, young massive stellar clusters are likely sources of the PeV (petaelectronvolt) regime cosmic rays (CRs). They can also be responsible for the cosmic ray composition, e.g., 22Ne/20Ne anomalous isotopic ratio in CRs. Efficient particle acceleration can be accompanied by super-adiabatic amplification of the fluctuating magnetic fields in the systems converting a part of kinetic power of the winds and supernovae into the magnetic energy through the CR-driven instabilities. The escape and CR propagation in the vicinity of the sources are affected by the non-linear CR feedback. These effects are expected to be important in starburst galaxies, which produce high-energy neutrinos and gamma-rays. We give a brief review of the theoretical models and observational data on high-energy particle acceleration and their radiation in star-forming regions with young stellar population.

This review aims at providing an up-to-date status and a general introduction to the subject of the numerical study of energetic particle acceleration and transport in turbulent astrophysical flows. The subject is also complemented by a short overview of recent progresses obtained in the domain of laser plasma experiments. We review the main physical processes at the heart of the production of a non-thermal distribution in both Newtonian and relativistic astrophysical flows, namely the first and second order Fermi acceleration processes. We also discuss shock drift and surfing acceleration, two processes important in the context of particle injection in shock acceleration. We analyze with some details the particle-in-cell (PIC) approach used to describe particle kinetics. We review the main results obtained with PIC simulations in the recent years concerning particle acceleration at shocks and in reconnection events. The review discusses the solution of Fokker-Planck problems with application to the study of particle acceleration at shocks but also in hot coronal plasmas surrounding compact objects. We continue by considering large scale physics. We describe recent developments in magnetohydrodynamic (MHD) simulations. We give a special emphasize on the way energetic particle dynamics can be coupled to MHD solutions either using a multi-fluid calculation or directly coupling kinetic and fluid calculations. This aspect is mandatory to investigate the acceleration of particles in the deep relativistic regimes to explain the highest Cosmic Ray energies.

In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $\gamma$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) {\gamma}-ray emitters. {\eta} Car is the most prominent member of this object class and is confirmed to emit phase-locked HE {\gamma} rays from hundreds of MeV to ~100 GeV energies. This work aims to search for and characterise the very-high-energy (VHE; E >100 GeV) {\gamma}-ray emission from {\eta} Car around the last periastron passage in 2014 with the ground-based High Energy Stereoscopic System (H.E.S.S.). Methods. The region around {\eta} Car was observed with H.E.S.S. between orbital phase p = 0.78 - 1.10, with a closer sampling at p {\approx} 0.95 and p {\approx} 1.10 (assuming a period of 2023 days). Optimised hardware settings as well as adjustments to the data reduction, reconstruction, and signal selection were needed to suppress and take into account the strong, extended, and inhomogeneous night sky background (NSB) in the {\eta} Car field of view. Tailored run-wise Monte-Carlo simulations (RWS) were required to accurately treat the additional noise from NSB photons in the instrument response functions. Results. H.E.S.S. detected VHE {\gamma}-ray emission from the direction of {\eta} Car shortly before and after the minimum in the X-ray light-curve close to periastron. Using the point spread function provided by RWS, the reconstructed signal is point-like and the spectrum is best described by a power law. The overall flux and spectral index in VHE {\gamma} rays agree within statistical and systematic errors before and after periastron. The {\gamma}-ray spectrum extends up to at least ~400 GeV. This implies a maximum magnetic field in a leptonic scenario in the emission region of 0.5 Gauss. No indication for phase-locked flux variations is detected in the H.E.S.S. data.

PSR B1259-63/LS 2883 is a gamma-ray binary system consisting of a pulsar in an eccentric orbit around a bright Oe stellar-type companion star that features a dense circumstellar disc. The high- and very-high-energy (HE, VHE) gamma-ray emission from PSR B1259-63/LS 2883 around the times of its periastron passage are characterised, in particular, at the time of the HE gamma-ray flares reported to have occurred in 2011, 2014, and 2017. Spectra and light curves were derived from observations conducted with the H.E.S.S.-II array in 2014 and 2017. A local double-peak profile with asymmetric peaks in the VHE light curve is measured, with a flux minimum at the time of periastron $t_p$ and two peaks coinciding with the times at which the neutron star crosses the companion's circumstellar disc ($\sim t_p \pm 16$ d). A high VHE gamma-ray flux is also observed at the times of the HE gamma-ray flares ($\sim t_p + 30$ d) and at phases before the first disc crossing ($\sim t_p - 35$ d). PSR B1259-63/LS 2883 displays periodic flux variability at VHE gamma-rays without clear signatures of super-orbital modulation in the time span covered by H.E.S.S. observations. In contrast, the photon index of the measured power-law spectra remains unchanged within uncertainties for about 200 d around periastron. Lower limits on exponential cut-off energies up to $\sim 40$ TeV are placed. At HE gamma-rays, PSR B1259-63/LS 2883 has now been detected also before and after periastron, close to the disc crossing times. Repetitive flares with distinct variability patterns are detected in this energy range. Such outbursts are not observed at VHEs, although a relatively high emission level is measured. The spectra obtained in both energy regimes displays a similar slope, although a common physical origin either in terms of a related particle population, emission mechanism, or emitter location is ruled out.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Gamma-ray bursts (GRBs) are brief flashes of gamma rays, considered to be the most energetic explosive phenomena in the Universe. The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow -- produced by the interaction between the ejected matter and the circumburst medium -- slows down, and a gradual decrease in brightness is observed. GRBs typically emit most of their energy via gamma-rays with energies in the kiloelectronvolt-to-megaelectronvolt range, but a few photons with energies of tens of gigaelectronvolts have been detected by space-based instruments. However, the origins of such high-energy (above one gigaelectronvolt) photons and the presence of very-high-energy (more than 100 gigaelectronvolts) emission have remained elussive. Here we report observations of very-high-energy emission in the bright GRB 180720B deep in the GRB afterglow -ten hours after the end of the prompt emission phase, when the X-ray flux had already decayed by four orders of magnitude. Two possible explanations exist for the observed radiation: inverse Compton emission and synchrotron emission of ultrarelativistic electrons. Our observations show that the energy fluxes in the X-ray and gamma-ray range and their photon indices remain comparable to each other throughout the afterglow. This discovery places distinct constraints on the GRB environment for both emission mechanisms, with the inverse Compton explanation alleviating the particle energy requirements for the emission observed at late times. The late timing of this detection has consequences for the future observations of GRBs at the highest energies.

Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes triggered ToO observations on February 18, 2015, with the goal of studying the gamma-ray emission in the VHE band. H.E.S.S. ToO observations were carried out during the nights of February 18, 19, 21, and 24, 2015. Together with Fermi-LAT, the multi-wavelength coverage of the flare includes Swift observations in soft-X-rays and optical/UV, and optical monitoring (photometry and spectro-polarimetry) by the Steward Observatory, the ATOM, the KAIT and the ASAS-SN telescope. VHE emission from PKS 0736+017 was detected with H.E.S.S. during the night of February 19, 2015, only. Fermi data indicate the presence of a gamma-ray flare, peaking at the time of the H.E.S.S. detection, with a flux doubling time-scale of around six hours. The gamma-ray flare was accompanied by at least a 1 mag brightening of the non-thermal optical continuum. No simultaneous observations at longer wavelengths are available for the night of the H.E.S.S. detection. The gamma-ray observations with H.E.S.S. and Fermi are used to put constraints on the location of the gamma-ray emitting region during the flare: it is constrained to be just outside the radius of the broad-line-region with a bulk Lorentz factor $\simeq 20$, or at the level of the radius of the dusty torus with Gamma > 60. PKS 0736+017 is the seventh FSRQ known to emit VHE photons and, at z=0.189, is the nearest so far. The location of the gamma-ray emitting region during the flare can be tightly constrained thanks to opacity, variability, and collimation arguments.

We aim to elucidate cosmic ray (CR) propagation in the weakly ionized environments of supernova remnants (SNRs) basing our analysis on the cosmic ray cloud (CRC) model developed by \citet{2013ApJ...768...73M} and \citet{2016MNRAS.461.3552N}. We solved two transport equations simultaneously: one for the CR pressure and one for the Alfv\'en wave energy density where CRs are initially confined in the SNR shock. Cosmic rays trigger a streaming instability and produce slab-type resonant Alfv\'en modes. The self-generated turbulence is damped by ion-neutral collisions and by noncorrelated interaction with Alfv\'en modes generated at large scales. We show that CRs leaking in cold dense phases such as those found in cold neutral medium (CNM) and diffuse molecular medium (DiM) can still be confined over distances of a few tens of parsecs from the CRC center for a few thousand years. At 10 TeV, CR diffusion can be suppressed by two or three orders of magnitude. This effect results from a reduced ion-neutral collision damping in the decoupled regime. We calculate the grammage of CRs in these environments. We find that in both single and multi-phase setups at 10 GeV, CNM and DiM media can produce grammage in the range 10-20 $\rm{g/cm^2}$ in the CNM and DiM phases. At 10 TeV, because of nonlinear propagation the grammage increases to values in the range 0.5-20 $\rm{g/cm^2}$ in these two phases. We also present preliminary calculations in inhomogeneous interstellar medium combining two or three different phases where we obtain the same trends.

Cosmic rays (CRs) are thought to play a dynamically important role in several key aspects of galaxy evolution, including the structure of the interstellar medium, the formation of galactic winds, and the non-thermal pressure support of halos. We introduce a numerical model solving for the CR streaming instability and acceleration of CRs at shocks with a fluid approach in the adaptive mesh refinement code ramses. CR streaming is solved with a diffusion approach and its anisotropic nature is naturally captured. We introduce a shock finder for the ramses code that automatically detects shock discontinuities in the flow. Shocks are the loci for CR injection, and their efficiency of CR acceleration is made dependent on the upstream magnetic obliquity according to the diffuse shock acceleration mechanism. We show that the shock finder accurately captures shock locations and estimates the shock Mach number for several problems. The obliquity-dependent injection of CRs in the Sedov solution leads to situations where the supernova bubble exhibits large polar caps (homogeneous background magnetic field), or a patchy structure of the CR distribution (inhomogeneous background magnetic field). Finally, we combine both accelerated CRs with streaming in a simple turbulent interstellar medium box, and show that the presence of CRs significantly modifies the structure of the gas.

Core-collapse supernovae produce fast shocks which expand into the dense circumstellar medium (CSM) of the stellar progenitor. Cosmic rays (CRs) accelerated at these shocks can induce the growth of electromagnetic fluctuations in the pre-shock medium. Using a self-similar description for the shock evolution, we calculate the growth time-scales of CR driven instabilities for SNe in general, and SN 1993J in particular. We find that extended SN shocks can trigger fast intra-day instabilities, strong magnetic field amplification, and CR acceleration. In particular, the non-resonant streaming instability can contribute to about 50 per cent of the magnetic field intensity deduced from radio data. This results in the acceleration of CR particles to energies of 1-10 PeV within a few days after the shock breakout.

The Crab nebula is one of the most studied cosmic particle accelerators, shining brightly across the entire electromagnetic spectrum up to very high-energy gamma rays. It is known from radio to gamma-ray observations that the nebula is powered by a pulsar, which converts most of its rotational energy losses into a highly relativistic outflow. This outflow powers a pulsar wind nebula (PWN), a region of up to 10~light-years across, filled with relativistic electrons and positrons. These particles emit synchrotron photons in the ambient magnetic field and produce very high-energy gamma rays by Compton up-scattering of ambient low-energy photons. While the synchrotron morphology of the nebula is well established, it was up to now not known in which region the very high-energy gamma rays are emitted. Here we report that the Crab nebula has an angular extension at gamma-ray energies of 52 arcseconds (assuming a Gaussian source width), significantly larger than at X-ray energies. This result closes a gap in the multi-wavelength coverage of the nebula, revealing the emission region of the highest energy gamma rays. These gamma rays are a new probe of a previously inaccessible electron and positron energy range. We find that simulations of the electromagnetic emission reproduce our new measurement, providing a non-trivial test of our understanding of particle acceleration in the Crab nebula.

The giant molecular cloud Sagittarius B2 (hereafter SgrB2) is the most massive region with ongoing high-mass star formation in the Galaxy. In the southern region of the 40-pc large envelope of SgrB2, we encounter the SgrB2(DS) region which hosts more than 60 high-mass protostellar cores distributed in an arc shape around an extended HII region. We use the Very Large Array in its CnB and D configurations, and in the frequency bands C (4--8 GHz) and X (8--12 GHz) to observe the whole SgrB2 complex. Continuum and radio recombination line maps are obtained. We detect radio continuum emission in SgrB2(DS) in a bubble-shaped structure. From 4 to 12 GHz, we derive a spectral index between -1.2 and -0.4, indicating the presence of non-thermal emission. We decompose the contribution from thermal and non-thermal emission, and find that the thermal component is clumpy and more concentrated, while the non-thermal component is more extended and diffuse. The radio recombination lines in the region are found to be not in local thermodynamic equilibrium (LTE) but stimulated by the non-thermal emission. The thermal free-free emission is likely tracing an HII region ionized by an O7 star, while the non-thermal emission can be generated by relativistic electrons created through first-order Fermi acceleration. We have developed a simple model of the SgrB2(DS) region and found that first-order Fermi acceleration can reproduce the observed flux density and spectral index.

We use particle-in-magnetohydrodynamics-cells to model particle acceleration and magnetic field amplification in a high Mach, parallel shock in three dimensions and compare the result to 2-D models. This allows us to determine whether 2-D simulations can be relied upon to yield accurate results in terms of particle acceleration, magnetic field amplification and the growth rate of instabilities. Our simulations show that the behaviour of the gas and the evolution of the instabilities are qualitatively similar for both the 2-D and 3-D models, with only minor quantitative differences that relate primarily to the growth speed of the instabilities. The main difference between 2-D and 3-D models can be found in the spectral energy distributions (SEDs) of the non-thermal particles. The 2-D simulations prove to be more efficient, accelerating a larger fraction of the particles and achieving higher velocities. We conclude that, while 2-D models are sufficient to investigate the instabilities in the gas, their results have to be treated with some caution when predicting the expected SED of a given shock.

Radio observations at metre-centimetre wavelengths shed light on the nature of the emission of HII regions. Usually this category of objects is dominated by thermal radiation produced by ionised hydrogen, namely protons and electrons. However, a number of observational studies have revealed the existence of HII regions with a mixture of thermal and non-thermal radiation. The latter represents a clue as to the presence of relativistic electrons. However, neither the interstellar cosmic-ray electron flux nor the flux of secondary electrons, produced by primary cosmic rays through ionisation processes, is high enough to explain the observed flux densities. We investigate the possibility of accelerating local thermal electrons up to relativistic energies in HII region shocks. We assumed that relativistic electrons can be accelerated through the first-order Fermi acceleration mechanism and we estimated the emerging electron fluxes, the corresponding flux densities, and the spectral indexes. We find flux densities of the same order of magnitude of those observed. In particular, we applied our model to the "deep south" (DS) region of Sagittarius B2 and we succeeded in reproducing the observed flux densities with an accuracy of less than 20% as well as the spectral indexes. The model also gives constraints on magnetic field strength ($0.3-4$ mG), density ($1-9\times10^4$ cm$^{-3}$), and flow velocity in the shock reference frame ($33-50$ km s$^{-1}$) expected in DS. We suggest a mechanism able to accelerate thermal electrons inside HII regions through the first-order Fermi acceleration. The existence of a local source of relativistic electrons can explain the origin of both the observed non-thermal emission and the corresponding spectral indexes.

Core-collapse supernovae exploding in dense winds are favorable sites for cosmic-ray (CR) acceleration to very high energies. We present our CR-radiation-hydrodynamics simulations of the explosion of a red supergiant. We study the evolution of the shock wave during the first day following core collapse, and estimate the time at which CR acceleration can start. We then calculate the maximum CR energy at the forward shock as a function of time, and show that it may already exceed 100 TeV only a few hours after shock breakout from the surface of the star.

The flat spectrum radio quasar 3C 279 is known to exhibit pronounced variability in the high-energy ($100\,$MeV$<E<100\,$GeV) $\gamma$-ray band, which is continuously monitored with Fermi-LAT. During two periods of high activity in April 2014 and June 2015 Target-of-Opportunity observations were undertaken with H.E.S.S. in the very-high-energy (VHE, $E>100\,$GeV) $\gamma$-ray domain. While the observation in 2014 provides an upper limit, the observation in 2015 results in a signal with $8.7\,\sigma$ significance above an energy threshold of $66\,$GeV. No VHE variability has been detected during the 2015 observations. The VHE photon spectrum is soft and described by a power-law index of $4.2\pm 0.3$. The H.E.S.S. data along with a detailed and contemporaneous multiwavelength data set provide constraints on the physical parameters of the emission region. The minimum distance of the emission region from the central black hole is estimated using two plausible geometries of the broad-line region and three potential intrinsic spectra. The emission region is confidently placed at $r\gtrsim 1.7\times10^{17}\,$cm from the black hole, i.e., beyond the assumed distance of the broad-line region. Time-dependent leptonic and lepto-hadronic one-zone models are used to describe the evolution of the 2015 flare. Neither model can fully reproduce the observations, despite testing various parameter sets. Furthermore, the H.E.S.S. data are used to derive constraints on Lorentz invariance violation given the large redshift of 3C 279.

This article aims at establishing new benchmark scenarios for Galactic cosmic-ray propagation in the GV-TV rigidity range, based on fits to the AMS-02 B/C data with the USINE v3.5 propagation code. We employ a new fitting procedure, cautiously taking into account data systematic error correlations in different rigidity bins and considering Solar modulation potential and leading nuclear cross-section as nuisance parameters. We delineate specific low, intermediate, and high-rigidity ranges that can be related to both features in the data and peculiar microphysics mechanisms resulting in spectral breaks. We single out a scenario which yields excellent fits to the data and includes all the presumably relevant complexity, the BIG model. This model has two limiting regimes: (i) the SLIM model, a minimal diffusion-only setup, and (ii) the QUAINT model, a convection-reacceleration model where transport is tuned by non-relativistic effects. All models lead to robust predictions in the high-energy regime ($\gtrsim10$GV), i.e. independent of the propagation scenario: at $1\sigma$, the diffusion slope $\delta$ is $[0.43-0.53]$, whereas $K_{10}$, the diffusion coefficient at 10GV, is $[0.26-0.36]$kpc$^2$Myr$^{-1}$; we confirm the robustness of the high-energy break, with a typical value $\Delta_h\sim 0.2$. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak ($\sim 4$GV) which might be related to some effective damping scale in the magnetic turbulence.

Young core-collapse supernovae with dense-wind progenitors may be able to accelerate cosmic-ray hadrons beyond the knee of the cosmic-ray spectrum, and this may result in measurable gamma-ray emission. We searched for gamma-ray emission from ten supernovae observed with the High Energy Stereoscopic System (H.E.S.S.) within a year of the supernova event. Nine supernovae were observed serendipitously in the H.E.S.S. data collected between December 2003 and December 2014, with exposure times ranging from 1.4 hours to 53 hours. In addition we observed SN 2016adj as a target of opportunity in February 2016 for 13 hours. No significant gamma-ray emission has been detected for any of the objects, and upper limits on the $>1$ TeV gamma-ray flux of the order of $\sim$10$^{-13}$ cm$^{-2}$s$^{-1}$ are established, corresponding to upper limits on the luminosities in the range $\sim$2 $\times$ 10$^{39}$ erg s$^{-1}$ to $\sim$1 $\times$ 10$^{42}$ erg s$^{-1}$. These values are used to place model-dependent constraints on the mass-loss rates of the progenitor stars, implying upper limits between $\sim$2 $\times 10^{-5}$ and $\sim$2 $\times 10^{-3}$M$_{\odot}$yr$^{-1}$ under reasonable assumptions on the particle acceleration parameters.

PKS 1830-211 is a known macrolensed quasar located at a redshift of z=2.5. Its high-energy gamma-ray emission has been detected with the Fermi-LAT instrument and evidence for lensing was obtained by several authors from its high-energy data. Observations of PKS 1830-211 were taken with the H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes in August 2014, following a flare alert by the Fermi- LAT collaboration. The H.E.S.S observations were aimed at detecting a gamma-ray flare delayed by 20-27 days from the alert flare, as expected from observations at other wavelengths. More than twelve hours of good quality data were taken with an analysis threshold of $\sim67$ GeV. The significance of a potential signal is computed as a function of the date as well as the average significance over the whole period. Data are compared to simultaneous observations by Fermi-LAT. No photon excess or significant signal is detected. An upper limit on PKS 1830-211 flux above 67 GeV is computed and compared to the extrapolation of the Fermi-LAT flare spectrum.

MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.

Latest precise cosmic-ray (CR) measurements and present gamma-ray observations have started challenging our understanding of CR transport and interaction in the Galaxy. Moreover, because the density of CRs is similar to the density of the magnetic field, gas, and starlight in the interstellar medium (ISM), CRs are expected to affect the ISM dynamics, including the physical and chemical processes that determine transport and star formation. In this context, observations of gamma-ray emission at MeV energies produced by the low-energy CRs are very important and urgent. A telescope covering the energy range between ~0.1 MeV and a few GeV with a sensitivity more than an order of magnitude better than previous instruments would allow for the first time to study in detail the low-energy CRs, providing information on their sources, their spectra throughout the Galaxy, their abundances, transport properties, and their role on the evolution of the Galaxy and star formation. Here we discuss the scientific prospects for studies of CRs, ISM (gas, interstellar photons, and magnetic fields) and associated gamma-ray emissions with such an instrument.

The 511 keV line from positron annihilation in the Galaxy was the first $\gamma$-ray line detected to originate from outside our solar system. Going into the fifth decade since the discovery, the source of positrons is still unconfirmed and remains one of the enduring mysteries in $\gamma$-ray astronomy. With a large flux of $\sim$10$^{-3}$ $\gamma$/cm$^{2}$/s, after 15 years in operation INTEGRAL/SPI has detected the 511 keV line at $>50\sigma$ and has performed high-resolution spectral studies which conclude that Galactic positrons predominantly annihilate at low energies in warm phases of the interstellar medium. The results from imaging are less certain, but show a spatial distribution with a strong concentration in the center of the Galaxy. The observed emission from the Galactic disk has low surface brightness and the scale height is poorly constrained, therefore, the shear number of annihilating positrons in our Galaxy is still not well know. Positrons produced in $\beta^+$-decay of nucleosynthesis products, such as $^{26}$Al, can account for some of the annihilation emission in the disk, but the observed spatial distribution, in particular the excess in the Galactic bulge, remains difficult to explain. Additionally, one of the largest uncertainties in these studies is the unknown distance that positrons propagate before annihilation. In this paper, we will summarize the current knowledge base of Galactic positrons, and discuss how next-generation instruments could finally provide the answers.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is galactic cosmic rays, supernova remnants, protostellar jets and superbubbles, and colliding wind binaries.

We study the problem of the escape and transport of Cosmic-Rays (CR) from a source embedded in a fully ionised, hot phase of the interstellar medium (HIM). In particular, we model the CR escape and their propagation in the source vicinity taking into account excitation of Alfv\'enic turbulence by CR streaming and mechanisms damping the self-excited turbulence itself. Our estimates of escape radii and times result in large values (100 pc, $2\times10^5$ yr) for particle energies $\lesssim20$ GeV and smaller values for particles with increasing energies (35 pc and 14 kyr at 1 TeV). These escape times and radii, when used as initial conditions for the CR propagation outside the source, result in relevant suppression of the diffusion coefficient (by a factor 5-10) on time-scales comparable with their (energy dependent) escape time-scale. The damping mechanisms are fast enough that even on shorter time scales the Alfv\'enic turbulence is efficiently damped, and the ratio between random and ordered component of the magnetic field is $\delta B/B_0\ll 1$, justifying the use of quasi-linear theory. In spite of the suppressed diffusion coefficient, and then the increased residence time in the vicinity (<200 pc) of their source, the grammage accumulated by CRs after their escape is found to be negligible (at all energies) as compared to the one accumulated while diffusing in the whole Galaxy, due to the low density of the HIM.

The blazar Mrk 501 (z=0.034) was observed at very-high-energy (VHE, $E\gtrsim 100$~GeV) gamma-ray wavelengths during a bright flare on the night of 2014 June 23-24 (MJD 56832) with the H.E.S.S. phase-II array of Cherenkov telescopes. Data taken that night by H.E.S.S. at large zenith angle reveal an exceptional number of gamma-ray photons at multi-TeV energies, with rapid flux variability and an energy coverage extending significantly up to 20 TeV. This data set is used to constrain Lorentz invariance violation (LIV) using two independent channels: a temporal approach considers the possibility of an energy dependence in the arrival time of gamma rays, whereas a spectral approach considers the possibility of modifications to the interaction of VHE gamma rays with extragalactic background light (EBL) photons. The non-detection of energy-dependent time delays and the non-observation of deviations between the measured spectrum and that of a supposed power-law intrinsic spectrum with standard EBL attenuation are used independently to derive strong constraints on the energy scale of LIV ($E_{\rm{QG}}$) in the subluminal scenario for linear and quadratic perturbations in the dispersion relation of photons. For the case of linear perturbations, the 95% confidence level limits obtained are $E_{\rm{QG},1} > 3.6 \times 10^{17} \ \rm{GeV} $ using the temporal approach and $E_{\rm{QG},1} > 2.6 \times 10^{19} \ \rm{GeV}$ using the spectral approach. For the case of quadratic perturbations, the limits obtained are $E_{\rm{QG},2} > 8.5 \times 10^{10} \ \rm{GeV} $ using the temporal approach and $E_{\rm{QG},2} > 7.8 \times 10^{11} \rm{ GeV}$ using the spectral approach.

Cosmic rays (CR) propagate through the galactic scales down to the smaller scales at which stars form. CRs are close to energy equipartition with the other components of the interstellar medium and can provide a support against gravity if pressure gradients develop. We study CRs propagation within a turbulent and magnetised bi-stable interstellar gas and identify the conditions necessary for CR trapping. We present a numerical study using 3D simulations of the evolution of a mixture of interstellar gas and CRs, in which turbulence is driven by stochastic forcing within a box of 40 pc. We explore a large parameter space (CR diffusion coefficient, magnetisation, driving scale and amplitude of the turbulence forcing, initial CR energy). We identify a clear transition in the interstellar dynamics for CR diffusion coefficient below a critical value, which depends on the characteristic length scale L as $D_{\rm crit}\simeq 3.1~10^{23}~{\rm~cm^2/s}~(L/{\rm 1~pc})^{q+1}$, where the exponent q relates the turbulent velocity dispersion to the length scale as $v\simeq L^q$. In our simulations this transition occurs around $D_{\rm crit}\simeq 10^{24-25}$ cm^2/s. The transition is recovered in all cases of our parameter study and is in very good agreement with our simple analytical estimate. In the trapped CR regime, the induced CR pressure gradients can modify the gas flow and provide a support against the thermal instability development. We discuss possible mechanisms that can reduce significantly the CR diffusion coefficients within the interstellar medium. CR pressure gradients can develop and modify the evolution of thermally bi-stable gas for diffusion coefficients $D\leq 10^{25}$ cm^2/s or in regions where the CR pressure exceeds the thermal one by factor>10. This study provides the basis of further works including more realistic CR diffusion coefficients, as well as local CR sources.

Aims: We present a detailed view of the pulsar wind nebula (PWN) HESS J1825-137. We aim to constrain the mechanisms dominating the particle transport within the nebula, accounting for its anomalously large size and spectral characteristics. Methods: The nebula is studied using a deep exposure from over 12 years of H.E.S.S. I operation, together with data from H.E.S.S. II improving the low energy sensitivity. Enhanced energy-dependent morphological and spatially-resolved spectral analyses probe the Very High Energy (VHE, E > 0.1 TeV) gamma-ray properties of the nebula. Results: The nebula emission is revealed to extend out to 1.5 degrees from the pulsar, ~1.5 times further than previously seen, making HESS J1825--137, with an intrinsic diameter of ~100 pc, potentially the largest gamma-ray PWN currently known. Characterisation of the nebula's strongly energy-dependent morphology enables the particle transport mechanisms to be constrained. A dependence of the nebula extent with energy of R $\propto$ E^\alpha with \alpha = -0.29 +/- 0.04 (stat) +/- 0.05 (sys) disfavours a pure diffusion scenario for particle transport within the nebula. The total gamma-ray flux of the nebula above 1~TeV is found to be (1.12 +/- 0.03 (stat) +/- 0.25 (sys)) $\times 10^{-11}$ cm$^{-2}$ s$^{-1}$, corresponding to ~64% of the flux of the Crab Nebula. Conclusions: HESS J1825-137 is a PWN with clear energy-dependent morphology at VHE gamma-ray energies. This source is used as a laboratory to investigate particle transport within middle-aged PWNe. Deep observations of this highly spatially-extended PWN enable a spectral map of the region to be produced, providing insights into the spectral variation within the nebula.

A hotspot at a position compatible with the BL Lac object 1ES 2322-409 was serendipitously detected with H.E.S.S. during observations performed in 2004 and 2006 on the blazar PKS 2316-423. Additional data on 1ES 2322-409 were taken in 2011 and 2012, leading to a total live-time of 22.3h. Point-like very-high-energy (VHE; E>100GeV) $\gamma$-ray emission is detected from a source centred on the 1ES 2322-409 position, with an excess of 116.7 events at a significance of 6.0$\sigma$. The average VHE $\gamma$-ray spectrum is well described with a power law with a photon index $\Gamma=3.40\pm0.66_{\text{stat}}\pm0.20_{\text{sys}}$ and an integral flux $\Phi(E>200GeV) = (3.11\pm0.71_{\rm stat}\pm0.62_{\rm sys})\times10^{-12} cm^{-2} s^{-1}$, which corresponds to 1.1$\%$ of the Crab nebula flux above 200 GeV. Multi-wavelength data obtained with Fermi LAT, Swift XRT and UVOT, RXTE PCA, ATOM, and additional data from WISE, GROND and Catalina, are also used to characterise the broad-band non-thermal emission of 1ES 2322-409. The multi-wavelength behaviour indicates day-scale variability. Swift UVOT and XRT data show strong variability at longer scales. A spectral energy distribution (SED) is built from contemporaneous observations obtained around a high state identified in Swift data. A modelling of the SED is performed with a stationary homogeneous one-zone synchrotron-self-Compton (SSC) leptonic model. The redshift of the source being unknown, two plausible values were tested for the modelling. A systematic scan of the model parameters space is performed, resulting in a well-constrained combination of values providing a good description of the broad-band behaviour of 1ES 2322-409.

Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM annihilation into mono-energetic gamma-rays and ii) DM in the form of pure WIMP multiplets that, annihilating into all electroweak bosons, produce a distinctive gamma-ray spectral shape with a high-energy peak at the DM mass and a lower-energy continuum. For case i), upper limits at 95\% confidence level of about $\langle \sigma v \rangle \lesssim 3 \times 10^{-25}$ cm$^3$ s$^{-1}$ are obtained in the mass range of 400 GeV to 1 TeV. For case ii), the full spectral shape of the models is used and several excluded regions are identified, but the thermal masses of the candidates are not robustly ruled out.

Using the simultaneous Infra-Red (IR) and X-ray light curves obtained by Kalamkar et al. (2016), we perform a Fourier analysis of the IR/X-ray timing correlations of the black hole X-ray binary (BHB) GX 339-4. The resulting IR vs X-ray Fourier coherence and lag spectra are similar to those obtained in previous studies of GX 339-4 using optical light curves. In particular, above 1 Hz, the lag spectrum features an approximately constant IR lag of about 100 ms. We model simultaneously the radio to IR Spectral Energy Distribution (SED), the IR Power Spectral Density (PSD), and the coherence and lag spectra using the jet internal shock model ISHEM assuming that the fluctuations of the jet Lorentz factor are driven by the accretion flow. It turns out that most of the spectral and timing features, including the 100 ms lag, are remarkably well reproduced by this model. The 100 ms time-scale is then associated with the travel time from the accretion flow to the IR emitting zone. Our exploration of the parameter space favours a jet which is at most mildly relativistic ($\bar{\Gamma}< 3$), and a linear and positive relation between the jet Lorentz factor and X-ray light curve i.e. $\Gamma(t)-1\propto L_{X}(t)$. The presence of a strong Low Frequency Quasi Periodic Oscillation (LFQPO) in the IR light curve could be caused by jet precession driven by Lense-Thirring precession of the jet-emitting accretion flow. Our simulations confirm that this mechanism can produce an IR LFQPO similar to that observed in GX 339-4.

We report on the measurement and investigation of pulsed high-energy $\gamma$-ray emission from the Vela pulsar, PSR B0833$-$45, based on 40.3 hours of observations with the largest telescope of H.E.S.S., CT5, in monoscopic mode, and on 8 years of data obtained with the Fermi-LAT. A dedicated very-low-threshold event reconstruction and analysis pipeline was developed and, together with the CT5 telescope response model, was validated using the Fermi-LAT data as reference. A pulsed $\gamma$-ray signal at a significance level of more than $15\sigma$ is detected from the P2 peak of the Vela pulsar light curve. Of a total of 15835 events, more than 6000 lie at an energy below 20 GeV, implying a significant overlap between H.E.S.S. II-CT5 and the Fermi-LAT. While the investigation of the pulsar light curve with the LAT confirms characteristics previously known up to 20 GeV, in the tens of GeV energy range, CT5 data show a change in the pulse morphology of P2, i.e., an extreme sharpening of its trailing edge, together with the possible onset of a new component at 3.4$\sigma$ significance level. Assuming a power-law model for the P2 spectrum, an excellent agreement is found for the photon indices ($\Gamma \simeq$ 4.1) obtained with the two telescopes above 10 GeV and an upper bound of 8% is derived on the relative offset between their energy scales. Using both instruments data, it is however shown that the spectrum of P2 in the 10-100 GeV has a pronounced curvature, i.e. a confirmation of the sub-exponential cutoff form found at lower energies with the LAT. This is further supported by the weak evidence for an emission above 100 GeV obtained with CT5. In contrast, converging indications are found from both CT5 and LAT data for the emergence of a hard component above 50 GeV in the leading wing (LW2) of P2, which possibly extends beyond 100 GeV.

Centaurus A (Cen A) is the nearest radio galaxy discovered as a very-high-energy (VHE; 100 GeV-100 TeV) $\gamma$-ray source by the High Energy Stereoscopic System (H.E.S.S.). It is a faint VHE $\gamma$-ray emitter, though its VHE flux exceeds both the extrapolation from early Fermi-LAT observations as well as expectations from a (misaligned) single-zone synchrotron-self Compton (SSC) description. The latter satisfactorily reproduces the emission from Cen A at lower energies up to a few GeV. New observations with H.E.S.S., comparable in exposure time to those previously reported, were performed and eight years of Fermi-LAT data were accumulated to clarify the spectral characteristics of the $\gamma$-ray emission from the core of Cen A. The results allow us for the first time to achieve the goal of constructing a representative, contemporaneous $\gamma$-ray core spectrum of Cen A over almost five orders of magnitude in energy. Advanced analysis methods, including the template fitting method, allow detection in the VHE range of the core with a statistical significance of 12$\sigma$ on the basis of 213 hours of total exposure time. The spectrum in the energy range of 250 GeV-6 TeV is compatible with a power-law function with a photon index $\Gamma=2.52\pm0.13_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$. An updated Fermi-LAT analysis provides evidence for spectral hardening by $\Delta\Gamma\simeq0.4\pm0.1$ at $\gamma$-ray energies above $2.8^{+1.0}_{-0.6}$ GeV at a level of $4.0\sigma$. The fact that the spectrum hardens at GeV energies and extends into the VHE regime disfavour a single-zone SSC interpretation for the overall spectral energy distribution (SED) of the core and is suggestive of a new $\gamma$-ray emitting component connecting the high-energy emission above the break energy to the one observed at VHE energies.

Core collapse supernovae (CCSNe) produce fast shocks which pervade the dense circum-stellar medium (CSM) of the stellar progenitor. Cosmic rays (CRs) if accelerated at these shocks can induce the growth of electromagnetic fluctuations in the foreshock medium. In this study, using a self-similar description of the shock evolution, we calculate the growth timescales of CR-driven instabilities. We select a sample of nearby core collapse radio supernova of type II and Ib/Ic. From radio data we infer the parameters which enter in the calculation of the instability growth times. We find that extended IIb SNe shocks can trigger fast intra day instabilities, strong magnetic field amplification and CR acceleration. In particular, the non-resonant streaming instability can contribute to about 50\% of the magnetic field intensity deduced from radio data. This results in the acceleration of CRs in the range 1-10 PeV within a few days after the shock breakout. In order to produce strong magnetic field amplification and CR acceleration a fast shocks pervading a dense CSM is necessary. In that aspect IIn supernov\ae~are also good candidates. But a detailed modeling of the blast wave dynamics coupled with particle acceleration is mandatory for this class of object before providing any firm conclusions. Finally, we find that the trans-relativistic object SN 2009bb even if it produces more modest magnetic field amplification can accelerate CRs up to 2-3 PeV within 20 days after the outburst.

(Abridged) Context. NGC 253 is one of only two starburst galaxies found to emit $\gamma$-rays from hundreds of MeV to multi-TeV energies. Accurate measurements of the very-high-energy (VHE) (E $>$ 100 GeV) and high-energy (HE) (E $>$ 60 MeV) spectra are crucial to study the underlying particle accelerators and cosmic-ray interaction and transport. Aims. The measurement of the VHE $\gamma$-ray emission of NGC 253 published in 2012 by H.E.S.S. was limited by large systematic uncertainties. Here, a measurement of the $\gamma$-ray spectrum of NGC 253 is investigated in both HE and VHE $\gamma$-rays. Methods. The data of H.E.S.S. observations are reanalysed using an updated calibration and analysis chain. The $Fermi$-LAT analysis employs more than 8 years of data processed using pass 8. The cosmic-ray particle population is evaluated from the combined HE--VHE $\gamma$-ray spectrum using NAIMA. Results. The VHE $\gamma$-ray energy spectrum is best fit by a power-law with a flux normalisation of $(1.34\,\pm\,0.14^{\mathrm{stat}}\,\pm\,0.27^{\mathrm{sys}}) \times 10^{-13} \mathrm{cm^{-2} s^{-1} TeV^{-1}}$ at 1 TeV -- about 40 \% above, but compatible with the value obtained in Abramowski et al. (2012). The spectral index $\Gamma = 2.39 \pm 0.14^{\mathrm{stat}} \pm 0.25^{\mathrm{sys}}$ is slightly softer than but consistent with the previous measurement. At energies above $\sim$3 GeV the HE spectrum is consistent with a power-law ranging into the VHE part of the spectrum measured by H.E.S.S. Conclusions. Two scenarios for the starburst nucleus are tested, in which the gas in the starburst nucleus acts as a target for hadronic cosmic rays. In these two models, the level to which NGC\,253 acts as a calorimeter is estimated to a range of $f_{\rm cal} = 0.1$ to $1$ while accounting for the measurement uncertainties.

Spectral lines are among the most powerful signatures for dark matter (DM) annihilation searches in very-high-energy $\gamma$-rays. The central region of the Milky Way halo is one of the most promising targets given its large amount of DM and proximity to Earth. We report on a search for a monoenergetic spectral line from self-annihilations of DM particles in the energy range from 300 GeV to 70 TeV using a two-dimensional maximum likelihood method taking advantage of both the spectral and spatial features of signal versus background. The analysis makes use of Galactic Center (GC) observations accumulated over ten years (2004 - 2014) with the H.E.S.S. array of ground-based Cherenkov telescopes. No significant $\gamma$-ray excess above the background is found. We derive upper limits on the annihilation cross section $\langle\sigma v\rangle$ for monoenergetic DM lines at the level of $\sim4\times10^{-28}$ cm$^{3}$s$^{-1}$ at 1 TeV, assuming an Einasto DM profile for the Milky Way halo. For a DM mass of 1 TeV, they improve over the previous ones by a factor of six. The present constraints are the strongest obtained so far for DM particles in the mass range 300 GeV - 70 TeV. Ground-based $\gamma$-ray observations have reached sufficient sensitivity to explore relevant velocity-averaged cross sections for DM annihilation into two $\gamma$-ray photons at the level expected from the thermal relic density for TeV DM particles.

A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A radio SNR candidate has been identified to be a counterpart to HESS J1534-571. The TeV source is therefore classified as a SNR. For the other two sources, HESS J1614-518 and HESS J1912+101, no identifying counterparts have been found, thus they remain SNR candidates for the time being. TeV-emitting SNRs are key objects in the context of identifying the accelerators of Galactic cosmic rays. The TeV emission of the relativistic particles in the new sources is examined in view of possible leptonic and hadronic emission scenarios, taking the current multiwavelength knowledge into account.

Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RX J0852.0$-$4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and spectrum of the whole RX J0852.0$-$4622 region. The H.E.S.S. data are combined with archival data from other wavebands and interpreted in the framework of leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the direct determination of the spectral characteristics of the parent particle population in leptonic and hadronic scenarios using only GeV-TeV data. Results. An updated analysis of the H.E.S.S. data shows that the spectrum of the entire SNR connects smoothly to the high-energy spectrum measured by Fermi-LAT. The increased data set makes it possible to demonstrate that the H.E.S.S. spectrum deviates significantly from a power law and is well described by both a curved power law and a power law with an exponential cutoff at an energy of $E_\mathrm{cut} = (6.7 \pm 1.2_\mathrm{stat} \pm 1.2_\mathrm{syst})$ TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous identification of the spectral shape as a power law with an exponential cutoff. No significant evidence is found for a variation of the spectral parameters across the SNR, suggesting similar conditions of particle acceleration across the remnant. A simple modeling using one particle population to model the SNR emission demonstrates that both leptonic and hadronic emission scenarios remain plausible. It is also shown that at least a part of the shell emission is likely due to the presence of a pulsar wind nebula around PSR J0855$-$4644.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

PKS 0625-354 (z=0.055) was observed with the four H.E.S.S. telescopes in 2012 during 5.5 hours. The source was detected above an energy threshold of 200 GeV at a significance level of 6.1$\sigma$. No significant variability is found in these observations. The source is well described with a power-law spectrum with photon index $\Gamma =2.84 \pm 0.50_{stat} \pm 0.10_{syst}$ and normalization (at $E_0$=1.0 TeV) $N_0(E_0)=(0.58 \pm 0.22_{stat} \pm 0.12_{syst})\times10^{-12}$ TeV$^{-1}$cm$^{-2}$s$^{-1}$. Multi-wavelength data collected with Fermi-LAT, Swift-XRT, Swift-UVOT, ATOM and WISE are also analysed. Significant variability is observed only in the Fermi-LAT $\gamma$-ray and Swift-XRT X-ray energy bands. Having a good multi-wavelength coverage from radio to very high energy, we performed a broadband modelling from two types of emission scenarios. The results from a one zone lepto-hadronic, and a multi-zone leptonic models are compared and discussed. On the grounds of energetics, our analysis favours a leptonic multi-zone model. Models associated to the X-ray variability constraint supports previous results suggesting a BL Lac nature of PKS 0625-354, with, however, a large-scale jet structure typical of a radio galaxy.

Shell-type supernova remnants (SNRs) are considered prime candidates for the acceleration of Galactic cosmic rays (CRs) up to the knee of the CR spectrum at $\mathrm{E} \approx \mathrm{3}\times \mathrm{10}^\mathrm{15}$ eV. Our Milky Way galaxy hosts more than 350 SNRs discovered at radio wavelengths and at high energies, of which 220 fall into the H.E.S.S. Galactic Plane Survey (HGPS) region. Of those, only 50 SNRs are coincident with a H.E.S.S source and in 8 cases the very high-energy (VHE) emission is firmly identified as an SNR. The H.E.S.S. GPS provides us with a legacy for SNR population study in VHE $\gamma$-rays and we use this rich data set to extract VHE flux upper limits from all undetected SNRs. Overall, the derived flux upper limits are not in contradiction with the canonical CR paradigm. Assuming this paradigm holds true, we can constrain typical ambient density values around shell-type SNRs to $n\leq 7~\textrm{cm}^\textrm{-3}$ and electron-to-proton energy fractions above 10~TeV to $\epsilon_\textrm{ep} \leq 5\times 10^{-3}$. Furthermore, comparisons of VHE with radio luminosities in non-interacting SNRs reveal a behaviour that is in agreement with the theory of magnetic field amplification at shell-type SNRs.

We give a brief review of the origin and acceleration of cosmic rays (CRs), emphasizing the production of CRs at different stages of supernova evolution by the first-order Fermi shock acceleration mechanism. We suggest that supernovae with trans-relativistic outflows, despite being rather rare, may accelerate CRs to energies above 10$^{18}$ eV over the first year of their evolution. Supernovae in young compact clusters of massive stars, and interaction powered superluminous supernovae, may accelerate CRs well above the PeV regime. We discuss the acceleration of the bulk of the galactic CRs in isolated supernova remnants and re-acceleration of escaped CRs by the multiple shocks present in superbubbles produced by associations of OB stars. The effects of magnetic field amplification by CR driven instabilities, as well as superdiffusive CR transport, are discussed for nonthermal radiation produced by nonlinear shocks of all speeds including trans-relativistic ones.

The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments. The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources. The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document.

Context. Recently, the high-energy (HE, 0.1-100 GeV) $\gamma$-ray emission from the object LMC P3 in the Large Magellanic Cloud (LMC) has been discovered to be modulated with a 10.3-day period, making it the first extra-galactic $\gamma$-ray binary. Aims. This work aims at the detection of very-high-energy (VHE, >100 GeV) $\gamma$-ray emission and the search for modulation of the VHE signal with the orbital period of the binary system. Methods. LMC P3 has been observed with the High Energy Stereoscopic System (H.E.S.S.); the acceptance-corrected exposure time is 100 h. The data set has been folded with the known orbital period of the system in order to test for variability of the emission. Energy spectra are obtained for the orbit-averaged data set, and for the orbital phase bin around the VHE maximum. Results. VHE $\gamma$-ray emission is detected with a statistical significance of 6.4 $\sigma$. The data clearly show variability which is phase-locked to the orbital period of the system. Periodicity cannot be deduced from the H.E.S.S. data set alone. The orbit-averaged luminosity in the $1-10$ TeV energy range is $(1.4 \pm 0.2) \times 10^{35}$ erg/s. A luminosity of $(5 \pm 1) \times 10^{35}$ erg/s is reached during 20% of the orbit. HE and VHE $\gamma$-ray emissions are anti-correlated. LMC P3 is the most luminous $\gamma$-ray binary known so far.

We search for high-energy gamma-ray emission from the binary neutron star merger GW170817 with the H.E.S.S. Imaging Air Cherenkov Telescopes. The observations presented here have been obtained starting only 5.3h after GW170817. The H.E.S.S. target selection identified regions of high probability to find a counterpart of the gravitational wave event. The first of these regions contained the counterpart SSS17a that has been identified in the optical range several hours after our observations. We can therefore present the first data obtained by a ground-based pointing instrument on this object. A subsequent monitoring campaign with the H.E.S.S. telescopes extended over several days, covering timescales from 0.22 to 5.2 days and energy ranges between $270\,\mathrm{GeV}$ to $8.55\,\mathrm{TeV}$. No significant gamma-ray emission has been found. The derived upper limits on the very-high-energy gamma-ray flux for the first time constrain non-thermal, high-energy emission following the merger of a confirmed binary neutron star system.