RESUMEN

The author reminisces about some of his experiences working with Monte Carlo techniques for Medical Physics applications.

Asunto(s)

RESUMEN

Introduction: EGSnrc software package is one of the computational packages for Monte Carlo simulation in radiation therapy and has several subset codes. Directional bremsstrahlung splitting (DBS) is a technique that applies braking radiations in interactions in this software. This study aimed to evaluate the effect of this technique on the simulation time, uncertainty, particle number of phase-space data, and photon beam spectrum resulting from a medical linear accelerator (LINAC). Materials and methods: The gantry of the accelerator, including the materials and geometries of different parts, was simulated using the BEAMnrc code (a subset code in the EGSnrc package). The phase-space data were recorded in different parts of the LINAC. The DBS values (1, 10, 100, and 1000) were changed, and their effects were evaluated on the simulation parameters and output spectra. Results: Increasing the DBS value from 1 to 1000 resulted in an increase in the simulation time from 1.778 to 11.310 hours, and increasing the number of particles in the phase-space plane (5 590 732-180 328 382). When the DBS had been picked up from 1 to 100, the simulation uncertainty decreased by about 1.29%. In addition, the DBS increment value from 100 to 1000 leads to an increase in uncertainty and simulation time of about 0.71% and 315%, respectively. Conclusion: Although using the DBS technique reduces the simulation time or uncertainty, increasing the DBS from a specific value, equal to 100 in our study, increases simulation uncertainties and times. Therefore, we propose considering a specific DBS value as we obtained for the Monte Carlo simulation of photon beams produced by linear accelerators.

RESUMEN

The availability of geometrical, physical, and initial beam parameters for Monte Carlo (MC) simulations of the Elekta Agility collimator head has become very difficult due to the proprietary nature of this data. This study presents strategies to independently determine the geometrical and physical properties of the components and initial beam parameters of the Agility collimator head for full beam simulations and postulates a benchmarking process using the EGSnrc MC toolkit. Target material of W (90%) and Re (10%) of 0.09 cm thickness, flattening filter of 1.77 cm thick stainless steel placed on 0.5 cm Al disc, and primary and secondary collimators of Tungsten alloy have been found to best fit the Agility head. The initial beam energy of 6.0 MeV with a radial distribution given as full-width half maxima (FWHM) of 0.301 cm (crossline) × 0.201 cm (inline) for 6 MV beam with a mean angular spread of 1.34° has been found best fitting the model. Variations of 0.29% and 0.59% have been noted in the measured and calculated values of TPR20,10 and D10 respectively. More than 90% dose points for all simulations passed the 2D gamma criteria of 3% DD, 3 mm DTA. This MC model of the Agility head can be used for dose calculation and validation of advanced treatment techniques.

Asunto(s)

RESUMEN

BACKGROUND: In this study, the dose distributions obtained by the algorithms used in Monaco treatment planning system (TPS) and Monte Carlo (MC) simulation were compared for small fields in the anthropomorphic RANDO phantom, and then, the results were analyzed using the gamma analysis method. MATERIALS AND METHODS: In the study, dose distributions obtained from the collapse cone algorithm, MC algorithm, and MC simulation were examined. The EGSnrc was utilized for MC simulation. RESULTS: In radiation fields smaller than 3 cm × 3 cm, the doses calculated by the CC algorithm are particularly high in the region of lung/soft-tissue interfaces. In the region of soft-tissue/vertebral interfaces, the doses calculated by the CC algorithm and the MC algorithm are compatible with the MC simulation. For each algorithm, the main reason for the non-overlapping dose curves in small fields compared to MC simulation is that the lateral electronic equilibrium loss is not taken into account by the algorithms. CONCLUSION: The doses calculated by the algorithms used in TPS may differ, especially in environments where density changes are sharp. Even if the radiation dose from different angles is calculated similarly in the target area by the algorithms, the calculated doses in the tissues in each radiation field path may be different. Therefore, to increase the quality of radiotherapy and to protect critical organs more accurately, the accuracy of the algorithms in TPS should be checked before treatment, especially in multi-field treatments such as stereotactic body radiation therapy and intensity-modulated radiotherapy for tumors in the abdominal region.

Asunto(s)

RESUMEN

Purpose. In this study, Monte Carlo (MC) simulations were done to relate the dose-response of the film to that in water. The effect of backscattering materials (PMMA, lead, polystyrene, and air) was investigated on its influence on film density for radionuclides including Am-241, Tc-99m, I-131, Cs-137.Methods. A BEAMnrc MC simulation was designed to score a phase-space file (PSF) below the container of the radionuclide under consideration to use as an input file for the subsequent DOSXYZnrc MC simulation. The geometry of the container holding the radionuclide was built using the component modules available in BEAMnrc. BEAMDP was used to investigate the container effect on the radionuclide spectrum as well as the fluence. The DOSXYZnrc simulation produced the absorbed dose in XR-QA2 and RT-QA2 GafchromicTMfilms. The DOSXYZnrc simulations were repeated for the GafchromicTMfilm now replaced with water to get the absorbed dose in water. From these results, conversion factors for the dose in water to the film dose for the different radionuclides, Am-241, Tc-99m, I-131, and Cs-137 were obtained. The actual film dose was calculated using the specific gamma exposure constant (Γ) at a distance of 50 cm for a point source approximation. From the BEAMnrc simulations, the particle fluence was extracted from PSFs to correct for the fluence at 0.1 cm below the sources from the fluence 50 cm away since the inverse square law will not apply to finite-size sources. The absorbed dose profiles in the film were compared to the absorbed dose profiles from the MC simulations.Results. A fitting function based on the neutron depletion model fits the optical density versus absorbed film dose data well and can be used as a calibration tool to obtain the film dose from its optical density. Lead as a backscatter material results in a higher optical density change but a lower absorbed dose. The XR-QA2 GafchromicTMfilm is more sensitive than the RT-QA2 GafchromicTMfilm, showing a more responsive optical density (OD) change in the energy range of radionuclides used in this study. Conversion factors were determined to convert the dose in water to the dose in GafchromicTMfilm. The Am-241 and I-131 simulated absorbed dose in the film to dose in water does not fluctuate as much as the simulated absorbed dose in film and water when using Tc-99m and Cs-137. Validation was shown for the comparison of the film and MC simulation absorbed dose profiles.Conclusions. MC BEAMnrc simulations are useful to simulate radionuclides and their containers. BEAMDP extracted energy spectra showed that the radionuclide containers produced a Compton effect on the energy spectra and added filtration on the lower spectral photon components. Extracted fluence ratios from PSFs were used to calculate the absorbed dose value at 0.1 cm distance from the source. By using the fit function, the dose in the film can be determined for known optical density values. The effect of the backscatter materials showed that the XR-QA2 GafchromicTMfilm results in higher optical density values than the RT-QA2 GafchromicTMfilm. The absorbed dose in both the films is comparable but not for a radionuclide such as Am-241 with an activity of 74MBq. The lead backscatter material showed to be the most prominent in optical density enhancement, and the air equivalent material was the least prominent. The XR-QA2 GafchromicTMfilm is the most sensitive and will be the best option if working with low energies. The absorbed dose in the XR-QA2 GafchromicTMfilm also showed a good comparison to the absorbed dose in water for the Am-241 radionuclide with an activity of 74MBq. The absorbed dose in the films compares well to the MC simulated doses.

Asunto(s)

RESUMEN

The commissioning and benchmark of a Monte Carlo (MC) model of the 6-MV Brainlab-Mitsubishi Vero4DRT linear accelerator for the purpose of quality assurance of clinical dynamic wave arc (DWA) treatment plans is reported. Open-source MC applications based on EGSnrc particle transport codes are used to simulate the medical linear accelerator head components. Complex radiotherapy irradiations can be simulated in a single MC run using a shared library format combined with BEAMnrc "source20." Electron energy tuning is achieved by comparing measured vs simulated percentage depth doses (PDDs) for MLC-defined field sizes in a water phantom. Electron spot size tuning is achieved by comparing measured and simulated inplane and crossplane beam profiles. DWA treatment plans generated from RayStation (RaySearch) treatment planning system (TPS) are simulated on voxelized (2.5 mm3 ) patient CT datasets. Planning target volume (PTV) and organs at risk (OAR) dose-volume histograms (DVHs) are compared to TPS-calculated doses for clinically deliverable dynamic volumetric modulated arc therapy (VMAT) trajectories. MC simulations with an electron beam energy of 5.9 MeV and spot size FWHM of 1.9 mm had the closest agreement with measurement. DWA beam deliveries simulated on patient CT datasets results in DVH agreement with TPS-calculated doses. PTV coverage agreed within 0.1% and OAR max doses (to 0.035 cc volume) agreed within 1 Gy. This MC model can be used as an independent dose calculation from the TPS and as a quality assurance tool for complex, dynamic radiotherapy treatment deliveries. Full patient CT treatment simulations are performed in a single Monte Carlo run in 23 min. Simulations are run in parallel using the Condor High-Throughput Computing software1 on a cluster of eight servers. Each server has two physical processors (Intel Xeon CPU E5-2650 0 @2.00 GHz), with 8 cores per CPU and two threads per core for 256 calculation nodes.

Asunto(s)

RESUMEN

OBJECTIVE: This study is to validate the utilization of Monte Carlo (MC) simulation to model the head of Primus linear accelerator, thereafter, using it to estimate the energy fluence distribution (EFD), the percentage depth dose (PDD), and beam profiles. MATERIALS AND METHODS: The BEAMNRC code that is based on the EGSNRC code has been used for modeling the linear accelerator head for 10 MeV electron beam with different applicator sizes (10 × 10, 15 × 15, and 20 × 20 cm2 ). The phase space was acquired from BEAMNRC at the end of each applicator and then used as an input file to DOSXYZNRC and BEAMDP to calculate the EFD, PDD, and beam profiles. RESULTS: There were a good consistency between the outcomes of the MC simulation and measured PDD and off-axis dose profiles that performed in a water phantom for all applicators. The PDD for the applicators proved to be favorable as a direct comparison of R100 , R90 , R80 , and R50 yielded results of < 2 mm, while it was 6 mm in R100 for the applicator 15 × 15 cm2 . The discrepancies in the surface doses (<3%) showed a quick decline in the build-up region and differences reached 0% within the first 2.4 mm. For the beam profiles comparison, the differences ranged from 2% (2 mm) to 3% (6 mm) for all applicators. CONCLUSION: Our examination demonstrated that the MC simulation by BEAMNRC code was accurate in modeling the Primus linear accelerator head.

Asunto(s)

RESUMEN

An empirical model for small circular electron fields was developed. This can be of great help in the treatment planning process for small circular electron fields. A complete dosimetric analysis of the circular fields defined by electron cutouts diameters (2 cm-9 cm) was done for nominal electron energies ranging between 6 MeV and 20 MeV using a 3D water phantom and a pin-point ion chamber. Properties studied included depth dose, in-air inverse-square fall-off, and beam profiles. The Varian Clianc 2100 C accelerator was modelled, benchmarked and Monte Carlo simulations were performed using the EGSnrc/BEAMnrc code for the small circular cutouts. A simple exponential model was found to accurately predict the very important therapeutic depth (90% of Dmax) for the small circular field size within an accuracy of better than 2 mm in most cases. The model has only two parameters (d90and 'b'). Also, the penumbra widths (90% of the off-axis profiles) of these small circular electron fields were studied and least square fitted to a simple quadratic model. Full dosimetric profiles of these small circular electron fields were further studied using the benchmarked Monte Carlo simulations. This study presents a simple model to predict the very important therapeutic depth (90% of Dmax) and a recipe to develop such an electron treatment model for any linear accelerator system. Such predictions can be extremely valuable and time saving prior to treatment planning involving not only small circular shaped electron fields but also irregularly shaped electron fields.

Asunto(s)

RESUMEN

PURPOSE: An EGSnrc based electron model was developed and validated for an Elekta Synergy® 160-leaf Agility™ linear accelerator. It was able to reproduce measured central axis (CAX) percentage depth dose (PDD) curves and off-axis profiles (OAPs) within 2%/2 mm, and relative output factors (ROFs) within 3%. METHODS: BEAMnrc component modules were used to model the accelerator accurately according to vendor supplied specifications. The electron beam focal spot size and input energy spectrum were determined through their effects on electron CAX PDDs and OAPs as benchmarked against water tank data. Phase space files were used as source input in DOSXYZnrc water phantom simulations. Dose distributions were calculated for six electron nominal energies, 11 field sizes and two source-to-surface distances. RESULTS: The full width at half maximum of the focal spot (assuming a Gaussian intensity distribution) was determined to be 1.50 mm. An asymmetrical input electron energy spectrum with a low-energy tail produced good agreement with measured data and solved the match in the build-up (BU) region for all electron energies used in this study. CONCLUSIONS: The improved input electron spectra for the electron model could predict dose distributions within 2%/2 mm of measured data. The model's success is embedded in the asymmetrical energy spectrum which provided a valuable free parameter which, by fine adjustment, improved the match in the BU region of dose distributions. Furthermore, focal spot parameters could be determined by means of iterative simulations.

Asunto(s)

RESUMEN

AIM: In this study, we investigated initial electron parameters of Siemens Artiste Linac with 6 MV photon beam using the Monte Carlo method. BACKGROUND: It is essential to define all the characteristics of initial electrons hitting the target, i.e. mean energy and full width of half maximum (FWHM) of the spatial distribution intensity, which is needed to run Monte Carlo simulations. The Monte Carlo is the most accurate method for simulation of radiotherapy treatments. MATERIALS AND METHODS: Linac head geometry was modeled using the BEAMnrc code. The phase space files were used as input file to DOSXYZnrc simulation to determine the dose distribution in a water phantom. We obtained percent depth dose curves and the lateral dose profile. All the results were obtained at 100 cm of SSD and for a 10 × 10 cm2 field. RESULTS: We concluded that there existed a good conformity between Monte Carlo simulation and measurement data when we used electron mean energy of 6.3 MeV and 0.30 cm FWHM value as initial parameters. We observed that FWHM values had very little effect on PDD and we found that the electron mean energy and FWHM values affected the lateral dose profile. However, these effects are between tolerance values. CONCLUSIONS: The initial parameters especially depend on components of a linac head. The phase space file which was obtained from Monte Carlo Simulation for a linac can be used as calculation of scattering, MLC leakage, to compare dose distribution on patients and in various studies.

RESUMEN

PURPOSE: A Monte Carlo model of a 6 MV medical linear accelerator (linac) unit built indigenously was developed using the BEAMnrc user code of the EGSnrc code system. The model was benchmarked against the measurements. Monte Carlo simulations were carried out for different incident electron beam parameters in the study. MATERIALS AND METHODS: Simulation of indigenously developed linac unit has been carried out using the Monte Carlo based BEAMnrc user-code of the EGSnrc code system. Using the model, percentage depth dose (PDD), and lateral dose profiles were studied using the DOSXYZnrc user code. To identify appropriate electron parameters, three different distributions of electron beam intensity were investigated. For each case, the kinetic energy of the incident electron was varied from 6 to 6.5 MeV (0.1 MeV increment). The calculated dose data were compared against the measurements using the PTW, Germany make RFA dosimetric system (water tank MP3-M and 0.125 cm3 ion chamber). RESULTS: The best fit of incident electron beam parameter was found for the combination of beam energy of 6.2 MeV and circular Gaussian distributed source in X and Y with FWHM of 1.0 mm. PDD and beam profiles (along both X and Y directions) were calculated for the field sizes from 5 cm × 5 cm to 25 cm × 25 cm. The dose difference between the calculated and measured PDD and profile values were under 1%, except for the penumbra region where the maximum deviation was found to be around 2%. CONCLUSIONS: A Monte Carlo model of indigenous linac (6 MV) has been developed and benchmarked against the measured data.

RESUMEN

OBJECTIVES: To evaluate the potential privileges of flattening filter-free (FFF) photon beams from Oncor® linac for 6 MV and 18 MV energies. METHODS: A Monte Carlo (MC) model of Oncor® linac was built using BEAMnrc MCCode and verified by the measured data using 6 MV and 18 MV energies. A comprehensive set of data was also characterized for MC model of Oncor® machine running with and without flattening filter (FF) for 6 MV and 18 MV beams in six field sizes. The investigated characteristics included mean energy, energy spectrum, photon spatial fluence, superficial dose, percent depth dose (PDD), dose output, and out-of-field dose with two indexes of lateral dose profile and isodose curve at three depths. RESULTS: Using FFF enhanced the energy uniformity 3.4±0.11% (6 MV) and 2.05±0.09% (18 MV) times and improved dose output by factor of 2.91 (6 MV) and 4.2 (18 MV) on the central axis, respectively. Using FFF also reduced the PDD dependencies by 9.1% (6 MV) and 5.57% (18 MV). In addition, using FFF had a lower out-of-field dose due to the reduced head scatter and softer spectra. CONCLUSIONS: The findings in this study suggested that using FFF, Oncor® machine could achieve better treatment results with lower dose toxicity and a shorter beam-on time.

Asunto(s)

RESUMEN

We have developed a new component module for the BEAMnrc software package, called SMLC, which models the tongue-and-groove structure of the Siemens Optifocus multileaf collimator. The ultimate goal is to perform accurate Monte Carlo simulations of the IMRT treatments carried out with Optifocus. SMLC has been validated by direct geometry checks and by comparing quantitatively the results of simulations performed with it and with the component module VARMLC. Measurements and Monte Carlo simulations of absorbed dose distributions of radiation fields sensitive to the tongue-and-groove effect have been performed to tune the free parameters of SMLC. The measurements cannot be accurately reproduced with VARMLC. Finally, simulations of a typical IMRT field showed that SMLC improves the agreement with experimental measurements with respect to VARMLC in clinically relevant cases. PACS NUMBER: 87.55. K.

Asunto(s)

RESUMEN

An automated Multi-Leaf Collimator (MLC) system has been developed as add-on for the cobalt-60 teletherapy machines available in India. The goal of the present computational study is to validate the MLC design using Monte Carlo (MC) modeling. The study was based on the Kirloskar-supplied Phoenix model machines that closely match the Atomic Energy of Canada Limited (AECL) theratron-80 machine. The MLC is a retrofit attachment to the collimator assembly, with 14 non-divergent leaf pairs of 40 mm thick, 7 mm wide, and 150 mm long tungsten alloy plates with rounded edges and 20 mm tongue and 2 mm groove in each leaf. In the present work, the source and collimator geometry has been investigated in detail to arrive at a model that best represents the measured dosimetric data. The authors have studied in detail the proto-I MLC built for cobalt-60. The MLC field sizes were MC simulated for 2 × 2 cm(2) to 14 × 14 cm(2) square fields as well as irregular fields, and the percent depth dose (PDD) and profile data were compared with ROPS(†) treatment planning system (TPS). In addition, measured profiles using the IMATRIXX system(‡) were also compared with the MC simulations. The proto-I MLC can define radiation fields up to 14 × 14 cm(2) within 3 mm accuracy. The maximum measured leakage through the leaf ends in closed condition was 3.4% and interleaf leakage observed was 7.3%. Good agreement between MC results, ROPS and IMATRIXX results has been observed. The investigation also supports the hypothesis that optical and radiation field coincidence exists for the square fields studied with the MLC. Plots of the percent depth dose (PDD) data and profile data for clinically significant irregular fields have also been presented. The MC model was also investigated to speed up the calculations to allow calculations of clinically relevant conformal beams. (†)Radiation Oncology Planning System (ROPS) is supplied by Tirumala Jyothi Computer Systems described at https://sites.google.com/site/tjcsrops/ (‡)IMATRIXX is supplied by IBA Dosimetry described at HYPERLINK http://www.iba-dosimetry.com.

RESUMEN

Dental prostheses made of high density material contribute to modify dose distribution in head and neck cancer treatment. Our objective is to quantify dose perturbation due to high density inhomogeneity with experimental measurements and Monte Carlo simulations. Firstly, measurements were carried in a phantom representing a human jaw with thermoluminescent detectors (GR200A) and EBT2 Gafchromic films in the vicinity of three samples: a healthy tooth, a tooth with amalgam and a Ni-Cr crown, irradiated in clinical configuration. Secondly, Monte Carlo simulations (BEAMnrc code) were assessed in an identical configuration. Experimental measurements and simulation results confirm the two well-known phenomena: firstly the passage from a low density medium to a high density medium induces backscattered electrons causing a dose increase at the interface, and secondly, the passage from a high density medium to a low density medium creates a dose decrease near the interface. So, the results show a 1.4% and 23.8% backscatter dose rise and attenuation after sample of 26.7% and 10.9% respectively for tooth with amalgam and crown compared to the healthy tooth. Although a tooth with amalgam has a density of about 12-13, the changes generated are not significant. However, the results for crown (density of 8) are very significant and the discordance observed may be due to calculation point size difference 0.8 mm and 0.25 mm respectively for TLD and Monte Carlo. The use of Monte Carlo simulations and experimental measurements provides objective evidence to evaluate treatment planning system results with metal dental prostheses.

Asunto(s)

RESUMEN

The Monte Carlo method is the most accurate method for simulation of radiation therapy equipment. The linear accelerators (linac) are currently the most widely used machines in radiation therapy centers. In this work, a Monte Carlo modeling of the Siemens ONCOR linear accelerator in 6 MV and 18 MV beams was performed. The results of simulation were validated by measurements in water by ionization chamber and extended dose range (EDR2) film in solid water. The linac's X-ray particular are so sensitive to the properties of primary electron beam. Square field size of 10 cm × 10 cm produced by the jaws was compared with ionization chamber and film measurements. Head simulation was performed with BEAMnrc and dose calculation with DOSXYZnrc for film measurements and 3ddose file produced by DOSXYZnrc analyzed used homemade MATLAB program. At 6 MV, the agreement between dose calculated by Monte Carlo modeling and direct measurement was obtained to the least restrictive of 1%, even in the build-up region. At 18 MV, the agreement was obtained 1%, except for in the build-up region. In the build-up region, the difference was 1% at 6 MV and 2% at 18 MV. The mean difference between measurements and Monte Carlo simulation is very small in both of ONCOR X-ray energy. The results are highly accurate and can be used for many applications such as patient dose calculation in treatment planning and in studies that model this linac with small field size like intensity-modulated radiation therapy technique.

RESUMEN

This investigation aims to design a practical multi-leaf collimator (MLC) system for the cobalt teletherapy machine and check its radiation properties using the Monte Carlo (MC) method. The cobalt machine was modeled using the BEAMnrc Omega-Beam MC system, which could be freely downloaded from the website of the National Research Council (NRC), Canada. Comparison with standard depth dose data tables and the theoretically modeled beam showed good agreement within 2%. An MLC design with low melting point alloy (LMPA) was tested for leakage properties of leaves. The LMPA leaves with a width of 7 mm and height of 6 cm, with tongue and groove of size 2 mm wide by 4 cm height, produced only 4% extra leakage compared to 10 cm height tungsten leaves. With finite (60)Co source size, the interleaf leakage was insignificant. This analysis helped to design a prototype MLC as an accessory mount on a cobalt machine. The complete details of the simulation process and analysis of results are discussed.