Commissioning of a medical accelerator photon beam Monte Carlo simulation using wide-field profiles

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

PURPOSE/OBJECTIVE

The 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/METHODS

Using 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.

RESULTS

Validation 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%.

CONCLUSION

SBRT treatment plans of Monaco agreed with MC results. It possible to use MC for treatment plans verifications as independent QC tool.

Spatial Mesh-Based Surface Source Model for the Electron Contamination of an 18 MV Photon Beams

Background

Source modeling is an approach to reduce computational burden in Monte Carlo simulations but at the cost of reduced accuracy. Although this method can be effective, one component of the source model that is exceptionally difficult to model is the electron contamination, a significant contributor to the skin and shallow dose.

Aims and Objectives

To improve the accuracy for the electron contamination component of the overall source model, we have generated a spatial mesh based surface source model.

Methods and Materials

The source model is located downstream from the flattening filter and mirror but upstream from the movable jaws. A typical phase space file uses around ten parameters per particle, but this method simplifies this number to five components. By using only the electron distance from the central axis, angles from the central axis and energy, the computational time and disk space required is greatly reduced.

Results and Conclusion

Despite the simplification in the source model, the electron contamination is still accurate to within 1.5%.

The comparison between 6 MV Primus LINAC simulation output using EGSnrc and commissioning data

Introduction

Monte Carlo calculation method is considered to be the most accurate method for dose calculation in radiotherapy. The purpose of this research is comparison between 6 MV Primus LINAC simulation output with commissioning data using EGSnrc and build a Monte Carlo geometry of 6 MV Primus LINAC as realistically as possible. The BEAMnrc and DOSXYZnrc (EGSnrc package) Monte Carlo model of the LINAC head was used as a benchmark.

Methods

In the first part, the BEAMnrc was used for the designing of the LINAC treatment head. In the second part, dose calculation and for the design of 3D dose file were produced by DOSXYZnrc. The simulated PDD and beam profile obtained were compared with that calculated using commissioning data. Good agreement was found between calculated PDD (1·1%) and beam profile using Monte Carlo simulation and commissioning data. After validation, TPR20,10, TMR and Spvalues were calculated in five different field.

Results

Good agreement was found between calculated values by using Monte Carlo simulation and commissioning data. Average differences for five field sizes in this approach is about 0·83% for Sp. for TPR20,10differences for field sizes 10×10 cm2is 0·29% and for TMR in five field sizes, the average value is ~1·6%.

Conclusion

In conclusion, the BEAMnrc and DOSXYZnrc codes package have very good accuracy in calculating dose distribution for 6 MV photon beam and it can be considered as a promising method for patient dose calculations and also the Monte Carlo model of primus linear accelerator built in this study can be used as method to calculate the dose distribution for cancer patients.

An investigation into the use of MMCTP to tune accelerator source parameters and testing its clinical application

This paper presents an alternative method to tune Monte Carlo electron beam parameters to match measured data using a minimal set of variables in order to reduce the model setup time prior to clinical implementation of the model. Monte Carlo calculations provide the possibility of a powerful treatment planning verification technique. The nonstandardized and nonautomated process of tuning the required accelerator model is one of the reasons for delays in the clinical implementation of Monte Carlo techniques. This work aims to establish and verify an alternative tuning method that can be carried out in a minimal amount of time, allowing it to be easily implemented in a clinical setting by personnel with minimal experience with Monte Carlo methods. This tuned model can then be incorporated into the MMCTP system to allow the system to be used as a second dose calculation check for IMRT plans. The technique proposed was used to establish the primary electron beam parameters for accelerator models for the Varian Clinac 2100 6 MV photon beam using the BEAMnrc Monte Carlo system. The method is intended to provide a clear, direct, and efficient process for tuning an accelerator model using readily available clinical quality assurance data. The tuning provides a refined model, which agrees with measured dose profile curves within 1.5% outside the penumbra or 3 mm in the penumbra, for square fields with sides of 3 cm up to 30 cm. These models can then be employed as the basis for Monte Carlo recalculations of dose distributions, using the MMCTP system, for clinical treatment plans, providing an invaluable assessment tool. This was tested on six IMRT plans and compared to the measurements performed for the pretreatment QA process. These Monte Carlo values for the average dose to the chamber volume agreed with measurements to within 0.6%.

PACS number: 87.55.km

Secondary neutron spectra from modern Varian, Siemens, and Elekta linacs with multileaf collimators

Neutrons are a by-product of high-energy x-ray radiation therapy (threshold for [gamma,n] reactions in high-Z material -7 MeV). Neutron production varies depending on photon beam energy as well as on the manufacturer of the accelerator. Neutron production from modern linear accelerators (linacs) has not been extensively compared, particularly in terms of the differences in the strategies that various manufacturers have used to implement multileaf collimators (MLCs) into their linac designs. However, such information is necessary to determine neutron dose equivalents for different linacs and to calculate vault shielding requirements. The purpose of the current study, therefore, was to measure the neutron spectra from the most up-to-date linacs from three manufacturers: Varian 21EX operating at 15, 18, and 20 MV, Siemens ONCOR operating at 15 and 18 MV, and Elekta Precise operating at 15 and 18 MV. Neutron production was measured by means of gold foil activation in Bonner spheres. Based on the measurements, the authors determined neutron spectra and calculated the average energy, total neutron fluence, ambient dose equivalent, and neutron source strength. The shapes of the neutron spectra did not change significantly between accelerators or even as a function of treatment energy. However, the neutron fluence, and therefore the ambient dose equivalent, did vary, increasing with increasing treatment energy. For a given nominal treatment energy, these values were always highest for the Varian linac. The current study thus offers medical physicists extensive information about the neutron production of MLC-equipped linacs currently in operation and provides them information vital for accurate comparison and prediction of neutron dose equivalents and calculation of vault shielding requirements.

Automatic determination of primary electron beam parameters in Monte Carlo simulation

In order to obtain realistic and reliable Monte Carlo simulations of medical linac photon beams, an accurate determination of the parameters that define the primary electron beam that hits the target is a fundamental step. In this work we propose a new methodology to commission photon beams in Monte Carlo simulations that ensures the reproducibility of a wide range of clinically useful fields. For such purpose accelerated Monte Carlo simulations of 2 x 2, 10 x 10, and 20 x 20 cm2 fields at SSD = 100 cm are carried out for several combinations of the primary electron beam mean energy and radial FWHM. Then, by performing a simultaneous comparison with the correspondent measurements for these same fields, the best combination is selected. This methodology has been employed to determine the characteristics of the primary electron beams that best reproduce a Siemens PRIMUS and a Varian 2100 CD machine in the Monte Carlo simulations. Excellent agreements were obtained between simulations and measurements for a wide range of field sizes. Because precalculated profiles are stored in databases, the whole commissioning process can be fully automated, avoiding manual fine-tunings. These databases can also be used to characterize any accelerators of the same model from different sites.

Estimation of the focal spot size and shape for a medical linear accelerator by Monte Carlo simulation

The focal spot size and shape of a medical linac are important parameters that determine the dose profiles, especially in the penumbral region. A relationship between the focal spot size and the dose profile penumbra has been studied and established from simulation results of the EGSnrc Monte Carlo code. A simple method is proposed to estimate the size and the shape of a linac's focal spot from the measured dose profile data.

Monte Carlo study of Siemens PRIMUS photoneutron production

Neutron production in radiotherapy facilities has been studied from the early days of modern linacs. Detailed studies are now possible using photoneutron capabilities of general-purpose Monte Carlo codes at energies of interest in medical physics. The present work studies the effects of modelling different accelerator head and room geometries on the neutron fluence and spectra predicted via Monte Carlo. The results from the simulation of a 15 MV Siemens PRIMUS linac show an 80% increase in the fluence scored at the isocentre when, besides modelling the components necessary for electron/photon simulations, other massive accelerator head components are included. Neutron fluence dependence on inner treatment room volume is analysed showing that thermal neutrons have a 'gaseous' behaviour and then a 1/V dependence. Neutron fluence maps for three energy ranges, fast (E > 0.1 MeV), epithermal (1 eV < E < 0.1 MeV) and thermal (E < 1 eV), are also presented and the influence of the head components on them is discussed.