Inclusion of thin target and source regions in alimentary and respiratory tract systems of mesh-type ICRP adult reference phantoms

Enhanced Internal Dosimetry for Alimentary Tract Organs in Nuclear Medicine based on the ICRP Mesh-Type Reference Phantoms

This study introduces a refined approach for more accurately estimating radiation doses to alimentary tract organs in nuclear medicine, by utilizing the ICRP pediatric and adult mesh-type reference computational phantoms (MRCPs) that improved the anatomical representation of these organs. Our initial step involved compiling a comprehensive dataset of electron Specific Absorbed Fractions (SAFs) for all source-target pairs of alimentary tract organs in both adult and pediatric phantoms, calculating SAFs for all cases in the present study only except those computed in the previous study for certain pediatric phantom cases. Subsequently, we determined S values for 1,252 radionuclides, facilitating dosimetry applications. The consistency of target and source masses for alimentary tract organs in the MRCPs with the reference values in ICRP Publication 89 led to noticeable differences in SAF, S values, and consequently, absorbed dose coefficients when compared to the stylized models in ICRP Publication 100. Notably, the S value ratios (MRCP/stylized) for selected radionuclides-11C, 18F, 68Ga, and 131I-ranged from 0.41 to 7.60. Particularly for therapeutic 131I-iodide in thyroid cancer, the use of MRCPs resulted in up to 1.49 times higher absorbed dose coefficients for the colon than those derived from stylized models, while the stomach dose coefficients decreased by a factor of 0.72. The application of our findings promises enhanced, more realistic dosimetry for alimentary tract organs, especially beneficial for radiopharmaceuticals likely to accumulate within these organs.

Radioiodine internal dose coefficients specific for Koreans

This study developed internal dose coefficients for radioiodine, tailored to the Korean population, by incorporating the Korean biokinetic model along with the Korean S values. The observed differences in dose coefficients for Koreans compared to the International Commission on Radiological Protection (ICRP) reference values noticeably varied depending on physical half-lives of iodine isotopes. For longer-lived isotopes such as I-125 and I-129, significant differences in thyroid dose coefficients were observed, with ratios (Korean/ICRP) from 0.30 to 0.55, indicating that actual doses for Koreans can be considerably lower than those evaluated based on the ICRP data. However, for short-lived iodine isotopes, such as I-131, the thyroid dose coefficients were comparable to the ICRP reference values (ratio=0.95-0.98). These comparable dose coefficients resulted from the lower thyroidal iodine uptake in the Korean model being almost entirely offset by the higher thyroid self-absorption S values in the Korean phantoms. Additionally, this study delves into the substantial differences in absorbed dose coefficients for non-thyroidal regions and effective dose coefficients, which arose not only from physiological/anatomical variability but also technical differences in phantom design. The use of Korean-specific dose coefficients is advisable particularly in scenarios predicting elevated doses, yielding a more precise and clinically relevant dose assessment.

Development of Respiratory Tract Organs for ICRP Pediatric Mesh-type Reference Computational Phantoms

ABSTRACT

As part of the activities of the International Commission on Radiological Protection (ICRP) Task Group 103, the present study developed a new set of respiratory tract organs consisting of the extrathoracic, bronchial, bronchiolar, and alveolar-interstitial regions for newborn, 1-, 5-, 10-, and 15-y-old males and females for use in pediatric mesh-type reference computational phantoms. The developed respiratory tract organs, while preserving the original topologies of those of the pediatric voxel-type reference computational phantoms of ICRP Publication 143, have improved anatomy and detailed structure and also include μm-thick target and source regions prescribed in ICRP Publication 66. The dosimetric impact of the developed respiratory tract organs was investigated by calculating the specific absorbed fraction for internal electron exposures, which were then compared with the ICRP Task Group 96 values. The results showed that except for the alveolar-interstitial region as a source region, the pediatric mesh phantoms showed larger specific absorbed fractions than the Task Group 96 values. The maximum difference was a factor of ~3.5 for the extrathoracic-2 basal cell and surface as target and source regions, respectively. These results reflect the differences in the target masses and geometry caused by the anatomical enhancement of the pediatric mesh phantoms. For the alveolar-interstitial region as a source region, the pediatric mesh phantoms showed larger values for low energy ranges and lower values with increasing energies, owing to the differences in the size and shape of the alveolar-interstitial region.

Stylized versus voxel phantoms: quantification of internal organ chord length distances

Dosimetric calculations, whether for radiation protection or nuclear medicine applications, are greatly influenced by the use of computational models of humans, called anthropomorphic phantoms. As anatomical models of phantoms have evolved and expanded, thus has the need for quantifying differences among each of these representations that yield variations in organ dose coefficients, whether from external radiation sources or internal emitters. This work represents an extension of previous efforts to quantify the differences in organ positioning within the body between a stylized and voxel phantom series. Where prior work focused on the organ depth distribution vis-à-vis the surface of the phantom models, the work described here quantifies the intra-organ and inter-organ distributions through calculation of the mean chord lengths. The revised Oak Ridge National Laboratory stylized phantom series and the University of Florida/National Cancer Institute voxel phantom series including a newborn, 1-, 5-, 10- and 15 year old, and adult phantoms were compared. Organ distances in the stylized phantoms were computed using a ray-tracing technique available through Monte Carlo radiation transport simulations in MCNP6. Organ distances in the voxel phantom were found using phantom matrix manipulation. Quantification of differences in organ chord lengths between the phantom series displayed that the organs of the stylized phantom series are typically situated farther away from one another than within the voxel phantom series. The impact of this work was to characterize the intra-organ and inter-organ distributions to explain the variations in updated internal dose coefficient quantities (i.e. specific absorbed fractions) while providing relevant data defining the spatial and volumetric organ distributions in the phantoms for use in subsequent internal dosimetric computations, with prospective relevance to patient-specific individualized dosimetry, as well as informing machine learning definition of organs using these reference models.

Development of alimentary tract organs for ICRP pediatric mesh-type reference computational phantoms

In line with the activities of Task Group 103 under the International Commission on Radiological Protection (ICRP), the present study was conducted to develop a new set of alimentary tract organs consisting of the oral cavity, oesophagus, stomach, small intestine, and colon for the newborn, 1 year-old, 5 year-old, 10 year-old, and 15 year-old males and females for use in the pediatric mesh-type reference computational phantoms (MRCPs). The developed alimentary tract organs of the pediatric MRCPs, while nearly preserving the original topology and shape of those of the pediatric voxel-type reference computational phantoms (VRCPs) of ICRPPublication 143, present considerable anatomical improvement and include all micrometre-scale target and source regions as prescribed in ICRPPublication 100. To investigate the dosimetric impact of the developed alimentary tract organs, organ doses and specific absorbed fractions were computed for certain external exposures to photons and electrons and internal exposures to electrons, respectively, which were then compared with the values computed using the current ICRP models (i.e. pediatric VRCPs and ICRP-100 stylised models). The results showed that for external exposures to penetrating radiations (i.e. photons >0.04 MeV), there was generally good agreement between the compared values, within a 10% difference, except for the oral mucosa. For external exposures to weakly penetrating radiations (i.e. low-energy photons and electrons), there were significant differences, up to a factor of ∼8300, owing to the geometric difference caused by the anatomical enhancement in the MRCPs. For internal exposures of electrons, there were significant differences, the maximum of which reached a factor of ∼73 000. This was attributed not only to the geometric difference but also to the target mass difference caused by the different luminal content mass and organ shape.

An update on computational anthropomorphic anatomical models

The prevalent availability of high-performance computing coupled with validated computerized simulation platforms as open-source packages have motivated progress in the development of realistic anthropomorphic computational models of the human anatomy. The main application of these advanced tools focused on imaging physics and computational internal/external radiation dosimetry research. This paper provides an updated review of state-of-the-art developments and recent advances in the design of sophisticated computational models of the human anatomy with a particular focus on their use in radiation dosimetry calculations. The consolidation of flexible and realistic computational models with biological data and accurate radiation transport modeling tools enables the capability to produce dosimetric data reflecting actual setup in clinical setting. These simulation methodologies and results are helpful resources for the medical physics and medical imaging communities and are expected to impact the fields of medical imaging and dosimetry calculations profoundly.

Development of paediatric mesh-type reference computational phantom series of International Commission on Radiological Protection

Very recently, Task Group 103 of the International Commission on Radiological Protection (ICRP) completed the development of the paediatric mesh-type reference computational phantoms (MRCPs) comprising ten phantoms (newborn, one year-old, five year-old, ten year-old, and fifteen year-old males and females). The paediatric MRCPs address the limitations of ICRPPublication 143's paediatric reference computational phantoms, which are in voxel format, stemming from the nature of the voxel geometry and the limited voxel resolutions. The paediatric MRCPs were constructed by converting the voxel-type reference phantoms to a high-quality mesh format with substantial enhancements in the detailed anatomy of the small and complex organs and tissues (e.g. bones, lymphatic nodes, and extra-thoracic region). Besides, the paediatric MRCPs were developed in consideration of the intra-organ blood contents and by modelling the micron-thick target and source regions of the skin, lens, urinary bladder, alimentary tract organs, and respiratory tract organs prescribed by the ICRP. For external idealised exposures, the paediatric MRCPs provide very similar effective dose coefficients (DC_Es) to those from the ICRP-143 phantoms but significantly different values for weakly penetrating radiations (e.g. the difference of ∼20 000 times for 10 keV electron beams). This paper introduces the developed paediatric MRCPs with a brief explanation of the construction process. Then, it discusses their computational performance in Geant4, PHITS, and MCNP6 in terms of memory usage and computation speed and their impact on dose calculations by comparing their calculated values of DC_Es for external exposures with those of the voxel-type reference phantoms.

Development of detailed pediatric eye models for lens dose calculations

The International Commission on Radiological Protection (ICRP) recently reduced the dose limit for the eye lens for occupational exposure from 150 mSv yr^-1to 20 mSv yr^-1, as averaged over defined periods of five years, with no annual dose in a single year exceeding 50 mSv, emphasizing the importance of the accurate estimation of lens dose. In the present study, for more accurate lens dosimetry, detailed eye models were developed for children and adolescents (newborns and 1, 5, 10, and 15 year olds), which were then incorporated into the pediatric mesh-type reference computational phantoms (MRCPs) and used to calculate lens dose coefficients (DCs) for photon and electron exposures. Finally, the calculated values were compared with those calculated with the adult MRCPs in order to determine the age dependence of the lens DCs. For photon exposures, the lens DCs of the pediatric MRCPs showed some sizable differences from those of the adult MRCPs at very low energies (10 and 15 keV), but the differences were all less than 35%, except for the posterior-anterior irradiation geometry, for which the lens dose is not of primary concern. For electron exposures, much larger differences were found. For the anterior-posterior (AP) and isotropic irradiation geometries, the largest differences between the lens DCs of the pediatric and adult phantoms were found in the energy range of 0.6-1 MeV, where the newborn lens DCs were larger by up to a factor of ∼5 than the adult. The lens DCs of the present study, which were calculated for the radiosensitive region of the lens, also were compared with those for the entire lens in the AP irradiation geometry. Our results showed that the DCs of the entire lens were similar to those of the radiosensitive region for 0.02-2 MeV photons and >2 MeV electrons, but that for the other energy ranges, significant differences were noticeable, i.e. 10%-40% for photons and up to a factor of ∼5 for electrons.

Development of skeletal systems for ICRP pediatric mesh-type reference computational phantoms

In 2016, the International Commission on Radiological Protection (ICRP) launched Task Group 103 (TG 103) for the explicit purpose of developing a new generation of adult and pediatric reference computational phantoms, named 'mesh-type reference computational phantoms (MRCPs)', that can overcome the limitations of voxel-type reference computational phantoms (VRCPs) of ICRPPublications 110and143due to their finite voxel resolutions and the nature of voxel geometry. After completing the development of the adult MRCPs, TG 103 has started the development of pediatric MRCPs comprising 10 phantoms (male and female versions of the reference newborn, 1-year-old, 5-year-old, 10-year-old, and 15-year-old). As part of the TG 103 project, within the present study, the skeletal systems, one of the most important and complex organ systems of the body, were developed for each phantom age and sex. The developed skeletal systems, while closely preserving the original bone topology of the pediatric VRCPs, present substantial improvements in the anatomy of complex and/or small bones. In order to investigate the dosimetric impact of the developed skeletons, the average absorbed doses and the specific absorbed fractions for radiosensitive skeletal tissues (i.e. active marrow and bone endosteum) were computed for some selected external and internal exposure cases, which were then compared with those calculated with the skeletons of pediatric VRCPs. The comparison result showed that the dose values of the pediatric MRCPs were generally similar to those of the pediatric VRCPs for highly penetrating radiations (e.g. photons >200 keV); however, for weakly penetrating radiations (e.g. photons ⩽200 keV and electrons), significant differences up to a factor of 140 were observed.

POLY2TET: a computer program for conversion of computational human phantoms from polygonal mesh to tetrahedral mesh

As a geometrical format for computational human phantoms, tetrahedral mesh (TM) is known to have significant advantages over polygonal mesh (PM), including higher compatibility with Monte Carlo radiation transport codes, higher computation speed, and the capability of modeling heterogeneous density variation in an organ of the phantom. In the present study, a computer program named POLY2TET was developed to convert the format of computational human phantoms from PM to TM and generate a sample source code or input file, as applicable, for the converted phantom to be used in some general-purpose Monte Carlo radiation transport codes (i.e. Geant4, PHITS, and MCNP6). The developed program was then tested using four existing high-fidelity PM phantoms. The computation speed, memory requirement, and initialisation time of the generated TM phantoms were also measured and compared with those of the original PM phantoms in Geant4. From the results of our test, it was concluded that the developed program successfully converts PM phantoms into the TM format. The organ doses calculated using the generated TM phantom for the three Monte Carlo codes all produced essentially identical dose values to those for the original PM phantoms in Geant4. The comparison of computation speed showed that compared to the original PM phantoms in Geant4, the TM phantoms in the three Monte Carlo codes were much faster in transporting the particles considered in the present study, i.e. by up to ∼2600 times for electron beams simulated in PHITS. The comparison of the memory requirement showed that the TM phantoms required more memory than the original PM phantoms, but, except for MCNP6, the memory required for the TM phantoms was still less than 12 GB, which typically is available in personal computers these days. For MCNP6, the required memory was much higher, i.e. 60-70 GB.

Neuroman: Voxel Phantoms from Surface Models of 300 Head Structures Including 12 Pairs of Cranial Nerves

For a precise simulation of electromagnetic radiation effects, voxel phantoms require detailed structures to approximate humans. The phantoms currently used still do not have sophisticated structures. This paper presents voxel and surface models of 300 head structures with cranial nerves and reports on a technique for voxel reconstruction of the cranial nerves having very thin and small structures. In real-color sectioned images of the head (voxel size: 0.1 mm), 300 structures were segmented using Photoshop. A surface reconstruction was performed automatically on Mimics. Voxel conversion was run on Voxel Studio. The abnormal shapes of the voxel models were found and classified into three types: thin cord, thin layers, and thin parts in the structures. The abnormal voxel models were amended using extended, filled, and manual voxelization methods devised for this study. Surface models in STL format and as PDF files of the 300 head structures were produced. The STL format has good scalability, so it can be used in most three-dimensional surface model software. The PDF file is very user friendly for students and researchers who want to learn the head anatomy. Voxel models of 300 head structures were produced (TXT format), and their voxel quantity and weight were measured. A voxel model is difficult to handle, and the surface model cannot use the radiation simulation. Consequently, the best method for making precise phantoms is one in which the flaws of the voxel and surface models complement each other, as in the present study.

Body-size-dependent phantom library constructed from ICRP mesh-type reference computational phantoms

Recently, ICRP Task Group 103 developed new mesh-type reference computational phantoms (MRCPs) for the adult male and female by converting the current voxel-type reference computational phantoms (VRCPs) of ICRP Publication 110 into a high-quality/fidelity mesh format. Utilizing the great deformability/flexibility of the MRCPs compared with the VRCPs, in the present study, we established a body-size-dependent phantom library by modifying the MRCPs. The established library includes 108 adult male and 104 adult female phantoms in different standing heights and body weights, covering most body sizes representative of Caucasian and Asian populations. Ten secondary anthropometric parameters with respect to standing height and body weight were derived from various anthropometric databases and applied in the construction of the phantom library. An in-house program for automatic phantom adjustment was developed and applied for practical construction of such a large number of phantoms in the library with minimized human intervention. Organ/tissue doses calculated with three male phantoms in different standing heights (165, 175, and 190 cm) with a fixed body weight of 80 kg for external exposures to broad parallel photon beams from 0.01 to 10⁴ MeV were compared, observing there are significant dose differences particularly for the photon energies <0.1 MeV in which the organ/tissue doses tended to increase with increasing standing height. In addition, the organ/tissue doses of three female phantoms in different body weights (45, 65, and 140 kg) with a fixed standing height of 165 cm were compared, showing a significant decreasing tendency with increasing body weight for the photon energies <10 MeV. For the higher energies, the opposite trend, interestingly, was observed; that is, the organ/tissue doses tended to increase with increasing body weight. The results, despite the limited number of exposure cases, suggest that the use of the body-size-dependent phantom library can improve the accuracy of individual dose estimates for many retrospective dosimetry studies by taking the body size of individuals into account.

Computation Speeds and Memory Requirements of Mesh-Type ICRP Reference Computational Phantoms in Geant4, MCNP6, and PHITS

Recently, Task Group 103 of the International Commission on Radiological Protection completed the development of new adult male and female mesh-type reference computational phantoms, which are planned for use in future International Commission on Radiological Protection dose coefficient calculations. In the present study, the performance of major Monte Carlo particle transport codes, i.e., Geant4, MCNP6, and PHITS, were investigated for the mesh-type reference computational phantoms by performing transport simulations of photons, electrons, neutrons, and helium ions for some external and internal exposures, and simultaneously measuring the memory usage, initialization time, and computation speed of the adult male mesh-type reference computational phantom in the codes. The measured results were then compared with the values measured with the current adult male voxel-type reference computational phantom in International Commission on Radiological Protection Publication 110 as well as five voxel phantoms produced from the adult male mesh-type reference computational phantom with different voxel resolutions, i.e., 0.1 × 0.1 × 0.1 mm, 0.6 × 0.6 × 0.6 mm, 1 × 1 × 1 mm, 2 × 2 × 2 mm, and 4 × 4 × 4 mm. From the results, it was found that in all of the codes, the memory usage of the mesh-type reference computational phantom is greater than that of the voxel-type reference computational phantom and the lowest resolution voxelized phantom, but it is sufficiently lower than the maximum memory, 64 GB, that can be installed in a personal computer. The required initialization time of the mesh-type reference computational phantom and of the voxel-type reference computational phantom and voxelized phantoms in resolutions lower than 0.6 × 0.6 × 0.6 mm was less than a few minutes in all of the codes. As for the computation speed among the codes, MCNP6 showed the worst performance for the mesh-type reference computational phantom, which was slower than that for the voxel-type reference computational phantom by up to ~50 times and slower than that for all of the voxelized phantoms by up to ~40 times. By contrast, PHITS showed the best performance for the mesh-type reference computational phantom, which was faster than that for the voxel-type reference computational phantom by up to ~3 times and faster than that for all of the voxelized phantoms by up to ~20 times. This high performance of PHITS is indeed encouraging considering that it is used nowadays by the International Commission on Radiological Protection for most dose coefficient calculations.

Mesh-type reference Korean phantoms (MRKPs) for adult male and female for use in radiation protection dosimetry

In the present study, to overcome the dosimetric limitations of the previous voxel-type reference Korean computational phantoms due to their limited voxel resolutions (i.e. on the order of millimeters) and the nature of voxel geometry, a pair of new reference Korean phantoms, called mesh-type reference Korean phantoms (MRKPs), were developed for the adult male and female in a high-quality/fidelity mesh format. The developed phantoms include all target and source regions required for effective dose calculation, even micrometer-scale target and source regions of the respiratory and alimentary tract organs, skin, urinary bladder, and eye lens. The developed phantoms, which are in either the polygon-mesh (PM) format or the tetrahedral-mesh (TM) format as necessary, can be directly used in several general-purpose Monte Carlo codes (e.g. Geant4, MCNP6, and PHITS) without voxelization. In order to understand the dosimetric impact of the new phantoms, the dose coefficients (=fluence-to-effective dose conversion coefficients) were calculated for photons and electrons with energies ranging from 10 keV to 10 GeV for the anterior-posterior (AP) irradiation geometry and compared with those of the previous voxel-type reference Korean phantoms. The results demonstrate that the effective dose coefficients of the MRKPs were generally similar to those of the previous voxel-type reference phantoms for photons; however, for electrons, significant differences were observed at energies lower than 1 MeV that were mainly due to the explicit definition of the 50 µm-thick radiosensitive target layer in the skin of the new mesh phantoms.

New mesh-type phantoms and their dosimetric applications, including emergencies

Committee 2 of the International Commission on Radiological Protection (ICRP) has constructed mesh-type adult reference computational phantoms by converting the voxel-type ICRP Publication 110 adult reference computational phantoms to a high-quality mesh format, and adding those tissues that were below the image resolution of the voxel phantoms and therefore not included in the Publication 110 phantoms. The new mesh phantoms include all the necessary source and target tissues for effective dose calculations, including the 8-40-µm-thick target layers of the alimentary and respiratory tract organs, thereby obviating the need for supplemental organ-specific stylised models (e.g. respiratory airways, alimentary tract organ walls and stem cell layers, lens of the eye, and skin basal layer). To see the impact of the new mesh-type reference phantoms, dose coefficients for some selected external and internal exposures were calculated and compared with the current reference values in ICRP Publications 116 and 133, which were calculated by employing the Publication 110 phantoms and the supplemental stylised models. The new mesh phantoms were also used to calculate dose coefficients for industrial radiography sources near the body, which can be used to estimate the organ doses of the worker who is accidentally exposed by an industrial radiography source; in these calculations, the mesh phantoms were deformed to reflect the size of the worker, and also to evaluate the effect of posture on dose coefficients.

Gamma and Beta Absorbed Dose Conversion Coefficients in the Range from 10 keV to 10 MeV for Accidental Exposures From Point Sources Placed in Clothing in Proximity to the Body

Retrospective dose assessment following acute radiation exposures during radiological incidents can be difficult and inaccurate due to the large uncertainties associated with dose estimation. However, rapid and accurate dose assessment is critical following an incident so that appropriate treatment can be provided to the patient as early as possible. Incident dose assessment relies heavily on biological dosimetry with corresponding large uncertainties for inhomogeneous exposures, resulting from the estimates of whole-body doses, while the assessment of absorbed doses to individual tissues might actually be more appropriate for acute radiation exposures. Incident exposure scenarios for orphan sources placed in a breast or back pants pocket were modeled using the International Commission on Radiological Protection computational reference male and female and the Monte Carlo N-particle code MCNP6 to compute absorbed dose conversion coefficients for organs of interest for monoenergetic photon and beta sources. The absorbed dose conversion coefficients are intended for use in conjunction with source information to rapidly estimate absorbed doses to organs of interest from radiological sources in one of the two pocket geometries. Absorbed dose conversion coefficients also have been calculated specifically for Co, Cs, and Ir. Those absorbed dose conversion coefficients were applied to data from a radiological incident in Yanango, Peru, for comparison with published dose assessments; the results agree within 20%. The conversion coefficients are expected to provide an accurate tool for assessing doses for the modeled geometries, provided uncertainties due to the exact source-body geometry and exposure time are considered.

Multi-threading performance of Geant4, MCNP6, and PHITS Monte Carlo codes for tetrahedral-mesh geometry

In this study, the multi-threading performance of the Geant4, MCNP6, and PHITS codes was evaluated as a function of the number of threads (N) and the complexity of the tetrahedral-mesh phantom. For this, three tetrahedral-mesh phantoms of varying complexity (simple, moderately complex, and highly complex) were prepared and implemented in the three different Monte Carlo codes, in photon and neutron transport simulations. Subsequently, for each case, the initialization time, calculation time, and memory usage were measured as a function of the number of threads used in the simulation. It was found that for all codes, the initialization time significantly increased with the complexity of the phantom, but not with the number of threads. Geant4 exhibited much longer initialization time than the other codes, especially for the complex phantom (MRCP). The improvement of computation speed due to the use of a multi-threaded code was calculated as the speed-up factor, the ratio of the computation speed on a multi-threaded code to the computation speed on a single-threaded code. Geant4 showed the best multi-threading performance among the codes considered in this study, with the speed-up factor almost linearly increasing with the number of threads, reaching ~30 when N  =  40. PHITS and MCNP6 showed a much smaller increase of the speed-up factor with the number of threads. For PHITS, the speed-up factors were low when N  =  40. For MCNP6, the increase of the speed-up factors was better, but they were still less than ~10 when N  =  40. As for memory usage, Geant4 was found to use more memory than the other codes. In addition, compared to that of the other codes, the memory usage of Geant4 more rapidly increased with the number of threads, reaching as high as ~74 GB when N  =  40 for the complex phantom (MRCP). It is notable that compared to that of the other codes, the memory usage of PHITS was much lower, regardless of both the complexity of the phantom and the number of threads, hardly increasing with the number of threads for the MRCP.

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.

Individualized adjustments to reference phantom internal organ dosimetry-scaling factors given knowledge of patient internal anatomy

Various computational tools are currently available that facilitate patient organ dosimetry in diagnostic nuclear medicine, yet they are typically restricted to reporting organ doses to ICRP-defined reference phantoms. The present study, while remaining computational phantom based, provides straightforward tools to adjust reference phantom organ dose for both internal photon and electron sources. A wide variety of monoenergetic specific absorbed fractions were computed using radiation transport simulations for tissue spheres of varying size and separation distance. Scaling methods were then constructed for both photon and electron self-dose and cross-dose, with data validation provided from patient-specific voxel phantom simulations, as well as via comparison to the scaling methodology given in MIRD Pamphlet No. 11. Photon and electron self-dose was found to be dependent on both radiation energy and sphere size. Photon cross-dose was found to be mostly independent of sphere size. Electron cross-dose was found to be dependent on sphere size when the spheres were in close proximity, owing to differences in electron range. The validation studies showed that this dataset was more effective than the MIRD 11 method at predicting patient-specific photon doses for at both high and low energies, but gave similar results at photon energies between 100 keV and 1 MeV. The MIRD 11 method for electron self-dose scaling was accurate for lower energies but began to break down at higher energies. The photon cross-dose scaling methodology developed in this study showed gains in accuracy of up to 9% for actual patient studies, and the electron cross-dose scaling methodology showed gains in accuracy up to 9% as well when only the bremsstrahlung component of the cross-dose was scaled. These dose scaling methods are readily available for incorporation into internal dosimetry software for diagnostic phantom-based organ dosimetry.