Structures in the cosmic ray energy spectra

All the components of Cosmic Rays (CR) have ‘structure’ in their energy spectra at some level, i.e. deviations from a simple power law, and their examination is relevant to the origin of the particles. Emphasis, here, is placed on the large-scale structures in the spectra of nuclei (the ‘knee’ at about 3 PeV), that of electrons plus positrons (a shallow ‘upturn’ at about 100 GeV) and the positron to electron plus positron ratio (an upturn starting at about 5 GeV).Fine structure is defined as deviations from the smooth spectra which already allow for the large-scale structure. Search for the fine structure has been performed in the precise data on positron to electron plus positron ratio measured by the AMS-02 experiment. Although no fine structure is indicated, it could in fact be present at the few percent level.     

Model approximation of cosmic ray spectrum

An analytical model which generalizes the equations describing the intensity of galactic cosmic rays (CR), including both processes, making it applicable in the inner heliosphere (where energy losses dominate) and outer heliosphere (influenced primarily by convection-diffusion processes) is derived. By a suitable choice of a parameter, the proposed model turns into two approximations: solution close to “force-field” model (describing the energy losses of CR in the inner heliosphere) and “convection-diffusion” equation (giving the reduction of CR intensity in the outer heliosphere). A mathematical relation between parameters in the proposed model and the modulation parameter Φ is derived.     

The role of the Galactic Halo and the Single Source in the formation of the cosmic ray anisotropy

The existence of the cosmic ray Halo in our Galaxy has been discussed for more than half a century. If it is real it could help to explain some puzzling features of the cosmic ray flux: its small radial gradient, nearly perfect isotropy and the low level of the fine structure in the energy spectra of the various particles. All these features could be understood if: (a) the Halo has a big size (b) cosmic rays in the Halo have a uniform spatial or radial distribution and (c) the cosmic ray density in the Halo is comparable or even higher than that in the Galactic Disk. The main topic of the paper concerns the present status of the anisotropy and a model for its formation. In our model the extremely small amplitude of the dipole anisotropy is due to the dilution of the anisotropy in the Disk by the dominating isotropic cosmic rays from the Halo. Some minor deviations from complete isotropy in the sub-PeV and PeV energy regions point out to the possible contribution of the Single Source with the phase of its first harmonic opposite to the phase produced by the Disk.     

Heavy quark production in ultra high energy cosmic ray interactions

In this paper we present a comprehensive study of the heavy quark production in ultra high energy cosmic ray interactions in the atmosphere considering that the primary cosmic ray can be either a photon, neutrino or a proton. The analysis is performed using a unified framework – the dipole formalism – and the saturation effects, associated to the physical process of parton recombination, are taken into account. We demonstrate that the contribution of heavy quarks for cosmic ray interactions is in general non-negligible and can be dominant depending of the process considered. Moreover, our results indicate that new dynamical mechanisms should be included in order to obtain reliable predictions for the heavy quark production in pp collisions at ultra high cosmic ray energies.     

From cosmic ray physics to cosmic ray astronomy: Bruno Rossi and the opening of new windows on the universe

Bruno Rossi is considered one of the fathers of modern physics, being also a pioneer in virtually every aspect of what is today called high-energy astrophysics. At the beginning of 1930s he was the pioneer of cosmic ray research in Italy, and, as one of the leading actors in the study of the nature and behavior of the cosmic radiation, he witnessed the birth of particle physics and was one of the main investigators in this fields for many years. While cosmic ray physics moved more and more towards astrophysics, Rossi continued to be one of the inspirers of this line of research. When outer space became a reality, he did not hesitate to leap into this new scientific dimension. Rossi’s intuition on the importance of exploiting new technological windows to look at the universe with new eyes, is a fundamental key to understand the profound unity which guided his scientific research path up to its culminating moments at the beginning of 1960s, when his group at MIT performed the first in situ measurements of the density, speed and direction of the solar wind at the boundary of Earth’s magnetosphere, and when he promoted the search for extra-solar sources of X rays. A visionary idea which eventually led to the breakthrough experiment which discovered Scorpius X-1 in 1962, and inaugurated X-ray astronomy.     

A model of the cosmic ray induced atmospheric neutron environment

In order to optimise the design of space instruments making use of detection materials with low atomic numbers, an understanding of the atmospheric neutron environment and its dependencies on time and position is needed. To produce a simple equation based model, Monte Carlo simulations were performed to obtain the atmospheric neutron fluxes produced by charged galactic cosmic ray interactions with the atmosphere. Based on the simulation results the omnidirectional neutron environment was parametrized including dependencies on altitude, magnetic latitude and solar activity. The upward- and downward-moving component of the atmospheric neutron flux are considered separately. The energy spectra calculated using these equations were found to be in good agreement with data from a purpose built balloon-borne neutron detector, high altitude aircraft data and previously published simulation based spectra.     

A composition dependent energy scale and the determination of the cosmic ray primary mass in the ankle region

At present there are still several open questions about the origin of the ultra high energy cosmic rays. However, great progress in this area has been made in recent years due to the data collected by the present generation of ground based detectors like the Pierre Auger Observatory and Telescope Array. In particular, it is believed that the study of the composition of the cosmic rays as a function of energy can play a fundamental role for the understanding of the origin of the cosmic rays.The observatories belonging to this generation are composed of arrays of surface detectors and fluorescence telescopes. The duty cycle of the fluorescence telescopes is ∼10% in contrast with the ∼100% of the surface detectors. Therefore, the energy calibration of the events observed by the surface detectors is performed by using a calibration curve obtained from a set of high quality events observed in coincidence by both types of detectors. The advantage of this method is that the reconstructed energy of the events observed by the surface detectors becomes almost independent of simulations of the showers because just a small part of the reconstructed energy (the missing energy), obtained from the fluorescence telescopes, comes from simulations. However, the calibration curve obtained in this way depends on the composition of the cosmic rays, which can introduce biases in composition analyses when parameters with a strong dependence on primary energy are considered. In this work we develop an analytical method to study these effects. We consider AMIGA (Auger Muons and Infill for the Ground Array), the low energy extension of the Pierre Auger Observatory corresponding to the surface detectors, to illustrate the use of the method. In particular, we study the biases introduced by an energy calibration dependent on composition on the determination of the mean value of the number of muons, at a given distance to the showers axis, which is one of the parameters most sensitive to primary mass and has an almost linear dependence with primary energy.     

Average anisotropy characteristics of high energy cosmic ray particles and geomagnetic disturbance index Ap

The average characteristics of the diurnal and semi-diurnal anisotropy of cosmic ray intensity at relativistic energies have been obtained by using data from the worldwide grid of neutron monitor for the period 1989 to 1996. The complex behaviour of the diurnal amplitudes and time of maxima (phase) and its association with the Ap index on a long-term and day-to-day basis have been studied. Even though the general characteristics, on a yearly average basis, have not changed significantly during this period, both the diurnal and semi-diurnal amplitudes and phases vary significantly, besides significant changes being observed for different interplanetary conditions on a short-term basis. It is found that the relationship between the Ap index and the diurnal vector is out of phase during the period 1991 to 1995. On a long-term basis, the correlation of diurnal variation with Ap index has been found to vary during the solar cycle. On a short-term basis, it has been observed that the high Ap days are usually associated with higher amplitudes with phase shifted to earlier hours.     

Searches for large-scale anisotropies of cosmic rays: Harmonic analysis and shuffling technique

The measurement of large scale anisotropies in cosmic ray arrival directions is generally performed through harmonic analyses of the right ascension distribution as a function of energy. These measurements are challenging due to the small expected anisotropies and meanwhile the relatively large modulations of observed counting rates due to experimental effects. In this paper, we present a procedure based on the shuffling technique to carry out these measurements, applicable to any cosmic ray detector without any additional corrections for the observed counting rates.     

Measurement of the radial density gradient of cosmic ray in the heliosphere by the GRAPES-3 experiment

A radial anisotropy in the flux of cosmic rays in heliosphere was theoretically predicted by Parker and others within the framework of the diffusion–convection mechanism. The solar wind is responsible for sweeping out the galactic cosmic rays, creating a radial density gradient within the heliosphere. This gradient coupled with the interplanetary magnetic field induces a flow of charged particles perpendicular to the ecliptic plane which was measured and correctly explained by Swinson, and is hereafter referred as ‘Swinson flow’. The large area GRAPES-3 tracking muon telescope offers a powerful probe to measure the Swinson flow and the underlying radial density gradient of the galactic cosmic rays at a relatively high rigidity of ∼100 GV. The GRAPES-3 data collected over a period of six years (2000–2005) were analyzed and the amplitude of the Swinson flow was estimated to be (0.0644 ± 0.0008)% of cosmic ray flux which was an ∼80σ effect. The phase of the maximum flow was at a sidereal time of (17.70 ± 0.05) h which was 18 min earlier than the expected value of 18 h. This small 18 min phase difference had a significance of ∼6σ indicating the inherent precision of the GRAPES-3 measurement. The radial density gradient of the galactic cosmic rays at a median rigidity of 77 GV was found to be 0.65% AU⁻¹.     

Prediction of the diffuse neutrino flux from cosmic ray interactions near supernova remnants

In this paper, we present high-energy neutrino spectra from 21 Galactic supernova remnants (SNRs), derived from gamma-ray measurements in the GeV–TeV range. We find that only the strongest sources, i.e. G40.5-0.5 in the north and Vela Junior in the south could be detected as single point sources by IceCube or KM3NeT, respectively. For the first time, it is also possible to derive a diffuse signal by applying the observed correlation between gamma-ray emission and radio signal. Radio data from 234 supernova remnants listed in Green’s catalog are used to show that the total diffuse neutrino flux is approximately a factor of 2.5 higher compared to the sources that are resolved so far. We show that the signal at above 10 TeV energies can actually become comparable to the diffuse neutrino flux component from interactions in the interstellar medium. Recently, the IceCube collaboration announced the detection of a first diffuse signal of astrophysical high-energy neutrinos. Directional information cannot unambiguously reveal the nature of the sources at this point due to low statistics. A number of events come from close to the Galactic center and one of the main questions is whether at least a part of the signal can be of Galactic nature. In this paper, we show that the diffuse flux from well-resolved SNRs is at least a factor of 20 below the observed flux.     

Search for point-like sources of cosmic rays with energies above 10 eV in the HiRes-I monocular data set

We report the results of a search for point-like deviations from isotropy in the arrival directions of ultra-high energy cosmic rays in the northern hemisphere. In the monocular data set collected by the High-Resolution Fly’s Eye, consisting of 1525 events with energy exceeding 10^18.5 eV, we find no evidence for point-like excesses. We place a 90% c.l. upper limit of 0.8 hadronic cosmic rays/km² yr on the flux from such sources for the northern hemisphere and place tighter limits as a function of position in the sky.     

Radio Cherenkov signals from the Moon: Neutrinos and cosmic rays

Neutrino production of radio Cherenkov signals in the Moon is the object of radio telescope observations. Depending on the energy range and detection parameters, the dominant contribution to the neutrino signal may come from interactions of the neutrino on the Moon facing the telescope, rather than neutrinos that have traversed a portion of the Moon. Using the approximate analytic expression of the effective lunar aperture from a recent paper by Gayley, Mutel and Jaeger, we evaluate the background from cosmic ray interactions in the lunar regolith. We also consider the modifications to the effective lunar aperture from generic non-standard model neutrino interactions. A background to neutrino signals are radio Cherenkov signals from cosmic ray interactions. For cosmogenic neutrino fluxes, neutrino signals will be difficult to observe because of low neutrino flux at the high energy end and large cosmic ray background in the lower energy range considered here. We show that lunar radio detection of neutrino interactions is best suited to constrain or measure neutrinos from astrophysical sources and probe non-standard neutrino-nucleon interactions such as microscopic black hole production.     

Numerical studies on the structure of the cosmic ray electron halo in starburst galaxies

The structure of the cosmic ray electron halo of a starburst galaxy depends strongly on the nature of galactic wind and the configuration of the magnetic field. We have investigated these dependencies by solving numerically the propagation of electrons originating in starburst galaxies, most likely in supernova remnants. The calculations are made for several models for the galactic winds and for the configuration of the magnetic fields for comparison with observations. Our simulation of a quasi-radio halo reproduces both the extended structure of ∼9 kpc and the subtle hollow structure near the polar region of the radio halo that are observed in the starburst galaxy NGC 253. These findings suggest the existence of strong galactic wind in NGC 253.      

Atmospheric ionization and cosmic rays: studies and measurements before 1912

The discovery of cosmic rays, a milestone in science, was based on the work by scientists in Europe and the New World and took place during a period characterized by nationalism and lack of communication. Many scientists that took part in this research a century ago were intrigued by the penetrating radiation and tried to understand the origin of it. Several important contributions to the discovery of the origin of cosmic rays have been forgotten; historical, political and personal facts might have contributed to their substantial disappearance from the history of science.     

The cosmic ray primary composition in the “knee” region through the EAS electromagnetic and muon measurements at EAS-TOP

The evolution of the cosmic ray primary composition in the energy range 10⁶–10⁷ GeV (i.e. the “knee” region) is studied by means of the e.m. and muon data of the Extensive Air Shower EAS-TOP array (Campo Imperatore, National Gran Sasso Laboratories). The measurement is performed through: (a) the correlated muon number (N_μ) and shower size (N_e) spectra, and (b) the evolution of the average muon numbers and their distributions as a function of the shower size. From analysis (a) the dominance of helium primaries at the knee, and therefore the possibility that the knee itself is due to a break in their energy spectrum (at E_k^He=(3.5±0.3)×10⁶ GeV) are deduced. Concerning analysis (b), the measurement accuracies allow the classification in terms of three mass groups: light (p,He), intermediate (CNO), and heavy (Fe). At primary energies E₀≈10⁶ GeV the results are consistent with the extrapolations of the data from direct experiments. In the knee region the obtained evolution of the energy spectra leads to: (i) an average steep spectrum of the light mass group (γ_p,He>3.1), (ii) a spectrum of the intermediate mass group harder than the one of the light component (γ_CNO2.75, possibly bending at E_k^CNO≈(6–7)×10⁶ GeV), (iii) a constant slope for the spectrum of the heavy primaries (γ_Fe2.3–2.7) consistent with the direct measurements. In the investigated energy range, the average primary mass increases from lnA=1.6–1.9 at E₀1.5×10⁶ GeV to lnA=2.8–3.1 at E₀1.5×10⁷ GeV. The result supports the standard acceleration and propagation models of galactic cosmic rays that predict rigidity dependent cut-offs for the primary spectra of the different nuclei. The uncertainties connected to the hadronic interaction model (QGSJET in CORSIKA) used for the interpretation are discussed.     

The time-dependent transport of cosmic rays in the heliosphere

The modulation of cosmic rays (CRs) in the heliosphere is a dynamic and therefore a highly time-dependent process. Numerical models with only a time-dependent neutral sheet prove to be successful when moderate to low solar activity occurs but fail to describe large and discrete steps in modulated CRs when solar activity is high. To explain this feature of heliospheric modulation, the concept of global merged interaction regions (GMIRs) is required. The combination of gradient, curvature and neutral sheet drifts with these GMIRs has so far been the most successful approach in explaining the 11-year and 22-year cycles in the long-term modulation of CRs.     

Contribution of nuclei accelerated by gamma-ray pulsars to cosmic rays in the Galaxy

We consider the galactic population of gamma-ray pulsars as possible sources of cosmic rays at and just above the “knee” in the observed cosmic ray spectrum at 10¹⁵–10¹⁶ eV. We suggest that iron nuclei may be accelerated in the outer gaps of pulsars, and then suffer partial photo-disintegration in the non-thermal radiation fields of the outer gaps. As a result, protons, neutrons, and surviving heavier nuclei are injected into the expanding supernova remnant. We compute the spectra of nuclei escaping from supernova remnants into the interstellar medium, taking into account the observed population of radio pulsars.

Our calculations, which include a realistic model for acceleration and propagation of nuclei in pulsar magnetospheres and supernova remnants, predict that heavy nuclei accelerated directly by gamma-ray pulsars could contribute about 20% of the observed cosmic rays in the knee region. Such a contribution of heavy nuclei to the cosmic ray spectrum at the knee can significantly increase the average value of lnA with increasing energy as is suggested by recent observations.