OT 081 is a well-known, luminous blazar that is remarkably variable in many energy bands. We present the first broadband study of the source, which includes very high energy (VHE, E > 100 GeV) gamma-ray data taken by the MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov telescopes) and H.E.S.S. (High Energy Stereoscopic System) imaging Cherenkov telescopes. The discovery of VHE gamma-ray emission happened during a high state of gamma-ray activity in July 2016, observed by many instruments from radio to VHE gamma-rays. We identify four states of activity of the source, one of which includes VHE gamma-ray emission. Variability in the VHE domain is found on daily time-scales. The intrinsic VHE spectrum can be described by a power law with index 3 . 27 +/- 0.44(stat) +/- 0.15(sys) (MAGIC) and 3 . 39 +/- 0.58(stat) +/- 0.64(sys) (H.E.S.S.) in the energy range of 55-300 and 120-500 GeV, respectively. The broadband emission cannot be successfully reproduced by a simple one-zone synchrotron self-Compton model. Instead, an additional external Compton component is required. We test a lepto-hadronic model that reproduces the data set well and a proton-synchrotron-dominated model that requires an extreme proton luminosity. Emission models that are able to successfully represent the data place the emitting region well outside of the broad-line region to a location at which the radiative environment is dominated by the infrared thermal radiation field of the dusty torus. In the scenario described by this flaring activity, the source appears to be a flat spectrum radio quasar (FSRQ), in contrast with past categorizations. This suggests that the source can be considered to be a transitional blazar, intermediate between BL Lac and FSRQ objects.

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.

Recent observations of a few young pulsar wind nebulae (PWNe) have revealed their morphologies in some detail. Given the availability of spatio-spectral-temporal data, we use our multi-zone (1D) leptonic emission code to model the PWNe associated with G29.7-0.3 (Kes 75) and G21.5-0.9 (G21.5) and obtain (by-eye) constraints on additional model parameters compared to spectral-only modelling. Kes 75 is a Galactic composite supernova remnant (SNR) with an embedded pulsar, PSR J1846-0258. X-ray studies reveal rapid expansion of Kes 75 over the past two decades. PWN G21.5 is also a composite SNR, powered by PSR J1833-1034. For Kes 75, we study a sudden plasma bulk speed increase that may be due to the magnetar-like outbursts of the central pulsar. An increase of a few percent in this speed does not result in any significant change in the model outputs. For G21.5, we investigate different diffusion coefficients and pulsar spin-down braking indices. We can reproduce the broadband spectra and X-ray surface brightness profiles for both PWNe, and the expansion rate, flux over different epochs, and X-ray photon index vs epoch and central radius for Kes 75 quite well. The latter three features are also investigated for G21.5. Despite obtaining reasonable fits overall, some discrepancies remain, pointing to further model revision. We find similar values to overlapping parameters between our 1D code and those of an independent 0D dynamical code (TIDE). Future work will incorporate spatial data from various energy wavebands to improve model constraints.

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.

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.

Context. To date, three pulsars have been firmly detected by imaging atmospheric Cherenkov telescopes (IACTs). Two of them reached the TeV energy range, challenging models of very high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the large-sized telescopes, which will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims. In this work we study the Crab pulsar emission with the LST-1, improving upon and complementing the results from other telescopes. Crab pulsar observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods. We analyzed a total of similar to 103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles of less than 50 degrees. To characterize the Crab pulsar emission over a broader energy range, a new analysis of the Fermi/LAT data, including similar to 14 years of observations, was also performed. Results. The Crab pulsar phaseogram, long-term light curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for the first peak and up to 700 GeV for the second peak The pulsed emission is detected with a significance level of 15.2 sigma. The two characteristic emission peaks of the Crab pulsar are clearly detected (> 10 sigma), as is the so-called bridge emission between them (5.7 sigma). We find that both peaks are described well by power laws, with spectral indices of similar to 3.44 and similar to 3.03, respectively. The joint analysis of Fermi/LAT and LST-1 data shows a good agreement between the two instruments in their overlapping energy range. The detailed results obtained from the first observations of the Crab pulsar with the LST-1 show the potential that CTAO will have to study this type of source.

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.

Context. Imaging atmospheric Cherenkov telescopes (IACTs) are used to observe very high-energy photons from the ground. Gamma rays are indirectly detected through the Cherenkov light emitted by the air showers they induce. The new generation of experiments, in particular the Cherenkov Telescope Array Observatory (CTAO), sets ambitious goals for discoveries of new gamma-ray sources and precise measurements of the already discovered ones. To achieve these goals, both hardware and data analysis must employ cuttingedge techniques. This also applies to the LST-1, the first IACT built for the CTAO, which is currently taking data on the Canary island of La Palma. Aims. This paper introduces a new event reconstruction technique for IACT data, aiming to improve the image reconstruction quality and the discrimination between the signal and the background from misidentified hadrons and electrons. Methods. The technique models the development of the extensive air shower signal, recorded as a waveform per pixel, seen by CTAO telescopes’ cameras. Model parameters are subsequently passed to random forest regressors and classifiers to extract information on the primary particle. Results. The new reconstruction was applied to simulated data and to data from observations of the Crab Nebula performed by the LST-1. The event reconstruction method presented here shows promising performance improvements. The angular and energy resolution, and the sensitivity, are improved by 10 to 20% over most of the energy range. At low energy, improvements reach up to 22%, 47%, and 50%, respectively. A future extension of the method to stereoscopic analysis for telescope arrays will be the next important step.

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.

The Tunka Advanced Instrument for Gamma- and cosmic-ray Astronomy (TAIGA) is a multicomponent experiment for the measurement of TeV to PeV gamma- and cosmic rays. Our goal is to establish a novel hybrid direct air shower technique, sufficient to access the energy domain of the long-sought Pevatrons. The hybrid air Cherenkov light detection technique combines the strengths of the HiSCORE shower front sampling array, and two similar to 4 m class, similar to 9.6 degrees field of view Imaging Air Cherenkov Telescopes (IACTs). The HiSCORE array provides good angular and shower core position resolution, while the IACTs provide the image shape and orientation for gamma-hadron separation. In future, an additional muon detector will be used for hadron tagging at >= 100 TeV energies. Here, only data from the first IACT of the TAIGA experiment are used. A randomforest algorithm was trained using Monte Carlo (MC) simulations and real data, and applied to 85 sigma of selected observational data tracking the Crab Nebula at a mean zenith angle of 33.5 degrees, resulting in a threshold energy of 6 TeV for this data set. The analysis was performed using the gammapy package. A total of 163.5 excess events were detected, with a statistical significance of 8.5 sigma. The observed spectrum of the Crab Nebula is best fit with a power law above 6 TeV with a flux normalization of (3.20 +/- 0.42) 10(-10) TeV-1 cm(-2) s(-1)at a reference energy of $13\, \mathrm{TeV}$ and a spectral index of -2.74 +/- 0.16.

The Large Magellanic Cloud (LMC) presents a unique environment for pulsar population studies due to its distinct star formation characteristics and proximity to the Milky Way. As part of the TRAPUM (TRAnsients and PUlsars with MeerKAT) Large Survey Project, we are using the core array of the MeerKAT radio telescope (MeerKAT) to conduct a targeted search of the LMC for radio pulsars at L-band frequencies, 856-1712 MHz. The excellent sensitivity of MeerKAT, coupled with a 2-h integration time, makes the survey 3 times more sensitive than previous LMC radio pulsar surveys. We report the results from the initial four survey pointings which have resulted in the discovery of seven new radio pulsars, increasing the LMC radio pulsar population by 30 per cent. The pulse periods of these new pulsars range from 278 to 1690 ms, and the highest dispersion measure is 254.20 pc cm(-3). We searched for, but did not find any significant pulsed radio emission in a beam centred on the SN 1987A remnant, establishing an upper limit of 6.3 mu Jy on its minimum flux density at 1400 MHz.

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems colocated with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to ≤11 known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with six undetected in radio. Overall, ≥236 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power E ̇ decreases to its observed minimum near 10 33 erg s −1 , approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties, and fit results to guide and be compared with modeling results.

LMC P3 is a gamma-ray binary comprising of an unconfirmed compact object and an O-star and is located in the Large Magellanic Cloud. Initially discovered in Fermi-LAT data, it shows an orbital period of 10.3 days. H.E.S.S. has reported the detected VHE gamma-ray emission during only 20% of the orbit, between orbital phases 0.2 and 0.4, which roughly corresponds to the inferior conjunction of the compact object. H.E.S.S. has continued the observations of this object since then. Here we will present new results obtained with a much deeper data set. The new data allow a more precise measurement of the location of the VHE gamma-ray peak along the orbit of the system to be made. We will interpret these results with respect to emission and absorption mechanisms in gamma-ray binary systems.

The Galactic Centre (GC) region is a highly interesting region for very high energy gamma-ray studies due to its proximity and diverse sources. It is also a unique place for Dark Matter (DM) searches, since we expect a large amount of DM in this region and it is nearby. Currently a new detector is under development to observe the GC region, called the Southern Wide field-of-view Gamma-ray Observatory (SWGO). This instrument will be the first water Cherenkov detector located in the Southern Hemisphere sensitive to >100 GeV gamma rays. In our work, we are going to present a simulation of the GC gamma-ray mission as seen by a SWGO-like observatory. To do so we are using the models published by the Cherenkov Telescope Array (CTA) Collaboration together with simulated Instrument Response Functions (IRFs) for testing purposes. This also allows us to predict the expected sensitivity of SWGO to WIMP DM annihilations using a template based method. A sensitivity study during the design phase of the observatory is important in order to choose the best detector design to get outstanding physics results.

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.

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.

We present a novel approach to the determination of the pointing of Imaging Atmospheric Cherenkov Telescopes (IACTs) using the trajectories of the stars in their camera’s field of view. The method starts with the reconstruction of the star positions from the Cherenkov camera data, taking into account the point spread function of the telescope, to achieve a satisfying reconstruction accuracy of the pointing position. A simultaneous fit of all reconstructed star trajectories is then performed with the orthogonal distance regression (ODR) method. ODR allows us to correctly include the star position uncertainties and use the time as an independent variable. Having the time as an independent variable in the fit makes it better suitable for various star trajectories. This method can be applied to any IACT and requires neither specific hardware nor interface or special data-taking mode. In this paper, we use the Large-Sized Telescope (LST) data to validate it as a useful tool to improve the determination of the pointing direction during regular data taking. The simulation studies show that the accuracy and precision of the method are comparable with the design requirements on the pointing accuracy of the LST (≤14″). With the typical LST event acquisition rate of 10 kHz, the method can achieve up to 50 Hz pointing monitoring rate, compared to 𝒪(1) Hz achievable with standard techniques. The application of the method to the LST prototype (LST-1) commissioning data shows the stable pointing performance of the telescope.

Fast Radio Bursts (FRBs) are highly energetic, extremely short-lived bursts of radio flashes. Despite extensive research, the exact cause of these outbursts remains a mystery. One of the most accredited models suggests that they originate from highly magnetized and rapidly spinning neutron stars known as magnetars. The high luminosity, short duration, and high dispersion measure of these events suggest they result from extreme, high-energy astrophysical processes of extragalactic origin. The number of detected FRBs, including repeating ones, has grown rapidly in recent years. Except for FRB20200428, that is associated to the galactic magnetar SGR1925+2154, no multi-wavelength counterparts to any FRB has been detected yet. The High Energy Stereoscopic System (H.E.S.S.) telescope has developed a program to uncover the nature of these mysterious events by searching for their gamma-ray counterparts. This contribution provides an overview of the searches for FRB sources conducted by H.E.S.S., including follow-up observations and simultaneous multi-wavelength campaigns with radio and X-ray observatories

The binary system Eta Carinae is a unique laboratory to study particle acceleration up to very high energies (VHE) under a wide range of conditions. Particles are thought to be accelerated at shocks forming in the wind collision region. Eta Carinae has been firmly established as a source of high energy gamma-rays in Fermi-LAT data over several orbits. With its highly eccentric orbit lasting 5.5 years, the periastron passage of the two stars is extremely close. This provides an opportunity to constrain the acceleration and absorption mechanisms of the system. Eta Carinae was detected above 200 GeV by the H.E.S.S. telescopes (H.E.S.S. Collaboration, 2020) based on data before and after the 2014 periastron. Unfortunately the 2014 periastron itself could not be observed by H.E.S.S. due to visibility constraints. Hence the 2020 periastron was the first periastron passage visible for the full 5 telescope H.E.S.S. array and was therefore monitored with an extensive observation campaign spanning the phase range from 0.97 to 1.05. We report the detection of a VHE signal from Eta Carinae during the 2020 periastron and describe its spectral properties together with simultaneous Fermi-LAT data. Together with previous and follow up observations, for the first time a VHE light curve spanning a full orbit is presented.