PeV Energies Research Articles

Constraining the Thickness of the Galactic Halo through Cosmic-Ray Anisotropy Using the Spatial-Dependent-Propagation Model

The spatial-dependent-propagation (SDP) model with a nearby source works well to reproduce the coevolving features of both cosmic-ray (CR)-nuclei spectra and anisotropy. However, it is well known that the Sun is actually deviating from the galactic disk. This will lead to a dominating anisotropy in the direction perpendicular to the galactic disk, which is discrepant with current observations. Thus, it is necessary to further investigate the effect of the solar offset on anisotropy. In this work, to the best of our knowledge, this is the first time that the combined studies of the solar offset, nuclei spectra, and anisotropy have been performed based on the SDP model. As a result, to reproduce CR spectra and anisotropy, the thickness of the inner halo (ZIH) needs to increase linearly with the displacement of the Sun. We also know that the PeV anisotropy could be used to estimate the value of the diffusion coefficient, thus breaking the degeneracy between the diffusion coefficient and halo thickness. Therefore, it is a good approach to constrain the halo thickness. Moreover, the anisotropy in the PeV energy region, as a new probe, might also shed new light on constraining the solar offset. It is hoped that the anisotropy of the energies from ∼TeV to PeV can be finely measured with the LHAASO experiment, leading to a better understanding of the thick halo.

Potential for Constraining Propagation Parameters of Galactic Cosmic Rays with the High Energy Cosmic-radiation Detection Facility on Board China’s Space Station

Precise measurements of the spectra of secondary and primary cosmic rays are crucial for understanding the origin and propagation of those energetic particles. The High Energy Cosmic-radiation Detection (HERD) facility on board China’s Space Station, which is expected to operate in 2027, will push the direct and precise measurements of cosmic-ray fluxes up to PeV energies. In this work, we investigate the potential of HERD for studying the propagation of cosmic rays using measurements of boron, carbon, and oxygen spectra. We find that, compared with the current results, the new HERD measurements can improve the accuracy of the propagation parameters by 8%–40%. The constraints on the injection spectra at high energies will also be improved.

Prediction of using LHAASO’s cosmic-ray electron measurements to constrain decaying heavy dark matter

Large High Altitude Air Shower Observatory (LHAASO) is an instrument designed for detecting cosmic rays (CRs) and gamma rays at TeV to PeV energies. The decays of heavy dark matter particles in the Galactic halo may produce high-energy electrons that can be detected by LHAASO. The main background for the LHAASO's CR electron measurements is the hadron residuals due to misidentification of the particle species. In this paper, we estimate the LHAASO's electron background using the known all-particle CR spectrum and the hadron rejection efficiency of LHAASO. With the estimated background, we predict the capability of LHAASO to constrain dark matter (DM) decay lifetime at 95% confidence level for various channels. We find that if neglecting systematic uncertainties, the CR electron measurement by LHAASO can improve the current best results by up to one order of magnitude for DM masses between 100--1000 TeV. However, indirect measurements of CR electrons by ground-based experiments suffer from sizable systematic uncertainties. With the systematic uncertainties included in the calculation, the projected constraints will be largely weakened. So, for using the CR electron observation of LHAASO to constrain the DM parameters, the key point is whether the systematic error can be effectively reduced.

PeV-scale leptogenesis, gravitational waves, and black holes from a SUSY-breaking phase transition

Supersymmetry is a highly motivated theoretical framework, whose scale of breaking may be at PeV energies, to explain null searches at the Large Hadron Collider. SUSY breaking through a first order phase transition may have occurred in the early universe, leading to potential gravitational wave signals. Constructing a realistic model for gauge-mediated supersymmetry breaking, we show that such a transition can also induce masses for heavy right-handed neutrinos and sneutrinos, whose CP-violating decays give leptogenesis at the PeV scale, and a novel mechanism of neutrino mass generation at one loop. For the same models we predict the possible gravity wave signals, and we study the possibility of production of primordial black holes during the phase transition.

Hint for a TeV neutrino emission from the Galactic Ridge with ANTARES

Interactions of cosmic ray protons, atomic nuclei, and electrons in the interstellar medium in the inner part of the Milky Way produce a γ-ray flux from the Galactic Ridge. If the γ-ray emission is dominated by proton and nuclei interactions, a neutrino flux comparable to the γ-ray flux is expected from the same sky region.Data collected by the ANTARES neutrino telescope are used to constrain the neutrino flux from the Galactic Ridge in the 1-100 TeV energy range. Neutrino events reconstructed both as tracks and showers are considered in the analysis and the selection is optimized for the search of an excess in the region |l|<30°, |b|<2°. The expected background in the search region is estimated using an off-zone region with similar sky coverage. Neutrino signal originating from a power-law spectrum with spectral index ranging from Γν=1 to 4 is simulated in both channels. The observed energy distributions are fitted to constrain the neutrino emission from the Ridge.The energy distributions in the signal region are inconsistent with the background expectation at ∼96% confidence level. The mild excess over the background is consistent with a neutrino flux with a power law with a spectral index 2.45−0.34+0.22 and a flux normalization dNνdEν=4.0−2.0+2.7×10−16 GeV−1cm−2s−1sr−1 at 40 TeV reference energy. Such flux is consistent with the expected neutrino signal if the bulk of the observed γ-ray flux from the Galactic Ridge originates from interactions of cosmic ray protons and nuclei with a power-law spectrum extending well into the PeV energy range.

Multiwavelength study of the galactic PeVatron candidate LHAASO J2108+5157

Context. Several new ultrahigh-energy (UHE) γ-ray sources have recently been discovered by the Large High Altitude Air Shower Observatory (LHAASO) collaboration. These represent a step forward in the search for the so-called Galactic PeVatrons, the enigmatic sources of the Galactic cosmic rays up to PeV energies. However, it has been shown that multi-TeV γ-ray emission does not necessarily prove the existence of a hadronic accelerator in the source; indeed this emission could also be explained as inverse Compton scattering from electrons in a radiation-dominated environment. A clear distinction between the two major emission mechanisms would only be made possible by taking into account multi-wavelength data and detailed morphology of the source. Aims. We aim to understand the nature of the unidentified source LHAASO J2108+5157, which is one of the few known UHE sources with no very high-energy (VHE) counterpart. Methods. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good-quality data. In addition, we analyzed 12 years of Fermi-LAT data, to better constrain emission of its high-energy (HE) counterpart 4FGL J2108.0+5155. We used naima and jetset software packages to examine the leptonic and hadronic scenario of the multi-wavelength emission of the source. Results. We found an excess (3.7σ) in the LST-1 data at energies E &gt; 3 TeV. Further analysis of the whole LST-1 energy range, assuming a point-like source, resulted in a hint (2.2σ) of hard emission, which can be described with a single power law with a photon index of Γ = 1.6 ± 0.2 the range of 0.3 − 100 TeV. We did not find any significant extended emission that could be related to a supernova remnant (SNR) or pulsar wind nebula (PWN) in the XMM-Newton data, which puts strong constraints on possible synchrotron emission of relativistic electrons. We revealed a new potential hard source in Fermi-LAT data with a significance of 4σ and a photon index of Γ = 1.9 ± 0.2, which is not spatially correlated with LHAASO J2108+5157, but including it in the source model we were able to improve spectral representation of the HE counterpart 4FGL J2108.0+5155. Conclusions. The LST-1 and LHAASO observations can be explained as inverse Compton-dominated leptonic emission of relativistic electrons with a cutoff energy of 100−30+70 TeV. The low magnetic field in the source imposed by the X-ray upper limits on synchrotron emission is compatible with a hypothesis of a PWN or a TeV halo. Furthermore, the spectral properties of the HE counterpart are consistent with a Geminga-like pulsar, which would be able to power the VHE-UHE emission. Nevertheless, the lack of a pulsar in the neighborhood of the UHE source is a challenge to the PWN/TeV-halo scenario. The UHE γ rays can also be explained as π0 decay-dominated hadronic emission due to interaction of relativistic protons with one of the two known molecular clouds in the direction of the source. Indeed, the hard spectrum in the LST-1 band is compatible with protons escaping a shock around a middle-aged SNR because of their high low-energy cut-off, but the origin of the HE γ-ray emission remains an open question.

Bump hunting in the diffuse flux of high-energy cosmic neutrinos

The origin of the bulk of the high-energy astrophysical neutrinos seen by IceCube, with TeV--PeV energies, is unknown. If they are made in photohadronic, i.e., proton-photon, interactions in astrophysical sources, this may manifest as a bump-like feature in their diffuse flux, centered around a characteristic energy. We search for evidence of this feature, allowing for variety in its shape and size, in 7.5 years of High-Energy Starting Events (HESE) collected by the IceCube neutrino telescope, and make forecasts using larger data samples from upcoming neutrino telescopes. Present-day data reveals no evidence of bump-like features, which allows us to constrain candidate populations of photohadronic neutrino sources. Near-future forecasts show promising potential for stringent constraints or decisive discovery of bump-like features. Our results provide new insight into the origins of high-energy astrophysical neutrinos, complementing those from point-source searches.

Multiwavelength analysis of Galactic Supernova Remnants

The origin of Galactic Cosmic Rays (CRs) and the possibility of Supernova Remnants (SNRs) being potential CR accelerators is still an open debate. The charged CRs can be detected indirectly by the γ-ray observatories through the π 0 production and consequent decay, leading to the generation of high-energy γ-rays. The goal of the study is to identify qualitative and quantitative trends in favour of hadronic scenario and search for SNRs which could be potential accelerators up to PeV energies (PeVatrons).We have performed a Multiwavelength (MWL) study using different radiative models to evaluate the hadronic contribution. The spectral energy distributions (SEDs) of selected SNRs are modeled using the Naima [1] package. Two different radiative scenarios are considered, pure leptonic and lepto-hadronic scenarios and different methods are used to evaluate their importance.This study shows that the lepto-hadronic scenario is favored for most SNRs. Two particular indicators of hadronic contribution come from the data around the π 0 production threshold and the data above a few TeV. The hard rise at the π 0 production threshold cannot be explained by leptonic processes. More data in this region would be valuable for these studies. For some SNRs, an important hadronic contribution is observed up to a few TeV, thus making them promising PeVatron candidates. In this high-energy region where the leptonic processes are expected to be suppressed, more data is required to help distinguish between the leptonic and hadronic origin of γ-ray emission. In the future, we intend to use the obtained model parameters to simulate data for CTA and assess its capability to identify PeVatrons.

Probing heavy sterile neutrinos at neutrino telescopes via the dipole portal

The question of how heavy a sterile neutrino can be probed in experiments leads us to investigate the Primakoff production of heavy sterile neutrinos from ultrahigh-energy neutrinos via the magnetic dipole portal. Despite the suppression from the small magnetic moment, the transition is significantly enhanced by tiny t-channel momentum transfers, similar to the resonant production of pions and axions in an external electromagnetic field. Based on the current IceCube measurement of astrophysical neutrinos up to PeV energies, strong constraints can already be derived on the transition magnetic moments of sterile neutrinos up to TeV masses. Moreover, we investigate the sensitivity of future tau neutrino telescopes, which are designed for EeV cosmogenic neutrino detection. We find that sterile neutrino masses as large as 30TeV can be probed at tau neutrino telescopes such as GRAND, POEMMA, and Trinity.

Searching Lorentz invariance violation from cosmic photon attenuation

Lorentz invariance violation (LIV) can change the threshold behavior predicted by special relativity and cause threshold anomalies which affect the propagation of cosmic photons. In this work, we focus on the threshold anomaly effect on cosmic photon attenuations by extragalactic background light (EBL) and discuss how to identify LIV from observations of very high energy (VHE) photons propagated from long distance in the universe. We point out that the Large High Altitude Air Shower Observatory (LHAASO), one of the most sensitive gamma-ray detector arrays currently operating at TeV and PeV energies, is an ideal facility for performing such LIV searching. We apply the proposed strategy to discuss the newly observed gamma-ray burst GRB 221009A to demonstrate the predictive ability of our suggestions.

X-Ray Studies of the Pulsar PSR J1420–6048 and Its TeV Pulsar Wind Nebula in the Kookaburra Region

We present a detailed analysis of broadband X-ray observations of the pulsar PSR J1420−6048 and its wind nebula (PWN) in the Kookaburra region with Chandra, XMM-Newton, and NuSTAR. Using the archival XMM-Newton and new NuSTAR data, we detected 68 ms pulsations of the pulsar and characterized its X-ray pulse profile, which exhibits a sharp spike and a broad bump separated by ∼0.5 in phase. A high-resolution Chandra image revealed a complex morphology of the PWN: a torus-jet structure, a few knots around the torus, one long (∼7′) and two short tails extending in the northwest direction, and a bright diffuse emission region to the south. Spatially integrated Chandra and NuSTAR spectra of the PWN out to 2.′5 are well-described by a power-law model with a photon index Γ ≈ 2. A spatially resolved spectroscopic study, as well as NuSTAR radial profiles of the 3–7 keV and 7–20 keV brightness, showed a hint of spectral softening with increasing distance from the pulsar. A multiwavelength spectral energy distribution (SED) of the source was then obtained by supplementing our X-ray measurements with published radio, Fermi-LAT, and H.E.S.S. data. The SED and radial variations of the X-ray spectrum were fit with a leptonic multizone emission model. Our detailed study of the PWN may be suggestive of (1) particle transport dominated by advection, (2) a low magnetic-field strength (B ∼ 5 μG), and (3) electron acceleration to ∼PeV energies.

A GEANT4 based simulation framework for the large area muon telescope of the GRAPES-3 experiment

The GRAPES-3 experiment located in Ooty, India, samples the electron and muon components in extensive air showers using an array of plastic scintillator detectors and a muon telescope (G3MT) consisting of proportional counters to study the composition of primary cosmic rays (PCRs) as well as γ-ray sources in the TeV–PeV energy range. The G3MT is designed with an appropriate mass absorber to shield the electromagnetic and hadronic components in the shower and to detect muons above 1 GeV×sec(θ) energy, incident from a zenith angle θ. We developed a simulation framework based on the GEANT4 toolkit to evaluate the response of shower particles such as muons, γ-rays, electrons and hadrons in the G3MT. We discuss the geometric modeling of the G3MT using GEANT4 starting with the proportional counter. We estimated the punch-through contribution of hadrons in the G3MT. We compare the simulated muon multiplicity distributions with the observed ones assuming PCR composition from a four population supernova remnant acceleration model namely H4a.

Limits on Neutrino Emission from GRB 221009A from MeV to PeV Using the IceCube Neutrino Observatory

Gamma-ray bursts (GRBs) have long been considered a possible source of high-energy neutrinos. While no correlations have yet been detected between high-energy neutrinos and GRBs, the recent observation of GRB 221009A—the brightest GRB observed by Fermi-GBM to date and the first one to be observed above an energy of 10 TeV—provides a unique opportunity to test for hadronic emission. In this paper, we leverage the wide energy range of the IceCube Neutrino Observatory to search for neutrinos from GRB 221009A. We find no significant deviation from background expectation across event samples ranging from MeV to PeV energies, placing stringent upper limits on the neutrino emission from this source.

MAGIC observations provide compelling evidence of hadronic multi-TeV emission from the putative PeVatron SNR G106.3+2.7

Context.Certain types of supernova remnants (SNRs) in our Galaxy are assumed to be PeVatrons, capable of accelerating cosmic rays (CRs) to ~ PeV energies. However, conclusive observational evidence for this has not yet been found. The SNR G106.3+2.7, detected at 1–100 TeV energies by different γ-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape, which can be divided into a head and a tail region with different physical conditions. However, in which region the 100 TeV emission is produced has not yet been identified because of the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear as to whether the origin of the γ-ray emission is leptonic or hadronic.Aims.With the better angular resolution provided by new MAGIC data compared to earlierγ-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1° resolution at TeV energies.Methods.We observed the SNR G106.3+2.7 using the MAGIC telescopes for 121.7 h in total – after quality cuts – between May 2017 and August 2019. The analysis energy threshold is ~0.2 TeV, and the angular resolution is 0.07−0.1°. We examined theγ-ray spectra of different parts of the emission, whilst benefitting from the unprecedented statistics and angular resolution at these energies provided by our new data. We also used measurements at other wavelengths such as radio, X-rays, GeVγ-rays, and 10 TeVγ-rays to model the emission mechanism precisely.Results.We detect extended γ-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detect a significantγ-ray emission with energies above 6.0 TeV from only the tail region suggests that the emissions above 10 TeV detected with air shower experiments (Milagro, HAWC, Tibet ASγand LHAASO) are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of ~1 PeV for that region. Such high-energy emission in this middle-aged SNR (4−10 kyr) can be explained by considering a scenario where protons escaping from the SNR in the past interact with surrounding dense gases at present.Conclusions.Theγ-ray emission region detected with the MAGIC telescopes in the SNR G106.3+2.7 is extended and spatially coincident with the radio continuum morphology. The multi-wavelength spectrum of the emission from the tail region suggests proton acceleration up to ~PeV, while the emission mechanism of the head region could either be hadronic or leptonic.

Probing LHAASO galactic PeVatrons through gamma-ray and neutrino correspondence

ABSTRACT Recently, Large High Altitude Air Shower Observatory (LHAASO) has detected several Galactic point sources of ultra high energy (UHE; Eγ &amp;gt; 100 TeV) gamma-rays. These gamma-rays are possibly created in leptonic or hadronic interactions of cosmic rays (CRs) of PeV energies. In the hadronic channel (p − p interaction), the gamma-rays are accompanied by neutrinos. The detection of neutrinos is therefore crucial in understanding CR acceleration in such objects. To estimate the neutrino flux, we adopt the two LHAASO sources (J2226+6057, J1908+0621) found to be spatially associated with the supernova remnants (SNR G106.3+2.7, SNR G40.5–0.5). For these two sources, the detected TeV-PeV gamma-ray spectra are found to be unusually hard (with spectral index ∼1.8). We develop a model of gamma-ray and neutrino emission based on the above two prototypes. The neutrino fluxes from these two sources are found to be below the IceCube sensitivity, but are detectable in upcoming IceCube-Gen2 and KM3NeT experiments. We further estimate the neutrino fluxes from similar other 10 LHAASO PeVatron sources and most of them are found to be detectable in IceCube-Gen2 and KM3NeT. Finally, we explore our model parameters, in particular the spectral power-law index and estimate the future potential of the neutrino detectors to probe CR acceleration in such Galactic sources.

Baikal-GVD neutrino telescope: Design reference 2022

Baikal-GVD is a large-scale Cherenkov detector designed to detect neutrinos in the TeV–PeV energy range with a goal to establish their sources. It is being deployed cluster by cluster and has been taking data since 2015. This contribution presents the telescope structure as of 2022.

Supernova connection of unidentified ultra-high-energy gamma-ray source LHAASO J2108+5157

ABSTRACT We present a simple phenomenological model of hadronic interaction between protons accelerated in an old supernova remnant (SNR) and cold protons situated within the associated molecular clouds (MCs). The accelerated protons from the old SNR escaped the SNR shock front, and got injected into the MCs at an earlier time, producing ultra-high-energy gamma-rays and neutrinos through inelastic proton–proton interaction. We also take into account the acceleration and subsequent escape of electrons from the SNR shock front. The escaped electrons produce gamma-rays through various radiative cooling mechanisms, after getting injected into the MCs. We use the model discussed in this letter to explain the multiwavelength spectral energy distribution of unidentified Galactic ultra-high-energy gamma-ray source LHAASO J2108+5157. We also discuss the feasibility of applying this model in other cases as well. Future observations can test the viability of the model discussed in this letter, which will in turn confirm that the SNRs can, in fact, accelerate particles up to PeV energies.

Where are Milky Way’s hadronic PeVatrons?

Observations of the Milky Way at TeV--PeV energies reveal a bright diffuse flux of hadronic cosmic rays and also bright point sources of gamma rays. If the gamma-ray sources are hadronic cosmic-ray accelerators, then they must also be neutrino sources. However, no neutrino sources have been detected. Where are they? We introduce a new population-based approach to probe Milky Way hadronic PeVatrons, demanding consistency between diffuse and point-source PeV-range data on cosmic rays, gamma rays, and neutrinos. For the PeVatrons, two extreme scenarios are allowed: (1) the hadronic cosmic-ray accelerators and the gamma-ray sources are the same objects, so that bright neutrino sources exist and improved telescopes can detect them, versus (2) the hadronic cosmic-ray accelerators and the gamma-ray sources are distinct, so that there are no detectable neutrino sources. The latter case is possible if hadronic accelerators have sufficiently thin column densities. We quantify present constraints and future prospects, showing how to reveal the nature of the hadronic PeVatrons.

Current Status of the HERD Space Mission Design

HERD (High Energy Radiation Detector) is a future experiment for space-based detection of cosmic rays and gamma ray astronomy onboard the Chinese Space Station. Its most innovative feature will be the event collection from 5 sides, and with its large size it will result in a one-order-of-magnitude jump in acceptance with respect to current largest calorimetric experiments. This will make possible to investigate cosmic ray spectra for each species up to the so-called knee region at PeV energies. The all-electron spectrum will be also measured up to about 10 TeV (depending on the actual flux) and the gamma sky will be studied with a large effective area from a few hundreds MeV up to few TeV, also allowing for the search for possible signatures from dark matter particles in our galaxy.

Upper Limits on the Isotropic Diffuse Flux of Cosmic PeV Photons from Carpet-2 Observations

Isotropic diffuse gamma-ray flux in the PeV energy band is an important tool for multimessenger tests of models of the origin of high-energy astrophysical neutrinos and for new-physics searches. So far, this flux has not yet been observed. Carpet-2 is an air-shower experiment capable of detecting astrophysical gamma rays with energies above 0.1 PeV. Here we report the upper limits on the isotropic gamma-ray flux from Carpet-2 data obtained in 1999-2011 and 2018-2022. These results, obtained with the new statistical method based on the shape of the muon-number distribution, summarize Carpet-2 observations as the upgraded installation, Carpet-3, starts its operation.