We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the HAWC and IceCube observatories in the Northern and Southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and holds a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map we determine the horizontal dipole components of the anisotropy δ 0h = 9.16 × 10 −4 and δ 6h = 7.25 × 10 −4 (±0.04 × 10 −4 ). In addition, we infer the direction (229.2 ± 3.5 • RA , 11.4 ± 3.0 • Dec.) of the interstellar magnetic field from the boundary between large scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large scale anisotropy to be δ N ∼ −3.97 +1.0 −2.0 × 10 −4 .

The results of 8-year R-band photopolarimetric data of blazar Mrk 421 collected from February 2008 to May 2016 are presented, along with extensive multiwavelength observations covering from radio to TeV γ-rays around the flares observed in May 2008, March 2010, and April 2013. The most important results are found in 2013 when the source displayed in the R-band a very high brightness state of 11.29 ± 0.03 mag (93.60 ± 1.53 mJy) on April 10th and a polarization degree of (11.00 ± 0.44)% on May 13th. The analysis of the optical data shows that the polarization variability is due to the superposition of two polarized components that might be produced in two distinct emitting regions. An intranight photopolarimetric variability study carried out over 7 nights after the 2013 April maximum found flux and polarization variations on the nights of April 14, 15, 16 and 19. In addition, the flux shows a minimum variability timescale of ∆ t =2.34±0.12 hours, and that the polarization degree presented variations of ∼ (1 -2) % in a timescale of ∆ t ∼ minutes. Also, a detailed analysis of the intranight data shows a coherence length of the large-scale magnetic field of l B 0.3 pc which is the same order of magnitude to the distance traveled by the relativistic shocks. This result suggests that there is a connection between the intranight polarimetric variations and spatial changes of the magnetic field. Analysis of the complete R-band data along with the historical optical light curve found for this object shows that Mrk 421 varies with a period of 16.26 ± 1.78 years.

The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy γ-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view, and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both data sets, the point-spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the γ-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. data set. Remaining differences in the overall γ-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the γ-ray sky between WCD and IACT techniques.

Gamma-ray bursts (GRBs) have been considered as potential very high energy photon emitters due to the large amount of energy released as well as the strong magnetic fields involved in their jets. However, the detection of teraelectronvolt photons is not expected from bursts beyond a redshift of z ≳ 0.1, due to their attenuation with the extragalactic background light (EBL). For these reasons, the recent observation of photons with energies of 18 and 251 TeV from GRB 221009A (z = 0.151) last 2022 October 9 has challenged what we know about the teraelectronvolt-emission mechanisms and the extragalactic background. In order to explain the teraelectronvolt observations, recent works exploring candidates of dark matter have started to appear. In this paper, we discuss the required conditions and limitations within the most plausible scenario, synchrotron self-Compton radiation in the GRB afterglow, to interpret the one 18 TeV photon observation besides the EBL. To avoid the Klein–Nishina effect, we find an improbable value of the microphysical magnetic parameter below 10−6 for a circumburst medium value >1 cm−3 (expected in the collapsar scenario). Therefore, we explore possible scenarios in terms of axion-like particles (ALPs) and dark photon mechanisms to interpret this highly energetic photon and we discuss the implications in the GRB energetics. We find that the ALPs and dark photon scenarios can explain the 18 teraelectronvolt photon but not the 251 teraelectronvolt photon.

Gamma-ray bursts (GRBs) are fascinating extragalactic objects. They represent a fantastic opportunity to investigate unique properties not exhibited in other sources. Multiwavelength afterglow observations from some short- and long-duration GRBs reveal an atypical long-lasting emission that evolves differently from the canonical afterglow light curves favoring the off-axis emission. We present an analytical synchrotron afterglow scenario and the hydrodynamical evolution of an off-axis top-hat jet decelerated in a stratified surrounding environment. The analytical synchrotron afterglow model is shown during the coasting, deceleration (off- and on-axis emission), and post–jet break decay phases, and the hydrodynamical evolution is computed by numerical simulations showing the time evolution of the Doppler factor, the half-opening angle, the bulk Lorentz factor, and the deceleration radius. We show that numerical simulations are in good agreement with those derived with our analytical approach. We apply the current synchrotron model and successfully describe the delayed nonthermal emission observed in a sample of long and short GRBs with evidence of off-axis emission. Furthermore, we provide constraints on the possible afterglow emission by requiring the multiwavelength upper limits derived for the closest Swift-detected GRBs and promising gravitational-wave events.