Monte Carlo calculation of specific absorbed fractions: variance reduction techniques

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.

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal photon sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction-fraction of particle energy emitted in the source tissue that is deposited in the target tissue-and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series-male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old-for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source-target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log-log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. The results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

Simple variance reduction in Monte Carlo calculations of specific absorbed fractions: Russian roulette and splitting at the source organ

PURPOSE

To investigate the capabilities of several variance reduction techniques in the calculation of specific absorbed fractions in cases where the source and the target organs are far away and/or the target organs have a small volume.

METHODS

The specific absorbed fractions have been calculated by using the Monte Carlo code PENELOPE and by assuming the thyroid gland as the source organ and the testicles, the urinary bladder, the uterus, and the ovaries as the target ones. A mathematical anthropomorphic phantom, similar to the MIRD-type phantoms, has been considered. Photons with initial energies of 50, 100 and 500 keV were emitted isotropically from the volume of the source organ. Simulations have been carried out by implementing the variance reduction techniques of splitting and Russian roulette at the source organ only and the interaction forcing at the target organs. The influence of the implementation details of those techniques have been investigated and optimal parameters have been determined. All simulations were run with a CPU time of 1.5 · 10⁵ s.

RESULTS

Specific absorbed fractions with relative uncertainties well below 10% have been obtained in most cases, agreeing with those used as reference. The best value for the factor defining the application of the Russian roulette technique was r = 0.3. The best value for the splitting number was between s = 3 and s = 10, depending on the specific energies and target organs.

CONCLUSIONS

The proposed strategy provides an effective method for computing specific absorbed fractions for the most unfavorable situations, with a computing effort that is considerably reduced with respect to other methodologies.

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.