Cancer risk estimation in Digital Breast Tomosynthesis using GEANT4 Monte Carlo simulations and voxel phantoms

The aim of this work was to estimate the risk of radiation induced cancer following the Portuguese breast screening recommendations for Digital Mammography (DM) when applied to Digital Breast Tomosynthesis (DBT) and to evaluate how the risk to induce cancer could influence the energy used in breast diagnostic exams. The organ doses were calculated by Monte Carlo simulations using a female voxel phantom and considering the acquisition of 25 projection images. Single organ cancer incidence risks were calculated in order to assess the total effective radiation induced cancer risk. The screening strategy techniques considered were: DBT in Cranio-Caudal (CC) view and two-view DM (CC and Mediolateral Oblique (MLO)).The risk of cancer incidence following the Portuguese screening guidelines (screening every two years in the age range of 50–80 years) was calculated by assuming a single CC DBT acquisition view as standalone screening strategy and compared with two-view DM. The difference in the total effective risk between DBT and DM is quite low. Nevertheless in DBT an increase of risk for the lung is observed with respect to DM. The lung is also the organ that is mainly affected when non-optimal beam energy (in terms of image quality and absorbed dose) is used instead of an optimal one. The use of non-optimal energies could increase the risk of lung cancer incidence by a factor of about 2.     

A Monte Carlo evaluation for effects of probable dimensional uncertainties of low dose rate brachytherapy seeds on dose

The aim of this study is to determine effects of size deviations of brachytherapy seeds on two dimensional dose distributions around the seed. Although many uncertainties are well known, the uncertainties which stem from geometric features of radiation sources are weakly considered and predicted. Neither TG-43 report which is not completely in common consensus, nor individual scientific MC and experimental studies include sufficient data for geometric uncertainties. Sizes of seed and its components can vary in a manufacturing deviation. This causes geometrical uncertainties, too. In this study, three seeds which have different geometrical properties were modeled using EGSnrc-Code Packages. Seeds were designed with all their details using the geometry package. 5% deviations of seed sizes were assumed. Modified seeds were derived from original seed by changing sizes by 5%. Normalizations of doses which were calculated from three kinds of brachytherapy seed and their derivations were found to be about 3%–20%. It was shown that manufacturing differences of brachytherapy seed cause considerable changes in dose distribution.     

A depth dose study between AAA and AXB algorithm against Monte Carlo simulation using AIP CT of a 4D dataset from a moving phantom

Aim

To identifying depth dose differences between the two versions of the algorithms using AIP CT of a 4D dataset.

Monte Carlo studies in Gold Nanoparticles enhanced radiotherapy: The impact of modelled parameters in dose enhancement

PurposeOver the last decades, Gold Nanoparticles (AuNPs) have been presented as an innovative approach in radiotherapy (RT) enhancement. Several studies have proven that the irradiation of tumors containing AuNPs could lead to more effective tumor control than irradiation alone. Studies with low kV photons and AuNPs conclude in encouraging results regarding the level of radioenhancement. However, experimental and theoretical studies with MV photons report controversial findings concerning the correlation between dose enhancement effect and tumor cell killing. The great variation in the experimental protocols and simulations complicates the comparison of their outcomes and depicts the need for limiting the variety of investigated parameters. Our purpose is to point out a possible direction for building realistic Monte Carlo (MC) models that could end up with promising results in MV photons RT enhancement.MethodsWe explored published in silico studies concerning AuNPs enhanced RT from 2010 to 2019. In this review, we discuss the different AuNPs and MV photon beams characteristics that have been reported and their effect in dose enhancement.ResultsAuNPs size, concentration, type of distribution along with photon beams energy and the presence of flattening filter in linear accelerators seem to be the major parameters that determine AuNPs radioenhancement in silico.ConclusionsPrior to AuNPs clinical translation in photon radiotherapy, in silico studies should emphasize on nanodosimetry and track structure codes than condensed history ones. Toxicity estimation and biological aspects should be implemented in MC simulations so as to achieve accurate and realistic modelling of AuNPs driven RT.     

Uncertainty analysis of penicillin V production using Monte Carlo simulation

Uncertainty and variability affect economic and environmental performance in the production of biotechnology and pharmaceutical products. However, commercial process simulation software typically provides analysis that assumes deterministic rather than stochastic process parameters and thus is not capable of dealing with the complexities created by variance that arise in the decision-making process. Using the production of penicillin V as a case study, this article shows how uncertainty can be quantified and evaluated. The first step is construction of a process model, as well as analysis of its cost structure and environmental impact. The second step is identification of uncertain variables and determination of their probability distributions based on available process and literature data. Finally, Monte Carlo simulations are run to see how these uncertainties propagate through the model and affect key economic and environmental outcomes. Thus, the overall variation of these objective functions are quantified, the technical, supply chain, and market parameters that contribute most to the existing variance are identified and the differences between economic and ecological evaluation are analyzed. In our case study analysis, we show that final penicillin and biomass concentrations in the fermenter have the highest contribution to variance for both unit production cost and environmental impact. The penicillin selling price dominates return on investment variance as well as the variance for other revenue-dependent parameters.     

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 calculations of radiotherapy dose in “homogeneous” anatomy

Given the substantial literature on the use of Monte Carlo (MC) simulations to verify treatment planning system (TPS) calculations of radiotherapy dose in heterogeneous regions, such as head and neck and lung, this study investigated the potential value of running MC simulations of radiotherapy treatments of nominally homogeneous pelvic anatomy. A pre-existing in-house MC job submission and analysis system, built around BEAMnrc and DOSXYZnrc, was used to evaluate the dosimetric accuracy of a sample of 12 pelvic volumetric arc therapy (VMAT) treatments, planned using the Varian Eclipse TPS, where dose was calculated with both the Analytical Anisotropic Algorithm (AAA) and the Acuros (AXB) algorithm. In-house TADA (Treatment And Dose Assessor) software was used to evaluate treatment plan complexity, in terms of the small aperture score (SAS), modulation index (MI) and a novel exposed leaf score (ELS/ELA). Results showed that the TPS generally achieved closer agreement with the MC dose distribution when treatments were planned for smaller (single-organ) targets rather than larger targets that included nodes or metastases. Analysis of these MC results with reference to the complexity metrics indicated that while AXB was useful for reducing dosimetric uncertainties associated with density heterogeneity, the residual TPS dose calculation uncertainties resulted from treatment plan complexity and TPS model simplicity. The results of this study demonstrate the value of using MC methods to recalculate and check the dose calculations provided by commercial radiotherapy TPSs, even when the treated anatomy is assumed to be comparatively homogeneous, such as in the pelvic region.     

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.     

Dual-energy contrast-enhanced digital mammography: Glandular dose estimation using a Monte Carlo code and voxel phantom

PurposeTo estimate the mean glandular dose of contrast enhanced digital mammography, using the EGSnrc Monte Carlo code and female adult voxel phantom.MethodsAutomatic exposure control of full field digital mammography system was used for the selection of the X-ray spectrum and the exposure settings for dual energy imaging. Measurements of the air-kerma and of the half value layers were performed and a Monte Carlo simulation of the digital mammography system was used to compute the mean glandular dose, for breast phantoms of various thicknesses, glandularities and for different X-ray spectra (low and high energy).ResultsFor breast phantoms of 2.0–8.0 cm thick and 0.1–100% glandular fraction, CC view acquisition, from AEC settings, can result in a mean glandular dose of 0.450 ± 0.022 mGy −2.575 ± 0.033 mGy for low energy images and 0.061 ± 0.021 mGy – 0.232 ± 0.033 mGy for high energy images. In MLO view acquisition mean glandular dose values ranged between 0.488 ± 0.007 mGy – 2.080 ± 0.021 mGy for low energy images and 0.065 ± 0.012 mGy – 0.215 ± 0.010 mGy for high energy images.ConclusionThe low kV part of contrast enhanced digital mammography is the main contributor to total mean glandular breast dose. The results of this study can be used to provide an estimated mean glandular dose for individual cases.     

Gold coated contact lens-type ocular in vivo dosimeter (CLOD) for monitoring of low dose in computed tomography: A Monte Carlo study

PurposeThis study reports a sensitivity enhancement of gold-coated contact lens-type ocular in vivo dosimeters (CLODs) for low-dose measurements in computed tomography (CT).MethodsMonte Carlo (MC) simulations were conducted to evaluate the dose enhancement from the gold (Au) layers on the CLODs. The human eye and CLODs were modeled, and the X-ray tube voltages were defined as 80, 120, and 140 kVp. The thickness of the Au layer attached to a CLOD ranged from 100 nm to 10 μm. The thickness of the active layer ranged from 20 to 140 μm. The dose ratio between the active layer of the Au-coated CLOD and a CLOD without a layer, i.e., the dose enhancement factor (DEF), was calculated.ResultsThe DEFs of the first 20-μm thick active layer of the 5-μm thick Au-coated CLOD were 18.4, 19.7, 20.2 at 80, 120, and 140 kVp, respectively. The DEFs decreased as the thickness of the active layer increased. The DEFs of 100-nm to 5-μm thick Au layers increased from 1.7 to 5.4 for 120-kVp X-ray tube voltage when the thickness of the active layer was 140 μm.ConclusionsThe MC results presented a higher sensitivity of Au-coated CLODs (∼20-times higher than that of CLODs without a gold layer). Au-coated CLODs can be applied to an evaluation of very low doses (a few cGy) delivered to patients during CT imaging.     

Monte Carlo dose verification of VMAT treatment plans using Elekta Agility 160-leaf MLC

In this study, we verified volumetric modulated arc therapy (VMAT) plans in an Elekta Synergy system with an integrated Agility 160-leaf multileaf collimator (MLC) by comparing them with Monte Carlo (MC)-calculated dose distributions using the AAPM TG-119 structure sets. The head configuration of the linear accelerator with the integrated MLC was simulated with the EGSnrc/BEAMnrc code. Firstly, the dosimetric properties of the MLC were evaluated with the MC technique and film measurements. Next, VMAT plans were created with the Pinnacle³ treatment planning system (TPS) for four regions in the AAPM TG-119 structures. They were then verified by comparing them with MC-calculated dose distributions using dose volume histograms (DVHs) and three-dimensional (3D) gamma analysis. The MC simulations for the Agility MLC dosimetric properties were in acceptable agreement with measurements. TPS-VMAT plans using TG-119 structure sets agreed with MC dose distributions within 2% in the comparison of D₉₅ in planning target volumes (PTVs) evaluated from DVHs. In contrast, higher dose regions such as D₂₀, D₁₀, and D₅ in PTVs for TPS tended to be smaller than MC values. This tendency was particularly noticeable for mock head and neck with complicated structures. In 3D gamma analysis, the passing rates with 3%/3mm criteria in PTVs were ≥99%, except for mock head and neck (89.5%). All passing rates for organs at risk (OARs) were in acceptable agreement of >96%. It is useful to verify dose distributions of PTVs and OARs in TPS-VMAT plans by using MC dose calculations and 3D gamma analysis.     

Mobile shielding evaluation on the fetal dose during a breast radiotherapy using Monte Carlo simulation

PurposeEvaluation of the out-of-field dose is an important aspect in radiotherapy. Due to the fetus radiosensitivity, this evaluation becomes even more conclusive when the patient is pregnant. In this work, a linear accelerator Varian Clinac 2100c operating at 6 MV, a pregnant anthropomorphic phantom (Maria), and different shields added above the abdominal region of the phantom were used for the analysis based on MCNPX. Methods: The simulations were performed for the medial and lateral projections, using either an open field collimation (10×16 cm²) or a multileaf collimator. The added shields (M1 and M2) were designed based on models proposed by Stovall et al. [1], intending to reduce the deposited dose on the fetus and related structures. Results: The presence of the shields showed to be effective in reducing the doses on the fetus, amniotic sac, and placenta, for example. A reduction of about 43% was found in the dose on the fetus when M2 was added, using the open field collimation, in comparison with the situation with no shield, being the lateral projection the main responsible for the dose. The use of MLC significatively reduced the doses in different structures, including on the fetus and amniotic sac, for example, in comparison to the open field situation. A slight increment on the dose in organs such as the eyes, thyroid and brain was found in both collimation systems, due to the presence of the shields. The contribution of the leakage radiation from the tube head of the linear accelerator was found to be in the order of µGy, being reduced by the presence of the M2 shield. Conclusion: Using the shields showed to be an essential feature in order to reduce the dose not only on the fetus, but also in important structures responsible to its development.     

Vacancy clustering and diffusion in germanium using kinetic lattice Monte Carlo simulations

We investigated vacancy-assisted self-diffusion in germanium by means of kinetic lattice Monte Carlo (KLMC) simulations below the melting temperature, for a vacancy concentration of 1 × 10¹⁸/cm³. At higher temperatures, fewer clusters formed, but there was less variation in the number of clusters than at lower temperatures as the time increased. Equilibrium diffusivities in the clustering region were 10² lower than those of free vacancies in the initial stage of KLMC simulations. They were expressed according to three temperature regimes: 6.5 × 10^? 4 exp(–0.35/k _B T) cm²/s at temperatures above 1100 K, 5.2 × 10⁵ exp(–2.32/k _B T) cm²/s at temperatures of 900–1100 K and 6.0 × 0^–7 exp(–0.19/k _B T) cm²/s at temperatures below 900 K. The effective mean migration energy, 1.1 eV, closely coincided with that of the 1.0–1.2 eV in experiments and was very different from the migration energy of the free vacancy.     

Validation of 4D Monte Carlo dose calculations using a programmable deformable lung phantom

PurposeTo validate the accuracy of 4D Monte Carlo (4DMC) simulations to calculate dose deliveries to a deforming anatomy in the presence of realistic respiratory motion traces. A previously developed deformable lung phantom comprising an elastic tumor was modified to enable programming of arbitrary motion profiles. 4D simulations of the dose delivered to the phantom were compared with the measurements.MethodsThe deformable lung phantom moving with irregular breathing patterns was irradiated using static and VMAT beam deliveries. Using the RADPOS 4D dosimetry system, point doses were measured inside and outside the tumor. Dose profiles were acquired using films along the motion path of the tumor (S-I). In addition to dose measurements, RADPOS was used to record the motion of the tumor during dose deliveries. Dose measurements were then compared against 4DMC simulations with EGSnrc/4DdefDOSXYZnrc using the recorded tumor motion.ResultsThe agreements between dose profiles from measurements and simulations were determined to be within 2%/2 mm. Point dose agreements were within 2σ of experimental and/or positional/dose reading uncertainties. 4DMC simulations were shown to accurately predict the sensitivity of delivered dose to the starting phase of breathing motions. We have demonstrated that our 4DMC method, combined with RADPOS, can accurately simulate realistic dose deliveries to a deforming anatomy moving with realistic breathing traces. This 4DMC tool has the potential to be used as a quality assurance tool to verify treatments involving respiratory motion. Adaptive treatment delivery is another area that may benefit from the potential of this 4DMC tool.     

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.     

Intercomparison of dose enhancement ratio and secondary electron spectra for gold nanoparticles irradiated by X-rays calculated using multiple Monte Carlo simulation codes

PurposeTargeted radiation therapy has seen an increased interest in the past decade. In vitro and in vivo experiments showed enhanced radiation doses due to gold nanoparticles (GNPs) to tumors in mice and demonstrated a high potential for clinical application. However, finding a functionalized molecular formulation for actively targeting GNPs in tumor cells is challenging. Furthermore, the enhanced energy deposition by secondary electrons around GNPs, particularly by short-ranged Auger electrons is difficult to measure. Computational models, such as Monte Carlo (MC) radiation transport codes, have been used to estimate the physical quantities and effects of GNPs. However, as these codes differ from one to another, the reliability of physical and dosimetric quantities needs to be established at cellular and molecular levels, so that the subsequent biological effects can be assessed quantitatively.MethodsIn this work, irradiation of single GNPs of 50 nm and 100 nm diameter by X-ray spectra generated by 50 and 100 peak kilovoltages was simulated for a defined geometry setup, by applying multiple MC codes in the EURADOS framework.ResultsThe mean dose enhancement ratio of the first 10 nm-thick water shell around a 100 nm GNP ranges from 400 for 100 kVp X-rays to 600 for 50 kVp X-rays with large uncertainty factors up to 2.3.ConclusionsIt is concluded that the absolute dose enhancement effects have large uncertainties and need an inter-code intercomparison for a high quality assurance; relative properties may be a better measure until more experimental data is available to constrain the models.