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.     

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.     

Ground-based simulations of galactic cosmic ray fragmentation and transport.

Since mean free paths for nuclear fragmentation are of the order of the ranges of primary Galactic Cosmic Ray (GCR) nuclei, determination of the radiation field produced by successive fragmentations of nuclei in material and tissue is essential to accurate assessment of GCR radiation risk to humans on long-duration missions outside the geomagnetosphere. We describe some recent measurements made at the Bevalac of heavy ion transport through materials, with representative results and examples of how they may be applied to aspects of the space radiation problem, including efforts to devise analytical tools for predicting biological effects and for designing spacecraft shielding.     

Analytical-HZETRN model for rapid assessment of active magnetic radiation shielding

The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA’s radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than ∼15 Tm and passive/structural shielding thicknesses greater than 40 g/cm² have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.     

Estimation of Galactic Cosmic Ray exposure inside and outside the Earth’s magnetosphere during the recent solar minimum between solar cycles 23 and 24

The evidently low solar activity observed between solar cycles 23 and 24 during the years 2008–2010 led to a substantial increase in the Galactic Cosmic Ray (GCR) intensity in comparison with preceding solar minima. As the GCRs consist of highly-ionizing charged particles having the potential to cause biological damage, they are a subject of concern for manned missions to space. With the enhanced particle fluxes observed between 2008 and 2010, it is reasonable to assume that the radiation exposure from GCR must have also increased to unusually high levels. In this paper, the GCR exposure outside and inside the Earth’s magnetosphere is numerically calculated for time periods starting from 1970 to the end of 2011 in order to investigate the increase in dose levels during the years 2008–2010 in comparison with the last three solar minima. The dose rates were calculated in a water sphere, used as a surrogate for the human body, either unshielded or surrounded by aluminium shielding of 0.3, 10 or 40 g/cm².     

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.     

Solar modulation and nuclear fragmentation effects in galactic cosmic ray transport through shielding.

Crews of manned interplanetary missions may accumulate significant radiation exposures from the galactic cosmic ray (GCR) environment in space. Estimates of how these dose levels are affected by the assumed temporal and spatial variations in the composition of the GCR environment, and by the effects of the spacecraft and body self-shielding on the transported radiation fields are presented. In this work, the physical processes through which shielding alters the transported radiation fields are described. We then present estimates of the effects on model calculations of (1) nuclear fragmentation model uncertainties, (2) solar modulation, (3) variations between solar cycles, and (4) proposed changes to the quality factors which relate dose equivalent to absorbed dose.     

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 space radiation shielding model of the Martian radiation environment experiment (MARIE)

The 2001 Mars Odyssey spacecraft was launched towards Mars on April 7, 2001. Onboard the spacecraft is the Martian radiation environment experiment (MARIE), which is designed to measure the background radiation environment due to galactic cosmic rays (GCR) and solar protons in the 20–500 MeV/n energy range. We present an approach for developing a space radiation-shielding model of the spacecraft that includes the MARIE instrument in the current mapping phase orientation. A discussion is presented describing the development and methodology used to construct the shielding model. For a given GCR model environment, using the current MARIE shielding model and the high-energy particle transport codes, dose rate values are compared with MARIE measurements during the early mapping phase in Mars orbit. The results show good agreement between the model calculations and the MARIE measurements as presented for the March 2002 dataset.     

Pion and electromagnetic contribution to dose: Comparisons of HZETRN to Monte Carlo results and ISS data

Recent work has indicated that pion production and the associated electromagnetic (EM) cascade may be an important contribution to the total astronaut exposure in space. Recent extensions to the deterministic space radiation transport code, HZETRN, allow the production and transport of pions, muons, electrons, positrons, and photons. In this paper, the extended code is compared to the Monte Carlo codes, Geant4, PHITS, and FLUKA, in slab geometries exposed to galactic cosmic ray (GCR) boundary conditions. While improvements in the HZETRN transport formalism for the new particles are needed, it is shown that reasonable agreement on dose is found at larger shielding thicknesses commonly found on the International Space Station (ISS). Finally, the extended code is compared to ISS data on a minute-by-minute basis over a seven day period in 2001. The impact of pion/EM production on exposure estimates and validation results is clearly shown. The Badhwar–O’Neill (BO) 2004 and 2010 models are used to generate the GCR boundary condition at each time-step allowing the impact of environmental model improvements on validation results to be quantified as well. It is found that the updated BO2010 model noticeably reduces overall exposure estimates from the BO2004 model, and the additional production mechanisms in HZETRN provide some compensation. It is shown that the overestimates provided by the BO2004 GCR model in previous validation studies led to deflated uncertainty estimates for environmental, physics, and transport models, and allowed an important physical interaction (π/EM) to be overlooked in model development. Despite the additional π/EM production mechanisms in HZETRN, a systematic under-prediction of total dose is observed in comparison to Monte Carlo results and measured data.     

航天员受银河宇宙线辐射的剂量计算

在近地空间(LEO)和深空探测中,航天员遭受的辐射风险主要来自于银河宇宙线(GCR)照射.银河宇宙线的辐射剂量是航天员辐射风险评价的基础.国际放射防护委员会(ICRP)于2013年提出了新的航天员空间辐射剂量估算方法,以更准确给出空间重离子辐射的剂量.基于此方法,开发了宇宙线粒子在物质中输运的蒙特卡罗程序,并在程序中实现用中国成年男性人体数字模型来仿真航天员.采用该程序计算了粒子(Z=1~92)各向同性照射航天员时器官的通量-器官剂量转换因数,并估算出航天员在近地轨道空间受银河宇宙线辐射的剂量.     

Energetic particle environment in near-Earth orbit.

The hazard of exposure to high doses of ionizing radiation is one of the primary concerns of extended manned space missions and a continuous threat for the numerous spacecraft in operation today. In the near-Earth environment the main sources of radiation are solar energetic particles (SEP), galactic cosmic rays (GCR), and geomagnetically trapped particles, predominantly protons and electrons. The intensity of the SEP and GCR source depends primarily on the phase of the solar cycle. Due to the shielding effect of the Earth's magnetic field, the observed intensity of SEP and GCR particles in a near-Earth orbit will also depend on the orbital parameters altitude and inclination. The magnetospheric source strength depends also on these orbital parameters because they determine the frequency and location of radiation belt passes. In this paper an overview of the various sources of radiation in the near-Earth orbit will be given and first results obtained with the Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX) will be discussed. SAMPEX was launched on 3 July 1992 into a near polar (inclination 82 degrees) low altitude (510 x 675 km) orbit. The SAMPEX payload contains four separate instruments of high sensitivity covering the energy range 0.5 to several hundred MeV/nucleon for ions and 0.4 to 30 MeV for electrons. This low altitude polar orbit with zenith-oriented instrumentation provides a new opportunity for a systematic study of the near-Earth energetic particle environment.     

Fragmentation of 1 GeV/nucleon iron ions in thick targets relevant for space exploration.

We have measured charged nuclear fragments produced by 1 GeV/nucleon 56Fe ions interacting with aluminium, polyethylene and lead. These materials are relevant for assessment of radiation risk for manned space flight. The data will be presented in a form suitable for comparison with models of nuclear fragmentation and transport, including linear energy transfer (LET) spectrum, fluence for iron and fragments, event-tack- and event-dose-averaged LET, total dose and iron contribution to dose.     

A model for predicting the radiation exposure for mission planning aboard the international space station

The International Space Station Cosmic Radiation Exposure Model (ISSCREM) has been developed as a possible tool for use in radiation mission planning as based on operational data collected with a tissue equivalent proportional counter (TEPC) aboard the ISS since 2000. It is able to reproduce the observed trapped radiation and galactic cosmic radiation (GCR) contributions to the total dose equivalent to within ±20% and ±10%, respectively, as would be measured by the onboard TEPC at the Zvezda Service Module panel 327 (SM-327). Furthermore, when these contributions are combined, the total dose equivalent that would be measured at this location is estimated to within ±10%. The models incorporated into ISSCREM correlate the GCR dose equivalent rate to the cutoff rigidity magnetic shielding parameter and the trapped radiation dose equivalent rate to atmospheric density inside the South Atlantic Anomaly. The GCR dose equivalent rate is found to vary minimally with altitude and TEPC module location however, due to the statistics and data available, the trapped radiation model could only be developed for the TEPC located at SM-327. Evidence of the variation in trapped radiation dose with detector orientation and the East–West asymmetry were observed at this location.     

Neutron yields from interactions of GCR-like beams in stopping targets.

In order to accurately determine the radiation risk to astronauts from GCR, the nature of the secondary radiation field created by the fragmentation of GCR in shielding and tissue must be understood. Due to the their high penetrabilities, neutrons are an important component of the secondary radiation field, especially for astronauts protected by thick shielding on lunar or Martian bases. Neutron yields from 435A MeV and 272A MeV Nb stopping in Nb and Al targets are presented, along with some preliminary analysis of neutron yields from 155A MeV C stopping in Al. Energy spectra and angular distributions are shown for neutron energies above 20 MeV. The data provides some information about the dependence of the neutron yield on projectile energy and target mass. Comparisons of the data with BUU calculations are also shown.     

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.     

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.     

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.