The Cherenkov Telescope Array Observatory (CTAO) will include telescopes of three different sizes, the smallest of which are the Small-Sized Telescopes (SSTs). In particular, the SSTs will be installed at the southern site of CTAO, on the Chilean Andes, and will cover the highest energy range of CTAO (up to similar to 300 TeV). The SSTs are developed by an international consortium of institutes that will provide them as an in-kind contribution to CTAO. The optical design of the SSTs is based on a Schwarzschild-Couder-like dual-mirror polynomial configuration, with a primary aperture of 4.3m diameter. They are equipped with a focal plane camera based on SiPM detectors covering a field of view of similar to 9 degrees. The preliminary design of the SST telescopes was evaluated and approved during the Product Review (PR) organised with CTAO in February 2023. The SST project is now going through a consolidation phase leading to the finalisation and submission of the final design to the Critical Design Review (CDR), expected to occur late 2024, after which the production and construction of the telescopes will begin leading to a delivery of the telescopes to CTAO southern site starting at the end of 2025-early 2026 onward. In this contribution we will present the progress of the SST programme, including the results of the PDR, the consolidation phase of the project and the plan up to the on-site integration of the telescopes.

SS 433 is a microquasar, a stellar binary system that launches collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.) and found an energy-dependent shift in the apparent position of the gamma-ray emission from the parsec-scale jets. These observations trace the energetic electron population and indicate that inverse Compton scattering is the emission mechanism of the gamma rays. Our modeling of the energy-dependent gamma-ray morphology constrains the location of particle acceleration and requires an abrupt deceleration of the jet flow. We infer the presence of shocks on either side of the binary system, at distances of 25 to 30 parsecs, and that self-collimation of the precessing jets forms the shocks, which then efficiently accelerate electrons.

Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the other two, but is significantly more distant. If the gamma-ray luminosities of the remnants of Tycho's and Kepler's SNe are similar, then the latter is expected to be one of the faintest gamma-ray sources within reach of the current generation IACT arrays. Here we report evidence at a statistical level of 4.6 sigma for a VHE signal from the remnant of Kepler's SN based on deep observations by the High Energy Stereoscopic System (H.E.S.S.) with an exposure of 152 hours. The measured integral flux above an energy of 226 GeV is ~0.3% of the flux of the Crab Nebula. The spectral energy distribution (SED) reveals a gamma-ray emitting component connecting the VHE emission observed with H.E.S.S. to the emission observed at GeV energies with Fermi-LAT. The overall SED is similar to that of the remnant of Tycho's SN, possibly indicating the same non-thermal emission processes acting in both these young remnants of thermonuclear SNe.

The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling.

The radio galaxy M87 is a variable very-high energy (VHE) gamma-ray source, exhibiting three major flares reported in 2005, 2008, and 2010. Despite extensive studies, the origin of the VHE gamma-ray emission is yet to be understood. In this study, we investigate the VHE gamma-ray spectrum of M87 during states of high gamma-ray activity, utilizing 20.2$\,$ hours the H.E.S.S. observations. Our findings indicate a preference for a curved spectrum, characterized by a log-parabola model with extra-galactic background light (EBL) model above 0.3$\,$TeV at the 4$\sigma$ level, compared to a power-law spectrum with EBL. We investigate the degeneracy between the absorption feature and the EBL normalization and derive upper limits on EBL models mainly sensitive in the wavelength range 12.4$\,$$\mu$m - 40$\,$$\mu$m.

Most gamma-ray detected active galactic nuclei are blazars with one of their relativistic jets pointing towards the Earth. Only a few objects belong to the class of radio galaxies or misaligned blazars. Here, we investigate the nature of the object PKS 0625-354, its gamma-ray flux and spectral variability and its broad-band spectral emission with observations from H.E.S.S., Fermi-LAT, Swift-XRT, and UVOT taken in November 2018. The H.E.S.S. light curve above 200 GeV shows an outburst in the first night of observations followed by a declining flux with a halving time scale of 5.9 h. The gamma gamma-opacity constrains the upper limit of the angle between the jet and the line of sight to similar to 10 degrees. The broad-band spectral energy distribution shows two humps and can be well fitted with a single-zone synchrotron self Compton emission model. We conclude that PKS 0625-354, as an object showing clear features of both blazars and radio galaxies, can be classified as an intermediate active galactic nuclei. Multi-wavelength studies of such intermediate objects exhibiting features of both blazars and radio galaxies are sparse but crucial for the understanding of the broad-band emission of gamma-ray detected active galactic nuclei in general.

Approximately one hundred sources of very-high-energy (VHE) gamma rays are known in the Milky Way. A survey of the entire Galactic Plane in the energy range from a few tens of GeV to a few hundred TeV has been proposed as a Key Science Project for the upcoming Cherenkov Telescope Array Observatory (CTAO). This article presents the status of the studies towards the Galactic Plane Survey (GPS). We build and make publicly available a sky model that combines data from observations of known gamma-ray emitters with state-of-the-art physically-driven models of synthetic populations of the main classes of established Galactic VHE sources, as well as of interstellar emission from cosmic-ray interactions in the Milky Way. We also perform an optimisation of the observation strategy. We use the improved sky model and observation strategy to simulate GPS data that are analysed using the methods and software tools under development for real data. We show that the GPS has the potential to increase the number of known Galactic VHE emitters by almost a factor of five. This corresponds to the detection of more than two hundred pulsar wind nebulae and a few tens of supernova remnants at average integral fluxes one order of magnitude lower than in the existing sample above 1 TeV, therefore opening the possibility to perform unprecedented population studies. The GPS also has the potential to provide new VHE detections of binary systems and pulsars, and to detect bright PeVatrons. Furthermore, the GPS will constitute a pathfinder for deeper follow-up observations of these source classes. Finally, we show that we can extract from GPS data an estimate of the contribution to diffuse emission from unresolved sources, and that there are good prospects of detecting interstellar emission and statistically distinguishing different scenarios. (Abridged)

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 G AMMAPY , 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.

PSR B1259–63/LS 2883 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 yr period, around an O9.5Ve star (LS 2883). At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission, for which the temporal and spectral properties of this emission are, for now, poorly understood. In this regard, very high-energy (VHE) emission is especially useful to study and constrain radiation processes and particle acceleration in the system. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of approximately 100 h of data taken over five months, from t p − 24 days to t p + 127 days around the system’s 2021 periastron passage (where t p is the time of periastron). We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent overall with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve – although sparsely sampled in this time period – we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index Γ = 2.65 ± 0.04 stat ± 0.04 sys , a value consistent with the spectral index obtained from the analysis of data collected with H.E.S.S. during the previous observations of the source. We report spectral variability with a difference of ΔΓ = 0.56 ± 0.18 stat ± 0.10 sys at 95% confidence intervals, between sub-periods of the 2021 dataset. We also detail our investigation into X-ray/TeV and GeV/TeV flux correlations in the 2021 periastron passage. We find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after t p + 25 days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from ∼ t p − 23 days to ∼ t p + 126 days. This suggests that the GeV and TeV emission originate from different electron populations.

Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $\alpha_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $\alpha_{\rm CRp}$ down to about $\Delta\alpha_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $\tau_{\chi}>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario.

Many models of dark matter (DM) are now widely considered and probed intensively with accelerators, underground detectors, and astrophysical experiments. Among the various approaches, high-energy astrophysical observations are extremely useful to complement laboratory searches for some DM candidates. In the near future, the Cherenkov Telescope Array (CTA) should enable us to access much heavier weakly interacting massive particles, as well as a broad range of other DM candidates. In this talk, we describe DM searches with CTA.

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce very high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy γ -ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the γ -ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The γ -ray luminosity above 0.5 TeV of both sources is (2–3) × 10 35 erg s −1 . This exceeds by more than a factor of 2 the luminosity of HESS J1646−458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the γ -ray emission from each source is extended with a significance of >3 σ and a Gaussian width of about 30 pc. For 30 Dor C, a connection between the γ -ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the γ -ray signal from R136 are discussed.

The measurement of very-high-energy cosmic-ray electrons is intrinsically difficult due to the very steep electron spectrum with low fluxes and an enormous background of hadronic cosmic rays. The large collection areas needed for such a measurement can be provided by ground-based imaging atmospheric Cherenkov telescopes. The High Energy Stereoscopic System (H.E.S.S.) has performed the first ground-based cosmic-ray electron measurement and thereby extended the measured range of the spectrum to several TeV. Here the H.E.S.S. measurement is presented, as well as an extension of the H.E.S.S. spectrum towards lower energies. At these energies, H.E.S.S. can probe recent ATIC measurements, which have been interpreted in terms of dark matter scenarios.

Context. Among the blazar class, extreme blazars have exceptionally hard intrinsic X-ray/TeV spectra, and extreme peak energies in their spectral energy distribution (SED). Observational evidence suggests that the non-thermal emission from extreme blazars is typically non-variable. All these unique features present a challenging case for blazar emission models, especially regarding those sources with hard TeV spectra. Aims. We aim to explore the X-ray and GeV observational features of a variety of extreme blazars, including extreme-TeV, extreme-synchrotron (extreme-Syn), and regular high-frequency-peaked BL Lac objects (HBLs). Furthermore, we aim to test the applicability of various blazar emission models that could explain the very hard TeV spectra. Methods. We conducted a detailed spectral analysis of X-ray data collected with AstroSat and Swift-XRT, along with quasi-simultaneous γ-ray data from Fermi-LAT, for five sources: 1ES 0120+340, RGB J0710+591, 1ES 1101−232, 1ES 1741+196, and 1ES 2322−409. We took three approaches to modelling the SEDs: (1) a steady-state one-zone synchrotron-self-Compton (SSC) code, (2) another leptonic scenario of co-accelerated electrons and protons on multiple shocks applied to the extreme-TeV sources only (e–p co-acceleration scenario), and (3) a one-zone hadro-leptonic (ONEHALE) code. The latter code is used twice to explain the γ-ray emission process: proton synchrotron and synchrotron emission of secondary pairs. Results. Our X-ray analysis provides well-constrained estimates of the synchrotron peak energies for both 1ES0120+340 and 1ES1741+196. These findings categorise these latter objects as extreme-synchrotron sources, as they consistently exhibit peak energies above 1 keV in different flux states. The multi-epoch X-ray and GeV data reveal spectral and flux variabilities in RGB J0710+591 and 1ES 1741+196, even on timescales of days to weeks. As anticipated, the one-zone SSC model adequately reproduces the SEDs of regular HBLs but encounters difficulties in explaining the hardest TeV emission. Hadronic models offer a reasonable fit to the hard TeV spectrum, though with the trade-off of requiring extreme jet powers. On the other hand, the lepto-hadronic scenario faces additional challenges in fitting the GeV spectra of extreme-TeV sources. Finally, the e–p co-acceleration scenario naturally accounts for the observed hard electron distributions and effectively matches the hardest TeV spectrum of RGB J0710+591 and 1ES 1101−232.

In 2021 July, PKS 1510−089 exhibited a significant flux drop in the high-energy γ-ray (by a factor 10) and optical (by a factor 5) bands and remained in this low state throughout 2022. Similarly, the optical polarization in the source vanished, resulting in the optical spectrum being fully explained through the steady flux of the accretion disk and the broad-line region. Unlike the aforementioned bands, the very-high-energy γ-ray and X-ray fluxes did not exhibit a significant flux drop from year to year. This suggests that the steady-state very-high-energy γ-ray and X-ray fluxes originate from a different emission region than the vanished parts of the high-energy γ-ray and optical jet fluxes. The latter component has disappeared through either a swing of the jet away from the line of sight or a significant drop in the photon production efficiency of the jet close to the black hole. Either change could become visible in high-resolution radio images.

The origin of the gamma-ray emission from M87 is currently a matter of debate. This work aims to localize the VHE (100 GeV-100 TeV) gamma-ray emission from M87 and probe a potential extended hadronic emission component in the inner Virgo Cluster. The search for a steady and extended gamma-ray signal around M87 can constrain the cosmic-ray energy density and the pressure exerted by the cosmic rays onto the intra-cluster medium, and allow us to investigate the role of the cosmic rays in the active galactic nucleus feedback as a heating mechanism in the Virgo Cluster. H.E.S.S. telescopes are sensitive to VHE gamma rays and have been utilized to observe M87 since 2004. We utilized a Bayesian block analysis to identify M87 emission states with H.E.S.S. observations from 2004 until 2021, dividing them into low, intermediate, and high states. Because of the causality argument, an extended ($\gtrsim$kpc) signal is allowed only in steady emission states. Hence, we fitted the morphology of the 120h low state data and found no significant gamma-ray extension. Therefore, we derived for the low state an upper limit of 58"(corresponding to $\approx$4.6kpc) in the extension of a single-component morphological model described by a rotationally symmetric 2D Gaussian model at 99.7% confidence level. Our results exclude the radio lobes ($\approx$30 kpc) as the principal component of the VHE gamma-ray emission from the low state of M87. The gamma-ray emission is compatible with a single emission region at the radio core of M87. These results, with the help of two multiple-component models, constrain the maximum cosmic-ray to thermal pressure ratio $X_{{CR,max.}}$$\lesssim$$0.32$ and the total energy in cosmic-ray protons (CRp) to $U_{CR}$$\lesssim$5$\times10^{58}$ erg in the inner 20kpc of the Virgo Cluster for an assumed CRp power-law distribution in momentum with spectral index $\alpha_{p}$=2.1.

We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2 degrees away from the best-fit position of the IceCube neutrino event 211208A. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution, and the gamma-ray data from Fermi-LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both measurements provide strong constraints on leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed spectral energy distribution (SED). The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models, and it would allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. A numerical lepto-hadronic model with external target photons reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power.

GRB221009A is the brightest gamma-ray burst ever detected. To probe the very-high-energy (VHE, $>$\!100 GeV) emission, the High Energy Stereoscopic System (H.E.S.S.) began observations 53 hours after the triggering event, when the brightness of the moonlight no longer precluded observations. We derive differential and integral upper limits using H.E.S.S. data from the third, fourth, and ninth nights after the initial GRB detection, after applying atmospheric corrections. The combined observations yield an integral energy flux upper limit of $\Phi_\mathrm{UL}^{95\%} = 9.7 \times 10^{-12}~\mathrm{erg\,cm^{-2}\,s^{-1}}$ above $E_\mathrm{thr} = 650$ GeV. The constraints derived from the H.E.S.S. observations complement the available multiwavelength data. The radio to X-ray data are consistent with synchrotron emission from a single electron population, with the peak in the SED occurring above the X-ray band. Compared to the VHE-bright GRB190829A, the upper limits for GRB221009A imply a smaller gamma-ray to X-ray flux ratio in the afterglow. Even in the absence of a detection, the H.E.S.S. upper limits thus contribute to the multiwavelength picture of GRB221009A, effectively ruling out an IC dominated scenario.

Context. HESS J1809$-$193 is an unassociated very-high-energy $\gamma$-ray source located on the Galactic plane. While it has been connected to the nebula of the energetic pulsar PSR J1809$-$1917, supernova remnants and molecular clouds present in the vicinity also constitute possible associations. Recently, the detection of $\gamma$-ray emission up to energies of $\sim$100 TeV with the HAWC observatory has led to renewed interest in HESS J1809$-$193. Aims. We aim to understand the origin of the $\gamma$-ray emission of HESS J1809$-$193. Methods. We analysed 93.2 h of data taken on HESS J1809$-$193 above 0.27 TeV with the High Energy Stereoscopic System (H.E.S.S.), using a multi-component, three-dimensional likelihood analysis. In addition, we provide a new analysis of 12.5 yr of Fermi-LAT data above 1 GeV within the region of HESS J1809$-$193. The obtained results are interpreted in a time-dependent modelling framework. Results. For the first time, we were able to resolve the emission detected with H.E.S.S. into two components: an extended component that exhibits a spectral cut-off at $\sim$13 TeV, and a compact component that is located close to PSR J1809$-$1917 and shows no clear spectral cut-off. The Fermi-LAT analysis also revealed extended $\gamma$-ray emission, on scales similar to that of the extended H.E.S.S. component. Conclusions. Our modelling indicates that based on its spectrum and spatial extent, the extended H.E.S.S. component is likely caused by inverse Compton emission from old electrons that form a halo around the pulsar wind nebula. The compact component could be connected to either the pulsar wind nebula or the supernova remnant and molecular clouds. Due to its comparatively steep spectrum, modelling the Fermi-LAT emission together with the H.E.S.S. components is not straightforward. (abridged)