Verification of a commercial implementation of the Macro-Monte-Carlo electron dose calculation algorithm using the virtual accelerator approach

Optimal values of the Electron Monte Carlo dose engine parameters

Background

The aim of this study was to indicate the most favorable - in terms of to the time of calculation and the uncertainty of determining the dose distribution - values of the parameters for the Electron Monte Carlo (eMC) algorithm in the Eclipse treatment planning system.

Materials and methods

Using the eMC algorithm and the variability of the values of its individual parameters, calculations of the electron dose distribution in the full-scattering virtual water phantom were performed, obtaining percentage depth doses, beam profiles, absolute dose values in points and calculation times. The reference data included water tank measurements such as relative dose distributions and absolute point doses.

Results

For 63 sets of calculation data created from selected values of the parameters for the eMC algorithm, calculation times were analyzed and the absolute calculated and measured doses were compared. Performing a statistical analysis made it possible to determine whether the differences in the values of deviations between the actual dose and the calculated dose in individual regions of the percentage depth dose curve and the beam profile are statistically significant between the analyzed sets of parameters.

Conclusions

Taking into account obtained results from the analysis of the discrepancy between the distribution of the calculated and measured dose, the correspondence of the absolute value of the calculated and measured dose and the duration of the calculation of the dose distribution, the optimal set of parameters was indicated for the eMC algorithm which allows obtaining the dose distribution and the number of monitor units in an acceptable time.

Stability analysis of a deterministic dose calculation for MRI-guided radiotherapy

Modern effort in radiotherapy to address the challenges of tumor localization and motion has led to the development of MRI guided radiotherapy technologies. Accurate dose calculations must properly account for the effects of the MRI magnetic fields. Previous work has investigated the accuracy of a deterministic linear Boltzmann transport equation (LBTE) solver that includes magnetic field, but not the stability of the iterative solution method. In this work, we perform a stability analysis of this deterministic algorithm including an investigation of the convergence rate dependencies on the magnetic field, material density, energy, and anisotropy expansion. The iterative convergence rate of the continuous and discretized LBTE including magnetic fields is determined by analyzing the spectral radius using Fourier analysis for the stationary source iteration (SI) scheme. The spectral radius is calculated when the magnetic field is included (1) as a part of the iteration source, and (2) inside the streaming-collision operator. The non-stationary Krylov subspace solver GMRES is also investigated as a potential method to accelerate the iterative convergence, and an angular parallel computing methodology is investigated as a method to enhance the efficiency of the calculation. SI is found to be unstable when the magnetic field is part of the iteration source, but unconditionally stable when the magnetic field is included in the streaming-collision operator. The discretized LBTE with magnetic fields using a space-angle upwind stabilized discontinuous finite element method (DFEM) was also found to be unconditionally stable, but the spectral radius rapidly reaches unity for very low-density media and increasing magnetic field strengths indicating arbitrarily slow convergence rates. However, GMRES is shown to significantly accelerate the DFEM convergence rate showing only a weak dependence on the magnetic field. In addition, the use of an angular parallel computing strategy is shown to potentially increase the efficiency of the dose calculation.

Full Monte Carlo and measurement-based overall performance assessment of improved clinical implementation of eMC algorithm with emphasis on lower energy range

New version 13.6.23 of the electron Monte Carlo (eMC) algorithm in Varian Eclipse™ treatment planning system has a model for 4MeV electron beam and some general improvements for dose calculation. This study provides the first overall accuracy assessment of this algorithm against full Monte Carlo (MC) simulations for electron beams from 4MeV to 16MeV with most emphasis on the lower energy range. Beams in a homogeneous water phantom and clinical treatment plans were investigated including measurements in the water phantom. Two different material sets were used with full MC: (1) the one applied in the eMC algorithm and (2) the one included in the Eclipse™ for other algorithms. The results of clinical treatment plans were also compared to those of the older eMC version 11.0.31. In the water phantom the dose differences against the full MC were mostly less than 3% with distance-to-agreement (DTA) values within 2mm. Larger discrepancies were obtained in build-up regions, at depths near the maximum electron ranges and with small apertures. For the clinical treatment plans the overall dose differences were mostly within 3% or 2mm with the first material set. Larger differences were observed for a large 4MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3mm. The discrepancies between the eMC and the full MC were generally larger for the second material set. The version 11.0.31 performed always inferiorly, when compared to the 13.6.23.

Monte Carlo N Particle code - Dose distribution of clinical electron beams in inhomogeneous phantoms

Electron dose distributions calculated using the currently available analytical methods can be associated with large uncertainties. The Monte Carlo method is the most accurate method for dose calculation in electron beams. Most of the clinical electron beam simulation studies have been performed using non- MCNP [Monte Carlo N Particle] codes. Given the differences between Monte Carlo codes, this work aims to evaluate the accuracy of MCNP4C-simulated electron dose distributions in a homogenous phantom and around inhomogeneities. Different types of phantoms ranging in complexity were used; namely, a homogeneous water phantom and phantoms made of polymethyl methacrylate slabs containing different-sized, low- and high-density inserts of heterogeneous materials. Electron beams with 8 and 15 MeV nominal energy generated by an Elekta Synergy linear accelerator were investigated. Measurements were performed for a 10 cm × 10 cm applicator at a source-to-surface distance of 100 cm. Individual parts of the beam-defining system were introduced into the simulation one at a time in order to show their effect on depth doses. In contrast to the first scattering foil, the secondary scattering foil, X and Y jaws and applicator provide up to 5% of the dose. A 2%/2 mm agreement between MCNP and measurements was found in the homogenous phantom, and in the presence of heterogeneities in the range of 1-3%, being generally within 2% of the measurements for both energies in a "complex" phantom. A full-component simulation is necessary in order to obtain a realistic model of the beam. The MCNP4C results agree well with the measured electron dose distributions.

Validation of an electron Monte Carlo dose calculation algorithm in the presence of heterogeneities using EGSnrc and radiochromic film measurements

The purpose of this study is to validate Eclipse's electron Monte Carlo algorithm (eMC) in heterogeneous phantoms using radiochromic films and EGSnrc as a reference Monte Carlo algorithm. Four heterogeneous phantoms are used in this study. Radiochromic films are inserted in these phantoms, including in heterogeneous media, and the measured relative dose distributions are compared to eMC calculations. Phantoms A, B, and C contain 1D heterogeneities, built with layers of lung‐ (phantom A) and bone‐ (phantoms B and C) equivalent materials sandwiched in Plastic Water. Phantom D is a thorax‐anthropomorphic phantom with 2D lung heterogeneities. Electron beams of 6, 9, 12 and 18 MeV from a Varian Clinac 2100 are delivered to these phantoms with a 10×10 cm2 applicator. Monte Carlo simulations with an independent algorithm (EGSnrc) are also used as a reference tool for two purposes: (1) as a second validation of the eMC dose calculations, and (2) to calculate the stopping power ratio between radiochromic films and bone medium, when dose is measured inside the heterogeneity. Percent depth dose (PDD) film measurements and eMC calculations agree within 2% or 3 mm for phantom A, and within 3% or 3 mm for phantoms B and C for almost all beam energies. One exception is observed with phantom B and the 6 MeV, where measured PDDs and those calculated with eMC differ by up to 4 mm. Gamma analysis of the measured and calculated 2D dose distributions in phantom D agree with criteria of 3%, 3 mm for 9, 12, and 18 MeV beams, and criteria of 5%, 3 mm for the 6 MeV beam. Dose calculations in heterogeneous media with eMC agree within 3% or 3 mm with radiochromic film measurements. Six (6) MeV beams are not modeled as accurately as other beam energies. The eMC algorithm is suitable for clinical dose calculations involving lung and bone.

PACS numbers: 87.10.Rt, 87.55.km

On the wall perturbation correction for a parallel-plate NACP-02 chamber in clinical electron beams

PURPOSE

In recent years, several Monte Carlo studies have been published concerning the perturbation corrections of a parallel-plate chamber in clinical electron beams. In these studies, a strong depth dependence of the relevant correction factors (p(wall) and P(cav)) for depth beyond the reference depth is recognized and it has been shown that the variation with depth is sensitive to the choice of the chamber's effective point of measurement. Recommendations concerning the positioning of parallel-plate ionization chambers in clinical electron beams are not the same for all current dosimetry protocols. The IAEA TRS-398 as well as the IPEM protocol and the German protocol DIN 6800-2 interpret the depth of measurement within the phantom as the water equivalent depth, i.e., the nonwater equivalence of the entrance window has to be accounted for by shifting the chamber by an amount deltaz. This positioning should ensure that the primary electrons traveling from the surface of the water phantom through the entrance window to the chamber's reference point sustain the same energy loss as the primary electrons in the undisturbed phantom. The objective of the present study is the determination of the shift deltaz for a NACP-02 chamber and the calculation of the resulting wall perturbation correction as a function of depth. Moreover, the contributions of the different chamber walls to the wall perturbation correction are identified.

METHODS

The dose and fluence within the NACP-02 chamber and a wall-less air cavity is calculated using the Monte Carlo code EGSnrc in a water phantom at different depths for different clinical electron beams. In order to determine the necessary shift to account for the nonwater equivalence of the entrance window, the chamber is shifted in steps deltaz around the depth of measurement. The optimal shift deltaz is determined from a comparison of the spectral fluence within the chamber and the bare cavity. The wall perturbation correction is calculated as the ratio between doses for the complete chamber and a wall-less air cavity.

RESULTS

The high energy part of the fluence spectra within the chamber strongly varies even with small chamber shifts, allowing the determination of deltaz within micrometers. For the NACP-02 chamber a shift deltaz = -0.058 cm results. This value is independent of the energy of the primary electrons as well as of the depth within the phantom and it is in good agreement with the value recommended in the German dosimetry protocol. Applying this shift, the calculated wall perturbation correction as a function of depth is varying less than 1% from zero up to the half value depth R50 for electron energies in the range of 6-21 MeV. The remaining depth dependence can mainly be attributed to the scatter properties of the entrance window. When neglecting the nonwater equivalence of the entrance window, the variation of p(wall) with depth is up to 10% and more, especially for low electron energies.

CONCLUSIONS

The variation of the wall perturbation correction for the NACP-02 chamber in clinical electron beams strongly depends on the positioning of the chamber. Applying a shift deltaz = -0.058 cm toward the focus ensures that the primary electron spectrum within the chamber bears the largest resemblance to the fluence of a wall-less cavity. Hence, the influence of the chamber walls on the perturbation correction can be separated out and the residual variation of p(wall) with depth is minimized.