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.

Quasi-periodic oscillations inferred during rare magnetar giant flare tails were initially interpreted as torsional oscillations of the neutron star (NS) crust, and have been more recently described as global core+crust perturbations. Similar frequencies are also present in high signal-to-noise magnetar short bursts. In magnetars, disturbances of the field are strongly coupled to the NS crust regardless of the triggering mechanism of short bursts. For low-altitude magnetospheric magnetar models of fast radio bursts (FRBs) associated with magnetar short bursts, such as the low-twist model, crustal oscillations may be associated with additional radio bursts in the encompassing short burst event (as recently suggested for SGR 1935+2154). Given the large extragalactic volume probed by wide-field radio transient facilities, this offers the prospect of studying NS crusts leveraging samples far more numerous than galactic high-energy magnetar bursts by studying statistics of sub-burst structure or clustered trains of FRBs. We sketch what may be learned and issues that may arise in such a program. Lower $l$-number eigenmodes (corresponding to FRB time intervals of $\sim5-50$ ms) are likely less susceptible than high-$l$ modes to confusion by systematic effects associated with the NS crust physics, magnetic field, and damping. They may be more promising in their utility, and also may corroborate models where FRBs arise from mature magnetars. Future observational characterization of such signals can also determine whether they can be employed as cosmological "standard oscillators" to constrain redshift, or can be used to constrain the mass of FRB-producing magnetars when reliable redshifts are available.

The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16 day period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a {\it rotational period}, an idea observationally motivated by the 6.7 hour period of the Galactic magnetar candidate, 1E 161348--5055. We explore three physical mechanisms that could lead to the creation of ultra long period magnetars: (i) enhanced spin-down due to episodic mass-loaded charged particle winds (e.g. as may accompany giant flares), (ii) angular momentum kicks from giant flares and (iii) fallback leading to long lasting accretion disks. We show that particle winds and fallback accretion can potentially lead to a sub-set of the magnetar population with ultra long periods, sufficiently long to accommodate FRB 180916 or 1E 161348--5055. If confirmed, such periods implicate magnetars in relatively mature states (ages $1-10$ kyr) and which possessed large internal magnetic fields at birth $B_{\rm int}\gtrsim 10^{16}$ G. In the low-twist magnetar model for FRBs, such long period magnetars may dominate FRB production for repeaters at lower isotropic-equivalent energies and broaden the energy distribution beyond that expected for a canonical population of magnetars which terminate their magnetic activity at shorter periods $P \lesssim 10$ s.

Fast Radio Bursts (FRBs) are millisecond-long radio pulses of extragalactic origin with peak luminosities far exceeding any Milky Way sources. The prevalent invocation for the FRB origin involves magnetars: young, magnetically powered neutron stars with the strongest magnetic fields in the Universe. A magnetar-defining signature is the emission of bright, hard X-ray bursts of sub-second duration. These occur in isolation or during a burst storm, when several hundred are observed within minutes to hours. On April 27th 2020, the Galactic magnetar SGR J1935+2154 entered an active period, emitting hundreds of X-ray bursts in a few hours. Remarkably, only one of these temporally coincided with an FRB-like radio burst. Here we report on the spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with NASA NICER and Fermi/GBM missions in their combined energy range. We demonstrate that the FRB-associated X-ray burst is very similar temporally, albeit strikingly different spectrally, from the 24 NICER/GBM bursts. If the FRB-associated burst were drawn from this magnetar burst population, its occurrence rate would be at most around 1 in 7000. This rarity combined with the unusual X-ray burst spectrum is perhaps indicative of an uncommon locale for the origin of the FRB-associated burst. We suggest that this unique event originated in quasi-polar open or closed magnetic field lines extending to high altitudes where radio emission can be generated, possibly from a collimated plasma flow.

We report on the timing and spectral properties of the soft X-ray emission from the magnetar 1E 2259+586 from January 2013, $\sim 8$ months after the detection of an anti-glitch, until September 2019, using the Neil Gehrels Swift and NICER observatories. During this time span, we detect two timing discontinuities. The first, occurring around 5 years after the April 2012 anti-glitch, is a relatively large spin-up glitch with a fractional amplitude $\Delta\nu/\nu=1.24(2)\times10^{-6}$. We find no evidence for flux enhancement or change in the spectral or pulse profile shape around the time of this glitch. This is consistent with the picture that a significant number of magnetar spin-up glitches are radiatively-quiet. Approximately 1.5 years later in April 2019, 1E 2259+586 exhibited an anti-glitch with spin-down of a fractional amplitude $\Delta\nu/\nu=-5.8(1)\times10^{-7}$; similar to the fractional change detected in 2012. We do not, however, detect any change to the pulse-profile shape or increase in the rms pulsed flux of the source, nor do we see any possible bursts from its direction around the time of the anti-glitch; all of which occurred during the 2012 event. Hence, similar to spin-up glitches, anti-glitches can occur silently. This may suggest that these phenomena originate in the neutron star interior, and that their locale and triggering mechanism do not necessarily have to be connected to the magnetosphere. Lastly, our observations suggest that the occurrence rate of spin-up and spin-down glitches is about the same in 1E 2259+586, with the former having a larger net fractional change.

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$.

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.

During a pointed 2018 NuSTAR observation, we detected a flare with a 2.2 hour duration from the magnetar 1RXS J170849.0$-$400910. The flare, which rose in $\sim25$ seconds to a maximum flux 6 times larger than the persistent emission, is highly pulsed with an rms pulsed fraction of $53\%$. The pulse profile shape consists of two peaks separated by half a rotational cycle, with a peak flux ratio of $\sim$2. The flare spectrum is thermal with an average temperature of 2.1 keV. Phase resolved spectroscopy show that the two peaks possess the same temperature, but differ in size. These observational results along with simple light curve modeling indicate that two identical antipodal spots, likely the magnetic poles, are heated simultaneously at the onset of the flare and for its full duration. Hence, the origin of the flare has to be connected to the global dipolar structure of the magnetar. This might best be achieved externally, via twists to closed magnetospheric dipolar field lines seeding bombardment of polar footpoint locales with energetic pairs. Approximately 1.86 hours following the onset of the flare, a short burst with its own 3-minute thermal tail occurred. The burst tail is also pulsating at the spin period of the source and phase-aligned with the flare profile, implying an intimate connection between the two phenomena. The burst may have been caused by a magnetic reconnection event in the same twisted dipolar field lines anchored to the surface hot spots, with subsequent return currents supplying extra heat to these polar caps.

We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). Motivated by the power-law fluence distributions of FRB 121102, we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol which directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar short bursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts (required in WT19) for individual repeaters implies the predicted burst fluence distribution is narrow, $\lesssim 3$ decades for yielding strains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and charge starvation, we obtain a death line for FRBs which segregates magnetars from the normal pulsar population, suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individual magnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-law luminosity function's low energy character depends on the population model, while the high energy index traces that of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalent energy/luminosity function peaks at $\sim10^{37}-10^{40}$ erg with a low-energy cutoff at $\sim 10^{37}$ erg. Lastly, we consider the local fluence distribution of FRBs, and find that it can constrain the subset of FRB-producing magnetar progenitors. Our model suggests that improvements in sensitivity may reveal flattening of the global FRB fluence distribution and saturation in FRB rates.

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 apply the standard radio pulsar rotating vector model to the white dwarf pulsar AR Sco's optical polarization position angle swings folded at the white dwarf's spin period as obtained by Buckley et al. (2017). Owing to the long duty cycle of spin pulsations with a good signal-to-noise ratio over the entire spin phase, in contrast to neutron star radio pulsars, we find well-constrained values for the magnetic obliquity $\alpha$ and observer viewing direction $\zeta$ with respect to the spin axis. We find $\cos\alpha=0.060^{+0.050}_{-0.053}$ and $\cos\zeta=0.49^{+0.09}_{-0.08}$, implying an orthogonal rotator with an observer angle $\zeta={60.4^\circ}^{+5.3^\circ}_{-6.0^\circ}$. This orthogonal nature of the rotator is consistent with the optical light curve consisting of two pulses per spin period, separated by $180^\circ$ in phase. Under the assumption that $\zeta \approx i$, where $i$ is the orbital inclination, and a Roche-lobe-filling companion M star, we obtain $m_{\rm WD} = 1.00^{+0.16}_{-0.10} M_\odot$ for the white dwarf mass. These polarization modeling results suggest the that non-thermal emission arises from a dipolar white dwarf magnetosphere and close to the star, with synchrotron radiation (if non-zero pitch angles can be maintained) being the plausible loss mechanism, marking AR Sco as an exceptional system for future theoretical and observational study.

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.

We show that the $\gamma$-ray pulsar observables, i.e., their total $\gamma$-ray luminosity, $L_{\gamma}$, spectral cut-off energy, $\epsilon_{\rm cut}$, stellar surface magnetic field, $B_{\star}$, and spin-down power $\dot{\mathcal{E}}$, obey a relation of the form $L_{\gamma}=f(\epsilon_{\rm cut},B_{\star},\dot{\mathcal{E}})$, which represents a 3D plane in their 4D log-space. Fitting the data of 88 pulsars of the second Fermi pulsar catalog, we show this relation to be $L_{\gamma}\propto \epsilon_{\rm cut}^{1.18\pm 0.24}B_{\star}^{0.17\pm 0.05}\dot{\mathcal{E}}^{0.41\pm 0.08}$, a pulsar fundamental plane (FP). We show that the observed FP is remarkably close to the theoretical relation $L_{\gamma}\propto \epsilon_{\rm cut}^{4/3}B_{\star}^{1/6}\dot{\mathcal{E}}^{5/12}$ obtained assuming that the pulsar $\gamma$-ray emission is due to curvature radiation by particles accelerated at the pulsar equatorial current sheet just outside the light cylinder. Interestingly, the FP seems incompatible with emission by synchrotron radiation. The corresponding scatter about the FP is $\sim 0.35$dex and can only partly be explained by the observational errors while the rest is probably due to the variation of the inclination and observer angles. We predict also that $\epsilon_{\rm cut}\propto \dot{\mathcal{E}}^{7/16}$ toward low $\dot{\mathcal{E}}$ for both young and millisecond pulsars implying that the observed death-line of $\gamma$-ray pulsars is due to $\epsilon_{\rm cut}$ dropping below the Fermi-band. Our results provide a comprehensive interpretation of the observations of $\gamma$-ray pulsars, setting requirement for successful theoretical modeling.

Marsh et al. detected radio and optical pulsations from the binary system AR Scorpii (AR Sco). This system, with an orbital period of 3.56 h, is composed of a cool, low-mass star and a white dwarf with a spin period of 1.95 min. Optical observations by Buckley et al. showed that the polarimetric emission from the white dwarf is strongly linearly polarised ( up to $\sim40\%$) with periodically changing intensities. This periodic non-thermal emission is thought to be powered by the highly magnetised ($ 5 \times 10^{8} $ G) white dwarf that is spinning down. The morphology of the polarisation signal, namely the position angle plotted against the phase angle, is similar to that seen in many radio pulsars. In this paper, we demonstrate that we can fit the traditional pulsar rotating vector model to the optical position angle. We used a Markov-chain-Monte-Carlo technique to find the best fit for the model yielding a magnetic inclination angle of $\alpha = (86.6^{+3.0}_{-2.8})^{\circ}$ and an observer angle of $\zeta = (60.5^{+5.3}_{-6.1})^{\circ}$. This modelling supports the scenario that the synchrotron emission originates above the polar caps of the white dwarf pulsar and that the latter is an orthogonal rotator.

We analyze the statistics of pulse arrival times in fast radio burst (FRB) 121102 and demonstrate that they are remarkably similar to statistics of magnetar high-energy short bursts. Motivated by this correspondence, we propose that repeating FRBs are generated during short bursts in the closed field line zone of magnetar magnetospheres via a pulsar-like emission mechanism. Crustal slippage events dislocate field line foot points, initiating intense particle acceleration and pair production, giving rise to coherent radio emission similar to that generated near pulsar polar caps. We argue that the energetics of FRB 121102 can be readily accounted for if the efficiency of the conversion of Poynting flux into coherent radio emission is $\sim10^{-4}-10^{-2}$, values consistent with empirical efficiencies of radio emission in pulsars and radio-loud magnetars. Such a mechanism could operate only in magnetars with preexisting low twist of the magnetosphere, so that the charge density in the closed zone is initially insufficient to screen the electric field provoked by the wiggling of magnetic field lines and is low enough to let $\sim 1$ GHz radio emission escape the magnetosphere, which can explain the absence of FRBs from known magnetars. The pair cascades crowd the closed flux tubes with plasma, screening the accelerating electric field, thus limiting the radio pulse duration to $\sim1$ ms. Within the framework of our model, the current dataset of the polarization angle variation in FRB 121102 suggests a magnetic obliquity $\alpha\lesssim40^\circ$ and viewing angle $\zeta$ with respect to the spin axis $\alpha<\zeta<180^\circ-\alpha$.

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.

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.

Over the last four decades, persistent and flaring emission of magnetars observed by various telescopes has provided us with a suite of light curves and spectra in soft and hard X-rays, with no emission yet detected above around 1 MeV. Attenuation of such high-energy photons by magnetic pair creation and photon splitting is expected to be active in the magnetospheres of magnetars, possibly accounting for the paucity of gamma-rays in their signals. This paper explores polarization-dependent opacities for these two QED processes in static vacuum dipole magnetospheres of highly-magnetized neutron stars, calculating attenuation lengths and determining escape energies, which are the maximum photon energies for transparency out to infinity. The numerical trajectory integral analysis in flat and curved spacetimes provides upper bounds of a few MeV or less to the visible energies for magnetars for locales proximate to the stellar surface. Photon splitting opacity alone puts constraints on the possible emission locales in their magnetospheres: regions within field loops of maximum altitudes 2-4 stellar radii are not commensurate with maximum detected energies of around 250 keV. These constraints apply not only to magnetar flares but also to their quiescent hard X-ray tail emission. An exploration of photon splitting attenuation in the context of a resonant inverse Compton scattering model for the hard X-ray tails derives distinctive phase-resolved spectroscopic and polarimetric signatures, of significant interest for future MeV-band missions such as AMEGO and e-ASTROGAM.

Enabled by the Fermi Large Area Telescope, we now know young and recycled pulsars fill the gamma-ray sky, and we are beginning to understand their emission mechanism and their distribution throughout the Galaxy. However, key questions remain: Is there a large population of pulsars near the Galactic center? Why do the most energetic pulsars shine so brightly in MeV gamma rays but not always at GeV energies? What is the source and nature of the pair plasma in pulsar magnetospheres, and what role does the polar cap accelerator play? Addressing these questions calls for a sensitive, wide-field MeV telescope, which can detect the population of MeV-peaked pulsars hinted at by Fermi and hard X-ray telescopes and characterize their spectral shape and polarization.

This White Paper explores advances in the study of Active Galaxies which will be enabled by new observing capabilities at MeV energies (hard X-rays to gamma-rays; 0.1-1000 MeV), with a focus on multi-wavelength synergies. This spectral window, covering four decades in energy, is one of the last frontiers for which we lack sensitive observations. Only the COMPTEL mission, which flew in the 1990s, has significantly probed this energy range, detecting a handful of AGN. In comparison, the currently active Fermi Gamma-ray Space Telescope, observing at the adjacent range of 0.1-100 GeV, is 100-1000 times more sensitive. This White Paper describes advances to be made in the study of sources as diverse as tidal disruption events, jetted AGN of all classes (blazars, compact steep-spectrum sources, radio galaxies and relics) as well as radio-quiet AGN, most of which would be detected for the first time in this energy regime. New and existing technologies will enable MeV observations at least 50-100 times more sensitive than COMPTEL, revealing new source populations and addressing several open questions, including the nature of the corona emission in non-jetted AGN, the precise level of the optical extragalactic background light, the accretion mode in low-luminosity AGN, and the structure and particle content of extragalactic jets.

The era of precision cosmology has revealed that about 85% of the matter in the universe is dark matter. Two well-motivated candidates are weakly interacting massive particles (WIMPs) and weakly interacting sub-eV particles (WISPs) (e.g. axions). Both WIMPs and WISPs possess distinct {\gamma}-ray signatures. Over the last decade, data taken between 50 MeV to >300 GeV by the Fermi Large Area Telescope (Fermi-LAT) have provided stringent constraints on both classes of dark matter models. Thus far, there are no conclusive detections. However, there is an intriguing {\gamma}-ray excess associated with the Galactic center that could be explained by WIMP annihilation. At lower energies, the poor angular resolution of the Fermi-LAT makes source identification challenging, inhibiting our ability to more sensitively probe both the Galactic center excess, as well as lower-mass WIMP and WISP models. Additionally, targeted WISP searches (e.g., those probing supernovae and blazars) would greatly benefit from enhanced energy resolution and polarization measurements in the MeV range. To address these issues, a new telescope that is optimized for MeV observations is needed. Such an instrument would allow us to explore new areas of dark matter parameter space and provide unprecedented access to its particle nature.

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 possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

The first multimessenger observation of a neutron star merger was independently detected in gamma-rays by Fermi-GBM and INTEGRAL SPI-ACS and gravitational waves by Advanced LIGO and Advanced Virgo. Gravitational waves are emitted from systems with accelerating quadrupole moments, and detectable sources are expected to be compact objects. Nearly all distant astrophysical gamma-ray sources are compact objects. Therefore, serendipitous observations of these two messengers will continue to uncover the sources of gravitational waves and gamma-rays, and enable multimessenger science across the Astro2020 thematic areas. This requires upgrades to the ground-based gravitational wave network and ~keV-MeV gamma-ray coverage for observations of neutron star mergers, and broadband coverage in both gravitational waves and gamma-rays to monitor other expected joint sources.

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.

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. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second 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 active galactic nuclei and gamma-ray bursts.

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.

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 number of low-mass millisecond pulsar (MSP) binaries in their rotation-powered state exhibit double-peaked X-ray orbital modulation centered at inferior pulsar conjunction. This state, which has been known to persist for years, has recently been interpreted as emission from a shock that enshrouds the pulsar. However, the pressure balance for such a configuration is a crucial unresolved issue. We consider two scenarios for pressure balance: a companion magnetosphere and stellar mass loss with gas dominance. It is found that the magnetospheric scenario requires several kilogauss poloidal fields for isobaric surfaces to enshroud the MSP as well as for the magnetosphere to remain stable if there is significant mass loss. For the gas-dominated scenario, it is necessary that the companion wind loses angular momentum prolifically as an advection or heating-dominated flow. Thermal bremsstrahlung cooling in the flow may be observable as a UV to soft X-ray component independent of orbital phase if the mass rate is high. We formulate the general requirements for shock stability against gravitational influences in the pulsar rotation-powered state for the gas-dominated scenario. We explore stabilizing mechanisms, principally irradiation feedback, which anticipates correlated shock emission and companion variability and predicts $ F_\gamma/F_{\rm X} \lesssim 14$ for the ratio of pulsar magnetospheric $\gamma$-ray to total shock soft-to-hard X-ray fluxes. This stability criterion implies an unbroken extension of X-ray power-law emission to hundreds of keV for some systems. We explore observational discriminants between the gas-dominated and magnetospheric scenarios, motivating contemporaneous radio through $\gamma$-ray monitoring of these systems.

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.

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.

(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.

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.

Hard, non-thermal, persistent pulsed X-ray emission extending between 10 keV and $\sim 150$ keV has been observed in nearly ten magnetars. For inner-magnetospheric models of such emission, resonant inverse Compton scattering of soft thermal photons by ultra-relativistic charges is the most efficient production mechanism. We present angle-dependent upscattering spectra and pulsed intensity maps for uncooled, relativistic electrons injected in inner regions of magnetar magnetospheres, calculated using collisional integrals over field loops. Our computations employ a new formulation of the QED Compton scattering cross section in strong magnetic fields that is physically correct for treating important spin-dependent effects in the cyclotron resonance, thereby producing correct photon spectra. The spectral cut-off energies are sensitive to the choices of observer viewing geometry, electron Lorentz factor, and scattering kinematics. We find that electrons with energies $\lesssim 15$ MeV will emit most of their radiation below 250 keV, consistent with inferred turnovers for magnetar hard X-ray tails. More energetic electrons still emit mostly below 1 MeV, except for viewing perspectives sampling field line tangents. Pulse profiles may be singly- or doubly-peaked dependent upon viewing geometry, emission locale, and observed energy band. Magnetic pair production and photon splitting will attenuate spectra to hard X-ray energies, suppressing signals in the Fermi-LAT band. The resonant Compton spectra are strongly polarized, suggesting that hard X-ray polarimetry instruments such as X-Calibur, or a future Compton telescope, can prove central to constraining model geometry and physics.

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.

The H.E.S.S. collaboration has discovered a new very high energy (VHE, E $>$ 0.1 TeV) $\gamma$-ray source, HESS J1741-302, located in the Galactic plane. Despite several attempts to constrain its nature, no plausible counterpart has been found so far at X-ray and MeV/GeV $\gamma$-ray energies, and the source remains unidentified. An analysis of 145-hour of observations of HESS J1741-302 at VHEs has revealed a steady and relatively weak TeV source ($\sim$1$\%$ of the Crab Nebula flux), with a spectral index of $\Gamma$ = 2.3 $\pm$ 0.2$_{\text{stat}}$ $\pm$ 0.2$_{\text{sys}}$, extending to energies up to 10 TeV without any clear signature of a cut-off. In a hadronic scenario, such a spectrum implies an object with particle acceleration up to energies of several hundred TeV. Contrary to most H.E.S.S. unidentified sources, the angular size of HESS J1741-302 is compatible with the H.E.S.S. point spread function at VHEs, with an extension constrained to be below 0.068$^{\circ}$ at a 99$\%$ confidence level. The $\gamma$-ray emission detected by H.E.S.S. can be explained both within a hadronic scenario, due to collisions of protons with energies of hundreds of TeV with dense molecular clouds, and in a leptonic scenario, as a relic pulsar wind nebula, possibly powered by the middle-aged (20 kyr) pulsar PSR B1737-30. A binary scenario, related to the compact radio source 1LC 358.266+0.038 found to be spatially coincident with the best fit position of HESS J1741-302, is also envisaged.

Pulsed non-thermal quiescent emission between 10 keV and around 150 keV has been observed in $\sim10$ magnetars. For inner magnetospheric models of such hard X-ray signals, resonant Compton upscattering of soft thermal photons from the neutron star surface is the most efficient radiative process. We present angle-dependent hard X-ray upscattering model spectra for uncooled monoenergetic relativistic electrons. The spectral cut-off energies are critically dependent on the observer viewing angles and electron Lorentz factor. We find that electrons with energies less than around 15 MeV will emit most of their radiation below 250 keV, consistent with the observed turnovers in magnetar hard X-ray tails. Moreover, electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulses phases. Our spectral computations use new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.

Non-thermal quiescent X-ray emission extending between 10 keV and around 150 keV has been seen in about 10 magnetars by RXTE, INTEGRAL, Suzaku, NuSTAR and Fermi-GBM. For inner magnetospheric models of such hard X-ray signals, inverse Compton scattering is anticipated to be the most efficient process for generating the continuum radiation, because the scattering cross section is resonant at the cyclotron frequency. We present hard X-ray upscattering spectra for uncooled monoenergetic relativistic electrons injected in inner regions of pulsar magnetospheres. These model spectra are integrated over bundles of closed field lines and obtained for different observing perspectives. The spectral turnover energies are critically dependent on the observer viewing angles and electron Lorentz factor. We find that electrons with energies less than around 15 MeV will emit most of their radiation below 250 keV, consistent with the turnovers inferred in magnetar hard X-ray tails. Electrons of higher energy still emit most of the radiation below around 1 MeV, except for quasi-equatorial emission locales for select pulse phases. Our spectral computations use a new state-of-the-art, spin-dependent formalism for the QED Compton scattering cross section in strong magnetic fields.

Index of H.E.S.S. conference proceedings to the 35th ICRC, Busan, Korea

Very high-energy gamma-rays (VHE, E>100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the universe to VHE gamma-rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extragalactic sources allows the EBL to be measured, which otherwise is very difficult to determine. An EBL-model independent measurement of the EBL SED with the H.E.S.S. array of Cherenkov telescopes is presented. It is obtained by extracting the EBL absorption signal from the reanalysis of high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is detected at a significance of 9.5 sigma, and the intensity of the EBL obtained in different spectral bands is presented together with the associated gamma-ray horizon.