In this article, we present Gammapy, an open-source Python package for the analysis of astronomical $\gamma$-ray data, and illustrate the functionalities of its first long-term-support release, version 1.0. Built on the modern Python scientific ecosystem, Gammapy provides a uniform platform for reducing and modeling data from different $\gamma$-ray instruments for many analysis scenarios. Gammapy complies with several well-established data conventions in high-energy astrophysics, providing serialized data products that are interoperable with other software packages. Starting from event lists and instrument response functions, Gammapy provides functionalities to reduce these data by binning them in energy and sky coordinates. Several techniques for background estimation are implemented in the package to handle the residual hadronic background affecting $\gamma$-ray instruments. After the data are binned, the flux and morphology of one or more $\gamma$-ray sources can be estimated using Poisson maximum likelihood fitting and assuming a variety of spectral, temporal, and spatial models. Estimation of flux points, likelihood profiles, and light curves is also supported. After describing the structure of the package, we show, using publicly available $\gamma$-ray data, the capabilities of Gammapy in multiple traditional and novel $\gamma$-ray analysis scenarios, such as spectral and spectro-morphological modeling and estimations of a spectral energy distribution and a light curve. Its flexibility and power are displayed in a final multi-instrument example, where datasets from different instruments, at different stages of data reduction, are simultaneously fitted with an astrophysical flux model.

The flat spectrum radio quasar (FSRQ) PKS 1510−089 is known for its complex multiwavelength behaviour and it is one of only a few FSRQs detected in very-high-energy (VHE,E > 100 GeV)γrays. The VHEγ-ray observations with H.E.S.S. and MAGIC in late May and early June 2016 resulted in the detection of an unprecedented flare, which revealed, for the first time, VHEγ-ray intranight variability for this source. While a common variability timescale of 1.5 h has been found, there is a significant deviation near the end of the flare, with a timescale of ∼20 min marking the cessation of the event. The peak flux is nearly two orders of magnitude above the low-level emission. For the first time, a curvature was detected in the VHEγ-ray spectrum of PKS 1510–089, which can be fully explained by the absorption on the part of the extragalactic background light. OpticalR-band observations with ATOM revealed a counterpart of theγ-ray flare, even though the detailed flux evolution differs from the VHEγ-ray light curve. Interestingly, a steep flux decrease was observed at the same time as the cessation of the VHEγ-ray flare. In the high-energy (HE,E > 100 MeV)γ-ray band, only a moderate flux increase was observed withFermi-LAT, while the HEγ-ray spectrum significantly hardens up to a photon index of 1.6. A search for broad-line region (BLR) absorption features in theγ-ray spectrum indicates that the emission region is located outside of the BLR. Radio very-long-baseline interferometry observations reveal a fast-moving knot interacting with a standing jet feature around the time of the flare. As the standing feature is located ∼50 pc from the black hole, the emission region of the flare may have been located at a significant distance from the black hole. If this is indeed a true correlation, the VHEγrays must have been produced far down in the jet, where turbulent plasma crosses a standing shock.

The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both datasets, the point spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the gamma-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. dataset. Remaining differences in the overall gamma-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the gamma-ray sky between WCD and IACT techniques.

Magnetar hyperflares are the most plausible explanation for fast radio bursts (FRBs)-enigmatic powerful radio pulses with durations of several milliseconds and high brightness temperatures. The first observational evidence for this scenario was obtained in 2020 April when an FRB was detected from the direction of the Galactic magnetar and soft gamma-ray repeater SGR 1935+2154. The FRB was preceded by two gamma-ray outburst alerts by the BAT instrument aboard the Swift satellite, which triggered follow-up observations by the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observed SGR 1935+2154 for 2 hr on 2020 April 28. The observations are coincident with X-ray bursts from the magnetar detected by INTEGRAL and Fermi-GBM, thus providing the first very high energy gamma-ray observations of a magnetar in a flaring state. High-quality data acquired during these follow-up observations allow us to perform a search for short-time transients. No significant signal at energies E > 0.6 TeV is found, and upper limits on the persistent and transient emission are derived. We here present the analysis of these observations and discuss the obtained results and prospects of the H.E.S.S. follow-up program for soft gamma-ray repeaters.

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

Theγ-ray pulsar PSR B1706 −44 and the adjacent supernova remnant (SNR) candidate G 343. 1−2. were observed by H.E.S.S. during a dedicated observation campaign in 2007. As a result of this observatio n campaign, a new source of very-high-energy (VHE; E > 100 GeV)γ-ray emission, HESS J1708 −443, was detected with a statistical significance of 7 σ, although no significant point-like emission was detected a t the position of the energetic pulsar itself. In this paper, the morphologic al and spectral analyses of the newly-discovered TeV source are presented. The centroid of HESS J1708 −443 is considerably o ffset from the pulsar and located near the apparent center of th e SNR, atαJ2000 = 17h08m11s ± 17s and δJ2000 = −44◦20′ ± 4′. The source is found to be significantly more extended than th e H.E.S.S. point spread function ( ∼0.1◦), with an intrinsic Gaussian width of 0.29 ◦ ± 0.04◦. Its integral flux between 1 and 10 TeV is ∼ 3.8× 10−12 ph cm−2 s−1, equivalent to 17% of the Crab Nebula flux in the same energy range. The measured energy spectrum is wel l-fit by a power law with a relatively hard photon index Γ = 2.0± 0.1stat± 0.2sys. Additional multi-wavelength data, including 330 MHz VLA ob servations, were used to investigate the VHE γ-ray source’s possible associations with the pulsar wind nebula of PSR B1706 −44 and/or with the complex radio structure of the partial shell-typ e SNR G 343.1 −2.3.

We search for an indirect signal of dark matter through very high-energy γ rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as γ rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dark matter dominated objects with well-measured kinematics and small uncertainties on their dark matter distribution profiles. In 2018, the five-telescopes of the high energy stereoscopic system observed the dwarf irregular galaxy WLM for 18 hours. We present the first analysis based on data obtained from an imaging atmospheric Cherenkov telescope for this subclass of dwarf galaxy. As we do not observe any significant excess in the direction of WLM, we interpret the result in terms of constraints on the velocity-weighted * Corresponding authors. cross section for dark matter pair annihilation h σ v i as a function of the dark matter particle mass for various continuum channels, as well as the prompt γγ emission. For the τ þ τ − channel, the limits reach a h σ v i value of about 4 × 10 − 22 cm 3 s − 1 for a dark matter particle mass of 1 TeV. For the prompt γγ channel, the upper limit reaches a h σ v i value of about 5 × 10 − 24 cm 3 s − 1 for a mass of 370 GeV. These limits represent an improvement of up to a factor 200, with respect to previous results for the dwarf irregular galaxies for TeV dark matter search. DOI: 10.1103/PhysRevD.103.102002

We search for an indirect signal of dark matter through very high-energy gamma rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as gamma rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dark matter dominated objects with well measured kinematics and small uncertainties on their dark matter distribution profiles. In 2018, the H.E.S.S. five-telescope array observed the dwarf irregular galaxy WLM for 18 hours. We present the first analysis based on data obtained from an imaging atmospheric Cherenkov telescope for this subclass of dwarf galaxy. As we do not observe any significant excess in the direction of WLM, we interpret the result in terms of constraints on the velocity-weighted cross section for dark matter pair annihilation as a function of the dark matter particle mass for various continuum channels as well as the prompt gamma-gamma emission. For the τ^+τ^- channel the limits reach a ⟨σ v ⟩ value of about 4× 10^-22 cm3s-1 for a dark matter particle mass of 1 TeV. For the prompt gamma-gamma channel, the upper limit reaches a ⟨σ v ⟩ value of about 5 ×10^-24 cm3s-1 for a mass of 370 GeV. These limits represent an improvement of up to a factor 200 with respect to previous results for the dwarf irregular galaxies for TeV dark matter search.

Dwarf spheroidal galaxy satellites of the Milky Way are prime targets for indirect detection of dark matter with gamma rays due to their proximity, high dark matter content and absence of non-thermal emission processes. Recently, the Dark Energy Survey (DES) revealed the existence of new ultra-faint dwarf spheroidal galaxies in the southern-hemisphere sky, therefore ideally located for ground-based observations with the imaging atmospheric Cherenkov telescope array H.E.S.S. We present a search for very-high-energy ($E\gtrsim100$ GeV) gamma-ray emission using H.E.S.S. observations carried out recently towards Reticulum II, Tucana II, Tucana III, Tucana IV and Grus II satellites. No significant very-high-energy gamma-ray excess is found from the observations on any individual object nor in the combined analysis of all the datasets. Using the most recent modeling of the dark matter distribution in the dwarf galaxy halo, we compute for the first time on DES satellites individual and combined constraints from Cherenkov telescope observations on the annihilation cross section of dark matter particles in the form of Weakly Interacting Massive Particles. The combined 95% C.L. observed upper limits reach $\langle \sigma v \rangle \simeq 1 \times 10^{-23}$ cm$^3$s$^{-1}$ in the $W^+W^-$ channel and $4 \times 10^{-26}$ cm$^3$s$^{-1}$ in the $\gamma\gamma$ channels for a dark matter mass of 1.5 TeV. The H.E.S.S. constraints well complement the results from Fermi-LAT, HAWC, MAGIC and VERITAS and are currently the most stringent in the $\gamma\gamma$ channels in the multi-GeV/multi-TeV mass range.

Here we report the results of the first ever contemporaneous multi-wavelength observation campaign on the BL Lac object PKS 2155−304 involvingSwift,NuSTAR,Fermi-LAT, and H.E.S.S. The use of these instruments allows us to cover a broad energy range, which is important for disentangling the different radiative mechanisms. The source, observed from June 2013 to October 2013, was found in a low flux state with respect to previous observations but exhibited highly significant flux variability in the X-rays. The high-energy end of the synchrotron spectrum can be traced up to 40 keV without significant contamination by high-energy emission. A one-zone synchrotron self-Compton model was used to reproduce the broadband flux of the source for all the observations presented here but failed for previous observations made in April 2013. A lepto-hadronic solution was then explored to explain these earlier observational results.

Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) γ-ray emitters. η Car is the most prominent member of this object class and is confirmed to emit phase-locked HE γ 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) γ-ray emission from η Car around the last periastron passage in 2014 with the ground-based High Energy Stereoscopic System (H.E.S.S.). Methods. The region around η Car was observed with H.E.S.S. between orbital phase p = 0.78−1.10, with a closer sampling at p ≈ 0.95 and p ≈ 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 η 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 γ-ray emission from the direction of η 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 γ rays agree within statistical and systematic errors before and after periastron. The γ-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.

The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short γ-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.

The unidentified very-high-energy (VHE;E> 0.1 TeV)γ-ray source, HESS J1826−130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steadyγ-ray flux from HESS J1826−130, which appears extended with a half-width of 0.21° ± 0.02stat° ± 0.05sys°. The source spectrum is best fit with either a power-law function with a spectral index Γ = 1.78 ± 0.10stat± 0.20sysand an exponential cut-off at 15.2−3.2+5.5TeV, or a broken power-law with Γ1= 1.96 ± 0.06stat± 0.20sys, Γ2= 3.59 ± 0.69stat± 0.20sysfor energies below and aboveEbr= 11.2 ± 2.7 TeV, respectively. The VHE flux from HESS J1826−130 is contaminated by the extended emission of the bright, nearby pulsar wind nebula, HESS J1825−137, particularly at the low end of the energy spectrum. Leptonic scenarios for the origin of HESS J1826−130 VHE emission related to PSR J1826−1256 are confronted by our spectral and morphological analysis. In a hadronic framework, taking into account the properties of dense gas regions surrounding HESS J1826−130, the source spectrum would imply an astrophysical object capable of accelerating the parent particle population up to ≳200 TeV. Our results are also discussed in a multiwavelength context, accounting for both the presence of nearby supernova remnants, molecular clouds, and counterparts detected in radio, X-rays, and TeV energies.