Search

This paper presents the yields of secondary light ions and neutrons produced from bombardments of thick targets of aluminum by protons, helium, and iron ions. In March 2016, ion beams of 400- and 800-MeV protons, 400-AMeV helium, and 400- and 800-AMeV iron bombarded aluminum targets of 20, 40, and 60 g/cm2 thickness at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL) over the course of 100 hours. An additional aluminum target of thickness 60 g/cm2 was placed 3.5 meters downstream to model the increased secondary particle yields in an enclosed space. Surveys of neutrons and light ions were taken with organic liquid scintillators at six angles between 10° and 135° off beam axis, with further light ion measurements taken with pairs of sodium iodide (NaI) detector arrays at two positions between 10° and 30°. Utilizing the organic scintillators, neutron events were isolated with pulse shape discrimination, while charged particles were discerned by comparing energy deposited in the detector and time of flight. Energy spectra were then obtained by using time-of-flight analysis techniques. Light ion spectra were generated using NaI detector array data by comparison of energy deposition between the NaI detectors pairs. These results will be compared to simulations of the experiment calculated by MCNP6. This experiment is a part of a multi-year project to supplement the amount of measured data available for validation and verification of radiation transport computer codes used in quantifying radiation exposure and assessing cancer risk incurred during manned space flight missions.

The importance of neutrons and light ions is considered when astronauts spend considerable time in thickly shielded regions of a spacecraft. This may be relevant for space missions both in and beyond low Earth orbit (LEO). In addition to heavy ion experiments at accelerators, it is suggested that an increased emphasis on experiments with lighter ions may be useful in reducing biological uncertainties.

The computationally efficient HZETRN code has been used in recent trade studies for lunar and Martian exploration and is currently being used in the engineering development of the next generation of space vehicles, habitats, and extra vehicular activity equipment. Code development has been based on a progression of approximations first assuming all particles are produced in the initiator direction of incidence (straight-ahead) later improved by treating neutrons produced in the backward hemisphere as moving straight-back (bi-directional). A new version (3DHZETRN) capable of transporting High charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation in transverse directions is developed. Herein, new algorithms for light ion and neutron propagation with well defined convergence criteria in 3D objects is developed and tested against Monte Carlo simulations of 3D effects.

Accurate galactic cosmic ray (GCR) models are required to assess crew exposure during long-duration missions to the Moon or Mars. Many of these models have been developed and compared to available measurements, with uncertainty estimates usually stated to be less than 15%. However, when the models are evaluated over a common epoch and propagated through to effective dose, relative differences exceeding 50% are observed. This indicates that the metrics used to communicate GCR model uncertainty can be better tied to exposure quantities of interest for shielding applications. This is the first of three papers focused on addressing this need. In this work, the focus is on quantifying the extent to which each GCR ion and energy group, prior to entering any shielding material or body tissue, contributes to effective dose behind shielding. Results can be used to more accurately calibrate model-free parameters and provide a mechanism for refocusing validation efforts on measurements taken over important energy regions. Results can also be used as references to guide future nuclear cross-section measurements and radiobiology experiments. It is found that GCR with Z>2 and boundary energies below 500 MeV/n induce less than 5% of the total effective dose behind shielding. This finding is important given that most of the GCR models are developed and validated against Advanced Composition Explorer/Cosmic Ray Isotope Spectrometer (ACE/CRIS) measurements taken below 500 MeV/n. It is therefore possible for two models to very accurately reproduce the ACE/CRIS data while inducing very different effective dose values behind shielding.

This is the last of three papers focused on quantifying the uncertainty associated with galactic cosmic rays (GCR) models used for space radiation shielding applications. In the first paper, it was found that GCR ions with Z>2 and boundary energy below 500 MeV/nucleon induce less than 5% of the total effective dose behind shielding. This is an important finding since GCR model development and validation have been heavily biased toward Advanced Composition Explorer/Cosmic Ray Isotope Spectrometer measurements below 500 MeV/nucleon. Weights were also developed that quantify the relative contribution of defined GCR energy and charge groups to effective dose behind shielding. In the second paper, it was shown that these weights could be used to efficiently propagate GCR model uncertainties into effective dose behind shielding. In this work, uncertainties are quantified for a few commonly used GCR models. A validation metric is developed that accounts for measurements uncertainty, and the metric is coupled to the fast uncertainty propagation method. For this work, the Badhwar-O'Neill (BON) 2010 and 2011 and the Matthia GCR models are compared to an extensive measurement database. It is shown that BON2011 systematically overestimates heavy ion fluxes in the range 0.5-4 GeV/nucleon. The BON2010 and BON2011 also show moderate and large errors in reproducing past solar activity near the 2000 solar maximum and 2010 solar minimum. It is found that all three models induce relative errors in effective dose in the interval [-20%, 20%] at a 68% confidence level. The BON2010 and Matthia models are found to have similar overall uncertainty estimates and are preferred for space radiation shielding applications.

In order to assess the astronaut exposure received within vehicles or habitats, accurate models of the ambient galactic cosmic ray (GCR) environment are required. Many models have been developed and compared to measurements, with uncertainty estimates often stated to be within 15%. However, intercode comparisons can lead to differences in effective dose exceeding 50%. This is the second of three papers focused on resolving this discrepancy. The first paper showed that GCR heavy ions with boundary energies below 500 MeV/n induce less than 5% of the total effective dose behind shielding. Yet, due to limitations on available data, model development and validation are heavily influenced by comparisons to measurements taken below 500 MeV/n. In the current work, the focus is on developing an efficient method for propagating uncertainties in the ambient GCR environment to effective dose values behind shielding. A simple approach utilizing sensitivity results from the first paper is described and shown to be equivalent to a computationally expensive Monte Carlo uncertainty propagation. The simple approach allows a full uncertainty propagation to be performed once GCR uncertainty distributions are established. This rapid analysis capability may be integrated into broader probabilistic radiation shielding analysis and also allows error bars (representing boundary condition uncertainty) to be placed around point estimates of effective dose.

Chromosome aberrations (CAs) are one of the effects of radiation exposure and can have implications for human health in the space environment, since they are related to cancer risk. In radiation research, chromosome aberrations are a convenient biomarker for carcinogenesis. To shed light on the formation and quality of chromosome aberrations in the space environment, many experiments and simulations have been performed using chromosome aberrations in human cells, induced by heavy ions, which are present in galactic cosmic rays (GCRs). In this work, the new simulation program, radiation-induced tracks, chromosome aberrations, repair and damage (RITCARD), is presented. This software program is based on the algorithm used in the NASA Radiation Track Image (NASARTI) model with some improvements. NASARTI and RITCARD are both comprised of four parts: a random walk (RW) algorithm for simulating chromosomes in a nucleus; a deoxyribonucleic acid (DNA) damage algorithm; a break repair process; and a function to assess and count chromosome aberrations. Prior to running RITCARD, the code, relativistic ion tracks (RITRACKS), is used to simulate detailed radiation track structure and calculate time-dependent differential voxel dose maps in a parallelepiped centered on a cell nucleus. The RITCARD program reads the pre-calculated voxel dose and locates the intersections between the voxels and the chromosomes that were simulated by random walk. Radiation-induced breaks occur strictly at these intersections with a probability that is a function of the voxel dose. When a break occurs in the random walk, the corresponding chromosome piece is cut into two fragments where each has a free end at the position of the break. RITCARD generates a collection of all fragments, free ends, and enlists free end pairs. In the next step, the algorithm simulates the time-dependent rejoining of free end pairs, using different probabilities for pairs originating from a given break (proper) or from different breaks (improper), which results in the formation of fragment sequences. By grouping these sequences, the program determines the number and types of aberrations, based on the same criteria used in our experiment. The new program is used to assess the yields of various types of chromosome aberrations in human fibroblast cells for several ions (1H+, 4He2+, 12C6+, 16O8+, 20Ne10+, 28Si14+, 48Ti22+ and 56Fe26+) with energies varying from 10 to 1,000 MeV/n. The results show linear and linear-quadratic dose dependence for most chromosome aberrations types. The calculation results were compared with those obtained by fluorescence in situ hybridization (FISH) experiments that were performed by our group. The simulations and experiments are in better agreement at lower LET. Regarding the simulation results, the coefficient of the linear part of the dose-dependence curve also peaks at an LET value of approximately 100 keV/lm, which evokes a relative biological effectiveness (RBE) peak found by other researchers.

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.

Safe and efficient mission operations in space require an accurate understanding of the physical interactions of space radiation. As the primary space radiation interacts with intervening materials, the composition and spectrum of the radiation environment changes. The production of secondary particles can make a significant contribution to radiation exposure. In this work, the NASA space radiation transport code, HZETRN, is extended to include the transport of electrons, positrons, and photons. The production of these particles is coupled to the initial cosmic ray radiation environment through the decay of neutral pions, which produce high energy photons, and through the decay of muons, which produce electrons and positrons. The photons, electrons, and positrons interact with materials producing more photons, electrons and positrons generating an electromagnetic cascade. The relevant cross sections, transport equation, and solution method are introduced. Electron and positron production in Earth’s atmosphere is investigated and compared to experimental balloon-flight measurements. Reasonable agreement is seen between HZETRN and data.

There is growing interest in flying biological experiments beyond low-Earth orbit (LEO) to measure biological responses potentially relevant to those expected during a human mission to Mars. Such experiments could be payloads onboard precursor missions, including unmanned private-public partnerships, as well as small low-cost spacecraft (satellites) designed specifically for biosentinel-type missions. It is the purpose of this paper to provide physical cosmic-ray interaction data and related information useful to biologists who may be planning such experiments. It is not the objective here to actually design such experiments or provide radiobiological response functions, which would be specific for each experiment and biological endpoint. Nuclide-specific flux and dose rates were calculated using OLTARIS and these results were used to determine particle traversal rates and doses in hypothetical biological targets. Comparisons are provided between GCR in interplanetary space and inside the ISS. Calculated probabilistic estimates of dose from solar particle events are also presented. Although the focus here is on biological experiments, the information provided may be useful for designing other payloads as well if the space radiation environment is a factor to be considered.

In March 2016 accelerator-based experiments colliding protons (0.4 and 0.8 GeV), helium (0.4 AGeV) and iron (0.4 and 0.8 AGeV) on thick aluminum targets with surface densities of 20, 40, and 60 g/cm2 were performed at the National Aeronautics and Space Administration Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory. Two targets were utilized in each experimental configuration. Hydrogen and helium ions were detected using organic liquid scintillators in conjunction with thin plastic scintillators at 10°, 30°, 45°, 60°, 80°, and 135° from beam axis. Time-of-flight techniques and pulse shape discrimination were used to identify light ion species in order to generate double differential energy spectra of the light ion fragments. Comparisons of these measured yields were compared with Monte Carlo calculations generated by MCNP6. These yields will be used to quantify uncertainty in radiation transport codes utilized in risk assessment for spaceflight missions with prolonged exposures to galactic cosmic rays.

The radiation environment at the Martian surface is, apart from occasional solar energetic particle events, dominated by galactic cosmic radiation, secondary particles produced in their interaction with the Martian atmosphere and albedo particles from the Martian regolith. The highly energetic primary cosmic radiation consists mainly of fully ionized nuclei creating a complex radiation field at the Martian surface. This complex field, its formation and its potential health risk posed to astronauts on future manned missions to Mars can only be fully understood using a combination of measurements and model calculations. In this work the outcome of a workshop held in June 2016 in Boulder, CO, USA is presented: experimental results from the Radiation Assessment Detector of the Mars Science Laboratory are compared to model results from GEANT4, HETC-HEDS, HZETRN, MCNP6, and PHITS. Charged and neutral particle spectra and dose rates measured between 15 November 2015 and 15 January 2016 and model results calculated for this time period are investigated.

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) flux model has been used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a comprehensive database of GCR measurements from various particle detectors to determine boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

Experimental work on light charged ion production from thick target shielding began this past May at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL). This paper presents the measured secondary light charged ion and neutron yields produced by 0.4- and 2.5-GeV protons and 0.4- and 1.0-AGeV iron ions striking a 30 g/cm2 aluminum target. Neutron and light charged ion (protons, deuterons, and tritons) measurements were taken with liquid scintillators and sodium iodide (NaI) detectors positioned at seven locations between 10 and 135 degrees off the beam axis to best cover the angular distributions of secondary particles, as determined by MCNPX simulations. In the liquid scintillators, neutron-gamma separation was achieved with pulse shape discrimination, and particle species were identified and isolated by analyzing the total charge deposited in the detector versus particle time of flight (TOF). After isolation, the TOF technique was utilized to produce energy spectra for protons, deuterons, tritons, and neutrons at various locations. Additionally, the stopping powers of light charged ions were compared in NaI detector pairs to generate energy spectra for protons, deuterons, and tritons. Preliminary results demonstrated promising agreement with MCNPX Monte Carlo transport code simulations for protons, deuterons, tritons, and neutrons, despite the lack of a full background characterization and optimization of detector settings. Results are expected to improve over the next three years with an increase in beam time, inclusion of specific liquid scintillator detection efficiencies, and an investigation of the physics model parameters in MCNPX. Future experiments will include the use of both forward and back targets composed of high-density polyethylene or aluminum with thicknesses of 20, 40, and 60 g/cm2. Furthermore, proton, helium, carbon, silicon, and iron projectiles will be utilized at energies of 0.4, 0.75, 1.5, and 2.5 AGeV. Ultimately, these measurements will be incorporated in the uncertainty analysis for the engineering codes that NASA uses to develop optimal shielding thicknesses for spacecraft and space habitat design.

Light ion improvements to the nuclear fragmentation model, NUCFRG, are reported. Improvements include the replacement of the simple light ion production model with a light ion coalescence model and an improved electromagnetic dissociation (EMD) formalism. Prior versions of the model provide reasonable overall agreement with measured data; however, those versions lack a physics-based description for coalescence and EMD. The version reported herein, NUCFRG3, has improved the theoretical descriptions of these mechanisms and offers additional benefits, such as the capability to calculate EMD cross-sections for single deuteron, triton, helion, and alpha particle emission. NUCFRG3 model evaluation and validation show that the predictive capability has been improved and strengthened by the light ion physics-based changes. Based on increased capability and better theoretical grounding, it is recommended that NUCFRG3 replace its predecessors for space radiation assessments and other applications.

During the design process, the configuration of space vehicles and habitats changes frequently and the merits of design changes must be evaluated. Methods for rapidly assessing astronaut exposure are therefore required. Typically, approximations are made to simplify the geometry and speed-up the evaluation of each design. In this work, the error associated with two common approximations used to simplify space radiation vehicle analyses, scaling into equivalent materials and material reordering, are investigated. Over thirty materials commonly found in spacesuits, vehicles, and human bodies are considered. Each material is placed in a material group (aluminum, polyethylene, or tissue), and the error associated with scaling and reordering was quantified for each material. Of the scaling methods investigated, range scaling is shown to be the superior method, especially for shields less than 30 g/cm2 exposed to a solar particle event. More complicated, realistic slabs are examined to quantify the separate and combined effects of using equivalent materials and reordering. The error associated with material reordering is shown to be at least comparable to, if not greater than, the error associated with range scaling. In general, scaling and reordering errors were found to grow with the difference between the average nuclear charge of the actual material and average nuclear charge of the equivalent material. Based on this result, a different set of equivalent materials (titanium, aluminum, and tissue) are substituted for the commonly used aluminum, polyethylene, and tissue. The realistic cases are scaled and reordered using the new equivalent materials, and the reduced error is shown.

The On-Line Tool for the Assessment of Radiation In Space (OLTARIS) is a World Wide Web based tool that assesses the effects of space radiation to humans in items such as spacecraft, habitats, rovers, and spacesuits. This document explains the basis behind the interface and framework used to input the data, perform the assessment, and output the results to the user as well as the physics, engineering, and computer science used to develop OLTARIS. The physics is based on the HZETRN2005 and NUCFRG2 research codes. The OLTARIS website is the successor to the SIREST website from the early 2000 s. Modifications have been made to the code to enable easy maintenance, additions, and configuration management along with a more modern web interface. Over all, the code has been verified, tested, and modified to enable faster and more accurate assessments. The next major areas of modification are more accurate transport algorithms, better uncertainty estimates, and electronic response functions. Improvements in the existing algorithms and data occur continuously and are logged in the change log section of the website.

The protection of astronauts and instrumentation from galactic cosmic rays and solar particle events is one of the primary constraints associated with mission planning in low earth orbit or deep space. To help satisfy this constraint, several computational tools have been developed to analyze the effectiveness of various shielding materials and structures exposed to space radiation. These tools are now being carefully scrutinized through a systematic effort of verification, validation, and uncertainty quantification. In this benchmark study, the deterministic transport code HZETRN is compared to the Monte Carlo transport codes HETC-HEDS and FLUKA for a 30 g/cm 2 water target protected by a 20 g/cm 2 aluminum shield exposed to a parameterization of the February 1956 solar particle event. Neutron and proton fluences as well as dose and dose equivalent are compared at various depths in the water target. The regions of agreement and disagreement between the three codes are quantified and discussed, and recommendations for future work are given.

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 2 aluminum slab followed by a 30 g/cm 2 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 2 H, 3 H, 3 He and 4 He fluxes are also examined.

Long term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures, and vehicles. In this paper, the numerical error associated with energy discretization in HZETRN is addressed. An inadequate numerical integration scheme in the transport algorithm is shown to produce large errors in the low energy portion of the neutron and light ion fluence spectra. It is further shown that the errors result from the narrow energy domain of the neutron elastic cross section spectral distributions, and that an extremely fine energy grid is required to resolve the problem under the current formulation. Two numerical methods are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. Convergence testing is completed by running the code for various environments and shielding materials with various energy grids to ensure stability of the newly implemented method.

To estimate astronaut health risk due to space radiation, one must have the ability to calculate various exposure-related quantities that are averaged over specific organs and tissue types. Such calculations require computational models of the ambient space radiation environment, particle transport, nuclear and atomic physics, and the human body. While significant efforts have been made to verify, validate, and quantify the uncertainties associated with many of these models and tools, relatively little work has focused on the uncertainties associated with the representation and utilization of the human phantoms. In this study, we first examine the anatomical properties of the Computerized Anatomical Man (CAM), Computerized Anatomical Female (CAF), Male Adult voXel (MAX), and Female Adult voXel (FAX) models by comparing the masses of various model tissues used to calculate effective dose to the reference values specified by the International Commission on Radiological Protection (ICRP). The MAX and FAX tissue masses are found to be in good agreement with the reference data, while major discrepancies are found between the CAM and CAF tissue masses and the reference data for almost all of the effective dose tissues. We next examine the distribution of target points used with the deterministic transport code HZETRN (High charge ( Z ) and Energy TRaNsport) to compute mass averaged exposure quantities. A numerical algorithm is presented and used to generate multiple point distributions of varying fidelity for many of the effective dose tissues identified in CAM, CAF, MAX, and FAX. The point distributions are used to compute mass averaged dose equivalent values under both a galactic cosmic ray (GCR) and solar particle event (SPE) environment impinging isotropically on three spherical aluminum shells with areal densities of 0.4 g/cm 2 , 2.0 g/cm 2 , and 10.0 g/cm 2 . The dose equivalent values are examined to identify a recommended set of target points for each of the tissues and to further assess the differences between CAM, CAF, MAX, and FAX. It is concluded that the previously published CAM and CAF point distributions were significantly under-sampled and that the set of point distributions presented here should be adequate for future studies involving CAM, CAF, MAX, or FAX. It is also found that the errors associated with the mass and location of certain tissues in CAM and CAF have a significant impact on the mass averaged dose equivalent values, and it is concluded that MAX and FAX are more accurate than CAM and CAF for space radiation analyses.

For the evaluation of organ dose and dose equivalent of astronauts on space shuttle and the International Space Station (ISS) missions, the CAMERA models of CAM (Computerized Anatomical Male) and CAF (Computerized Anatomical Female) of human tissue shielding have been implemented and used in radiation transport model calculations at NASA. One of new human geometry models to meet the “reference person” of International Commission on Radiological Protection (ICRP) is based on detailed Voxel (volumetric and pixel) phantom models denoted for male and female as MAX (Male Adult voXel) and FAX (Female Adult voXel), respectively. We compared the CAM model predictions of organ doses to those of MAX model, since the MAX model represents the male adult body with much higher fidelity than the CAM model currently used at NASA. Directional body-shielding mass was evaluated for over 1500 target points of MAX for specified organs considered to be sensitive to the induction of stochastic effects. Radiation exposures to solar particle event (SPE), trapped protons, and galactic cosmic ray (GCR) were assessed at the specific sites in the MAX phantom by coupling space radiation transport models with the relevant body-shielding mass. The development of multiple-point body-shielding distributions at each organ made it possible to estimate the mean and variance of organ doses at the specific organ. For the estimate of doses to the blood forming organs (BFOs), data on active marrow distributions in adult were used to weight the bone marrow sites over the human body. The discrete number of target points of MAX organs resulted in a reduced organ dose and dose equivalent compared to the results of CAM organs especially for SPE, and should be further investigated. Differences of effective doses between the two approaches were found to be small (

Exposure estimates inside space vehicles, surface habitats, and high altitude aircrafts exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETC-HEDS, FLUKA, and MCNPX, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light particle transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.

Exposure estimates inside space vehicles, surface habitats, and high altitude aircraft exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETCHEDS and FLUKA, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light ion (A

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.

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.

Context: The Radiation Assessment Detector (RAD) on the Mars Science Laboratory (MSL) has been measuring the radiation environment on the surface of Mars since August 6th 2012. MSL-RAD is the first instrument to provide detailed information about charged and neutral particle spectra and dose rates on the Martian surface, and one of the primary objectives of the RAD investigation is to help improve and validate current radiation transport models. Aims: Applying different numerical transport models with boundary conditions derived from the MSL-RAD environment the goal of this work was to both provide predictions for the particle spectra and the radiation exposure on the Martian surface complementing the RAD sensitive range and, at the same time, validate the results with the experimental data, where applicable. Such validated models can be used to predict dose rates for future manned missions as well as for performing shield optimization studies. Methods: Several particle transport models (GEANT4, PHITS, HZETRN/OLTARIS) were used to predict the particle flux and the corresponding radiation environment caused by galactic cosmic radiation on Mars. From the calculated particle spectra the dose rates on the surface are estimated. Results: Calculations of particle spectra and dose rates induced by galactic cosmic radiation on the Martian surface are presented. Although good agreement is found in many cases for the different transport codes, GEANT4, PHITS, and HZETRN/OLTARIS, some models still show large, sometimes order of magnitude discrepancies in certain particle spectra. We have found that RAD data is helping to make better choices of input parameters and physical models. Elements of these validated models can be applied to more detailed studies on how the radiation environment is influenced by solar modulation, Martian atmosphere and soil, and changes due to the Martian seasonal pressure cycle. By extending the range of the calculated particle spectra with respect to the experimental data additional information about the radiation environment is gained, and the contribution of different particle species to the dose is estimated.

NASA currently uses one-dimensional deterministic transport to generate values of the organ dose equivalent needed to calculate stochastic radiation risk following crew space exposures. In this study, organ absorbed doses and dose equivalents are calculated for 50th percentile male and female astronaut phantoms using both the NASA High Charge and Energy Transport Code to perform one-dimensional deterministic transport and the Particle and Heavy Ion Transport Code System to perform three-dimensional Monte Carlo transport. Two measures of radiation risk, effective dose and risk of exposure-induced death (REID) are calculated using the organ dose equivalents resulting from the two methods of radiation transport. For the space radiation environments and simplified shielding configurations considered, small differences (

In March 2016, secondary neutron production from thick-target shielding experiments were conducted at the National Aeronautics and Space Administration’s (NASA) Space Radiation Laboratory at Brookhaven National Laboratory. Ion beams of proton, helium, and iron projectiles were aimed at aluminum targets with areal densities of 20, 40, and 60 g/cm2 . The ion beams were extracted at energies of 400 and 800 AMeV and neutron yields were measured with liquid scintillators at 10°, 30°, 45°, 60°, 80°, and 135° off the beam axis. A second 60 g/cm2 aluminum target was placed 3.5 m downstream from the middle of front target to study backscattered neutrons. Double differential thick-target neutron yields for various combinations of projectile, projectile energy, target material, target thickness, and detector location were produced using the time-of-flight technique. These measurements will help NASA perform uncertainty analyses on their transport codes and contribute to shielding design studies for future space applications.

The near-Earth space radiation environment includes energetic galactic cosmic rays (GCR), high intensity proton and electron belts, and the potential for solar particle events (SPE). These sources may penetrate shielding materials and deposit significant energy in sensitive electronic devices on board spacecraft and satellites. Material and design optimization methods may be used to reduce the exposure and extend the operational lifetime of individual components and systems. Since laboratory experiments are expensive and may not cover the range of particles and energies relevant for space applications, such optimization may be done computationally with efficient algorithms that include the various constraints placed on the component, system, or mission. In the present work, the web-based tool OLTARIS (On-Line Tool for the Assessment of Radiation in Space) is presented, and the applicability of the tool for rapidly analyzing exposure levels within either complicated shielding geometries or user-defined material slabs exposed to space radiation is demonstrated. An example approach for material optimization is also presented. Slabs of various advanced multifunctional materials are defined and exposed to several space radiation environments. The materials and thicknesses defining each layer in the slab are then systematically adjusted to arrive at an optimal slab configuration.

The radiation environment on the Moon includes albedo neutrons produced by primary particles interacting with the lunar surface. In this work, HZETRN2010 is used to calculate the albedo neutron contribution to effective dose as a function of shielding thickness for four different space radiation environments and to determine to what extent various factors affect such estimates. First, albedo neutron spectra computed with HZETRN2010 are compared to Monte Carlo results in various radiation environments. Next, the impact of lunar regolith composition on the albedo neutron spectrum is examined, and the variation on effective dose caused by neutron fluence-to-effective dose conversion coefficients is studied. A methodology for computing effective dose in detailed human phantoms using HZETRN2010 is also discussed and compared. Finally, the combined variation caused by environmental models, shielding materials, shielding thickness, regolith composition and conversion coefficients on the albedo neutron contribution to effective dose is determined. It is shown that a single percentage number for characterizing the albedo neutron contribution to effective dose can be misleading. In general, the albedo neutron contribution to effective dose is found to vary between 1-32%, with the environmental model, shielding material and shielding thickness being the driving factors that determine the exact contribution. It is also shown that polyethylene or other hydrogen-rich materials may be used to mitigate the albedo neutron exposure.

The On-Line Tool for the Assessment of Radiation in Space (OLTARIS, https://oltaris.nasa.gov) is a web-based set of tools and models that allows engineers and scientists to assess the effects of space radiation on spacecraft, habitats, rovers, and spacesuits. The site is intended to be a design tool for those studying the effects of space radiation for current and future missions as well as a research tool for those developing advanced material and shielding concepts. The tools and models are built around the HZETRN radiation transport code and are primarily focused on human- and electronic-related responses. The focus of this paper is to highlight new capabilities that have been added to support deep space (outside Low Earth Orbit) missions. Specifically, the electron, proton, and heavy ion design environments for the Europa mission have been incorporated along with an efficient coupled electron-photon transport capability to enable the analysis of complicated geometries and slabs exposed to these environments. In addition, a neutron albedo lunar surface environment was also added, that will be of value for the analysis of surface habitats. These updates will be discussed in terms of their implementation and on how OLTARIS can be used by instrument vendors, mission designers, and researchers to analyze their specific requirements.12

Risk to astronauts due to ionizing radiation exposure is a primary concern for missions beyond Low Earth Orbit (LEO) and will drive mission architecture requirements, mission timelines, and operational practices. For short missions, radiation risk is dominated by the possibility of a large Solar Particle Event (SPE). Longer duration missions have both SPE and Galactic Cosmic Ray (GCR) risks. SPE exposure can contribute significantly toward cancer induction in combination with GCR. As mission duration increases, mitigation strategies must address the combined risks from SPE and GCR exposure. In this paper, full mission exposure assessments were performed for the proposed long duration lunar surface mission scenarios. In order to accomplish these assessments, previously developed radiation shielding models for a proposed lunar habitat and rover were utilized. End-to-End mission exposure assessments were performed by first calculating exposure rates for locations in the habitat, rover, and during Extra-Vehicular Activities (EVA). Subsequently, total mission exposures were evaluated for the proposed timelines. Mission exposure results, assessed in terms of effective dose, are presented for the proposed timelines and recommendations are made for improved astronaut shielding and safer operational practices.