PHITS simulations of the Matroshka experiment

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

Simulations of absorbed dose on the phantom surface of MATROSHKA-R experiment at the ISS

The health risks associated with exposure to various components of space radiation are of great concern when planning manned long-term interplanetary missions, such as future missions to Mars. Since it is not possible to measure the radiation environment inside of human organs in deep space, simulations based on radiation transport/interaction codes coupled to phantoms of tissue equivalent materials are used. However, the calculated results depend on the models used in the codes, and it is therefore necessary to verify their validity by comparison with measured data. The goal of this paper is to compare absorbed doses obtained in the MATROSHKA-R experiment performed at the International Space Station (ISS) with simulations performed with the three-dimensional Monte Carlo Particle and Heavy-Ion Transport code System (PHITS). The absorbed dose was measured using passive detectors (packages of thermoluminescent and plastic nuclear track detectors) placed on the surface of the spherical tissue equivalent phantom MATROSHKA-R, which was exposed aboard the ISS in the Service Zvezda Module from December 2005 to September 2006. The data calculated by PHITS assuming an ISS shielding of 3 g/cm² and 5 g/cm² aluminum mass thickness were in good agreement with the measurements. Using a simplified geometrical model of the ISS, the influence of variations in altitude and wall mass thickness of the ISS on the calculated absorbed dose was estimated. The uncertainties of the calculated data are also discussed; the relative expanded uncertainty of absorbed dose in phantom was estimated to be 44% at a 95% confidence level.     

PHITS simulations of the Protective curtain experiment onboard the Service module of ISS: Comparison with absorbed doses measured with TLDs

“Protective curtain” was the physical experiment onboard the International Space Station (ISS) aimed on radiation measurement of the dose – reducing effect of the additional shielding made of hygienic water-soaked wipes and towels placed on the wall in the crew cabin of the Service module Zvezda. The measurements were performed with 12 detector packages composed of thermoluminescent detectors (TLDs) and plastic nuclear track detectors (PNTDs) placed at the Protective curtain, so that they created pairs of shielded and unshielded detectors.     

An update about recent developments of the PHITS code

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

IVIDIL experiment onboard the ISS

The experiment IVIDIL (Influence of Vibrations on Diffusion in Liquids) is scheduled to be performed in forthcoming fall 2009 onboard the ISS, inside the SODI instrument mounted in the Glovebox on the ESA Columbus module. It is planned to carry out 39 experimental runs with each of them lasting 18 h. The objective of the experiment is threefold.     

Space Shuttle drops down the SAA doses on ISS

Long-term analysis of data from two radiation detection instruments on the International Space Station (ISS) shows that the docking of the Space Shuttle drops down the measured dose rates in the region of the South Atlantic Anomaly (SAA) by a factor of 1.5–3. Measurements either by the R3DE detector, which is outside the ISS at the EuTEF facility on the Columbus module behind a shielding of less than 0.45 g cm⁻², and by the three detectors of the Liulin-5 particle telescope, which is inside the Russian PEARS module in the spherical tissue equivalent phantom behind much heavier shielding demonstrate that effect. Simultaneously the estimated averaged incident energies of the incoming protons rise up from about 30 to 45 MeV. The effect is explained by the additional shielding against the SAA 30–150 MeV protons, provided by the 78 tons Shuttle to the instruments inside and outside of the ISS. An additional reason is the ISS attitude change (performed for the Shuttle docking) leading to decreasing of dose rates in two of Liulin-5 detectors because of the East–West proton fluxes asymmetry in SAA. The Galactic Cosmic Rays dose rates are practically not affected.     

Statistical validation of HZETRN as a function of vertical cutoff rigidity using ISS measurements

Measurements taken in Low Earth Orbit (LEO) onboard the International Space Station (ISS) and transit vehicles have been extensively used to validate radiation transport models. Primarily, such comparisons were done by integrating measured data over mission or trajectory segments so that individual comparisons to model results could be made. This approach has yielded considerable information but is limited in its ability to rigorously quantify and differentiate specific model errors or uncertainties. Further, as exploration moves beyond LEO and measured data become sparse, the uncertainty estimates derived from these validation cases will no longer be applicable. Recent improvements in the underlying numerical methods used in HZETRN have resulted in significant decreases in code run time. Therefore, the large number of comparisons required to express error as a function of a physical quantity, like cutoff rigidity, are now possible. Validation can be looked at in detail over any portion of a flight trajectory (e.g. minute by minute) such that a statistically significant number of comparisons can be made. This more rigorous approach to code validation will allow the errors caused by uncertainties in the geometry models, environmental models, and nuclear physics models to be differentiated and quantified. It will also give much better guidance for future model development. More importantly, it will allow a quantitative means of extrapolating uncertainties in LEO to free space. In this work, measured data taken onboard the ISS during solar maximum are compared to results obtained with the particle transport code HZETRN. Comparisons are made at a large number (∼77,000) of discrete time intervals, allowing error estimates to be given as a function of cutoff rigidity. It is shown that HZETRN systematically underestimates exposure quantities at high cutoff rigidity. The errors are likely associated with increased angular variation in the geomagnetic field near the equator, the lack of pion production in HZETRN, and errors in high energy nuclear physics models, and will be the focus of future work.     

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.     

Comparison of the mean quality factors for astronauts calculated using the Q-functions proposed by ICRP,ICRU, and NASA

For the estimation of the radiation risk for astronauts, not only the organ absorbed doses but also their mean quality factors must be evaluated. Three functions have been proposed by different organizations for expressing the radiation quality, including the Q(L), Q(y), and Q_NASA(Z, E) relationships as defined in International Committee of Radiological Protection (ICRP) Publication 60, International Commission on Radiation Units and Measurements (ICRU) Report 40, and National Aeronautics and Space Administration (NASA) TP-2011-216155, respectively. The Q(L) relationship is the most simple and widely used for space dosimetry, but the use of the latter two functions enables consideration of the difference in the track structure of various charged particles during the risk estimation.     

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.     

An improved neutron transport algorithm for HZETRN

Long-term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures and vehicles. It is necessary to expose the numerical tools commonly used in radiation analyses to extensive verification, validation and uncertainty quantification. 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. Since adding a sufficient number of energy points will render the code computationally inefficient, we revisit the light ion and neutron transport theory developed for HZETRN and focus on neutron elastic interactions. Two numerical methods (average value and collocation) are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. An energy grid convergence study is conducted to demonstrate the improved stability of the new methods. Based on the results of the convergence study and the ease of implementation, the average value method with a 100 point energy grid is found to be suitable for future use in HZETRN.     

Preparation of the Biochip experiment on the EXPOSE-R2 mission outside the International Space Station

Biochips might be suited for planetary exploration. Indeed, they present great potential for the search for biomarkers – molecules that are the sign of past or present life in space – thanks to their size (miniaturized devices) and sensitivity. Their detection principle is based on the recognition of a target molecule by affinity receptors fixed on a solid surface. Consequently, one of the main concerns when developing such a system is the behavior of the biological receptors in a space environment. In this paper, we describe the preparation of an experiment planned to be part of the EXPOSE-R2 mission, which will be conducted on the EXPOSE-R facility, outside the International Space Station (ISS), in order to study the resistance of biochip models to space constraints (especially cosmic radiation and thermal cycling). This experiment overcomes the limits of ground tests which do not reproduce exactly the space parameters. Indeed, contrary to ground experiments where constraints are applied individually and in a limited time, the biochip models on the ISS will be exposed to cumulated constraints during several months. Finally, this ISS experiment is a necessary step towards planetary exploration as it will help assessing whether a biochip can be used for future exploration missions.     

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.     

Characterization of the near Earth radiation environment by Liulin type spectrometers

Comprehensive study of the dose, flux and deposited energy spectra shape data obtained by Liulin type spectrometers on spacecraft (five different experiments) and aircraft since 2001 is performed with the aim of understanding how well these parameters can characterize the type of predominant particles and their energy in the near Earth radiation environment. Three different methods for characterisation of the incoming radiation from Liulin spectrometers are described. The results revealed that the most informative one is by the shape of the deposited energy spectra. Spectra generated by Galactic Cosmic Rays (GCR) protons and their secondaries are with linear falling shape in the coordinates deposited energy/deposited per channel dose rate. The position of the maximum of the deposited energy spectra inside the South Atlantic Anomaly (SAA) region depends on the incident energy of the incoming protons. Spectra generated by relativistic electrons in the outer radiation belt have a maximum in the first channels. For higher energy depositions these spectra are similar to the GCR spectra. Mixed radiation by protons and electrons and/or bremsstrahlung is characterized by spectra with 2 maxima. All type of spectra has a knee close to 6.2 MeV deposited energy, which correspond to the stopping energy of protons in the detector. Dose to flux ratio known also as specific dose is another high information parameter, which is given by experimentally obtained formulae [Heffner, J. Nuclear radiation and safety in space. M. Atomizdat. 115, 1971 (in Russian)] connecting the dose to flux ratio and the incident energy of the particles.     

卫星舱内两类剂量计测量结果的比较

本文利用1994年和1996年两次返回式卫星的搭载条件对舱内辐射剂量进行了对比测量．通过对比测量,研究不同掺杂、不同厚度LiF剂量计测量空间辐射剂量的特点;研究GM计数管计数和LiF剂量间的转换系数以及转换系数随屏蔽状况的变化;由剂量和GM计数研究粒子平均碰撞阻止本领的估计方法．结果表明,不同掺杂、不同厚度的LiF剂量计测量结果间无显著差异,而转换系数几乎不受舱内位置和屏蔽状态的影响．不同厚度LiF剂量计,不同屏蔽状态的GM计数管计数和剂量—计数转换系数的比较研究以及对粒子平均碰撞阻止本领的估计表明,舱内辐射剂量起决定作用的是高能粒子成分,其平均碰撞阻止本领估计约为5MeV／cm．     

Depth dose measurements with the Liulin-5 experiment inside the spherical phantom of the MATROSHKA-R project onboard the International Space Station

The Liulin-5 experiment is a part of the international project MATROSHKA-R on the Russian segment of the ISS, which uses a tissue-equivalent spherical phantom equipped with a set of radiation detectors. The objective of the MATROSHKA-R project is to provide depth dose distribution of the radiation field inside the sphere in order to get more information on the distribution of dose in a human body. Liulin-5 is a charged particle telescope using three silicon detectors. It measures time resolved energy deposition spectra, linear energy transfer (LET) spectra, particle flux, and absorbed doses of electrons, protons and heavy ions, simultaneously at three depths along the radius of the phantom. Measurements during the minimum of the solar activity in cycle 23 show that the average absorbed daily doses at 40 mm depth in the phantom are between 180 μGy/day and 220 μGy/day. The absorbed doses at 165 mm depth in the phantom decrease by a factor of 1.6–1.8 compared to the doses at 40 mm depth due to the self-shielding of the phantom from trapped protons. The average dose equivalent at 40 mm depth is 590 ± 32 μSV/day and the galactic cosmic rays (GCR) contribute at least 70% of the total dose equivalent at that depth. Shown is that due to the South Atlantic Anomaly (SAA) trapped protons asymmetry and the direction of Liulin-5 lowest shielding zone the dose rates on ascending and descending nodes in SAA are different. The data obtained are compared to data from other radiation detectors on ISS.     

Utilization of CAM,CAF, MAX,and FAX for space radiation analyses using HZETRN

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.0 g/cm², and 10.0 g/cm². 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.