• Energy Research

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.

Context. 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 bright broad-band emission observed at phases close to periastron offers a unique opportunity to study particle acceleration and radiation processes in binary systems. Observations at gamma-ray energies constrain these processes through variability and spectral characterisation studies. Aims. 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. Short-term and average emission characteristics of PSR B1259–63/LS 2883 are determined. Super-orbital variability is searched for in order to investigate possible cycle-to-cycle VHE flux changes due to different properties of the companion star’s circumstellar disc and/or the conditions under which the HE gamma-ray flares develop. Methods. Spectra and light curves were derived from observations conducted with the H.E.S.S-II array in 2014 and 2017. Phase-folded light curves are compared with the results obtained in 2004, 2007, and 2011. Fermi-LAT observations from 2010/11, 2014, and 2017 are analysed. Results. A local double-peak profile with asymmetric peaks in the VHE light curve is measured, with a flux minimum at the time of periastron tp and two peaks coinciding with the times at which the neutron star crosses the companion’s circumstellar disc (~tp ± 16 d). A high VHE gamma-ray flux is also observed at the times of the HE gamma-ray flares (~tp + 30 d) and at phases before the first disc crossing (~tp − 35 d). The spectral energy range now extends to below 200 GeV and up to ~45 TeV. Conclusions. 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 the monitoring of the source with the H.E.S.S. telescopes. This flux variability is most probably caused by the changing environmental conditions, particularly at times close to periastron passage at which the neutron star is thought to cross the circumstellar disc of the companion star twice. In contrast, the photon index remains unchanged within uncertainties for about 200 d around periastron. 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.

Context.Traditionally, TeV-γ-ray astronomy has been conducted by experiments employing proprietary data and analysis software. However, the next generation ofγ-ray instruments, such as the Cherenkov Telescope Array Observatory (CTAO), will be operated as open observatories. Alongside the data, they will also make the associated software tools available to a wider community. This necessity prompted the development of open, high-level, astronomical software customized for high-energy astrophysics.Aims.In this article, we present Gammapy, an open-source Python package for the analysis of astronomicalγ-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γ-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.Methods.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 γ-ray instruments. After the data are binned, the flux and morphology of one or moreγ-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.Results.After describing the structure of the package, we show, using publicly available gamma-ray data, the capabilities of Gammapy in multiple traditional and novelγ-ray analysis scenarios, such as spectral and spectro-morphological modeling and estimations of a spectral energy distribution and a light curve. Its flexibility and its 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.