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.     

High energy cosmic ray spectroscopy. IV. The evidence from direct observations at lower energies and directional anisotropies

In previous papers (A.D. Erlykin, A.W. Wolfendale, Astropart. Phys. 7 (1997) 1, 203; 8 (1998) 265; J. Phys. G 23 (1997) 979), we presented evidence for structure in the size spectrum of cosmic ray air showers which we interpreted as due to the presence of oxygen and iron nuclei from a local, recent, supernova remnant. Although the energies in question are 3 × 10¹⁵ eV and 1.2 × 10¹⁶ eV, well above those where direct measurements are possible, the direct measurements are, in fact, relevant. We find that the direct measurements are quite consistent with an extrapolation back of our spectra. Indeed, taken alone, the direct measurements themselves provide strong evidence for the existence of an extra, single source contribution to the total energy spectrum. The paper also includes a discussion of the high energy electron spectrum, anisotropies and the likely site of the local SNR.     

Changes of the cosmic-ray mass composition in the 10–10 eV energy range

Data taken with ten Cosmic Ray Tracking (CRT) detectors and the HEGRA air-shower array on La Palma, Canary Islands, have been analysed to investigate changes of the cosmic ay mass composition at the ‘knee’ of the cosmic-ray flux spectrum near 10¹⁵ eV energy. The analysis is based on the angular distributions of particles in air showers. HEGRA data provided the shower size, direction, and core position and CRT data the particle track information. It is shown that the angular distribution of muons in air showers is sensitive to the composition over a wide range of shower sizes and, thus, primary cosmic-ray energies with little systematic uncertainties. Results can be easily expressed in terms of ln A of primary cosmic rays. In the lower part of the energy range covered, we have considerable overlap with direct composition measurements by the JACEE collaboration and find compatible results in the observed rise of ln A. Above about 10¹⁵ eV energy we find no or at most a slow further rise of ln A. Simple cosmic-ray composition models are presented which are fully consistent with our results as well as the JACEE flux and composition measurements and the flux measurements of the Tibet ASγ collaboration. Minimal three-parameter composition models defined by the same power-law slope of all elements below the knee and a common change in slope at a fixed rigidity are inconsistent with these data.     

Flux limits on ultra high energy neutrinos with AMANDA-B10

Data taken during 1997 with the AMANDA-B10 detector are searched for a diffuse flux of neutrinos of all flavors with energies above 10¹⁶ eV. At these energies the Earth is opaque to neutrinos, and thus neutrino induced events are concentrated at the horizon. The background are large muon bundles from down-going atmospheric air shower events. No excess events above the background expectation are observed and a neutrino flux following E⁻², with an equal mix of all flavors, is limited to E²Φ(10¹⁵ eV < E < 3 × 10¹⁸ eV)  0.99 × 10⁻⁶ GeV cm⁻² s⁻¹ sr⁻¹ at 90% confidence level. This is the most restrictive experimental bound placed by any neutrino detector at these energies. Bounds to specific extraterrestrial neutrino flux predictions are also presented.     

The mass composition of cosmic rays near 10 eV as deduced from measurements made at Volcano Ranch

Linsley used the Volcano Ranch array to collect data on the lateral distribution of showers produced by cosmic rays at energies above 10¹⁷ eV. Very precise measurements of the steepness of the lateral distribution function were made on 366 events. The current availability of sophisticated hadronic interaction models has prompted an interpretation of the measurements. In this analysis we use the Monte Carlo code to generate showers, together with to simulate the detector response to ground particles. The results show that, with the assumption of a bi-modal proton and iron mix, iron is the dominant component of cosmic rays near 10¹⁸ eV, assuming that hadronic interactions are well-described by at this energy range. The Volcano Ranch data set, as available to us, does not allow a straightforward assignment of energy for each event. It is thus not possible to give the energy dependence of the mass composition.     

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.     

Angular correlation of ultra-high energy cosmic rays with compact radio-loud quasars

Angular correlations of ultra-high energy cosmic rays with cosmologically distant sources may provide clues to these mysterious events. We compare cosmic ray tracks with energies above 10²⁰ eV to a compilation of radio-loud compact QSO positions. The statistical method emphasizes invariant quantities and a test of statistical independence of track and source distributions. Statistical independence is ruled out by several independent statistics at confidence levels of orders 97–99%, indicating that track directions and QSO source positions are correlated at a significant level.     

An auroral F-region study using in situ measurements by the Atmosphere Explorer-C satellite

On 14 July 1974 the Atmosphere Explorer-C satellite flew through an aurora at F-region altitudes just after local midnight. The effects of the particle influx are clearly evident in the ion densities, the 6300 Å airglow, and the electron and ion temperatures. This event provided an opportunity to study the agreement between the observed ion densities and those calculated from photochemical theory using in situ measurements of such atmospheric parameters as the neutral densities and the differential electron energy spectra obtained along the satellite track. Good agreement is obtained for the ions O₂⁺, NO⁺ and N₂⁺ using photochemical theory and measured rate constants and electron impact cross sections. Atomic nitrogen densities are calculated from the observed [NO⁺]/[O₂⁺] ratio. In the region of most intense electron fluxes (20 erg cm⁻² sec⁻¹) at 280 km, the N density is found to be between 2 and 7 × 10⁷ cm⁻³. The resulting N densities are found to account for approx. 60% of the production of N⁺ through electron impact on N and the resonant charge exchange of O⁺(²P) with N(⁴S). This reaction also provides a significant source of O(¹S) in the aurora at F-region altitudes. In the region of intense fast electron influx, the reaction with atomic nitrogen is found to be the main loss of O⁺(²P).     

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.     

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.     

Simulations of a giant hybrid air shower detector

We have used Monte Carlo simulations to investigate the capabilities of a giant air shower observatory designed to detect showers initiated by cosmic rays with energies exceeding 10¹⁹ eV. The observatory is to consist of an array of detectors that will characterise the air shower at ground level, and optical detectors to measure the fluorescence light emitted by the shower in the atmosphere. Using these detectors together in a ‘hybrid’ configuration, we find that precise geometrical reconstruction of the shower axis is possible, leading to excellent resolution in energy. The technique is also shown to provide very good reconstruction below 10¹⁹ eV, at energies where the ground array is not fully efficient.     

Reanalysis of GU miniarray data using CORSIKA

The GU miniarray is a ultra high energy cosmic ray (UHECR) detector consisting of eight plastic scintillators of carpet area 2 m², each viewed by fast PMTs. It is used to detect Giant EAS by the method of time spread measurement of secondary particles produced in the atmosphere. The energies of the air showers have been reestimated using CORSIKA program. As in the original analysis the Cosmic Ray energy was determined via its relation to the ground level parameter N_s, the shower size. This relation was obtained previously through a best fit relation in agreement with QGS model and Yakutsk data. In this work we use CORSIKA code with QGSJET model of high energy hadronic interactions to simulate miniarray data leading to a modified relation between primary energy and shower size. A revised energy spectrum is reported for 10¹⁷–10¹⁹ eV primary energy.     

Comparison of AGASA data with CORSIKA simulation

An interpretation of Akeno giant air shower array (AGASA) data by comparing the experimental results with the simulated ones by cosmic ray simulation for KASCADE (CORSIKA) has been made. General features of the electromagnetic component and low energy muons observed by AGASA can be well reproduced by CORSIKA. The form of the lateral distribution of charged particles agrees well with the experimental one between a few hundred metres and 2000 m from the core, irrespective of the hadronic interaction model studied and the primary composition (proton or iron). It does not depend on the primary energy between 10^17.5 and 10²⁰ eV as the experiment shows. If we evaluate the particle density measured by scintillators of 5 cm thickness at 600 m from the core S₀(600), suffix 0 denotes the vertically incident shower) by taking into account the similar conditions as in the experiment, the conversion relation from S₀(600) to the primary energy is expressed as E (eV)=2.15×10¹⁷S₀(600)^1.015 within 10% uncertainty among the models and composition used, which suggests the present AGASA conversion factor is the lower limit. Although the form of the muon lateral distribution fits well to the experiment within 1000 m from the core, the absolute values change with hadronic interaction model and primary composition. The slope of the ρ_μ(600) (muon density above 1 GeV at 600 m from the core) vs. S₀(600) relation in experiment is flatter than that in simulation of any hadronic model and primary composition. As the experimental slope is constant from 10¹⁵ to 10¹⁹ eV, we need to study this relation in a wide primary energy range to infer the rate of change of chemical composition with energy.     

A possible gas for solar neutrino spectroscopy

We discuss the possibility of using pure CF₄ to fill a 2000 m³ Time Projection Chamber in order to detect the solar neutrinos through the elastic scattering v_ee⁻ → v_ee⁻, with the threshold of 100 keV on the kinetic energy of the recoiling electron. In a volume of 2000 m³ of CF₄ at normal pressure and room temperature, which corresponds to a mass of 7.4 ton, we expect ~ 3300 of such events per year. The detector can give the spectrum of the low energy neutrinos from the Sun and it can identify solar neutrinos of different origin: pp, ⁷Be, and, eventually, ⁸B. We find that ¹⁴C is a possible severe source of background: it is necessary to have a ratio ¹⁴C/¹²C lower than 10⁻¹⁹ in order to be able to identify the pp neutrinos.     

Space–time foam and cosmic-ray interactions

It has been proposed that propagation of cosmic-rays at extreme-energy may be sensitive to Lorentz-violating metric fluctuations (“foam”). We investigate the changes in interaction thresholds for cosmic-rays and gamma-rays interacting on the CMB and IR backgrounds, for a class of stochastic models of space–time foam. The strength of the foam is characterized by the factor (E/M_P)^a, where a is a phenomenological suppression parameter. We find that there exists a critical value of a (dependent on the particular reaction: a_crit3 for cosmic-rays, 1 for gamma-rays), below which the threshold energy can only be lowered, and above which the threshold energy may be raised, but at most by a factor of two. Thus, it does not appear possible in this class of models to extend cosmic-ray spectra significantly beyond their classical absorption energies. However, the lower thresholds resulting from foam may have signatures in the cosmic-ray spectrum. In the context of this foam model, we find that cosmic-ray energies cannot exceed the fundamental Planck scale, and so set a lower bound of 10⁸ TeV for the scale of gravity. We also find that suppression of p→pπ⁰ and γ→e⁻e⁺ “decays” favors values aa_crit. Finally, we comment on the apparent non-conservation of particle energy–momentum, and speculate on its re-emergence as dark energy in the foamy vacuum.     

A comparison of cosmic ray composition measurements at the highest energies

In recent years the Fly's Eye and Akeno groups have presented analyses of the cosmic ray mass composition at energies above 10¹⁷ eV. While the analysis of the Fly's Eye group points to a likely change in mass composition from heavy to light at energies above 10¹⁸ eV, the Akeno analysis favours an unchanging composition. However, the two groups base their conclusions on simulations using quite different hadronic models. Here we present a comparison of the experiments using the same hadronic model and find that the agreement between the experiments is much improved. Under this model, both experiments measure a composition rich in iron around 10¹⁷ eV which becomes lighter at higher energies. However, the agreement is not complete, which indicates scope for improvement of the interaction model, or perhaps the need for a re-examination of the experimental results.     

The anisotropy of cosmic ray arrival directions around 10 eV

Anisotropy in the arrival directions of cosmic rays with energies above 10¹⁷ eV is studied using data from the Akeno 20 km² array and the Akeno Giant Air Shower Array (AGASA), using a total of about 114 000 showers observed over 11 years. In the first harmonic analysis, we have found a strong anisotropy of 4% around 10¹⁸ eV, corresponding to a chance probability of 0.2% after taking the number of independent trials into account. with two-dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region.     

Deactivation of N2A Σumolecules in the aurora

Recent rocket observations of the N₂ V-K (Vegard-Kaplan) system in the aurora have been reinterpreted using an atmospheric model based on mass spectrometer measurements in an aurora of similar intensity at the same time of year. In contrast to the original interpretation, we find that population by cascade from the C³Π_u and B³Π_g states in the A³Σ_u⁺v=0,1 levels, as calculated using recently measured electron excitation cross sections, accurately accounts for the observed relative emission rates (IV-K/12PG_0.0). In addition there is no need to change the production rate of A ³ Σ _u⁺ molecules relative to that of C³Π_uv=0 as a function of altitude in order to fit the profile of the deactivation probability to the atmospheric model. Quenching of A ³ Σ _u⁺ molecules at high altitudes is dominated by atomic oxygen. The rate constants for the v=0 and v=1 levels are 8 × 10⁻¹¹ cm³ sec⁻¹ and 1.7 × 10⁻¹⁰ cm³ sec⁻¹ respectively, as determined using the model atmosphere mentioned above. Recent observations with a helium cooled mass spectrometer suggest that conventional mass spectrometer measurements tend to underestimate the atomic oxygen relative concentration. The rate coefficients may therefore be too large by as much as a factor of 3. Below 130 Km we find that it is possible to account for the deactivation in bright auroras by invoking large nitric oxide concentrations, similar to those recently observed mass spectrometrically and using a rate constant of 8 × 10⁻¹¹ cm³ sec⁻¹ for both the v=1 levels. This rate constant is very nearly the same as that measured in the laboratory (7 × 10⁻¹¹ cm³ sec⁻¹). Molecular oxygen appears not to play a significant role in deactivating the lower A ³ Σ _u⁺ levels.     

Properties of the infrared dark cloud G79.2+0.38

The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the ¹²CO and its isotopic molecules ¹³CO and ¹⁸CO. The dense molecular cores defined by the C¹⁸O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H₂) (5 – 12) × 10²² cm⁻² and the mean number density n (3 ± 1) × 10⁴ cm⁻³. Their sizes are 1.7 and 1.2 pc in ¹³CO(1-0) line, 1.2 and 0.6 pc in C¹⁸O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 10² to 2 × 10³ M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p^{−0.34±0.02}. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules ¹³CO and C¹⁸O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X_13/18 between ¹³CO and C¹⁸O with the column density.     

Resonance scattering by N2

A study is made of the intensity distribution among the bands of the Meinel and first negative system of N₂⁺ due to resonance scattering of sunlight. Absolute transition probabilities are used to calculate the relative populations among the ion states under resonance scattering conditions; the mean lifetime for deactivation is the parameter which determines the amount of resonance scattering. Photon scattering rates are calculated for most of the ion bands and it is suggested that an appropriate value for the 3914 Å band would be 0·050 photons/ sec per ion. Observations of the Δυ = −1 sequence of the first negative system in the twilight spectrum are reported. Extended vibrational development is detected which indicates that only about 80 per cent of the emission is resonance scattered. Sunlit auroral spectra of N₂⁺, however, which have been generally considered to be due predominantly to resonance scattering, indicates only about 40 per cent of the emission is due to resonance scattering. Measurable effects resulting from a charge-transfer ion source (O⁺(²D)) are predicted.