Auto-commissioning of a Monte Carlo electron beam model with application to photon MLC shaped electron fields

Objective.Presently electron beam treatments are delivered using dedicated applicators. An alternative is the usage of the already installed photon multileaf collimator (pMLC) enabling efficient electron treatments. Currently, the commissioning of beam models is a manual and time-consuming process. In this work an auto-commissioning procedure for the Monte Carlo (MC) beam model part representing the beam above the pMLC is developed for TrueBeam systems with electron energies from 6 to 22 MeV.Approach.The analytical part of the electron beam model includes a main source representing the primary beam and a jaw source representing the head scatter contribution each consisting of an electron and a photon component, while MC radiation transport is performed for the pMLC. The auto-commissioning of this analytical part relies on information pre-determined from MC simulations, in-air dose profiles and absolute dose measurements in water for different field sizes and source to surface distances (SSDs). For validation calculated and measured dose distributions in water were compared for different field sizes, SSDs and beam energies for eight TrueBeam systems. Furthermore, a sternum case in an anthropomorphic phantom was considered and calculated and measured dose distributions were compared at different SSDs.Main results.Instead of the manual commissioning taking up to several days of calculation time and several hours of user time, the auto-commissioning is carried out in a few minutes. Measured and calculated dose distributions agree generally within 3% of maximum dose or 2 mm. The gamma passing rates for the sternum case ranged from 96% to 99% (3% (global)/2 mm criteria, 10% threshold).Significance.The auto-commissioning procedure was successfully implemented and applied to eight TrueBeam systems. The newly developed user-friendly auto-commissioning procedure allows an efficient commissioning of an MC electron beam model and eases the usage of advanced electron radiotherapy utilizing the pMLC for beam shaping.

Technical Note: Parallel implementation of the EGSnrc Monte Carlo simulation of ionizing radiation transport using OpenMP

PURPOSE

To present the implementation of a new option for parallel processing of the EGSnrc Monte Carlo system using the OpenMP API, as an alternative to the provided method based on the use of a batch queuing system (BQS).

METHODS

The parallel solution presented, called OMP_EGS, makes use of OpenMP features to control the workload distribution between the compute units. These features were inserted into the original EGSnrc source code through properly defined macros. In order to validate the platform, the possibility of producing results in exact agreement with the serial implementation was assessed. The performance of OMP_EGS was evaluated against the BQS method, in terms of parallel speedup and efficiency.

RESULTS

As the OpenMP features can be activated or deactivated depending on the compilation options, the implementation of the platform allowed the direct recovery of the original serial implementation. The validation tests showed that OMP_EGS was able to reproduce the exact same results as the serial implementation. The performance and scalability tests showed that OMP_EGS is a better alternative than the EGSnrc BQS parallel implementation, both in terms of runtime and parallel efficiency.

CONCLUSIONS

The presented solution has several advantages over the BQS-based parallel implementation available for the EGSnrc system. One of the main advantages is that, in contrast to the BQS alternative, it can be implemented using different compilers and operative systems, which turns it into a compact and portable solution that can be used on a wide range of working environments. It does not introduce artifacts on the simulated distributions, as it only handles the distribution of work among the available computing resources and it proved to have a better performance.

MCTP system model based on linear programming optimization of apertures obtained from sequencing patient image data maps

PURPOSE

The authors present a hybrid direct multileaf collimator (MLC) aperture optimization model exclusively based on sequencing of patient imaging data to be implemented on a Monte Carlo treatment planning system (MC-TPS) to allow the explicit radiation transport simulation of advanced radiotherapy treatments with optimal results in efficient times for clinical practice.

METHODS

The planning system (called CARMEN) is a full MC-TPS, controlled through aMATLAB interface, which is based on the sequencing of a novel map, called "biophysical" map, which is generated from enhanced image data of patients to achieve a set of segments actually deliverable. In order to reduce the required computation time, the conventional fluence map has been replaced by the biophysical map which is sequenced to provide direct apertures that will later be weighted by means of an optimization algorithm based on linear programming. A ray-casting algorithm throughout the patient CT assembles information about the found structures, the mass thickness crossed, as well as PET values. Data are recorded to generate a biophysical map for each gantry angle. These maps are the input files for a home-made sequencer developed to take into account the interactions of photons and electrons with the MLC. For each linac (Axesse of Elekta and Primus of Siemens) and energy beam studied (6, 9, 12, 15 MeV and 6 MV), phase space files were simulated with the EGSnrc/BEAMnrc code. The dose calculation in patient was carried out with the BEAMDOSE code. This code is a modified version of EGSnrc/DOSXYZnrc able to calculate the beamlet dose in order to combine them with different weights during the optimization process.

RESULTS

Three complex radiotherapy treatments were selected to check the reliability of CARMEN in situations where the MC calculation can offer an added value: A head-and-neck case (Case I) with three targets delineated on PET/CT images and a demanding dose-escalation; a partial breast irradiation case (Case II) solved with photon and electron modulated beams (IMRT + MERT); and a prostatic bed case (Case III) with a pronounced concave-shaped PTV by using volumetric modulated arc therapy. In the three cases, the required target prescription doses and constraints on organs at risk were fulfilled in a short enough time to allow routine clinical implementation. The quality assurance protocol followed to check CARMEN system showed a high agreement with the experimental measurements.

CONCLUSIONS

A Monte Carlo treatment planning model exclusively based on maps performed from patient imaging data has been presented. The sequencing of these maps allows obtaining deliverable apertures which are weighted for modulation under a linear programming formulation. The model is able to solve complex radiotherapy treatments with high accuracy in an efficient computation time.

Intensity- and energy-modulated electron radiotherapy by means of an xMLC for head and neck shallow tumors

The purpose of this paper is to assess the feasibility of delivering intensity- and energy-modulated electron radiation treatment (MERT) by a photon multileaf collimator (xMLC) and to evaluate the improvements obtained in shallow head and neck (HN) tumors. Four HN patient cases covering different clinical situations were planned by MERT, which used an in-house treatment planning system that utilized Monte Carlo dose calculation. The cases included one oronasal, two parotid and one middle ear tumors. The resulting dose-volume histograms were compared with those obtained from conventional photon and electron treatment techniques in our clinic, which included IMRT, electron beam and mixed beams, most of them using fixed-thickness bolus. Experimental verification was performed with plane-parallel ionization chambers for absolute dose verification, and a PTW ionization chamber array and radiochromic film for relative dosimetry. A MC-based treatment planning system for target with compromised volumes in depth and laterally has been validated. A quality assurance protocol for individual MERT plans was launched. Relative MC dose distributions showed a high agreement with film measurements and absolute ion chamber dose measurements performed at a reference point agreed with MC calculations within 2% in all cases. Clinically acceptable PTV coverage and organ-at-risk sparing were achieved by using the proposed MERT approach. MERT treatment plans, based on delivery of intensity-modulated electron beam using the xMLC, for superficial head and neck tumors, demonstrated comparable or improved PTV dose homogeneity with significantly lower dose to normal tissues. The clinical implementation of this technique will be able to offer a viable alternative for the treatment of shallow head and neck tumors.

Modulated electron radiotherapy treatment planning using a photon multileaf collimator for post-mastectomized chest walls

BACKGROUND AND PURPOSE

To evaluate the feasibility of using a photon MLC (xMLC) for modulated electron radiotherapy treatment (MERT) as an alternative to conventional post-mastectomy chest wall (CW) irradiation. A Monte Carlo (MC) based planning system was developed to overcome the inaccuracy of the 'pencil beam' algorithm. MC techniques are known to accurately calculate the dose distributions of electron beams, allowing the explicit simulation of electron interactions within the MLC.

MATERIALS AND METHODS

Four real clinical CW cases were planned using MERT which were compared with the conventional electron treatments based on blocks and by a straightforward approach using the MLC, and not the blocks (as an intermediate step to MERT) to shape the same segments with SSD between 60 and 70 cm depending on PTV size. MC calculations were verified with an array of ionization chambers and radiochromic films in a solid water phantom.

RESULTS

Tests based on gamma analysis between MC dose distributions and radiochromic film measurements showed an excellent agreement. Differences in the absolute dose measured with a plane-parallel chamber at a reference point were below 3% for all cases. MERT solution showed a better PTV coverage and a significant reduction of the doses to the organs at risk (OARs).

CONCLUSION

MERT can effectively improve the current electron treatments by obtaining a better PTV coverage and sparing healthy tissues. More directly, block-shaped treatments could be replaced by MLC-shaped non-modulated segments providing similar results.

Micro ionization chamber dosimetry in IMRT verification: Clinical implications of dosimetric errors in the PTV

BACKGROUND AND PURPOSE

Absolute dose measurements for Intensity Modulated Radiotherapy (IMRT) beamlets is difficult due to the lack of lateral electron equilibrium. Recently we found that the absolute dosimetry in the penumbra region of the IMRT beamlet, can suffer from significant errors (Capote et al., Med Phys 31 (2004) 2416-2422). This work has the goal to estimate the error made when measuring the Planning Target Volume's (PTV) absolute dose by a micro ion chamber (microIC) in typical IMRT treatment. The dose error comes from the assumption that the dosimetric parameters determining the absolute dose are the same as for the reference conditions.

MATERIALS AND METHODS

Two IMRT treatment plans for common prostate carcinoma case, derived by forward and inverse optimisation, were considered. Detailed geometrical simulation of the microIC and the dose verification set-up was performed. The Monte Carlo (MC) simulation allows us to calculate the delivered dose to water and the dose delivered to the active volume of the ion chamber. However, the measured dose in water is usually derived from chamber readings assuming reference conditions. The MC simulation provides needed correction factors for ion chamber dosimetry in non reference conditions.

RESULTS

Dose calculations were carried out for some representative beamlets, a combination of segments and for the delivered IMRT treatments. We observe that the largest dose errors (i.e. the largest correction factors) correspond to the smaller contribution of the corresponding IMRT beamlets to the total dose delivered in the ionization chamber within PTV.

CONCLUSION

The clinical impact of the calculated dose error in PTV measured dose was found to be negligible for studied IMRT treatments.

Microionization chamber for reference dosimetry in IMRT verification: clinical implications on OAR dosimetric errors

Intensity modulated radiotherapy (IMRT) has become a treatment of choice in many oncological institutions. Small fields or beamlets with sizes of 1 to 5 cm2 are now routinely used in IMRT delivery. Therefore small ionization chambers (IC) with sensitive volumes 0.1 cm3 are generally used for dose verification of an IMRT treatment. The measurement conditions during verification may be quite different from reference conditions normally encountered in clinical beam calibration, so dosimetry of these narrow photon beams pertains to the so-called non-reference conditions for beam calibration. This work aims at estimating the error made when measuring the organ at risk's (OAR) absolute dose by a micro ion chamber (microIC) in a typical IMRT treatment. The dose error comes from the assumption that the dosimetric parameters determining the absolute dose are the same as for the reference conditions. We have selected two clinical cases, treated by IMRT, for our dose error evaluations. Detailed geometrical simulation of the microIC and the dose verification set-up was performed. The Monte Carlo (MC) simulation allows us to calculate the dose measured by the chamber as a dose averaged over the air cavity within the ion-chamber active volume (D(air)). The absorbed dose to water (D(water)) is derived as the dose deposited inside the same volume, in the same geometrical position, filled and surrounded by water in the absence of the ion chamber. Therefore, the D(water)/D(air) dose ratio is the MC estimator of the total correction factor needed to convert the absorbed dose in air into the absorbed dose in water. The dose ratio was calculated for the microIC located at the isocentre within the OARs for both clinical cases. The clinical impact of the calculated dose error was found to be negligible for the studied IMRT treatments.