Rafal Moderski

Context . Very high-energy (VHE) γ -ray emission detected from the central molecular zone (CMZ) hints at diffusive propagation of cosmic rays (CRs) injected continuously near the Galactic center (GC). Aims . Using H.E.S.S. VHE γ -ray observations, we aim to construct a multicomponent description of the region in order to derive the spatial and spectral distributions of CRs and discuss their potential origin. Methods . We rely on a spectro-morphological analysis, in which spectral and spatial parameters of different components are fit simultaneously, and on a three-dimensional description of the CMZ and CR distributions. Results . The GC ridge emission is well reproduced with the assumed gas distribution model weighted by a 1 /r ^α CR density profile with α = 1.10 ± 0.05 _stat ± 0.10 _syst . We detect no significant spectral variations across the CMZ, which supports a continuous injection from the GC. We report a significant (>3 σ ) curvature in the GC ridge γ -ray spectrum, implying a parent proton break (or cutoff energy) of E _b ∼ 15–50 TeV ( E _cut ∼ 60–160 TeV). We also derive the best-fit spectral parameters of the other sources in the region and report a significant cutoff at E ∼ 5 TeV in the spectrum of the pulsar wind nebula G0.9+0.1. Finally we constrain the position of the injection site to be close to the GC (ruling out the Arches and Quintuplet clusters as major contributors to the emission). We discuss possible origins for the GC ridge emission, potential mechanisms for causing the observed curvature, and a possible contribution of the gas content to the emission of HESS J1745–290, and derive qualitative constraints on the line-of-sight positions of the clouds shaping the CMZ.

The propagation of very-high-energy (VHE; E _γ ≥ 100 GeV) gamma rays from extragalactic sources is affected by interactions with photons of the extragalactic background light (EBL), resulting in pair production that attenuates the intrinsic gamma-ray flux. This interaction renders the Universe increasingly opaque to VHE photons at high energies and redshifts. New physics scenarios involving axion-like particles (ALPs) could modify this expected optical depth. In particular, ALPs with masses m _a ∼ 10 eV can decay into two photons over cosmological timescales, thereby contributing to the diffuse EBL. If such ALPs constitute a significant fraction of the dark matter density, their decay would enhance the EBL intensity and consequently increase the gamma-ray optical depth. In this study, we investigate this scenario using a large sample of gamma-ray spectra observed with the High Energy Stereoscopic System (H.E.S.S.) and the Fermi Large Area Telescope. We model the contribution of decaying ALPs to the EBL and assess their impact on the spectra of blazars across redshifts. By comparing these observations with standard EBL models, we place constraints on the properties of heavy ALPs — specifically their mass and photon coupling — and evaluate their viability as a dark matter candidate capable of modifying the gamma-ray transparency of the Universe. From the combined analysis, and under the assumption that ALPs constitute the entire dark matter density, we derive 95% confidence exclusion limits on the photon-ALP coupling down to m _aγ ∼ 7×10 ^-12 GeV ^-1 for masses m _a ∼ 15 eV. These constraints are competitive with existing astrophysical bounds and provide complementary sensitivity to other techniques, closing a previously unconstrained region of parameter space in the m _a ∼ 2.5–20 eV range.

The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB 170817A, heralded the birth of multimessenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), which is ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV–TeV GW follow-up strategies for CTAO using a multistep simulation pipeline and to estimate the expected rate of joint GW–GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV–TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star–black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope.

At the center of our Galaxy lies the bright γ -ray point-like source HESS J1745-290, which is compatible in position with Sgr A ^⋆ , although an association between the two remains uncertain. Using data obtained between 2004 and 2019 with the High Energy Stereoscopic System (H.E.S.S.) on the Galactic center region, we studied the variability of HESS J1745-290 over 353 hours of observations collected over 16 years, representing the largest dataset gathered yet on this region at TeV energies. We performed a 3D maximum-likelihood analysis of the central source and the diffuse γ -ray emission in the Galactic center region. This analysis allowed us to extract the spectral and morphological intrinsic behavior of the two components. By performing this analysis on an annual basis, we derived the light curve of HESS J1745-290 and the diffuse emission over the past 16 years. The 3D maximum-likelihood analysis method allowed us to separate the central source from the overlapping diffuse emission, enabling a recalibration of the former by the latter and alleviating some of the systematic effects. We find no long-term or yearly variability. We also provide an estimate of the sensitivity of H.E.S.S. to variation of this specific source over 16 years. We rule out any yearly γ -ray flux variation of this source larger than 30%, as well as any linear flux variation exceeding 30% over this time period.

Context. Recent observational efforts using imaging atmospheric Cherenkov telescopes (IACTs) have led to firm detections of very-high-energy (VHE) signals from bright gamma-ray bursts (GRBs), often at moderate redshifts. Aims. This work presents 15 years of H.E.S.S. GRB observations and examines their implications through population comparisons and selected modelling cases. Methods. GRBs are a key science target of the High Energy Stereoscopic System (H.E.S.S.). With a low-energy threshold (≲100 GeV) and rapid repointing capabilities, H.E.S.S. can begin follow-up observations within tens of seconds after a GRB trigger, covering the late prompt or early afterglow phases. Results. We report GRB follow-up observations with H.E.S.S. from 2004 to 2019, which resulted in no significant VHE signals (aside from the detections of GRB 180720B and GRB 190829A). The resulting upper limits comprise the largest set available for GRBs at VHE. Conclusions. A subset of bursts with favourable conditions were selected for X-ray analysis and emission modelling. Population studies were performed to compare detected and non-detected GRBs. The results indicate that VHE-detected GRBs are not a distinct population, but tend to feature luminous X-ray emission and favourable redshift and observing conditions. This highlights the potential of next-generation IACTs such as the Cherenkov Telescope Array Observatory (CTAO), whose lower energy threshold will enhance the detection of fainter and more distant GRBs.

Dwarf spheroidal galaxies (dSphs) are excellent targets for indirect dark matter (DM) searches using gamma-ray telescopes because they are thought to have high DM content and a low astrophysical background. The sensitivity of these searches is improved by combining the observations of dSphs made by different gamma-ray telescopes. We present the results of a combined search by the most sensitive currently operating gamma-ray telescopes, namely: the satellite-borne Fermi -LAT telescope; the ground-based imaging atmospheric Cherenkov telescope arrays H.E.S.S., MAGIC, and VERITAS; and the HAWC water Cherenkov detector. Individual datasets were analyzed using a common statistical approach. Results were subsequently combined via a global joint likelihood analysis. We obtain constraints on the velocity-weighted cross section 〈 σv 〉 for DM self-annihilation as a function of the DM particle mass. This five-instrument combination allows the derivation of up to 2-3 times more constraining upper limits on 〈 σv 〉 than the individual results over a wide mass range spanning from 5 GeV to 100 TeV. Depending on the DM content modeling, the 95% confidence level observed limits reach 1.5×10 ^-24 cm ³ s ^-1 and 3.2×10 ^-25 cm ³ s ^-1 , respectively, in the τ ⁺ τ ^- annihilation channel for a DM mass of 2 TeV.

HESS J0852-463 is a TeV γ-ray source located in the Galactic plane. The region consists of a supernova remnant (SNR, RX J0852.0-4622) with a shell-like morphology, commonly referred to as Vela Junior, and a pulsar denoted PSR J0855-4644. Pulsars are among the most efficient leptonic accelerators in our Galaxy, making this region particularly interesting to study. We utilise the most recent data taken by the High Energy Stereoscopic System (H.E.S.S.), to investigate any γ-ray emission associated with the pulsar in this region, PSR J0855-4644. We applied a full forward folding method on the H.E.S.S. data. Utilising 3D modelling techniques, we evaluated the TeV γ-ray emission towards the various components of this complex system. The distinct energy-dependent morphology observed in our data motivates further investigation of this source. We resolved the emission in the Vela Junior region into various components, several of which correspond to the SNR itself. In particular, we find a new extended component which is coincident with the position of PSR J0855-4644. The spectrum follows a power-law distribution with a best-fit index of ΓE = 1.81 \pm 0.07stat which differs from the properties of the surrounding γ-ray emission of the Vela Junior SNR. A one-zone leptonic joint fit between the X-rays (from XMM-Newton) and γ-rays (from H.E.S.S.) leads to a lower limit on the magnetic field of 1.6μG and a spectral index of α = 1.88 \pm 0.01, in line with expectations of pulsar wind nebulae (PWNe). In this paper, we report the first detection of the PWN of PSR J0855-4644 at TeV energies with the H.E.S.S. experiment at a significance of 12.2σ. This is attributed to the advanced techniques of the 3D analysis. Based on the pulsar's characteristics, its PWN is consistent with the known TeV PWNe population in the Galaxy.

Context. PKS 0346-27 is a low synchrotron peaked blazar at redshift 0.991. The very high energy (VHE; E > 100 GeV) spectra of blazars are always affected by γγ absorption by the extragalactic background light (EBL), and subsequently no blazars have been detected in VHE γ -rays at redshifts exceeding 1. Aims. This is the goal of a target-of-opportunity (ToO) programme by H.E.S.S.: to observe flaring high-redshift ( z ≳ 1) blazars. Importantly, extending the redshift range of VHE-detected blazars to z ≳ 1 will yield insights into the cosmological evolution of both the VHE blazar population and the EBL. Methods. We report H.E.S.S. ToO and multi-wavelength observations of the blazar PKS 0346−27. We analysed and modelled the H.E.S.S. data together with simultaneous data from Fermi -LAT, Swift (XRT and UVOT), using single-zone leptonic and hadronic models. Results. PKS 0346-27 was detected by H.E.S.S. at a significance of 6.3 σ during one night on 3 November 2021, while for other nights before and after this day, upper limits on the VHE flux have been determined. No evidence for intra-night γ -ray variability has been found. A flare in high-energy ( E > 100 MeV) γ -rays detected by Fermi -LAT preceded the H.E.S.S. detection by 2 days. A fit with a single-zone emission model to the contemporaneous spectral energy distribution during the detection night was possible with a proton-synchrotron-dominated hadronic model, requiring a proton-kinetic-energy-dominated jet power temporarily exceeding the source’s Eddington limit, although alternative (e.g. multi-zone) models cannot be ruled out. A one-zone leptonic model is, in principle, also able to fit the flare-state spectral energy distribution. However, it requires implausible parameter choices, in particular, extreme Doppler and bulk Lorentz factors of ≳80.

Context. Microquasars have emerged as promising candidates to explain the cosmic-ray flux at petaelectronvolt energies. LHAASO observations revealed V4641 Sgr as the most extreme example so far. Its gamma-ray spectrum extends up to 800 TeV, which requires particles with multi-PeV energy. The TeV emission is highly extended, which challenges expectations given the reported low-inclination angle of the V4641 Sgr jets. Aims. We spatially and spectrally resolved the gamma-ray emission from V4641 Sgr and investigated the particle acceleration in the system. Methods. Using ≈100 h of H.E.S.S. data, we performed a spectro-morphological study of the gamma-ray emission around V4641 Sgr. We employed HI and dedicated CO observations of the region to infer the target material for cosmic-ray interactions. Results. We detected multi-TeV emission around V4641 Sgr with a high significance. The emission region is elongated, and its major and minor axes are 0.34° ±0.01 _syst ^{±0.04 _stat} and 0.06° ±0.01 _syst ^{±0.01 _stat} , respectively. We found a power-law spectrum with an index ≈1.8, and together with results from other gamma-ray instruments, this reveals a spectral energy distribution (SED) that peaks at energies of ≈100 TeV for the first time. We found indications (3 σ ) of a two-component morphology, with indistinguishable spectral properties. The position of V4641 Sgr is inconsistent with the best-fit position of the single-component model and with the dip between the two components. We found no significant evidence of an energy-dependent morphology. No dense gas was found at any distance towards V4641 Sgr, which places an upper limit of n _gas ≲ 0.2 cm ⁻³ within the gamma-ray emission region. Conclusions. The peak of the SED at ≈100 TeV identifies V4641 Sgr as a candidate cosmic-ray accelerator beyond the so-called knee. The absence of dense target gas places stringent energetic constraints on hadronic interpretations, however. The H.E.S.S. measurement requires an unusually hard (≈1.5) spectral index for the protons. A leptonic scenario faces fewer obstacles if the particle transport is fast enough to avoid losses and to reproduce the observed energy-independent morphology. The absence of bright X-ray emission across the gamma-ray emission region requires a magnetic field strength ≲3 μG, however. Our findings favour a leptonic origin of the gamma-ray emission. This conclusion does not exclude hadron acceleration in the V4641 Sgr system.