Monte Carlo reference data sets for imaging research: Executive summary of the report of AAPM Research Committee Task Group 195.

The use of Monte Carlo simulations in diagnostic medical imaging research is widespread due to its flexibility and ability to estimate quantities that are challenging to measure empirically. However, any new Monte Carlo simulation code needs to be validated before it can be used reliably. The type and degree of validation required depends on the goals of the research project, but, typically, such validation involves either comparison of simulation results to physical measurements or to previously published results obtained with established Monte Carlo codes. The former is complicated due to nuances of experimental conditions and uncertainty, while the latter is challenging due to typical graphical presentation and lack of simulation details in previous publications. In addition, entering the field of Monte Carlo simulations in general involves a steep learning curve. It is not a simple task to learn how to program and interpret a Monte Carlo simulation, even when using one of the publicly available code packages. This Task Group report provides a common reference for benchmarking Monte Carlo simulations across a range of Monte Carlo codes and simulation scenarios. In the report, all simulation conditions are provided for six different Monte Carlo simulation cases that involve common x-ray based imaging research areas. The results obtained for the six cases using four publicly available Monte Carlo software packages are included in tabular form. In addition to a full description of all simulation conditions and results, a discussion and comparison of results among the Monte Carlo packages and the lessons learned during the compilation of these results are included. This abridged version of the report includes only an introductory description of the six cases and a brief example of the results of one of the cases. This work provides an investigator the necessary information to benchmark his/her Monte Carlo simulation software against the reference cases included here before performing his/her own novel research. In addition, an investigator entering the field of Monte Carlo simulations can use these descriptions and results as a self-teaching tool to ensure that he/she is able to perform a specific simulation correctly. Finally, educators can assign these cases as learning projects as part of course objectives or training programs.

Peak skin and eye lens radiation dose from brain perfusion CT based on Monte Carlo simulation.

OBJECTIVE: The purpose of our study was to accurately estimate the radiation dose to skin and the eye lens from clinical CT brain perfusion studies, investigate how well scanner output (expressed as volume CT dose index [CTDI(vol)]) matches these estimated doses, and investigate the efficacy of eye lens dose reduction techniques. MATERIALS AND METHODS: Peak skin dose and eye lens dose were estimated using Monte Carlo simulation methods on a voxelized patient model and 64-MDCT scanners from four major manufacturers. A range of clinical protocols was evaluated. CTDI(vol) for each scanner was obtained from the scanner console. Dose reduction to the eye lens was evaluated for various gantry tilt angles as well as scan locations. RESULTS: Peak skin dose and eye lens dose ranged from 81 mGy to 348 mGy, depending on the scanner and protocol used. Peak skin dose and eye lens dose were observed to be 66-79% and 59-63%, respectively, of the CTDI(vol) values reported by the scanners. The eye lens dose was significantly reduced when the eye lenses were not directly irradiated. CONCLUSION: CTDI(vol) should not be interpreted as patient dose; this study has shown it to overestimate dose to the skin or eye lens. These results may be used to provide more accurate estimates of actual dose to ensure that protocols are operated safely below thresholds. Tilting the gantry or moving the scanning region further away from the eyes are effective for reducing lens dose in clinical practice. These actions should be considered when they are consistent with the clinical task and patient anatomy.

The feasibility of patient size-corrected, scanner-independent organ dose estimates for abdominal CT exams.

PURPOSE: A recent work has demonstrated the feasibility of estimating the dose to individual organs from multidetector CT exams using patient-specific, scanner-independent CTDIvol-to-organ-dose conversion coefficients. However, the previous study only investigated organ dose to a single patient model from a full-body helical CT scan. The purpose of this work was to extend the validity of this dose estimation technique to patients of any size undergoing a common clinical exam. This was done by determining the influence of patient size on organ dose conversion coefficients generated for typical abdominal CT exams. METHODS: Monte Carlo simulations of abdominal exams were performed using models of 64-slice MDCT scanners from each of the four major manufacturers to obtain dose to radiosensitive organs for eight patient models of varying size, age, and gender. The scanner-specific organ doses were normalized by corresponding CTDIvol values and averaged across scanners to obtain scanner-independent CTDIvol-to-organ-dose conversion coefficients for each patient model. In order to obtain a metric for patient size, the outer perimeter of each patient was measured at the central slice of the abdominal scan region. Then, the relationship between CTDIvol-to-organ-dose conversion coefficients and patient perimeter was investigated for organs that were directly irradiated by the abdominal scan. These included organs that were either completely ("fully irradiated") or partly ("partially irradiated") contained within the abdominal exam region. Finally, dose to organs that were not at all contained within the scan region ("nonirradiated") were compared to the doses delivered to fully irradiated organs. RESULTS: CTDIvol-to-organ-dose conversion coefficients for fully irradiated abdominal organs had a strong exponential correlation with patient perimeter. Conversely, partially irradiated organs did not have a strong dependence on patient perimeter. In almost all cases, the doses delivered to nonirradiated organs were less than 5%, on average across patient models, of the mean dose of the fully irradiated organs. CONCLUSIONS: This work demonstrates the feasibility of calculating patient-specific, scanner-independent CTDIvol-to-organ-dose conversion coefficients for fully irradiated organs in patients undergoing typical abdominal CT exams. A method to calculate patient-specific, scanner-specific, and exam-specific organ dose estimates that requires only knowledge of the CTDIvol for the scan protocol and the patient's perimeter is thus possible. This method will have to be extended in future studies to include organs that are partially irradiated. Finally, it was shown that, in most cases, the doses to nonirradiated organs were small compared to the dose to fully irradiated organs.

The feasibility of a scanner-independent technique to estimate organ dose from MDCT scans: using CTDIvol to account for differences between scanners.

PURPOSE: Monte Carlo radiation transport techniques have made it possible to accurately estimate the radiation dose to radiosensitive organs in patient models from scans performed with modern multidetector row computed tomography (MDCT) scanners. However, there is considerable variation in organ doses across scanners, even when similar acquisition conditions are used. The purpose of this study was to investigate the feasibility of a technique to estimate organ doses that would be scanner independent. This was accomplished by assessing the ability of CTDIvol measurements to account for differences in MDCT scanners that lead to organ dose differences. METHODS: Monte Carlo simulations of 64-slice MDCT scanners from each of the four major manufacturers were performed. An adult female patient model from the GSF family of voxelized phantoms was used in which all ICRP Publication 103 radiosensitive organs were identified. A 120 kVp, full-body helical scan with a pitch of 1 was simulated for each scanner using similar scan protocols across scanners. From each simulated scan, the radiation dose to each organ was obtained on a per mA s basis (mGy/mA s). In addition, CTDIvol values were obtained from each scanner for the selected scan parameters. Then, to demonstrate the feasibility of generating organ dose estimates from scanner-independent coefficients, the simulated organ dose values resulting from each scanner were normalized by the CTDIvol value for those acquisition conditions. RESULTS: CTDIvol values across scanners showed considerable variation as the coefficient of variation (CoV) across scanners was 34.1%. The simulated patient scans also demonstrated considerable differences in organ dose values, which varied by up to a factor of approximately 2 between some of the scanners. The CoV across scanners for the simulated organ doses ranged from 26.7% (for the adrenals) to 37.7% (for the thyroid), with a mean CoV of 31.5% across all organs. However, when organ doses are normalized by CTDIvoI values, the differences across scanners become very small. For the CTDIvol, normalized dose values the CoVs across scanners for different organs ranged from a minimum of 2.4% (for skin tissue) to a maximum of 8.5% (for the adrenals) with a mean of 5.2%. CONCLUSIONS: This work has revealed that there is considerable variation among modern MDCT scanners in both CTDIvol and organ dose values. Because these variations are similar, CTDIvol can be used as a normalization factor with excellent results. This demonstrates the feasibility of establishing scanner-independent organ dose estimates by using CTDIvol to account for the differences between scanners.

A method to generate equivalent energy spectra and filtration models based on measurement for multidetector CT Monte Carlo dosimetry simulations.

The purpose of this study was to present a method for generating x-ray source models for performing Monte Carlo (MC) radiation dosimetry simulations of multidetector row CT (MDCT) scanners. These so-called "equivalent" source models consist of an energy spectrum and filtration description that are generated based wholly on the measured values and can be used in place of proprietary manufacturer's data for scanner-specific MDCT MC simulations. Required measurements include the half value layers (HVL1 and HVL2) and the bowtie profile (exposure values across the fan beam) for the MDCT scanner of interest. Using these measured values, a method was described (a) to numerically construct a spectrum with the calculated HVLs approximately equal to those measured (equivalent spectrum) and then (b) to determine a filtration scheme (equivalent filter) that attenuates the equivalent spectrum in a similar fashion as the actual filtration attenuates the actual x-ray beam, as measured by the bowtie profile measurements. Using this method, two types of equivalent source models were generated: One using a spectrum based on both HVL1 and HVL2 measurements and its corresponding filtration scheme and the second consisting of a spectrum based only on the measured HVL1 and its corresponding filtration scheme. Finally, a third type of source model was built based on the spectrum and filtration data provided by the scanner's manufacturer. MC simulations using each of these three source model types were evaluated by comparing the accuracy of multiple CT dose index (CTDI) simulations to measured CTDI values for 64-slice scanners from the four major MDCT manufacturers. Comprehensive evaluations were carried out for each scanner using each kVp and bowtie filter combination available. CTDI experiments were performed for both head (16 cm in diameter) and body (32 cm in diameter) CTDI phantoms using both central and peripheral measurement positions. Both equivalent source model types result in simulations with an average root mean square (RMS) error between the measured and simulated values of approximately 5% across all scanner and bowtie filter combinations, all kVps, both phantom sizes, and both measurement positions, while data provided from the manufacturers gave an average RMS error of approximately 12% pooled across all conditions. While there was no statistically significant difference between the two types of equivalent source models, both of these model types were shown to be statistically significantly different from the source model based on manufacturer's data. These results demonstrate that an equivalent source model based only on measured values can be used in place of manufacturer's data for Monte Carlo simulations for MDCT dosimetry.

Variability of surface and center position radiation dose in MDCT: Monte Carlo simulations using CTDI and anthropomorphic phantoms.

The larger coverage afforded by wider z-axis beams in multidetector CT (MDCT) creates larger cone angles and greater beam divergence, which results in substantial surface dose variation for helical and contiguous axial scans. This study evaluates the variation of absorbed radiation dose in both cylindrical and anthropomorphic phantoms when performing helical or contiguous axial scans. The approach used here was to perform Monte Carlo simulations of a 64 slice MDCT. Simulations were performed with different radiation profiles (simulated beam widths) for a given collimation setting (nominal beam width) and for different pitch values and tube start angles. The magnitude of variation at the surface was evaluated under four different conditions: (a) a homogeneous CTDI phantom with different combinations of pitch and simulated beam widths, (b) a heterogeneous anthropomorphic phantom with one measured beam collimation and various pitch values, (c) a homogeneous CTDI phantom with fixed beam collimation and pitch, but with different tube start angles, and (d) pitch values that should minimize variations of surface dose-evaluated for both homogeneous and heterogeneous phantoms. For the CTDI phantom simulations, peripheral dose patterns showed variation with percent ripple as high as 65% when pitch is 1.5 and simulated beam width is equal to the nominal collimation. For the anterior surface dose on an anthropomorphic phantom, the percent ripple was as high as 40% when the pitch is 1.5 and simulated beam width is equal to the measured beam width. Low pitch values were shown to cause beam overlaps which created new peaks. Different x-ray tube start angles create shifts of the peripheral dose profiles. The start angle simulations showed that for a given table position, the surface dose could vary dramatically with minimum values that were 40% of the peak when all conditions are held constant except for the start angle. The last group of simulations showed that an "ideal" pitch value can be determined which reduces surface dose variations, but this pitch value must take into account the measured beam width. These results reveal the complexity of estimating surface dose and demonstrate a range of dose variability at surface positions for both homogeneous cylindrical and heterogeneous anthropomorphic phantoms. These findings have potential implications for small-sized dosimeter measurements in phantoms, such as with TLDs or small Farmer chambers.

Reducing radiation dose to selected organs by selecting the tube start angle in MDCT helical scans: a Monte Carlo based study.

PURPOSE: Previous work has demonstrated that there are significant dose variations with a sinusoidal pattern on the peripheral of a CTDI 32 cm phantom or on the surface of an anthropomorphic phantom when helical CT scanning is performed, resulting in the creation of "hot" spots or "cold" spots. The purpose of this work was to perform preliminary investigations into the feasibility of exploiting these variations to reduce dose to selected radiosensitive organs solely by varying the tube start angle in CT scans. METHODS: Radiation dose to several radiosensitive organs (including breasts, thyroid, uterus, gonads, and eye lenses) resulting from MDCT scans were estimated using Monte Carlo simulation methods on voxelized patient models, including GSF's Baby, Child, and Irene. Dose to fetus was also estimated using four pregnant female models based on CT images of the pregnant patients. Whole-body scans were simulated using 120 kVp, 300 mAs, both 28.8 and 40 mm nominal collimations, and pitch values of 1.5, 1.0, and 0.75 under a wide range of start angles (0 degree-340 degrees in 20 degrees increments). The relationship between tube start angle and organ dose was examined for each organ, and the potential dose reduction was calculated. RESULTS: Some organs exhibit a strong dose variation, depending on the tube start angle. For small peripheral organs (e.g., the eye lenses of the Baby phantom at pitch 1.5 with 40 mm collimation), the minimum dose can be 41% lower than the maximum dose, depending on the tube start angle. In general, larger dose reductions occur for smaller peripheral organs in smaller patients when wider collimation is used. Pitch 1.5 and pitch 0.75 have different mechanisms of dose reduction. For pitch 1.5 scans, the dose is usually lowest when the tube start angle is such that the x-ray tube is posterior to the patient when it passes the longitudinal location of the organ. For pitch 0.75 scans, the dose is lowest when the tube start angle is such that the x-ray tube is anterior to the patient when it passes the longitudinal location of the organ. CONCLUSIONS: Helical MDCT scanning at pitch 1.5 and pitch 0.75 results in "cold spots" and "hot spots" that are created both at surface and in-depth locations within patients. For organs that have a relatively small longitudinal extent, dose can vary considerably with different start angles. While current MDCT systems do not provide the user with the ability to control the tube start angle, these results indicate that in these specific situations (pitch 1.5 or pitch 0.75, small organs and especially small patients), there could be significant dose savings to organs if that functionality would be provided.