Amir A. Bahadori

Response functions for computing absorbed dose to skeletal tissues from photon irradiation--an update.

A comprehensive set of photon fluence-to-dose response functions (DRFs) is presented for two radiosensitive skeletal tissues-active and total shallow marrow-within 15 and 32 bone sites, respectively, of the ICRP reference adult male. The functions were developed using fractional skeletal masses and associated electron-absorbed fractions as reported for the UF hybrid adult male phantom, which in turn is based upon micro-CT images of trabecular spongiosa taken from a 40 year male cadaver. The new DRFs expand upon both the original set of seven functions produced in 1985, and a 2007 update calculated under the assumption of secondary electron escape from spongiosa. In this study, it is assumed that photon irradiation of the skeleton will yield charged particle equilibrium across all spongiosa regions at energies exceeding 200 keV. Kerma coefficients for active marrow, inactive marrow, trabecular bone and spongiosa at higher energies are calculated using the DRF algorithm setting the electron-absorbed fraction for self-irradiation to unity. By comparing kerma coefficients and DRF functions, dose enhancement factors and mass energy-absorption coefficient (MEAC) ratios for active marrow to spongiosa were derived. These MEAC ratios compared well with those provided by the NIST Physical Reference Data Library (mean difference of 0.8%), and the dose enhancement factors for active marrow compared favorably with values calculated in the well-known study published by King and Spiers (1985 Br. J. Radiol. 58 345-56) (mean absolute difference of 1.9 percentage points). Additionally, dose enhancement factors for active marrow were shown to correlate well with the shallow marrow volume fraction (R(2) = 0.91). Dose enhancement factors for the total shallow marrow were also calculated for 32 bone sites representing the first such derivation for this target tissue.

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.

A semiconductor radiation imaging pixel detector for space radiation dosimetry

Progress in the development of high-performance semiconductor radiation imaging pixel detectors based on technologies developed for use in high-energy physics applications has enabled the development of a completely new generation of compact low-power active dosimeters and area monitors for use in space radiation environments. Such detectors can provide real-time information concerning radiation exposure, along with detailed analysis of the individual particles incident on the active medium. Recent results from the deployment of detectors based on the Timepix from the CERN-based Medipix2 Collaboration on the International Space Station (ISS) are reviewed, along with a glimpse of developments to come. Preliminary results from Orion MPCV Exploration Flight Test 1 are also presented.

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.

The effect of anatomical modeling on space radiation dose estimates: a comparison of doses for NASA phantoms and the 5th, 50th, and 95th percentile male and female astronauts

The National Aeronautics and Space Administration (NASA) performs organ dosimetry and risk assessment for astronauts using model-normalized measurements of the radiation fields encountered in space. To determine the radiation fields in an organ or tissue of interest, particle transport calculations are performed using self-shielding distributions generated with the computer program CAMERA to represent the human body. CAMERA mathematically traces linear rays (or path lengths) through the computerized anatomical man (CAM) phantom, a computational stylized model developed in the early 1970s with organ and body profiles modeled using solid shapes and scaled to represent the body morphometry of the 1950 50th percentile (PCTL) Air Force male. With the increasing use of voxel phantoms in medical and health physics, a conversion from a mathematical-based to a voxel-based ray-tracing algorithm is warranted. In this study, the voxel-based ray tracer (VoBRaT) is introduced to ray trace voxel phantoms using a modified version of the algorithm first proposed by Siddon (1985 Med. Phys. 12 252-5). After validation, VoBRAT is used to evaluate variations in body self-shielding distributions for NASA phantoms and six University of Florida (UF) hybrid phantoms, scaled to represent the 5th, 50th, and 95th PCTL male and female astronaut body morphometries, which have changed considerably since the inception of CAM. These body self-shielding distributions are used to generate organ dose equivalents and effective doses for five commonly evaluated space radiation environments. It is found that dosimetric differences among the phantoms are greatest for soft radiation spectra and light vehicular shielding.

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.

Calculation of Organ Doses for a Large Number of Patients Undergoing CT Examinations

OBJECTIVE The objective of our study was to develop an automated calculation method to provide organ dose assessment for a large cohort of pediatric and adult patients undergoing CT examinations. MATERIALS AND METHODS We adopted two dose libraries that were previously published: the volume CT dose index-normalized organ dose library and the tube current-exposure time product (100 mAs)-normalized weighted CT dose index library. We developed an algorithm to calculate organ doses using the two dose libraries and the CT parameters available from DICOM data. We calculated organ doses for pediatric (n = 2499) and adult (n = 2043) CT examinations randomly selected from four health care systems in the United States and compared the adult organ doses with the values calculated from the ImPACT calculator. RESULTS The median brain dose was 20 mGy (pediatric) and 24 mGy (adult), and the brain dose was greater than 40 mGy for 11% (pediatric) and 18% (adult) of the head CT studies. Both the National Cancer Institute (NCI) and ImPACT methods provided similar organ doses (median discrepancy < 20%) for all organs except the organs located close to the scanning boundaries. The visual comparisons of scanning coverage and phantom anatomies revealed that the NCI method, which is based on realistic computational phantoms, provides more accurate organ doses than the ImPACT method. CONCLUSION The automated organ dose calculation method developed in this study reduces the time needed to calculate doses for a large number of patients. We have successfully used this method for a variety of CT-related studies including retrospective epidemiologic studies and CT dose trend analysis studies.

Response functions for computing absorbed dose to skeletal tissues from neutron irradiation

Spongiosa in the adult human skeleton consists of three tissues-active marrow (AM), inactive marrow (IM) and trabecularized mineral bone (TB). AM is considered to be the target tissue for assessment of both long-term leukemia risk and acute marrow toxicity following radiation exposure. The total shallow marrow (TM(50)), defined as all tissues lying within the first 50 µm of the bone surfaces, is considered to be the radiation target tissue of relevance for radiogenic bone cancer induction. For irradiation by sources external to the body, kerma to homogeneous spongiosa has been used as a surrogate for absorbed dose to both of these tissues, as direct dose calculations are not possible using computational phantoms with homogenized spongiosa. Recent micro-CT imaging of a 40 year old male cadaver has allowed for the accurate modeling of the fine microscopic structure of spongiosa in many regions of the adult skeleton (Hough et al 2011 Phys. Med. Biol. 56 2309-46). This microstructure, along with associated masses and tissue compositions, was used to compute specific absorbed fraction (SAF) values for protons originating in axial and appendicular bone sites (Jokisch et al 2011 Phys. Med. Biol. 56 6857-72). These proton SAFs, bone masses, tissue compositions and proton production cross sections, were subsequently used to construct neutron dose-response functions (DRFs) for both AM and TM(50) targets in each bone of the reference adult male. Kerma conditions were assumed for other resultant charged particles. For comparison, AM, TM(50) and spongiosa kerma coefficients were also calculated. At low incident neutron energies, AM kerma coefficients for neutrons correlate well with values of the AM DRF, while total marrow (TM) kerma coefficients correlate well with values of the TM(50) DRF. At high incident neutron energies, all kerma coefficients and DRFs tend to converge as charged-particle equilibrium is established across the bone site. In the range of 10 eV to 100 MeV, substantial differences are observed among the kerma coefficients and DRF. As a result, it is recommended that the AM kerma coefficient be used to estimate the AM DRF, and that the TM kerma coefficient be used to estimate the TM(50) DRF below 10 eV. Between 10 eV and 100 MeV, the appropriate DRF should be used as presented in this study. Above 100 MeV, spongiosa kerma coefficients apply well for estimating skeletal tissue doses. DRF values for each bone site as a function of energy are provided in an electronic annex to this article available at http://stacks.iop.org/0031-9155/56/6873/mmedia.

A comparative study of space radiation organ doses and associated cancer risks using PHITS and HZETRN

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 (<8%) in the effective dose and REID are found. However, for the galactic cosmic ray (GCR) boundary condition, compensating errors are observed, indicating that comparisons between the integral measurements of complex radiation environments and code calculations can be misleading. Code-to-code benchmarks allow for the comparison of differential quantities, such as secondary particle differential fluence, to provide insight into differences observed in integral quantities for particular components of the GCR spectrum.

Dosimetric impacts of microgravity: an analysis of 5th, 50th and 95th percentile male and female astronauts

Computational phantoms serve an important role in organ dosimetry and risk assessment performed at the National Aeronautics and Space Administration (NASA). A previous study investigated the impact on organ dose equivalents and effective doses from the use of the University of Florida hybrid adult male (UFHADM) and adult female (UFHADF) phantoms at differing height and weight percentiles versus those given by the two existing NASA phantoms, the computerized anatomical man (CAM) and female (CAF) (Bahadori et al 2011 Phys. Med. Biol. 56 1671-94). In the present study, the UFHADM and UFHADF phantoms of different body sizes were further altered to incorporate the effects of microgravity. Body self-shielding distributions are generated using the voxel-based ray tracer (VoBRaT), and the results are combined with depth dose data from the NASA codes BRYNTRN and HZETRN to yield organ dose equivalents and their rates for a variety of space radiation environments. It is found that while organ dose equivalents are indeed altered by the physiological effects ofmicrogravity, the magnitude of the change in overall risk (indicated by the effective dose) is minimal for the spectra and simplified shielding configurations considered. The results also indicate, however, that UFHADMand UFHADF could be useful in designing dose reduction strategies through optimized positioning of an astronaut during encounters with solar particle events.