A pilot search for MeV gamma-ray emission from five galaxy clusters using archival COMPTEL data
Abstract
We search for MeV gamma-ray emission between 0.75-30 MeV from five galaxy clusters, viz. Coma, VIRGO, SPT-CL J2012-5649, Bullet, and El Gordo, using archival data from the COMPTEL telescope. For this purpose we use three search templates: point source, radial disk and radial Gaussian. We do not detect any signals from Coma, SPT-CL J2012-5649, Bullet and El Gordo clusters with the 95% c.l. photon energy flux limit . For VIRGO, we detect a non-zero signal between 0.75 to 1.50 MeV having marginal significance of about , with the observed energy flux equal to . However, we do not confirm the previously reported evidence in literature for a gamma-ray line from Coma and VIRGO clusters between 5-7 MeV.
I Introduction
Galaxy clusters are the largest gravitationally bound and virialized structures in the Universe and act as wonderful laboratories for Cosmology, galaxy evolution Kravtsov and Borgani (2012); Allen et al. (2011); Vikhlinin et al. (2014) and fundamental Physics Böhringer and Chon (2016); Desai (2018); Bora and Desai (2021a, b); Mendonça et al. (2021); Bora et al. (2022). They have been observed over a broad range of the electromagnetic spectrum from radio waves Feretti et al. (2012) to hard X-rays with energies up to 30 keV Wik et al. (2014); Rojas Bolivar et al. (2023). In the last two decades a large number of new galaxy clusters have been discovered through dedicated surveys at optical, X-ray, and millimeter wavelengths. At higher energies, the detection of gamma-rays from galaxy clusters (at GeV energies) is still a hotly debated issue, with only the Coma cluster showing evidence for gamma-ray emission from the intra-cluster medium Keshet and Reiss (2018); Xi et al. (2018); Baghmanyan et al. (2022). Gamma-rays could be produced in galaxy clusters due to hadronic interactions or inverse Compton scattering processes in the intra-cluster medium Pinzke and Pfrommer (2010). Another possible source of gamma-rays in clusters could be from dark matter annihilations if the dark matter is a WIMP Colafrancesco et al. (2006). All possible mechanisms of gamma-ray production in clusters have been recently reviewed in Manna and Desai (2024).
Most recently, we carried out a systematic search for gamma-ray emission between 1-300 GeV using 15 years of Fermi-LAT data. Manna and Desai (2024), from a sample of 300 galaxy clusters selected from the South Pole Telescope 2500 square degree cluster survey Bleem et al. (2015); Bocquet et al. (2019), whose redshifts were estimated from optical/NIR follow-ups Desai et al. (2012); Song et al. (2012). The clusters in the SPT-SZ survey have been detected using the Sunyaev-Zeldovich effect Sunyaev and Zeldovich (1972) and hence represent a mass limited sample. From this catalog, we were able to detect gamma-ray emission from one galaxy cluster with significance , viz. SPT-CLJ 2012-5649 with six other clusters having significance between Manna and Desai (2024). However we could not discern if this emission is from the intra-cluster medium or because of radio galaxies in the cluster.
However, there is a paucity of studies on clusters in the soft gamma-ray energy band between hard X-rays and gamma-rays at MeV energies . One reason is that there has only been one space-based telescope that has imaged the sky in this energy range from 0.75-30 MeV, which is the Imaging Compton Telescope (COMPTEL). COMPTEL was one of the three telescopes onboard the Compton Gamma Ray Observatory (CGRO), which operated from 1991-2000 Schoenfelder et al. (1993). No other space-based gamma-ray telescope launched within the last two decades has covered the above energy range. The Fermi-LAT and AGILE telescopes are sensitive at higher energies ( 30 MeV), whereas INTEGRAL is sensitive up to 10 MeV. Therefore, there has only been one paper in literature which has searched for gamma-ray emission from galaxy clusters at energies between between 0.75-30 MeV using COMPTEL Iyudin et al. (2004) (I04, hereafter). This work looked for MeV emission from galaxy clusters near the North Galactic pole. I04 found evidence for a gamma-ray line between 5.6-7.6 MeV with statistical significance of around from Coma cluster as well as from two elliptical galaxies M87 and M49 in the VIRGO cluster, which they argued is due to the de-excitation gamma-ray line emission from , which could be produced from spallation interactions of cosmic rays with ions in the intra-cluster medium. The observed -ray line flux was equal to photons//sec and photons//sec, from VIRGO and Coma, respectively. We also note that M87 is a FR-I radio galaxy with a misaligned jet and has been detected at very high energies ( GeV) by Fermi-LAT, H.E.S.S., MAGIC, and VERITAS Molero et al. (2023).
Although the COMPTEL data is still archived at NASA’s HEASARC, the legacy software used to process the COMPTEL data was decommissioned shortly after the end of CGRO mission and was not easily accessible to the wider scientific community. Therefore, there were only a handful of results from the COMPTEL mission in the last twenty odd years (Strong and Collmar (2019) and references therein). Recently, a new plugin as well as custom Python scripts have been added to the widely used GammaLib and ctools libraries Knödlseder et al. (2016) to analyze the COMPTEL data, facilitating a seamless scientific analysis Knödlseder et al. (2022). Given that the energy range probed by COMPTEL is an uncharted territory with no approved experiment in the near future probing this energy range, it behooves us to analyze this legacy data to look for gamma-ray emissions from galaxy clusters.
Given some of the intriguing hints for GeV gamma-ray emission from galaxy clusters Manna and Desai (2024), it would be interesting to look for detections at MeV energies from the same clusters. Therefore, as a pilot search, we analyze the archival COMPTEL data to look for detections at MeV energies . The targets we consider for our analysis include VIRGO, Coma, SPT-CLJ 2012-5649, Bullet (SPT-CL J0658-5556), and El Gordo (SPT-CL J0102-4915) clusters. As mentioned above, I04 had detected MeV emission from both Coma and the VIRGO clusters. Most recently, we detected emission from SPT-CL J2012-5649 (also known as Abell 3667) between 1-10 GeV using the Fermi-LAT telescope. Therefore, it would be interesting to check if the signal seen in SPT-CL J2012-5649 at GeV energies persists at lower energies. Both the Bullet and El Gordo clusters are massive merging clusters, which have proved to be a very good laboratories for testing CDM and modified gravity theories Clowe et al. (2006); Asencio et al. (2023). Therefore they would also be very good targets for gamma-ray searches.
II Data Analysis
The COMPTEL telescope, onboard the Compton Gamma Ray Observatory (CGRO) launched in 1991, was a pioneering instrument that revolutionized our understanding of the high-energy universe. It employed a unique imaging technique called Compton scattering, where gamma rays interact with electrons within the detector, changing their direction. By analyzing these scattered photons, it could reconstruct the location and energy of the original gamma-ray source. The telescope field of view was about one steradian. It was sensitive to gamma-rays between 0.75- 30 MeV, where the energy and angular resolution ranged between 5-8% and , respectively depending on the photon energy. During its 9.7 years of operation, COMPTEL did about 340 distinct pointings (with a duration of two weeks), where each pointing had a field of view radius of around Strong and Collmar (2019). More details on the observing specifications, detector calibration, and performance of COMPTEL can be found in Schoenfelder et al. (1993). A large number of astrophysical sources have been detected by COMPTEL such as pulsars, AGNs, X-ray binaries, gamma-ray bursts, solar flares, 26Al, supernova remnants, extragalactic diffuse gamma-ray background, etc Schoenfelder et al. (1993); Schönfelder et al. (2000).
We comprehensively analysed the COMPTEL data to investigate the potential presence of MeV gamma-ray emission from five galaxy clusters: the Coma Cluster, SPT-CL J2012-5649 Cluster, the Bullet Cluster,the El Gordo Cluster and the Virgo Cluster. We carried out a systematic search within a radius surrounding each target cluster utilizing the revamped COMPTEL analysis framework implemented in the ctools software package Knödlseder et al. (2022). We employed the comobsselect tool to extract the pertinent viewing periods for the clusters under investigation, applying precise filtering based on their coordinates and a specified radius. The selection process effectively removes unwanted background noise. Using the comobsbin tool, we implemented energy binning to optimize the signal-to-noise ratio. We used 16 logarithmically-spaced energy bins spanning the 0.75-30 MeV range. With this spacing, we could detect the gamma-ray lines which have been reported in I04. For SPT-CL J2012-5649, El Gordo, VIRGO and Bullet Cluster, we combined the data from two or more distinct viewing periods using the comobsadd tool with 80 bins each in the and directions, where and are the Compton scattering directions. Unlike the other clusters, for Coma cluster, we got data from only one viewing period in our specified radius (). The effective exposure for all our clusters can be found in Table 1.
Our search analysis was done using three different templates, viz. point source, radial disk and radial Gaussian templates for our analysis. We used the same radial model definitions used in Fermi-LAT analyses Wood et al. (2017) 111available at https://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/xml_model_defs.html. For the spectral analysis, we utilized the comobsmodel tool to fit source emission models to the COMPTEL data. We then employed the comlixfit tool to perform maximum likelihood fitting of the COMPTEL data utilizing the iterative SRCLIX algorithm. Furthermore, the instrumental background Model was modelled using the BGDLIXE background computation model. By accounting for the instrumental background components, BGDLIXE enables a more accurate and reliable characterization of the true MeV gamma-ray signal originating from the target galaxy clusters. We leveraged the comlixmap tool to construct Test Statistics (TS) maps (cf. Eq. 1) by applying the SRCLIX on every test source position. The test statistic provides a quantitative indication of the evidence against a null hypothesis. Between a model with both signal and background components () and a model with only the background (), the test statistic (TS) can be expressed as Mattox et al. (1996):
(1) |
where denotes the likelihood when both the source model () and the background model () are jointly fitted to the data. It reflects how well the combined model aligns with the observed data, while corresponds to the log-likelihood, when only the background model is considered. Crucially, under the assumption that the background model adequately captures the data, the TS follows a distribution with degrees of freedom, where denotes the number of free parameters within the source model component. The same statistics is also used in IceCube and Fermi-LAT analysis Pasumarti and Desai (2022); Manna and Desai (2024).
For each cluster we also constructed the Spectral Energy Distribution (SED) plots using the csspec tool and obtained upper limits (in case of null results) with the ctulimit tool. These SED plots, constructed within the energy range of 0.75-30 MeV using 16 logarithmic energy bins and the BINS spectrum generation method, visually represent the upper limits on the flux density of potential gamma-ray sources within the clusters. The upper limits for the differential photon flux were determined using a reference energy of 30 MeV, chosen to maximize the sensitivity to potential spectral variations within the MeV range. We used BINS method to run the csspec script, the spectral model is replaced by a binned function that fits all data. We calculate integral photon and energy flux upper limits over the full energy range of 0.75 MeV to 30 MeV, encompassing the analyzed spectrum to provide constraints for clusters with null detections. The energy flux signifies the cumulative photon energy deposited within a specified energy range.
III Results
We now report the results from our searches from the five clusters between the energy range of 0.75-30 MeV, spanning the full COMPTEL observing period from 1991 to 2000. Table 1 summarizes the maximum TS value for each cluster for all the three templates along with the effective exposure. For Coma, SPT-CL J2012-5649, Bullet, and El Gordo, the TS values are less than 3.0. For VIRGO we see TS values greater than 3.0 for all the three input templates, with the largest value of 5.6 obtained for the radial disk template. Consequently, the corresponding SED plots only show upper limits for all the clusters except VIRGO. For illustration, we show the corresponding SED plots for a point source template for Coma, SPT-CL J2012-5649, Bullet, and El Gordo clusters in Figures 1, 2, 3, 4, respectively. The upper limits for the other templates are of similar values. For the VIRGO cluster, we show the SED plots for all the three templates. We obtain a non-zero flux in the first three energy bins between to MeV, having marginal significance between and null detections at higher energies. For the VIRGO cluster, we first used the fixed position of RA= and Dec=. We also redid the fit by keeping the position a free parameter and we get a best position consistent with our initial position and with the same SED as before. For the Coma as well as VIRGO cluster we also redid our analyses with different binnings including the one used in I04, but it does not change our results. Therefore, we cannot reproduce the results for line emission from the Coma cluster reported in I04. Although we detect a signal from the VIRGO cluster, it is seen at a lower energy than what was found in I04.
We summarize our results for all the clusters in Table 2. We report the 95% confidence level upper limits for the differential photon flux (at 30 MeV), the integral photon flux, and photon energy flux. For Coma, SPT-CL J2012-5649, Bullet, and El Gordo, we report the differential flux upper limit at 30 MeV, along with integral flux and energy flux upper limits. For Coma cluster, the integral photon flux limit we obtain of about is comparable to the reported flux in I04. All the remaining clusters with null results also show similar values for the differential and integral photon flux as well as the energy flux. For the VIRGO cluster, we report the aforementioned measured quantities along with error bars, which can be found in the last three rows of Table 2. For the point source template, the total energy flux equals erg/cm2/s. When we use the radial disk and radial Gaussian templates, we get total photon energy fluxes of erg/cm2/s and erg/cm2/s respectively. Once again, the photon flux is comparable to the line flux reported in I04, albeit measured at lower energies. We note that the AGN in M87 could be the cause of the observed MeV emission,s ince M87 has been detected at very high energies.
Cluster | Template | TS Values |
Exposure
(cm2/s) |
---|---|---|---|
Point | 0.02 | - | |
Coma | Radial Disk | 0.06 | |
Radial Gaussian | 0.01 | ||
Point | 1.44 | - | |
SPT-CL J2012-5649 | Radial Disk | 2.74 | |
Radial Gaussian | 1.9 | ||
Point | 0 | - | |
Bullet | Radial Disk | 0 | |
Radial Gaussian | 0 | ||
Point | 1.12 | - | |
El Gordo | Radial Disk | 2.40 | |
Radial Gaussian | 1.00 | ||
Point | 3.0 | - | |
VIRGO | Radial Disk | 5.6 | |
Radial Gaussian | 3.6 |
Cluster | Template |
Differential Flux at 30 MeV
(ph/cm2/s/MeV) |
Integral Flux
(ph/cm2/s) |
Energy Flux
(erg/cm2/s) |
---|---|---|---|---|
Point | ||||
Coma | Radial Disk | |||
Radial Gaussian | ||||
Point | ||||
SPT-CL J2012-5649 | Radial Disk | |||
Radial Gaussian | ||||
Point | ||||
Bullet | Radial Disk | |||
Radial Gaussian | ||||
Point | ||||
El Gordo | Radial Disk | |||
Radial Gaussian | ||||
Point | ||||
VIRGO | Radial Disk | |||
Radial Gaussian |
IV Conclusions
There has only been one previous search (reported in I04) for gamma-ray emission from galaxy clusters in the energy range between 0.75-30 MeV using the COMPTEL telescope. I04 had found a gamma-ray line between 5-7 MeV from two galaxy clusters namely, VIRGO and Coma. Although, there have been a plethora of sources detected at GeV energies from the Fermi-LAT satellite, no other systematic search for emission from galaxy clusters has been carried in the MeV energy region. The main reason is that there is no other telescope apart from COMPTEL, which has mapped the universe in MeV gamma-rays. Even the COMPTEL telescope was decommissioned in 2000 and the software used to analyze the data was available only to the core institutes involved in the COMPTEL collaboration.
In the last two years, the software used to analyze the archival COMPTEL data has been revamped and made publicly available, facilitating a seamless analysis of this data Knödlseder et al. (2022). We leverage on this to carry out a pilot search for MeV gamma-ray emission using COMPTEL from five galaxy clusters, namely Coma, VIRGO, Bullet, El Gordo and SPT-CL J2012-5649. We used three search templates for the same: point source, radial disk and radial Gaussian. Our results are summarized in Table 2. We report null results from Coma, Bullet, El Gordo, and SPT-CLJ 2012-5649. We calculate the 95% c.l. upper flux limits on the differential photon flux (at 30 MeV), integral photon flux, and energy flux. The energy flux limits are erg//sec for all the clusters with no detections (cf. Table 2). Therefore, we do not confirm the gamma-ray line seen for the Coma cluster in I04. For the VIRGO cluster, we detect a signal between 0.75 and 1.50 MeV with significance between and null detections at higher energies up to 30 MeV. The SED for the VIRGO cluster using all the three search templates can be found in Figs. 5, 6, 7. The observed energy flux is between , depending on the template used. The observed signal is at a lower energy than that observed in I04. Therefore, although we see a signal at low significance around 1 MeV, we do not corroborate the gamma-ray line observed for VIRGO in I04. We note that since M87 (in the center VIRGO cluster) contains a misaligned jet from which VHE emissions have been detected, this AGN in M87 could also give rise to the observed MeV emission.
Therefore, this pilot search for MeV gamma-ray emission with COMPTEL from a few selected galaxy clusters is a proof of principles application of the revamped COMPTEL software in the analysis of galaxy clusters and paves the way for extension to other clusters. We shall also extend this analysis to X-ray and SZ selected samples. We shall also discuss theoretical implications of our results, including constraints on inverse-Compton emission and possible limits on magnetic fields for each of the clusters in a future work. More detailed studies of galaxy clusters at MeV energies should soon be possible with future missions such as COSI Zoglauer et al. (2021) or GECCO Orlando et al. (2022), which are sensitive up to 10 MeV. Beyond 10 MeV, a whole bunch of experiments have been proposed but not yet scheduled for launch at the time of writing, such as GAMMA-400 Galper et al. (2014); Egorov et al. (2020), Advanced Compton Telescope Boggs (2006), Advanced Energetic Pair Telescope (AdEPT) Hunter et al. (2014), PANGU Wu et al. (2014, 2015) , GRAMS Aramaki et al. (2020a, b), MAST Dzhatdoev and Podlesnyi (2019), All-Sky-ASTROGAM Tatischeff et al. (2019). All of these missions should provide more insights into gamma-ray emission from galaxy clusters at energies greater than 10 MeV.
In the spirit of open science, we have made our analyses templates and data used for making the SED plots publicly available at https://github.com/siddhantmannaiith/Comptel-Analysis-
Acknowledgements.
SM extend his sincere gratitude to the Government of India, Ministry of Education (MOE) for their continuous support through the stipend, which has played a crucial role in the successful completion of our research. We are also very grateful to Jürgen Knödlseder for thoroughly explaining the nuances of the COMPTEL analysis software and helping us with the analysis. We also acknowledge the anonymous referee for several useful comments on the manuscript.References
- Kravtsov and Borgani (2012) A. V. Kravtsov and S. Borgani, Ann. Rev. Astron. Astrophys. 50, 353 (2012), eprint 1205.5556.
- Allen et al. (2011) S. W. Allen, A. E. Evrard, and A. B. Mantz, Ann. Rev. Astron. Astrophys. 49, 409 (2011), eprint 1103.4829.
- Vikhlinin et al. (2014) A. A. Vikhlinin, A. V. Kravtsov, M. L. Markevich, R. A. Sunyaev, and E. M. Churazov, Physics Uspekhi 57, 317-341 (2014).
- Böhringer and Chon (2016) H. Böhringer and G. Chon, Modern Physics Letters A 31, 1640008 (2016), eprint 1610.02855.
- Desai (2018) S. Desai, Physics Letters B 778, 325 (2018), eprint 1708.06502.
- Bora and Desai (2021a) K. Bora and S. Desai, JCAP 2021, 012 (2021a), eprint 2008.10541.
- Bora and Desai (2021b) K. Bora and S. Desai, JCAP 2021, 052 (2021b), eprint 2104.00974.
- Mendonça et al. (2021) I. E. C. R. Mendonça, K. Bora, R. F. L. Holanda, S. Desai, and S. H. Pereira, JCAP 2021, 034 (2021), eprint 2109.14512.
- Bora et al. (2022) K. Bora, R. F. L. Holanda, S. Desai, and S. H. Pereira, European Physical Journal C 82, 17 (2022), eprint 2106.15805.
- Feretti et al. (2012) L. Feretti, G. Giovannini, F. Govoni, and M. Murgia, Astronomy and Astrophysics Review 20, 54 (2012), eprint 1205.1919.
- Wik et al. (2014) D. R. Wik, A. Hornstrup, S. Molendi, G. Madejski, F. A. Harrison, A. Zoglauer, B. W. Grefenstette, F. Gastaldello, K. K. Madsen, N. J. Westergaard, et al., Astrophys. J. 792, 48 (2014), eprint 1403.2722.
- Rojas Bolivar et al. (2023) R. A. Rojas Bolivar, D. R. Wik, A. Tümer, F. Gastaldello, J. Hlavacek-Larrondo, P. Nulsen, V. Vacca, G. Madejski, M. Sun, C. L. Sarazin, et al., Astrophys. J. 954, 76 (2023), eprint 2308.00969.
- Keshet and Reiss (2018) U. Keshet and I. Reiss, Astrophys. J. 869, 53 (2018), eprint 1709.07442.
- Xi et al. (2018) S.-Q. Xi, X.-Y. Wang, Y.-F. Liang, F.-K. Peng, R.-Z. Yang, and R.-Y. Liu, Phys. Rev. D 98, 063006 (2018).
- Baghmanyan et al. (2022) V. Baghmanyan, D. Zargaryan, F. Aharonian, R. Yang, S. Casanova, and J. Mackey, MNRAS 516, 562 (2022), eprint 2110.00309.
- Pinzke and Pfrommer (2010) A. Pinzke and C. Pfrommer, MNRAS 409, 449 (2010), eprint 1001.5023.
- Colafrancesco et al. (2006) S. Colafrancesco, S. Profumo, and P. Ullio, Astron. & Astrophys. 455, 21 (2006), eprint astro-ph/0507575.
- Manna and Desai (2024) S. Manna and S. Desai, JCAP 2024, 017 (2024), eprint 2310.07519.
- Bleem et al. (2015) L. E. Bleem, B. Stalder, T. de Haan, K. A. Aird, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, B. A. Benson, et al., Astrophys. J. Suppl. Ser. 216, 27 (2015), eprint 1409.0850.
- Bocquet et al. (2019) S. Bocquet, J. P. Dietrich, T. Schrabback, L. E. Bleem, M. Klein, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, et al., Astrophys. J. 878, 55 (2019), eprint 1812.01679.
- Desai et al. (2012) S. Desai, R. Armstrong, J. J. Mohr, D. R. Semler, J. Liu, E. Bertin, S. S. Allam, W. A. Barkhouse, G. Bazin, E. J. Buckley-Geer, et al., Astrophys. J. 757, 83 (2012), eprint 1204.1210.
- Song et al. (2012) J. Song, A. Zenteno, B. Stalder, S. Desai, L. E. Bleem, K. A. Aird, R. Armstrong, M. L. N. Ashby, M. Bayliss, G. Bazin, et al., Astrophys. J. 761, 22 (2012), eprint 1207.4369.
- Sunyaev and Zeldovich (1972) R. A. Sunyaev and Y. B. Zeldovich, Comments on Astrophysics and Space Physics 4, 173 (1972).
- Schoenfelder et al. (1993) V. Schoenfelder, H. Aarts, K. Bennett, H. de Boer, J. Clear, W. Collmar, A. Connors, A. Deerenberg, R. Diehl, A. von Dordrecht, et al., Astrophys. J. Suppl. Ser. 86, 657 (1993).
- Iyudin et al. (2004) A. F. Iyudin, H. Böhringer, V. Dogiel, and G. Morfill, Astron. & Astrophys. 413, 817 (2004).
- Molero et al. (2023) M. Molero, P. Batista, L. Fortson, M. Nievas Rosillo, E. Pueschel, D. Ribeiro, and M. Vázquez Acosta, arXiv e-prints arXiv:2309.14159 (2023), eprint 2309.14159.
- Strong and Collmar (2019) A. Strong and W. Collmar, Memorie della Societa Astronomica Italiana 90, 297 (2019), eprint 1907.07454.
- Knödlseder et al. (2016) J. Knödlseder, M. Mayer, C. Deil, J. B. Cayrou, E. Owen, N. Kelley-Hoskins, C. C. Lu, R. Buehler, F. Forest, T. Louge, et al., Astron. & Astrophys. 593, A1 (2016), eprint 1606.00393.
- Knödlseder et al. (2022) J. Knödlseder, W. Collmar, M. Jarry, and M. McConnell, Astron. & Astrophys. 665, A84 (2022), eprint 2207.13404.
- Clowe et al. (2006) D. Clowe, M. Bradač, A. H. Gonzalez, M. Markevitch, S. W. Randall, C. Jones, and D. Zaritsky, Astrophys. J. Lett. 648, L109 (2006), eprint astro-ph/0608407.
- Asencio et al. (2023) E. Asencio, I. Banik, and P. Kroupa, Astrophys. J. 954, 162 (2023), eprint 2308.00744.
- Schönfelder et al. (2000) V. Schönfelder, K. Bennett, J. J. Blom, H. Bloemen, W. Collmar, A. Connors, R. Diehl, W. Hermsen, A. Iyudin, R. M. Kippen, et al., Astronomy and Astrophysics Supplement Series 143, 145 (2000), eprint astro-ph/0002366.
- Wood et al. (2017) M. Wood, R. Caputo, E. Charles, M. Di Mauro, J. Magill, J. S. Perkins, and Fermi-LAT Collaboration, in 35th International Cosmic Ray Conference (ICRC2017) (2017), vol. 301 of International Cosmic Ray Conference, p. 824, eprint 1707.09551.
- Mattox et al. (1996) J. R. Mattox, D. L. Bertsch, J. Chiang, B. L. Dingus, S. W. Digel, J. A. Esposito, J. M. Fierro, R. C. Hartman, S. D. Hunter, G. Kanbach, et al., Astrophys. J. 461, 396 (1996).
- Pasumarti and Desai (2022) V. Pasumarti and S. Desai, JCAP 2022, 002 (2022), eprint 2210.12804.
- Zoglauer et al. (2021) A. Zoglauer, T. Siegert, A. Lowell, B. Mochizuki, C. Kierans, C. Sleator, D. H. Hartmann, H. Lazar, H. Gulick, J. Beechert, et al., arXiv e-prints arXiv:2102.13158 (2021), eprint 2102.13158.
- Orlando et al. (2022) E. Orlando, E. Bottacini, A. A. Moiseev, A. Bodaghee, W. Collmar, T. Ensslin, I. V. Moskalenko, M. Negro, S. Profumo, S. W. Digel, et al., JCAP 2022, 036 (2022), eprint 2112.07190.
- Galper et al. (2014) A. M. Galper, V. Bonvicini, N. P. Topchiev, O. Adriani, R. L. Aptekar, I. V. Arkhangelskaja, A. I. Arkhangelskiy, L. Bergstrom, E. Berti, G. Bigongiari, et al., arXiv e-prints arXiv:1412.4239 (2014), eprint 1412.4239.
- Egorov et al. (2020) A. E. Egorov, N. P. Topchiev, A. M. Galper, O. D. Dalkarov, A. A. Leonov, S. I. Suchkov, and Y. T. Yurkin, JCAP 2020, 049 (2020), eprint 2005.09032.
- Boggs (2006) S. E. Boggs, New Astronomy Reviews 50, 604 (2006), eprint astro-ph/0608532.
- Hunter et al. (2014) S. D. Hunter, P. F. Bloser, G. O. Depaola, M. P. Dion, G. A. DeNolfo, A. Hanu, M. Iparraguirre, J. Legere, F. Longo, M. L. McConnell, et al., Astroparticle Physics 59, 18 (2014), eprint 1311.2059.
- Wu et al. (2014) X. Wu, M. Su, A. Bravar, J. Chang, Y. Fan, M. Pohl, and R. Walter, in Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray, edited by T. Takahashi, J.-W. A. den Herder, and M. Bautz (2014), vol. 9144 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, p. 91440F, eprint 1407.0710.
- Wu et al. (2015) X. Wu, J. Chang, R. Walter, M. Su, G. Ambrosi, M. Böttcher, M. Chernyakova, Y. Fan, C. Farnier, F. Gargano, et al., in 34th International Cosmic Ray Conference (ICRC2015) (2015), vol. 34 of International Cosmic Ray Conference, p. 964.
- Aramaki et al. (2020a) T. Aramaki, P. O. H. Adrian, G. Karagiorgi, and H. Odaka, Astroparticle Physics 114, 107 (2020a), eprint 1901.03430.
- Aramaki et al. (2020b) T. Aramaki, J. Asaadi, Y. Ichinohe, Y. Inoue, G. Karagiorgi, J. Leyva, R. Mukherjee, H. Odaka, K. Perez, W. Seligman, et al., arXiv e-prints arXiv:2009.03754 (2020b), eprint 2009.03754.
- Dzhatdoev and Podlesnyi (2019) T. Dzhatdoev and E. Podlesnyi, Astroparticle Physics 112, 1 (2019), eprint 1902.01491.
- Tatischeff et al. (2019) V. Tatischeff, A. De Angelis, M. Tavani, U. Oberlack, R. Walter, G. Ambrosi, A. Argan, P. von Ballmoos, S. Brandt, A. Bulgarelli, et al., Memorie della Societa Astronomica Italiana 90, 137 (2019), eprint 1905.07806.