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 HAWC Observatory collected 6 yr of extensive data, providing an ideal platform for long-term monitoring of blazars in the very high energy (VHE) band, without bias toward specific flux states. HAWC continuously monitors blazar activity at TeV energies, focusing on sources with a redshift of z ≤ 0.3, based on the Third Fermi-LAT Catalog of High-Energy sources. We specifically focused our analysis on Mrk 421 and Mrk 501, as they are the brightest blazars observed by the HAWC Observatory. With a data set of 2143 days, this work significantly extends the monitoring previously published, which was based on 511 days of observation. By utilizing HAWC data for the VHE γ -ray emission in the 300 GeV–100 TeV energy range, in conjunction with Swift-XRT data for the 0.3–10 keV X-ray emission, we aim to explore potential correlations between these two bands. For Mrk 501, we found evidence of a long-term correlation. Additionally, we identified a period in the light curve where the flux was very low for more than 2 yr. On the other hand, our analysis of Mrk 421 measured a strong linear correlation for quasi-simultaneous observations collected by HAWC and Swift-XRT. This result is consistent with a linear dependence and a multiple-zone synchrotron self-Compton model to explain the X-ray and γ -ray emission. Finally, as suggested by previous findings, we confirm a harder-when-brighter behavior in the spectral evolution of the flux properties for Mrk 421. These findings contribute to the understanding of blazar emissions and their underlying mechanisms.

We report the total energy spectrum of cosmic rays in the energy interval from 10 TeV to 1 PeV, which lies in the energy region where both direct and indirect cosmic ray experiments overlap. The all-particle spectrum was obtained from an unfolding analysis of 5.3 years of data collected with the High Altitude Water Cherenkov (HAWC) observatory for zenith angles smaller than or equal to 35 degrees. The study was carried out in the framework of the QGSJET-II-04 hadronic interaction model. The measured spectrum confirms the presence of a knee-like feature at tens of TeV. In our analysis, the position of this softening is found at 40.2 +/- 1.0(stat.)+6.2 -6.4(sys.) TeV. The measured spectral indices before and after the break are gamma 1 = -2.53 +/- 0.01(stat.)+0.04 -0.05(sys.) and gamma 2 = -2.71 +/- 0.01(stat.)+0.03-0.04(sys.), respectively.

Supernova remnants are one potential source class considered a PeVatron (i.e., capable of accelerating cosmic rays above PeV energies). The shock fronts produced after the explosion of the supernova are ideal regions for particle acceleration. IC 443 is a supernova remnant that has been studied extensively at different wavelengths. Using 2966 days of gamma-ray data from the High Altitude Water Cherenkov (HAWC) observatory, we study the emission of IC 443 with the objective of finding signatures of cosmic-ray acceleration at the PeV scale. Using a maximum likelihood method, we find a point source located at (alpha = 94 .degrees 42, delta = 22 .degrees 35) that we associate with IC 443. The measured spectrum is a simple power law with an index of -3.14 +/- 0.18, which is consistent with previous TeV observations. Although we cannot confirm that IC 443 is a hadronic PeVatron, we do not find any sign that the spectrum has a cutoff at tens of TeV energies, with the spectrum extending to similar to 30 TeV. Furthermore, we also find a new extended component in the region whose emission is described by a simple power law with an index of -2.49 +/- 0.08 and which we call HAWC J0615+2213. While we show evidence that this new source might be a new TeV halo, we defer a detailed analysis of this new source to another publication.

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.

Axion-like particles (ALPs) are pseudo-Nambu–Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 h of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the radio galaxy in the center of the cluster: NGC 1275. By modeling the magnetic field surrounding this target, we searched for spectral indications of ALP presence. Despite finding no statistical evidence of ALP signatures, we were able to exclude ALP models in the sub-micro electronvolt range. Our analysis improved upon previous work by calculating the full likelihood and statistical coverage for all considered models across the parameter space. Consequently, we achieved the most stringent limits to date for ALP masses around 50 neV, with cross sections down to gaγ=3×10−12 GeV−1.

Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C is modelled together with multi-wavelength data within a synchrotron and synchrotron-self Compton (SSC) scenario. We find that SSC can explain the broadband data well from the optical to the very-high-energy band. For the late-time radio data, a different component is needed to account for the observed emission. Differently from previous GRBs detected in the very-high-energy range, the model for GRB~201216C strongly favors a wind-like medium. The model parameters have values similar to those found in past studies of the afterglows of GRBs detected up to GeV energies.

We have developed a new set of Application-Specific Integrated Circuits (ASICs) of the TARGET family (CTC and CT5TEA), designed for the readout of signals from photosensors in cameras of Imaging Atmospheric Cherenkov Telescopes (IACTs) for ground-based gamma-ray astronomy. We present the performance and design details. Both ASICs feature 16 channels, with CTC being a Switched-Capacitor Array (SCA) sampler at 0.5 to 1 GSa/s with a 16,384 sample deep storage buffer, including the functionality to digitize full waveforms at arbitrary times. CT5TEA is its companion trigger ASIC (though may be used on its own), which provides trigger information for the analog sum of four (and 16) adjacent channels. Since sampling and triggering takes place in two separate ASICs, the noise due to interference from the SCA is suppressed, and allows a minimal trigger threshold of <= 2.5 mV (0.74 photo electrons (p.e.)) with a trigger noise of <= 0.5 mV (0.15 p.e.). For CTC, a maximal input voltage range from -0.5 V up to 1.7 V is achieved with an effective bit range of > 11.6 bits and a baseline noise of 0.7 mV. The cross-talk improved to <= 1% over the whole -3 dB bandwidth of 220 MHz and even down to 0.2% for 1.5 V pulses of 10 ns width. Not only is the performance presented, but a temperature- stable calibration routine for pulse mode operation is introduced and validated. The resolution is found to be similar to 2.5% at 33.7 mV (10 p.e.) and <= 0.3% at 337 mV (100 p.e.) with an integrated non-linearity of < 1.6 mV. Developed for the Small-Sized Telescope (SST) and Schwarzschild-Couder Telescope (SCT) cameras of the Cherenkov Telescope Array Observatory (CTAO), CTC and CT5TEA are deployed for both prototypes and shall be integrated into the final versions.

The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $\nu_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $\nu_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $\nu_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model.

The High-Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory, located on the side of the Sierra Negra volcano in Mexico, has been fully operational since 2015. The HAWC collaboration has recently significantly improved their extensive air shower reconstruction algorithms, which has notably advanced the observatory performance. The energy resolution for primary gamma rays with energies below 1 TeV was improved by including a noise-suppression algorithm. Corrections have also been made to systematic errors in direction fitting related to the detector and shower plane inclinations, O(0.degrees 1) biases in highly inclined showers, and enhancements to the core reconstruction. The angular resolution for gamma rays approaching the HAWC array from large zenith angles (>37 degrees) has improved by a factor of 4 at the highest energies (>70 TeV) as compared to previous reconstructions. The inclusion of a lateral distribution function fit to the extensive air shower footprint on the array to separate gamma-ray primaries from cosmic-ray ones based on the resulting chi(2) values improved the background rejection performance at all inclinations. At large zenith angles, the improvement in significance is a factor of 4 compared to previous HAWC publications. These enhancements have been verified by observing the Crab Nebula, which is an overhead source for the HAWC Observatory. We show that the sensitivity to Crab-like point sources (E-2.63) with locations overhead to 30 degrees zenith is comparable to or less than 10% of the Crab Nebula's flux between 2 and 50 TeV. Thanks to these improvements, HAWC can now detect more sources, including the Galactic center.

The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

The decay or annihilation of dark matter particles may produce a steady flux of very-high-energy gamma rays detectable above the diffuse background. Nearby clusters of galaxies provide excellent targets to search for the signatures of particle dark matter interactions. In particular, the Virgo cluster spans several degrees across the sky and can be efficiently probed with a wide field-of-view instrument. The High Altitude Water Cherenkov (HAWC) observatory, due to its wide field of view and sensitivity to gamma rays at an energy scale of 300 GeV--100 TeV is well-suited for this search. Using 2141 days of data, we search for gamma-ray emission from the Virgo cluster, assuming well-motivated dark matter sub-structure models. Our results provide some of the strongest constraints on the decay lifetime of dark matter for masses above 10 TeV.

PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required.

This paper investigates the origin of the gamma-ray emission from MGRO J1908+06 in the GeV-TeV energy band. By analyzing the data collected by the Fermi Large Area Telescope, the Very Energetic Radiation Imaging Telescope Array System, and High Altitude Water Cherenkov, with the addition of spectral data previously reported by LHAASO, a multiwavelength study of the morphological and spectral features of MGRO J1908+06 provides insight into the origin of the gamma-ray emission. The mechanism behind the bright TeV emission is studied by constraining the magnetic field strength, the source age, and the distance through detailed broadband modeling. Both spectral shape and energy-dependent morphology support the scenario that inverse Compton emission of an evolved pulsar wind nebula associated with PSR J1907+0602 is responsible for the MGRO J1908+06 gamma-ray emission with a best-fit true age of T = 22 +/- 9 kyr and a magnetic field of B = 5.4 +/- 0.8 mu G, assuming the distance to the pulsar d(PSR) = 3.2 kpc.

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)

We report an observation of ultrahigh-energy (UHE) gamma rays from the Galactic center (GC) region, using 7 yr of data collected by the High-Altitude Water Cherenkov (HAWC) Observatory. The HAWC data are best described as a point-like source (HAWC J1746-2856) with a power-law spectrum (dN dE = phi(E/26 TeV)(gamma)), where gamma = -2.88 +/- 0.15(stat) - 0.1(sys) and phi = 1.5 x 10(-15) (TeV cm(2) s)(-1) +/- 0.3(stat) (+0.08sys)(-0.13sys) extending from 6 to 114 TeV. We find no evidence of a spectral cutoff up to 100 TeV using HAWC data. Two known point-like gamma-ray sources are spatially coincident with the HAWC gamma-ray excess: Sgr A* (HESS J1745-290) and the Arc (HESS J1746-285). We subtract the known flux contribution of these point sources from the measured flux of HAWC J1746-2856 to exclude their contamination and show that the excess observed by HAWC remains significant (>5 sigma), with the spectrum extending to >100 TeV. Our result supports that these detected UHE gamma rays can originate via hadronic interaction of PeV cosmic-ray protons with the dense ambient gas and confirms the presence of a proton PeVatron at the GC.

Very-high-energy (0.1-100 TeV) gamma-ray emission was observed in HAWC data from the lobes of the microquasar SS 433, making them the first set of astrophysical jets that were resolved at TeV energies. In this work, we update the analysis of SS 433 using 2,565 days of data from the High Altitude Water Cherenkov (HAWC) observatory. Our analysis reports the detection of a point-like source in the east lobe at a significance of $6.6\,\sigma$ and in the west lobe at a significance of $8.2\,\sigma$. For each jet lobe, we localize the gamma-ray emission and identify a best-fit position. The locations are close to the X-ray emission sites "e1" and "w1" for the east and west lobes, respectively. We analyze the spectral energy distributions and find that the energy spectra of the lobes are consistent with a simple power-law $\text{d}N/\text{d}E\propto E^{\alpha}$ with $\alpha = -2.44^{+0.13+0.04}_{-0.12-0.04}$ and $\alpha = -2.35^{+0.12+0.03}_{-0.11-0.03}$ for the east and west lobes, respectively. The maximum energy of photons from the east and west lobes reaches 56 TeV and 123 TeV, respectively. We compare our observations to various models and conclude that the very-high-energy gamma-ray emission can be produced by a population of electrons that were efficiently accelerated.

Galactic cosmic rays (GCRs) are charged particles that reach the heliosphere almost isotropically in a wide energy range. In the inner heliosphere, the GCR flux is modulated by solar activity so that only energetic GCRs reach the lower layers of the solar atmosphere. In this work, we propose that high-energy GCRs can be used to explore the solar magnetic fields at low coronal altitudes. We used GCR data collected by the High-Altitude Water Cherenkov observatory to construct maps of GCR flux coming from the Sun’s sky direction and studied the observed GCR deficit, known as Sun shadow (SS), over a 6 yr period (2016–2021) with a time cadence of 27.3 days. We confirm that the SS is correlated with sunspot number, but we focus on the relationship between the photospheric solar magnetic field measured at different heliolatitudes and the relative GCR deficit at different energies. We found a linear relationship between the relative deficit of GCRs represented by the depth of the SS and the solar magnetic field. This relationship is evident in the observed energy range of 2.5–226 TeV, but is strongest in the range of 12.4 33.4 TeV, which implies that this is the best energy range to study the evolution of magnetic fields in the low solar atmosphere.

HESS J1809-193 is an unidentified TeV source, first detected by the High Energy Stereoscopic System (H.E.S.S.) collaboration. The emission originates in a source-rich region that includes several supernova remnants (SNRs) and pulsars including SNR G11.1+0.1, SNR G11.0-0.0, and the young radio pulsar PSR J1809-1917. Originally classified as a pulsar wind nebula candidate, recent studies show the peak of the TeV region overlapping with a system of molecular clouds. This resulted in the revision of the original leptonic scenario to look for alternate hadronic scenarios. Marked as a potential PeVatron candidate, this region has been studied extensively by H.E.S.S. due to its emission extending up to several tens of TeV. In this work, we use 2398 days of data from the High Altitude Water Cherenkov (HAWC) observatory to carry out a systematic source search of the HESS J1809-193 region. We were able to resolve emission detected as an extended component (modelled as a symmetric Gaussian with a 1 σ radius of 0.°21) with no clear cutoff at high energies and emitting photons up to 210 TeV. We model the multiwavelength observations for the region around HESS J1809-193 using a time-dependent leptonic model and a lepto-hadronic model. Our model indicates that both scenarios could explain the observed data within the region of HESS J1809-193.

The first TeV gamma-ray source with no lower energy counterparts, TeV J2032+4130, was discovered by HEGRA. It appears in the third HAWC catalog as 3HWC J2031+415 and it is a bright TeV gamma-ray source whose emission has previously been resolved as two sources: HAWC J2031+415 and HAWC J2030+409. While HAWC J2030+409 has since been associated with the Fermi Large Area Telescope Cygnus Cocoon, no such association for HAWC J2031+415 has yet been found. In this work, we investigate the spectrum and energy-dependent morphology of HAWC J2031+415. We associate HAWC J2031+415 with a gamma-ray binary system containing the pulsar PSR J2032+4127 and its companion MT91 213. We study HAWC data to observe their periastron in 2017. Additionally, we perform a combined multiwavelength analysis using radio, X-ray, and gamma-ray emission. We conclude that HAWC J2031+415 and, by extension, TeV J2032+4130 are most probably a pulsar wind nebula powered by PSR J2032+4127.