Search For Gamma Rays Research Articles

Searching for dark matter annihilation with IceCube and P-ONE

We present a new search for weakly interacting massive particles utilizing ten years of public IceCube data, setting more stringent bounds than previous IceCube analysis on massive dark matter to neutrino annihilation. We also predict the future potential of the new neutrino observatory, P-ONE, showing that it may even exceed the sensitivities of Fermi-LAT gamma-ray searches by about 1–2 orders of magnitude in 1–10 TeV regions. This analysis considers the diffuse dark matter self-annihilation to neutrinos via direct and indirect channels, from the galactic dark matter halo and extra-galactic sources. We also predict that P-ONE will be capable of pushing these bounds further than IceCube, even reaching the thermal relic abundance utilizing a galactic center search for extended run-time.

Deep learning techniques for imaging air Cherenkov telescopes

Very-high-energy (VHE) gamma rays and charged cosmic rays (CCRs) provide an observational window into the acceleration mechanisms of extreme astrophysical environments. One of the major challenges at imaging air Cherenkov telescopes (IACTs) designed to look for VHE gamma rays, is the separation of air showers initiated by CCRs which form a background to gamma-ray searches. Two other less well-studied problems at IACTs are (a) the classification of different primary nuclei among the CCR events, and (b) identification of anomalous events initiated by beyond-the-Standard-Model (BSM) particles that could give rise to shower signatures which differ from the standard images of either gamma rays or CCR showers. The problems of categorizing the primary particle that initiates a shower image, or the problem of tagging anomalous shower events in a model-independent way, are problems that are well suited to a machine learning approach. Traditional studies that have explored gamma-ray/CCR separation have used a multivariate analysis based on derived shower properties, which contains significantly reduced information about the shower. In our work, we address the problems outlined above by using machine learning architectures trained on full simulated shower images, as opposed to training on just a few derived shower properties. We illustrate the techniques of binary and multicategory classification using convolutional neural networks, and we also pioneer the use of autoencoders for anomaly detection at VHE gamma-ray experiments. The latter technique has been studied previously in the context of collider physics, to tag anomalous BSM candidates in a model-independent way. In this study, for the first time, we demonstrate the efficacy of these techniques in the domain of VHE gamma-ray experiments. As a case study, we apply our techniques to the High Energy Stereoscopic System experiment. However, the real strength of the techniques that we broach here in the context of VHE gamma-ray observatories, is that these methods can be applied broadly to any other IACT---such as the upcoming Cherenkov Telescope Array---or can even be suitably adapted to CCR experiments.

Deep Search for Gamma-Ray Emission from the Accreting X-Ray Pulsar 1A 0535+262

Binary systems are a well-established subclass of gamma-ray sources. The high mass X-ray binary pulsar 1A 0535+262 has been considered to be a possible gamma-ray emitter for a long time, although former gamma-ray searches using the Fermi Large Area Telescope (LAT) and VERITAS data resulted in upper limits only. We aim at a deep search for gamma-ray emission and pulsations from 1A 0535+262 using more than 13 yr of Fermi-LAT data. The analysis was performed for both the whole Fermi-LAT data set, as well as for the X-ray outbursts that 1A 0535+262 has experienced since the launch of Fermi. Various X‐ray observations have been used to generate the ephemeris for the pulsation search. We also investigate the long-term gamma-ray flux variability and perform orbital phase-resolved analysis for the outbursts. We did not detect any steady or pulsed gamma-ray emission from 1A 0535+262 during the whole Fermi-LAT mission span or its X-ray outbursts. We thus derived the deepest gamma-ray luminosity upper limits to date at the 95% confidence level to be around (2.3–4.7) × 1032 erg s−1 depending on different spectral indices assumed, which results in a ratio of L γ to L X (2–150 keV) being (1.9–3.9) × 10−6.

New approaches for gamma-hadron separation at the IceCube Neutrino Observatory

The IceCube Neutrino Observatory located at the geographic South Pole is composed of two detectors. One is the in-ice optical array, which measures high-energy muons from air showers and charged particles produced by the interaction of high-energy neutrinos in the ice. The other is an array of ice-Cherenkov tanks at the surface, called IceTop, which is used both as veto for the in-ice neutrino measurements and for detecting cosmic-ray air showers. In the next decade, the IceCube-Gen2 extension will increase the surface coverage including surface radio antennas and scintillator panels on the footprint of an extended optical array in the ice. The combination of the current surface and in-ice detectors can be exploited for the study of cosmic rays and the search for PeV gamma rays. The in-ice detector measures the high-energy muonic component of air showers, whereas the signal in IceTop is dominated by the electromagnetic component. The relative size of the muonic and electromagnetic components is different for gamma-and hadron-induced air showers. Thus, the gamma-hadron separation of cosmic rays is attempted using machine learning techniques including deep learning. Here, different approaches are presented. Finally, the prospects for the detection of PeV photons with IceCube-Gen2 will be discussed.

Indirect detection of low-mass dark matter through the π0 and η windows

We consider the search for gamma-rays produced by the annihilation or decay of low-mass dark matter which couples to quarks. In this scenario, most of the photons are produced from the decays of $\pi^0$ or $\eta$ mesons. These decays produce distinctly different photon signatures due to the difference in meson mass. We assess the ability of the future MeV-range observatories to constrain the hadronic final states produced by dark matter annihilation or decay from the shape of the resulting photon spectrum. We then comment on how this information can be used to determine properties of the dark matter coupling to the quark current, based on the approximate symmetries of low-energy QCD.

Searching for velocity-dependent dark matter annihilation signals from extragalactic halos

We consider gamma-ray signals of dark matter annihilation in extragalactic halos in the case where dark matter annihilates from a p-wave or d-wave state. In these scenarios, signals from extragalactic halos are enhanced relative to other targets, such as the Galactic Center or dwarf spheroidal galaxies, because the typical relative speed of the dark matter is larger in extragalactic halos. We perform a mock data analysis of gamma rays produced by dark matter annihilation in halos detected by the Sloan Digital Sky Survey. We include a model for uncorrelated galactic and extragalactic gamma ray backgrounds, as well as a simple model for backgrounds due to astrophysical processes in the extragalactic halos detected by the survey. We find that, for models which are still allowed by other gamma ray searches, searches of extragalactic halos with the current Fermi exposure can produce evidence for dark matter annihilation, though it is difficult to distinguish the p-wave and d-wave scenarios. With a factor 10×larger exposure, though, discrimination of the velocity-dependence is possible.

Neutrino point source searches for dark matter spikes

Any dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation, leading to a potentially detectable neutrino signal. In this paper we examine 10 - 105 M ⊙ black holes associated with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today, in light of neutrino data from the ANTARES [1] and IceCube [2] detectors. In various regions of the sky, we determine the minimum distance away from the solar system that a dark matter spike must be in order to have not been detected as a neutrino point source for a variety of representative dark matter annihilation channels. Given these constraints on the distribution of dark matter spikes in the Galaxy, we place significant limits on the formation of the first generation of stars in early minihalos — stronger than previous limits from gamma-ray searches in Fermi Gamma-Ray Space Telescope data. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted; thus neutrino observations may be used to constrain the properties of Dark Stars. The limits are particularly strong for heavier WIMPs. For WIMP masses ∼ 5TeV, we show that ≲ 10 % of minihalos can host first stars that collapse into BHs larger than 103 M ⊙.

Spatial extension of dark subhalos as seen by Fermi -LAT and the implications for WIMP constraints

Spatial extension has been hailed as a "smoking gun" in the gamma-ray search of dark galactic subhalos, which would appear as unidentified sources for gamma-ray telescopes. In this work, we study the sensitivity of the Fermi-LAT to extended subhalos using simulated data based on a realistic sky model. We simulate spatial templates for a set of representative subhalos, whose parameters were derived from our previous work with N-body cosmological simulation data. We find that detecting an extended subhalo and finding an unequivocal signal of angular extension requires, respectively, a flux 2 to 10 times larger than in the case of a point-like source. By studying a large grid of models, where parameters such as the WIMP mass, annihilation channel or subhalo model are varied significantly, we obtain the response of the LAT as a function of the product of annihilation cross section times the J-factor. Indeed, we show that spatial extension can be used as an additional filter to reject subhalos candidates among the pool of unidentified LAT sources, as well as a "smoking gun" for positive identification. For instance, typical angular extensions of a few tenths of degree are expected for the considered scenarios. Finally, we also study the impact of the obtained LAT sensitivity to such extended subhalos on the achievable dark matter constraints, which are a few times less constraining than comparable point-source limits.

The MUSE-Faint survey. II. The dark matter-density profile of the ultra-faint dwarf galaxy Eridanus 2

Aims. We use stellar line-of-sight velocities to constrain the dark-matter density profile of Eridanus 2, an ultra-faint dwarf galaxy with an absolute V-band magnitude MV = −7.1 that corresponds to a stellar mass M* ≈ 9 × 104 M⊙. We furthermore derive constraints on fundamental properties of self-interacting and fuzzy dark matter scenarios. Methods. We present new observations of Eridanus 2 from MUSE-Faint, a survey of ultra-faint dwarf galaxies with the Multi Unit Spectroscopic Explorer on the Very Large Telescope, and determine line-of-sight velocities for stars inside the half-light radius. Combined with literature data, we have 92 stellar tracers out to twice the half-light radius. With these tracers we constrain models of cold dark matter, self-interacting dark matter, and fuzzy dark matter, using CJAM and pyGravSphere for the dynamical analysis. The models of self-interacting and fuzzy dark matter relate the density profile to the self-interaction coefficient and the dark-matter particle mass, respectively. Results. We find substantial evidence (Bayes factor ∼10−0.6) for cold dark matter (a cuspy halo) over self-interacting dark matter (a cored halo) and weak evidence (Bayes factor ∼10−0.4) for fuzzy dark matter over cold dark matter. We find a virial mass M200 ∼ 108 M⊙ and astrophysical factors J(αcJ) ~ 1011 M⊙2 kpc−5 and D(αcD) ~ 102 − 102.5 M⊙ kpc−2 (proportional to dark-matter annihilation and decay signals, respectively), the exact values of which depend on the density profile model. The mass-to-light ratio within the half-light radius is consistent with the literature. We do not resolve a core (rc < 47 pc, 68% confidence level) or a soliton (rsol < 7.2 pc, 68% confidence level). These limits are equivalent to an effective self-interaction coefficient fΓ < 2.2 × 10−29 cm3 s−1 eV−1 c2 and a fuzzy-dark-matter particle mass ma > 4.0 × 10−20 eV c−2. The constraint on self-interaction is complementary to those from gamma-ray searches. The constraint on fuzzy-dark-matter particle mass is inconsistent with those obtained for larger dwarf galaxies, suggesting that the flattened density profiles of those galaxies are not caused by fuzzy dark matter.

Search for gamma rays from SNe with a variable-size sliding-time-window analysis of the Fermi-LAT data

ABSTRACT We present a systematic search for γ-ray emission from supernovae (SNe) in the Fermi Large Area Telescope (LAT) Pass 8 data. The sample of targets consists of 55,880 candidates from the Open Supernova Catalogue. We searched for γ rays from SNe by means of a variable-size sliding-time-window analysis. Our results confirm the presence of transient γ-ray emission from the sources of non-AGN classes, including transitional pulsars, solar flares, γ-ray bursts, novae, and the Crab Nebula, which are projected near some of these SN’s positions, and also strengthen support to the variable signal in the direction of SN iPTF14hls. The analysis is successful in finding both short (e.g. solar flares) and long (e.g. transitional pulsars) high flux states. Our search reveals two new γ-ray transient signals occurred in 2019 in the directions of optical transients that are SN candidates, AT2019bvr and AT2018iwp, with their flux increases within 6 months after the dates of SN’s discoveries. These signals are bright and their variability is at a higher statistical level than that of iPTF14hls. An exploration of archival multiwavelength observations towards their positions is necessary to establish their association with SNe or other classes of sources. Our analysis, in addition, shows a bright transient γ-ray signal at low Galactic latitudes in the direction of PSR J0205+6449. In addition, we report the results of an all-sky search for γ-ray transient sources. This provided two additional candidates to gamma-ray transient sources.

Tunka-Grande and TAIGA-Muon scintillation arrays: status and prospects

The Tunka-Grande and TAIGA-Muon arrays are the part of a single experimental complex, which also includes the Tunka-133 and TAIGA-HiSCORE (High Sensitivity COsmic Rays and gamma Explorer) wide-angle Cherenkov arrays, TAIGA-IACT array (Imaging Atmospheric Cherenkov Telescope) and Tunka-Rex radio antennas array (Tunka Radio Extension). This complex is located in the Tunka Valley (Buryatia Republic, Russia), 50 km from Lake Baikal. It is aimed at investigating the energy spectrum and mass composition of charged cosmic rays in the energy range 100 TeV - 1000 PeV, searching for diffuse gamma rays above 100 TeV and studying local sources of gamma rays with energies above 30 TeV.This report outlines 3 key points. The first is a description of the Tunka-Grande and TAIGA-Muon scintillation arrays. The second part presents preliminary results of the search for diffuse gamma rays with energies above 50 PeV according to the Tunka-Grande data. The third part is devoted to the prospects of the search for diffuse gamma rays with energies above 100 TeV using the TAIGA-Muon array.

Carpet—2 Search for Gamma Rays above 100 TeV in Coincidence with HAWC and IceCube Alerts

We report on the search of astrophysical gamma rays with energies in the 100 TeV to several PeV range arriving in directional and temporal coincidence with public alerts from HAWC (TeV gamma rays) and IceCube (neutrinos above ∼ 100TeV). The observations have been performed with the Carpet-2 air-shower detector at the Baksan Neutrino Observatory, working in the “photon-friendly” mode since 2018. Photon candidate showers are selected by their low muon content. No significant excess of the photon candidates have been observed, and upper limits on gamma-ray fluences associated with the alerts are obtained. For events with good viewing conditions, the Carpet-2 effective area for photons is of the order of the IceCube effective area for neutrinos of the same energy, so the constraints start to probe the production of neutrinos in fast flares of Galactic sources.

Possibilities of the Tunka-Grande and TAIGA-Muon scintillation arrays with the TAIGA-HiSCORE Cherenkov array joint operation in the research of cosmic and gamma rays

The Tunka-Grande scintillation array is part of the TAIGA Gamma Observatory. It is intended for investigation of energy spectrum and mass composition of primary cosmic rays in the energy range 10 PeV–10 EeV and the search for diffuse cosmic gamma rays. The TAIGA-HiSCORE Cherenkov array aims at observing gamma-rays with the energy from 1 TeV. TAIGA-Muon low-threshold scintillation detector array is a network of surface and underground detectors for registration charge particles of EAS. Currently, 3 clusters have been deployed. The first cluster is running in test mode. It is planed that in the future the total area of the TAIGA-Muon will be about 2000 sq. m. and it will search astrophysical gamma-rays in the energy range from 100 TeV together with the Tunka-Grande scintillation array and the Cherenkov experiments of the TAIGA Gamma Observatory. To evaluate the possibility of join operation Tunka-Grande, TAIGA-Muon and TAIGA-HiSCORE, a simulation was performed using the CORSIKA and Geant4 software packages. The status of model-based studies is presented and assessed the prospects for joint operation of the arrays.

Search for PeV Gamma-Ray Emission from the Southern Hemisphere with 5 Yr of Data from the IceCube Observatory

Abstract The measurement of diffuse PeV gamma-ray emission from the Galactic plane would provide information about the energy spectrum and propagation of Galactic cosmic rays, and the detection of a pointlike source of PeV gamma-rays would be strong evidence for a Galactic source capable of accelerating cosmic rays up to at least a few PeV. This paper presents several unbinned maximum-likelihood searches for PeV gamma-rays in the Southern Hemisphere using 5 yr of data from the IceTop air shower surface detector and the in-ice array of the IceCube Observatory. The combination of both detectors takes advantage of the low muon content and deep shower maximum of gamma-ray air showers and provides excellent sensitivity to gamma-rays between ∼0.6 and 100 PeV. Our measurements of pointlike and diffuse Galactic emission of PeV gamma-rays are consistent with the background, so we constrain the angle-integrated diffuse gamma-ray flux from the Galactic plane at 2 PeV to 2.61 × 10−19 cm−2 s−1 TeV−1 at 90% confidence, assuming an E −3 spectrum, and we estimate 90% upper limits on pointlike emission at 2 PeV between 10−21 and 10−20 cm−2 s−1 TeV−1 for an E −2 spectrum, depending on decl. Furthermore, we exclude unbroken power-law emission up to 2 PeV for several TeV gamma-ray sources observed by the High Energy Spectroscopic System and calculate upper limits on the energy cutoffs of these sources at 90% confidence. We also find no PeV gamma-rays correlated with neutrinos from IceCube’s high-energy starting event sample. These are currently the strongest constraints on PeV gamma-ray emission.

Search for cosmic gamma rays of ultra-high energies with the Baksan air shower array

Recently two collaborations (Tibet ASγ and HAWC) almost simultaneously reported about their measurements of ultra-high-energy (> 100 TeV) gamma rays from the Crab Nebula. The Baksan air shower array (BASA) recorded a burst from this very source as long ago as more than 20 years. At the moment the new experiment on gamma ray astronomy (Carpet-3 project) started at Baksan, and it is appropriate to recall old results of BASA to be compared with the new results, already published and expected.

CARPET-2 SEARCH FOR GAMMA RAYS ABOVE 100 TEV IN COINCIDENCE WITH HAWC AND ICECUBE ALERTS

CARPET-2 SEARCH FOR GAMMA RAYS ABOVE 100 TEV IN COINCIDENCE WITH HAWC AND ICECUBE ALERTS

Neutrinos and gamma rays from long-lived mediator decays in the Sun

We investigate a scenario where dark matter (DM) particles can be captured and accumulate in the Sun, and subsequently annihilate into a pair of long-lived mediators. These mediators can decay further out in the Sun or outside of the Sun. Compared to the standard scenario where DM particles annihilate directly into Standard Model particles close to the solar core, here we also obtain fluxes of gamma rays and charged cosmic rays. We simulate this scenario using a full three-dimensional model of the Sun, and include interactions and neutrino oscillations. In particular, we perform a model-independent study of the complementarity between neutrino and gamma ray fluxes by comparing the recent searches from IceCube, Super-Kamiokande, Fermi-LAT, ARGO and HAWC. We find that the resulting neutrino fluxes are significantly higher at high energy when the mediators decay further out in the Sun. We also find that gamma ray searches place stronger constraints than neutrino searches on these models even in cases where the mediators decay mainly inside the Sun, except in the approximately inner 10% of the Sun where neutrino searches are more powerful. We present our results in a model-independent manner and release a new version of the \ extsf{WimpSim} code that can be used to simulate this scenario for arbitrary mediator models.

Left–right symmetric models with a mixture of keV–TeV dark matter

We discuss the possibility of realising a multi-component dark matter (DM) scenario with widely separated DM masses: one having keV scale mass and the other with GeV–TeV scale mass, within the framework of left right symmetric models. Due to gauge interactions, both the DM candidates are produced thermally in the early Universe but overproducing the keV mass candidate. We consider one of the right handed neutrinos to be decaying at late epochs, just before the big bang nucleosynthesis, in order to dilute the thermally overproduced keV DM. We constrain the parameter space from the requirement of producing sub-dominant keV–TeV DM, satisfying indirect detection constraints from gamma ray searches and producing the tantalising 3.55 keV monochromatic x-ray line, reported by several groups to be present in galaxy and galaxy cluster data, from the decay of a 7.1 keV DM on cosmological scales. We find that these requirements can keep the right sector gauge boson masses around a few TeV while requiring some of the right handed neutrinos in the sub GeV regime.

The Search for Diffuse Gamma Rays Using Data from the Tunka-Grande Experiment

The Tunka-Grande array is part of an experimental complex located in the Tunka Valley (Republic of Buryatia, Russia) about 50 km from Lake Baikal. This complex also contains the Tunka-133 and Tunka-Rex arrays. The aim of this complex is to study the primary cosmic ray energy spectrum and mass composition in the energy range of 1016–1018 eV, and to search for diffuse gamma rays in the energy range of 5 × 1016–5 × 1017 eV. The design of the Tunka-Grande array and the procedure for reconstructing the parameters of extensive air showers (EASes) are described, and preliminary results are presented from the search for diffuse gamma rays with energies of more than 5 × 1016 eV.

Dependence of accessible dark matter annihilation cross sections on the density profiles of dwarf spheroidal galaxies with the Cherenkov Telescope Array

Dwarf spheroidal galaxies are excellent targets in gamma-ray searches for dark matter. We consider dark matter searches in dwarf spheroidal galaxies (dSphs) with the Cherenkov Telescope Array (CTA). The aim of this work is to reveal a quantitative and precise dependence of the accessible dark matter annihilation cross-sections on the dark matter density profiles of dSphs and on the distance to them. In most data analyses, researchers have assumed point-like signals from dSphs because it is difficult to resolve the expected emission profiles with current gamma-ray observatories. In future however, CTA will be able to resolve the peak emission profiles in dSphs. We take several variations of the dark matter density profile of Draco dSph as examples and analyze the simulated observations of with CTA. We derive the accessible region of the dark matter annihilation cross-section with each dark matter density profile. The accessible region of the annihilation cross-section can differ by a factor of 10 among plausible profiles. We also examine the dependence on the distance to the target dSphs by assuming the same profiles of dSphs at different distances. Closer targets are better due to the higher J-factor, while their spatial extension significantly degrades our reach to the annihilation cross-section compared to the value expected from a simple distance-scaling of the J-factor. Spatial extension of the source affects the probable parameter region in energy-dependent ways. In some gamma-ray energy ranges, this behaviour becomes moderately dependent on the properties of the observation facility.