Gamma-ray Observations Research Articles

Axion-Photon Conversion of LHAASO Multi-TeV and PeV Photons

Abstract The Large High Altitude Air Shower Observatory (LHAASO) has reported the detection of a large number of multi-TeV-scale photon events including also several PeV-scale gamma-ray-photon events with energy as high as 1.4 PeV. The possibility that some of these events may have extragalactic origins is not yet excluded. Here we propose a mechanism for the traveling of very-high-energy (VHE) and ultra-high-energy (UHE) photons based upon the axion-photon conversion scenario, which allows extragalactic above-threshold photons to be detected by observers on the Earth. We show that the axion-photon conversation can serve as an alternative mechanism, besides the threshold anomaly due to Lorentz invariance violation, for the very-high-energy features of the newly observed gamma ray burst GRB 221009A.

Axion-Photon Conversion of LHAASO Multi-TeV and PeV Photons

The Large High Altitude Air Shower Observatory (LHAASO) has reported the detection of a large number of multi-TeV-scale photon events also including several PeV-scale gamma-ray-photon events with energy as high as 1.4 PeV. The possibility that some of these events may have extragalactic origins is not yet excluded. Here we propose a mechanism for the traveling of very-high-energy and ultra-high-energy photons based upon the axion-photon conversion scenario, which allows extragalactic above-threshold photons to be detected by observers on the Earth. We show that the axion-photon conversation can serve as an alternative mechanism, besides the threshold anomaly due to Lorentz invariance violation, for the very-high-energy features of the newly observed gamma ray burst GRB 221009A.

Evaluation of the potential of a gamma-ray observatory to detect astrophysical neutrinos through inclined showers

We assess the capabilities of a ground-based gamma-ray observatory to detect astrophysical neutrinos with energies in the $100\,{\rm TeV}$ to $100\,{\rm PeV}$ range. The identification of these events would be done through the measurement of very inclined extensive air showers induced by downward-going and upward-going neutrinos. The discrimination of neutrino-induced showers in the overwhelming cosmic-ray background is achieved by analysing the balance of the total electromagnetic and muonic signals of the shower at the ground. We demonstrate that a ${\rm km^2}$-scale wide field-of-view ground-based gamma-ray observatory could detect a couple of Very-High to Ultra-High energy (VHE-UHE) neutrino events per year with a reasonable pointing accuracy, making it an interesting facility for multi-messenger studies with both photons and neutrinos.

Validation of standardized data formats and tools for ground-level particle-based gamma-ray observatories

Context. Ground-based γ-ray astronomy is still a rather young field of research, with strong historical connections to particle physics. This is why most observations are conducted by experiments with proprietary data and analysis software, as is usual in the particle physics field. However, in recent years, this paradigm has been slowly shifting toward the development and use of open-source data formats and tools, driven by upcoming observatories such as the Cherenkov Telescope Array (CTA). In this context, a community-driven, shared data format (the gamma-astro-data-format, or GADF) and analysis tools such as Gammapy and ctools have been developed. So far, these efforts have been led by the Imaging Atmospheric Cherenkov Telescope community, leaving out other types of ground-based γ-ray instruments. Aims. We aim to show that the data from ground particle arrays, such as the High-Altitude Water Cherenkov (HAWC) observatory, are also compatible with the GADF and can thus be fully analyzed using the related tools, in this case, Gammapy. Methods. We reproduced several published HAWC results using Gammapy and data products compliant with GADF standard. We also illustrate the capabilities of the shared format and tools by producing a joint fit of the Crab spectrum including data from six different γ-ray experiments. Results. We find excellent agreement with the reference results, a powerful confirmation of both the published results and the tools involved. Conclusions. The data from particle detector arrays such as the HAWC observatory can be adapted to the GADF and thus analyzed with Gammapy. A common data format and shared analysis tools allow multi-instrument joint analysis and effective data sharing. To emphasize this, a sample of Crab nebula event lists is made public with this paper. Because of the complementary nature of pointing and wide-field instruments, this synergy will be distinctly beneficial for the joint scientific exploitation of future observatories such as the Southern Wide-field Gamma-ray Observatory and CTA.

Observe Gamma-Rays and Neutrinos Associated with Ultra-High Energy Cosmic Rays

IceCube measures a diffuse neutrino flux comparable to the Waxman-Bahcall bound, which suggests the possibility that the ultra-high energy cosmic rays (UHECRs) have a common origin with diffuse high energy neutrinos. We propose high energy gamma-ray and/or neutrino observations toward the arrival directions of UHECRs to search for the sources and test this possibility. We calculate the detection probability of gamma-ray/neutrino sources, and find that the average probability per UHECR of >10 EeV is ∼10% if the sensitivity of the gamma-ray or neutrino telescope is ∼10−12 erg cm−2 s−1 and the source number density is ∼10−5 Mpc−3. Future gamma-ray and neutrino observations toward UHECRs, e.g., by LHAASO-WCDA, CTA, IceCube/Gen2, are encouraged to constrain the density of UHECR sources or even identify the sources of UHECRs.

First constraints on axionlike particles from Galactic sub-PeV gamma rays

Experimental refinements and technical innovations in the field of extensive air shower telescopes have enabled measurements of Galactic cosmic-ray interactions in the sub-PeV range, providing new avenues for the search for new physics and dark matter. For the first time, we exploit sub-PeV (1 TeV -- 1 PeV) observations of Galactic diffuse gamma rays by HAWC and Tibet AS$\gamma$ to search for an axion-like-particle (ALP) induced gamma-ray signal directly linked to the origin of the IceCube extragalactic high-energy neutrino flux. Indeed, the production of high-energy neutrinos in extragalactic sources implies the concomitant production of gamma rays at comparable energies. Within the magnetic field of the neutrino emitting sources, gamma rays may efficiently convert into ALPs, escape their host galaxy un-attenuated, propagate through intergalactic space, and reconvert into gamma rays in the magnetic field of the Milky Way. Such a scenario creates an all-sky diffuse high-energy gamma-ray signal in the sub-PeV range. Accounting for the guaranteed Galactic astrophysical gamma-ray contributions from cosmic-ray interactions with gas and radiation and from sub-threshold sources, we set competitive upper limits on the photon-ALP coupling constant $g_{a\gamma\gamma}$. We find $g_{a\gamma\gamma} < 2.1\times10^{-11}$ GeV$^{-1}$ for ALP masses $m_a \leq 2\times10^{-7}$ eV at a 95\% confidence level. Our results are comparable to previous limits on ALPs derived from the TeV gamma-ray domain and progressively close the mass gap towards ADMX limits. The code and data to reproduce the results of this study are available on GitHub \url{https://github.com/ceckner/subPeVALPs}.

Probing particle acceleration at trans-relativistic shocks with off-axis gamma-ray burst afterglows

ABSTRACT Particle acceleration is expected to be different between relativistic and non-relativistic collisionless shocks. We show that electromagnetic counterparts to gravitational waves (GWs), gamma-ray burst (GRB) afterglows, are ideal targets for observing trans-relativistic evolution of accelerated electron distribution because the GWs spot nearby GRBs with off-axis jets, otherwise missed in gamma-ray observations. We find that the relativistic spectral slope begins to change steeply near the peak time of the light curve and approaches the non-relativistic limit in about 10 times the peak time. The trans-relativistic evolution of the afterglow synchrotron spectrum is consistent with GRB 170817A observations within errors, and will be measurable in similar but more distant events at a GW horizon ∼200 Mpc in a denser environment. We roughly estimate that such events represent a fraction of 10–50 per cent of the GRB 170817A-like off-axis short GRBs. We also find that the spectral evolution does not depend on the jet structure if their light curves are similar to each other.

The Mercedes water Cherenkov detector

The concept of a small, single-layer water Cherenkov detector, with three photomultiplier tubes (PMTs), placed at its bottom in a 120^{circ } star configuration (Mercedes Water Cherenkov Detector) is presented. The PMTs are placed near the lateral walls of the stations with an adjustable inclination and may be installed inside or outside the water volume. To illustrate the technical viability of this concept and obtain a first-order estimation of its cost, an engineering design was elaborated. The sensitivity of these stations to low energy Extensive Air Shower (EAS) electrons, photons and muons is discussed, both in compact and sparse array configurations. It is shown that the analysis of the intensity and time patterns of the PMT signals, using machine learning techniques, enables the tagging of muons, achieving an excellent gamma/hadron discrimination for TeV showers. This concept minimises the station production and maintenance costs, allowing for a highly flexible and fast installation. Mercedes Water Cherenkov Detectors (WCDs) are thus well-suited for use in high-altitude large gamma-ray observatories covering an extended energy range from the low energies, closing the gap between satellite and ground-based measurements, to very high energy regions, beyond the PeV scale.

Gamma/hadron discrimination at high energies through the azimuthal fluctuations of air shower particle distributions at the ground

Wide field-of-view gamma-ray observatories must fight the overwhelming cosmic ray background to identify very-high-energy astrophysical gamma-ray events.This work introduces a novel gamma/hadron discriminating variable, LCm, which quantifies the azimuthal non-uniformity of the particle distributions at the ground. This non-uniformity, due to the presence of hadronic sub-showers, is higher in proton-induced showers than in gamma showers. The discrimination power of this new variable is then discussed, as a function of the air shower array fill factor, in the energy range 10TeV to 1PeV, and compared to the classical gamma/hadron discriminator based on the measurement of the number of muons at the ground. The results obtained are extremely encouraging, paving the way for the use of the proposed quantity in present and future large ground-array gamma-ray observatories.

Classification of gamma-ray targets for velocity-dependent and subhalo-boosted dark-matter annihilation

Gamma-ray observations have long been used to constrain the properties of dark matter (DM), with a strong focus on weakly interacting massive particles annihilating through velocity-independent processes. However, in the absence of clear-cut observational evidence for the simplest candidates, the interest of the community in more complex DM scenarios involving a velocity-dependent cross-section has been growing steadily over the past few years. We present the first systematic study of velocity-dependent DM annihilation (in particular p-wave annihilation and Sommerfeld enhancement) in a variety of astrophysical objects, not only including the well-studied Milky Way dwarf satellite galaxies, but nearby dwarf irregular galaxies and local galaxy clusters as well. Particular attention is given to the interplay between velocity dependence and DM halo substructure. Uncertainties related to halo mass, phase-space and substructure modelling are also discussed in this velocity-dependent context. We show that, for s-wave annihilation, extremely large subhalo boost factors are to be expected, up to 1011 in clusters and up to 106–107 in dwarf galaxies where subhalos are usually assumed not to play an important role. Boost factors for p-wave annihilation are smaller but can still reach 103 in clusters. The angular extension of the DM signal is also significantly impacted, with e.g. the cluster typical emission radius increasing by a factor of order 10 in the s-wave case. We also compute the signal contrast of the objects in our sample with respect to annihilation happening in the Milky Way halo. Overall, we find that the hierarchy between the brightest considered targets depends on the specific details of the assumed particle-physics model.

An updated view and perspectives on high-energy gamma-ray emissionfrom SGR J1935+2154 and its environment

SGR J1935+2154 was discovered in 2016 and is currently one of the most burst-active Soft Gamma-ray Repeaters (SGR), having emitted many X-ray bursts in recent years. In one of our previous articles, we investigated the contribution to high-energy and very high-energy gamma-ray emission (VHE, E > 100 GeV) due to cosmic-ray acceleration of SNR G57.2+0.8 hosting SGR J1935+2154 using the GALPROP propagation code. However, follow-up observations of SGR 1935+2154 were made for 2 hours on April 28, 2020, using the High Energy Stereoscopic System (H.E.S.S.). The observations coincide with X-ray bursts detected by INTEGRAL and Fermi/Gamma-ray Burst Monitor (GBM). These are the first high-energy gamma-ray observations of an SGR in a flaring state, and upper limits on sustained and transient emission have been derived. Now that new H.E.S.S. observations have been made, it is interesting to update our model with respect to these new upper limits. We extend our previous results to a more general situation using the new version of GALPROP. We obtain a hadronic model that confirms the results discussed by H.E.S.S. . This leads to an optimistic prospect that cosmic ray gamma rays from SGR J1935+2154 can contribute to the overall gamma energy density distribution and in particular to the diffusion gamma rays from the Galactic center.

Geminga contribution to the cosmic-ray positron excess according to the gamma-ray observations

We attempt to interpret the cosmic-ray positron excess by injection from the nearby pulsar Geminga, assuming a two-zone diffusion scenario and an injection spectrum with a low energy cutoff. Since the high energy positrons and electrons from Geminga can induce γ rays via inverse Compton scattering, we take into account the extended γ-ray observations around Geminga from HAWC for ∼10 TeV and from Fermi-LAT for GeV. According to the extended γ-ray observation claimed by an analysis of Fermi-LAT data, we find that Geminga could explain the positron excess for a 30% energy conversion efficiency into positrons and electrons. However, based on the constraint on the extended γ rays given by another Fermi-LAT analysis, positrons from Geminga would be insufficient to account for the positron excess. Further robust analysis of Fermi-LAT data for the extended γ rays would be crucial to determine whether Geminga can explain the positron excess or not.

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

ABSTRACT MAXI J1820+070 is a low-mass X-ray binary with a black hole (BH) as a compact object. This binary underwent an exceptionally bright X-ray outburst from 2018 March to October, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 h of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy (HE) gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ∼500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential HE and very-HE gamma-ray emitting region should be located at a distance from the BH ranging between 1011 and 1013 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA.

Constraining Einstein-dilaton models using joint gravitational-wave and electromagnetic observations

With the advent of gravitational-wave astronomy it has now been possible to constrain modified theories of gravity that were invoked to explain the dark energy. In a class of dilaton models, distances to cosmic sources inferred from electromagnetic and gravitational wave observations would differ due to the presence of a friction term. In such theories, the ratio of the Newton's constant to the fine structure constant varies with time. In this paper we explore the degree to which it will be possible to test such models. If collocated sources (e.g. supernovae and binary neutron star mergers), but not necessarily multimessengers, can be identified by electromagnetic telescopes and gravitational-wave detectors one can probe if light and gravitational radiation are subject to the same laws of propagation over cosmological distances. This helps in constraining the variation of Newton's constant relative to fine-structure constant. The next generation of gravitational wave detectors, such as the Cosmic Explorer and Einstein Telescope, in tandem with the Vera Rubin Observatory and gamma ray observatories such as the Fermi Space Observatory will be able to detect or constrain such variations at the level of a few parts in 100. We apply this method to GW170817 with distances inferred by the LIGO and Virgo detectors and the observed Kilonova.

Constraints on the Very High Energy Gamma-Ray Emission from Short GRBs with HAWC

Many gamma-ray bursts (GRBs) have been observed from radio wavelengths, and a few at very high energies (VHEs, >100 GeV). The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is well suited to study transient phenomena at VHEs owing to its large field of view and duty cycle. These features allow for searches of VHE emission and can probe different model assumptions of duration and spectra. In this paper, we use data collected by HAWC between 2014 December and 2020 May to search for emission in the energy range from 80 to 800 GeV coming from a sample of 47 short GRBs that triggered the Fermi, Swift, and Konus satellites during this period. This analysis is optimized to search for delayed and extended VHE emission within the first 20 s of each burst. We find no evidence of VHE emission, either simultaneous or delayed, with respect to the prompt emission. Upper limits (90% confidence level) derived on the GRB fluence are used to constrain the synchrotron self-Compton forward-shock model. Constraints for the interstellar density as low as 10−2 cm−3 are obtained when assuming z = 0.3 for bursts with the highest keV fluences such as GRB 170206A and GRB 181222841. Such a low density makes observing VHE emission mainly from the fast-cooling regime challenging.

CRPropa 3.2 — an advanced framework for high-energy particle propagation in extragalactic and galactic spaces

The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, largely due to the inflow of data collected by large-scale cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe.This version represents a major leap forward, significantly expanding the simulation framework and opening up the possibility for many more astrophysical applications. This includes, among others: efficient simulation of high-energy particles in diffusion-dominated domains, self-consistent and fast modelling of electromagnetic cascades with an extended set of channels for photon production, and studies of cosmic-ray diffusion tensors based on updated coherent and turbulent magnetic-field models. Furthermore, several technical updates and improvements are introduced with the new version, such as: enhanced interpolation, targeted emission of sources, and a new propagation algorithm (Boris push). The detailed description of all novel features is accompanied by a discussion and a selected number of example applications.

COMPTEL data analysis using GammaLib and ctools

More than 20 yr after the end of NASA’s Compton Gamma-Ray Observatory mission, the data collected by its Imaging Compton Telescope (COMPTEL) still provide the most comprehensive and deepest view of our Universe in MeV gamma rays. While most of the COMPTEL data are archived at NASA’s High Energy Astrophysics Science Archive Research Center (HEASARC), the absence of any publicly available software for their analysis means the data cannot benefit from the scientific advances made in the field of gamma-ray astronomy at higher energies. To make this unique treasure again accessible for science, we developed open source software that enables a comprehensive and modern analysis of the archived COMPTEL telescope data. Our software is based on a dedicated plug-in to the GammaLib library, a community-developed toolbox for the analysis of astronomical gamma-ray data. We implemented high-level scripts for building science analysis workflows in ctools, a community-developed gamma-ray astronomy science analysis software framework. We describe the implementation of our software and provide the underlying algorithms. Using data from the HEASARC archive, we demonstrate that our software reproduces derived data products that were obtained in the past using the proprietary COMPTEL software. We furthermore demonstrate that our software reproduces COMPTEL science results published in the literature. This brings the COMPTEL telescope data back into life, allowing them to benefit from recent advances in gamma-ray astronomy, and gives the community a means to unveil its still hidden treasures.

TeV Dark Matter Searches in the Extragalactic Gamma-ray Sky

High-energetic gamma rays from astrophysical targets constitute a unique probe for annihilation or decay of heavy particle dark matter (DM). After several decades, diverse null detections have resulted in strong constraints for DM particle masses up to the TeV scale. While the gamma-ray signature is expected to be universal from various targets, uncertainties of astrophysical origin strongly affect and weaken the limits. At the same time, spurious signals may originate from non-DM related processes. The many gamma-ray targets in the extragalactic sky being searched for DM play a crucial role to keep these uncertainties under control and to ultimately achieve an unambiguous DM detection. Lately, a large progress has been made in combined analyses of TeV DM candidates towards different targets by using data from various instruments and over a wide range of gamma-ray energies. These approaches not only resulted in an optimal exploitation of existing data and an improved sensitivity, but also helped to level out target- and instrument-related uncertainties. This review gathers all searches in the extragalactic sky performed so far with the space-borne Fermi-Large Area Telescope, the ground-based imaging atmospheric Cherenkov telescopes, and the High-Altitude Water Cherenkov Gamma-Ray Observatory (HAWC). We discuss the different target classes and provide a complete list of all analyses so far.

Evidence for PeV Proton Acceleration from Fermi-LAT Observations of SNR G106.3+2.7.

The existence of a "knee" at energy ∼1 PeV in the cosmic-ray spectrum suggests the presence of Galactic PeV proton accelerators called "PeVatrons." Supernova remnant (SNR) G106.3+2.7 is a prime candidate for one of these. The recent detection of very high energy (0.1-100TeV) gamma rays from G106.3+2.7 may be explained either by the decay of neutral pions or inverse Compton scattering by relativistic electrons. We report an analysis of 12years of Fermi-LAT gamma-ray data that shows that the GeV-TeV gamma-ray spectrum is much harder and requires a different total electron energy than the radio and x-ray spectra, suggesting it has a distinct, hadronic origin. The nondetection of gamma rays below 10GeV implies additional constraints on the relativistic electron spectrum. A hadronic interpretation of the observed gamma rays is strongly supported. This observation confirms the long-sought connection between Galactic PeVatrons and SNRs. Moreover, it suggests that G106.3+2.7 could be the brightest member of a new population of SNRs whose gamma-ray energy flux peaks at TeV energies. Such a population may contribute to the cosmic-ray knee and be revealed by future very high energy gamma-ray detectors.

Sensitivity of future gamma-ray telescopes to primordial black holes

The strongest existing constraints on primordial black holes with masses in the range of $m_{\rm BH} \sim 10^{15}-10^{17} \, {\rm g}$ have been derived from measurements of the local cosmic-ray electron-positron flux by Voyager 1, and MeV-scale gamma-ray observations of the Inner Galaxy by COMPTEL and INTEGRAL. In this paper, we evaluate the sensitivity of future MeV-scale gamma-ray telescopes such as e-ASTROGAM or AMEGO to Hawking radiation. We show that such an instrument would be able to provide the strongest constraints on black holes in the mass range of $m_{\rm BH} \sim (0.6-20) \times 10^{16} \, {\rm g}$, typically exceeding current constraints by approximately two orders of magnitude. In scenarios in which the observed 511 keV excess is the result of Hawking radiation, we find that e-ASTROGAM or AMEGO would not only be able to detect the Hawking radiation from the Inner Galaxy, but could precisely measure the abundance and mass distribution of the black holes responsible for this signal.