Search

Purpose: To develop, implement, and validate a full 1.5 T/7 MV magnetic resonance (MR)-Linac accelerator head and cryostat model in EGSnrc for high precision dose calculations accounting for magnetic field effects that are independent from the vendor treatment planning system., Methods: Primary electron beam parameters for the implemented model were adapted to be in accordance with measured dose profiles of the Elekta Unity (Elekta AB, Stockholm, Sweden). Parameters to be investigated were the mean electron energy as well as the Gaussian radial intensity and energy distributions. Energy tuning was done comparing depth dose profiles simulated with monoenergetic beams of varying energies to measurements. The optimum radial intensity distribution was found by varying the radial full width at half maximum (FWHM) and comparing simulated and measured lateral profiles. The influence of the energy distribution was investigated by comparing simulated lateral and depth dose profiles with varying energy spreads to measured data. Comparison of simulations and measurements was performed by calculating average and maximum local dose deviations. The model was validated recalculating a clinical intensity-modulated radiation therapy plan for the MR-Linac and comparing the resulting dose distribution with simulations from the commercial treatment planning system Monaco using the gamma criterion., Results: Comparison of simulated and measured data showed that the optimum initial electron beam for MR-Linac simulations was monoenergetic with an electron energy of (7.4 ± 0.2) MeV. The optimum Gaussian radial intensity distribution has a FWHM of (2.2 ± 0.3) mm. The average relative deviations were smaller than 1% for all simulated profiles with optimum electron parameters, whereas the largest maximum deviation of 2.07% was found for the 22 × 22 cm 2 cross-plane profile. Profiles were insensitive to energy spread variations. The IMRT plan recalculated with the final MR-Linac model with optimized initial electron beam parameters showed a gamma pass rate of 99.83 % using a gamma criterion of 3%/3 mm., Conclusions: The EGSnrc MR-Linac model developed in this study showed good accordance with measurements and was successfully used to recalculate a first full clinical IMRT treatment plan. Thus, it shows the general possibility for future secondary dose calculations of full IMRT plans with EGSnrc, which needs further detailed investigations before clinical use., (© 2019 The Authors. Medical Physics published by Wiley Periodicals, Inc. on behalf of American Association of Physicists in Medicine.)

In this study, a Monte Carlo (MC)-based beam model for an ELEKTA linear accelerator was established. The beam model is based on the EGSnrc Monte Carlo code, whereby electron beams with nominal energies of 10, 12 and 15 MeV were considered. For collimation of the electron beam, only the integrated photon multi-leaf-collimators (MLCs) were used. No additional secondary or tertiary add-ons like applicators, cutouts or dedicated electron MLCs were included. The source parameters of the initial electron beam were derived semi-automatically from measurements of depth-dose curves and lateral profiles in a water phantom. A routine to determine the initial electron energy spectra was developed which fits a Gaussian spectrum to the most prominent features of depth-dose curves. The comparisons of calculated and measured depth-dose curves demonstrated agreement within 1%/1 mm. The source divergence angle of initial electrons was fitted to lateral dose profiles beyond the range of electrons, where the imparted dose is mainly due to bremsstrahlung produced in the scattering foils. For accurate modelling of narrow beam segments, the influence of air density on dose calculation was studied. The air density for simulations was adjusted to local values (433 m above sea level) and compared with the standard air supplied by the ICRU data set. The results indicate that the air density is an influential parameter for dose calculations. Furthermore, the default value of the BEAMnrc parameter 'skin depth' for the boundary crossing algorithm was found to be inadequate for the modelling of small electron fields. A higher value for this parameter eliminated discrepancies in too broad dose profiles and an increased dose along the central axis. The beam model was validated with measurements, whereby an agreement mostly within 3%/3 mm was found.