Explore the vast range of titles available.

GRAPES-3 is a mid-altitude (2200 m) and near-equatorial (11.°4N) air shower array, overlapping in its field of view for cosmic-ray observations with experiments that are located in the Northern and Southern Hemispheres. We analyze a sample of 3.7 × 109 cosmic-ray events collected by the GRAPES-3 experiment between 2013 January 1 and 2016 December 31 with a median energy of ∼16 TeV for study of small-scale (<60°) angular-scale anisotropies. We observed two structures, labeled A and B, that deviate from the expected isotropic distribution of cosmic rays in a statistically significant manner. Structure A spans 50°–80° in R.A. and from −15° to 30° in decl. The relative excess observed in structure A is at the level of (6.5 ± 1.3) × 10−4 with a statistical significance of 6.8 standard deviations. Structure B is observed in the R.A. range 110°–140° and at decl. from −10° to 30°. The relative excess observed in this region is at the level of (4.9 ± 1.4) × 10−4 with a statistical significance of 4.7 standard deviations. These structures are consistent with those reported by Milagro, ARGO-YBJ, and HAWC. These observations could provide a better understanding of the sources of cosmic rays, their propagation, and the magnetic structures in our Galaxy.

The past decade has brought impressive advances in the astrophysics of cosmic rays (CRs) and multiwavelength astronomy, thanks to the new instrumentation launched into space and built on the ground. Modern technologies employed by those instruments provide measurements with unmatched precision, enabling searches for subtle signatures of dark matter and new physics. Understanding the astrophysical backgrounds to better precision than the observed data is vital in moving to this new territory. A state-of-the-art CR propagation code, called GalProp, is designed to address exactly this challenge. Having 25 yr of development behind it, the GalProp framework has become a de facto standard in the astrophysics of CRs, diffuse photon emissions (radio to γ-rays), and searches for new physics. GalProp uses information from astronomy, particle physics, and nuclear physics to predict CRs and their associated emissions self-consistently, providing a unifying modeling framework. The range of its physical validity covers 18 orders of magnitude in energy, from sub-keV to PeV energies for particles and from μeV to PeV energies for photons. The framework and the data sets are public and are extensively used by many experimental collaborations and by thousands of individual researchers worldwide for interpretation of their data and for making predictions. This paper details the latest release of the GalProp framework and updated cross sections, further developments of its initially auxiliary data sets for models of the interstellar medium that grew into independent studies of the Galactic structure—distributions of gas, dust, radiation, and magnetic fields—as well as the extension of its modeling capabilities. Example applications included with the distribution illustrating usage of the new features are also described.

We present a suite of models of the coherent magnetic field of the Galaxy based on new divergence-free parametric functions describing the global structure of the field. The model parameters are fit to the latest full-sky Faraday rotation measures (RMs) of extragalactic sources and polarized synchrotron intensity (PI) maps from the Wilkinson Microwave Anisotropy Probe and Planck. We employ multiple models for the density of thermal and cosmic-ray electrons in the Galaxy, needed to predict the sky maps of RMs and PI for a given Galactic magnetic field (GMF) model. The robustness of the inferred properties of the GMF is gauged by studying many combinations of parametric field models and electron density models. We determine the pitch angle of the local magnetic field (11° ± 1°), explore the evidence for a grand-design spiral coherent magnetic field (inconclusive), determine the strength of the toroidal and poloidal magnetic halo fields below and above the disk (magnitudes the same for both hemispheres within ≈10%), set constraints on the half-height of the cosmic-ray diffusion volume (≥2.9 kpc), investigate the compatibility of RM- and PI-derived magnetic field strengths (compatible under certain assumptions), and check if the toroidal halo field could be created by the shear of the poloidal halo field due to the differential rotation of the Galaxy (possibly). A set of eight models is identified to help quantify the present uncertainties in the coherent GMF spanning different functional forms, data products, and auxiliary input. We present the corresponding sky maps of rates for axion–photon conversion in the Galaxy and deflections of ultrahigh-energy cosmic rays.

Neutron monitors (NMs) are large ground-based detectors of atmospheric secondary particles, mostly neutrons, from primary cosmic rays. Their sky direction and rigidity imply a well-defined incoming (asymptotic) direction in space. From 2015 December 16 to 2017 January 8, 6 of the 18 NM counters had been transferred from McMurdo to Jang Bogo, both in Antarctica, so data from similar detectors were recorded simultaneously at these two nearby NM stations. Autocorrelations of these NM count rates are well fit as the sum of three components: an exponential function and a cosine with a period of 1 day, both centered at zero lag, plus a constant. Fitting the cross correlation of the two count rates, the functions are no longer centered at zero lag. The best-fit cosine phase is at time lag −160.22 ± 0.12 minutes. Calculating cosmic-ray trajectories in Earth's magnetic field throughout the time interval, the mean difference in response-weighted asymptotic longitudes corresponds to time lag −169.41 ± 0.31 minutes, in close agreement with the observed lag. Thus, the cosine term is consistent with and provides a technique to cleanly measure the cosmic-ray anisotropy. In contrast, the peak term shows a time lag of –14.55 minutes, much closer to the –9.60 minutes lag in rotation due to the difference in geographic longitude. We find a similar behavior in the correlations between other pairs of stations. We propose that rapid fluctuations in the counting rate may be primarily due to cosmic-ray particles of very high energy.

Long-term observations of the Galactic center by Fermi and HESS have revealed a novel phenomenon: the high-energy gamma-ray spectrum from the gamma-ray source HESS J1745-290 exhibits a double power-law structure. In this study, we propose a new explanation for this phenomenon. We suggest that the low-energy (GeV) power-law spectrum originates from interactions between trapped background “sea” cosmic ray particles and the dense gaseous environment near the Galactic center. In contrast, the bubble-like structure in the high-energy (TeV) spectrum is produced by protons accelerated during active phases of the Galactic center, through the same physical process. Based on this framework, we first calculate the gamma-ray emission generated by cosmic ray protons accelerated in the Galactic center. Then, using a spatially dependent cosmic ray propagation model, we compute the energy spectrum of background “sea” cosmic ray protons and their associated diffuse gamma-ray emission in the Galactic center region. The results closely reproduce the observations from Fermi-LAT and HESS, suggesting that their long-term data support this picture: high-energy cosmic rays in the local region originate from nearby cosmic ray sources, while low-energy cosmic rays are a unified contribution from distant cosmic ray sources. We anticipate that this double power-law structure may be widely present in the halo of a Galactic cosmic-ray source or a slow-diffusion region. We hope that future observations will detect more such sources, allowing us to further test and validate our model.

Supernova remnants (SNRs) have long been suspected to be the primary sources of Galactic cosmic rays. Over the past decades, great strides have been made in the modeling of particle acceleration, magnetic field amplification, and escape from SNRs. Yet while many SNRs have been observed in nonthermal emission in radio, X-rays, and gamma rays, there is no evidence for any individual object contributing to the locally observed flux. Here, we propose a particular spectral signature from individual remnants that is due to the energy-dependent escape from SNRs. For young and nearby sources, we predict fluxes enhanced by tens of percent in narrow rigidity intervals; given the percent-level flux uncertainties of contemporary cosmic-ray data, such features should be readily detectable. We model the spatial and temporal distribution of sources and the resulting distribution of fluxes with a Monte Carlo approach. The decision tree that we have trained on simulated data is able to discriminate with very high significance between the null hypothesis of a smooth distribution of sources and the scenario with a stochastic distribution of individual sources. We suggest that this cosmic-ray energy-dependent injection time (CREDIT) scenario be considered in experimental searches to identify individual SNRs as cosmic-ray sources.

The detection of cosmic neutrinos with energies above a teraelectronvolt (TeV) offers a unique exploration into astrophysical phenomena 1 , 2 – 3 . Electrically neutral and interacting only by means of the weak interaction, neutrinos are not deflected by magnetic fields and are rarely absorbed by interstellar matter: their direction indicates that their cosmic origin might be from the farthest reaches of the Universe. High-energy neutrinos can be produced when ultra-relativistic cosmic-ray protons or nuclei interact with other matter or photons, and their observation could be a signature of these processes. Here we report an exceptionally high-energy event observed by KM3NeT, the deep-sea neutrino telescope in the Mediterranean Sea 4 , which we associate with a cosmic neutrino detection. We detect a muon with an estimated energy of 12 0 − 60 + 110 petaelectronvolts (PeV). In light of its enormous energy and near-horizontal direction, the muon most probably originated from the interaction of a neutrino of even higher energy in the vicinity of the detector. The cosmic neutrino energy spectrum measured up to now 5 , 6 – 7 falls steeply with energy. However, the energy of this event is much larger than that of any neutrino detected so far. This suggests that the neutrino may have originated in a different cosmic accelerator than the lower-energy neutrinos, or this may be the first detection of a cosmogenic neutrino 8 , resulting from the interactions of ultra-high-energy cosmic rays with background photons in the Universe. A very high-energy muon observed by the KM3NeT experiment in the Mediterranean Sea is evidence for the interaction of an exceptionally high-energy neutrino of cosmic origin.  

Recent studies have shown that the anisotropy is of great value to decipher cosmic rays’ origin and propagation. We have built a unified scenario to describe the observations of the energy spectra and the large-scale anisotropy and called attention to their synchronous evolution with energy. In this work, the impact of the local regular magnetic field (LRMF) and corresponding anisotropic diffusion on large-scale anisotropy have been investigated. When the perpendicular diffusion coefficient is much smaller than the parallel one, the dipole anisotropy points to the LRMF and the observational phase below 100 TeV could be reproduced. Moreover, we find that the dipole phase above 100 TeV strongly depends on the evolution of local diffusion. But the current measurements at that energy are still scarce. We suggest that more precise measurements at that energy could be carried out to unveil the local diffusion and further the local turbulence.

Neutron monitors (NMs) are ground-based devices designed to measure cosmic-ray count rates by monitoring atmospheric neutrons from cosmic-ray showers. We present results from new electronics that have recorded cross-counter time delay histograms for the Princess Sirindhorn Neutron Monitor (PSNM) at the summit of Doi Inthanon, Thailand. From these histograms, we have extracted the cross-counter leader fraction (L) and corrected it for atmospheric effects. For large counter separation, we measure nearly constant L ≈ 0.997, implying that 0.3% of counts in one counter are temporally associated with later counts on a given distant counter. Monte Carlo simulations confirm that individual secondary particles cannot account for the associated counts at large counter separation, which instead requires a contribution from multiple secondary particles in the same cosmic-ray shower that is apparently independent of distance over 3 to 7.5 m. We infer that ≈4.5% of PSNM counts are associated with a later count in at least one of its 18 counters from a different secondary particle in the same shower. Monte Carlo simulations of atmospheric showers and NM yield functions can be validated using our measurements of neutron multiplicity across counters and the contributions of single and multiple secondary particles. These measurements also improve understanding of the single-counter L, which has been used for precise tracking of cosmic-ray spectral variations and extending the range of NM observations to higher energies.

A dipole anisotropy in ultra–high-energy cosmic ray (UHECR) arrival directions, of extragalactic origin, is now firmly established at energies E > 8 EeV. Furthermore, the UHECR angular power spectrum shows no power at smaller angular scales than the dipole, apart from hints of possible individual hot or warm spots for energy thresholds ≳40 EeV. Here we exploit the magnitude of the dipole and the limits on smaller-scale anisotropies to place constraints on two quantities: the extragalactic magnetic field (EGMF) and the number density of UHECR sources or the volumetric event rate if UHECR sources are transient. We also vary the bias between the extragalactic matter and the UHECR source densities, reflecting whether UHECR sources are preferentially found in over- or underdense regions, and find that little or no bias is favored. We follow Ding et al. (2021) in using the CosmicFlows-2 density distribution of the local universe as our baseline distribution of UHECR sources, but we improve and extend that work by employing an accurate and self-consistent treatment of interactions and energy losses during propagation. Deflections in the Galactic magnetic field are treated using either the full JF12 magnetic field model, with both random and coherent components, or just the coherent part, to bracket the impact of the GMF on the dipole anisotropy. This large-scale structure model gives good agreement with both the direction and magnitude of the measured dipole anisotropy and forms the basis for simulations of discrete sources and the inclusion of EGMF effects.