Search for: All records

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

Abstract Optical emissions associated with Terrestrial Gamma ray Flashes (TGFs) have recently become important subjects in space‐based and ground‐based observations as they can help us understand how TGFs are produced during thunderstorms. In this paper, we present the first time‐resolved leader spectra of the optical component associated with a downward TGF. The TGF was observed by the Telescope Array Surface Detector (TASD) simultaneously with other lightning detectors, including a Lightning Mapping Array (LMA), an INTerFerometer (INTF), a Fast Antenna (FA), and a spectroscopic system. The spectroscopic system recorded leader spectra at 29,900 frames per second (33.44 s time resolution), covering a spectral range from 400 to 900 nm, with 2.1 nm per pixel. The recordings of the leader spectra began 11.7 ms before the kA return stroke and at a height of 2.37 km above the ground. These spectra reveal that optical emissions of singly ionized nitrogen and oxygen occur between 167 s before and 267 s after the TGF detection, while optical emissions of neutrals (H I, 656 nm; N I, 744 nm, and O I, 777 nm) occur right at the moment of the detection. The time‐dependent spectra reveal differences in the optical emissions of lightning leaders with and without downward TGFs. 

Abstract On 11 September 2021, two small thunderstorms developed over the Telescope Array Surface Detector (TASD) that produced an unprecedented number of six downward terrestrial gamma ray flashes (TGFs) within one‐hour timeframe. The TGFs occurred during the initial stage of negative cloud‐to‐ground flashes whose return strokes had increasingly large peak currents up to 223 kA, 147 GeV energy deposit in up to 25 1.2 km‐spaced surface detectors, and intermittent bursts of gamma‐rays with total durations up to 717 s. The analyses are based on observations recorded by the TASD network, complemented by data from a 3D lightning mapping array, broadband VHF interferometer, fast electric field change sensor, high‐speed video camera, and the National Lightning Detection Network. The TGFs of the final two flashes had gamma fluences of and 8, logarithmically bridging the gap between previous TASD and satellite‐based detections. The observations further emphasize the similarity between upward and downward TGF varieties, suggesting a common mechanism for their production. 

ABSTRACT The dwarf spheroidal galaxies (dSphs) orbiting the Milky Way are widely regarded as systems supported by velocity dispersion against self-gravity, and as prime targets for the search for indirect dark matter (DM) signatures in the GeV-to-TeV $$\gamma$$-ray range owing to their lack of astrophysical $$\gamma$$-ray background. We present forecasts of the sensitivity of the forthcoming Cherenkov Telescope Array Observatory (CTAO) to annihilating or decaying DM signals in these targets. An original selection of candidates is performed from the current catalogue of known objects, including both classical and ultrafaint dSphs. For each, the expected DM content is derived using the most comprehensive photometric and spectroscopic data available, within a consistent framework of analysis. This approach enables the derivation of novel astrophysical factor profiles for indirect DM searches, which are compared with results from the literature. From an initial sample of 64 dSphs, eight promising targets are identified – Draco I, Coma Berenices, Ursa Major II, Ursa Minor, and Willman 1 in the North, Reticulum II, Sculptor, and Sagittarius II in the South – for which different DM density models yield consistent expectations, leading to robust predictions. CTAO is expected to provide the strongest limits above $$\sim$$10 TeV, reaching velocity-averaged annihilation cross sections of $$\sim 5\times 10^{-25}$$ cm$^3$ s$$^{-1}$$ and decay lifetimes up to $$\sim 10^{26}$$ s for combined limits. The dominant uncertainties arise from the imprecise determination of the DM content, particularly for ultrafaint dSphs. Observation strategies are proposed that optimize either deep exposures of the best candidates or diversified target selections. 

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

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 (∼3σ) positive correlation between high-energy (HE;E > 100 MeV) gamma rays and the ultraviolet band. We detected marginally significant (∼3σ) 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 recentIXPEobservations. 

We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array (TA) experiment with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion Letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extragalactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields. The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E>100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as 2 ×10−5 Mpc−3, that is the conservative lower limit for the source number density. 

We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ∼ 10 EeV that becomes lighter up to E ∼ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. 

ABSTRACT 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\gt $$ 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\pm 0.44_{\rm stat}\pm 0.15_{\rm sys}$$ (MAGIC) and $$3.39\pm 0.58_{\rm stat}\pm 0.64_{\rm 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 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 theFermi-LAT telescopes, yielded a conservative 95% confidence upper limit ofz ≤ 0.244 for the redshift of this blazar.