Inference of the optimal pretarget electron beam parameters in a Monte Carlo virtual linac model through simulated annealing

Reconstruction of the electron source intensity distribution of a clinical linear accelerator using in-air measurements and a genetic algorithm

Background and purpose

The electron source intensity distribution of a clinical linear accelerator has a great influence on the calculation of output factors for small radiation fields where source occlusion by the collimating devices takes place. The purpose of this study was to present a new method for the electron source reconstruction problem.

Materials and methods

The measurements were performed in-air using diode and 6 MV 1 × 1 cm² photon field in flattening filter-free mode. In Monte Carlo simulation, an electron target area was divided into a number of square subsources. Then, the in-air doses in 2D silicon chip array were calculated individually from each subsource. A genetic algorithm search was applied in order to determine the optimal weight factors for all subsources that provide the best agreement between simulated and measured doses.

Results

It was found that the reconstructed electron source intensity from a clinical linear accelerator has the two-dimensional elliptical double Gaussian distribution. The source intensity distribution consisted of two intensity components along the in-plane (x) and cross-plane (y) directions characterized by full width half-maximum (FWHM): FWHM_x1 = 0.27 cm, FWHM_x2 = 0.08 cm, FWHM_y1 = 0.24 cm, FWHM_y2 = 0.06 cm, where broader components are 81% and 53% of the total intensity along × and y axis respectively.

Conclusions

The obtained results demonstrated an elliptical double Gaussian intensity distribution of the incident electron source. We anticipate that the proposed method has universal applications independent of the type of linear accelerator, modality or energy.

Monte Carlo and analytic modeling of an Elekta Infinity linac with Agility MLC: Investigating the significance of accurate model parameters for small radiation fields

PURPOSE

To explain the deviation observed between measured and Monaco calculated dose profiles for a small field (i.e., alternating open-closed MLC pattern). A Monte Carlo (MC) model of an Elekta Infinity linac with Agility MLC was created and validated against measurements. In addition, an analytic model which predicts the fluence at the isocenter plane was used to study the impact of multiple beam parameters on the accuracy of dose calculations for small fields.

METHODS

A detailed MC model of a 6 MV Elekta Infinity linac with Agility MLC was created in EGSnrc/BEAMnrc and validated against measurements. An analytic model using primary and secondary virtual photon sources was created and benchmarked against the MC simulations and the impact of multiple beam parameters on the accuracy of the model for a small field was investigated. Both models were used to explain discrepancies observed between measured/EGSnrc simulated and Monaco calculated dose profiles for alternating open-closed MLC leaves.

RESULTS

MC-simulated dose profiles (PDDs, cross- and in-line profiles, etc.) were found to be in very good agreements with measurements. The best fit for the leaf bank rotation was found to be 9 mrad to model the defocusing of Agility MLC. Moreover, a very good agreement was observed between results from the analytic model and MC simulations for a small field. Modifying the radial size of the incident electron beam in the BEAMnrc model improved the agreement between Monaco and EGSnrc calculated dose profiles by approximately 16% and 30% in the position of maxima and minima, respectively.

CONCLUSION

Accurate modeling of the full-width-half-maximum (FWHM) of the primary photon source as well as the MLC leaf design (leaf bank rotation, etc.) is essential for accurate calculations of dose delivered by small radiation fields when using virtual source or MC models of the beam.

Monte Carlo Investigation of Photon Beam Characteristics and its Variation with Incident Electron Beam Parameters for Indigenous Medical Linear Accelerator

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

Comparative evaluation of modern dosimetry techniques near low- and high-density heterogeneities

The purpose of this study is to compare performance of several dosimetric methods in heterogeneous phantoms irradiated by 6 and 18 MV beams. Monte Carlo (MC) calculations were used, along with two versions of Acuros XB, anisotropic analytical algorithm (AAA), EBT2 film, and MOSkin dosimeters. Percent depth doses (PDD) were calculated and measured in three heterogeneous phantoms. The first two phantoms were a 30×30×30 cm3 solid‐water slab that had an air‐gap of 20×2.5×2.35 cm3. The third phantom consisted of 30×30×5 cm3 solid water slabs, two 30×30×5 cm3 slabs of lung, and one 30×30×1 cm3 solid water slab. Acuros XB, AAA, and MC calculations were within 1% in the regions with particle equilibrium. At media interfaces and buildup regions, differences between Acuros XB and MC were in the range of +4.4% to −12.8%. MOSkin and EBT2 measurements agreed to MC calculations within ∼2.5%, except for the first centimeter of buildup where differences of 4.5% were observed. AAA did not predict the backscatter dose from the high‐density heterogeneity. For the third, multilayer lung phantom, 6 MV beam PDDs calculated by all TPS algorithms were within 2% of MC. 18 MV PDDs calculated by two versions of Acuros XB and AAA differed from MC by up to 2.8%, 3.2%, and 6.8%, respectively. MOSkin and EBT2 each differed from MC by up to 2.9% and 2.5% for the 6 MV, and by −3.1% and ∼2% for the 18 MV beams. All dosimetric techniques, except AAA, agreed within 3% in the regions with particle equilibrium. Differences between the dosimetric techniques were larger for the 18 MV than the 6 MV beam. MOSkin and EBT2 measurements were in a better agreement with MC than Acuros XB calculations at the interfaces, and they were in a better agreement to each other than to MC. The latter is due to their thinner detection layers compared to MC voxel sizes.

PACS numbers: 87.55.K‐, 87.55.kd, 87.55.km, 87.53.Bn, 87.55.k

Monte Carlo investigation into feasibility and dosimetry of flat flattening filter free beams

Flattening filter free (FFF) beams due to their non-uniformity, are sub-optimal for larger field sizes. The purpose of this study was to investigate the incident electron beam distributions that would produce flat FFF (F4) beams without the use of a flattening filter (FF). Monte Carlo (MC) simulations with BEAMnrc and DOSXYZnrc codes have been performed to evaluate feasibility of this approach. The dose distributions in water for open 6 MV beams were simulated using the Varian 21EX linac head model, which will be called the FF model. The FF was then removed from the FF model, and MC simulations were performed using (1) 6 MeV electrons incident on the target and (2) a 6 MeV electron beam with electron angular distributions optimized to provide as flat dose profiles as possible. Configuration (1) represents FFF beam while configuration (2) allowed producing a F4 beam. Optimizations have also been performed to produce flattest profiles for a set of dose rates (DRs) in the range from 1.25 to 2.4 of the DR of FF beam. Profiles and percentage depth doses (PDDs) from 6 MV F4 beams have been calculated and compared to those from the FF beam. Calculated profiles demonstrated improved flatness of the FFF beams. In fact, up to field sizes within the circle of 35 cm diameter the flatness of F4 beam at dmax was better or comparable to that of FF beam. At 20 cm off-axis the dose increased from 52% for FFF to 92% for F4 beam. Also, profiles of F4 beams did not change considerably with depth. PDDs from F4 beams were similar to those of the FFF beam. The DR for the largest modeled (44 cm diameter) F4 beam was higher than the DR from FF beam by a factor of 1.25. It was shown that the DR can be increased while maintaining beam flatness, but at the cost of reduced field size.

On Monte Carlo modeling of megavoltage photon beams: a revisited study on the sensitivity of beam parameters

PURPOSE

To commission Monte Carlo beam models for five Varian megavoltage photon beams (4, 6, 10, 15, and 18 MV). The goal is to closely match measured dose distributions in water for a wide range of field sizes (from 2 x 2 to 35 x 35 cm2). The second objective is to reinvestigate the sensitivity of the calculated dose distributions to variations in the primary electron beam parameters.

METHODS

The GEPTS Monte Carlo code is used for photon beam simulations and dose calculations. The linear accelerator geometric models are based on (i) manufacturer specifications, (ii) corrections made by Chibani and Ma ["On the discrepancies between Monte Carlo dose calculations and measurements for the 18 MV Varian photon beam," Med. Phys. 34, 1206-1216 (2007)], and (iii) more recent drawings. Measurements were performed using pinpoint and Farmer ionization chambers, depending on the field size. Phase space calculations for small fields were performed with and without angle-based photon splitting. In addition to the three commonly used primary electron beam parameters (E(AV) is the mean energy, FWHM is the energy spectrum broadening, and R is the beam radius), the angular divergence (theta) of primary electrons is also considered.

RESULTS

The calculated and measured dose distributions agreed to within 1% local difference at any depth beyond 1 cm for different energies and for field sizes varying from 2 x 2 to 35 x 35 cm2. In the penumbra regions, the distance to agreement is better than 0.5 mm, except for 15 MV (0.4-1 mm). The measured and calculated output factors agreed to within 1.2%. The 6, 10, and 18 MV beam models use theta = 0 degrees, while the 4 and 15 MV beam models require theta = 0.5 degrees and 0.6 degrees, respectively. The parameter sensitivity study shows that varying the beam parameters around the solution can lead to 5% differences with measurements for small (e.g., 2 x 2 cm2) and large (e.g., 35 x 35 cm2) fields, while a perfect agreement is maintained for the 10 x 10 cm2 field. The influence of R on the central-axis depth dose and the strong influence of theta on the lateral dose profiles are demonstrated.

CONCLUSIONS

Dose distributions for very small and very large fields were proved to be more sensitive to variations in E(AV), R, and theta in comparison with the 10 x 10 cm2 field. Monte Carlo beam models need to be validated for a wide range of field sizes including small field sizes (e.g., 2 x 2 cm2).