All about MAX: a male adult voxel phantom for Monte Carlo calculations in radiation protection dosimetry

Monte Carlo modeling and simulation of a new 3D printed phantom for WBC calibration with ballistic gel as a tissue substitute

Whole-body counter (WBC) systems are used for in vivo monitoring in occupational internal dosimetry, typically calibrated using physical anthropomorphic phantoms. Our research group previously 3D-printed the Reference Female Phantom for Internal Dosimetry (RFPID) without internal organs specifically designed for WBC calibration. The RFPID and it is intended to fill it homogenously with ballistic gel, which is commonly used as a tissue equivalent in ballistic studies. However, comprehensive characterization of its physicochemical properties and radiological behavior as a tissue surrogate for dosimetry is limited. This study aims to evaluate the suitability of ballistic gel as a tissue substitute for physical phantoms in WBC system calibration and to analyze the RFPID as a model for WBC calibration. Ballistic gel tests determined its density and attenuation coefficients, comparing it to muscle, water, and PMMA. The RFPID was modeled and simulated using MCNP6.2 code and placed in an in vivo monitoring system using an 8"x4″ NaI(Tl) scintillator detector previously validated. The simulations were repeated with the RCP_AF of ICRP-110. Results indicate that ballistic gel has a density approximately 6% different from muscle and shows similar linear attenuation coefficients to muscle at intermediate and high energy levels (186-2200 keV). Simulations revealed a disparity of less than 9% in counting efficiency between RFPID and RCP_AF for energies from 100 to 3000 keV, confirming the phantom's suitability for WBC calibration and ballistic gel's viability as a tissue surrogate in internal dosimetry.

Considering clonal hematopoiesis of indeterminate potential in space radiation risk analysis for hematologic cancers and cardiovascular disease

BACKGROUND

Expanding human presence in space through long-duration exploration missions and commercial space operations warrants improvements in approaches for quantifying crew space radiation health risks. Currently, risk assessment models for radiogenic cancer and cardiovascular disease consider age, sex, and tobacco use, but do not incorporate other modifiable (e.g., body weight, physical activity, diet, environment) and non-modifiable individual risk factors (e.g., genetics, medical history, race/ethnicity, family history) that may greatly influence crew health both in-mission and long-term. For example, clonal hematopoiesis of indeterminate potential (CHIP) is a relatively common age-related condition that is an emerging risk factor for a variety of diseases including cardiovascular disease and cancer. CHIP carrier status may therefore exacerbate health risks associated with space radiation exposure.

METHODS

In the present study, published CHIP hazard ratios were used to modify background hazard rates for coronary heart disease, stroke, and hematologic cancers in the National Aeronautics and Space Administration space radiation risk assessment model. The risk of radiation exposure-induced death for these endpoints was projected for a future Mars exploration mission scenario.

RESULTS

Here we show appreciable increases in the lifetime risk of exposure-induced death for hematologic malignancies, coronary heart disease, and stroke, which are observed as a function of age after radiation exposure for male and female crew members that are directly attributable to the elevated health risks for CHIP carriers.

CONCLUSIONS

We discuss the importance of evaluating individual risk factors such as CHIP as part of a comprehensive space radiation risk assessment strategy aimed at effective risk communication and disease surveillance for astronauts embarking on future exploration missions.

Skeletal dosimetry in a microCT image-based rat model for external photon irradiation

Human skeletal dosimetry has experienced great developments in radiation protection in recent years by using the heterogeneous skeletal model. While for the rats experimentally used in radiation medicine, the investigation on skeletal dosimetry were mainly based on the homogeneous skeletal model, leading to inaccurate assessments of dose to radiosensitive tissues of red bone marrow (RBM) and bone surface. The purpose of this study is to develop a rat model with heterogeneous skeletal system and to investigate the dose difference in bone tissues for external photon irradiation. The high resolution of microCT images of a rat weighing 335 g were segmented into bone cortical, bone trabecular, bone marrow as well as other organs to construct the rat model. The absorbed dose to bone cortical, bone trabecular and bone marrow were calculated respectively by using Monte Carlo simulation for 22 external monoenergetic photon beams between 10 keV and 10 MeV under four different irradiation geometries conditions (left lateral [LL], right lateral [RL], dorsal-ventral [DV], ventral-dorsal [VD]). The calculated absorbed dose data were expressed as dose conversion coefficients and presented in this article, and the effect of irradiation conditions, photon energies and bone tissues density on the skeletal dose was discussed. The results showed that the dose conversion coefficients varying the photon energy for bone cortical, bone trabecular and bone marrow exhibit different trends and have the same sensitivity to irradiation conditions. The dose difference in bone tissues indicated that bone cortical and bone trabecular have significant attenuation effect on the energy deposition in bone marrow and bone surface for photon energies below 0.2 MeV. The set of dose conversion coefficients in this work can be used to determine the absorbed dose to skeletal system for external photon irradiation and to supplement the rat skeletal dosimetry.

Voxel model of a rabbit: assessment of absorbed doses in organs after CT examination performed by two different protocols

The objective of this work was to assess absorbed doses in organs and tissues of a rabbit, following computed tomography (CT) examinations, using a dedicated 3D voxel model. Absorbed doses in relevant organs were calculated using the MCNP5 Monte Carlo software. Calculations were perfomed for two standard CT protocols, using tube voltages of 110 kVp and 130 kVp. Absorbed doses were calculated in 11 organs and tissues, i.e., skin, bones, brain, muscles, heart, lungs, liver, spleen, kidney, testicles, and fat tissue. The doses ranged from 15.3 to 28.3 mGy, and from 40.2 to 74.3 mGy, in the two investigated protocols. The organs that received the highest dose were bones and kidneys. In contrast, brain and spleen were organs that received the smallest doses. Doses in organs which are stretched along the body did not change significantly with distance. On the other hand, doses in organs which are localized in the body showed maximums and minimums. Using the voxel model, it is possible to calculate the dose distribution in the rabbit's body after CT scans, and study the potential biological effects of CT doses in certain organs. The voxel model presented in this work can be used to calculated doses in all radiation experiments in which rabbits are used as experimental animals.

A tool for precise calculation of organ doses in voxelised geometries using GAMOS/Geant4 with a graphical user interface

Introduction: The limit of the method of calculating organ doses using voxelised phantoms with a Monte Carlo simulation code is that dose calculation errors in the boundaries of the organs are especially relevant for thin, small or complex geometries. In this report, we describe a tool that helps overcome this problem, accurately calculating organ doses by applying the “parallel geometry” utility feature of Geant4 through the GAMOS framework.

Methods and methods: We have tried to simplify the use of this tool by automatically processing the different DICOM image modalities (CT, PT, ST, NM), and by including the automatic conversion of the structures found in a DICOM RTSTRUCT file into Geant4 volumes that build the parallel geometry. For Nuclear Medicine applications, the DICOM PT, ST or NM images are converted into probabilities of generation of primary particles in each voxel, and the DICOM CT images into materials and material densities. For radiotherapy treatments, the DICOM RTPlan or RTIonPlan may also be used, hence the user only needs to describe the accelerator geometry. We also provide a Graphical User Interface for ease of use by for inexperienced users in Monte Carlo.

Results: We have tested the functionality of the tool with an I-131 thyroid cancer treatment, and obtained the expected energy deposition and dose differences, given that the particle source, geometry and structures are defined.

Conclusions: In summary, we provide an easy-to-use tool to calculate, with high accuracy, organ doses, taking into account their exact geometry as painted by the medical personnel on a voxelised phantom.

CRaTER observations and permissible mission duration for human operations in deep space

Prolonged exposure to the galactic cosmic ray (GCR) environment is a potentially limiting factor for manned missions in deep space. Evaluating the risk associated with the expected GCR environment is an essential step in planning a deep space mission. This requires an understanding of how the local interstellar spectrum is modulated by the heliospheric magnetic field (HMF) and how observed solar activity is manifested in the HMF over time. While current GCR models agree reasonably well with measured observations of GCR flux on the first matter, they must rely on imperfect or loose correlations to describe the latter. It is more accurate to use dose rates directly measured by instruments in deep space to quantify the GCR condition for a given period of time. In this work, dose rates observed by the Cosmic Ray Telescope for the Effects of Radiation (CRaTER) instrument are used to obtain the local GCR intensity and composition as a function of time. A response function is constructed that relates observed dose rates to solar modulation potential using a series of Monte Carlo radiation transport calculations. The record of observed solar modulation potential vs. time is then used to calculate a recent historical record of permissible mission duration (PMD) according to NASA's permissible exposure limits (PEL). Tables are provided for extreme values of PMD. Additional tables include risk of exposure-induced death (at upper 95% confidence interval) accrual rates and NASA effective dose rates as a function of solar modulation potential, astronaut age, sex, and shielding thickness. The significance of the PMD values reported in relation to likely transit duration requirements for future exploration missions is discussed. There is general agreement between CRaTER observations and the prescription of solar modulation vs. time given by the Badhwar-O'Neill 2014 GCR model. However, CRaTER observations do capture the effects of significant heliospheric transients, among other features, that are missing from the prescription of solar modulation potential vs. time.

NASA's first ground-based Galactic Cosmic Ray Simulator: Enabling a new era in space radiobiology research

With exciting new NASA plans for a sustainable return to the moon, astronauts will once again leave Earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic radiation (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPEs with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include carcinogenesis, central nervous system (CNS) effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including circulatory and heart disease, as well as potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy-ion exposures and the validation of countermeasures in a relevant space environment. Historically, most research on understanding space radiation-induced health risks has been performed using acute exposures of monoenergetic single-ion beams. However, the space radiation environment consists of a wide variety of ion species over a broad energy range. Using the fast beam switching and controls systems technology recently developed at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory, a new era in radiobiological research is possible. NASA has developed the “GCR Simulator” to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (approximately 65%–75%) and helium ions (approximately 10%–20%) with heavier ions (Z ≥ 3) contributing the remainder. The GCR simulator exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the 5 heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader system is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low-energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires approximately 75 minutes. To more closely simulate the low-dose rates found in space, sequential field exposures can be divided into daily fractions over 2 to 6 weeks, with individual beam fractions as low as 0.1 to 0.2 mGy. In the large beam configuration (60 × 60 cm²), 54 special housing cages can accommodate 2 to 3 mice each for an approximately 75 min duration or 15 individually housed rats. On June 15, 2018, the NSRL made a significant achievement by completing the first operational run using the new GCR simulator. This paper discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NSRL to accelerate our understanding and mitigation of health risks faced by astronauts. Ultimately, the GCR simulator will require validation across multiple radiogenic risks, endpoints, doses, and dose rates.

This study describes how NASA’s new earth-based galactic cosmic ray simulator is being used to accelerate our understanding of the effects of space radiation exposure on astronauts and to validate countermeasures for exploration missions. For the first time, research teams can study mixed field ion and dose rate effects in a simulated space environment.

A methodology for investigating the impact of medical countermeasures on the risk of exposure induced death

The space radiation environment is composed of ionizing particles that may pose health risks to crew members during Low Earth Orbit (LEO) and deep space missions. NASA has established astronaut career radiation limits for cancer of 3% Risk of Exposure Induced Death (REID) at the 95% confidence level. The REID is the increased lifetime risk of death from cancer due to radiation exposure in comparison to an unexposed background population and has been traditionally mitigated by passive shielding design concepts and limiting safe days in space. Additional reduction in radiation exposure risk may be achieved with Medical Countermeasures (MCM). Recent meta-analyses have demonstrated the efficacy of aspirin in the reduction of the background colorectal cancer incidence and mortality rates for specific cohorts. Additional studies of warfarin in patients greater than 50 years of age have indicated statistically significant decreases in stomach, bladder, brain, prostate, and lung cancer incidence as compared to control groups. While ultimate selection of suitable countermeasures will be the responsibility of flight surgeons, this paper presents a general methodology for incorporating MCM into the NASA Space Radiation Cancer Risk model and includes modifications of the background mortality rates (hazard rates) and the radiation risk coefficients to numerically quantify the benefits of MCM. As examples of the method, aspirin and warfarin will be employed as MCM in a sensitivity analysis to compute the REID for astronauts embarking on a one-year deep space mission scenario.

Whole body effective dose equivalent dataset for MAX and FAX shielded with Common Aerospace Materials in deep space

Materials have a primary purpose in the design of space vehicles, such as fuels, walls, racks, windows, etc. Additionally, each will also affect space radiation protection. Using the On-Line Tool for the Assessment of Radiation in Space (OLTARIS), version 3.5, analysis package, this article includes the whole body effective dose equivalent (E_D) data from human phantoms being shielded by 59 aerospace materials for deep space travel. To represent the average anatomy of an astronaut, the Female Adult voXel (FAX), 2005 version, and the Male Adult voXel (MAX), 2005 version, human phantoms are used. A simple spherical geometry, which is composed of a spherical shell with the human phantom placed in the center, is also used. Eighteen shielding thicknesses ranging from 0.01 to 1000 g per centimetres squared are evaluated and the ray distribution used in this study is the 1002 geodesic. All aerospace materials are categorized into four groups: metals, polymers, composites, and fuels, hydrides, and liquid gases. These materials include common fuels and propellants used in space travel, engineered materials developed to significantly increase the absorption of secondary radiation, and materials in the early stages of development for the purpose of meeting both shielding and structural needs of future spacecraft missions. The data in this article is used for the paper, "Evaluating the Effectiveness of Common Aerospace Materials at Lowering the Whole Body Effective Dose Equivalent in Deep Space," [13].

Anthropomorphic left ventricular mesh phantom: a framework to investigate the accuracy of SQUEEZ using Coherent Point Drift for the detection of regional wall motion abnormalities

We present an anthropomorphically accurate left ventricular (LV) phantom derived from human computed tomography (CT) data to serve as the ground truth for the optimization and the spatial resolution quantification of a CT-derived regional strain metric (SQUEEZ) for the detection of regional wall motion abnormalities. Displacements were applied to the mesh points of a clinically derived end-diastolic LV mesh to create analytical end-systolic poses with physiologically accurate endocardial strains. Normal function and regional dysfunction of four sizes [1, 2/3, 1/2, and 1/3 American Heart Association (AHA) segments as core diameter], each exhibiting hypokinesia (70% reduction in strain) and subtle hypokinesia (40% reduction in strain), were simulated. Regional shortening (RS CT ) estimates were obtained by registering the end-diastolic mesh to each simulated end-systolic mesh condition using a nonrigid registration algorithm. Ground-truth models of normal function and of hypokinesia were used to identify the optimal parameters in the registration algorithm and to measure the accuracy of detecting regional dysfunction of varying sizes and severities. For normal LV function,RS CT values in all 16 AHA segments were accurate to within ± 5 % . For cases with regional dysfunction, the errors inRS CT around the dysfunctional region increased with decreasing size of dysfunctional tissue.

Advances in Computational Human Phantoms and Their Applications in Biomedical Engineering - A Topical Review

Over the past decades, significant improvements have been made in the field of computational human phantoms (CHPs) and their applications in biomedical engineering. Their sophistication has dramatically increased. The very first CHPs were composed of simple geometric volumes, e.g., cylinders and spheres, while current CHPs have a high resolution, cover a substantial range of the patient population, have high anatomical accuracy, are poseable, morphable, and are augmented with various details to perform functionalized computations. Advances in imaging techniques and semi-automated segmentation tools allow fast and personalized development of CHPs. These advances open the door to quickly develop personalized CHPs, inherently including the disease of the patient. Because many of these CHPs are increasingly providing data for regulatory submissions of various medical devices, the validity, anatomical accuracy, and availability to cover the entire patient population is of utmost importance. The article is organized into two main sections: the first section reviews the different modeling techniques used to create CHPs, whereas the second section discusses various applications of CHPs in biomedical engineering. Each topic gives an overview, a brief history, recent developments, and an outlook into the future.

Solar particle event storm shelter requirements for missions beyond low Earth orbit

Protecting spacecraft crews from energetic space radiations that pose both chronic and acute health risks is a critical issue for future missions beyond low Earth orbit (LEO). Chronic health risks are possible from both galactic cosmic ray and solar energetic particle event (SPE) exposures. However, SPE exposures also can pose significant short term risks including, if dose levels are high enough, acute radiation syndrome effects that can be mission- or life-threatening. In order to address the reduction of short term risks to spaceflight crews from SPEs, we have developed recommendations to NASA for a design-standard SPE to be used as the basis for evaluating the adequacy of proposed radiation shelters for cislunar missions beyond LEO. Four SPE protection requirements for habitats are proposed: (1) a blood-forming-organ limit of 250 mGy-equivalent for the design SPE; (2) a design reference SPE environment equivalent to the sum of the proton spectra during the October 1989 event series; (3) any necessary assembly of the protection system must be completed within 30 min of event onset; and (4) space protection systems must be designed to ensure that astronaut radiation exposures follow the ALARA (As Low As Reasonably Achievable) principle.

Modeling "Textured" Bones in Virtual Human Phantoms

The purpose of this study was to develop detailed and realistic models of the cortical and trabecular bones in the spine, ribs, and sternum and incorporate them into the library of virtual human phantoms (XCAT). Cortical bone was modeled by 3D morphological erosion of XCAT homogenously defined bones with an average thickness measured from the CT dataset upon which each individual XCAT phantom was based. The trabecular texture was modeled using a power law synthesis algorithm where the parameters were tuned using high-resolution anatomical images of the Human Visible Female. The synthesized bone textures were added into the XCAT phantoms. To qualitatively evaluate the improved realism of the bone modeling, CT simulations of the XCAT phantoms were acquired with and without the textured bone modeling. The 3D power spectrum of the anatomical images exhibited a power law behavior (R² = 0.84), as expected in fractal and porous textures. The proposed texture synthesis algorithm was able to synthesize textures emulating real anatomical images, with the simulated CT images with the prototyped bones were more realistic than those simulated with the original XCAT models. Incorporating intra-organ structures, the "textured" phantoms are envisioned to be used to conduct virtual clinical trials in the context of medical imaging in cases where the actual trials are infeasible due to the lack of ground truth, cost, or potential risks to the patients.

Computational phantoms, ICRP/ICRU, and further developments

Phantoms simulating the human body play a central role in radiation dosimetry. The first computational body phantoms were based upon mathematical expressions describing idealised body organs. With the advent of more powerful computers in the 1980s, voxel phantoms have been developed. Being based on three-dimensional images of individuals, they offer a more realistic anatomy. Hence, the International Commission on Radiological Protection (ICRP) decided to construct voxel phantoms representative of the adult Reference Male and Reference Female for the update of organ dose coefficients. Further work on phantom development has focused on phantoms that combine the realism of patient-based voxel phantoms with the flexibility of mathematical phantoms, so-called 'boundary representation' (BREP) phantoms. This phantom type has been chosen for the ICRP family of paediatric reference phantoms. Due to the limited voxel resolution of the adult reference computational phantoms, smaller tissues, such as the lens of the eye, skin, and micron-thick target tissues in the respiratory and alimentary tract regions, could not be segmented properly. In this context, ICRP Committee 2 initiated a research project with the goal of producing replicas of the ICRP Publication 110 phantoms in polygon mesh format, including all source and target regions, even those with micron resolution. BREP phantoms of the fetus and the pregnant female at various stages of gestation complete the phantoms available for radiation protection computations.

MCNPX CALCULATIONS OF SPECIFIC ABSORBED FRACTIONS IN SOME ORGANS OF THE HUMAN BODY DUE TO APPLICATION OF 133Xe, 99mTc and 81mKr RADIONUCLIDES

Monte Carlo simulations were performed to evaluate treatment doses with wide spread used radionuclides 133Xe, 99mTc and 81mKr. These different radionuclides are used in perfusion or ventilation examinations in nuclear medicine and as indicators for cardiovascular and pulmonary diseases. The objective of this work was to estimate the specific absorbed fractions in surrounding organs and tissues, when these radionuclides are incorporated in the lungs. For this purpose a voxel thorax model has been developed and compared with the ORNL phantom. All calculations and simulations were performed by means of the MCNP5/X code.

Development of a paediatric head voxel model database for dosimetric applications

OBJECTIVE

To develop a database of paediatric head voxel models intended for Monte Carlo (MC) dosimetric applications.

METHODS

Seventeen head and neck CT image data sets were retrieved from the picture archiving and communicating system of our hospital and were reformed into voxel models. 22 organs were segmented at each data set. The segmented organ masses were compared to the respective age- and gender-specific ICRP reference mass value. Adjustments were made such that segmented and reference mass values coincide within a tolerance of 10%. A dental cone beam CT cleft palate simulation study was set up to demonstrate the applicability of our database to MC frameworks and to investigate the need for age- and gender-specific paediatric models.

RESULTS

The designed database covers the age range from 2 months to 14 years old. Each model represents a reference head voxel phantom for its corresponding age and gender category. The simulation study revealed absorbed organ dose differences larger than 50% among the 5, 8 and 12 years old models when exposed to identical conditions.

CONCLUSION

Children cannot be represented by one average phantom covering the entire age range like adults due to the fact that their organs change rapidly in size and shape. A database of paediatric head voxel models was designed to enable dose calculations via MC simulations. Advances in knowledge: The application of each model of the database to MC frameworks provides age- and gender-specific organ dose estimations from medical exposures in the head and neck region.

COMPARISON OF ORGAN DOSES IN HUMAN PHANTOMS: VARIATIONS DUE TO BODY SIZE AND POSTURE

Organ dose calculations performed using human phantoms can provide estimates of astronauts' health risks due to cosmic radiation. However, the characteristics of such phantoms strongly affect the estimation precision. To investigate organ dose variations with body size and posture in human phantoms, a non-uniform rational B-spline boundary surfaces model was constructed based on cryosection images. This model was used to establish four phantoms with different body size and posture parameters, whose organs parameters were changed simultaneously and which were voxelised with 4 × 4 × 4 mm3 resolution. Then, using Monte Carlo transport code, the organ doses caused by ≤500 MeV isotropic incident protons were calculated. The dose variations due to body size differences within a certain range were negligible, and the doses received in crouching and standing-up postures were similar. Therefore, a standard Chinese phantom could be established, and posture changes cannot effectively protect astronauts during solar particle events.

Converting computed tomography images into photon interaction coefficients by using stoichiometric calibration and parametric fit models

PURPOSE

X ray and γ-ray are widely applied in radiology, radiotherapy, and nuclear medicine. Linear attenuation coefficients and linear energy absorption coefficients are essential for dose calculation and image correction. In this study, a method that entails combining the stoichiometric calibration and parametric physical models was developed to convert computed tomography (CT) images into the linear attenuation coefficients and linear energy absorption coefficients.

METHODS

A calibration scan was performed using standard tissue-equivalent materials to obtain the characteristics of the x-ray energy spectrum. Subsequently, relationships between CT numbers and tissue parameters were established using standard soft tissue and bone tissue data adopted from the literature. The linear attenuation coefficient and linear energy absorption coefficient were calculated using the parametric fit model.

RESULTS

The results showed a linear relationship between CT numbers and tissue parameters. The tissue-equivalent materials differed from real human tissues, leading to considerable errors in estimation of mass attenuation coefficients when the photon energy was lower than 50 keV. Mass attenuation coefficients and mass energy transfer coefficients of five tissues were calculated and validated using clinical CT images. The error was less than ± 5% and ± 8%, compared with the values of the International Commission on Radiation Units (ICRU) 46 report.

CONCLUSIONS

The probability of photon interaction with tissues and physical characteristics of tissues can be accurately evaluated by using the proposed method and applied in various clinical applications.

Rotating and translating anthropomorphic head voxel models to establish an horizontal Frankfort plane for dental CBCT Monte Carlo simulations: a dose comparison study

In order to carry out Monte Carlo (MC) dosimetry studies, voxel phantoms, modeling human anatomy, and organ-based segmentation of CT image data sets are applied to simulation frameworks. The resulting voxel phantoms preserve patient CT acquisition geometry; in the case of head voxel models built upon head CT images, the head support with which CT scanners are equipped introduces an inclination to the head, and hence to the head voxel model. In dental cone beam CT (CBCT) imaging, patients are always positioned in such a way that the Frankfort line is horizontal, implying that there is no head inclination. The orientation of the head is important, as it influences the distance of critical radiosensitive organs like the thyroid and the esophagus from the x-ray tube. This work aims to propose a procedure to adjust head voxel phantom orientation, and to investigate the impact of head inclination on organ doses in dental CBCT MC dosimetry studies. The female adult ICRP, and three in-house-built paediatric voxel phantoms were in this study. An EGSnrc MC framework was employed to simulate two commonly used protocols; a Morita Accuitomo 170 dental CBCT scanner (FOVs: 60  ×  60 mm² and 80  ×  80 mm², standard resolution), and a 3D Teeth protocol (FOV: 100  ×  90 mm²) in a Planmeca Promax 3D MAX scanner. Result analysis revealed large absorbed organ dose differences in radiosensitive organs between the original and the geometrically corrected voxel models of this study, ranging from  -45.6% to 39.3%. Therefore, accurate dental CBCT MC dose calculations require geometrical adjustments to be applied to head voxel models.

MIDA: A Multimodal Imaging-Based Detailed Anatomical Model of the Human Head and Neck

Computational modeling and simulations are increasingly being used to complement experimental testing for analysis of safety and efficacy of medical devices. Multiple voxel- and surface-based whole- and partial-body models have been proposed in the literature, typically with spatial resolution in the range of 1-2 mm and with 10-50 different tissue types resolved. We have developed a multimodal imaging-based detailed anatomical model of the human head and neck, named "MIDA". The model was obtained by integrating three different magnetic resonance imaging (MRI) modalities, the parameters of which were tailored to enhance the signals of specific tissues: i) structural T1- and T2-weighted MRIs; a specific heavily T2-weighted MRI slab with high nerve contrast optimized to enhance the structures of the ear and eye; ii) magnetic resonance angiography (MRA) data to image the vasculature, and iii) diffusion tensor imaging (DTI) to obtain information on anisotropy and fiber orientation. The unique multimodal high-resolution approach allowed resolving 153 structures, including several distinct muscles, bones and skull layers, arteries and veins, nerves, as well as salivary glands. The model offers also a detailed characterization of eyes, ears, and deep brain structures. A special automatic atlas-based segmentation procedure was adopted to include a detailed map of the nuclei of the thalamus and midbrain into the head model. The suitability of the model to simulations involving different numerical methods, discretization approaches, as well as DTI-based tensorial electrical conductivity, was examined in a case-study, in which the electric field was generated by transcranial alternating current stimulation. The voxel- and the surface-based versions of the models are freely available to the scientific community.

Fluence-to-effective dose conversion coefficients from a Saudi population based phantom for monoenergetic photon beams from 10 keV to 20 MeV

Fluence-to-dose conversion coefficients are important quantities for radiation protection, derived from Monte Carlo simulations of the radiation particles through a stylised phantom or voxel based phantoms. The voxel phantoms have been developed for many ethnic groups for their accurate reflection of the anatomy. In this study, we used the Monte Carlo code MCNPX to calculate the photon fluence-to-effective dose conversion coefficients with a voxel phantom based on the Saudi Arabian male population. Six irradiation geometries, anterior-posterior (AP), posterior-anterior (PA), left lateral (LLAT), right lateral (RLAT), rotational (ROT) and isotropic (ISO) were simulated for monoenergetic photon beams from 10 keV to 20 MeV. We compared the coefficients with the reference values in ICRP Publication 116. The coefficients in the AP and PA geometries match the reference values to 9% and 12% on average as measured by root mean square while those in the LLAT, RLAT ROT and ISO geometries differ, mostly below, from the reference by 23, 22, 15 and 16%, respectively. The torso of the Saudi phantom is wider than the ICRP reference male phantom and likely to cause more attenuation to the lateral beam. The ICRP reference coefficients serve well for the Saudi male population as conservative estimations for the purpose of radiation protection.

Monte Carlo simulation of a whole-body counter using IGOR phantoms

Whole-body counting is a technique of choice for assessing the intake of gamma-emitting radionuclides. An appropriate calibration is necessary, which is done either by experimental measurement or by Monte Carlo (MC) calculation. The aim of this work was to validate a MC model for calibrating whole-body counters (WBCs) by comparing the results of computations with measurements performed on an anthropomorphic phantom and to investigate the effect of a change in phantom's position on the WBC counting sensitivity. GEANT MC code was used for the calculations, and an IGOR phantom loaded with several types of radionuclides was used for the experimental measurements. The results show a reasonable agreement between measurements and MC computation. A 1-cm error in phantom positioning changes the activity estimation by >2%. Considering that a 5-cm deviation of the positioning of the phantom may occur in a realistic counting scenario, this implies that the uncertainty of the activity measured by a WBC is ∼10-20%.

An exponential growth of computational phantom research in radiation protection, imaging, and radiotherapy: a review of the fifty-year history

Radiation dose calculation using models of the human anatomy has been a subject of great interest to radiation protection, medical imaging, and radiotherapy. However, early pioneers of this field did not foresee the exponential growth of research activity as observed today. This review article walks the reader through the history of the research and development in this field of study which started some 50 years ago. This review identifies a clear progression of computational phantom complexity which can be denoted by three distinct generations. The first generation of stylized phantoms, representing a grouping of less than dozen models, was initially developed in the 1960s at Oak Ridge National Laboratory to calculate internal doses from nuclear medicine procedures. Despite their anatomical simplicity, these computational phantoms were the best tools available at the time for internal/external dosimetry, image evaluation, and treatment dose evaluations. A second generation of a large number of voxelized phantoms arose rapidly in the late 1980s as a result of the increased availability of tomographic medical imaging and computers. Surprisingly, the last decade saw the emergence of the third generation of phantoms which are based on advanced geometries called boundary representation (BREP) in the form of Non-Uniform Rational B-Splines (NURBS) or polygonal meshes. This new class of phantoms now consists of over 287 models including those used for non-ionizing radiation applications. This review article aims to provide the reader with a general understanding of how the field of computational phantoms came about and the technical challenges it faced at different times. This goal is achieved by defining basic geometry modeling techniques and by analyzing selected phantoms in terms of geometrical features and dosimetric problems to be solved. The rich historical information is summarized in four tables that are aided by highlights in the text on how some of the most well-known phantoms were developed and used in practice. Some of the information covered in this review has not been previously reported, for example, the CAM and CAF phantoms developed in 1970s for space radiation applications. The author also clarifies confusion about 'population-average' prospective dosimetry needed for radiological protection under the current ICRP radiation protection system and 'individualized' retrospective dosimetry often performed for medical physics studies. To illustrate the impact of computational phantoms, a section of this article is devoted to examples from the author's own research group. Finally the author explains an unexpected finding during the course of preparing for this article that the phantoms from the past 50 years followed a pattern of exponential growth. The review ends on a brief discussion of future research needs (a supplementary file '3DPhantoms.pdf' to figure 15 is available for download that will allow a reader to interactively visualize the phantoms in 3D).

Study of the counting efficiency of a WBC setup by using a computational 3D human body library in sitting position based on polygonal mesh surfaces

A realistic computational 3D human body library, called MaMP and FeMP (Male and Female Mesh Phantoms), based on polygonal mesh surface geometry, has been created to be used for numerical calibration of the whole body counter (WBC) system of the nuclear power plant (NPP) in Doel, Belgium. The main objective was to create flexible computational models varying in gender, body height, and mass for studying the morphology-induced variation of the detector counting efficiency (CE) and reducing the measurement uncertainties. First, the counting room and an HPGe detector were modeled using MCNPX (Monte Carlo radiation transport code). The validation of the model was carried out for different sample-detector geometries with point sources and a physical phantom. Second, CE values were calculated for a total of 36 different mesh phantoms in a seated position using the validated Monte Carlo model. This paper reports on the validation process of the in vivo whole body system and the CE calculated for different body heights and weights. The results reveal that the CE is strongly dependent on the individual body shape, size, and gender and may vary by a factor of 1.5 to 3 depending on the morphology aspects of the individual to be measured.

[Development of a digital chest phantom for studies on energy subtraction techniques]

Digital chest phantoms continue to play a significant role in optimizing imaging parameters for chest X-ray examinations. The purpose of this study was to develop a digital chest phantom for studies on energy subtraction techniques under ideal conditions without image noise. Computed tomography (CT) images from the LIDC (Lung Image Database Consortium) were employed to develop a digital chest phantom. The method consisted of the following four steps: 1) segmentation of the lung and bone regions on CT images; 2) creation of simulated nodules; 3) transformation to attenuation coefficient maps from the segmented images; and 4) projection from attenuation coefficient maps. To evaluate the usefulness of digital chest phantoms, we determined the contrast of the simulated nodules in projection images of the digital chest phantom using high and low X-ray energies, soft tissue images obtained by energy subtraction, and "gold standard" images of the soft tissues. Using our method, the lung and bone regions were segmented on the original CT images. The contrast of simulated nodules in soft tissue images obtained by energy subtraction closely matched that obtained using the gold standard images. We thus conclude that it is possible to carry out simulation studies based on energy subtraction techniques using the created digital chest phantoms. Our method is potentially useful for performing simulation studies for optimizing the imaging parameters in chest X-ray examinations.

Monte Carlo modeling for dose assessment in cone beam CT for oral and maxillofacial applications

PURPOSE

To adjust Monte Carlo modeling for dose assessment in dedicated cone beam computed tomography (CBCT) of the oral and maxillofacial region.

METHODS

Two different CBCT systems with different fields of view (FOVs), beam qualities, and scan geometries were modeled using Monte Carlo simulation. Dose calculations for typical CBCT examinations were performed with the head and neck part of four computational anatomical phantoms.

RESULTS

Simulation results compared favorably to values acquired experimentally using physical phantoms (in the literature). For a given phantom scanned with 90 kV and 60 × 60 mm FOV, effective dose per mAs was on average the same for the two different systems. Exposing the four phantoms under identical settings for the same CBCT system resulted in variations in organ doses of greater than 100%, leading to differences in effective dose of 30%. For one system, the dose dependence on the operating tube potential can be described with a quadratic polynomial function. Dose distributions over the axial plane were presented as contour plots.

CONCLUSIONS

Monte Carlo modeling is an efficient and accurate means of evaluating dose distributions for dedicated cone beam oral and maxillofacial CT. Results suggest large dose differences among patients undergoing the same examination on a given system, supporting approaches toward patient specific dosimetry. A dedicated and standardized computational phantom for head and neck dosimetry should be established.

Development of newborn and 1-year-old reference phantoms based on polygon mesh surfaces

The purpose of this study is the development of paediatric reference phantoms for newborn and 1-year-old infants to be used for the calculation of organ and tissue equivalent doses in radiation protection. The study proposes a method for developing anatomically highly sophisticated paediatric phantoms without using medical images. The newborn and 1-year-old hermaphrodite phantoms presented here were developed using three-dimensional (3D) modelling software applied to anatomical information taken from atlases, textbooks and images provided by the Department of Anatomy of the Federal University of Pernambuco, Brazil. The method uses polygon mesh surfaces to model body contours, the shape of organs as well as their positions and orientations in the human body. Organ and tissue masses agree with corresponding data given by the International Commission on Radiological Protection for newborn and 1-year-old reference children. Bones were segmented into cortical bone, spongiosa, medullary marrow and cartilage to allow for the use of μCT images of trabecular bone for skeletal dosimetry. Anatomical results show 3D images of the phantoms' surfaces, organs and skeletons, as well as tables with organ and tissue masses or skeletal tissue volumes. Dosimetric results present comparisons of organ and tissue absorbed doses or specific absorbed fractions between the newborn and 1-year-old phantoms and corresponding data for other paediatric stylised or voxel phantoms. Most differences were found to be below 10%.

The properties of human body phantoms used in calculations of electromagnetic fields exposure by wireless communication handsets or hand-operated industrial devices

According to international guidelines, the assessment of biophysical effects of exposure to electromagnetic fields (EMF) generated by hand-operated sources needs the evaluation of induced electric field (E(in)) or specific energy absorption rate (SAR) caused by EMF inside a worker's body and is usually done by the numerical simulations with different protocols applied to these two exposure cases. The crucial element of these simulations is the numerical phantom of the human body. Procedures of E(in) and SAR evaluation due to compliance analysis with exposure limits have been defined in Institute of Electrical and Electronics Engineers standards and International Commission on Non-Ionizing Radiation Protection guidelines, but a detailed specification of human body phantoms has not been described. An analysis of the properties of over 30 human body numerical phantoms was performed which has been used in recently published investigations related to the assessment of EMF exposure by various sources. The differences in applicability of these phantoms in the evaluation of E(in) and SAR while operating industrial devices and SAR while using mobile communication handsets are discussed. The whole human body numerical phantom dimensions, posture, spatial resolution and electric contact with the ground constitute the key parameters in modeling the exposure related to industrial devices, while modeling the exposure from mobile communication handsets, which needs only to represent the exposed part of the human body nearest to the handset, mainly depends on spatial resolution of the phantom. The specification and standardization of these parameters of numerical human body phantoms are key requirements to achieve comparable and reliable results from numerical simulations carried out for compliance analysis against exposure limits or within the exposure assessment in EMF-related epidemiological studies.

RADAR reference adult, pediatric, and pregnant female phantom series for internal and external dosimetry

A new generation of reference computational phantoms, based on image-based models tied to the reference masses defined by the International Commission on Radiological Protection (ICRP) for dose calculations, is presented.

METHODS

Anatomic models based on nonuniform rational b-spline modeling techniques were used to define reference male and female adults, 15-y-olds, 10-y-olds, 5-y-olds, 1-y-olds, newborns, and pregnant women at 3 stages of gestation, using the defined reference organ masses in ICRP publication 89. Absorbed fractions and specific absorbed fractions for internal emitters were derived using standard Monte Carlo radiation transport simulation codes.

RESULTS

Differences were notable between many pairs of organs in specific absorbed fractions because of the improved realism of the models, with adjacent organs usually closer and sometimes touching. Final estimates of absorbed dose for radiopharmaceuticals, for example, were only slightly different overall, as many of the differences were small and most pronounced at low radiation energies. Some new important organs were defined (salivary glands, prostate, eyes, and esophagus), and the identity of a few gastrointestinal tract organs changed.

CONCLUSION

A new generation of reference models for standardized internal and external dose calculations has been defined. The models will be implemented in standardized software for internal dose calculations and be used to produce new standardized dose estimates for radiopharmaceuticals and other applications.

Construction of an extended library of adult male 3D models: rationale and results

In order to best cover the possible extent of heights and weights of male adults the construction of 25 whole body 3D models has been undertaken. Such a library is thought to be useful to specify the uncertainties and relevance of dosimetry calculations carried out with models representing individuals of average body heights and weights. Representative 3D models of Caucasian body types are selected in a commercial database according to their height and weight, and 3D models of the skeleton and internal organs are designed using another commercial dataset. A review of the literature enabled one to fix volume or mass target values for the skeleton, soft organs, skin and fat content of the selected individuals. The composition of the remainder tissue is fixed so that the weight of the voxel models equals the weight of the selected individuals. After mesh and NURBS modelling, volume adjustment of the selected body shapes and additional voxel-based work, 25 voxel models with 109 identified organs or tissue are obtained. Radiation transport calculations are carried out with some of the developed models to illustrate potential uses. The following points are discussed throughout this paper: justification of the fixed or obtained models’ features regarding available and relevant literature data; workflow and strategy for major modelling steps; advantages and drawbacks of the obtained library as compared with other works. The construction hypotheses are explained and justified in detail since future calculation results obtained with this library will depend on them.

Hybrid computational phantoms representing the reference adult male and adult female: construction and applications for retrospective dosimetry

Currently, two classes of computational phantoms have been developed for dosimetry calculation: (1) stylized (or mathematical) and (2) voxel (or tomographic) phantoms describing human anatomy through mathematical surface equations and 3D voxel matrices, respectively. Mathematical surface equations in stylized phantoms are flexible, but the resulting anatomy is not as realistic. Voxel phantoms display far better anatomical realism, but they are limited in terms of their ability to alter organ shape, position, and depth, as well as body posture. A new class of computational phantoms called hybrid phantoms takes advantage of the best features of stylized and voxel phantoms-flexibility and anatomical realism, respectively. In the current study, hybrid computational phantoms representing the adult male and female reference anatomy and anthropometry are presented. These phantoms serve as the starting framework for creating patient or worker sculpted whole-body phantoms for retrospective dose reconstruction. Contours of major organs and tissues were converted or segmented from computed tomography images of a 36-y-old Korean volunteer and a 25-y-old U.S. female patient, respectively, with supplemental high-resolution CT images of the cranium. Polygon mesh models for the major organs and tissues were reconstructed and imported into Rhinoceros™ for non-uniform rational B-spline (NURBS) surface modeling. The resulting NURBS/polygon mesh models representing body contour and internal anatomy were matched to anthropometric data and reference organ mass data provided by Centers for Disease Control and Prevention and International Commission on Radiation Protection, respectively. Finally, two hybrid adult male and female phantoms were completed where a total of eight anthropometric data categories were matched to standard values within 4% and organ volumes matched to ICRP data within 1% with the exception of total skin. The hybrid phantoms were voxelized from the NURBS phantoms at resolutions of 0.158 × 0.158 × 0.158 cm and 0.126 × 0.126 × 0.126 cm for the male and female, respectively. To highlight the flexibility of the hybrid phantoms, graphical displays are given of (1) underweight and overweight adult male phantoms, (2) a sitting position for the adult female phantom, and (3) extraction and higher-resolution voxelization of the small intestine for localized dosimetry of mucosal and stem cell layers. These phantoms are used to model radioactively contaminated individuals and to then assess time-dependent detector count rate thresholds corresponding to 50, 250, and 500 mSv effective dose, as might be needed during in-field radiological triage by first responders or first receivers.

Selected organ dose conversion coefficients for external photons calculated using ICRP adult voxel phantoms and Monte Carlo code FLUKA

The adult reference male and female computational voxel phantoms recommended by ICRP are adapted into the Monte Carlo transport code FLUKA. The FLUKA code is then utilised for computation of dose conversion coefficients (DCCs) expressed in absorbed dose per air kerma free-in-air for colon, lungs, stomach wall, breast, gonads, urinary bladder, oesophagus, liver and thyroid due to a broad parallel beam of mono-energetic photons impinging in anterior-posterior and posterior-anterior directions in the energy range of 15 keV-10 MeV. The computed DCCs of colon, lungs, stomach wall and breast are found to be in good agreement with the results published in ICRP publication 110. The present work thus validates the use of FLUKA code in computation of organ DCCs for photons using ICRP adult voxel phantoms. Further, the DCCs for gonads, urinary bladder, oesophagus, liver and thyroid are evaluated and compared with results published in ICRP 74 in the above-mentioned energy range and geometries. Significant differences in DCCs are observed for breast, testis and thyroid above 1 MeV, and for most of the organs at energies below 60 keV in comparison with the results published in ICRP 74. The DCCs of female voxel phantom were found to be higher in comparison with male phantom for almost all organs in both the geometries.

Development of 5- and 10-year-old pediatric phantoms based on polygon mesh surfaces

PURPOSE

The purpose of this study is the development of reference pediatric phantoms for 5- and 10-year-old children to be used for the calculation of organ and tissue equivalent doses in radiation protection.

METHODS

The study proposes a method for developing anatomically highly sophisticated pediatric phantoms without using medical images. The 5- and 10-year-old male and female phantoms presented here were developed using 3D modeling software applied to anatomical information taken from atlases and textbooks. The method uses polygon mesh surfaces to model body contours, the shape of organs as well as their positions, and orientations in the human body. Organ and tissue masses comply with the corresponding data given by the International Commission on Radiological Protection (ICRP) for the 5- and 10-year-old reference children. Bones were segmented into cortical bone, spongiosa, medullary marrow, and cartilage to allow for the use of micro computer tomographic (microCT) images of trabecular bone for skeletal dosimetry.

RESULTS

The four phantoms, a male and a female for each age, and their organs are presented in 3D images and their organ and tissue masses in tables which show the compliance of the ICRP reference values. Dosimetric data, calculated for the reference pediatric phantoms by Monte Carlo methods were compared with corresponding data from adult mesh phantoms and pediatric stylized phantoms. The comparisons show reasonable agreement if the anatomical differences between the phantoms are properly taken into account.

CONCLUSIONS

Pediatric phantoms were developed without using medical images of patients or volunteers for the first time. The models are reference phantoms, suitable for regulatory dosimetry, however, the 3D modeling method can also be applied to medical images to develop patient-specific phantoms.

Variation in lunar neutron dose estimates

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.

Organ-specific external dose coefficients and protective apron transmission factors for historical dose reconstruction for medical personnel

While radiation absorbed dose (Gy) to the skin or other organs is sometimes estimated for patients from diagnostic radiologic examinations or therapeutic procedures, rarely is occupationally-received radiation absorbed dose to individual organs/tissues estimated for medical personnel; e.g., radiologic technologists or radiologists. Generally, for medical personnel, equivalent or effective radiation doses are estimated for compliance purposes. In the very few cases when organ doses to medical personnel are reconstructed, the data is usually for the purpose of epidemiologic studies; e.g., a study of historical doses and risks to a cohort of about 110,000 radiologic technologists presently underway at the U.S. National Cancer Institute. While ICRP and ICRU have published organ-specific external dose conversion coefficients (DCCs) (i.e., absorbed dose to organs and tissues per unit air kerma and dose equivalent per unit air kerma), those factors have been published primarily for mono-energetic photons at selected energies. This presents two related problems for historical dose reconstruction, both of which are addressed here. It is necessary to derive conversion factor values for (1) continuous distributions of energy typical of diagnostic medical x-rays (bremsstrahlung radiation), and (2) energies of particular radioisotopes used in medical procedures, neither of which are presented in published tables. For derivation of DCCs for bremsstrahlung radiation, combinations of x-ray tube potentials and filtrations were derived for different time periods based on a review of relevant literature. Three peak tube potentials (70 kV, 80 kV, and 90 kV) with four different amounts of beam filtration were determined to be applicable for historic dose reconstruction. The probabilities of these machine settings were assigned to each of the four time periods (earlier than 1949, 1949-1954, 1955-1968, and after 1968). Continuous functions were fit to each set of discrete values of the ICRP/ICRU mono-energetic DCCs and the functions integrated over the air-kerma weighted photon fluence of the 12 defined x-ray spectra. The air kerma-weighted DCCs in this work were developed specifically for an irradiation geometry of anterior to posterior (AP) and for the following tissues: thyroid, breast, ovary, lens of eye, lung, colon, testes, heart, skin (anterior side only), red bone marrow (RBM), and brain. In addition, a series of functional relationships to predict DT Ka-1 values for RBM dependent on body mass index [BMI (kg m-2) ≡ weight per height] and average photon energy were derived from a published analysis. Factors to account for attenuation of radiation by protective lead aprons were also developed. Because lead protective aprons often worn by radiology personnel not only reduce the intensity of x-ray exposure but also appreciably harden the transmitted fluence of bremsstrahlung x-rays, DCCs were separately calculated for organs possibly protected by lead aprons by considering three cases: no apron, 0.25 mm Pb apron, and 0.5 mm Pb apron. For estimation of organ doses from conducting procedures with radioisotopes, continuous functions of the reported mono-energetic values were developed, and DCCs were derived by estimation of the function at relevant energies. By considering the temporal changes in primary exposure-related parameters (e.g., energy distribution), the derived DCCs and transmission factors presented here allow for more realistic historical dose reconstructions for medical personnel when monitoring badge readings are the primary data on which estimation of an individual's organ doses are based.

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.