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

Report of AAPM Task Group 155: Megavoltage photon beam dosimetry in small fields and non-equilibrium conditions

Small-field dosimetry used in advance treatment technologies poses challenges due to loss of lateral charged particle equilibrium (LCPE), occlusion of the primary photon source, and the limited choice of suitable radiation detectors. These challenges greatly influence dosimetric accuracy. Many high-profile radiation incidents have demonstrated a poor understanding of appropriate methodology for small-field dosimetry. These incidents are a cause for concern because the use of small fields in various specialized radiation treatment techniques continues to grow rapidly. Reference and relative dosimetry in small and composite fields are the subject of the International Atomic Energy Agency (IAEA) dosimetry code of practice that has been published as TRS-483 and an AAPM summary publication (IAEA TRS 483; Dosimetry of small static fields used in external beam radiotherapy: An IAEA/AAPM International Code of Practice for reference and relative dose determination, Technical Report Series No. 483; Palmans et al., Med Phys 45(11):e1123, 2018). The charge of AAPM task group 155 (TG-155) is to summarize current knowledge on small-field dosimetry and to provide recommendations of best practices for relative dose determination in small megavoltage photon beams. An overview of the issue of LCPE and the changes in photon beam perturbations with decreasing field size is provided. Recommendations are included on appropriate detector systems and measurement methodologies. Existing published data on dosimetric parameters in small photon fields (e.g., percentage depth dose, tissue phantom ratio/tissue maximum ratio, off-axis ratios, and field output factors) together with the necessary perturbation corrections for various detectors are reviewed. A discussion on errors and an uncertainty analysis in measurements is provided. The design of beam models in treatment planning systems to simulate small fields necessitates special attention on the influence of the primary beam source and collimating devices in the computation of energy fluence and dose. The general requirements for fluence and dose calculation engines suitable for modeling dose in small fields are reviewed. Implementations in commercial treatment planning systems vary widely, and the aims of this report are to provide insight for the medical physicist and guidance to developers of beams models for radiotherapy treatment planning systems.

Dosimetric verification of lung phantom calculated by collapsed cone convolution: A Monte Carlo and experimental evaluation

OBJECTIVE

To evaluate the dose calculation accuracy in the Prowess Panther treatment planning system (TPS) using the collapsed cone convolution (CCC) algorithm.

METHODS

The BEAMnrc Monte Carlo (MC) package was used to predict the dose distribution of photon beams produced by the Oncor® linear accelerator (linac). The MC model of an 18 MV photon beam was verified by measurement using a p-type diode dosimeter. Percent depth dose (PDD) and dose profiles were used for comparison based on three field sizes: 5×5, 10×10, and 20×20cm2. The accuracy of the CCC dosimetry was also evaluated using a plan composed of a simple parallel-opposed field (11×16cm2) in a lung phantom comprised of four tissue simulating media namely, lung, soft tissue, bone and spinal cord. The CCC dose calculation accuracy was evaluated by MC simulation and measurements according to the dose difference and 3D gamma analysis. Gamma analysis was carried out through comparison of the Monte Carlo simulation and the TPS calculated dose.

RESULTS

Compared to the dosimetric results measured by the Farmer chamber, the CCC algorithm underestimated dose in the planning target volume (PTV), right lung and lung-tissue interface regions by about -0.11%, -1.6 %, and -2.9%, respectively. Moreover, the CCC algorithm underestimated the dose at the PTV, right lung and lung-tissue interface regions in the order of -0.34%, -0.4% and -3.5%, respectively, when compared to the MC simulation. Gamma analysis results showed that the passing rates within the PTV and heterogeneous region were above 59% and 76%. For the right lung and spinal cord, the passing rates were above 80% for all gamma criteria.

CONCLUSIONS

This study demonstrates that the CCC algorithm has potential to calculate dose with sufficient accuracy for 3D conformal radiotherapy within the thorax where a significant amount of tissue heterogeneity exists.

Experimental evaluation of a GPU-based Monte Carlo dose calculation algorithm in the Monaco treatment planning system

A new GPU-based Monte Carlo dose calculation algorithm (GPUMCD), devel-oped by the vendor Elekta for the Monaco treatment planning system (TPS), is capable of modeling dose for both a standard linear accelerator and an Elekta MRI linear accelerator. We have experimentally evaluated this algorithm for a standard Elekta Agility linear accelerator. A beam model was developed in the Monaco TPS (research version 5.09.06) using the commissioned beam data for a 6 MV Agility linac. A heterogeneous phantom representing several scenarios - tumor-in-lung, lung, and bone-in-tissue - was designed and built. Dose calculations in Monaco were done using both the current clinical Monte Carlo algorithm, XVMC, and the new GPUMCD algorithm. Dose calculations in a Pinnacle TPS were also produced using the collapsed cone convolution (CCC) algorithm with heterogeneity correc-tion. Calculations were compared with the measured doses using an ionization chamber (A1SL) and Gafchromic EBT3 films for 2 × 2 cm2, 5 × 5 cm2, and 10 × 10 cm2 field sizes. The percentage depth doses (PDDs) calculated by XVMC and GPUMCD in a homogeneous solid water phantom were within 2%/2 mm of film measurements and within 1% of ion chamber measurements. For the tumor-in-lung phantom, the calculated doses were within 2.5%/2.5 mm of film measurements for GPUMCD. For the lung phantom, doses calculated by all of the algorithms were within 3%/3 mm of film measurements, except for the 2 × 2 cm2 field size where the CCC algorithm underestimated the depth dose by ~ 5% in a larger extent of the lung region. For the bone phantom, all of the algorithms were equivalent and calculated dose to within 2%/2 mm of film measurements, except at the interfaces. Both GPUMCD and XVMC showed interface effects, which were more pronounced for GPUMCD and were comparable to film measurements, whereas the CCC algorithm showed these effects poorly.

New approach in lung cancer radiotherapy offers better normal tissue sparing

PURPOSE

Medical images are more than pictures. They contain additional quantitative information which can be interrogated, quantified, and utilized. Besides anatomical information computed tomography (CT) imaging data provide electron density information. Radiotherapy use of this density information is limited to its application only in dose calculations. The direct product of dose, density, and volume forms a quantity called integral dose. The integral dose delivered to a volume of interest is the total energy deposited in that volume. Here it is hypothesized that minimization of the integral dose is advantageous in radiotherapy planning. The purpose of this work is to study the incorporation of quantitative imaging information in radiotherapy inverse optimization through total energy minimization (Energy hereafter).

DESIGN

Twenty lung patient plans were studied. For each patient density was quantified on voxel-by-voxel basis through image gray value-to-density conversion curves. Energy-based objective function was used for inverse radiotherapy plan optimization. The obtained plans were evaluated in the light of current standard of care, based on dose-volume (DVH) optimization approach.

RESULTS

The statistical significance analyses of the results indicated that the doses to normal tissue were between 14% and 45% lower, when Energy-based optimization was used instead of DVH-based optimization.

CONCLUSION

Incorporation of quantitative imaging information, through CT derived density, in the optimization cost function allows reduction of dose to normal tissue for NSCLC cases. Energy-based radiotherapy plans result in lower normal tissue dose and potentially lower complication rates compared to standard of care.

Monte Carlo as a tool to evaluate the effect of different lung densities on radiotherapy dose distribution

This study aims at evaluating the effects of different lung densities on dose distribution after irradiation at different field sizes, by comparing experimental measurements, GEANT4 Monte Carlo (MC) simulations and two TPS calculation algorithms on ad hoc phantoms. Irradiations were performed with a Varian Clinac 2100 C/D with a nominal energy of 6 MV. Dosimetric experimental measurements were obtained with radiochromic films. A model based on GEANT4 MC code was developed to simulate both the accelerator and the phantoms. Results of dose distribution show an acceptable agreement between MC simulations and experimental measurements, both in the tumour-equivalent region and in the normal tissue-equivalent ones. On the opposite, results vary among the TPS algorithms, especially in regions of lung-equivalent material at low density, but also at the interface between lung- and tumour-equivalent materials.

An institutional experience of quality assurance of a treatment planning system on photon beam

AIM

The purpose of the present study is to show the application of the IAEA TRS-430 QA procedures of Eclipse™v7.5 TPS for photon energies. In addition, the trends of the deviations found in the conducted tests were determined.

BACKGROUND

In the past, the lack of complete TPS QA procedures led to some serious accidents. So, QA in the radiotherapy treatment planning process is essential for determination of accuracy in the radiotherapy process and avoidance of treatment errors.

MATERIALS AND METHODS

The calculations of TPS and measurements of irradiations of the treatment device were compared in the study. As a result, the local dose deviation values (δ 1: central beam axis, δ 2: penumbra and build up region, δ 3: inside field, δ 4: outside beam edges, δ 50-90: beam fringe, RW50: radiological width) and their confidence limit values (including systematic and random errors) were obtained.

RESULTS

The confidence limit values of δ 4 were detected to increase with expanding field size. The values of δ 1 and δ 3 of hard wedge were larger than open fields. The values of δ 2 and δ 50-90 of the inhomogeneity effect test were larger, especially than other tests of this study. The average deviation was showed to increase with the rise of the wedge angle. The values of δ 3 and δ 4 of lung irradiation were outside tolerance.

CONCLUSIONS

The QA of TPS was done and it was found that there were no reservations in its use in patient treatment. The trend of the deviations is shown.

Independent absorbed-dose calculation using the Monte Carlo algorithm in volumetric modulated arc therapy

Purpose

To report the result of independent absorbed-dose calculations based on a Monte Carlo (MC) algorithm in volumetric modulated arc therapy (VMAT) for various treatment sites.

Methods and materials

All treatment plans were created by the superposition/convolution (SC) algorithm of SmartArc (Pinnacle V9.2, Philips). The beam information was converted into the format of the Monaco V3.3 (Elekta), which uses the X-ray voxel-based MC (XVMC) algorithm. The dose distribution was independently recalculated in the Monaco. The dose for the planning target volume (PTV) and the organ at risk (OAR) were analyzed via comparisons with those of the treatment plan.

Before performing an independent absorbed-dose calculation, the validation was conducted via irradiation from 3 different gantry angles with a 10- × 10-cm² field. For the independent absorbed-dose calculation, 15 patients with cancer (prostate, 5; lung, 5; head and neck, 3; rectal, 1; and esophageal, 1) who were treated with single-arc VMAT were selected. To classify the cause of the dose difference between the Pinnacle and Monaco TPSs, their calculations were also compared with the measurement data.

Result

In validation, the dose in Pinnacle agreed with that in Monaco within 1.5%. The agreement in VMAT calculations between Pinnacle and Monaco using phantoms was exceptional; at the isocenter, the difference was less than 1.5% for all the patients. For independent absorbed-dose calculations, the agreement was also extremely good. For the mean dose for the PTV in particular, the agreement was within 2.0% in all the patients; specifically, no large difference was observed for high-dose regions. Conversely, a significant difference was observed in the mean dose for the OAR. For patients with prostate cancer, the mean rectal dose calculated in Monaco was significantly smaller than that calculated in Pinnacle.

Conclusions

There was no remarkable difference between the SC and XVMC calculations in the high-dose regions. The difference observed in the low-dose regions may have arisen from various causes such as the intrinsic dose deviation in the MC calculation, modeling accuracy, and CT-to-density table used in each planning system It is useful to perform independent absorbed-dose calculations with the MC algorithm in intensity-modulated radiation therapy commissioning.

Validation of GPU based TomoTherapy dose calculation engine

PURPOSE

The graphic processing unit (GPU) based TomoTherapy convolution/superposition(C/S) dose engine (GPU dose engine) achieves a dramatic performance improvement over the traditional CPU-cluster based TomoTherapy dose engine (CPU dose engine). Besides the architecture difference between the GPU and CPU, there are several algorithm changes from the CPU dose engine to the GPU dose engine. These changes made the GPU dose slightly different from the CPU-cluster dose. In order for the commercial release of the GPU dose engine, its accuracy has to be validated.

METHODS

Thirty eight TomoTherapy phantom plans and 19 patient plans were calculated with both dose engines to evaluate the equivalency between the two dose engines. Gamma indices (Γ) were used for the equivalency evaluation. The GPU dose was further verified with the absolute point dose measurement with ion chamber and film measurements for phantom plans. Monte Carlo calculation was used as a reference for both dose engines in the accuracy evaluation in heterogeneous phantom and actual patients.

RESULTS

The GPU dose engine showed excellent agreement with the current CPU dose engine. The majority of cases had over 99.99% of voxels with Γ(1%, 1 mm) < 1. The worst case observed in the phantom had 0.22% voxels violating the criterion. In patient cases, the worst percentage of voxels violating the criterion was 0.57%. For absolute point dose verification, all cases agreed with measurement to within ±3% with average error magnitude within 1%. All cases passed the acceptance criterion that more than 95% of the pixels have Γ(3%, 3 mm) < 1 in film measurement, and the average passing pixel percentage is 98.5%-99%. The GPU dose engine also showed similar degree of accuracy in heterogeneous media as the current TomoTherapy dose engine.

CONCLUSIONS

It is verified and validated that the ultrafast TomoTherapy GPU dose engine can safely replace the existing TomoTherapy cluster based dose engine without degradation in dose accuracy.

Influence of dose calculation algorithms on the predicted dose distribution and NTCP values for NSCLC patients

PURPOSE

To investigate differences in calculated doses and normal tissue complication probability (NTCP) values between different dose algorithms.

METHODS

Six dose algorithms from four different treatment planning systems were investigated: Eclipse AAA, Oncentra MasterPlan Collapsed Cone and Pencil Beam, Pinnacle Collapsed Cone and XiO Multigrid Superposition, and Fast Fourier Transform Convolution. Twenty NSCLC patients treated in the period 2001-2006 at the same accelerator were included and the accelerator used for treatments were modeled in the different systems. The treatment plans were recalculated with the same number of monitor units and beam arrangements across the dose algorithms. Dose volume histograms of the GTV, PTV, combined lungs (excluding the GTV), and heart were exported and evaluated. NTCP values for heart and lungs were calculated using the relative seriality model and the LKB model, respectively. Furthermore, NTCP for the lungs were calculated from two different model parameter sets. Calculations and evaluations were performed both including and excluding density corrections.

RESULTS

There are found statistical significant differences between the calculated dose to heart, lung, and targets across the algorithms. Mean lung dose and V20 are not very sensitive to change between the investigated dose calculation algorithms. However, the different dose levels for the PTV averaged over the patient population are varying up to 11%. The predicted NTCP values for pneumonitis vary between 0.20 and 0.24 or 0.35 and 0.48 across the investigated dose algorithms depending on the chosen model parameter set. The influence of the use of density correction in the dose calculation on the predicted NTCP values depends on the specific dose calculation algorithm and the model parameter set. For fixed values of these, the changes in NTCP can be up to 45%.

CONCLUSIONS

Calculated NTCP values for pneumonitis are more sensitive to the choice of algorithm than mean lung dose and V20 which are also commonly used for plan evaluation. The NTCP values for heart complication are, in this study, not very sensitive to the choice of algorithm. Dose calculations based on density corrections result in quite different NTCP values than calculations without density corrections. It is therefore important when working with NTCP planning to use NTCP parameter values based on calculations and treatments similar to those for which the NTCP is of interest.

Dosimetry of oblique tangential photon beams calculated by superposition/convolution algorithms: a Monte Carlo evaluation

Although there are many works on evaluating dose calculations of the anisotropic analytical algorithm (AAA) using various homogeneous and heterogeneous phantoms, related work concerning dosimetry due to tangential photon beam is lacking. In this study, dosimetry predicted by the AAA and collapsed cone convolution (CCC) algorithm was evaluated using the tangential photon beam and phantom geometry. The photon beams of 6 and 15 MV with field sizes of 4 × 4 (or 7 × 7), 10 × 10 and 20 × 20 cm², produced by a Varian 21 EX linear accelerator, were used to test performances of the AAA and CCC using Monte Carlo (MC) simulation (EGSnrc-based code) as a benchmark. Horizontal dose profiles at different depths, phantom skin profiles (i.e., vertical dose profiles at a distance of 2 mm from the phantom lateral surface), gamma dose distributions, and dose-volume histograms (DVHs) of skin slab were determined. For dose profiles at different depths, the CCC agreed better with doses in the air-phantom region, while both the AAA and CCC agreed well with doses in the penumbra region, when compared to the MC. Gamma evaluations between the AAA/CCC and MC showed that deviations of 2D dose distribution occurred in both beam edges in the phantom and air-phantom interface. Moreover, the gamma dose deviation is less significant in the air-phantom interface than the penumbra. DVHs of skin slab showed that both the AAA and CCC underestimated the width of the dose drop-off region for both the 6 and 15 MV photon beams. When the gantry angle was 0°, it was found that both the AAA and CCC overestimated doses in the phantom skin profiles compared to the MC, with various photon beam energies and field sizes. The mean dose differences with doses normalized to the prescription point for the AAA and CCC were respectively: 7.6% ± 2.6% and 2.1% ± 1.3% for a 10 × 10 cm2 field, 6 MV; 16.3%± 2.1% and 6.7% ± 2.1% for a 20 × 20 cm2 field, 6 MV; 5.5% ± 1.2% and 1.7% ± 1.4% for a 10 × 10 cm2, 15 MV; 18.0% ± 1.3% and 8.3% ± 1.8% for a 20 × 20 cm², 15 MV. However, underestimations of doses in the phantom skin profile were found with small fields of 4 × 4 and 7 × 7 cm² for the 6 and 15 MV photon beams, respectively, when the gantry was turned 5° anticlockwise. As surface dose with tangential photon beam geometry is important in some radiation treatment sites such as breast, chest wall and sarcoma, it is found that neither of the treatment planning system algorithms can predict the dose well at depths shallower than 2 mm. The dosimetry data and beam and phantom geometry in this study provide a better knowledge of a dose calculation algorithm in tangential-like irradiation.

Lung dose for minimally moving thoracic lesions treated with respiration gating

PURPOSE

To evaluate incidental doses to benign lung tissue for patients with minimally moving lung lesions treated with respiratory gating.

METHODS AND MATERIALS

Seventeen lung patient plans were studied retrospectively. Tumor motion was less than 5 mm in all cases. For each patient, mid-ventilation (MidVen) and mid-inhalation (MidInh) breathing phases were reconstructed. The MidInh phase was centered on the end-of-inhale (EOI) phase within a 30% gating window. Planning target volumes, heart, and spinal cord were delineated on the MidVen phase and transferred to the MidInh phase. Lungs were contoured separately on each phase. Intensity-modulated radiotherapy plans were generated on the MidVen phases. The plans were transferred to the MidInh phase, and doses were recomputed. The evaluation metric was based on dose indices, volume indices, generalized equivalent uniform doses, and mass indices for targets and critical structures. Statistical tests were used to establish the significance of the differences between the reference (MidVen) and compared (MidInh) dose distributions.

RESULTS

Statistical tests demonstrated that the indices evaluated for targets, cord, and heart differed by within 2.3%. The index differences in the lungs, however, are in excess of 6%, indicating the potentially achievable lung sparing and/or dose escalation.

CONCLUSIONS

Respiratory gating is a clinical option for patients with minimally moving lung lesions treated at EOI. Gating will be more beneficial for larger tumors, since dose escalation in those cases will result in a larger increase in the tumor control probability.

Determination of the optimal statistical uncertainty to perform electron-beam Monte Carlo absorbed dose estimation in the target volume

PURPOSE OF STUDY

Monte Carlo based treatment planning system are known to be more accurate than analytical methods for performing absorbed dose estimation, particularly in and near heterogeneities. However, the required computation time can still be an issue. The present study focused on the determination of the optimum statistical uncertainty in order to minimise computation time while keeping the reliability of the absorbed dose estimation in treatments planned with electron-beams.

MATERIALS AND METHODS

Three radiotherapy plans (medulloblastoma, breast and gynaecological) were used to investigate the influence of the statistical uncertainty of the absorbed dose on the target volume dose-volume histograms (spinal cord, intramammary nodes and pelvic lymph nodes, respectively).

RESULTS

The study of the dose-volume histograms showed that for statistical uncertainty levels (1 S.D.) above 2 to 3%, the standard deviation of the mean dose in the target volume calculated from the dose-volume histograms increases by at least 6%, reflecting the gradual flattening of the dose-volume histograms.

CONCLUSIONS

This work suggests that, in clinical context, Monte Carlo based absorbed dose estimations should be performed with a maximum statistical uncertainty of 2 to 3%.

Monte Carlo vs. pencil beam based optimization of stereotactic lung IMRT

Background

The purpose of the present study is to compare finite size pencil beam (fsPB) and Monte Carlo (MC) based optimization of lung intensity-modulated stereotactic radiotherapy (lung IMSRT).

Materials and methods

A fsPB and a MC algorithm as implemented in a biological IMRT planning system were validated by film measurements in a static lung phantom. Then, they were applied for static lung IMSRT planning based on three different geometrical patient models (one phase static CT, density overwrite one phase static CT, average CT) of the same patient. Both 6 and 15 MV beam energies were used. The resulting treatment plans were compared by how well they fulfilled the prescribed optimization constraints both for the dose distributions calculated on the static patient models and for the accumulated dose, recalculated with MC on each of 8 CTs of a 4DCT set.

Results

In the phantom measurements, the MC dose engine showed discrepancies < 2%, while the fsPB dose engine showed discrepancies of up to 8% in the presence of lateral electron disequilibrium in the target. In the patient plan optimization, this translates into violations of organ at risk constraints and unpredictable target doses for the fsPB optimized plans. For the 4D MC recalculated dose distribution, MC optimized plans always underestimate the target doses, but the organ at risk doses were comparable. The results depend on the static patient model, and the smallest discrepancy was found for the MC optimized plan on the density overwrite one phase static CT model.

Conclusions

It is feasible to employ the MC dose engine for optimization of lung IMSRT and the plans are superior to fsPB. Use of static patient models introduces a bias in the MC dose distribution compared to the 4D MC recalculated dose, but this bias is predictable and therefore MC based optimization on static patient models is considered safe.

Endo-rectal balloon cavity dosimetry in a phantom: performance under IMRT and helical tomotherapy beams

BACKGROUND AND PURPOSE

The use of endo-rectal balloons as immobilisation devices in external beam radiotherapy for prostate cancer has led to improved target position reproducibility and a decrease in rectal toxicity. The air cavity created by an endo-rectal balloon in photon radiotherapy perturbs the dose distribution. In this study, the effect of the balloon cavity on the dose distribution and the accuracy to which two treatment planning systems calculate the dose distribution were investigated.

MATERIALS AND METHODS

Single beams as well as 3D conformal, conventional IMRT and helical tomotherapy treatment plans were investigated using a specifically constructed phantom. Radiochromic film was used to measure the cavity wall doses and cavity wall DVHs.

RESULTS

For a 70 Gy prescription dose both the Pinnacle and TomoTherapy TPSs over-predicted the anterior cavity wall dose by 1.43 Gy, 3.92 Gy and 2.67 Gy for 3D conformal, conventional IMRT and helical tomotherapy, respectively. The posterior cavity wall dose was under-predicted by 2.62 Gy, 2.01 Gy and 4.79 Gy for 3D conformal, conventional IMRT and helical tomotherapy, respectively. An over-prediction by the Pinnacle RTPS of the V50, V60, V65 and V70 values for the cavity wall DVH was measured for the 3D conformal and conventional IMRT cases. These reductions may lead to a less than expected rectal toxicity. The TomoTherapy RTPS under-predicted the V50, V60, V65 and V70 values which may lead to higher rectal toxicity than predicted.

CONCLUSION

Calculation of dose around an air cavity created by an endo-rectal balloon provides a challenge for radiotherapy planning systems. Various electronic disequilibrium situations exist due to the cavity, which can lead to a lower anterior rectal wall and higher posterior rectal wall dose than that calculated by planning systems. This has consequences for comparisons of dose volume constraints between different modalities.

Commissioning a fast Monte Carlo dose calculation algorithm for lung cancer treatment planning

A commercial Monte Carlo simulation package, NXEGS 1.12 (NumeriX LLC, New York, NY), was commissioned for photon‐beam dose calculations. The same sets of measured data from 6‐MV and 18‐MV beams were used to commission NXEGS and Pinnacle 6.2b (Philips Medical Systems, Andover, MA). Accuracy and efficiency were compared against the collapsed cone convolution algorithm implemented in Pinnacle 6.2b, together with BEAM simulation (BEAMnrc 2001: National Research Council of Canada, Ottawa, ON). We investigated a number of options in NXEGS: the accuracy of fast Monte Carlo, the re‐implementation of EGS4, post‐processing technique (dose de‐noising algorithm), and dose calculation time. Dose distributions were calculated with NXEGS, Pinnacle, and BEAM in water, lung‐slab, and air‐cylinder phantoms and in a lung patient plan. We compared the dose distributions calculated by NXEGS, Pinnacle, and BEAM. In a selected region of interest (7725 voxels) in the lung phantom, all but 1 voxel had a γ (3% and 3 mm thresholds) of 1 or less for the dose difference between the NXEGS re‐implementation of EGS4 and BEAM, and 99% of the voxels had a γ of 1 or less for the dose difference between NXEGS fast Monte Carlo and BEAM. Fast Monte Carlo with post‐processing was up to 100 times faster than the NXEGS re‐implementation of EGS4, while maintaining ±2% statistical uncertainty. With air inhomogeneities larger than 1 cm, post‐processing preserves the dose perturbations from the air cylinder. When 3 or more beams were used, fast Monte Carlo with post‐processing was comparable to or faster than Pinnacle 6.2b collapsed cone convolution.

PACS numbers: 87.18.Bb, 87.53.Wz

An evaluation of planning techniques for stereotactic body radiation therapy in lung tumors

PURPOSE

To evaluate four planning techniques for stereotactic body radiation therapy (SBRT) in lung tumors.

METHODS AND MATERIALS

Four SBRT plans were performed for 12 patients with stage I/II non-small-cell lung cancer under the following conditions: (1) conventional margins on free-breathing CT (plan 1), (2) generation of an internal target volume (ITV) using 4DCT with beam delivery under free-breathing conditions (plan 2), (3) gating at end-exhale (plan 3), and (4) gating at end-inhale (plan 4). Planning was performed following the RTOG 0236 protocol with a prescription dose of 54 Gy (3 fractions). For each plan 4D dose was calculated using deformable-image registration.

RESULTS

There was no significant difference in tumor dose delivered by the 4 plans. However, compared with plan 1, plans 2-4 reduced total lung BED by 1.9+/-1.2, 3.1+/-1.6 and 3.5+/-2.1 Gy, reduced mean lung dose by 0.8+/-0.5, 1.5+/-0.8, and 1.6+/-1.0 Gy, reduced V20 by 1.5+/-1.0%, 2.7+/-1.4%, and 2.8+/-1.8%, respectively, with p<0.01. Compared with plan 2, plans 3-4 reduced lung BED by 1.2+/-1.0 and 1.6+/-1.5 Gy, reduced mean lung dose by 0.6+/-0.5 and 0.8+/-0.7 Gy, reduced V20 by 1.2+/-1.1% and 1.3+/-1.5%, respectively, with p<0.01. The differences in lung BED, mean dose and V20 of plan 4 compared with plan 3 were insignificant.

CONCLUSIONS

Tumor dose coverage was statistically insignificant between all plans. However, compared with plan 1, plans 2-4 significantly reduced lung doses. Compared with plan 2, plan 3-4 also reduced lung toxicity. The difference in lung doses between plan 3 and plan 4 was not significant.

Hybrid IMRT for treatment of cancers of the lung and esophagus

PURPOSE

To report on a hybrid intensity-modulated radiation therapy (IMRT; static plus IMRT beams treated concurrently) technique for lung and esophageal patients to reduce the volume of lung treated to low doses while delivering a conformal dose distribution.

METHODS

Treatment plans were analyzed for 18 patients (12 lung and 6 esophageal). Patients were treated with a hybrid technique that concurrently combines static (approximately two-thirds dose) and IMRT (approximately one-third dose) beams. These plans were compared with conventional three-dimensional (3D; non-IMRT) plans and all IMRT plans using custom four- and five-field arrangements and nine equally spaced coplanar beams. Plans were optimized to reduce V13 and V5 values. Dose-volume histograms were calculated for the planning target volume, heart, and the ipsilateral, contralateral, and total lung. Lung volumes V5, V13, V20, V30; mean lung dose (MLD); and the generalized equivalent uniform dose (gEUD) were calculated for each plan.

RESULTS

Hybrid plans treated significantly smaller total and contralateral lung volumes with low doses than nine-field IMRT plans. Largest reductions were for contralateral lung V5, V13, and V20 values for lung (-11%, -15%, -7%) and esophageal (-16%, -20%, -7%) patients. Smaller reductions were found also for 3D and four- and five-field IMRT plans. MLD and gEUDs were similar for all plan types. The 3D plans treated much larger extra planning target volumes to prescribed dose levels.

CONCLUSIONS

Hybrid IMRT demonstrated advantages for reduction of low-dose lung volumes in the thorax for reducing low dose to lung while also reducing the potential magnitude of dose deviations due to intrafraction motion and small field calculation accuracy.

Multibeam tomotherapy: a new treatment unit devised for multileaf collimation, intensity-modulated radiation therapy

A fully integrated system for treatment planning, application, and verification for automated multileaf collimator (MLC) based, intensity-modulated, image-guided, and adaptive radiation therapy (IMRT, IGRT and ART, respectively) is proposed. Patient comfort, which was the major development goal, will be achieved through a new unit design and short treatment times. Our device for photon beam therapy will consist of a new dual energy linac with five fixed treatment heads positioned evenly along one plane but one electron beam generator only. A minimum of moving parts increases technical reliability and reduces motion times to a minimum. Motion is allowed solely for the MLCs, the robotic patient table, and the small angle gantry rotation of +/- 36 degrees. Besides sophisticated electron beam guidance, this compact setup can be built using existing modules. The flattening-filter-free treatment heads are characterized by reduced beam-on time and contain apertures restricted in one dimension to the area of maximum primary fluence output. In the case of longer targets, this leads to a topographic intensity modulation, thanks to the combination of "step and shoot" MLC delivery and discrete patient couch motion. Owing to the limited number of beam directions, this multislice cone beam serial tomotherapy is referred to as "multibeam tomotherapy." Every patient slice is irradiated by one treatment head at any given moment but for one subfield only. The electron beam is then guided to the next head ready for delivery, while the other heads are preparing their leaves for the next segment. The "Multifocal MLC-positioning" algorithm was programmed to enable treatment planning and optimize treatment time. We developed an overlap strategy for the longitudinally adjacent fields of every beam direction, in doing so minimizing the field match problem and the effects of possible table step errors. Clinical case studies show for the same or better planning target volume coverage, better organ-at-risk sparing, and comparable mean integral dose to the normal tissue a reduction in treatment time by more than 50% to only a few minutes in comparison to high-quality 3-D conformal and IMRT treatments. As a result, it will be possible to incorporate features for better patient positioning and image guidance, while sustaining reasonable overall treatment times at the same time. The virtual multibeam tomotherapy design study TOM'5-CT contains a dedicated electron beam CT (TOM'AGE) and an objective optical topometric patient positioning system (TOPOS). Thanks to the wide gantry bore of 120 cm and slim gantry depths of 70 cm, patients can be treated very comfortably, in all cases tumor-isocentrically, as well as with noncoplanar beam arrangements as in stereotactic radiosurgery with a couch rotation of up to +/- 54 degrees. The TOM'5 treatment unit on which this theoretical concept is based has a stand-alone depth of 40 cm and an outer diameter of 245 cm; the focus-isocenter distance of the heads is 100 cm with a field size of 40 cm x 7 cm and 0.5 cm leaves, which operate perpendicular to the axis of table motion.

On the dosimetric behaviour of photon dose calculation algorithms in the presence of simple geometric heterogeneities: comparison with Monte Carlo calculations

A comparative study was performed to reveal differences and relative figures of merit of seven different calculation algorithms for photon beams when applied to inhomogeneous media. The following algorithms were investigated: Varian Eclipse: the anisotropic analytical algorithm, and the pencil beam with modified Batho correction; Nucletron Helax-TMS: the collapsed cone and the pencil beam with equivalent path length correction; CMS XiO: the multigrid superposition and the fast Fourier transform convolution; Philips Pinnacle: the collapsed cone. Monte Carlo simulations (MC) performed with the EGSnrc codes BEAMnrc and DOSxyznrc from NRCC in Ottawa were used as a benchmark. The study was carried out in simple geometrical water phantoms (rho = 1.00 g cm(-3)) with inserts of different densities simulating light lung tissue (rho = 0.035 g cm(-3)), normal lung (rho = 0.20 g cm(-3)) and cortical bone tissue (rho = 1.80 g cm(-3)). Experiments were performed for low- and high-energy photon beams (6 and 15 MV) and for square (13 x 13 cm2) and elongated rectangular (2.8 x 13 cm2) fields. Analysis was carried out on the basis of depth dose curves and transverse profiles at several depths. Assuming the MC data as reference, gamma index analysis was carried out distinguishing between regions inside the non-water inserts or inside the uniform water. For this study, a distance to agreement was set to 3 mm while the dose difference varied from 2% to 10%. In general all algorithms based on pencil-beam convolutions showed a systematic deficiency in managing the presence of heterogeneous media. In contrast, complicated patterns were observed for the advanced algorithms with significant discrepancies observed between algorithms in the lighter materials (rho = 0.035 g cm(-3)), enhanced for the most energetic beam. For denser, and more clinical, densities a better agreement among the sophisticated algorithms with respect to MC was observed.

Dosimetric verification of the anisotropic analytical algorithm for radiotherapy treatment planning

BACKGROUND AND PURPOSE

To investigate the accuracy of photon dose calculations performed by the Anisotropic Analytical Algorithm, in homogeneous and inhomogeneous media and in simulated treatment plans.

MATERIALS AND METHODS

Predicted dose distributions were compared with ionisation chamber and film measurements for a series of increasingly complex situations. Initially, simple and complex fields in a homogeneous medium were studied. The effect of inhomogeneities was investigated using a range of phantoms constructed of water, bone and lung substitute materials. Simulated treatment plans were then produced using a semi-anthropomorphic phantom and the delivered doses compared to the doses predicted by the Anisotropic Analytical Algorithm.

RESULTS

In a homogeneous medium, agreement was found to be within 2% dose or 2mm dta in most instances. In the presence of heterogeneities, agreement was generally to within 2.5%. The simulated treatment plan measurements agreed to within 2.5% or 2mm.

CONCLUSIONS

The accuracy of the algorithm was found to be satisfactory at 6 and 10MV both in homogeneous and inhomogeneous situations and in the simulated treatment plans. The algorithm was more accurate than the Pencil Beam Convolution model, particularly in the presence of low density heterogeneities.

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-TMS's collapsed cone model Helax-CC) and one pencil beam algorithm (Helax-TMS's 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 patient's 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.

Comparison of dose calculation algorithms for treatment planning in external photon beam therapy for clinical situations

A study of the performance of five commercial radiotherapy treatment planning systems (TPSs) for common treatment sites regarding their ability to model heterogeneities and scattered photons has been performed. The comparison was based on CT information for prostate, head and neck, breast and lung cancer cases. The TPSs were installed locally at different institutions and commissioned for clinical use based on local procedures. For the evaluation, beam qualities as identical as possible were used: low energy (6 MV) and high energy (15 or 18 MV) x-rays. All relevant anatomical structures were outlined and simple treatment plans were set up. Images, structures and plans were exported, anonymized and distributed to the participating institutions using the DICOM protocol. The plans were then re-calculated locally and exported back for evaluation. The TPSs cover dose calculation techniques from correction-based equivalent path length algorithms to model-based algorithms. These were divided into two groups based on how changes in electron transport are accounted for ((a) not considered and (b) considered). Increasing the complexity from the relatively homogeneous pelvic region to the very inhomogeneous lung region resulted in less accurate dose distributions. Improvements in the calculated dose have been shown when models consider volume scatter and changes in electron transport, especially when the extension of the irradiated volume was limited and when low densities were present in or adjacent to the fields. A Monte Carlo calculated algorithm input data set and a benchmark set for a virtual linear accelerator have been produced which have facilitated the analysis and interpretation of the results. The more sophisticated models in the type b group exhibit changes in both absorbed dose and its distribution which are congruent with the simulations performed by Monte Carlo-based virtual accelerator.