MPPD: A User-Friendly Posture Deformation Program for Mesh-Type Computational Phantoms.

ABSTRACT: Recently, the International Commission on Radiological Protection (ICRP) released adult Mesh-type Reference Computational Phantoms (MRCPs), which have great advantage in high deformability. Previous studies have exploited their high deformability to investigate the dosimetric influence of varying statures and postures, demonstrating significant variations in radiation doses. However, the previous studies are constrained by their inability to consider both stature and posture concurrently and by the limited range of postures analyzed. In the present study, a computer program named MPPD (Mesh-type Phantom Posture Deformer) was developed, a user-friendly graphical user interface that enables users to adjust the posture of adult MRCPs and corresponding library phantoms. The MPPD program was applied to deform five adult male phantoms of different statures into sitting and kneeling postures, showcasing its rapid computational speed and minimal RAM usage. The effectiveness of the MPPD program for dose calculation was also investigated by computing the detriment-weighted doses for MPPD-deformed adult male MRCPs, which showed good agreement with dose values for existing posture-deformed phantoms of the previous study. Furthermore, as an application of the MPPD program, the combined dosimetric impact of stature and posture was investigated, which is the inaugural effort to estimate doses by considering these factors concurrently. The result showed that the impact of stature and posture on radiation doses could largely vary depending on the radiation source, highlighting the importance of simultaneous consideration of stature and posture for accurate dose estimation.

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.

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.

Influence of Body Posture on Internal Organ Dosimetry: Radiocesium Exposure Modeling Using Novel Posture-dependent Mesh Computational Phantoms.

ABSTRACT: Current practice in reference internal dosimetry assumes a fixed upright standing posture is maintained throughout the dose-integration period. Recently, the mesh-type ICRP adult reference computational phantoms were transformed into different body postures (e.g., sitting, squatting) for use in occupational dose reconstruction applications. Here, for the first time, we apply this phantom series to the study of organ dose estimates following radionuclide intake. We consider the specific cases of 137 Cs and 134 Cs ingestion (accidental/occupational intake) with attention to variability in absorbed dose as a function of posture. The ICRP Publication 137 systemic biokinetic model for soluble cesium ingestion was used to compute organ-level time-integrated activity coefficients for reference adults, over a 50-y dose-integration period, for 134 Cs and 137 Cs (and its radioactive progeny 137m Ba). Mean posture time-allocations (h d -1 for standing, sitting, and lying) were taken from published survey data. In accord with modern dosimetry formalisms (e.g., MIRD, ICRP), a posture weighting factor was introduced that accounts for the fraction of time spent within each independent posture. Absorbed dose coefficients were computed using PHITS Monte Carlo simulations. ICRP 103 tissue weighting factors were applied along with the posture weighting factors to obtain committed effective dose per unit intake (Sv Bq -1 ). For 137 Cs ingestion, most organ absorbed dose coefficients were negligibly to marginally higher (< ~3%) for sitting or crouched (lying fetal/semi-fetal) postures maintained over the dose commitment period, relative to the upright standing posture. The committed effective dose coefficients were 1.3 × 10 -8 Sv Bq -1 137 Cs for standing, sitting, or crouched postures; thus, the posture-weighted committed effective dose was not significantly different than the committed effective dose for a maintained upright standing posture. For 134 Cs ingestion, most organ absorbed dose coefficients for the sitting and crouched postures were significantly larger than the standing posture, but the differences were still considered minor (< ~8% for most organs). The committed effective dose coefficients were 1.2 × 10 -8 Sv Bq -1 134 Cs for the standing posture and 1.3 × 10 -8 Sv Bq -1 134 Cs for the sitting or crouched posture. The posture-weighted committed effective dose was 1.3 × 10 -8 Sv Bq -1 134 Cs. Body posture has minor influence on organ-level absorbed dose coefficients and committed effective dose for ingestion of soluble 137 Cs or 134 Cs.

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.

Incorporation of micro-CT-based detailed bone models into ICRP adult mesh-type reference computational phantoms.

Objective.The red bone marrow (RBM) and bone endosteum (BE), which are required for effective dose calculation, are macroscopically modeled in the reference phantoms of the international commission on radiological protection (ICRP) due to their microscopic and complex histology. In the present study, the detailed bone models were developed to simplify the dose calculation process for skeletal dosimetry.Approach.The detailed bone models were developed based on the bone models developed at the University of Florida. A new method was used to update the definition of BE region by storing the BE location indices using virtual sub-voxels. The detailed bone models were then installed in the spongiosa regions of the ICRP mesh-type reference computational phantoms (MRCPs) via the parallel geometry feature of the Geant4 code.Main results.Comparing the results between the detailed-bone-installed MRCPs and the original MRCPs with the absorbed dose to spongiosa and fluence-to-dose response function (DRF)-based methods, the DRF-based method showed much smaller but still significant differences. Compared with the values given in ICRPPublications116 and 133, the differences were very large (i.e. several orders of magnitudes), due mainly to the anatomical improvement of the skeletal system in the MRCPs; that is, spongiosa and medullary cavity are fully enclosed by cortical bone in the MRCPs but not in the ICRP-110 phantoms.Significance.The detailed bone models enable the direct calculation of the absorbed doses to the RBM and BE, simplifying the dose calculation process and potentially improving the consistency and accuracy of skeletal dosimetry.

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 (DCEs) 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 DCEs 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.

Patient Size-Dependent Dosimetry Methodology Applied to 18F-FDG Using New ICRP Mesh Phantoms.

Despite the known influence of anatomic variability on internal dosimetry, dosimetry for 18F-FDG and other diagnostic radiopharmaceuticals is routinely derived using reference phantoms, which embody population-averaged morphometry for a given age and sex. Moreover, phantom format affects dosimetry estimates to varying extent. Here, we applied newly developed mesh format reference phantoms and a patient-dependent phantom library to assess the impact of height, weight, and body contour variation on dosimetry of 18F-FDG. We compared the mesh reference phantom dosimetry estimates with corresponding estimates from common software to identify differences related to phantom format or software implementation. Our study serves as an example of how more precise patient size-dependent dosimetry methodology could be performed. Methods: Absorbed dose coefficients were computed for the adult mesh reference phantoms and derivative patient-dependent phantom series by Monte Carlo simulation using the PHITS radiation transport code within PARaDIM software. The dose coefficients were compared with reference absorbed dose coefficients obtained from ICRP Publication 128, or generated using software including OLINDA 2.1, OLINDA 1.1, and IDAC-dose 2.1. Results: Differences in dosimetry arising from anatomical variations were shown to be significant, with detriment-weighted dose coefficients for the percentile-specific phantoms varying by up to ±40% relative to the corresponding reference phantom effective dose coefficients, irrespective of phantom format. Similar variations were seen in the individual organ absorbed dose coefficients for the percentile-specific phantoms relative to the reference phantoms. The effective dose coefficient for the mesh reference adult was 0.017 mSv/MBq, which was 5% higher than estimated by a corresponding voxel phantom, and 10% lower than estimated by the stylized phantom format. Conclusion: We observed notable variability in 18F-FDG dosimetry across morphometrically different patients, supporting the use of patient-dependent phantoms for more accurate dosimetric estimations relative to standard reference dosimetry. These data may help in optimizing imaging protocols and research studies, in particular when longer-lived isotopes are employed.

Detailed tooth models for ICRP mesh-type reference computational phantoms.

For use in electron paramagnetic resonance dosimetry with tooth enamel, in the present study, very detailed mesh-type tooth models composed of 198 individual tooth models (i.e. newborn: 20; 1 year: 28; 5 years: 48; 10 years: 38; 15 years: 32; and adult: 32) were developed for each sex. The developed tooth models were then implanted in the International Commission on Radiological Protection pediatric and adult mesh-type reference computational phantoms and used to calculate tooth enamel doses, by Monte Carlo simulations with Geant4, for external photon exposures in several idealized irradiation geometries. The calculated dose values were then compared to investigate the dependency of the enamel dose on the age and sex of the phantom and the sites of the teeth. The results of the present study generally show that, if the photon energy is low (i.e. <0.1 MeV), the enamel dose is significantly affected by the age and sex of the phantom and also the sites of the teeth used for dose calculation; the differences are frequently greater than a few times or even orders of magnitude. However, with a few exceptions, the enamel dose was hardly affected by these parameters for energies between 0.1 and 3 MeV. For energies >3 MeV, moderate differences were observed (i.e., up to a factor of two), due to the existence of dose build-up in the head of the phantom for high-energy photons. The calculated dose values were also compared with those of the previous studies where voxel and mathematical models were used to calculate the enamel doses. The results again show significant differences at low energies, e.g., up to ∼3500 times at 0.015 MeV, which are mainly due to the differences in the level of tooth-modeling detailedness.

Dose conversion coefficients for neutron external exposures with five postures: walking, sitting, bending, kneeling, and squatting.

In a previous study, posture-dependent dose coefficients (DCs) for photon external exposures were calculated using the adult male and female mesh-type reference computational phantoms (MRCPs) of the International Commission on Radiological Protection (ICRP) that had been transformed into five non-standing postures (i.e. walking, sitting, bending, kneeling, and squatting). As an extension, the present study was conducted to establish another DC dataset for external exposures to neutrons by performing Monte Carlo radiation transport simulations with the adult male and female MRCPs in the five non-standing postures. The resulting dataset included the DCs for absorbed doses (i.e., organ/tissue-averaged absorbed doses) delivered to 29 individual organs/tissues, and for effective doses for neutron energies ranging from 10-9 to 104 MeV in six irradiation geometries: antero-posterior (AP), posteroanterior (PA), left-lateral (LLAT), right-lateral (RLAT), rotational (ROT), and isotropic (ISO) geometries. The comparison of DCs for the non-standing MRCPs with those of the standing MRCPs showed significant differences. In the lateral irradiation geometries, for example, the standing MRCPs overestimate the breast DCs of the squatting MRCPs by up to a factor of 4 due to the different arm positions but underestimate the gonad DCs by up to about 17 times due to the different leg positions. The impact of different postures on effective doses was generally less than that on organ doses but still significant; for example, the standing MRCPs overestimate the effective doses of the bending MRCPs only by 20% in the AP geometry at neutron energies less than 50 MeV, but underestimate those of the kneeling MRCPs by up to 40% in the lateral geometries at energies less than 0.1 MeV.

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.

Dosimetric considerations of 99mTc-MDP uptake within the epiphyseal plates of the long bones of pediatric patients.

Skeletal scintigraphy is most performed in pediatric patients using the radiopharmaceutical 99mTc labelled methylene diphosphonate (99mTc-MDP). Reference biokinetic models for 99mTc-MDP indicate 50% of the administered activity is uniformly localized to the interior bone surfaces (trabecular and cortical regions), yet imaging data clearly show some preferential uptake to the epiphyseal growth plates of the long bones. To explore the dosimetric consequences of these regional activity concentrations, we have modified mesh-type computational phantoms of the International Commission on Radiological Protection (ICRP) reference pediatric series to explicitly include geometric models of the epiphyseal growth plates (2 mm in thickness) within the left/right, distal/proximal ends of the humeri, radii, ulnae, femora, tibia, and fibulae. Bone mineral activity from the ICRP Publication 128 biokinetic model for 99mTc-MDP (ICRP 2015) was then partitioned to the growth plates at values of 0.5%, 4.4%, 8.3%, 12.2%, 16.1%, and 20%. Radiation transport simulations were performed to compute 99mTc S-values and organ dose coefficients to the soft tissues and to bone site-specific regions of spongiosa. As the percentage of bone activity assigned to the growth plates was increased (from 0.5% to 20%), absorbed doses to the soft tissue organs, active bone marrow, bone endosteum (BE), as well as the detriment-weighted dose, were shown to decrease from their nominal values (no substantial growth plate activity), while epiphyseal plate self-doses increased. In the 15 year old male phantom, moving from 0.5% to 20% relative bone activity within the epiphyseal plates resulted in a 15% reduction in active marrow (AM) and BE dose, a 10% reduction in mean soft tissue and detriment-weighted dose, and a 6.3-fold increase in epiphyseal plate self-dose. In the newborn female phantom, we observed a 18% decrease in AM and BE dose, a 10% decrease in mean soft tissue dose, a 15% decrease in detriment-weighted dose, and 12.8-fold increase in epiphyseal plate self-dose. Increases (to 3 mm) and decreases (to 1 mm) in the assumed growth plate thickness of our models were shown to impact only the growth plate self-dose. Future work in differential quantification of 99mTc-MDP activity-growth plates versus other bone surfaces-is required to provide clinically realistic data on activity partitioning as a function of patient age, and perhaps skeletal site. The phantom series presented here may be used to develop more optimized age-related guidance on 99mTc-MDP administered activities to children.

Body-size-dependent Iodine-131Svalues.

In a recent epidemiologic risk assessment on late health effects of patients treated with radioactive iodine (RAI), organ/tissue doses of the patients were estimated based on iodine-131Svalues derived from the reference computational phantoms of the International Commission on Radiological Protection (ICRP). However, the use of theSvalues based on the reference phantoms may lead to significant biases in the estimated doses of patients whose body sizes (height and weight) are significantly different from the reference body sizes. To fill this critical gap, we established a comprehensive dataset of body-size-dependent iodine-131Svalues (rT← thyroid) for 30 radiosensitive target organs/tissues by performing Monte Carlo dose calculations coupled with a total of 212 adult male and female computational phantoms in different heights and weights. We observed that theSvalues tend to decrease with increasing body height; for example, theSvalue (gonads ← thyroid) of the 160 cm male phantom is about 3 times higher than that of the 190 cm male phantom at the 70 kg weight. We also observed that theSvalues tend to decrease with increasing body weight for some organs/tissues; for example, theSvalue (skin ← thyroid) of the 45 kg female phantom is about two times higher than that of the 130 kg female phantom at the 160 cm height. For other organs/tissues, which are relatively far from the thyroid, in contrast, theSvalues tend to increase with increasing body weight; for example, theSvalue (bladder ← thyroid) of the 45 kg female phantom is about 2 times lower than that of the 130 kg female phantom. Overall, the majority of the body-size-dependentSvalues deviated to within 25% from those of the reference phantoms. We believe that the use of body-size-dependentSvalues in dose reconstructions should help quantify the dosimetric uncertainty in epidemiologic investigations of RAI-treated patients.

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.

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 104 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.

Dose coefficients of percentile-specific computational phantoms for photon external exposures.

The use of dose coefficients (DCs) based on the reference phantoms recommended by the International Commission on Radiological Protection (ICRP) with a fixed body size may produce errors to the estimated organ/tissue doses to be used, for example, for epidemiologic studies depending on the body size of cohort members. A set of percentile-specific computational phantoms that represent 10th, 50th, and 90th percentile standing heights and body masses in adult male and female Caucasian populations were recently developed by modifying the mesh-type ICRP reference computational phantoms (MRCPs). In the present study, these percentile-specific phantoms were used to calculate a comprehensive dataset of body-size-dependent DCs for photon external exposures by performing Monte Carlo dose calculations with the Geant4 code. The dataset includes the DCs of absorbed doses for 29 individual organs/tissues from 0.01 to 104 MeV photon energy, in the antero-posterior, postero-anterior, right lateral, left lateral, rotational, and isotropic geometries. The body-size-dependent DCs were compared with the DCs of the MRCPs in the reference body size, showing that the DCs of the MRCPs are generally similar to those of the 50th percentile standing height and body mass phantoms over the entire photon energy region except for low energies (≤ 0.03 MeV); the differences are mostly less than 10%. In contrast, there are significant differences in the DCs between the MRCPs and the 10th and 90th percentile standing height and body mass phantoms (i.e., H10M10 and H90M90). At energies of less than about 10 MeV, the MRCPs tended to under- and over-estimate the organ/tissue doses of the H10M10 and H90M90 phantoms, respectively. This tendency was revised at higher energies. The DCs of the percentile-specific phantoms were also compared with the previously published values of another phantom sets with similar body sizes, showing significant differences particularly at energies below about 0.1 MeV, which is mainly due to the different locations and depths of organs/tissues between the different phantom libraries. The DCs established in the present study should be useful to improve the dosimetric accuracy in the reconstructions of organ/tissue doses for individuals in risk assessment for epidemiologic investigations taking body sizes into account.

DOSE COEFFICIENTS OF MESH-TYPE ICRP REFERENCE COMPUTATIONAL PHANTOMS FOR EXTERNAL EXPOSURES OF NEUTRONS, PROTONS, AND HELIUM IONS.

To overcome inherent limitations of the Voxel-type Reference Computational Phantoms (VRCPs) due to the limited voxel resolutions and the nature of voxel geometry, the International Commission on Radiological Protection (ICRP) has developed the adult male and female Mesh-type Reference Computational Phantoms (MRCPs). We previously used the MRCPs to calculate a complete set of dose coefficients (DCs) for idealized external exposures of photons and electrons (Yeom et al. NET in press). In the present study, we extended the previous study to include additional radiation particles (neutrons, protons, and helium ions) into the DC library by conducing Monte Carlo radiation transport simulations with the Geant4 code. The MRPC-based DCs were compared with the existing reference DCs of ICRP Publication 116 which are based on the ICRP VRCPs to investigate impact of the new mesh-type reference phantoms on the DC values. We found that the MRCPs generally provide DCs of organ/tissue doses and effective doses similar to those from the VRCPs for penetrating radiations (uncharged particles), whereas significant DC differences were observed for weakly penetrating radiations (charged particles) mainly due to the improved representation of the detailed anatomical structures in the MRCPs over the VRCPs.

PARaDIM: A PHITS-Based Monte Carlo Tool for Internal Dosimetry with Tetrahedral Mesh Computational Phantoms.

Mesh-type and voxel-based computational phantoms comprise the current state of the art for internal dose assessment via Monte Carlo simulations but excel in different aspects, with mesh-type phantoms offering advantages over their voxel counterparts in terms of their flexibility and realistic representation of detailed patient- or subject-specific anatomy. We have developed PARaDIM (pronounced "paradigm": Particle and Heavy Ion Transport Code System-Based Application for Radionuclide Dosimetry in Meshes), a freeware application for implementing tetrahedral mesh-type phantoms in absorbed dose calculations. It considers all medically relevant radionuclides, including α, ß, γ, positron, and Auger/conversion electron emitters, and handles calculation of mean dose to individual regions, as well as 3-dimensional dose distributions for visualization and analysis in a variety of medical imaging software. This work describes the development of PARaDIM, documents the measures taken to test and validate its performance, and presents examples of its uses. Methods: Human, small-animal, and cell-level dose calculations were performed with PARaDIM and the results compared with those of widely accepted dosimetry programs and literature data. Several tetrahedral phantoms were developed or adapted using computer-aided modeling techniques for these comparisons. Results: For human dose calculations, agreement of PARaDIM with OLINDA 2.0 was good-within 10%-20% for most organs-despite geometric differences among the phantoms tested. Agreement with MIRDcell for cell-level S value calculations was within 5% in most cases. Conclusion: PARaDIM extends the use of Monte Carlo dose calculations to the broader community in nuclear medicine by providing a user-friendly graphical user interface for calculation setup and execution. PARaDIM leverages the enhanced anatomic realism provided by advanced computational reference phantoms or bespoke image-derived phantoms to enable improved assessments of radiation doses in a variety of radiopharmaceutical use cases, research, and preclinical development. PARaDIM can be downloaded freely at www.paradim-dose.org.

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.