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.     

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.     

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.     

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.     

Gamma ray spectrum of the crab nebula in the multi TeV region

The THEMISTOCLE array of 18 Cherenkov detectors which has a 3 TeV gamma energy threshold, has detected a signal from the Crab nebula at a 5.8 standard deviation level. Information on the energy spectrum is obtained in the range 3–15 TeV. The integrated flux can be fitted with the form, Φ (> E) = (3.7 ± 0.5) × 10^-12 (E/5)-^{−1.5 ± 0.20} cm⁻² s⁻¹ (E in TeV) compatible with the extrapolation of results at lower energies. The Crab signal is used to measure the angular resolution of the multi-telescope technique. The value obtained is 2.3 mr (0.15°) in agreement with the results of simulations, and confirms the interest of this new method for multi-TeV gamma-ray detection.     

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.     

Dynamics of the O H2O charge-transfer reaction and its implications for space-borne measurements

O⁺(⁴S) + H₂O charge-transfer cross-section and product-ion time-of-flight measurements are presented over a center-of-mass collision energy range of 0.2–15 eV. The results are obtained with a newly constructed guided-ion beam experiment. The measured energy dependence of the charge-transfer crosssection agrees well with previous measurements and ADO model predictions of Bates [Chem. Phys. Lett. 82, 396 (1981) Proc. R. Soc. Lond. A384, 289 (1982) Chem. Phys. Lett. 111, 428 (1984)] at low collision energies, but exhibits significant differences at collision energies above 2 eV. The product kinetic energy analysis shows that most of the product ions are produced with near-thermal kinetic energies. At low collision energies, a significant backscattered channel is observed that is associated with complex formation. This channel exhibits a cross-section of approximately 1.4 Ųat a collision energy of 2.65 eV, corresponding to a laboratory ion energy of 5 eV.     

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.     

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.     

The HNCO ring in the Sgr B2 region

The large-scale structure associated with the 2′N HNCO peak in Sgr B2 [Minh, Y.C., Haikala, L., Hjalmarson, Å., Irvine, W.M., 1998. ApJ 498, 261 (Paper I)] has been investigated. A ring-like morphology of the HNCO emission has been found; this structure may be colliding with the Principal Cloud of Sgr B2. This “HNCO Ring” appears to be centered at (l,b) = (0.7°,−0.07°), with a radius of 5 pc and a total mass of 1.0 × 10⁵ to 1.6 × 10⁶ M. The expansion velocity of the Ring is estimated to be 30–40 km s⁻¹, which gives an expansion time scale of 1.5 × 10⁵ year. The morphology suggests that collision between the Ring and the Principal Cloud may be triggering the massive star formation in the Sgr B2 cloud sequentially, with the latest star formation taking place at the 2′N position. The chemistry related to HNCO is not certain yet, but if it forms mainly via reaction with the evaporated OCN⁻ from icy grain mantles, the observed enhancement of the HNCO abundance can be understood as resulting from shocks associated with the collision between the Principal Cloud and the expanding HNCO Ring.     

Analysis of the variation in inclination of Intercosmos 13 Rocket, 1975-22B

The orbit of Intercosmos 13 rocket (1975-22B) has been determined at 103 epochs between 30 April 1975 and 10 April 1980 from almost 7000 observations. One hundred and three values of inclination have been determined and corrections incoporated for the effects due to zonal harmonic, lunisolar and tesseral harmonic perturbations, precession, and solid Earth tides. The modified data have been analysed to yield values of the atmospheric rotation rate, Λ rev day⁻¹, viz. Λ = 0.94 ± 0.10 at an average height of 322 ± 6 km and Λ = 1.27 ± 0.02 at 288 km. Analysis of the inclination near 14th-order resonance has indicated lumped harmonic values 10^{9 1.0}_1.4 = − 76.13 ± 12.47, 10^{9 1,0}₁₄ = − 29.89 ± 32.64, 10^{9 −1.2}₁₄ = − 63.11 ± 15.44 10^{9 −1.2}₁₄ = − 32.52 ± 26.96, for inclination 82.952°.     

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.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.     

Gamma-ray bursters as sources of cosmic rays

From the little we know of the physical conditions in γ-ray bursters, it seems that they are potentially effective in the acceleration of high-energy cosmic rays (CRs), especially if the bursters are at cosmological distances. We find that, with the observed statistics and fluxes of γ-ray bursts, cosmological bursters may be an important source of cosmic rays in two regions of the observed spectrum: (1) At the very-high-energy end (E > 10¹⁹ eV), where CRs must be of extragalactic origin. (2) Around and above the spectral feature that has been described as a bump and/or a knee, which occurs around 10¹⁵ eV. The occasional bursters that occur inside the Galaxy — about once in a few hundred thousand years if burst emission is isotropic; more often, if it is beamed — could maintain the density of galactic cosmic rays at the observed level in this range. These two energy ranges might correspond to two typical CR energy scales characteristic of bursters: one pertinent to CR acceleration due to interaction of a magnetized-fireball front with an ambient medium; the other to acceleration in the fireball itself (e.g. shock acceleration).     

The potential of the ground based arrays of imaging atmospheric Cherenkov telescopes. II. Gamma ray flux sensitivities

We discuss the capability of ‘100 GeV’ class imaging atmospheric Cherenkov telescope (IACT) arrays as future powerful instruments of ground-based gamma-ray astronomy. It is assumed that the array is gathered from individually triggered quadrangular 4-IACT ‘cells’ with a linear size of about 100 m. The multi-cell concept allows coverage of large detection areas economically, and at the same time the effective exploitation of the stereoscopic approach of determination of the shower parameters using information obtained by several IACTs simultaneously. Determination of arrival directions of γ-ray primaries on an event-by-event basis with accuracy δθ ≤ 0.1° combined with high suppression efficiency (at both the hardware and software levels) of the background hadronic showers by a factor of ≈ 10³, and large, up to 1 km² collection areas, can provide minimum detectable energy fluxes of ≥ 100 GeV γ-rays from point sources down to 10⁻¹³ erg/cm² s which is about 3 orders of magnitude lower than the current sensitivities achieved by the satellite-borne detectors at MeV and GeV energies. High sensitivities of multi-IACT arrays would partially compensate the limited efficiency of the technique for all-sky surveys, as well as allow study of moderately extended (≤ 1°) γ-ray sources. IACT arrays with minimum detectable fluence of ≥ 100 GeV γ-rays S_γ < 10⁻⁸ erg/cm² are well suited for effective exploration of highly sporadic nonthermal phenomena from different classes of astrophysical objects on time-scales from ≤ 1 s to several minutes.     

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.     

Primary proton spectrum around the knee deduced from the emulsion-chamber data obtained at Mts. Fuji and Kanbala

We have done extensive Monte Carlo simulations using the new simulation codes of CORSIKA and COSMOS to compare with the gamma-family data obtained at Mts. Fuji (3750 m above sea level) and Kanbala (5500 m above sea level). Then, we estimated the primary proton and helium spectra around the knee energy region using a multiple-layered feed-forward neural network as a classifier of primary particle kind. The selection efficiency of proton-induced family events is estimated to be 82%. The flux value of protons at 2×10¹⁵ eV is (5.5±1.5)×10⁻¹⁴ (m⁻² s⁻¹ sr⁻¹ GeV⁻¹). The result suggests heavy-enriched primary composition around the knee region.     

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.