High Energy Cosmic Rays Research Articles

Cosmic rays from annihilation of heavy dark matter particles

The origin of the ultra high energy cosmic rays via annihilation of heavy stable, fermions “f”, of the cosmological dark matter (DM) is studied. The particles in question are supposed to be created by the scalaron decays in R2 modified gravity. The novel part of our approach is the assumption that the mass of these carriers of DM is slightly below than a half of the scalaron mass. In such a case the phase space volume becomes tiny. It leads to sufficiently low probability of “f” production, so their average cosmological energy density could be equal to the observed energy density of dark matter. Several regions of the universe, where the annihilation could take place, are studied. They include the whole universe under the assumption of homogeneous energy density, the high density DM clump in the galactic center, the cloud of DM in the Galaxy with realistic density distribution, and high density clumps of DM in the Galaxy. Possible resonance annihilation of ff¯ into energetic light particles is considered. We have shown that the proposed scenario can successfully explain the origin of the ultrahigh energy flux of cosmic rays where canonical astrophysical mechanisms are not operative.

Hadron-ion collisions in Pythia and the vector-meson dominance model for photoproduction

We present an extension to the Pythia Monte Carlo event generator that enables simulations of collisions between a generic hadron beam on a nuclear target with energy variation in event-by-event basis. This builds upon Pythia ’s module for heavy ions, Angantyr, as well as previous work on simulating hadron-proton collisions. As such, the extensions in this work are largely technical, except for a rudimentary model for hadronic fluctuations. With hadron-ion simulations, we implement an explicit vector-meson dominance (VMD) model that can be used to simulate interactions of hadronic component of real photons in photo-nuclear collisions. Such processes can be studied in ultra-peripheral heavy-ion collisions and in the future also with the upcoming Electron-Ion Collider. Our work also has applications to hadronic showers, e.g. air showers initiated by high-energy cosmic rays. We first validate the VMD model by comparing to HERA photoproduction data on proton target. Then we apply this to generate events for ultra-peripheral heavy-ion collisions at the LHC and present the results corresponding to the event-selection criteria matching to a recent ATLAS analysis. We find that single-particle multiplicity and rapidity distributions are well in line with the measured ones. We also construct the Fourier coefficients from two-particle correlations for the simulated events and study whether the resulting azimuthal anisotropies are consistent with the ATLAS results.

Search for Joint Multimessenger Signals from Potential Galactic Cosmic-Ray Accelerators with HAWC and IceCube

The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

Large-scale Cosmic-ray Anisotropies with 19 yr of Data from the Pierre Auger Observatory

Results are presented for the measurement of large-scale anisotropies in the arrival directions of ultra–high-energy cosmic rays detected at the Pierre Auger Observatory during 19 yr of operation, prior to AugerPrime, the upgrade of the observatory. The 3D dipole amplitude and direction are reconstructed above 4 EeV in four energy bins. Besides the established dipolar anisotropy in R.A. above 8 EeV, the Fourier amplitude of the 8–16 EeV energy bin is now also above the 5σ discovery level. No time variation of the dipole moment above 8 EeV is found, setting an upper limit to the rate of change of such variations of 0.3% yr−1 at the 95% confidence level. Additionally, the results for the angular power spectrum are shown, demonstrating no other statistically significant multipoles. The results for the equatorial dipole component down to 0.03 EeV are presented, using for the first time a data set obtained with a trigger that has been optimized for lower energies. Finally, model predictions are discussed and compared with observations, based on two source emission scenarios obtained in the combined fit of spectrum and composition above 0.6 EeV.

Global Analysis of the Extended Decreases in Cosmic Rays Observed with Worldwide Networks of Neutron Monitors and Muon Detectors: Temporal Variation of the Rigidity Spectrum and Its Implication

This paper presents the global analysis of two extended decreases in the galactic cosmic-ray intensity observed by worldwide networks of ground-based detectors in 2012. This analysis is capable of separately deriving the cosmic-ray density (or omnidirectional intensity) and anisotropy, each as a function of time and rigidity. A simple diffusion model along the spiral field line between Earth and a cosmic-ray barrier indicates the long duration of these events, resulting from about 190° eastern extent of a barrier such as an interplanetary shock followed by the sheath region and/or the corotating interaction region (CIR). It is suggested that the coronal mass ejection merging with and compressing the preexisting CIR at its flank can produce such an extended barrier. The derived rigidity spectra of the density and anisotropy both vary in time during each event period. In particular we find that the temporal feature of the “phantom Forbush decrease (FD)” reported in an analyzed period is dependent on rigidity, and looks quite different at different rigidities. From these rigidity spectra of the density and anisotropy, we derive the rigidity spectrum of the average parallel mean free path of pitch angle scattering along the spiral field line and infer the power spectrum of the magnetic fluctuation and its temporal variation. The possible physical cause of the strong rigidity dependence of the phantom FD is also discussed. These results demonstrate the high-energy cosmic rays observed at Earth responding to remote space weather.

Lateral Distribution Function Estimation of Electrons and Muons using Nishimura-Kamata-Greisen Function

Abstract The density of charged particles in extensive air showers reaching the surface of the Earth was calculated by estimating the lateral distribution function (LDF) of different primary particles at high energies. LDF simulation was performed using the AIRES system (version 19.04.10), a set of programs and subroutines designed to simulate ultra-high-energy air showers resulting from the interaction of cosmic rays with the Earth’s atmosphere. This system includes algorithms for fast simulation and output data management, and can simulate particle showers realistically and manage the related data efficiently. Its aim is to contribute to research on high-energy cosmic ray interactions for two charged particles, like muons and electrons, at very high energies (1016, 1018, and 1019 eV) and taking into account the effect of particles, such as protons, helium nuclei, and iron nuclei, and primary energies and zenith angles (0°, 20°, and 40°). The LDF was also calculated using the Nishimura–Kamata–Greisen function, and good agreement was found with the results produced by the AIRES system for high-energy muons and electrons created by primary particles.

Impact of dark matter scattering on the trajectory of high-energy cosmic rays

Impact of dark matter scattering on the trajectory of high-energy cosmic rays

Cosmic rays, radiocarbon and the beginning of high-energy particle physics: Marietta Blau and the “disintegration stars”

Abstract This brief review describes the discovery of cosmic rays by Victor Hess from the Vienna Radium Institute and the later contribution of Marietta Blau through her observation of “disintegration stars” in photographic emulsion plates exposed to cosmic-ray bombardment. Marietta Blau, a nearly forgotten cosmic-ray pioneer from the Vienna Radium Institute, developed the nuclear emulsion technique for studying nuclear reactions, eventually discovering the disintegration of nuclei through high-energy cosmic rays. Blau survived the Holocaust by escaping to Mexico City from 1939 to 1944. Starting in 1948 at Columbia University, later as a staff member of Brookhaven National Laboratory and then University of Miami, she performed fundamental and original research with nuclear emulsions exposed to 3-GeV protons at the Brookhaven Cosmotron and to 6-GeV protons at the Berkeley Bevatron. Blau returned to Vienna in 1960, where at the Radium Institute a classical β-decay counting facility for radiocarbon dating had been installed, which was finally superseded by the Vienna Environmental Research Accelerator (VERA), a modern versatile accelerator mass spectrometry facility.

Charging and Discharging Modeling of Inertial Sensors Based on Ultraviolet Charge Management

Inertial sensors act as inertial references in space gravitational wave detection missions, necessitating that their internal test mass (TM) maintains minimal residual acceleration noise. High-energy particles and cosmic rays in space, along with ion pumps in ground-based torsion pendulum experiments, can cause charge accumulation on the TM surface, leading to increased residual acceleration noise. Consequently, a charge management system was introduced to control the TM’s charge. The complex light path propagation within the electrode housing (EH) makes the TM’s charging and discharging process difficult to theoretically calculate and fully simulate. To address this issue, we propose a simulation method for charging and discharging inertial sensors within ultraviolet (UV) charge management systems. This method innovatively considers the impact of photoelectron emission angle and the TM’s position offset from symmetry on performance. The method also simulates charging and discharging rates over time under conditions of symmetry and preliminarily examines factors affecting the TM’s equilibrium potential. Simulation results indicate that this method effectively models the charge management system’s operation, providing a valuable reference for system design.

Multi-Messenger Connection in High-Energy Neutrino Astronomy

Low fluxes of astrophysical neutrinos at TeV energies, and the overwhelming background of atmospheric neutrinos below that, render the current paradigm of neutrino astronomy a severely statistics-limited one. While many hints have emerged, all the evidence gathered by IceCube and ANTARES, over the course of almost a decade and a half of operation, has fallen short of providing any conclusive answer to the puzzle of the origin of high-energy cosmic rays and neutrinos. The advancement of the field is thus closely associated with not only the neutrino observatories coming online in the next few years, but also on the coordinated efforts of the EM, GW and cosmic ray communities to develop dedicated channels and infrastructure that allow for the swift and comprehensive multi-messenger follow-up of relevant events detected in any of these sectors. This paper highlights the strides that have been already taken in that direction and the fruits that they have borne, as well as the challenges that lie ahead.

Anomalous events in the TUS orbital detector data

The main goal of the TUS experiment was to search for and study ultra high-energy cosmic rays with energies E > 70 EeV. At the same time, the TUS detector registered a number of unusual events, the origin of which was unclear. Events that are unique and not similar to EAS are the subject of the study presented in this paper. Events such as gamma-ray bursts (GRBs), out-of-aperture upward going EASs accompanied by lightning flashes, as well as terrestrial gamma-ray flashes (TGFs) are considered as their possible sources.

Constraining the hadronic properties of star-forming galaxies above 1 GeV with 15-years Fermi-LAT data

Star-forming and starburst galaxies (SFGs and SBGs) are considered to be powerful emitters of non-thermal γ-rays and neutrinos, due to their intense phases of star-formation activity, which should confine high-energy Cosmic-Rays (CRs) inside their environments. On this regard, the Fermi-LAT collaboration has found a correlation between the γ-ray and infrared luminosities for a sample of local sources. Yet, the physics behind these non-thermal emission is still under debate. We provide novel constraints on the tight relation between γ-rays and star formation rate (SFR) exploiting 15 years of public Fermi-LAT data. Thus, we probe the calorimetric fraction Fcal of high-energy protons in SFGs and SBGs, namely, the fraction of high-energy protons actually producing high-energy γ-rays and neutrinos. Further, we extrapolate this information to their diffuse γ-ray and neutrino emissions constraining their contribution to the extra-galactic gamma-ray background (EGB) and the diffuse neutrino flux. Using the publicly-available fermitools, we analyse 15.3 years of γ-ray between 1-1000 GeV data for 70 sources, 56 of which were not previously detected. We relate this emission to a theoretical model for SBGs in order to constrain Fcal for each source and then study its correlation with the star formation rate of the sources. Firstly, we find at 4σ level an indication of γ-ray emission for other two SBGs, namely M 83 and NGC 1365. By contrast, we find that, even with the new description of background, the significance for the γ-ray emission of M 33 (initially reported as discovered) still stands at ~ 4σ (as already reported by previous works).Along with previous findings, the flux of each detected source is consistent with a ~ E -2.3/2.4 spectrum, compatible with the injected CR flux inferred in the Milky-Way. We also notice that the correlation between Fcal and the SFR is in accordance with the expected scaling relation for CR escape dominated by advection. We remark that undiscovered sources strongly constrain Fcal at 95% CL, providing fundamental information when we interpret the results as common properties of SFGs and SBGs. Finally, we find that these sources might contribute (12 ± 3)% to the EGB, while the corresponding diffuse neutrino flux strongly depends on the spectral index distribution along the source class.

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

ABSTRACT The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

Velocity decorrelation functions of high-energy cosmic rays propagating in magnetic fields

Velocity decorrelation functions of high-energy cosmic rays propagating in magnetic fields

Very-high-energy γ-Ray Emission from Young Massive Star Clusters in the Large Magellanic Cloud

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce very high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy γ-ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the γ-ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The γ-ray luminosity above 0.5 TeV of both sources is (2–3) × 1035 erg s−1. This exceeds by more than a factor of 2 the luminosity of HESS J1646−458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the γ-ray emission from each source is extended with a significance of >3σ and a Gaussian width of about 30 pc. For 30 Dor C, a connection between the γ-ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the γ-ray signal from R136 are discussed.

Anisotropies of diffusive ultra-high energy cosmic rays in [formula omitted] gravity theory

Understanding the anisotropy of ultra high-energy cosmic rays (UHECRs) is crucial for unravelling the origins and propagation mechanisms of these enigmatic particles. In this work, we studied the dipolar anisotropy of UHECRs in the diffusive regime by considering three cosmological models: the standard ΛCDM model, f(R) gravity power-law model and the Starobinsky model. This work aims to see the role of the f(R) gravity theory in understanding the anisotropy of UHECRs without condoning the standard cosmology. We found that the amplitude of the dipolar anisotropy is sensitive to these cosmological models, with the f(R) power-law model predicting the largest amplitude, while the ΛCDM model predicting the smallest amplitude at most of the energies in the range considered. The predicted amplitude of the Starobinsky model lies within the range of the ΛCDM one. This work not only provides a way for exploration of UHECRs anisotropy within different cosmological contexts but also may pave the way for new avenues of research at the intersection of high-energy astrophysics.

The HERO Project (High Energy Cosmic Ray Observatory): Objectives and Design Layout

The HERO Project (High Energy Cosmic Ray Observatory): Objectives and Design Layout

The cosmic-ray spectrum in the PeV to EeV energy range

Cosmic rays around the so-called knee in the spectrum at around PeV primary energy are generally galactic in origin. Observations on the form of their energy spectrum and their mass composition are fundamental tools to understand the origin, acceleration and propagation mechanism of high-energy cosmic rays. In addition, it is required to find signatures to clarify the transition from galactic to extragalactic sources, which are believed to be responsible for the highest-energy cosmic rays above EeV. This brief review focuses on recent experimental results around the knee of the all-particle energy spectrum and composition in the energy range of the knee up to EeV energies.

The Charged Particle Detection System of the SRD Module onboard the China Space Station

The radiation field environment outside the Earth varies greatly with space location and the cycle of solar activity, and the radiation environment near the space station’s orbit is much more complex than that on the Earth’s surface. Various high-energy cosmic ray particles and secondary particles produced by them and the bulkhead of the space station greatly impact the health of astronauts and the working conditions of instruments. Each aircraft is always equipped with some dedicated radiation monitoring instruments to assess the exposure risk for astronauts, and to analyze the causes of instrument failures. The Space Radiation Detector Module (SRDM) is working in the China Space Station (CSS) to measure the radiation environment inside, including two parts: the Charged Particle Detection System (CPDS) and the Neutron Detection System (NDS). The CPDS, which is the main content of this paper, contains a detector unit, that consists of three silicon detectors a BGO calorimeter, and the corresponding readout electronics unit of the detector unit. Ground test results show that the detection system can detect various charged particles from hydrogen to nitrogen ions with an energy resolution of less than 15%. The actual measurement results for a period in orbit show that the main types of charged particles in the cabin are protons and α particles, with measured energies ranging from 0.8 to 265.2 MeV for protons and from 1.0 to 61.6 MeV/A for α particles(where A is the mass number), and linear energy density (LET) range mainly from 1.0 to 612.4 keV/μm. The radiation environment data measured in CSS can provide an important reference value for the exposure risk, life science experiments, and the status of instruments on board.

Radiation tolerance test and performance verification of pnCCD at high temperatures for future satellite mission HiZ-GUNDAM

A wide-field X-ray survey in the soft X-ray band is crucial for future satellite missions in the astronomical observations. HiZ-GUNDAM, currently under development, is a proposed satellite designed to observe soft X-ray transients including gamma-ray bursts. This satellite employs wide-field X-ray monitors consisting of lobster-eye optics and focal-plane pixel sensors in the soft X-ray band of 0.4–4 keV. A pnCCD is a candidate for focal-plane Si pixel detectors, featuring a back-illuminated X-ray CCD, large pixel size (70–100 µm), and a large active image area of approximately 55 × 55 mm2 for the flight model. Here, we investigated the basic characteristics and performance of the small-size pnCCD with 128 × 256 pixels, a pixel size of 75 µm, and a detector thickness of 450 µm. High-energy cosmic rays such as protons can degrade the performance of pnCCDs by increasing dark current and charge transfer inefficiency due to ionizing and displacement damage. These factors may affect soft X-ray observations, potentially causing the degradation of lower-detectable energy thresholds and an increase in the number of hot pixels. Therefore, we conducted a radiation tolerance test at room temperature using a proton beam. After irradiating the pnCCD with 10-MeV protons equivalent to three years of nominal operation for HiZ-GUNDAM, we found that the operation temperatures of the pnCCD should be lower than -35 °C. This requirement will be incorporated into the design of the mission operation system.