Comparison of space radiation GCR models to AMS heavy ion data.

Galactic cosmic rays (GCR) are a constant source of radiation that constitutes one of the major hazards during deep space exploration missions for both astronauts and hardware. In this work, GCR models commonly used by the space radiation protection community are compared with recently published high-precision, high-resolution measurements of cosmic ray lithium, beryllium, boron, carbon, nitrogen, and oxygen fluxes along with their ratios (Li/B, Li/C, Li/O, Be/B, Be/C, Be/O, B/C, B/O, C/O, N/B, N/O) from the Alpha Magnetic Spectrometer (AMS). All of the models were developed and calibrated prior to the publication of this AMS data, therefore this is an opportunity to validate the models against an independent data set. This paper is a compliment to the previously published comparison of GCR models with AMS hydrogen, helium, and the boron-to-carbon ratio (Norbury et al., 2018).

Advances in space radiation physics and transport at NASA.

The space radiation environment is a complex mixture of particle types and energies originating from sources inside and outside of the galaxy. These environments may be modified by the heliospheric and geomagnetic conditions as well as planetary bodies and vehicle or habitat mass shielding. In low Earth orbit (LEO), the geomagnetic field deflects a portion of the galactic cosmic rays (GCR) and all but the most intense solar particle events (SPE). There are also dynamic belts of trapped electrons and protons with low to medium energy and intense particle count rates. In deep space, the GCR exposure is more severe than in LEO and varies inversely with solar activity. Unpredictable solar storms also present an acute risk to astronauts if adequate shielding is not provided. Near planetary surfaces such as the Earth, moon or Mars, secondary particles are produced when the ambient deep space radiation environment interacts with these surfaces and/or atmospheres. These secondary particles further complicate the local radiation environment and modify the associated health risks. Characterizing the radiation fields in this vast array of scenarios and environments is a challenging task and is currently accomplished with a combination of computational models and dosimetry. The computational tools include models for the ambient space radiation environment, mass shielding geometry, and atomic and nuclear interaction parameters. These models are then coupled to a radiation transport code to describe the radiation field at the location of interest within a vehicle or habitat. Many new advances in these models have been made in the last decade, and the present review article focuses on the progress and contributions made by workers and collaborators at NASA Langley Research Center in the same time frame. Although great progress has been made, and models continue to improve, significant gaps remain and are discussed in the context of planned future missions. Of particular interest is the juxtaposition of various review committee findings regarding the accuracy and gaps of combined space radiation environment, physics, and transport models with the progress achieved over the past decade. While current models are now fully capable of characterizing radiation environments in the broad range of forecasted mission scenarios, it should be remembered that uncertainties still remain and need to be addressed.

Comparison of space radiation GCR models to recent AMS data.

This paper is the third in a series of comparisons of American (NASA) and Russian (ROSCOSMOS) space radiation calculations. The present work focuses on calculation of fluxes of galactic cosmic rays (GCR), which are a constant source of radiation that constitutes one of the major hazards during deep space exploration missions for both astronauts/cosmonauts and hardware. In this work, commonly used GCR models are compared with recently published measurements of cosmic ray Hydrogen, Helium, and the Boron-to-Carbon ratio from the Alpha Magnetic Spectrometer (AMS). All of the models were developed and calibrated prior to the publication of the AMS data; therefore this an opportunity to validate the models against an independent data set.

Solar proton exposure of an ICRU sphere within a complex structure part II: Ray-trace geometry.

A computationally efficient 3DHZETRN code with enhanced neutron and light ion (Z ≤ 2) propagation was recently developed for complex, inhomogeneous shield geometry described by combinatorial objects. Comparisons were made between 3DHZETRN results and Monte Carlo (MC) simulations at locations within the combinatorial geometry, and it was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in ray-trace geometry. This latest extension enables the code to be used within current engineering design practices utilizing fully detailed vehicle and habitat geometries. Through convergence testing, it is shown that fidelity in an actual shield geometry can be maintained in the discrete ray-trace description by systematically increasing the number of discrete rays used. It is also shown that this fidelity is carried into transport procedures and resulting exposure quantities without sacrificing computational efficiency.

Solar proton exposure of an ICRU sphere within a complex structure Part I: Combinatorial geometry.

The 3DHZETRN code, with improved neutron and light ion (Z≤2) transport procedures, was recently developed and compared to Monte Carlo (MC) simulations using simplified spherical geometries. It was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in general combinatorial geometry. A more complex shielding structure with internal parts surrounding a tissue sphere is considered and compared against MC simulations. It is shown that even in the more complex geometry, 3DHZETRN agrees well with the MC codes and maintains a high degree of computational efficiency.

Galactic cosmic ray simulation at the NASA Space Radiation Laboratory.

Most accelerator-based space radiation experiments have been performed with single ion beams at fixed energies. However, the space radiation environment consists of a wide variety of ion species with a continuous range of energies. Due to recent developments in beam switching technology implemented at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL), it is now possible to rapidly switch ion species and energies, allowing for the possibility to more realistically simulate the actual radiation environment found in space. The present paper discusses a variety of issues related to implementation of galactic cosmic ray (GCR) simulation at NSRL, especially for experiments in radiobiology. Advantages and disadvantages of different approaches to developing a GCR simulator are presented. In addition, issues common to both GCR simulation and single beam experiments are compared to issues unique to GCR simulation studies. A set of conclusions is presented as well as a discussion of the technical implementation of GCR simulation.

3DHZETRN: Neutron leakage in finite objects.

The 3DHZETRN formalism was recently developed as an extension to HZETRN with an emphasis on 3D corrections for neutrons and light ions. Comparisons to Monte Carlo (MC) simulations were used to verify the 3DHZETRN methodology in slab and spherical geometry, and it was shown that 3DHZETRN agrees with MC codes to the degree that various MC codes agree among themselves. One limitation of such comparisons is that all of the codes (3DHZETRN and three MC codes) utilize different nuclear models/databases; additionally, using a common nuclear model is impractical due to the complexity of the software. It is therefore difficult to ascertain if observed discrepancies are caused by transport code approximations or nuclear model differences. In particular, an important simplification in the 3DHZETRN formalism assumes that neutron production cross sections can be represented as the sum of forward and isotropic components, where the forward component is subsequently solved within the straight-ahead approximation. In the present report, previous transport model results in specific geometries are combined with additional results in related geometries to study neutron leakage using the Webber 1956 solar particle event as a source boundary condition. A ratio is defined to quantify the fractional neutron leakage at a point in a finite object relative to a semi-infinite slab geometry. Using the leakage ratio removes some of the dependence on the magnitude of the neutron production and clarifies the effects of angular scattering and absorption with regard to differences between the models. Discussion is given regarding observed differences between the MC codes and conclusions drawn about the need for further code development. Although the current version of 3DHZETRN is reasonably accurate compared to MC simulations, this study shows that improved leakage estimates can be obtained by replacing the isotropic/straight-ahead approximation with more detailed descriptions.

3DHZETRN: Shielded ICRU spherical phantom.

A computationally efficient 3DHZETRN code capable of simulating High (H) Charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation was recently developed for a simple homogeneous shield object. Monte Carlo benchmarks were used to verify the methodology in slab and spherical geometry, and the 3D corrections were shown to provide significant improvement over the straight-ahead approximation in some cases. In the present report, the new algorithms with well-defined convergence criteria are extended to inhomogeneous media within a shielded tissue slab and a shielded tissue sphere and tested against Monte Carlo simulation to verify the solution methods. The 3D corrections are again found to more accurately describe the neutron and light ion fluence spectra as compared to the straight-ahead approximation. These computationally efficient methods provide a basis for software capable of space shield analysis and optimization.

Low Earth orbit assessment of proton anisotropy using AP8 and AP9 trapped proton models.

The completion of the International Space Station (ISS) in 2011 has provided the space research community with an ideal evaluation and testing facility for future long duration human activities in space. Ionized and secondary neutral particles radiation measurements inside ISS form the ideal tool for validation of radiation environmental models, nuclear reaction cross sections and transport codes. Studies using thermo-luminescent detectors (TLD), tissue equivalent proportional counter (TPEC), and computer aided design (CAD) models of early ISS configurations confirmed that, as input, computational dosimetry at low Earth orbit (LEO) requires an environmental model with directional (anisotropic) capability to properly describe the exposure of trapped protons within ISS. At LEO, ISS encounters exposure from trapped electrons, protons and geomagnetically attenuated galactic cosmic rays (GCR). For short duration studies at LEO, one can ignore trapped electrons and ever present GCR exposure contributions during quiet times. However, within the trapped proton field, a challenge arises from properly estimating the amount of proton exposure acquired. There exist a number of models to define the intensity of trapped particles. Among the established trapped models are the historic AE8/AP8, dating back to the 1980s and the recently released AE9/AP9/SPM. Since at LEO electrons have minimal exposure contribution to ISS, this work ignores the AE8 and AE9 components of the models and couples a measurement derived anisotropic trapped proton formalism to omnidirectional output from the AP8 and AP9 models, allowing the assessment of the differences between the two proton models. The assessment is done at a target point within the ISS-11A configuration (circa 2003) crew quarter (CQ) of Russian Zvezda service module (SM), during its ascending and descending nodes passes through the south Atlantic anomaly (SAA). The anisotropic formalism incorporates the contributions of proton narrow pitch angle (PA) and east-west (EW) effects. Within SAA, the EW anisotropy results in different level of exposure to each side of the ISS Zvezda SM, allowing angular evaluation of the anisotropic proton spectrum. While the combined magnitude of PA and EW effects at LEO depends on a multitude of factors such as trapped proton energy, orientation and altitude of the spacecraft along the velocity vector, this paper draws quantitative conclusions on the combined anisotropic magnitude differences within ISS SM target point between AP8 and AP9 models.

Implementation of Gy-Eq for deterministic effects limitation in shield design.

The NCRP has recently defined RBE values and a new quantity (Gy-Eq) for use in estimation of deterministic effects in space shielding and operations. The NCRP's RBE for neutrons is left ambiguous and not fully defined. In the present report we will suggest a complete definition of neutron RBE consistent with the NCRP recommendations and evaluate attenuation properties of deterministic effects (Gy-Eq) in comparison with other dosimetric quantities.