Leen Paelinck

Accuracy of patient dose calculation for lung IMRT: A comparison of Monte Carlo, convolution/superposition, and pencil beam computations

The accuracy of dose computation within the lungs depends strongly on the performance of the calculation algorithm in regions of electronic disequilibrium that arise near tissue inhomogeneities with large density variations. There is a lack of data evaluating the performance of highly developed analytical dose calculation algorithms compared to Monte Carlo computations in a clinical setting. We compared full Monte Carlo calculations (performed by our Monte Carlo dose engine MCDE) with two different commercial convolution/superposition (CS) implementations (Pinnacle-CS and Helax-TMSs collapsed cone model Helax-CC) and one pencil beam algorithm (Helax-TMSs pencil beam model Helax-PB) for 10 intensity modulated radiation therapy (IMRT) lung cancer patients. Treatment plans were created for two photon beam qualities (6 and 18 MV). For each dose calculation algorithm, patient, and beam quality, the following set of clinically relevant dose-volume values was reported: (i) minimal, median, and maximal dose (Dmin, D50, and Dmax) for the gross tumor and planning target volumes (GTV and PTV); (ii) the volume of the lungs (excluding the GTV) receiving at least 20 and 30 Gy (V20 and V30) and the mean lung dose; (iii) the 33rd percentile dose (D33) and Dmax delivered to the heart and the expanded esophagus; and (iv) Dmax for the expanded spinal cord. Statistical analysis was performed by means of one-way analysis of variance for repeated measurements and Tukey pairwise comparison of means. Pinnacle-CS showed an excellent agreement with MCDE within the target structures, whereas the best correspondence for the organs at risk (OARs) was found between Helax-CC and MCDE. Results from Helax-PB were unsatisfying for both targets and OARs. Additionally, individual patient results were analyzed. Within the target structures, deviations above 5% were found in one patient for the comparison of MCDE and Helax-CC, while all differences between MCDE and Pinnacle-CS were below 5%. For both Pinnacle-CS and Helax-CC, deviations from MCDE above 5% were found within the OARs: within the lungs for two (6 MV) and six (18 MV) patients for Pinnacle-CS, and within other OARs for two patients for Helax-CC (for Dmax of the heart and D33 of the expanded esophagus) but only for 6 MV. For one patient, all four algorithms were used to recompute the dose after replacing all computed tomography voxels within the patients skin contour by water. This made all differences above 5% between MCDE and the other dose calculation algorithms disappear. Thus, the observed deviations mainly arose from differences in particle transport modeling within the lungs, and the commissioning of the algorithms was adequately performed (or the commissioning was less important for this type of treatment). In conclusion, not one pair of the dose calculation algorithms we investigated could provide results that were consistent within 5% for all 10 patients for the set of clinically relevant dose-volume indices studied. As the results from both CS algorithms differed significantly, care should be taken when evaluating treatment plans as the choice of dose calculation algorithm may influence clinical results. Full Monte Carlo provides a great benchmarking tool for evaluating the performance of other algorithms for patient dose computations.

The influence of small field sizes, penumbra, spot size and measurement depth on perturbation factors for microionization chambers

The purpose of this study was the investigation of perturbation factors for microionization chambers in small field dosimetry and the influence of penumbra for different spot sizes. To this purpose, correlated sampling was implemented in the EGSnrc Monte Carlo (MC) user code cavity: CScavity. CScavity was first benchmarked against results in the literature for an NE2571 chamber. An efficiency increase of 17 was attained for the calculation of a realistic chamber perturbation factor in a water phantom. Calculations have been performed for microionization chambers of type PinPoint 31006 and PinPoint 31016 in full BEAMnrc linac simulations. Investigating the physical backgrounds of the differences for these small field settings, perturbation factors have been split up into (1) central electrode perturbation, (2) wall perturbation, (3) air-to-water perturbation (chamber volume air-to-water) and (4) water volume perturbation (water chamber volume to 1 mm(3) voxel). The influence of different spot sizes, position in penumbra, measuring depth and detector geometry on these perturbation factors has been investigated, in a 0.8 x 0.8 cm(2) field setting. p(cel) for the PP31006 steel electrode shows a variation of up to 1% in the lateral position, but only 0.4% for the PP31016 with an Al electrode. The air-to-water perturbation in the optimal scanning direction for both profiles and depth is most influenced by the radiation field, and only to a small extent the chamber geometry. The PP31016 geometry (shorter, larger radius) requires less total perturbation within the central axis of the field, but results in slightly larger variations off axis in the optimal scanning direction. Smaller spot sizes (0.6 mm FWHM) and sharper penumbras, compared to larger spot sizes (2 mm FWHM), result in larger perturbation starting in the penumbra. The longer geometries of the PP31006/14/15 exhibit in the non-optimal scanning direction large variations in total perturbation (p(tot) 1.201(4) (0.6 mm spot, 3 mm off axis, type A MC uncertainty) to 0.803(4) (5 mm off axis)) mainly due to volume perturbation. Therefore in IMRT settings, when the detector is not always in the optimal scanning direction, the PP31016 geometry requires less extreme perturbation (max p(tot) 1.130(3)) and shows less variation. However, these results suggest that small variations in positioning, spot size or MLC result in large differences in perturbation factors. Therefore even these 0.016 cm(3) ionization chambers are limited in their use for a field setting of 0.8 x 0.8 cm(2), as used in this investigation.

Experimental verification of lung dose with radiochromic film: comparison with Monte Carlo simulations and commercially available treatment planning systems

The purpose of this study was to assess the absorbed dose in and around lung tissue by performing radiochromic film measurements, Monte Carlo simulations and calculations with superposition convolution algorithms. We considered a layered polystyrene phantom of 12 x 12 x 12 cm3 containing a central cavity of 6 x 6 x 6 cm3 filled with Gammex RMI lung-equivalent material. Two field configurations were investigated, a small 1 x 10 cm2 field and a larger 10 x 10 cm2 field. First, we performed Monte Carlo simulations to investigate the influence of radiochromic film itself on the measured dose distribution when the film intersects a lung-equivalent region and is oriented parallel to the central beam axis. To that end, the film and the lung-equivalent materials were modelled in detail, taking into account their specific composition. Next, measurements were performed with the film oriented both parallel and perpendicular to the central beam axis to verify the results of our Monte Carlo simulations. Finally, we digitized the phantom in two commercially available treatment planning systems, Helax-TMS version 6.1A and Pinnacle version 6.2b, and calculated the absorbed dose in the phantom with their incorporated superposition convolution algorithms to compare with the Monte Carlo simulations. Comparing Monte Carlo simulations with measurements reveals that radiochromic film is a reliable dosimeter in and around lung-equivalent regions when the film is positioned perpendicular to the central beam axis. Radiochromic film is also able to predict the absorbed dose accurately when the film is positioned parallel to the central beam axis through the lung-equivalent region. However, attention must be paid when the film is not positioned along the central beam axis, in which case the film gradually attenuates the beam and decreases the dose measured behind the cavity. This underdosage disappears by offsetting the film a few centimetres. We find deviations of about 3.6% between Monte Carlo and the superposition convolution algorithm of Pinnacle behind the lung region, for both field configurations. Pinnacle is quite accurate in the lung region. Deviations up to 5.6% for the small field are found in the lung region between Monte Carlo and the superposition convolution algorithm of Helax-TMS. Behind the lung region, Helax-TMS is in better agreement with Monte Carlo. Radiochromic film measurements or Monte Carlo simulations are reliable methods to establish the dose in and around lung tissue.

Intensity-Modulated Arc Therapy with Simultaneous Integrated Boost in the Treatment of Primary Irresectable Cervical Cancer

Purpose:To report on the planning procedure, quality control, and clinical implementation of intensity-modulated arc therapy (IMAT) delivering a simultaneous integrated boost (SIB) in patients with primary irresectable cervix carcinoma.Patients and Methods:Six patients underwent PET-CT (positron emission tomography-computed tomography) and MRI (magnetic resonance imaging) before treatment planning. Prescription (25 fractions) was(1) a median dose (D50) of 62, 58 and 56 Gy to the primary tumor (GTV_cervix), primary clinical target volume (CTV_cervix) and its planning target volume (PTV_cervix), respectively;(2) a D50 of 60 Gy to the PET-positive lymph nodes (GTV_nodes);(3) a minimal dose (D98) of 45 Gy to the planning target volume of the elective lymph nodes (PTV_nodes).IMAT plans were generated using an anatomy-based exclusion tool with the aid of weight and leaf position optimization. The dosimetric delivery of IMAT was validated preclinically using radiochromic film dosimetry.Results:Five to nine arcs were needed to create valid IMAT plans. Dose constraints on D50 were not met in two patients (both GTV_cervix: 1 Gy and 3 Gy less). D98 for PTV_nodes was not met in three patients (1 Gy each). Film dosimetry showed excellent gamma evaluation. There were no treatment interruptions.Conclusion:IMAT allows delivering an SIB to the macroscopic tumor without compromising the dose to the elective lymph nodes or the organs at risk. The clinical implementation is feasible.Ziel:Evaluation einer intensitätsmodulierten Rotationsbestrahlung (IMAT) mit Applikation eines simultanen integrierten Boosts (SIB) zur primären Behandlung des fortgeschrittenen Zervixkarzinoms.Patienten und Methodik:Sechs Patientinnen mit einem fortgeschrittenen Zervixkarzinom wurden einer MRT- (Magnetresonanztomographie) und PET-CT-basierten (Positronenemissionstomographie-Computertomographie) Bestrahlungsplanung für eine IMAT unterzogen und bestrahlt. Das Dosis-Zeit-Muster wurde, bezogen auf die entsprechenden Zielvolumina für 25 Fraktionen, wie folgt festgelegt:1. Eine mediane Dosis (D50) von 62 Gy, 58 Gy und 56 Gy wurde für das makroskopische Zervixkarzinom (GTV), das klinische Zielvolumen (CTV) und das Planungszielvolumen (PTV) verschrieben.2. Eine mediane Dosis von 60 Gy wurde für die PET-positiven regionären Lymphknoten festgelegt.3. Elektiv zu bestrahlende regionäre Lymphknoten sollten eine minimale Dosis (D98) von 45 Gy erhalten.Die IMAT-Pläne wurden mit Hilfe eines anatomiebasierten Ausschlussalgorithmus durch Optimierung und Wichtung von Leafpositionen erzeugt. Die präklinische Dosimetrie erfolgte mittels Filmdosimetrie.Ergebnisse:Insgesamt fünf bis neun Rotationsfelder waren zur Erzeugung geeigneter IMAT-Pläne erforderlich. Bei zwei Patientinnen war die angestrebte Dosis für das makroskopische Zervixkarzinom (GTV) 1 Gy und 3 Gy zu niedrig. In drei Fällen wurde die minimale Dosis (D98) an den elektiv zu behandelnden Lymphknoten um je 1 Gy unterschritten. Die Daten der Filmdosimetrie zeigten eine ausgezeichnete Gammabewertung. Die Bestrahlung konnte in allen Fällen ohne Unterbrechung appliziert werden.Schlussfolgerung:Die klinische Umsetzung der IMAT mit SIB des Zervixkarzinoms ist ohne Dosiskompromisse an elektiven Lymphknotenstationen und Risikoorganen möglich.

Application of monomer/polymer gel dosimetry to study the effects of tissue inhomogeneities on intensity-modulated radiation therapy (IMRT) dose distributions

BACKGROUND AND PURPOSE When planning an intensity-modulated radiation therapy (IMRT) treatment in a heterogeneous region (e.g. the thorax), the dose computation algorithm of a treatment planning system may need to account for these inhomogeneities in order to obtain a reliable prediction of the dose distribution. An accurate dose verification technique such as monomer/polymer gel dosimetry is suggested to verify the outcome of the planning system. MATERIALS AND METHODS The effects of low-density structures: (a) on narrow high-energy (18 MV) photon beams; and (b) on a class-solution IMRT treatment delivered to a thorax phantom have been examined using gel dosimetry. The used phantom contained air cavities that could be filled with water to simulate a homogeneous or heterogeneous configuration. The IMRT treatment for centrally located lung tumors was delivered on both cases, and gel derived dose maps were compared with computations by both the GRATIS and Helax-TMS planning system. RESULTS Dose rebuildup due to electronic disequilibrium in a narrow photon beam is demonstrated. The gel measurements showed good agreement with diamond detector measurements. Agreement between measured IMRT dose maps and dose computations was demonstrated by several quantitative techniques. An underdosage of the planning target volume (PTV) was revealed. The homogeneity of the phantom had only a minor influence on the dose distribution in the PTV. An expansion of low-level isodoses in the lung volume was predicted by collapsed cone computations in the heterogeneous case. CONCLUSIONS For the class-solution described, the dose in centrally located mediastinal tumors can be computed with sufficient accuracy, even when neglecting the lower lung density. Polymer gel dosimetry proved to be a valuable technique to verify dose calculation algorithms for IMRT in 3D in heterogeneous configurations.

Monte Carlo modeling of the ModuLeaf miniature MLC for small field dosimetry and quality assurance of the clinical treatment planning system

The purpose of this investigation was the verification of both the measured data and quality of the implementation of the add-on ModuLeaf miniature multileaf collimator (ML mMLC) into the clinical treatment planning system for conformal stereotactic radiosurgery treatment. To this end the treatment head with ML mMLC was modeled in the BEAMnrc Monte Carlo (MC) code. The 6 MV photon beams used in the setup were first benchmarked with a set of measurements. A total ML mMLC transmission of 1.13% of the 10 x 10 cm2 open field dose was measured and reproduced with the BEAMnrc/DOSXYZnrc code. Correspondence between calculated and measured output factors (OFs) was within 2%. Correspondence between MC and measured profiles was within 2% dose and 2 mm distance, only for the smallest 0.5 x 0.5 cm2 field the results were within 3% dose. In the next step, the MC model was compared with Gafchromic film measurements and Pinnacle(3) 7.4 f (convolution superposition algorithm) calculated dose distributions, using a gamma evaluation comparison, for a multi-beam patient setup delivered to a Lucytrade mark phantom. The gamma evaluation of the MC versus Gafchromic film resulted in 3.4% of points not fulfilling gamma <or= 1 for a 2%/2 mm criterion, the Pinnacle(3) 7.4 f versus Gafchromic results 3.8% and Pinnacle versus MC less than 1%. For specific patients with lesions of 8 cc and 0.2 cc, Monte Carlo and Pinnacle simulations of the plans were performed and compared using DVH evaluation. DVHs corresponded within 2% dose and 2% volume.

MCDE: a new Monte Carlo dose engine for IMRT

A new accurate Monte Carlo code for IMRT dose computations, MCDE (Monte Carlo dose engine), is introduced. MCDE is based on BEAMnrc/DOSXYZnrc and consequently the accurate EGSnrc electron transport. DOSXYZnrc is reprogrammed as a component module for BEAMnrc. In this way both codes are interconnected elegantly, while maintaining the BEAM structure and only minimal changes to BEAMnrc.mortran are necessary. The treatment head of the Elekta SLiplus linear accelerator is modelled in detail. CT grids consisting of up to 200 slices of 512 x 512 voxels can be introduced and up to 100 beams can be handled simultaneously. The beams and CT data are imported from the treatment planning system GRATIS via a DICOM interface. To enable the handling of up to 50 x 10(6) voxels the system was programmed in Fortran95 to enable dynamic memory management. All region-dependent arrays (dose, statistics, transport arrays) were redefined. A scoring grid was introduced and superimposed on the geometry grid, to be able to limit the number of scoring voxels. The whole system uses approximately 200 MB of RAM and runs on a PC cluster consisting of 38 1.0 GHz processors. A set of in-house made scripts handle the parallellization and the centralization of the Monte Carlo calculations on a server. As an illustration of MCDE, a clinical example is discussed and compared with collapsed cone convolution calculations. At present, the system is still rather slow and is intended to be a tool for reliable verification of IMRT treatment planning in the case of the presence of tissue inhomogeneities such as air cavities.

The importance of accurate linear accelerator head modelling for IMRT Monte Carlo calculations

Two Monte Carlo dose engines for radiotherapy treatment planning, namely a beta release of Peregrine and MCDE (Monte Carlo dose engine), were compared with Helax-TMS (collapsed cone superposition convolution) for a head and neck patient for the Elekta SLi plus linear accelerator. Deviations between the beta release of Peregrine and MCDE up to 10% were obtained in the dose volume histogram of the optical chiasm. It was illustrated that the differences are not caused by the particle transport in the patient, but by the modelling of the Elekta SLi plus accelerator head and more specifically the multileaf collimator (MLC). In MCDE two MLC modules (MLCQ and MLCE) were introduced to study the influence of the tongue-and-groove geometry, leaf bank tilt and leakage on the actual dose volume histograms. Differences in integral dose in the optical chiasm up to 3% between the two modules have been obtained. For single small offset beams though the FWHM of lateral profiles obtained with MLCE can differ by more than 1.5 mm from profiles obtained with MLCQ. Therefore, and because the recent version of MLCE is as fast as MLCQ, we advise to use MLCE for modelling the Elekta MLC. Nevertheless there still remains a large difference (up to 10%) between Peregrine and MCDE. By studying small offset beams we have shown that the profiles obtained with Peregrine are shifted, too wide and too flat compared with MCDE and phantom measurements. The overestimated integral doses for small beam segments explain the deviations observed in the dose volume histograms. The Helax-TMS results are in better agreement with MCDE, although deviations exceeding 5% have been observed in the optical chiasm. Monte Carlo dose deviations of more than 10% as found with Peregrine are unacceptable as an influence on the clinical outcome is possible and as the purpose of Monte Carlo treatment planning is to obtain an accuracy of 2%. We would like to emphasize that only the Elekta MLC has been tested in this work, so it is certainly possible that alpha releases of Peregrine provide more accurate results for other accelerators.

Comparison of build-up dose between Elekta and Varian linear accelerators for high-energy photon beams using radiochromic film and clinical implications for IMRT head and neck treatments

Skin toxicity has been reported for IMRT of head and neck cancer. The purpose of this study was to investigate the dose in the build-up region delivered by a 6 MV treatment plan for which important skin toxicity was observed. We also investigated if the different designs of the treatment head of an Elekta and a Varian linear accelerator, especially the lower position of the Varian multi-leaf collimator, give rise to different build-up doses. For regular square open beams, the build-up dose along the central beam axis is higher for the Varian machine than for the Elekta machine, both for 6 MV and 18 MV. At the Elekta machine at 18 MV, the superficial dose of a diamond shaped 10 x 10 cm2 field is 3.6% lower than the superficial dose of a regular 10 x 10 cm2 field. This effect is not seen at 6 MV. At the Varian machine, the superficial dose of the diamond shaped field is respectively 3.5 and 14.2% higher than the superficial dose of the regular 10 x 10 cm2 field for 6 MV and 18 MV. Despite the differences measured in build-up dose for single beams between the Elekta and the Varian linear accelerator, there were no measurable differences in superficial dose when a typical IMRT dose plan of 6 MV for a head and neck tumour is executed at the two machines.

Adaptive radiotherapy for locally advanced non-small cell lung cancer, can we predict when and for whom?

ABSTRACT Background. Adaptive radiotherapy (ART) could be a tool to reduce toxicity and to facilitate dose escalation in stage III NSCLC. Our aim was to identify the most appropriate time and potential benefit of ART. Material and methods. We analyzed volume reduction and dosimetric consequences of 41 patients who were treated with concurrent (cCRT) (n = 21) or sequential (sCRT) chemoradiotherapy to a median dose of 70 Gy, 2 Gy/F. At every treatment fraction a cone-beam CT (CBCT) was performed. The gross tumor volume (GTV-T) was adapted (exclusion of lymph nodes) to create the GTV-T-F1. Every fifth fraction (F5–F30), the GTV-T-F1 was adapted on the CBCT to create a GTV-T-Fx. Dose volume histograms were recalculated for every GTV-T-Fx, enabling to create lookup tables to predict the theoretical dosimetric advantage on common lung dose constraints. Results. The average GTV reduction was 42.1% (range 4.0–69.3%); 50.1% and 33.7% for the cCRT and sCRT patients, respectively. A linear relationship between GTV-T-F1 volume and absolute volume decrease was found for both groups. The mean V5, V20, V30 and mean lung dose increased by 0.8, 3.1, 5.2 and 3.4%, respectively. A larger increase (p < 0.05) was observed for peripheral tumors and cCRT. Lookup tables were generated. Conclusion. ART offers the most beneficial dosimetric effects when performed around fraction 15, especially for patients with a large initial GTV-T treated by cCRT.