Nick Reynaert

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.

On the accuracy of monomer/polymer gel dosimetry in the proximity of a high-dose-rate 192Ir source.

The aim of this work was to investigate the applicability of MR-based polymer gel dosimetry to measure the absorbed dose distribution at short distance from an iridium-192 brachytherapy point source. In this paper, different methodological problems that may result in significant errors in the measured dose distribution are discussed. First of all the extent to which physicochemical mechanisms alter the dose response is discussed. The most important among these are the oxygen permeability of the catheter material and monomer-diffusion-related effects during irradiation. The effect of oxygen on the dose-R2 curve has been determined quantitatively and an oxygen map is performed using a well-defined external irradiation beam. The effect of diffusion of monomer during irradiation has been simulated. Another contribution of errors is related to magnetic susceptibility differences between the catheter and the gel during scanning the irradiated gel. The magnetic field distortion has been mapped by using both an experimental MRI technique and by simulation. Moreover, in constructing a dose-versus-distance curve by polar averaging, the sensitivity to the exact distance between source and point of measurement and to partial volume effects is illustrated. An optimization program is proposed to determine the location of the source on a sub-pixel scale.

Underdosage of the upper-airway mucosa for small fields as used in intensity-modulated radiation therapy: A comparison between radiochromic film measurements, Monte Carlo simulations, and collapsed cone convolution calculations

Head-and-neck tumors are often situated at an air-tissue interface what may result in an underdosage of part of the tumor in radiotherapy treatments using megavoltage photons, especially for small fields. In addition to effects of transient electronic disequilibrium, for these small fields, an increased lateral electron range in air will result in an important extra reduction of the central axis dose beyond the cavity. Therefore dose calculation algorithms need to model electron transport accurately. We simulated the trachea by a 2 cm diameter cylindrical air cavity with the rim situated 2 cm beneath the phantom surface. A 6 MV photon beam from an Elekta SLiplus linear accelerator, equipped with the standard multileaf collimator (MLC), was assessed. A 10 x 2 cm2 and a 10 x 1 cm2 field, both widthwise collimated by the MLC, were applied with their long side parallel to the cylinder axis. Central axis dose rebuild-up was studied. Radiochromic film measurements were performed in an in-house manufactured polystyrene phantom with the films oriented either along or perpendicular to the beam axis. Monte Carlo simulations were performed with BEAM and EGSnrc. Calculations were also performed using the pencil beam (PB) algorithm and the collapsed cone convolution (CCC) algorithm of Helax-TMS (MDS Nordion, Kanata, Cahada) version 6.0.2 and using the CCC algorithm of Pinnacle (ADAC Laboratories, Milpitas, CA, USA) version 4.2. A very good agreement between the film measurements and the Monte Carlo simulations was found. The CCC algorithms were not able to predict the interface dose accurately when lateral electronic disequilibrium occurs, but were shown to be a considerable improvement compared to the PB algorithm. The CCC algorithms overestimate the dose in the rebuild-up region. The interface dose was overestimated by a maximum of 31% or 54%, depending on the implementation of the CCC algorithm. At a depth of 1 mm, the maximum dose overestimation was 14% or 24%.

Monte Carlo evaluation of the AAA treatment planning algorithm in a heterogeneous multilayer phantom and IMRT clinical treatments for an Elekta SL25 linear accelerator

The Anisotropic Analytical Algorithm (AAA) is a new pencil beam convolution/superposition algorithm proposed by Varian for photon dose calculations. The configuration of AAA depends on linear accelerator design and specifications. The purpose of this study was to investigate the accuracy of AAA for an Elekta SL25 linear accelerator for small fields and intensity modulated radiation therapy (IMRT) treatments in inhomogeneous media. The accuracy of AAA was evaluated in two studies. First, AAA was compared both with Monte Carlo (MC) and the measurements in an inhomogeneous phantom simulating lung equivalent tissues and bone ribs. The algorithm was tested under lateral electronic disequilibrium conditions, using small fields (2 x 2 cm(2)). Good agreement was generally achieved for depth dose and profiles, with deviations generally below 3% in lung inhomogeneities and below 5% at interfaces. However, the effects of attenuation and scattering close to the bone ribs were not fully taken into account by AAA, and small inhomogeneities may lead to planning errors. Second, AAA and MC were compared for IMRT plans in clinical conditions, i.e., dose calculations in a computed tomography scan of a patient. One ethmoid tumor, one orophaxynx and two lung tumors are presented in this paper. Small differences were found between the dose volume histograms. For instance, a 1.7% difference for the mean planning target volume dose was obtained for the ethmoid case. Since better agreement was achieved for the same plans but in homogeneous conditions, these differences must be attributed to the handling of inhomogeneities by AAA. Therefore, inherent assumptions of the algorithm, principally the assumption of independent depth and lateral directions in the scaling of the kernels, were slightly influencing AAAs validity in inhomogeneities. However, AAA showed a good accuracy overall and a great ability to handle small fields in inhomogeneous media compared to other pencil beam convolution algorithms.

Conversion of CT numbers into tissue parameters for Monte Carlo dose calculations: a multi-centre study.

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

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.

Parameter dependence of the MCNP electron transport in determining dose distributions

In this paper, a detailed study of the electron transport in MCNP is performed, separating the effects of the energy binning technique on the energy loss rate, the scattering angles, and the sub-step length as a function of energy. As this problem is already well known, in this paper we focus on the explanation as to why the default mode of MCNP can lead to large deviations. The resolution dependence was investigated as well. An error in the MCNP code in the energy binning technique in the default mode (DBCN 18 card = 0) was revealed, more specific in the updating of cross sections when a sub-step is performed corresponding to a high-energy loss. This updating error is not present in the ITS mode (DBCN 18 card = 1) and leads to a systematically lower dose deposition rate in the default mode. The effect is present for all energies studied (0.5-10 MeV) and depends on the geometrical resolution of the scoring regions and the energy grid resolution. The effect of the energy binning technique is of the same order of that of the updating error for energies below 2 MeV, and becomes less important for higher energies. For a 1 MeV point source surrounded by homogeneous water, the deviation of the default MCNP results at short distances attains 9% and remains approximately the same for all energies. This effect could be corrected by removing the completion of an energy step each time an electron changes from an energy bin during a sub-step. Another solution consists of performing all calculations in the ITS mode. Another problem is the resolution dependence, even in the ITS mode. The higher the resolution is chosen (the smaller the scoring regions) the faster the energy is deposited along the electron track. It is proven that this is caused by starting a new energy step when crossing a surface. The resolution effect should be investigated for every specific case when calculating dose distributions around beta sources. The resolution should not be higher than 0.85*(1-EFAC)*CSDA, where EFAC is the energy loss per energy step and CSDA a continuous slowing down approximation range. This effect could as well be removed by determining the cross sections for energy loss and multiple scattering at the average energy of an energy step and by sampling the cross sections for each sub-step. Overall, we conclude that MCNP cannot be used without a caution due to possible errors in the electron transport. When care is taken, it is possible to obtain correct results that are in agreement with other Monte Carlo codes.

Dosimetric comparison of different treatment modalities for stereotactic radiosurgery of arteriovenous malformations and acoustic neuromas

PURPOSE We investigated the influence of beam modulation on treatment planning by comparing four available stereotactic radiosurgery (SRS) modalities: Gamma-Knife-Perfexion, Novalis-Tx Dynamic-Conformal-Arc (DCA) and Dynamic-Multileaf-Collimation-Intensity-Modulated-radiotherapy (DMLC-IMRT), and Cyberknife. MATERIAL AND METHODS Patients with arteriovenous malformation (n = 10) or acoustic neuromas (n = 5) were planned with different treatment modalities. Paddick conformity index (CI), dose heterogeneity (DH), gradient index (GI) and beam-on time were used as dosimetric indices. RESULTS Gamma-Knife-Perfexion can achieve high degree of conformity (CI = 0.77 ± 0.04) with limited low-doses (GI = 2.59 ± 0.10) surrounding the inhomogeneous dose distribution (D(H) = 0.84 ± 0.05) at the cost of treatment time (68.1 min ± 27.5). Novalis-Tx-DCA improved this inhomogeneity (D(H) = 0.30 ± 0.03) and treatment time (16.8 min ± 2.2) at the cost of conformity (CI = 0.66 ± 0.04) and Novalis-TX-DMLC-IMRT improved the DCA CI (CI = 0.68 ± 0.04) and inhomogeneity (D(H) = 0.18 ± 0.05) at the cost of low-doses (GI = 3.94 ± 0.92) and treatment time (21.7 min ± 3.4) (p<0.01). Cyberknife achieved comparable conformity (CI = 0.77 ± 0.06) at the cost of low-doses (GI = 3.48 ± 0.47) surrounding the homogeneous (D(H) = 0.22 ± 0.02) dose distribution and treatment time (28.4min±8.1) (p<0.01). CONCLUSIONS Gamma-Knife-Perfexion will comply with all SRS constraints (high conformity while minimizing low-dose spread). Multiple focal entries (Gamma-Knife-Perfexion and Cyberknife) will achieve better conformity than High-Definition-MLC of Novalis-Tx at the cost of treatment time. Non-isocentric beams (Cyberknife) or IMRT-beams (Novalis-Tx-DMLC-IMRT) will spread more low-dose than multiple isocenters (Gamma-Knife-Perfexion) or dynamic arcs (Novalis-Tx-DCA). Inverse planning and modulated fluences (Novalis-Tx-DMLC-IMRT and CyberKnife) will deliver the most homogeneous treatment. Furthermore, Linac-based systems (Novalis and Cyberknife) can perform image verification at the time of treatment delivery.

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.