Propagation of cosmic ray nuclei above 50 GeV/amu

It is now well-known that the relative abundance of secondary cosmic ray nuclei which are produced by spallation in interstellar space depends on energies. The mechanism responsible for this energy dependence are not well understood because of the scarcity of data on the cosmic ray composition at energy above 50 GeV/amu. The University of Chicago experiment on the Spacelab-2 flown by the Space Shuttle in 1985 is aimed at the measurement of the elemental composition in this unexplored energy regime. Our latest results on the energy dependence of the secondary to primary abundance ratio using the pairs of elements B,C and N,O, are discussed.     

Measurement of cosmic ray elemental composition from the CAKE balloon experiment

CAKE (Cosmic Abundances below Knee Energies) was a prototype balloon experiment for the determination of the charge spectra and abundances of the primary cosmic rays (CR) with Z > 10. It was a passive instrument made of layers of CR39® and Lexan®/Makrofol® nuclear track detectors; it had a geometric acceptance of ∼0.7 m² sr for Fe nuclei. Here, the scanning and analysis strategies, the algorithms used for the off-line filtering and for the tracking in automated mode of the primary cosmic rays are presented, together with the resulting CR charge distribution and their abundances.     

Precise measurements of the cosmic ray antiproton spectrum with BESS including the effects of solar modulation

The Balloon Borne Experiment with a Superconducting Spectrometer (BESS) has measured the energy spectrum of cosmic-ray antiprotons between 0.18 and 4.20 GeV in eight flights between 1993 and 2002. Above about 1 GeV, models in which antiprotons are secondary products of the interactions of primary cosmic rays with the interstellar gas agree with the BESS antiproton spectrum. Below 1 GeV, the data show a possible excess antiproton flux compared to secondary model predictions, suggesting the presence of an additional source of antiprotons. The antiproton/proton ratios measured between 1993 and 1999, during the Sun’s positive-polarity phase, are consistent with simple models of solar modulation. However, results from the 2000 flight, following the solar magnetic field reversal, show a sudden increase in the antiproton/proton ratio and tend to favor a charge-sign-dependent drift model. To extend BESS measurements to lower energies, an evolutionary instrument, BESS-Polar, is under construction for polar flight in 2004.     

Galactic cosmic ray composition and energy spectra.

Galactic cosmic ray nuclei represent a significant risk to long-duration spaceflight outside the magnetosphere. We review briefly existing measurements of the composition and energy spectra of heavy cosmic ray nuclei, pointing out which species and energy ranges are most critical to assessing cosmic ray risks for spaceflight. Key data sets are identified and a table of cosmic ray abundances is presented for elements from H to Ni (Z = 1 to 28). Because of the 22-year nature of the solar modulation cycle, data from the approaching 1998 solar minimum is especially important to reducing uncertainties in the cosmic ray radiation hazard. It is recommended that efforts to model this hazard take advantage of approaches that have been developed to model the astrophysical aspects of cosmic rays.     

Galactic cosmic ray flux simulation and prediction.

A dynamic galactic cosmic ray model is proposed to quantitatively describe the z=1-28 ions and electrons of E=10-10(5) MeV/nucleon and their particle flux variations around the Earth's orbit and beyond the Earth's magnetosphere due to diverse large-scale variations of solar activity factors. The variations of large-scale heliospheric magnetic fields and the galactic cosmic ray flux variation time delays relative to solar activity variations are simulated. The lag characteristics and sunspot number predictions having been determined in detail, the model can be used to predict galactic cosmic ray flux levels.     

Fluxes and nuclear abundances of cosmic rays inside the magnetosphere using a transmission function approach

At 1 AU and outside the Earth’s magnetosphere, the relative abundances to protons for He (He/p), C (C/p) and Fe (Fe/p) nuclei were calculated using the observation data of AMS-01 (for p and He) and HEAO-3 (for C and Fe) above 0.8 GeV/nucleon. In addition, the transmission function (TF) for the GCR propagation inside the magnetosphere was evaluated using the IGRF and T96 (introduced by Tsyganenko and Stern) models to obtain permitted and forbidden trajectories inside the magnetosphere. The TF allowed one to derive the primary He-nuclei fluxes in the same geomagnetic regions of AMS-01 observations. These fluxes were found in good agreement with the observation data. Furthermore inside the magnetosphere in addition to the flux of helium, it allowed one to obtain those of the primary p, C, and Fe nuclei and the relative abundances of He, C and Fe nuclei to protons from the same observation data of AMS-01 and HEAO-3 above ≈0.8 GeV/nucleon. Up to a geomagnetic latitude of ≈45.84°, the relative isotopic abundances were found to depend on the mass number I_isot and, on average, range from a factor ≈2.31 up to ≈3.35 larger than those outside the magnetosphere at 1 AU. Thus, the magnetospheric isotopic/nuclear relative abundances differ from those inside the solar cavity and those in the interstellar space. The usage of the TF approach can allow one to determine the nuclear abundances in the magnetosphere at any geomagnetic latitude and, thus, any orbit, provided that the CR spectra are determined at 1 AU.     

Observations with the SMM gamma-ray spectrometer

The Gamma Ray Spectrometer on the SMM satellite has observed solar cosmic energetic photon transients since 17 February 1980. Using the data available through 1981, new results have been obtained on ion acceleration phenomena in solar flares. It now is evident that both ion and electron acceleration can take place impulsively, simultaneously or within seconds of one another. That the impulsive acceleration process can produce ions with energies as high as GeV/nucleon is directly shown by observations of neutrons at the Earth with energies of several hundred MeV. These two facts and the relative timing of hard X-ray emissions provide new constraints on solar flare particle acceleration theory. New flare spectra have also been observed showing new nuclear γ-ray lines not previously observed from ²⁴Mg, ²⁰Ne and ⁵⁶Fe as well as from other elements. These spectral observations provide new information on the relative abundances of the accelerated and target nuclei. Following a review of the solar data and implications for flare theories we will also give a brief review of the results obtained on nonsolar γ-ray bursts. Most such bursts have photon spectra extending to MeV energies but with little, if any, evidence for spectral features.     

Normalized ionization yield function for various nuclei obtained with full Monte Carlo simulations

Several recent results important for production of ion pairs in the Earth atmosphere by various primary cosmic ray nuclei are presented. The direct ionization by various primary cosmic ray nuclei is explicitly obtained. The longitudinal profile of atmospheric cascades is sensitive to the energy and mass (charge) of the primary particle. In this study different cosmic ray nuclei are considered as primaries, namely Helium, Oxygen and Iron nuclei. The cosmic ray induced ionization is obtained on the basis of CORSIKA 6.52 code simulations using FLUKA 2006 and QGSJET II hadronic interaction models. The energy of the primary particles is normalized to GeV per nucleon. In addition, the ionization yield function Y is normalized as ion pair production per nucleon. The obtained ionization yield functions Y for various primaries are compared. The presented results and their application are discussed.     

Effects of solar modulation on the cosmic ray positron fraction

We implemented a 2D Monte Carlo model to simulate the solar modulation of galactic cosmic rays. The model is based on the Parker’s transport equation which contains diffusion, convection, particle drift and energy loss. Following the evolution in time of the solar activity, we are able to modulate a local interstellar spectrum (LIS), that we assumed isotropic beyond the termination shock, down to the Earth position inside the heliosphere. In this work we focused our attention to the cosmic ray positron fraction at energy below ∼10 GeV, showing how the particle drift processes could explain different results for AMS-01 and PAMELA. We compare our modulated spectra with observations at Earth, and then make a prediction of the cosmic ray positron fraction for the AMS-02 experiment.     

The extension of the INCL model for simulation of shielding in space

Radiation hazard for space missions is mainly due to cosmic ray protons, helium nuclei and light ions, whose energy spectrum is maximum around 1 GeV per nucleon but remains non-negligible for energies up to 15 GeV per nucleon. Nuclear reactions induced by high energy protons are often described by intranuclear cascade plus evaporation models. The attention is focused here on the Liège Intranuclear Cascade model (INCL), which has been shown to reproduce fairly well a great deal of experimental data for nucleon-induced reactions in the 200 MeV to 2 GeV range, when coupled with the ABLA evaporation-fission code. In order to extend the model to other conditions relevant for space radiation, three improvements of INCL are under development. They are reported on here. First, the reaction model has been extended to nucleon–nucleus reactions at incident energies up to 15 GeV, mainly by the inclusion of additional pion production channels in nucleon–nucleon collisions during the cascade. Second, a coalescence mechanism for the emission of light charged particles has been implemented recently. Finally, the model has been modified in order to accommodate light ions as projectiles. First results are shown and compared with illustrative experimental data. Implications for issues concerning radiation protection in space are discussed.     

Effect of random nature of cosmic ray sources – Supernova remnants – on cosmic ray intensity fluctuations,anisotropy, and electron energy spectrum

The random nature of sources (the supernova remnants) leads to the fluctuations of cosmic ray intensity in space and time. We calculate the expected fluctuations in a flat-halo diffusion model for particles with energies from 0.1 to 10³ TeV. The data on energy spectra and anisotropy of very high energy protons, nuclei and electrons, and the astronomical data on supernova remnants, the potential sources of cosmic rays, are used to constrain the value of the cosmic-ray diffusion coefficient and its dependence on energy.     

Simulation experiments of the cosmic ray isotope telescope on ETS VI

A computer simulation was carried out to evaluate the basic characteristics of a Δ E×E cosmic ray telescope consisting of 23 solid state detectors including 3 position sensitive detectors with large effective area. Based on the simulation, the geometric factor of the telescope is deduced to be as large as 22.5 cm²sr, almost independent of charge and energy concerned. The energy ranges to be covered by the telescope are, for example, 18–98 MeV/n for Li and 56–339 MeV/n for Fe. By analyzing simulated data, the mass resolution for iron in the overall energy range covered by the telescope is estimated as about 0.22 amu in standard deviation. The expected counting rates and mass-histograms are simulated for Galactic cosmic rays and solar energetic particles.     

Research by the Tasmanian cosmic ray group during the International Geophysical Year

Systematic recording of the cosmic radiation commenced in Hobart in 1946 and at Mawson in Antarctica in 1955, making these two of the longest running cosmic ray observatories in the world. For the IGY, observations were also made at a sub-Antarctic island and near the equator, and an airborne survey of the nucleonic component was made from Geomagnetic Latitude −60°, south of Australia, to Japan and back. At Hobart there were neutron monitors, vertical and inclined muon telescopes, an ionization chamber, and two muon telescopes at ∼40 m of water equivalent underground. The research based on these and other observations determined the energy dependence of the Forbush and 11-year variations and concentrated, in particular, on understanding the anisotropic nature of galactic cosmic rays up to 150 GeV; the anisotropies in the onset phase of Forbush decreases; and the anisotropies in solar cosmic ray events. An investigation was initiated to calculate the trajectories and cutoff rigidities of cosmic rays in a high order simulation of the geomagnetic field. This was completed in 1959–60.     

Modeling the acceleration and modulation of anomalous cosmic ray oxygen

After the solar wind termination shock crossings of the Voyager spacecraft, the acceleration of anomalous cosmic rays has become a very contentious subject. In this paper we examine several topics pertinent to anomalous cosmic ray oxygen acceleration and transport using a numerical cosmic ray modulation model. These include the effects of drifts on a purely Fermi I accelerated spectra, the effects of introducing higher charge states of oxygen into the modulation model, examining the viability of momentum diffusion as a re-acceleration process in the heliosheath and examining energy spectra, and intensity gradients, in the inner heliosphere during consecutive drift cycles.     

Recent results form measurements of the energy spectrum of cosmic-ray induced neutrons aboard an ER-2 airplane and on the ground.

Crews of future high-altitude commercial aircraft may be significantly exposed to atmospheric cosmic radiation from galactic cosmic rays (GCR). To help determine such exposures, the Atmospheric Ionizing Radiation Project, an international collaboration of 15 laboratories, made simultaneous radiation measurements with 14 instruments on a NASA ER-2 high-altitude aircraft. The primary instrument was a sensitive extended-energy multisphere neutron spectrometer, which was also used to make measurements on the ground. Its detector responses were calculated for neutrons and charged hadrons at energies up to 100 GeV using the radiation transport code MCNPX. We have now recalculated the detector responses including the effects of the airplane structure. We are also using new FLUKA calculations of GCR-induced hadron spectra in the atmosphere to correct for spectrometer counts produced by charged hadrons. Neutron spectra are unfolded from the corrected measured count rates using the MAXED code. Results for the measured cosmic-ray neutron spectrum (thermal to >10 GeV), total neutron fluence rate, and neutron dose equivalent and effective dose rates, and their dependence on altitude and geomagnetic cutoff generally agree well with results from recent calculations of GCR-induced neutron spectra.     

Time-dependent cosmic ray modulation

Time-dependent cosmic ray modulation is calculated over multiple solar cycles using our well established two-dimensional time-dependent modulation model. Results are compared to Voyager 1, Ulysses and IMP cosmic ray observations to establish compatibility. A time-dependence in the diffusion and drift coefficients, implicitly contained in recent expressions derived by , ,  and , is incorporated into the cosmic ray modulation model. This results in calculations which are compatible with spacecraft observations on a global scale over consecutive solar cycles. This approach compares well to the successful compound approach of Ferreira and Potgieter (2004). For both these approaches the magnetic field magnitude, variance of the field and current sheet tilt angle values observed at Earth are transported time-dependently into the outer heliosphere. However, when results are compared to observations for extreme solar maximum, the computed step-like modulation is not as pronounced as observed. This indicates that some additional merging of these structures into more pronounced modulation barriers along the way is needed.     

Trapped antiprotons produced by cosmic rays in the Earth's magnetosphere.

The existence of significant fluxes of antiparticles in the Earth magnetosphere has been predicted on theoretical considerations in this article. These antiparticles (positrons or antiprotons) at several hundred kilometers of altitudes, we believe are not of direct extraterrestrial origin, but are the natural products of nuclear reactions of the high energy primary cosmic rays (CR) and trapped protons (TP) confined in the terrestrial radiation belt, with the constituents of terrestrial atmosphere. Extraterrestrial positrons and antiprotons born in nuclear reactions of the same CR particles passing through only 5-7 g/cm2 of interstellar matter, exhibit lower fluxes compared to the antiprotons born at hundreds of g/cm2 in the atmosphere, which when confined in the magnetic field of the Earth (in any other planet), get accumulated. We present the results of the computations of the antiproton fluxes at 10 MeV to several GeV energies due to CR particle interactions with the matter in the interstellar space, and also with the residual atmosphere at altitudes of approximately 1000 km over the Earth's surface. The estimates show that the magnetospheric antiproton fluxes are greater by two orders of magnitude compared to the extraterrestrial fluxes measured at energies <1-2 GeV.     

A search for the signature of microquasars in the cosmic ray iron spectrum measured by TIGER

In a recent paper Heinz and Sunyaev suggest that relativistic jets observed in microquasars might result in narrow features in the energy spectra of heavy cosmic rays with ≈1 to ≈10 GeV/nuc. They further argue that such features might be observable if there has been one or more microquasars nearby within the last few million years. We report preliminary results of a search for evidence of such features using data from a 32-day balloon flight of the Trans-Iron Galactic Element Recorder (TIGER). Although this flight took place near solar maximum, calculations of the broadening effects of solar modulation indicate that a narrow feature of sufficient intensity should still be observable. An energy spectrum for iron with high statistical significance has been derived from ≈100,000 Fe events in the energy range from about 2.5 to 10 GeV/nuc. Although our preliminary results do not reveal any obvious features, we will discuss the possibility of observing such features with TIGER and other instruments.     

A model of galactic cosmic rays for use in calculating linear energy transfer spectra.

The galactic cosmic rays (GCR) contain fully stripped nuclei, from Hydrogen to beyond the Iron group, accelerated to high energies and are a major component of the background radiation encountered by satellites and interplanetary spacecraft. This paper presents a GCR model which is based upon our current understanding of the astrophysics of GCR transport through interstellar and interplanetary space. The model can be used to predict the energy spectra for all stable and long-lived radioactive species from H to Ni over an energy range from 50 to 50,000 MeV/nucleon as a function of a single parameter, the solar modulation level phi. The details of this model are summarized, phi is derived for the period 1974 to present, and results from this model during the 1990/1991 CRRES mission are presented.