The simulation of extensive air showers and particle cascades in general is a cornerstone of modern astroparticle physics. For more than two decades, CORSIKA, currently in version 7, has been one of the most widely used tools for this purpose. However, its architecture reflects design constraints of an earlier computing era, as well as increasingly limiting extensibility, maintainability, and adaptability to modern experimental requirements. CORSIKA 8 is a complete redesign of the original CORSIKA code, implemented in modern C++ and based on contemporary software engineering principles. It introduces a modular and extensible simulation framework with explicit handling of units, flexible geometry, and environment descriptions. In this paper, we present the design philosophy and core architecture of CORSIKA 8, describe the implementation of electromagnetic and hadronic shower physics, and validate air shower simulations against CORSIKA 7. The results demonstrate good agreement at the similar to 5-10% level for key observables, confirming the physics fidelity of CORSIKA 8. The observed differences reflect the current level of intrinsic systematic uncertainties in state-of-the-art air-shower simulations. We also showcase new use cases that were beyond the capabilities of version 7, such as the simulation of cross-media showers and particle cascades in ice, including radio-signal propagation.

Massive brane fluctuations, called branons, behave as weakly interacting massive particles, which is one of the most favored class of candidates to fulfill the role of the dark matter (DM), an elusive kind of matter beyond the Standard Model. We present a multi-target search in dwarf spheroidal galaxies for branon DM annihilation signatures with a total exposure of 354 hours with the ground-based gamma-ray telescope system MAGIC. This search led to the most constraining limits on branon DM in the sub-TeV and multi-TeV DM mass range. Our most stringent limit on the thermally-averaged annihilation cross-section (at 95% confidence level) corresponds to ⟨ σv ⟩ ≃ 1.9 × 10 -24 cm 3 s -1 at a branon mass of ∼ 1.5 TeV.

The BL Lacertae object VER J0521+211 underwent a notable flaring episode in February 2020. A short-term monitoring campaign, led by the MAGIC (Major Atmospheric Gamma Imaging Cherenkov) collaboration, covering a wide energy range from radio to very high-energy (VHE, 100 GeV < E < 100 TeV) gamma rays was organised to study its evolution. These observations resulted in a consistent detection of the source over six consecutive nights in the VHE gamma-ray domain. Combining these nightly observations with an extensive set of multi-wavelength data made modelling of the blazar's spectral energy distribution (SED) possible during the flare. This modelling was performed with a focus on two plausible emission mechanisms: (i) a leptonic two-zone synchrotron-self-Compton scenario, and (ii) a lepto-hadronic one-zone scenario. Both models effectively replicated the observed SED from radio to the VHE gamma-ray band. Furthermore, by introducing a set of evolving parameters, both models were successful in reproducing the evolution of the fluxes measured in different bands throughout the observing campaign. Notably, the lepto-hadronic model predicts enhanced photon and neutrino fluxes at ultra-high energies (E > 100 TeV). While the photon component, generated via decay of neutral pions, is not directly observable as it is subject to intense pair production (and therefore extinction) through interactions with the cosmic microwave background photons, neutrino detectors (e.g. IceCube) can probe the predicted neutrino component. Finally, the analysis of the gamma-ray spectra, observed by MAGIC and the Fermi-LAT telescopes, yielded a conservative 95% confidence upper limit of z <= 0.244 for the redshift of this blazar.

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.

CORSIKA 8 is a new framework for simulations of particle cascades in air and dense media implemented in modern C++17, based on past experience with existing codes, in particular CORSIKA 7. The flexible and modular structure of the project allows the development of independent modules that can produce a fully customizable particle shower simulation. The radio module in particular is designed to treat the electric field calculation and its propagation through complex media to each observer location in an autonomous and flexible way. It already allows for the simultaneous simulation of the radio emission calculated with two independent time-domain formalisms, the "Endpoint formalism" as previously implemented in CoREAS and the "ZHS" algorithm as ported from ZHAireS. The design acts as the baseline interface for current and future development for the simulation of radio emission from particle showers in standard and complex scenarios, such as cross-media showers penetrating from air into ice. In this work, we present the design and implementation of the radio module in CORSIKA 8, along with validation studies and a direct comparison of the radio emission from air showers simulated with CORSIKA 8, CORSIKA 7 and ZHAireS. We also present the impact of simulation details such as the step size of simulated particle tracks on radio-emission simulations and perform a direct comparison of the "Endpoints" and "ZHS" formalisms for the same underlying air showers. Finally, we present an in-depth comparison of CORSIKA 8 and CORSIKA 7 for optimum simulation settings and discuss the relevance of observed differences in light of reconstruction efforts for the energy and mass of cosmic rays.

Aims. Mrk 421 was in its most active state around early 2010, which led to the highest TeV gamma-ray flux ever recorded from any active galactic nuclei (AGN). We aim to characterize the multiwavelength behavior during this exceptional year for Mrk 421, and evaluate whether it is consistent with the picture derived with data from other less exceptional years. Methods. We investigated the period from November 5, 2009, (MJD 55140) until July 3, 2010, (MJD 55380) with extensive coverage from very-high-energy (VHE; E>100 GeV) gamma rays to radio with MAGIC, VERITAS, Fermi-LAT, RXTE, Swift, GASP-WEBT, VLBA, and a variety of additional optical and radio telescopes. We characterized the variability by deriving fractional variabilities as well as power spectral densities (PSDs). In addition, we investigated images of the jet taken with VLBA and the correlation behavior among different energy bands. Results. Mrk 421 was in widely different states of activity throughout the campaign, ranging from a low-emission state to its highest VHE flux ever recorded. We find the strongest variability in X-rays and VHE gamma rays, and PSDs compatible with power-law functions with indices around 1.5. We observe strong correlations between X-rays and VHE gamma rays at zero time lag with varying characteristics depending on the exact energy band. We also report a marginally significant (similar to 3 sigma) positive correlation between high-energy (HE; E>100 MeV) gamma rays and the ultraviolet band. We detected marginally significant (similar to 3 sigma) correlations between the HE and VHE gamma rays, and between HE gamma rays and the X-ray, that disappear when the large flare in February 2010 is excluded from the correlation study, hence indicating the exceptionality of this flaring event in comparison with the rest of the campaign. The 2010 violent activity of Mrk 421 also yielded the first ejection of features in the VLBA images of the jet of Mrk 421. Yet the large uncertainties in the ejection times of these unprecedented radio features prevent us from firmly associating them to the specific flares recorded during the 2010 campaign. We also show that the collected multi-instrument data are consistent with a scenario where the emission is dominated by two regions, a compact and extended zone, which could be considered as a simplified implementation of an energy-stratified jet as suggested by recent IXPE observations.

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.

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 current generation of Imaging Atmospheric Cherenkov Telescopes (IACTs), comprised of major installations such as the MAGIC, H.E.S.S. and VERITAS telescopes, is sometimes called the 3rd generation of such instruments. These telescopes use multipixel cameras composed of hundreds up to thousands photomultiplier tubes (PMTs). The supremacy of PMTs is currently being challenged by photon sensors, rapidly spreading in popularity: the silicon photomultipliers (SiPMs), that are becoming a valid alternative thanks to their high PDE, low operating voltage and flexibility in installation. In this report, we investigate the performance of an existing 3rd generation IACT array (taking as a working example MAGIC), in which PMTs would be replaced with SiPMs, applying generalized simulations, not tuned for a specific hardware solution. We find an energy-dependent improvement in sensitivity, reaching a factor of three at the current trigger threshold energy. Interestingly, we also find that the stronger sensitivity of SiPMs in the red part of the spectrum, a source of background for IACTs, does not affect this performance.

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.

The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA.

ABSTRACT By studying the variability of blazars across the electromagnetic spectrum, it is possible to resolve the underlying processes responsible for rapid flux increases, so-called flares. We report on an extremely bright X-ray flare in the high-peaked BL Lacertae object Markarian 421 (Mrk 421) that occurred simultaneously with enhanced γ-ray activity detected at very high energies by First G-APD Cherenkov Telescope on 2019 June 9. We triggered an observation with XMM–Newton, which observed the source quasi-continuously for 25 h. We find that the source was in the brightest state ever observed using XMM–Newton, reaching a flux of 2.8 × 10−9 $\mathrm{erg\, cm^{-2}\, s^{-1}}$ over an energy range of 0.3–10 keV. We perform a spectral and timing analysis to reveal the mechanisms of particle acceleration and to search for the shortest source-intrinsic time-scales. Mrk 421 exhibits the typical harder-when-brighter behaviour throughout the observation and shows a clock-wise hysteresis pattern, which indicates that the cooling dominates over the acceleration process. While the X-ray emission in different sub-bands is highly correlated, we can exclude large time lags as the computed z-transformed discrete correlation functions are consistent with a zero lag. We find rapid variability on time-scales of 1 ks for the 0.3–10 keV band and down to 300 s in the hard X-ray band (4–10 keV). Taking these time-scales into account, we discuss different models to explain the observed X-ray flare, and find that a plasmoid-dominated magnetic reconnection process is able to describe our observation best.

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.

With the motivation to improve experimental gains and precision, established astroparticle experiments are currently undergoing massive upgrades. In addition, several new experiments are being built or planned. With the resulting gain in observational quality, the amount and accuracy of simulated data required for the analysis is also rising. In order to meet the increasing requirements and complexity due to the experiments' growth and to provide a unified software ecosystem, it was decided to re-develop the de facto standard extensive air shower simulation CORSIKA completely in C++ based on the original Fortran code. Since one of the largest runtime consumers is the propagation of millions of optical Cherenkov and fluorescence photons, and many experiments are starting to use them for measurements, it was decided to develop hardware-accelerated code to speed up the simulation. Specific methods have been developed to propagate photons on deep learning acceleration hardware similar to classical GPUs to take additional advantage of the current and future growth of the deep learning sector. In particular, Nvidia accelerators were tested.

Recently, several gamma-ray bursts (GRBs) have been detected in the very-high-energy (VHE) gamma-ray energy range by ground-based gamma-ray experiments such as MAGIC, H.E.S.S., and LHASSO. For some GRBs, the VHE emission is consistent with synchrotron self-Compton (SSC) emission from high-energy electrons accelerated in the forward shock of the relativistic jet. However, more statistics are needed to further constrain the emission models. GRB 201216C is a long bright GRB detected in a broad energy range from radio to VHE. The redshift is estimated to be 1.1, making the GRB the most distant source detected in the VHE energy range. MAGIC started the observation 56 seconds after the GRB was triggered by the Swift-BAT telescope. We performed a detailed analysis and detected the signal with about 6 sigma in the first 20 minutes. MAGIC continued the observation for 2.2 hours on the same night and 4 hours on the next night. No signal was detected later than 40 minutes after the GRB trigger. We have performed modelling of the multi-wavelength emission using the MAGIC data. We analysed simultaneous optical data from Liverpool Telescope with MAGIC and included the results in the modelling. The sub-TeV emission is consistent with the single-zone SSC model in the forward shock. In this presentation, we show the final results of the MAGIC data analysis of GRB 201216C and discuss the emission mechanism of the multi-wavelength data.