A comparison of total reaction cross section models used in particle and heavy ion transport codes

Understanding the interactions and propagations of high energy protons and heavy ions are essential when trying to estimate the biological effects of Galactic Cosmic Rays (GCR) and Solar Particle Events (SPE) on personnel in space. To be able to calculate the shielding properties of different materials and radiation risks, particle and heavy ion transport codes are needed. In all particle and heavy ion transport codes, the probability function that a projectile particle will collide within a certain distance x in the matter depends on the total reaction cross sections, and the calculated partial fragmentation cross sections scale with the total reaction cross sections. It is therefore crucial that accurate total reaction cross section models are used in the transport calculations. In this paper, different models for calculating nucleon–nucleus and nucleus–nucleus total reaction cross sections are compared with each other and with measurements. The uncertainties in the calculations with the different models are discussed, as well as their overall performances with respect to the available experimental data. Finally, a new compilation of experimental data is briefly presented.     

Theoretical model of HZE particle fragmentation by hydrogen targets.

The fragmenting of high energy, heavy ions (HZE particles) by hydrogen targets is an important, physical process in several areas of space radiation research. In this work quantum mechanical optical model methods for estimating cross sections for HZE particle fragmentation by hydrogen targets are presented. The cross sections are calculated using a modified abrasion-ablation collision formalism adapted from a nucleus-nucleus collision model. Elemental and isotopic production cross sections are estimated and compared with report measurements for the breakup of neon, sulphur, and iron, nuclei at incident energies between 400 and 910 MeV/nucleon. Good agreement between theory and experiment is obtained.     

An update about recent developments of the PHITS code

PHITS (Particle and Heavy-Ion Transport code System) is a general-purpose three-dimensional Monte Carlo code, developed and maintained by RIST, JAEA and KEK in Japan together with Sihver et al. at Chalmers in Sweden. PHITS can deal with the transports of all varieties of hadrons and heavy ions with energies up to around 100 GeV/nucleon, and in this paper the current status of PHITS is presented. We introduce a relativistically covariant version of JQMD, called R-JQMD, that features an improved ground state initialization algorithm, and we will present the introduction of electron and photon transport in PHITS using EGS5, which have increased the energy region for the photon and energy transport from up to around 3 GeV to up to several hundred GeV depending on the atomic number of the target. We show how the accuracy in dose and fluence calculations can be improved by using tabulated cross sections. Benchmarking of shielding and irradiation effects of high energy protons in different materials relevant for shielding of accelerator facilities is also presented. In particular, we show that PHITS can be used for estimating the dose received by aircrews and personnel in space. In recent years, many countries have issued regulations or recommendations to set annual dose limitations for aircrews. Since estimation of cosmic-ray spectra in the atmosphere is an essential issue for the evaluation of aviation doses, we have calculated these spectra using PHITS. The accuracy of the atmospheric propagation simulation of cosmic-ray performed by PHITS has been well verified by experimental cosmic-ray spectra taken under various conditions. Based on a comprehensive analysis of the simulation results, an analytical model called “PARMA” has been proposed for instantaneously estimating the atmospheric cosmic-ray spectra below the altitude of 20 km. We have also performed preliminary simulations of long-term dose distribution measurements at the ISS performed with the joint ESA-FSA experiment MATROSHKA-R (MTR-R) led by the Russian Federation Institute of Biomedical Problems (IMBP) and the ESA supported experiment MATROSHKA (MTR), led by the German Aerospace Center (DLR). For the purpose of examining the applicability of PHITS to the shielding design in space, the absorbed doses in a tissue equivalent water phantom inside an imaginary space vessel has been estimated for different shielding materials of different thicknesses. The results confirm previous results which indicate that PHITS is a suitable tool when performing shielding design studies of spacecrafts.     

Depth dependency of neutron density produced by cosmic rays in the lunar subsurface

Depth dependency of neutrons produced by cosmic rays (CRs) in the lunar subsurface was estimated using the three-dimensional Monte Carlo particle and heavy ion transport simulation code, PHITS, incorporating the latest high energy nuclear data, JENDL/HE-2007. The PHITS simulations of equilibrium neutron density profiles in the lunar subsurface were compared with the measurement by Apollo 17 Lunar Neutron Probe Experiment (LNPE). Our calculations reproduced the LNPE data except for the 350–400 mg/cm² region under the improved condition using the CR spectra model based on the latest observations, well-tested nuclear interaction models with systematic cross section data, and JENDL/HE-2007.     

Nuclear fragmentation database for GCR transport code development

A critical need for NASA is the ability to accurately model the transport of heavy ions in the Galactic Cosmic Rays (GCR) through matter, including spacecraft walls, equipment racks, etc. Nuclear interactions are of great importance in the GCR transport problem, as they can cause fragmentation of the incoming ion into lighter ions. Since the radiation dose delivered by a particle is proportional to the square of (charge/velocity), fragmentation reduces the dose delivered by incident ions. The other mechanism by which dose can be reduced is ionization energy loss, which can lead to some particles stopping in the shielding. This is the conventional notion of shielding, but it is not applicable to human spaceflight since the particles in the GCR tend to be too energetic to be stopped in the relatively thin shielding that is possible within payload mass constraints. Our group has measured a large number of fragmentation cross sections, intended to be used as input to, or for validation of, NASA’s radiation transport models. A database containing over 200 charge-changing cross sections and over 2000 fragment production cross sections has been compiled. In this report, we examine in detail the contrast between fragment measurements at large acceptance and small acceptance. We use output from the PHITS Monte Carlo code to test our assumptions using as an example ⁴⁰Ar data (and simulated data) at a beam energy of 650 MeV/nucleon. We also present preliminary analysis in which isotopic resolution was attained for beryllium fragments produced by beams of ¹⁰B and ¹¹B. Future work on the experimental data set will focus on extracting and interpreting production cross sections for light fragments.     

Estimates of HZE particle contributions to SPE radiation exposures on interplanetary missions.

Estimates of radiation doses resulting from possible HZE (high energy heavy ion) components of solar particle events (SPEs) are presented for crews of manned interplanetary missions. The calculations assume a model spectrum obtained by folding measured solar flare HZE particle abundances with the measured energy spectra of SPE alpha particles. These hypothetical spectra are then transported through aluminum spacecraft shielding. The results, presented as estimates of absorbed dose and dose equivalent, indicate that HZE components by themselves are not a major concern for crew protection but should be included in any overall risk assessment. The predictions are found to be sensitive to the assumed spectral hardness parameters.     

A program to measure new energetic particle nuclear interaction cross sections.

The Transport Collaboration, consisting of researchers from institutions in France, Germany, Italy and the USA, has established a program to make new measurements of nuclear interaction cross sections for heavy projectiles (Z > or = 2) in targets of liquid H2, He and heavier materials. Such cross sections directly affect calculations of galactic and solar cosmic ray transport through matter and are needed for accurate radiation hazard assessment. To date, the collaboration has obtained data using the LBL Bevalac HISS facility with 20 projectiles from 4He to 58Ni in the energy range 393-910 MeV/nucleon. Preliminary results from the analysis of these data are presented here and compared to other measurements and to cross section prediction formulae.     

Early and late damages induced by heavy charged particle irradiation in embryonic tissue of Arabidopsis seeds.

Early and late effects of accelerated heavy ions (HZE) on the embryonic tissue of Arabidopsis thaliana seeds were investigated seeing that initial cells of the plant eumeristems resemble the original cells of animal and human tissues with continuous cell proliferation. The endpoints measured were lethality and tumorization in the M1-generation for early effects and embryonic lethality in the M2-generation for late effects. The biological endpoints are plotted as functions of the physical parameters of the irradiation i.e. ion fluence (p/cm2), dose (Gray), charge Z and linear energy transfer (LET). The results presented contribute to the estimation of the principles of biological HZE effects and thus may help to develop a unified theory which could explain the whole sequence from physical and chemical reactions to biological responses connected with heavy ion radiation. Additionally, the data of this paper may be used for the discussion of the quality factor for heavy ion irradiation needed for space missions and for HZE-application in radio-therapy by use of accelerators (UNILAC, (SIS/ESR), BEVALAC).     

PHITS simulations of the Matroshka experiment

The radiation environment in space is very different from the one encountered on Earth. In addition to the sparsely ionizing radiation, there are particles of different Z with energies ranging from keV up to hundreds of GeV which can cause severe damage to both electronics and humans. It is therefore important to understand the interactions of these highly ionizing particles with different materials such as the hull of space vehicles, human organs and electronics. We have used the Particle and Heavy-Ion Transport code System (PHITS), which is a three-dimensional Monte Carlo code able to calculate interactions and transport of particles and heavy ions with energies up to 100 GeV/nucleon in most matter. PHITS is developed and maintained by a collaboration between RIST (Research Organization for Information Science & Technology), JAEA (Japan Atomic Energy Agency), KEK (High Energy Accelerator Research Organization), Japan and Chalmers University of Technology, Sweden. For the purpose of examining the applicability of PHITS to the shielding design we have simulated the ESA facility Matroshka (MTR) designed and lead by the German Aerospace Center (DLR). Preliminary results are presented and discussed in this paper.     

Results on Artemia cysts, lettuce and tobacco seeds in the Biobloc 4 experiment flown aboard the Soviet Biosatellite Cosmos 1129.

Artemia cysts, lettuce and tobacco seeds were flown aboard the Cosmos 1129 for 19 days. A correlative method was used in order to determine the passage of cosmic heavy ions (HZE particles) through the biological test objects. This space flight resulted in a decrease on hatchability, nucleic acid and protein synthesis in hydrated Artemia cysts. HZE particle effects on plant cellular chromosomes are confirmed. In tobacco seeds, a stimulating effect on germination rate and a higher frequency of abnormalities were observed. Dormant biological objects are a very suitable material to study cosmic ray effects: these objects can be arranged in monolayers and sandwiched between visual track detectors in order to determine the passage of the cosmic heavy ions (HZE particles). On the other hand this method allows us to study effects of microgravity and those of the protonic component of cosmic rays in the objects not hit by the HZE articles.     

Overview of nuclear fragmentation models and needs.

It has been known for some time that adequate assessment of spacecraft shield requirements and concomitant estimates of astronauts radiation exposures from galactic cosmic radiation requires accurate, quantitative methods for characterizing these radiation fields as they pass through thick absorbers. The main nuclear interaction processes involved are (1) nuclear elastic and inelastic collisions, and (2) nuclear breakup (fragmentation) and electromagnetic dissociation (EMD). Nuclear fragmentation and EMD are important because they alter the elemental and isotopic composition of the transported radiation fields. At present, there is no suitably accurate theory for predicting nuclear fragmentation cross sections for all collision pairs and energies of interest in space radiation protection. Typical cross-section differences between theory and experiment range from about 25 percent to a factor of two. The resulting errors in transported flux, for high linear energy transfer (LET) particles, are comparble to these cross-section errors. In this overview, theoretical models of heavy ion fragmentation currently used to generate input data bases for cosmic-ray transport and shielding codes are reviewed. Their shortcomings are discussed. Further actions needed to improve their accuracy and generality are presented.     

Assessment and requirements of nuclear reaction databases for GCR transport in the atmosphere and structures

The transport properties of galactic cosmic rays (GCR) in the atmosphere, material structures, and human body (self-shielding) are of interest in risk assessment for supersonic and subsonic aircraft and for space travel in low-Earth orbit and on interplanetary missions. Nuclear reactions, such as knockout and fragmentation, present large modifications of particle type and energies of the galactic cosmic rays in penetrating materials. We make an assessment of the current nuclear reaction models and improvements in these model for developing required transport code data bases. A new fragmentation data base (QMSFRG) based on microscopic models is compared to the NUCFRG2 model and implications for shield assessment made using the HZETRN radiation transport code. For deep penetration problems, the build-up of light particles, such as nucleons, light clusters and mesons from nuclear reactions in conjunction with the absorption of the heavy ions, leads to the dominance of the charge Z = 0, 1, and 2 hadrons in the exposures at large penetration depths. Light particles are produced through nuclear or cluster knockout and in evaporation events with characteristically distinct spectra which play unique roles in the build-up of secondary radiation's in shielding. We describe models of light particle production in nucleon and heavy ion induced reactions and make an assessment of the importance of light particle multiplicity and spectral parameters in these exposures.     

Ground-based measurements of galactic cosmic ray fragmentation in shielding.

The mean free path for nuclear interactions of galactic cosmic-rays is comparable to shielding and tissue thicknesses present in human interplanetary exploration, resulting in a significant fraction of nuclear reaction products at depth. In order to characterize the radiation field, the energy spectrum, the angular distribution, and the multiplicity of each type of secondary particles must also be known as a function of depth. Reactions can take place anywhere in a thick absorber; therefore, it is necessary to know these quantities as a function of particle energy for all particles produced. HZE transport methods are used to predict the radiation field; they are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The most cost-effective method of validation is a comparison with ground-based experimental measurements. A research program to provide such validation measurements using neon, iron and other accelerated heavy ion beams will be discussed and illustrated using results from ongoing experiments and their comparison with current transport codes. The extent to which physical measurements yield radiobiological predictions will be discussed.     

In vitro neurotoxic effects of 1 GeV/n iron particles assessed in retinal explants.

The heavy ion component of the cosmic radiation remains problematic to the assessment of risk in manned space flight. The biological effectiveness of HZE particles has yet to be established, particularly with regard to nervous tissue. Using heavy ions accelerated at the AGS of Brookhaven National Laboratory, we study the neurotoxic effects of iron particles. We exposed retinal explants, taken from chick embryos, to determine the dose response relationships for neurite outgrowth. Morphometric techniques were used to evaluate the in vitro effects of 1 GeV/a iron particles (LET 148 keV/micrometer). Iron particles produced a dose-dependent reduction of neurite outgrowth with a maximal effect achieved with a dose of 100 cGy. Doses as low as 10-50 cGy were able to induce reductions of the neurite outgrowth as compared to the control group. Neurite generation is a more sensitive parameter than neurite elongation, suggesting different mechanism of radiation damage in our model. These results showed that low doses/fluences of iron particles could impair the retinal ganglion cells' capacity to generate neurites indicating the highly neurotoxic capability of this heavy charged particle.     

Microlesions: theory and reality.

Efforts to assess radiation risk in space have been complicated by the considerable unknowns regarding the biological effects of the heavy ion component (HZE particles) of the cosmic rays. The attention has focused primarily on the assignation of a quality factor (Q) which would take into account the greater effectiveness of heavy ions vis-a-vis other forms of ionizing radiation. If however, as the so-called "Microlesion Theory" allows, the passage of HZE particles through living tissue produces unique biological damage, the traditional use of Q becomes meaningless. Therefore, it is critical to determine if microlesions, in fact, do exist. While the concept does not necessarily require detectable morphological damage, "tunnel-lesions" or holes in ocular tissues have been cited as evidence of microlesions. These data, however, are open to reinterpretation. On-going light, scanning and transmission electron microscopic studies of the corneas, lenses and retinas of rat eyes exposed to 450 MeV/amu 56Fe ions thus far have not revealed tunnel-lesion damage. The morphological effects of the heavy ions have been found to be qualitatively similar to the changes following other kinds of ionizing radiation.     

Liquid drop model considerations in HZE particle fragmentation by hydrogen

The fragmenting of high energy, heavy ions by hydrogen targets is an important physical process in several areas of space radiation protection research. Quantum mechanical, optical model methods for calculating cross sections for particle fragmentation by hydrogen have been developed from a modified abrasion-ablation collision formalism. The abrasion stage is treated as a knockout process which leaves the residual prefragment in an excited state. In the ablation stage the prefragment deexcites to produce the final fragment. The prefragment excitation energies are estimated from a combination of liquid drop model considerations and frictional-spectator interaction processes. Estimates of elemental and isotopic production cross sections are in good agreement with published cross section measurements.     

An improved Green's function for ion beam transport.

Ion beam transport theory allows testing of material transmission properties in the laboratory environment generated by particle accelerators. This is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are the main emphasis of the present work. In consequence, an analytic solution of the linear Boltzmann equation is pursued in the form of a Green's function allowing flexibility in application to a broad range of boundary value problems. It has been established that simple solutions can be found for high charge and energy (HZE) ions by ignoring nuclear energy downshifts and dispersion. Such solutions were found to be supported by experimental evidence with HZE ion beams when multiple scattering was added. Lacking from the prior solutions were range and energy straggling and energy downshift with dispersion associated with nuclear events. Recently, we have found global solutions including these effects providing a broader class of HZE ion solutions.