Cosmic rays are often modeled as charged particles. This allows their non-ballistic propagation in magnetized structures to be captured. In certain situations, a neutral cosmic ray component can arise. For example, cosmic ray neutrons are produced in considerable numbers through hadronic pp and p$γ$ interactions. At ultrahigh energies, the decay timescales of these neutrons is dilated, allowing them to traverse distances on the scale of galactic and cosmological structures. Unlike charged cosmic rays, neutrons are not deflected by magnetic fields. They propagate ballistically at the speed of light in straight lines. The presence of a neutral baryonic cosmic ray component formed in galaxies, clusters and cosmological filaments can facilitate the escape and leakage of cosmic rays from magnetic structures that would otherwise confine them. We show that, by allowing confinement breaking, the formation of cosmic-ray neutrons by high-energy hadronic interactions in large scale astrophysical structures can modify the exchange of ultra high-energy particles across magnetic interfaces between galaxies, clusters, cosmological filaments and voids.
The IceCube Observatory comprises a cubic-kilometer particle detector deep in the Antarctic ice and the cosmic-ray air-shower array IceTop at the surface above. Previous analyses of the cosmic-ray composition have used coincident events with IceTop detecting the electromagnetic shower footprint as well as GeV muons, while the sensors submerged in the ice measure the TeV muons from the same events. The energy range of previous composition analyses, however, has been limited to 3 PeV primary energy and above, whereas the IceTop all-particle energy spectrum has been extended down to 250 TeV. This contribution presents a method to reconstruct the combined spectrum of cosmic-ray protons and helium nuclei, starting at 200 TeV primary energy. The resulting H+He spectrum closes the gap in the measurements of light cosmic rays between IceCube as well as KASCADE and experiments measuring in the TeV energy range, such as DAMPE and HAWC.
This study presents preliminary results from the analysis of cosmic-ray anisotropy using air showers detected by the IceTop surface array between 2011 and 2022. With improved statistical precision and updated Monte Carlo simulation events compared to previous IceTop reports, we investigate anisotropy patterns across four energy ranges spanning from 300 TeV to 6.9 PeV. This work extends the measurement of cosmic-ray anisotropy in the southern hemisphere to higher energies than previously achieved with IceTop. Our results provide a foundation for exploring potential connections between the observed anisotropy, the energy spectrum, and the mass composition of the cosmic-ray flux.
This article is the rapporteur's summary of the cosmic ray indirect sessions of the 38th International Cosmic Ray Conference in Nagoya, Japan. The rapporteur highlights cosmic ray indirect observatories around the world, and reviews a selection of the latest results regarding the cosmic ray energy spectrum, mass composition, anisotropy, hadronic interaction models, theory, geophysics, interdisciplinary research, and future projects.
The hypothesis that the entire cosmic ray spectrum, from $\lesssim1\,{\rm GeV}$ to $\gtrsim100\,{\rm EeV}$ energy, can be accounted for by diffusive shock acceleration on increasingly large scales is critically examined. Specifically, it is conjectured that Galactic cosmic rays, up to $\sim3\,{\rm PeV}$, are mostly produced by local supernova remnants, from which they escape upstream. These cosmic rays initiate a powerful magnetocentrifugal wind, removing disk mass and angular momentum before passing through the Galactic Wind Termination Shock at a radius $\sim200\,{\rm kpc}$, where they can be re-accelerated to account for observed cosmic rays up to $\sim30\,{\rm PeV}$. The cosmic rays transmitted downstream from more powerful termination shocks associated with other galaxies can be further accelerated at Intergalactic Accretion Shocks to the highest energies observed. In this interpretation, the highest rigidity observed particles are protons; the highest energy particles are heavy nuclei, such as iron. A universal "bootstrap" prescription, coupling the energy density of the magnetic turbulence to that of the resonant cosmic rays, is proposed, initially for the highest energy par
HAWC is a ground-based observatory consisting of 300 water Cherenkov detectors, which observes the extensive air showers induced by cosmic rays from some TeV to a few PeV and, in particular, gamma rays from 300 GeV to more than 100 TeV. One of the crucial features required for a detector of extensive air showers is the estimation of the primary energy of the events to study the spectra of cosmic and gamma rays. For HAWC there are currently two gamma-ray energy estimators: one relies on a ground density parameter, while the other utilizes an artificial neural network. For the cosmic ray energy estimation, there is only one estimator based on maximum likelihood procedures and measurements of the lateral charge distribution of the events. It is worthwhile to update the cosmic-ray energy estimator due to recent improvements of the extensive air shower offline-reconstruction techniques in HAWC. Therefore, we implemented an artificial neural network to reconstruct the primary energy of hadronic events trained with several observables that characterize the air showers. We trained several models and evaluated their performance against the existing cosmic ray energy estimator. In this work,
Cosmic rays are the only agent able to penetrate into the interior of dense molecular clouds. Depositing (part of) their energy through ionisation, cosmic rays play an essential role in determining the physical and chemical evolution of star-forming regions. To a first approximation their effect can be quantified by the cosmic-ray induced ionization rate. Interestingly, theoretical estimates of the ionization rate assuming the cosmic-ray spectra observed in the local interstellar medium result in an ionization rate that is one to two orders of magnitude below the values inferred from observations. However, due to the discrete nature of sources, the local spectra of MeV cosmic rays are in general not representative for the spectra elsewhere in the Galaxy. Such stochasticity effects have the potential of reconciling modelled ionization rates with measured ones. Here, we model the distribution of low-energy cosmic-ray spectra expected from a statistical population of supernova remnants in the Milky Way. The corresponding distribution for the ionization rate is derived and confronted with data. We find that the stochastic uncertainty helps with explaining the surprisingly high ionizati
Anisotropic diffusion is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and the interplay of cosmic rays with the Galactic magnetic field. This diffusion term contributes to the highly stiff nature of the cosmic ray transport equation. To conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators (namely, Crank-Nicolson (CN)) have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods to treat the linear anisotropc diffusion equation in the presence of advection and time-independent and time-dependent sources. These methods allow us to take even larger step sizes that can substantially speed-up the simulations whilst generating highly accurate solutions. In or subsequent work, we will use these exponential solvers in the Picard code to study anisotropic cosmic ray diffusion and we will consider additional physical processes such as continuous momentum losses and reacceleration.
The Global Cosmic-ray Observatory (GCOS) is a proposed large-scale observatory for studying ultra-high-energy cosmic particles, including ultra-high-energy cosmic rays (UHECRs), photons, and neutrinos. Its primary goal is to characterise the properties of the highest-energy particles in Nature with unprecedented accuracy, and to identify their elusive sources. With an aperture at least a ten-fold larger than existing observatories, this next-generation facility should start operating after 2030, when present-day detectors will gradually cease their activities. Here we briefly review the scientific case motivating GCOS. We present the status of the project, preliminary ideas for its design, and some estimates of its capabilities.
Recent balloon-borne and satellite experiments have established new features in the behavior of the spectra of cosmic rays. An analysis of all the data showed that hardening of most abundant primary cosmic ray nuclei spectra with increasing rigidity is observed in $>200$ GV region. At the same time, rigidity dependencies of secondary cosmic rays are distinctly different. The AMS-02 data show that above 200 GV the secondary cosmic rays Li, Be, B harden more than the primary He, C, O. In this paper we discuss a new scenario that self-consistently describes these new features of the cosmic ray spectra. We demonstrate that the measured by the AMS-02 changes in the slope of the cosmic ray spectra caused by the transition from the contribution of multiple distant Galactic sources, including the nuclear interactions of the particles accelerated by these sources with an interstellar medium during their wandering in the Galaxy, to the contribution of mainly local ones. We also found that the spectral universality is observed at rigidity $R> 10^5$ GV.
The description of the transport of cosmic rays in magnetized media is central to both acceleration and propagation of these particles in our Galaxy and outside. The investigation of the process of particle acceleration, especially at shock waves, has already emphasized that non-linear effects such as self-generation of waves and dynamical reaction of cosmic rays on the background plasmas, are crucial if to achieve a physical understanding of the origin of cosmic rays. Here we discuss how similar non-linear effects on Galactic scales may affect the propagation of cosmic rays, not only through the excitation of plasma waves important for particle scattering, but also by inducing the motion of the interstellar medium in the direction opposite to the gravitational pull exerted by matter in the Galaxy, thereby resulting in the launching of a wind. The recent discovery of several unexpected features in cosmic ray spectra (discrepant hardening, spectral breaks in the H and He spectra, rising positron fraction and unexpectedly hard antiproton spectrum) raises the question of whether at least some of these effects may be attributed to poorly understood aspects of cosmic ray transport.
Plasma outflow from a gravitational potential well with cosmic rays and self-excited Alfvén waves with cooling and wave damping is studied in the hydrodynamics regime. We study outflows in the presence of cosmic ray and Alfvén waves including the effect of cooling and wave damping. We seek physically allowable steady-state subsonic-supersonic transonic solutions. We adopted a multi-fluid hydrodynamical model for the cosmic ray plasma system. Thermal plasma, cosmic rays, and self-excited Alfvén waves are treated as fluids. Interactions such as cosmic-ray streaming instability, cooling, and wave damping were fully taken into account. We considered one-dimensional geometry and explored steady-state solutions. The model is reduced to a set of ordinary differential equations, which we solved for subsonic-supersonic transonic solutions with given boundary conditions at the base of the gravitational potential well. We find that physically allowable subsonic-supersonic transonic solutions exist for a wide range of parameters. We studied the three-fluid system (considering only forward-propagating Alfvén waves) in detail. We examined the cases with and without cosmic ray diffusion separatel
In the standard picture of cosmic ray transport the propagation of charged cosmic rays through turbulent magnetic fields is described as a random walk with cosmic rays scattering on magnetic field turbulence. This is in good agreement with the highly isotropic cosmic ray arrival directions as this diffusion process effectively isotropizes the cosmic ray distribution. High-statistics observatories like IceCube and HAWC have however observed significant deviations from isotropy down to very small angular scales. This is in strong tension with this standard picture of cosmic ray propagation. While large scale multipoles arise naturally, for example due to the earth's motion relative to the isotropic cosmic ray distribution, there is no intuitive mechanism to account for the observed anisotropies at smaller angular scales. By relaxing one of the standard assumptions of quasi linear theory and treating correlations between fluxes of cosmic rays from different directions explicitly we show that higher multipoles also are to be expected from particle propagation through turbulent magnetic fields. We present a first analytical calculation of the angular power spectrum assuming a physically
Cosmic rays are a fundamental source of ionization for molecular and diffuse clouds, influencing their chemical, thermal, and dynamical evolution. The amount of cosmic rays inside a cloud also determines the $γ$-ray flux produced by hadronic collisions between cosmic rays and cloud material. We study the spectrum of cosmic rays inside and outside of a diffuse cloud, by solving the stationary transport equation for cosmic rays including diffusion, advection and energy losses due to ionization of neutral hydrogen atoms. We found that the cosmic ray spectrum inside a diffuse cloud differs from the one in the interstellar medium (ISM) for energies smaller than $E_{br}\approx 100$ MeV, irrespective of the model details. Below $E_{br}$, the spectrum is harder (softer) than that in the ISM if the latter is a power law $\propto p^{-s}$ with $s$ larger (smaller) than $\sim0.42$. As a consequence also the ionization rate due to CRs is strongly affected. Assuming an average Galactic spectrum similar to the one inferred from AMS-2 and Voyager 1 data, we discuss the resulting ionization rate in a typical diffuse cloud.
Our purpose is to evaluate the rate of the maximum energy and the acceleration rate that cosmic rays acquire in the non-relativistic diffusive shock acceleration as it could apply during their lifetime in various astrophysical sites, where highly oblique shocks exist. We examine numerically (using Monte Carlo simulations) the effect of the diffusion coefficients on the energy gain and the acceleration rate, by testing the role between the obliquity of the magnetic field to the shock normal, and the significance of both perpendicular cross-field diffusion and parallel diffusion coefficients to the acceleration rate. We find (and justify previous analytical work - Jokipii 1987) that in highly oblique shocks the smaller the perpendicular diffusion gets compared to the parallel diffusion coefficient values, the greater the energy gain of the cosmic rays to be obtained. An explanation of the cosmic ray spectrum in high energies, between $10^{15}$eV and about $10^{18}$eV is claimed, as we estimate the upper limit of energy that cosmic rays could gain in plausible astrophysical regimes; interpreted by the scenario of cosmic rays which are injected by three different kind of sources, (a) s
We discuss the influence of large scale cosmic magnetic fields on the propagation of hadronic cosmic rays above 10^19 eV based on large scale structure simulations. Our simulations suggest that rather substantial deflection up to several tens of degrees at 10^20 eV are possible for nucleon primaries. Further, spectra and composition of cosmic rays from individual sources can depend on magnetic fields surrounding these sources in intrinsically unpredictable ways. This is true even if deflection from such individual sources is small. We conclude that the influence of large scale cosmic magnetic fields on ultra-high energy cosmic ray propagation is currently hard to quantify. We discuss possible reasons for discrepant results of simulations by Dolag et al. which predict deflections of at most a few degrees for nucleons. We finally point out that even in these latter simulations a possible heavy component would in general suffer substantial deflection.
Diffusive shock acceleration is the prime candidate for efficient acceleration of cosmic rays. Galactic cosmic rays are believed to originate predominantly from this process in supernova remnant shock waves. Confinement of the cosmic rays in the shock region is key in making the mechanism effective. It has been known that on small scales (smaller than the typical gyroradius) high-amplitude non-resonant instabilities arise due to cosmic ray streaming ahead of the shock. For the efficiency of scattering of the highest energy cosmic rays it is of interest to determine the type of instabilities that act on longer length scales, i.e. larger than the cosmic ray gyroradius. We will present the results of our analysis of an instability that acts in this regime and will discuss its driving mechanism and typical growth times.
The High Altitude Water Cherenkov (HAWC) observatory is a ground-based air-shower detector designed to study the TeV gamma and cosmic ray windows. The observatory is composed of a densely packed array of $300$ water Cherenkov tanks, $4.5$ m deep and $7.3$ diameter with $4$ photomultipliers (PMT) each, distributed on a $22,000 \, \mbox{m}^2$ surface. The instrument registers the number of hit PMT's as well as the timing information and the total charge seen by the photomultipliers during the shower event. From the analysis of these data, shower observables such as the arrival direction, the core position at ground, the age and the primary energy are estimated, from which the energy spectrum of cosmic rays and its composition can be studied. In this work, we will describe the methodologies of HAWC cosmic ray analyses, including the Bayesian spectral unfolding procedure used to determine the all-particle and the light component energy spectra of cosmic rays. We will see that the distribution of the shower age vs the fraction of hit PMT's contains information about the composition of cosmic rays.
Cosmic rays and related measurements: IceTop overview; IceTop calibration; Cosmic ray composition; Air shower measurements; Studies of cosmic ray muons; Studies of PeV gamma rays; Cosmic ray anisotropy; Forbush decrease; Submitted papers to the 32nd International Cosmic Ray Conference, Beijing 2011.
The theory of Galactic Winds, driven by the cosmic-ray pressure gradient, is reviewed both on the magnetohydrodynamic and on the kinetic level. In this picture the magnetic field of the Galaxy above the dense gas disk is assumed to have a flux tube geometry, the flux tubes rising locally perpendicular out of the disk to become radially directed at large distances, with the cosmic-ray sources located deep within the Galactic disk. At least above the gas disk, the magnetic fluctuations which resonantly scatter the cosmic rays are selfconsistently excited as Alf{`e}n waves by the escaping cosmic rays. The fluctuation amplitudes remain finite through nonlinear wave dissipation. The spatially increasing speed of the resulting outflow results in a diffusion-convection boundary whose position depends on particle momentum. It replaces the escape boundary of static diffusion models. New effects like overall Galactic mass and angular momentum loss as well as gas heating beyond the disk appear. Also particle re-acceleration in the distant wind halo suggests itself. The resulting magnetohydrodynamic flow properties and the cosmic-ray transport properties are compared with observations. On the