Simulations of an accelerator-based shielding experiment using the particle and heavy-ion transport code system PHITS.

In order to estimate the biological effects of HZE particles, an accurate knowledge of the physics of interaction of HZE particles is necessary. Since the heavy ion transport problem is a complex one, there is a need for both experimental and theoretical studies to develop accurate transport models. RIST and JAERI (Japan), GSI (Germany) and Chalmers (Sweden) are therefore currently developing and bench marking the General-Purpose Particle and Heavy-Ion Transport code System (PHITS), which is based on the NMTC and MCNP for nucleon/meson and neutron transport respectively, and the JAM hadron cascade model. PHITS uses JAERI Quantum Molecular Dynamics (JQMD) and the Generalized Evaporation Model (GEM) for calculations of fission and evaporation processes, a model developed at NASA Langley for calculation of total reaction cross sections, and the SPAR model for stopping power calculations. The future development of PHITS includes better parameterization in the JQMD model used for the nucleus-nucleus reactions, and improvement of the models used for calculating total reaction cross sections, and addition of routines for calculating elastic scattering of heavy ions, and inclusion of radioactivity and burn up processes. As a part of an extensive bench marking of PHITS, we have compared energy spectra of secondary neutrons created by reactions of HZE particles with different targets, with thicknesses ranging from <1 to 200 cm. We have also compared simulated and measured spatial, fluence and depth-dose distributions from different high energy heavy ion reactions. In this paper, we report simulations of an accelerator-based shielding experiment, in which a beam of 1 GeV/n Fe-ions has passed through thin slabs of polyethylene, Al, and Pb at an acceptance angle up to 4 degrees.     

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.     

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.     

The physics of the FLUKA code: Recent developments

FLUKA is a Monte-Carlo code able to simulate interaction and transport of hadrons, heavy ions and electromagnetic particles from few keV (or thermal neutron) to cosmic ray energies in whichever material. The highest priority in the design and development of the code has always been the implementation and improvement of sound and modern physical models. A summary of the FLUKA physical models is given, while recent developments are described in detail: among the others, extensions of the intermediate energy hadronic interaction generator, refinements in photon cross sections and interaction models, analytical on-line evolution of radio-activation and remnant dose. In particular, new developments in the nucleus–nucleus interaction models are discussed. Comparisons with experimental data and examples of applications of relevance for space radiation are also provided.     

Transport of light ions in matter

A recent set of light ion experiments are analyzed using the Green's function method of solving the Boltzmann equation for ions of high charge and energy (the GRNTRN transport code) and the NUCFRG2 fragmentation database generator code. Although the NUCFRG2 code reasonably represents the fragmentation of heavy ions, the effects of light ion fragmentation requires a more detailed nuclear model including shell structure and short range correlations appearing as tightly bound clusters in the light ion nucleus. The most recent NUCFRG2 code is augmented with a quasielastic alpha knockout model and semiempirical adjustments (up to 30 percent in charge removal) in the fragmentation process allowing reasonable agreement with the experiments to be obtained. A final resolution of the appropriate cross sections must await the full development of a coupled channel reaction model in which shell structure and clustering can be accurately evaluated.     

The FLUKA code for space applications: recent developments.

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions.     

Mutation induction in yeast by very heavy ions.

Resistance to canavanine was studied in haploid yeast after exposure to heavy ions (argon to uranium) of energies between 1 and 10 MeV/u covering a LET-range up to about 10000 keV/micrometer. Mutations were found in all instances but the induction cross sections increased with ion energy. This is taken to mean that the contribution of penumbra electrons plays an important role. The probability to recover surviving mutants is highest if the cell is not directly hit by the particle. The experiments demonstrate that the geometrical dimensions of the target cell nucleus as well as its sensitivity in terms of survival have a critical influence on mutation induction with very heavy ions.     

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.     

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.     

Comparisons of several transport models in their predictions in typical space radiation environments

We have used several transport codes to calculate dose and dose equivalent values as well as the particle spectra behind a slab or inside a spherical shell shielding in typical space radiation environments. Two deterministic codes, HZETRN and UPROP, and two Monte Carlo codes, FLUKA and Geant4, are included. A soft solar particle event, a hard solar particle event, and a solar minimum galactic cosmic rays environment are considered; and the shielding material is either aluminum or polyethylene. We find that the dose values and particle spectra from HZETRN are in general rather consistent with Geant4 except for neutrons. The dose equivalent values from HZETRN and Geant4 are not far from each other, but the HZETRN values behind shielding are often lower than the Geant4 values. Results from FLUKA and Geant4 are mostly consistent for considered cases. However, results from the legacy code UPROP are often quite different from the other transport codes, partly due to its non-consideration of neutrons. Comparisons for the spherical shell geometry exhibit the same qualitative features as for the slab geometry. In addition, results from both deterministic and Monte Carlo transport codes show that the dose equivalent inside the spherical shell decreases from the center to the inner surface and this decrease is large for solar particle events; consistent with an earlier study based on deterministic radiation transport results. This study demonstrates both the consistency and inconsistency among these transport models in their typical space radiation predictions; further studies will be required to pinpoint the exact physics modules in these models that cause the differences and thus may be improved.     

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.     

Comparison of the transport codes HZETRN,HETC and FLUKA for galactic cosmic rays

The HZETRN deterministic radiation code is one of several tools developed to analyze the effects of harmful galactic cosmic rays (GCR) and solar particle events on mission planning and shielding for astronauts and instrumentation. This paper is a comparison study involving the two Monte Carlo transport codes, HETC–HEDS and FLUKA and the deterministic transport code, HZETRN. Each code is used to transport an ion from the 1977 solar minimum GCR spectrum impinging upon a 20 g/cm² aluminum slab followed by a 30 g/cm² water slab. This research is part of a systematic effort of verification and validation to quantify the accuracy of HZETRN and determine areas where it can be improved. Comparisons of dose and dose equivalent values at various depths in the water slab are presented in this report. This is followed by a comparison of the proton and forward, backward and total neutron flux at various depths in the water slab. Comparisons of the secondary light ion ²H, ³H, ³He and ⁴He fluxes are also examined.     

Heavy ion induced DNA double strand breaks in cells of E. coli.

Vegetative cells of E. coli differing in their radiosensitivity have been used in heavy ion irradiation experiment. Besides inactivation measurements also the induction of DNA double strand breaks (DSB) have been measured using the method of pulse-field gel electrophoresis. This method allows to separate linear DNA with length up to 8 Mio base pairs. After irradiation with heavy ions we find a higher amount of low molecular weight fragments when compared to sparsely ionizing radiation. This agrees with the idea that heavy ions as a structured radiation have a high probability to induce more than one strand break in a DNA molecule if the particle hits the DNA. The amount of intact DNA remaining in the agarose plugs decreases exponentially for increasing radiation doses or particle fluences. From these curves cross sections for the induction of DSB after heavy ion irradiation have been determined. These results will be discussed in comparison to the results for cell survival.     

Deterministic pion and muon transport in Earth’s atmosphere

An accurate understanding of the physical interactions and transport of space radiation is important for safe and efficient space operations. Secondary particles produced by primary particle interactions with intervening materials are an important contribution to radiation risk. Pions are copiously produced in the nuclear interactions typical of space radiations and can therefore be an important contribution to radiation exposure. Charged pions decay almost exclusively to muons. As a consequence, muons must also be considered in space radiation exposure studies. In this work, the NASA space radiation transport code HZETRN has been extended to include the transport of charged pions and muons. The relevant transport equation, solution method, and implemented cross sections are reviewed. Muon production in the Earth’s upper atmosphere is then investigated, and comparisons with recent balloon flight measurements of differential muon flux are presented. Muon production from the updated version of HZETRN is found to match the experimental data well.     

Validation of the new trapped environment AE9/AP9/SPM at low Earth orbit

The completion of the international space station (ISS) in 2011 has provided the space research community an ideal proving ground for future long duration human activities in space. Ionizing radiation measurements in ISS form the ideal tool for the validation of radiation environmental models, nuclear transport codes and nuclear reaction cross sections. Indeed, prior measurements on the space transportation system (STS; shuttle) provided vital information impacting both the environmental models and the nuclear transport code developments by indicating the need for an improved dynamic model of the low Earth orbit (LEO) trapped environment. Additional studies using thermo-luminescent detector (TLD), tissue equivalent proportional counter (TEPC) area monitors, and computer aided design (CAD) model of earlier ISS configurations, confirmed STS observations that, as input, computational dosimetry requires an environmental model with dynamic and directional (anisotropic) behavior, as well as an accurate six degree of freedom (DOF) definition of the vehicle attitude and orientation along the orbit of ISS.