Monte Carlo simulation of the dose distribution of ICRP adult reference computational phantoms for acquisitions with a 320 detector-row cone-beam CT scanner

PurposeThe purpose of this study was to develop and validate a Monte Carlo (MC) simulation tool for patient dose assessment for a 320 detector-row CT scanner, based on the recommendations of International Commission on Radiological Protection (ICRP). Additionally, the simulation was applied on four clinical acquisition protocols, with and without automatic tube current modulation (TCM).MethodsThe MC simulation was based on EGS4 code and was developed specifically for a 320 detector-row cone-beam CT scanner. The ICRP adult reference phantoms were used as patient models. Dose measurements were performed free-in-air and also in four CTDI phantoms: 150 mm and 350 mm long CT head and CT body phantoms. The MC program was validated by comparing simulations results with these actual measurements acquired under the same conditions. The measurements agreed with the simulations across all conditions within 5%. Patient dose assessment was performed for four clinical axial acquisitions using the ICRP adult reference phantoms, one of them using TCM.ResultsThe results were nearly always lower than those obtained from other dose calculator tools or published in other studies, which were obtained using mathematical phantoms in different CT systems. For the protocol with TCM organ doses were reduced by between 28 and 36%, compared to the results obtained using a fixed mA value.ConclusionsThe developed simulation program provides a useful tool for assessing doses in a 320 detector-row cone-beam CT scanner using ICRP adult reference computational phantoms and is ready to be applied to more complex protocols.     

A Monte Carlo simulation for the estimation of patient dose in rest and stress cardiac computed tomography with a 320-detector row CT scanner

PurposeTo estimate organ dose and effective dose for patients for cardiac CT as applied in an international multicenter study (CORE320) with a 320-Detector row CT scanner using Monte Carlo (MC) simulations and voxelized phantoms. The effect of positioning of the arms, off-centering the patient and heart rate on patient dose was analyzed.MethodsA MC code was tailored to simulate the geometry and characteristics of the CT scanner. The phantoms representing the adult reference male and female were implemented according to ICRP 110. Effective dose and organ doses were obtained for CT acquisition protocols for calcium scoring, coronary angiography and myocardial perfusion.ResultsFor low heart rate, the normalized effective dose (E) for cardiac CT was higher for female (5.6 mSv/100 mAs) compared to male (2.2 mSv/100 mAs) due to the contribution of female breast tissue. Averaged E for female and male was 11.3 mSv for the comprehensive cardiac protocol consisting of calcium scoring (1.9 mSv); coronary angiography including rest cardiac perfusion (5.1 mSv) and stress cardiac perfusion (4.3 mSv). These values almost doubled at higher heart rates (20.1 mSv). Excluding the arms increased effective dose by 6–8%, centering the patient showed no significant effect. The k-factor (0.028 mSv/mGy.cm) derived from this study leads to effective doses up to 2–3 times higher than the values obtained using now outdated methodologies.ConclusionMC modeling of cardiac CT examinations on realistic voxelized phantoms allowed us to assess patient doses accurately and we derived k-factors that are well above those published previously.     

Fast Monte Carlo codes for occupational dosimetry in interventional radiology

PurposeInterventional radiology techniques cause radiation exposure both to patient and personnel. The radiation dose to the operator is usually measured with dosimeters located at specific points above or below the lead aprons. The aim of this study is to develop and validate two fast Monte Carlo (MC) codes for radiation transport in order to improve the assessment of individual doses in interventional radiology. The proposed methodology reduces the number of required dosemeters and provides immediate dose results.MethodsTwo fast MC simulation codes, PENELOPE/penEasyIR and MCGPU-IR, have been developed. Both codes have been validated by comparing fast MC calculations with the multipurpose PENELOPE MC code and with measurements during a realistic interventional procedure.ResultsThe new codes were tested with a computation time of about 120 s to estimate operator doses while a standard simulation needs several days to obtain similar uncertainties. When compared with the standard calculation in simple set-ups, MCGPU-IR tends to underestimate doses (up to 5%), while PENELOPE/penEasyIR overestimates them (up to 18%). When comparing both fast MC codes with experimental values in realistic set-ups, differences are within 25%. These differences are within accepted uncertainties in individual monitoring.ConclusionThe study highlights the fact that computational dosimetry based on the use of fast MC codes can provide good estimates of the personal dose equivalent and overcome some of the limitations of occupational monitoring in interventional radiology. Notably, MCGPU-IR calculates both organ doses and effective dose, providing a better estimate of radiation risk.     

Monte Carlo simulation of eye lens dose reduction from CT scan using organ based tube current modulation

PurposeTo investigate lens dose reduction with organ based tube current modulation (TCM) using the Monte Carlo method.MethodsTo calculate lens dose with organ based TCM, 36 pairs of X-ray sources with bowtie filters were placed around the patient head using a projection angle interval of 10° for one rotation of Computed Tomography (CT). Each projection was simulated respectively. Both voxelized and stylized eye models and Chinese reference male phantoms were used in the simulation, and tube voltages 80, 100, 120 and 140 kVp were used.ResultsDose differences between two eye models were less than 20%, but large variations were observed among dose results from different projections of all tube voltages investigated. Dose results from 0° (AP) directions were 60 times greater than those from 180° (PA) directions, which enables organ based TCM reduce lens doses by more than 47%.ConclusionsOrgan based TCM may be used to reduce lens doses. Stylized eye models are more anatomically realistic compared with voxelized eye models and are more reliable for dose evaluation.     

Shielding effectiveness of X-ray protective garment

PurposeCertification of the X-ray shielding garment is based on attenuation testing on flat material samples. We investigated the difference of shielding effectiveness compared to realistic use when the garment is worn on the body of a staff person.MethodsAttenuation factors of X-ray protective aprons have been evaluated for several clinical scenarios with Monte Carlo (MC) calculations based on the ICRP female reference model and an experimental setup. The MC calculated attenuation factors refer to the effective dose E, whereas the measured attenuation factors refer to the personal dose equivalent Hp(10). The calculated/measured factors were compared to the attenuation factors of the identical materials measured under the conditions of the standard IEC 61331-1 that is currently in use for the type testing of X-ray protective aprons.ResultsAs a result, for example, at a common tube voltage of 80 kV, the real attenuation factors of a 0.35 mm Pb apron worn by a 3-dimensional body were 38% to 76% higher than when measured under IEC conditions on flat samples. The MC-calculated organ doses show the maximum contribution to E being within the operatoŕs abdomen/pelvis region.ConclusionsWith our findings, personal X-ray protective garments could be improved in effectiveness.     

Monte Carlo as quality control tool of stereotactic body radiation therapy treatment plans

Purpose/objectiveThe objective of this study was to verify the accuracy of treatment plans of stereotactic body radiation therapy (SBRT) and to verify the feasibility of the use of Monte Carlo (MC) as quality control (QC) on a daily basis.Material/methodsUsing EGSnrc, a MC model of Agility™ linear accelerator was created. Various measurements (Percentage depth dose (PDD), Profiles and Output factors) were done for different fields sizes from 1x1 up to 40x40 (cm²). An iterative model optimization was performed to achieve adequate parameters of MC simulation. 40 SBRT patient’s dosimetry plans were calculated by Monaco™ 3.1.1. CT images, RT-STRUCT and RT-PLAN files from Monaco™ being used as input for Moderato MC code. Finally, dose volume histogram (DVH) and paired t-tests for each contour were used for dosimetry comparison of the Monaco™ and MC.ResultsValidation of MC model was successful, as <2% difference comparing to measurements for all field’s sizes. The main energy of electron source incident on the target was 5.8 MeV, and the full width at half maximum (FWHM) of Gaussian electron source were 0.09 and 0.2 (cm) in X and Y directions, respectively. For 40 treatment plan comparisons, the minimum absolute difference of mean dose of planning treatment planning (PTV) was 0.1% while the maximum was 6.3%. The minimum absolute difference of Max dose of PTV was 0.2% while the maximum was 8.1%.ConclusionSBRT treatment plans of Monaco agreed with MC results. It possible to use MC for treatment plans verifications as independent QC tool.     

A new analytical formula for neutron capture gamma dose calculations in double-bend mazes in radiation therapy

BackgroundPhotoneutrons are produced in radiation therapy with high energy photons. Also, capture gamma rays are the byproduct of neutrons interactions with wall material of radiotherapy rooms.AimIn the current study an analytical formula was proposed for capture gamma dose calculations in double bend mazes in radiation therapy rooms.Materials and methodsA total of 40 different layouts with double-bend mazes and a 18 MeV photon beam of Varian 2100 Clinac were simulated using MCNPX Monte Carlo (MC) code. Neutron capture gamma ray dose equivalent was calculated by the MC method along the maze and at the maze entrance door of all the simulated rooms. Then, all MC resulted data were fitted to an empirical formula for capture gamma dose calculations. Wu–McGinley analytical formula for capture gamma dose equivalent at the maze entrance door in single-bend mazes was also used for comparison purposes.ResultsFor capture gamma dose equivalents at the maze entrance door, the difference of 2–11% was seen between MC and the derived equation, while the difference of 36–87% was found between MC and the Wu–McGinley methods.ConclusionOur results showed that the derived formula results were consistent with the MC results for all of 40 different geometries. However, as a new formula, further evaluations are required to validate its use in practical situations. Finally, its application is recommend for capture gamma dose calculations in double-bend mazes to improve shielding calculations.     

Monte Carlo dose in a prosthesis phantom based on exact geometry vs streak artefact contaminated CT data as benchmarked against Gafchromic film measurements

PurposeIn radiotherapy, accurate calculation of patient radiation dose is very important for good clinical outcome. In the presence of metallic implants, the dose calculation accuracy could be compromised by metal artefacts generated in computed tomography (CT) images of patients. This study investigates the influence of metal-induced CT artefacts on MC dose calculations in a pelvic prosthesis phantom.MethodsA pelvic phantom containing unilateral Ti prosthesis was CT-scanned and accurate Hounsfield unit (HU) values were assigned to known materials of the phantom as opposed to HU values produced through the artefact CT images of the phantom. Using the DOSXYZnrc MC code, dose calculations were computed in the phantom model constructed from the original CT images containing the artefacts and artefact-free images made from the exact geometry of the phantom with known materials. The dose calculations were benchmarked against Gafchromic EBT3 film measurements using 15 MeV electron and 10 MV photon beams.ResultsThe average deviations between film and MC dose data decreased from 3 ± 2% to 1 ± 1% and from about 6 ± 2% to 3 ± 1% for the artefact and artefact-free phantom models against film data for the electron and photon fields, respectively.ConclusionsFor the Ti prosthesis phantom, the presence of metal-induced CT artefacts could cause dose inaccuracies of about 3%. Construction of an artefact-free phantom model made from the exact geometry of the phantom with known materials to overcome the effect of artefacts is advantageous compared to using CT data directly of which the exact tissue composition is not well-known.     

Validation of irtGPUMCD,a GPU-based Monte Carlo internal dosimetry framework for radionuclide therapy

PurposeMonte Carlo (MC) simulations are highly desirable for dose treatment planning and evaluation in radiation oncology. This is true also in emerging nuclear medicine applications such as internal radiotherapy with radionuclides. The purpose of this study is the validation of irtGPUMCD, a GPU-based MC code for dose calculations in internal radiotherapy.MethodsThe female and male phantoms of the International Commission on Radiological Protection (ICRP 110) were used as benchmarking geometries for this study focused on ¹⁷⁷Lu and including ^99mTc and ¹³¹I. Dose calculations were also conducted for a real patient. For phantoms, twelve anatomical structures were considered as target/source organs. The S-values were evaluated with irtGPUMCD simulations (10⁸ photons), with gamma branching ratios of ICRP 107 publication. The ¹⁷⁷Lu electrons S-values were calculated for source organs only, based on local deposition of dose in irtGPUMCD. The S-value relative difference between irtGPUMCD and IDAC-DOSE were evaluated for all targets/sources considered. A DVHs comparison with GATE was conducted. An exponential track length estimator was introduced in irtGPUMCD to increase computational efficiency.ResultsThe relative S-value differences between irtGPUMCD and IDAC-DOSE were <5% while this comparison with GATE was <1%. The DVHs dosimetric indices comparison between GATE and irtGPUMCD for the patient led to an excellent agreement (<2%). The time required for the simulation of 10⁸ photons was 1.5 min for the female phantom, and one minute for the real patient (<1% uncertainty). These results are promising and let envision the use of irtGPUMCD for internal dosimetry in clinical applications.     

An automated Monte Carlo QC system for volumetric modulated arc therapy: Possibilities and challenges

PurposeTo develop and implement an automated Monte Carlo (MC) system for patient specific VMAT quality control in a patient geometry that generates treatment planning system (TPS) compliant DICOM objects and includes a module for 3D analysis of dose deviations. Also, the aims were to recommend diagnose specific tolerance criteria and an evaluation procedure.MethodsThe EGSnrc code package formed the basis for development of the MC system. The workflow consists of a number of modules connected to a TPS by means of manual DICOM exports and imports which were executed sequentially without user interaction. DVH comparison was performed in the TPS. In addition, MC- and TPS dose distributions were analysed by applying the normalized dose difference (NDD) formalism. NDD failure maps and a pass rate for a certain threshold were obtained. 170 clinical plans (prostate, thorax, head-and-neck and gynecological) were selected for analysis.ResultsAgreement within 1.5% was found between clinical- and MC data for the mean dose to the target volumes and within 3% for parameters more sensitive to the shape of the DVH e.g. D_98% PTV. Regarding the NDD analysis, tolerance criteria 2%/3 mm were established for prostate plans and 3%/3 mm for the rest of the cases.ConclusionsAn automated MC system was developed and implemented. Evaluation procedure is recommended with NDD-analysis as a first step. For pass rate < 95%, the evaluation continues with comparison of DVH parameters. For deviations larger than 2%, a visual inspection of the clinical- and MC dose distributions is performed.     

Simplified patient-specific renal dosimetry in 177Lu therapy: a proof of concept

PurposeThe aim of this proof-of-concept study is to propose a simplified personalized kidney dosimetry procedure in ¹⁷⁷Lu peptide receptor radionuclide therapy (PRRT) for neuroendocrine tumors and metastatic prostate cancer. It relies on a single quantitative SPECT/CT acquisition and multiple radiometric measurements executed with a collimated external probe, properly directed on kidneys.MethodsWe conducted a phantom study involving external count-rate measurements in an abdominal phantom setup filled with activity concentrations of ^99mTc, reproducing patient-relevant organ effective half-lives occurring in ¹⁷⁷Lu PRRT. GATE Monte Carlo (MC) simulations of the experiment, using ^99mTc and ¹⁷⁷Lu as sources, were performed. Furthermore, we tested this method via MC on a clinical case of ¹⁷⁷Lu-DOTATATE PRRT with SPECT/CT images at three time points (2, 20 and 70 hrs), comparing a simplified kidney dosimetry, employing a single SPECT/CT and probe measurements at three time points, with the complete MC dosimetry.ResultsThe experimentally estimated kidney half-life with background subtraction applied was compatible within 3% with the expected value. The MC simulations of the phantom study, both with ^99mTc and ¹⁷⁷Lu, confirmed a similar level of accuracy. Concerning the clinical case, the simplified dosimetric method led to a kidney dose estimation compatible with the complete MC dosimetry within 6%, 12% and 2%, using respectively the SPECT/CT at 2, 20 and 70 hrs.ConclusionsThe proposed simplified procedure provided a satisfactory accuracy and would reduce the imaging required to derive the kidney absorbed dose to a unique quantitative SPECT/CT, with consequent benefits in terms of clinic workflows and patient comfort.     

Fetal radiation dose in three common CT examinations during pregnancy – Monte Carlo study

PurposeTo determine fetal doses in different stages of pregnancy in three common computed tomography (CT) examinations: pulmonary CT angiography, abdomino-pelvic and trauma scan with Monte Carlo (MC) simulations.MethodsAn adult female anthropomorphic phantom was scanned with a 64-slice CT using pulmonary angiography, abdomino-pelvic and trauma CT scan protocols. Three different sized gelatin boluses placed on the phantom’s abdomen simulated different stages of pregnancy. Intrauterine dose was used as a surrogate to a dose absorbed to the fetus. MC simulations were performed to estimate uterine doses. The simulation dose levels were calibrated with volumetric CT dose index (CTDI_vol) measurements and MC simulations in a cylindrical CTDI body phantom and compared with ten point doses measured with metal-oxide-semiconductor field-effect-transistor dosimeters. Intrauterine volumes and uterine walls were segmented and the respective dose volume histograms were calculated.ResultsThe mean intrauterine doses in different stages of pregnancy varied from 0.04 to 1.04 mGy, from 4.8 to 5.8 mGy, and from 9.8 to 12.6 mGy in the CT scans for pulmonary angiography, abdomino-pelvic and trauma CT scans, respectively. MC simulations showed good correlation with the MOSFET measurement at the measured locations.ConclusionsThe three studied examinations provided highly varying fetal doses increasing from sub-mGy level in pulmonary CT angiography to notably higher levels in abdomino-pelvic and trauma scans where the fetus is in the primary exposure range. Volumetric dose distribution offered by MC simulations in an appropriate anthropomorphic phantom provides a comprehensive dose assessment when applied in adjunct to point-dose measurements.     

Total body irradiation for standard treatment rooms: a robust sweeping beam technique with respect to the body shape

BackgroundThe purpose of this work is to improve a sweeping beam technique for total body irradiation (TBI) on a low flat couch using a varying patient thickness model. We designed a flat couch for total body irradiation in supine and prone position. Three generic arcs with rectangular segments for a patient torso thickness of 16, 22 and 28 cm were generated with respect to varying patient thickness of four particular parts of the body: head, torso, thighs and calves.Materials and methodsLongitudinal and transversal dose profiles were measured using an ionization chamber and the EBT3 gafchromic film in a solid water slab phantom. The robustness of the method was examined in phantoms of different thicknesses.ResultsMeasured dose homogeneity stays within ±10% of prescribed dose for all of the three patient thickness models. The robustness of the method was evaluated as the increase in dose in the phantom center of 0.7% per 1 cm reduction in phantom thickness.ConclusionThe method is applicable for the broad range of patient sizes, comfortable for patients, robust and suitable for standard treatment rooms with a standard linear accelerator. It requires minimal investments into equipment.     

Generating and using patient-specific whole-body models for organ dose estimates in CT with increased accuracy: Feasibility and validation

The estimation of patient dose using Monte Carlo (MC) simulations based on the available patient CT images is limited to the length of the scan. Software tools for dose estimation based on standard computational phantoms overcome this problem; however, they are limited with respect to taking individual patient anatomy into account. The purpose of this study was to generate whole-body patient models in order to take scattered radiation and over-scanning effects into account. Thorax examinations were performed on three physical anthropomorphic phantoms at tube voltages of 80 kV and 120 kV; absorbed dose was measured using thermoluminescence dosimeters (TLD). Whole-body voxel models were built as a combination of the acquired CT images appended by data taken from widely used anthropomorphic voxel phantoms. MC simulations were performed both for the CT image volumes alone and for the whole-body models. Measured and calculated dose distributions were compared for each TLD chip position; additionally, organ doses were determined.MC simulations based only on CT data underestimated dose by 8%–15% on average depending on patient size with highest underestimation values of 37% for the adult phantom at the caudal border of the image volume. The use of whole-body models substantially reduced these errors; measured and simulated results consistently agreed to better than 10%.This study demonstrates that combined whole-body models can provide three-dimensional dose distributions with improved accuracy. Using the presented concept should be of high interest for research studies which demand high accuracy, e.g. for dose optimization efforts.     

Evaluation of patient organ doses from kilovoltage cone-beam CT imaging in radiation therapy

BackgroundCurrently, CBCT system is an indispensable component of radiation therapy units. Because of that, it is important in treatment planning and diagnosis. CBCT is also an crucial tool for patient positioning and verification in image-guided radiation therapy (IGRT). Therefore, it is critical to investigate the patient organ doses arising from CBCT imaging. The purpose of this study is to evaluate patient organ doses and effective dose to patients from three different protocols of Elekta Synergy XVI system for kV CBCT imaging examinations in image guided radiation therapy.Materials and methodsOrgan dose measurements were done with thermoluminescent dosimeters in Alderson RA NDO male phantom for head & neck (H&N), chest and pelvis protocols of the Elekta Synergy XVI kV CBCT system. From the measured organ dose, effective dose to patients were calculated according to the International Commission on Radiological Protection 103 report recommendations.ResultsFor H&N, chest and pelvis scans, the organ doses were in the range of 0.03–3.43 mGy, 6.04–22.94 mGy and 2.5–25.28 mGy, respectively. The calculated effective doses were 0.25 mSv, 5.56 mSv and 4.72 mSv, respectively.ConclusionThe obtained results were consistent with the most published studies in the literature. Although the doses to patient organs from the kV CBCT system were relatively low when compared with the prescribed treatment dose, the amount of delivered dose should be monitored and recorded carefully in order to avoid secondary cancer risk, especially in pediatric examinations.     

Dose distribution comparison in volumetric-modulated arc therapy plans for head and neck cancers with and without an external body contour extended technique

AimThis study compared volumetric-modulated arc therapy (VMAT) plans for head and neck cancers with and without an external body contour extended technique (EBCT).BackgroundDose calculation algorisms for VMAT have limitations in the buildup region.Materials and methodsThree VMAT plans were enrolled, with one case having a metal artifact from an artificial tooth. The proper dose was calculated using Eclipse version 11.0. The body contours were extended 2 cm outward from the skin surface in three-dimensional space, and the dose was recalculated with an anisotropic analytical algorithm (AAA) and Acuros XB (AXB). Monitor units (MUs) were set, and the dose distributions in the planning target volume (PTV), clinical target volume, and organ at risk (OAR) and conformity index (CI) with and without an EBCT were compared. The influence of a metal artifact outside of the thermoplastic head mask was also compared.ResultsThe coverage of PTV by the 95% dose line near the patient’s skin was increased drastically by using an EBCT. Plan renormalization had a negligible impact on MUs and doses delivered to OARs. CI of PTV with a 6-MV photon beam was closer to 1 than that with a 10-MV photon beam when both AAA and AXB were used in all cases. Metal artifacts outside the head mask had no effect on dose distribution.ConclusionsAn EBCT is needed to estimate the proper dose at object volumes near the patient’s skin and can improve the accuracy of the calculated dose at target volumes.     

Multiple-source models for electron beams of a medical linear accelerator using BEAMDP computer code

AimThe aim of this work was to develop multiple-source models for electron beams of the NEPTUN 10PC medical linear accelerator using the BEAMDP computer code.BackgroundOne of the most accurate techniques of radiotherapy dose calculation is the Monte Carlo (MC) simulation of radiation transport, which requires detailed information of the beam in the form of a phase-space file. The computing time required to simulate the beam data and obtain phase-space files from a clinical accelerator is significant. Calculation of dose distributions using multiple-source models is an alternative method to phase-space data as direct input to the dose calculation system.Materials and methodsMonte Carlo simulation of accelerator head was done in which a record was kept of the particle phase-space regarding the details of the particle history. Multiple-source models were built from the phase-space files of Monte Carlo simulations. These simplified beam models were used to generate Monte Carlo dose calculations and to compare those calculations with phase-space data for electron beams.ResultsComparison of the measured and calculated dose distributions using the phase-space files and multiple-source models for three electron beam energies showed that the measured and calculated values match well each other throughout the curves.ConclusionIt was found that dose distributions calculated using both the multiple-source models and the phase-space data agree within 1.3%, demonstrating that the models can be used for dosimetry research purposes and dose calculations in radiotherapy.