Determining optimization of the initial parameters in Monte Carlo simulation for linear accelerator radiotherapy

Monte Carlo Model Validation of 6MV Beam of OMID, the First Iranian Linear Accelerator

Monte Carlo (MC) techniques are regarded as an accurate method to simulate the dose calculation in radiotherapy for many years. The present paper aims to validate the simulated model of the 6-MV beam of OMID linear accelerator (BEHYAAR Company) by EGSnrc codes system and also investigate the effects of initial electron beam parameters (energy, radial full width at half maximum, and mean angular spread) on dose distributions. For this purpose, the comparison between the calculated and measured percentage depth dose (PDD) and lateral dose profiles was done by gamma index (GI) with 1%-1 mm acceptance criteria. MC model validating was done for 3 cm × 3 cm, 5 cm × 5 cm, 8 cm × 8 cm, 10 cm × 10 cm, and 20 cm × 20 cm field sizes. To study the sensitivity of model to beam parameters, the field size was selected as 10 cm × 10 cm and 30 cm × 30 cm. All lateral dose profiles were obtained at 10 cm. Excellent agreement was achieved with a 99.2% GI passing percentage for PDD curves and at least 93.8% GI for lateral dose profiles for investigated field sizes. Our investigation confirmed that the lateral dose profile severely depends on the considered source parameters in this study. PDD only considerably depends on the initial electron beam energy. Therefore, source parameters should not be specified independently. These results indicate that the current model of OMID 6-MV Linac is well established, and the accuracy of the simulation is high enough to be used in various applications.

Comparison of Monaco treatment planning system algorithms and Monte Carlo simulation for small fields in anthropomorphic RANDO phantom: The esophagus case

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.

Simplified sigmoidal curve fitting for a 6 MV FFF photon beam of the Halcyon to determine the field size for beam commissioning and quality assurance

Background

An O-ring gantry-type linear accelerator (LINAC) with a 6-MV flattening filter-free (FFF) photon beam, Halcyon, includes a reference beam that contains representative information such as the percent depth dose, profile and output factor for commissioning and quality assurance. However, because it does not provide information about the field size, we proposed a method to determine all field sizes according to all depths for radiation therapy using simplified sigmoidal curve fitting (SCF).

Methods

After mathematical definition of the SCF using four coefficients, the defined curves were fitted to both the reference data (RD) and the measured data (MD). For good agreement between the fitting curve and the profiles in each data set, the field sizes were determined by identifying the maximum point along the third derivative of the fitting curve. The curve fitting included the field sizes for beam profiles of 2 × 2, 4 × 4, 6 × 6, 8 × 8, 10 × 10, 20 × 20 and 28 × 28 cm² as a function of depth (at 1.3, 5, 10 and 20 cm). The field size results from the RD were compared with the results from the MD using the same condition.

Results

All fitting curves show goodness of fit, R², values that are greater than 0.99. The differences in field size between the RD and the MD were within the range of 0 to 0.2 cm. The smallest difference in the field sizes at a depth of 10 cm, which is a surface-to-axis distance, was reported.

Conclusion

Application of the SCF method has been proven to accurately capture the field size of the preconfigured RD and the measured FFF photon beam data for the Halcyon system. The current work can be useful for beam commissioning as a countercheck methodology to determine the field size from RD in the treatment planning system of a newly installed Halcyon system and for routine quality assurance to ascertain the correctness of field sizes for clinical use of the Halcyon system.

Monte Carlo commissioning of radiotherapy LINAC-Introducing an improved methodology

Purpose

Monte Carlo (MC) commissioning of medical linear accelerator (LINAC) is a time-consuming process involving a comparison between measured and simulated cross beam/lateral profiles and percentage depth doses (PDDs) for various field sizes. An agreement between these two data sets is sought by trial and error method while varying the incident electron beam parameters, such as electron beam energy or width, etc. This study aims to improve the efficiency of MC commissioning of a LINAC by assessing the feasibility of using a limited number of simulated PDDs.

Materials and methods

Using EGSnrc codes, a Varian Clinac 2100 unit has been commissioned for 6 MV photon beam, and a methodology has been proposed to identify the incident electron beam parameters in a speedier fashion. Impact of voxel size in 3-dimensions and cost functions used for comparison of the measured and simulated data have been investigated along with the role of interpolation.

Results

A voxel size of 1 × 1×0.5 cm³ has been identified as suitable for accurate and fast commissioning of the LIANC. The optimum number of simulated PDDs (required for further interpolation) has been found to be five.

Conclusion

The present study suggests that PDDs alone at times can be insufficient for an unambiguous commissioning process and should be supported by including the lateral beam profiles in the process.

Determination of initial electron parameters by means of Monte Carlo simulations for the Siemens Artiste Linac 6 MV photon beam

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 cm² 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.

Benchmarking of Monte Carlo model of Siemens Oncor® linear accelerator for 18MV photon beam: Determination of initial electron beam parameters

OBJECTIVE

This study aims to benchmark a Monte Carlo (MC) model of the 18 MV photon beam produced by the Siemens Oncor® linac using the BEAMnrc and DOSXYZnrc codes.

METHODS

By matching the percentage depth doses and beam profiles calculated by MC simulations with measurements, the initial electron beam parameters including electron energy, full width at half maximum (spatial FWHM), and mean angular spread were derived for the 10×10 cm2 and 20×20 cm2 field sizes. The MC model of the 18 MV photon beam was then validated against the measurements for different field sizes (5×5, 30×30 and 40×40 cm2) by gamma index analysis.

RESULTS

The optimum values for electron energy, spatial FWHM and mean angular spread were 14.2 MeV, 0.08 cm and 0.8 degree, respectively. The MC simulations yielded the comparable measurement results of these optimum parameters. The gamma passing rates (with acceptance criteria of 1% /1 mm) for percentage depth doses were found to be 100% for all field sizes. For cross-line profiles, the gamma passing rates were 100%, 97%, 95%, 96% and 95% for 5×5, 10×10, 20×20, 30×30 and 40×40 cm2 field sizes, respectively.

CONCLUSIONS

By validation of the MC model of Siemens Oncor® linac using various field sizes, it was found that both dose profiles of small and large field sizes were very sensitive to the changes in spatial FWHM and mean angular spread of the primary electron beam from the bending magnet. Hence, it is recommended that both small and large field sizes of the 18 MV photon beams should be considered in the Monte Carlo linac modeling.

Monte Carlo simulation of linac irradiation with dynamic wedges

This study aims to simulate the dose distributions of LINAC with dynamic wedges (DWs) under various field sizes and wedge angles by the BEAMnrc code with DYNJAWS component module. These were compared with those calculated by the treatment planning system (TPS) and the measured data. All percentage depth doses (PDDs) were found to be in good agreement between TPS, Monte Carlo (MC) and measurements made in open fields and fields with DWs. For dose profiles, compared with the MC and the measurements, TPS gives reliable results for large field sizes (>10 × 10 cm(2)) but results in significant errors in small field sizes (5 × 5 cm(2)). The entrance surface doses calculated by TPS were found to be significantly overestimated. For depths deeper than 0.5 cm, TPS yields PDDs in agreement with MC simulations.