• Energy Research

Abstract Extreme high-frequency-peaked BL Lac objects (EHBLs) are blazars that exhibit extremely energetic synchrotron emission. They also feature nonthermal gamma-ray emission whose peak lies in the very high-energy (VHE, E > 100 GeV) range, and in some sources exceeds 1 TeV: this is the case for hard-TeV EHBLs such as 1ES 0229+200. With the aim of increasing the EHBL population, 10 targets were observed with the MAGIC telescopes from 2010 to 2017, for a total of 265 hr of good-quality data. The data were complemented by coordinated Swift observations. The X-ray data analysis confirms that all but two sources are EHBLs. The sources show only a modest variability and a harder-when-brighter behavior, typical for this class of objects. At VHE gamma-rays, three new sources were detected and a hint of a signal was found for another new source. In each case, the intrinsic spectrum is compatible with the hypothesis of a hard-TeV nature of these EHBLs. The broadband spectral energy distributions (SEDs) of all sources are built and modeled in the framework of a single-zone, purely leptonic model. The VHE gamma-ray-detected sources were also interpreted with a spine–layer model and a proton synchrotron model. The three models provide a good description of the SEDs. However, the resulting parameters differ substantially in the three scenarios, in particular the magnetization parameter. This work presents the first mini catalog of VHE gamma-ray and multiwavelength observations of EHBLs.

ABSTRACT We investigate the physical nature and origin of the gamma-ray emission from the extended source HESS J1841−055 observed at TeV and GeV energies. We observed HESS J1841−055 at TeV energies for a total effective time of 43 h with the MAGIC telescopes, in 2012 and 2013. Additionally, we analysed the GeV counterpart making use of about 10 yr of Fermi-LAT data. Using both Fermi-LAT and MAGIC, we study both the spectral and energy-dependent morphology of the source for almost four decades of energy. The origin of the gamma-ray emission from this region is investigated using multiwaveband information on sources present in this region, suggested to be associated with this unidentified gamma-ray source. We find that the extended emission at GeV–TeV energies is best described by more than one source model. We also perform the first energy-dependent analysis of the HESS J1841−055 region at GeV–TeV. We find that the emission at lower energies comes from a diffuse or extended component, while the major contribution of gamma rays above 1 TeV arises from the southern part of the source. Moreover, we find that a significant curvature is present in the combined observed spectrum of MAGIC and Fermi-LAT. The first multiwavelength spectral energy distribution of this unidentified source shows that the emission at GeV–TeV energies can be well explained with both leptonic and hadronic models. For the leptonic scenario, bremsstrahlung is the dominant emission compared to inverse Compton. On the other hand, for the hadronic model, gamma-ray resulting from the decay of neutral pions (π0) can explain the observed spectrum. The presence of dense molecular clouds overlapping with HESS J1841−055 makes both bremsstrahlung and π0-decay processes the dominant emission mechanisms for the source.

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.