Measurement of the cosmic ray composition at the knee with the SPASE-2/AMANDA-B10 detectors

The mass composition of high-energy cosmic rays at energies above 10¹⁵ eV can provide crucial information for the understanding of their origin. Air showers were measured simultaneously with the SPASE-2 air shower array and the AMANDA-B10 Cherenkov telescope at the South Pole. This combination has the advantage to sample almost all high-energy shower muons and is thus a new approach to the determination of the cosmic ray composition. The change in the cosmic ray mass composition was measured versus existing data from direct measurements at low energies. Our data show an increase of the mean log atomic mass lnA by about 0.8 between 500 TeV and 5 PeV. This trend of an increasing mass through the “knee” region is robust against a variety of systematic effects.     

Ultra high energy cosmic rays and the Galactic halo

An analysis has been made of the fraction of ultra high energy cosmic rays (above 10¹⁸ eV) which could be due to processes involved in two possible ‘Models’. The first is the Giant Magnetic Halo Model and the second is the Dark Matter Halo Model. We find that the former, in which heavy nuclei are trapped in a giant halo, fails for energies above about 3 × 10¹⁹ eV. For the Dark Matter Halo Model, in which relic particles follow the “conventional” dark matter and whose decays give ultra high energy cosmic rays, the predicted anisotropies are much higher than those observed. The lack of observation of a finite flux from the Andromeda Galaxy means that the conclusion is insensitive to the spatial scale size of the assumed halo distribution. It is concluded that less than 10% of the ultra high energy cosmic rays come from relic particles in the Galactic halo.     

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.     

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.     

Implications of particle acceleration in active galactic nuclei for cosmic rays and high energy neutrino astronomy

We consider the production of high energy neutrinos and cosmic rays in radio-quiet active galactic nuclei (AGN) or in the central regions of radio-loud AGN. We use a model in which acceleration of protons takes place at a shock in an accretion flow onto a supermassive black hole, and follow the cascade that results from interactions of the accelerated protons in the AGN environment. We use our results to estimate the diffuse high energy neutrino intensity and cosmic ray intensity due to AGN. We discuss our results in the context of high energy neutrino telescopes under construction, and measurements of the cosmic ray composition in the region of the “knee” in the energy spectrum at 10⁷ GeV.     

Identification of “small” dust impacts in the Ulysses dust detector data

Since October 1990, 3 weeks after the launch of the Ulysses spacecraft, the dust detector onboard recorded impacts of cosmic dust particles. Besides dust impacts, the detector recorded noise from a variety of sources. So far, a very rigid scheme had been applied to eliminate noise from impact data. The data labeled “big” dust impacts previously led to the identification of interstellar dust and of dust streams from Jupiter. The analysis presented here is concerned with data of signals of small amplitudes which are strongly contaminated by noise. Impacts identified in this data set are called “small” impacts. It is shown that dust impacts can be clearly distinguished from noise for most of the events due to the multi-coincidence characteristics of the instrument. 516 “small” impacts have been identified. For an additional 119 events, strong arguments can be given that they are probably small dust impacts. Thereby, the total number of dust impacts increases from 333 to 968 in the time period from 28 October 1990 to 31 December 1992. This increase permits a better statistical analysis, especially of the Jupiter dust streams which consist mostly of small and fast particles. Additional dust streams have been identified between the already known streams before and after Jupiter flyby. The dependence of the deflection from the Jupiter direction, the stream intensity and width on Jupiter distance support the assertion that they have been emitted from the Jovian system. The masses of the 635 “small” dust particles range from 6 × 10⁻¹⁷ to 3 × 10⁻¹⁰ g with a mean value of 1 × 10⁻¹² g, which compares to a range from 1 × 10⁻¹⁶ to 4 × 10⁻⁹ g with a mean value of 2 × 10⁻¹¹ g for the previously identified 333 “big” dust particles.     

The cosmic ray energy spectrum between 10 and 10 eV

The energy spectrum of cosmic rays with primary energies between 10¹⁴ eV and 10¹⁶ eV has been studied with the CASA-MIA air shower array. The measured differential energy spectrum is a power law (dj/dE ∝ E^−y) with spectral indices γ of 2.66±0.02 below approximately 10¹⁵ eV and 3.00±0.05 above. A new method is used for measuring primary energy derived from ground-based data in a compositionally insensitive way. In contrast with some previous reports, the “knee” of the energy spectrum does not appear sharp, but rather a smooth transition over energies from 10¹⁵ eV to 3.0 × 10¹⁵ eV.     

Extragalactic UHE proton spectrum and prediction for iron-nuclei flux at 10–10 GeV

We investigate the problem of transition from galactic cosmic rays to extragalactic ultra-high energy cosmic rays. Using the model for extragalactic ultra-high energy cosmic rays and observed all-particle cosmic ray spectrum, we calculate the galactic spectrum of iron nuclei in the energy range 10⁸–10⁹ GeV. The flux and spectrum predicted at lower energies agree well with the KASCADE data. The transition from galactic to extragalactic cosmic rays is distinctly seen in spectra of protons and iron nuclei, when they are measured separately. The shape of the predicted iron spectrum agrees with the Hall diffusion.     

Decaying nuclei and the age of cosmic rays in the galaxy

The propagation of radioactive nuclei of cosmic rays in a flat diffusion galactic model (sources and the main gaseous mass are concentrated in the galactic disc) is considered. The corresponding results are not reducible to the results of a simple homogeneous model. It is shown that the recent data on the Be¹⁰ nuclei abundance in cosmic rays do not contradict the occurrence of a large cosmic ray halo.     

Transient perturbations and their effects in the heliosphere, the geo-magnetosphere, and the Earth’s atmosphere: Space weather perspective

We discuss the effects of certain dynamic features of space environment in the heliosphere, the geo-magnetosphere, and the earth’s atmosphere. In particular, transient perturbations in solar wind plasma, interplanetary magnetic field, and energetic charged particle (cosmic ray) fluxes near 1 AU in the heliosphere have been discussed. Transient variations in magnetic activity in geo-magnetosphere and solar modulation effects in the heliosphere have also been studied. Emphasis is on certain features of transient perturbations related to space weather effects. Relationships between geomagnetic storms and transient modulations in cosmic ray intensity (Forbush decreases), especially those caused by shock-associated interplanetary disturbances, have been studied in detail. We have analysed the cosmic ray, geomagnetic and interplanetary plasma/field data to understand the physical mechanisms of two phenomena namely, Forbush decrease and geomagnetic storms, and to search for precursors to Forbush decrease (and geomagnetic storms) that can be used as a signature to forecast space weather. It is shown that the use of cosmic ray records has practical application for space weather predictions. Enhanced diurnal anisotropy and intensity deficit of cosmic rays have been identified as precursors to Forbush decreases in cosmic ray intensity. It is found that precursor to smaller (less than 5%) amplitude Forbush decrease due to weaker interplanetary shock is enhanced diurnal anisotropy. However, larger amplitude (greater than 5%) Forbush decrease due to stronger interplanetary shock shows loss cone type intensity deficit as precursor in ground based intensity record. These precursors can be used as inputs for space weather forecast.     

The measurement of the decorrelation function in underground muon pairs as a probe of primary cosmic ray interactions

We discuss the measurement of the average angle within pairs of muons detected in underground experiments versus their relative distance, as a tool for studying primary Cosmic Ray interactions. We propose to call Φ(d) “decorrelation” function. Under simplifying approximation, the decorrelation function can be computed analytically to show its dependence on details of cosmic ray interactions as well as on the propagation in the rock. We argue that this new measurement is useful in a real experiment, complementing traditional ones like the decoherence function.     

HI imaging the low red-shift cosmic web

Only in recent years has the realization emerged that galaxies do not dominate the universal baryon budget but are merely the brightest pearls of an underlying cosmic web. Although the gas in these inter-galactic filaments is moderately to highly ionized, QSO absorption lines have shown that the surface area increases dramatically in going down to lower HI column densities. The first image of the cosmic web in HI emission has just been made of the Local Group filament connecting M31 and M33. The corresponding HI distribution function is in very good agreement with that of the QSO absorption lines, confirming the 30-fold increase in surface area expected between 10¹⁹ and 10¹⁷ cm⁻². The critical observational challenge is crossing the “HI desert”, the range of log(N_HI) from about 19.5 down to 18, over which photo-ionization by the intergalactic radiation field produces an exponential decline in the neutral fraction from essentially unity down to a few percent. Nature is kinder again to the HI observer below log(N_HI) = 18, where the neutral fraction decreases only very slowly with log(N_HI). With the SKA, we can begin the systematic study of the cosmic web beyond the Local Group. With moderate integration times, the necessary resolution and sensitivity can be achieved out to distances beyond the Virgo cluster. When combined with targeted optical and UV absorption line observations, the total baryonic masses and enrichment histories of the cosmic web will be determined over the complete range of environmental over-densities.     

Chemical composition of primary cosmic rays with energies from 10 to 10 eV

The chemical composition of primary cosmic rays with energies from 10¹⁵ to 10^16.5 eV, so called “knee” region, is examined. We have observed the time structures of air Čerenkov light associated with air showers at Mt. Chacaltaya, Bolivia, since 1995. The distribution of a parameter that characterizes the observed time structures is compared with that calculated with a Monte Carlo technique for various chemical compositions. Then the energy dependence of the average logarithmic mass numbers ln A of the primary cosmic rays is determined. The present result at 10^15.3 eV is almost consistent with the result of JACEE (A12) and shows gradual increase in ln A as a function of the primary energy (A24 at 10¹⁶ eV). Form the comparison of the observational results with several theoretical models, we conclude that the supernova explosion of massive stars is a plausible candidate for the origin of cosmic rays around the “knee” region.     

Cosmic ray composition and energy spectrum from 1–30 PeV using the 40-string configuration of IceTop and IceCube

The mass composition of high energy cosmic rays depends on their production, acceleration, and propagation. The study of cosmic ray composition can therefore reveal hints of the origin of these particles. At the South Pole, the IceCube Neutrino Observatory is capable of measuring two components of cosmic ray air showers in coincidence: the electromagnetic component at high altitude (2835 m) using the IceTop surface array, and the muonic component above ∼1 TeV using the IceCube array. This unique detector arrangement provides an opportunity for precision measurements of the cosmic ray energy spectrum and composition in the region of the knee and beyond. We present the results of a neural network analysis technique to study the cosmic ray composition and the energy spectrum from 1 PeV to 30 PeV using data recorded using the 40-string/40-station configuration of the IceCube Neutrino Observatory.     

Power spectra of cosmic ray intensity variations at the ground and atmospheric pressure oscillations

Power spectrum estimates are computed for cosmic ray and pressure variations in the frequency range of 1.6 × 10⁻⁶ to 4.15 × 10⁻³ c/sec for three data sets each of 27 days length (2 min interval) recorded at Chacaltaya (Bolivia) during 1965–1966. The general trend of these spectra showy ƒ⁻ fit having exponent values from −1.5 to −1.9 for L.F. side. From the semi-diurnal peaks both in cosmic rays and pressure, the average value of pressure coefficient is found to be 0.3 per cent/mm of Hg. For the same three sets of data, a detailed analysis of cross-spectrum, coherence and phase relationship between cosmic rays and pressure is carried out in the frequency range of 1–12 c/day. Besides semi-diurnal peak, variations of 3, 4, 7 and 8 c/day are found to be common both in cosmic rays and pressure. The value of pressure coefficient and residual amplitude of cosmic rays for these particular frequencies are calculated.     

ISOCAM Long Wave Detector Glitches: Data and Ray-Tracing Simulations

The affected pixel number distribution of the ISOCAM Long-Wave (LW) detector in cosmic ray-induced glitchesis calculated. The methods employed are Monte Carlo ray-tracing techniques and the taxi metric, which allows direct calculation of the affected pixels based on the knowledge of entry and exit points of the ray. The simulation results are compared with long-term experimental data. Based on the simulations combined with the CREME96 cosmic ray model for solar quiet period, it is estimated that the detector is on average traversed by 0.3 cosmic ray protons per second. For the experimental data, a corresponding minimum bound of 0.14 cosmic ray protons per second is obtained. Deviations in the simulated and measured pixel number distributions are discussed.     

Constraints on structure formation models from the Sunyaev–Zel'dovich effect

In the context of cold dark matter (CDM) cosmological models, we have simulated images of the brightness temperature fluctuations in the cosmic microwave background (CMB) sky owing to the Sunyaev–Zel'dovich (S–Z) effect in a cosmological distribution of clusters. We compare the image statistics with recent ATCA limits on arcmin-scale CMB anisotropy. The S–Z effect produces a generically non-Gaussian field and we compute the variance in the simulated temperature-anisotropy images, after convolution with the ATCA beam pattern, for different cosmological models. All the models are normalized to the 4-yr COBE data. We find an increase in the simulated-sky temperature variance with increase in the cosmological density parameter Ω₀. A comparison with the upper limits on the sky variance set by the ATCA appears to rule out our closed-universe model: low-Ω₀ open-universe models are preferred. The result is independent of any present day observations of σ ₈.     

Structure in the EAS size spectrum: analysis of the CASA-MIA results

We present an analysis of the very recent data from the CASA-MIA extensive air shower array (M.A.K. Glasmacher, Ph.D. Thesis, University of Michigan (1998) in terms of a search for the structure in the shower size spectrum claimed by us from an analysis of “the world's data”. Our earlier claim is found to be supported to the extent that there is strong evidence for the existence of structure in the spectrum which cannot obviously be explained by the conventional Galactic Modulation Model. There is modest evidence for the structure being of “our” form and strong support for “our” mass composition when “corrected” to the interaction model advocated by us. None of the results are inconsistent with there having been a recent, nearby, single supernova.     

The energy scale of inflation: is the hunt for the primordial B-mode a waste of time?

Recent theoretical results indicate that the detection of primordial gravity waves from inflation may be a hopeless task. First, foregrounds from lensing put a strict lower limit on the detectability of the B-mode polarization signal in the cosmic microwave background, the “smoking gun” for tensor (gravity wave) fluctuations. Meanwhile, widely accepted theoretical arguments indicate that the amplitude of gravity waves produced in inflation will be below this limit. I argue that failure is not inevitable, and that the effort to detect the primordial signal in the B-mode, whether it succeeds or fails, will yield crucial information about the nature of inflation.     

Origin of cosmic rays: Galactic models with halo

Many years ago physical and radio-astronomical arguments and data led to the assumption that cosmic rays in the Galaxy (and probably in other galaxies) fill a more or less extended halo, but are not concentrated in the disk. It was not so long ago, however, that the existence of a radio-halo was discovered, in which the effective dimensions increase with a decrease in frequency. The frequency decrease occurs when relativistic electrons diffuse from the disk, losing energy due to bremsstrahlung and Compton scattering.Meanwhile, some ambiguity on the question of the existence of a radio-halo, and other reasons, have led to a rather wide use of disk models, particularly those in which cosmic rays are present in the Galaxy only for a periodT _cr,d3×10⁶ yr. The authors have repeatedly stated the inadmissibility of such models and, generally, a homogeneous (leaky box) model for the origin of cosmic rays. The new data concerning the amount of radioactive¹⁰Be nuclei in cosmic rays near the Earth in no way contradict the halo models in which the lifetime of cosmic rays isT _cr,h10⁸ yr. In connection with the continuing controversy, the present paper is devoted to a detailed consideration of the difference between the homogeneous and diffusion models. Within the latter models some calculations on the chemical composition of cosmic rays have been carried out, which concern not only stable but also radioactive isotopes.