• Energy Research

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.

The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy γ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula’s γ-ray emission between 1 GeV and ∼100 TeV, that is, over five orders of magnitude in energy. Using the open-source software package GAMMAPY, we combined 11.4 yr of data from the Fermi Large Area Telescope and 80 h of High Energy Stereoscopic System (H.E.S.S.) data at the event level and provide a measurement of the spatial extension of the nebula and its energy spectrum. We find evidence for a shrinking of the nebula with increasing γ-ray energy. Furthermore, we fitted several phenomenological models to the measured data, finding that none of them can fully describe the spatial extension and the spectral energy distribution at the same time. Especially the extension measured at TeV energies appears too large when compared to the X-ray emission. Our measurements probe the structure of the magnetic field between the pulsar wind termination shock and the dust torus, and we conclude that the magnetic field strength decreases with increasing distance from the pulsar. We complement our study with a careful assessment of systematic uncertainties.

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.