- Open Access
Monte Carlo vs. Pencil Beam based optimization of stereotactic lung IMRT
© Sikora et al; licensee BioMed Central Ltd. 2009
Received: 8 October 2009
Accepted: 12 December 2009
Published: 12 December 2009
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.
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.
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.
Optimization of stereotactic lung intensity-modulated radiotherapy (lung IMSRT) is challenging for two reasons. First, because of uncertainties in dose calculation in the presence of tissue interfaces (between lung and tumour), and the commonly small fields. Second, because of the uncertainties in dose calculation when optimizing on a static patient model representing a moving target.
Although it has been well documented that conventional (Pencil beam (PB), superposition/convolution) algorithms fail to some degree when calculating dose to lung [1–12], most treatment planning systems (TPSs) use these algorithms both for optimization of lung IMSRT and for final dose calculation, and MC only as a benchmarking tool. Fraass et al  have stated that the use of MC dose calculation algorithms for clinical planning improves the dose accuracy in heterogeneous regions of lung and bony anatomy, in particular when applying very small field sizes which exhibit lateral electron disequilibrium effects. Ideally, a TPS comprises a MC dose engine, since an inaccurate dose algorithm will not only introduce dose errors, but will also lead to wrongly optimized treatment plans.
The main argument for conventional algorithms has been that the MC dose calculation algorithms are too slow, especially for use in optimization of IMSRT . However, with the computational power offered by modern computers combined with an efficient MC system [15, 16], it should be feasible to use MC not only for recalculation of IMSRT treatment plans but also for optimization itself. The biological IMRT MC-TPS HYPERION[17, 18] has been used clinically since 2002. The current version allows for optimization based on both fsPB and MC. In case of MC, it uses a Virtual Source Model (VSM) of the accelerator head [15, 16] together with the XVMC dose engine for simulation in the patient .
In this study, we first compare our conventional fsPB and MC algorithms to measurements in a lung phantom in order to benchmark their overall accuracy. Then, we apply both algorithms in optimization of lung IMSRT to investigate if direct optimization with MC can provide an advantage over optimization with a conventional algorithm and recalculation with MC.
The comparison of MC with conventional algorithms is not independent of the patient model. Commonly, lung IMSRT is performed on a static model of the moving target and different approaches are used to represent the distribution of densities within the Planning Target Volume (PTV). In this study, we therefore also investigate the influence of various patient models of the lung patient geometry on the results of optimization and compare these results to the real 4D accumulated dose, recalculated with MC on each of 8 computer tomography scans (CTs) representing 8 breathing phases .
Materials and methods
Verification of the MC and fsPB dose algorithms in a lung phantom
The MC and fsPB algorithms investigated in this study have been presented in detail elsewhere (see [19, 16] and [21, 22], respectively). Both algorithms calculate dose-to-tissue, since the fsPB is commissioned based on the result of simulation with the MC algorithm.
The lung phantom
The lung phantom dose distribution calculation, measurements and 2D dose comparisons
For each tumour size, a conformal plan consisting of one vertical beam was created using 6 MV and 15 MV nominal beam energies, and the dose delivered to the phantom was calculated both with MC and fsPB. The maximum dose in the tumour was below 1 Gy to avoid saturation of the radiographic films.
A fine dose grid of 1 mm3 was used and the statistical uncertainty of the MC simulations was set to 0.5%. The presence of the dosimetric film was taken into account in the dose calculations by introducing a 1 mm thick layer filled with a density of ρ = 1.1 at the position of the film in the CT scan of the lung phantom (as the CT acquisition was made without the film).
Radiographic Kodak X-Omat V films were used for the measurements. These were cut to fit the dimensions of the phantom and placed into the central plane of the cork cube between the hemispheres of the tumour perpendicular to the beam direction (Figure 1). The cork cube was then tightly fixed to the container and the phantom was positioned with the help of the positioning markers. An uncertainty in the position of the film plane of up to 2 degrees relatively to the beam direction was impossible to avoid because of how the slit in the cork cube was cut (Figure 1). An Elekta Synergy S accelerator was used to irradiate the phantoms. All films were developed at the same time and scanned with a Radlink LaserPro 16 scanner. The optical densities of the films were converted to absolute dose by using a calibration curve created for solid water.
The MC and fsPB dose distributions were compared to the 2D film measurements by using the γ-index method , with an acceptance criteria of 3% dose difference and 3 mm distance-to-agreement (3%/3 mm). Also the tighter 2%/2 mm criteria was investigated. The artefacts from the positioning markers of the phantom were excluded from the analysis and no attempt was made to correct for the slight uncertainty in the positioning of the film. The in-plane and cross-plane profiles through the isocentre where extracted from the 2D film measurements and compared to the MC and fsPB dose distributions. Additionally, depth dose curves were compared between MC and fsPB. The depth dose curves were not measured due to unprecise alignment and a relatively high density contribution of the film when its plane is located along the beam direction.
The treatment planning optimization tool
a poisson cell-kill EUD (equivalent uniform dose) model to maximize tumour control ,
and a maximum overdosage constraint to avoid target hotspots and to control maximum organ at risk doses in terms of the root mean square (rms) overdosage above a given threshold dose.
Plan optimization is a two-stage process. In stage I, the fsPB computes dose distributions for a large number of beamlets that constitute the fluence profiles. Then, the weights of these beamlets are optimized to yield the idealized fluence profiles. Stage II starts with an initial segmentation of the fluence profiles into a sequence of deliverable MLC segments. Here, the segment doses can either be obtained from a concatenation of their beamlet doses or by a MC calculation. The weights and shapes of these sequences are optimized as described in . In case a segment changes shape, the MC calculation is re-run. An approach via MC-calculated beamlets meets two obstacles: firstly, for a given incident history density, the point-wise dose uncertainty of a beamlet is greater because the dose is smeared out over a greater volume, which in turn causes instabilities in the optimization. Secondly, by definition a beamlet dose cannot include MLC effects which requires a re-computation of the full segment doses in stage II anyway and leads to wrong (usually too steep) field penumbra in stage I.
Optimization of lung IMSRT plans with MC and PB
Patient data and models
One example patient which had previously been treated with lung IMSRT at the University Hospital of Tübingen was used in the study. The patient had a tumour with a diameter of 2.4 cm which was located posteriorly in the lower right lung. The breathing excursion was 2.9 cm mainly in the cranio-cadual direction. A respiratory-correlated CT (RCCT) dataset was acquired with a Siemens Sensation Open scanner reconstructed with 1 × 1 × 3 mm3 voxel size. The CT dataset was grouped into eight CT sets (0/25/50/75% inhale and 100/75/50/25% exhale), where the 0% inhale CT was used as the planning geometry (the static exhale (planning) CT). Contours of the clinical target volumes (CTVs) from all breathing phases as well as of the organs at risk (OARs) were defined and approved by a radiologist. The planning target volume (PTV) was defined as the internal target volume (ITV), i.e. the union of the CTVs from all breathing phases, expanded by 2 mm in order to account for setup uncertainties. This resulted in a PTV of 47.1 cm3, around three times larger than the clinical target volume (CTV).
one phase static CT : this model uses the PTV and the exhale planning CT without any attempt to correct for density variations,
minimum density overwrite, one phase static CT : in this model, a density less then 0.4 (an empirical value to avoid erratic fluence modulation caused by density-related underdosage in the target volume) was raised to this value within the PTV on the planning CT to mimic the effective density during free breathing irradiation,
average CT: this model uses the PTV and a superposition of all RCCTs, such that each voxel has a density equal to the weighted average of the Hounsfield values from all breathing phases.
Static lung IMSRT planning
Prescribed and resulting isoeffects/EUDs for all patient models plans optimized with fsPB and MC and recalculated with 4DMC; the recalculated MC dose on the static patient models of the fsPB optimized plans are shown in brackets.
model 1, one phase static CT
target* poisson EUD
target* rms overdosage
lung R, mean dose
lung R, mean damage
skin, rms overdosage
model 2, minimum density overwrite one phase static CT
target* poisson EUD
target* rms overdosage
lung R, mean dose
lung R, mean damage
skin, rms overdosage
model 3, average CT
target* poisson EUD
target* rms overdosage
lung R, mean dose
lung R, mean damage
skin, rms overdosage
A beamlet size of 4 × 2 mm2 was used and segments smaller than 0.64 cm2 were not allowed. We used a 2.5 × 2.5 × 2.5 mm3 dose calculation grid size and a 3% statistical uncertainty per MLC segment for the MC calculation. After optimization, the fsPB plans were recalculated in the static geometries with the MC dose engine. Additionally, all plans were recalculated by 4DMC.
For 4D MC plan calculation (4DMC), the dose was computed in each of the eight static geometries of the RCCT dataset. These MC doses per instance were weighted according to their share of the breathing cycle and accumulated in a common reference geometry (here, the exhale planning CT) by dose warping derived from deformable registration. Details of the 4DMC recalculation method are described in details in . It is important to notice that the PTV encloses a volume where the tumor can be found with certain probability which depends on the breathing pattern. Therefore, it is not possible to evaluate the real dose to the CTV without employing the 4DMC recalculation where the accumulated dose in the CTV is calculated. In the following, we use a common term - target volume, which denotes the PTV for static planning and the CTV for 4DMC recalculation, respectively.
Experimental verification of MC and fsPB dose engines
Agreement between the 2D dose distributions as measured with film and calculated with MC and fsPB for the 3%/3 mm and 2%/2 mm γ acceptance criteria; the acceptance (in %) is listed for the whole region (all), the tumour region (tumour) and the region outside the tumour (lung) and for all tumour sizes.
MC 3%/3 mm
fsPB 3%/3 mm
MC 2%/2 mm
fsPB 2%/2 mm
Tumour diameter [cm]
Static IMSRT dose distributions
The MC recalculation on the static patient models of the fsPB optimized plans shows that most of the prescribed OAR constraints were actually violated. In addition, we found that the fsPB dose calculation algorithm underestimated the target-EUD in case of 6 MV and overestimated the target-EUD in case of 15 MV for all patient models.
In this study we have shown how dose calculation of small fields in the presence of tissue heterogeneities and static modelling of a moving target influence fsPB and MC dose calculation and optimization for lung IMSRT. The results of the lung IMSRT planning are presented for one extreme case (small tumor, large movements, density inhomogeneities) in order to emphasize possible dose calculation and patient model uncertainties and less extreme cases will produce smaller problems.
While the MC dose calculation performed well for a static lung phantom, the fsPB algorithm always underestimated the target dose compared to film measurements. This confirms what has been found by others for various superposition/convolution and collapsed cone algorithms [10, 12, 28, 8, 3, 5, 29]. The MC can safely be regarded as the superior algorithm for dose calculation on a static geometry.
When looking at real patient data of lung IMSRT, the picture gets more complicated and it is difficult to interpret the results from fsPB planning. Vanderstraeten et al.  investigated dose calculation with various commercially available conventional dose engines, and showed overestimation or underestimation of the target dose when comparing to the MC calculated dose, depending on the algorithm. Our results showed that our fsPB algorithm both underestimated and overestimated the target-EUD depending on the energy (6 or 15 MV). This is related to uncertainties in penumbra widening factors (especially parameter fu 1(ρ)) of the fit functions implemented in the fsPB for low densities (Figure 4, Jeleń et al ). The dose to the ipsilateral lung calculated by the fsPB was underestimated for all static plans while MC static plans agreed very well with the 4DMC recalculations. Thus, for this patient example lung IMSRT treatment planning with the fsPB results in higher complication probability then MC-based planning.
In addition to an unpredictable target-EUD, optimization with fsPB led to violation of most OAR constraints for both low and high energy photon beams when compared to MC recalculated dose on both the static and 4D patient models. This means that optimization with fsPB and recalculation with MC is impractical, because the violated OAR constraints would have to be readjusted and the whole optimization procedure rerun. Bearing this in mind, optimization with fsPB and recalculation with MC may not necessarily be more efficient than MC based optimization which needed less than 1.5 hours in our case. Although MC is superior to fsPB on a static geometry, the real accumulated dose to the patient is influenced by the breathing motion during delivery. In our study, the 4D MC recalculations of the MC optimized treatment plans showed that the target doses calculated on the static patient models were always lower than the 4D accumulated target dose.
The differences between static and 4D recalculated MC dose depended on the patient model. The largest difference, as expected, was found for the one phase static CT model, where a large part of the target was occupied by low density lung tissue to where the optimizer tried to deliver a higher fluence in order to achieve a homogeneous target dose. The minimum density overwrite one phase static CT model worked well for our patient, especially for the 6 MV plan. The advantage of this patient model is that the optimizer is not miss-leaded by low densities in the target. In contrast, the average CT model still contains a significant fraction of low density lung tissue (like the one phase static CT model) to where the optimizer also tries to boost the dose but which is rarely visited by the tumor.
For the fsPB optimized plans, 4D MC recalculation displayed a complex picture where dose calculation uncertainties and geometrical uncertainties either added up or cancelled out. Because only geometrical uncertainties were present with the MC optimized plans, these were more predictable. Although always underestimating the target dose compared to 4D MC, the MC optimized plans showed only small deviations in the OAR isoeffects/EUDs, especially for relatively static organs. The 4DMC dose were higher because the dose "follows" the higher density of the tumor during breathing. We therefore consider MC based optimization on static patient models to be a safe method for lung IMSRT planning for the minimum density or average CT models tested in this study.
4D MC optimization was not considered in the present paper. It eliminates the uncertainties in the geometrical model and provides the optimal dose for free breathing IMSRT treatment with a negligible increase of optimization and calculation time at the price of a greater dependency on the consistency between predicted and realized patient breathing motion .
The MC dose engine was superior to fsPB in presence of lateral electron disequilibrium on static geometries. With an efficient MC system, MC based optimization of lung IMSRT is feasible, also given clinical time constraints. MC dose optimization on static patient models always underestimated the CTV dose compared to 4DMC recalculations, but OAR differences were very small. Therefore, MC optimization of lung IMSRT using a static patient model is recommended.
We are grateful to Prof. Dirk Verellen from the University Hospital Brussels for kindly lending the lung phantom. The author is thankful to Liv Bolstad Hysing for helpful discussions and remarks concerning the manuscript. This work was supported by the European Community under a Marie Curie Host Fellowship for Early Stage Researchers Training and Elekta AB, Sweden.
- Solberg TD, DeMarco JJ, Holly FE, Smathers JB, DeSalles AA: Monte Carlo treatment planning for stereotactic radiosurgery. Radiother Oncol 1998, 49: 73-84. 10.1016/S0167-8140(98)00065-6View ArticlePubMedGoogle Scholar
- Ma CM, Mok E, Kapur A, Pawlicki T, Findley D, Brain S, Forster K, Boyer AL: Clinical implementation of a Monte Carlo treatment planning system. Med Phys 1999,26(10):2133-43. 10.1118/1.598729View ArticlePubMedGoogle Scholar
- Paelinck L, Reynaert N, Thierens H, De Neve W, De Wagter C: Experimental verification of lung dose with radiochromic film: comparison with Monte Carlo simulations and commercially available treatment planning systems. Phys Med Biol 2005,50(9):2055-69. 10.1088/0031-9155/50/9/009View ArticlePubMedGoogle Scholar
- Arnfield MR, Siantar CH, Siebers J, Garmon P, Cox L, Mohan R: The impact of electron transport on the accuracy of computed dose. Med Phys 2000,27(6):1266-74. 10.1118/1.599004View ArticlePubMedGoogle Scholar
- Solberg TD, Holly FE, De Salles AA, Wallace RE, Smathers JB: Implications of tissue heterogeneity for radiosurgery in head and neck tumors. Int J Radiat Oncol Biol Phys 1995, 32: 235-9.View ArticlePubMedGoogle Scholar
- Aspradakis MM, Morrison RH, Richmond ND, Steele A: Experimental verification of convolution/superposition photon dose calculations for radiotherapy treatment planning. Phys Med Biol 2003,48(17):2873-93. 10.1088/0031-9155/48/17/309View ArticlePubMedGoogle Scholar
- McDermott PN, He T, DeYoung A: Dose calculation accuracy of lung planning with a commercial IMRT treatment planning system. J Appl Clin Med Phys 2003,4(4):341-51. 10.1120/1.1623172View ArticlePubMedGoogle Scholar
- Tsiakalos MF, Theodorou K, Kappas C, Zefkili S, Rosenwold JC: Analysis of the penumbra enlargement in lung versus the quality index of photon beams: a methodology to check the dose calculation algorithm. Med Phys 2004,31(4):943-9. 10.1118/1.1669085View ArticlePubMedGoogle Scholar
- Ding W, Johnston PN, Wong TPY, Bubb IF: Investigation of photon beam models in heterogeneous media of modern radiotherapy. Australas Phys Eng Sci Med 2004,27(2):39-48. 10.1007/BF03178375View ArticlePubMedGoogle Scholar
- Dobler B, Walter C, Knopf A, Fabri D, Loeschel R, Polednik M, Schneider F, Wenz F, Lohr F: Optimization of extracranial stereotactic radiation therapy of small lung lesions using accurate dose calculation algorithms. Radiat Oncol 2006, 1: 45. 10.1186/1748-717X-1-45PubMed CentralView ArticlePubMedGoogle Scholar
- Carrasco P, Jornet N, Duch MA, Weber L, Ginjaume M, Eudaldo T, Jurado D, Ruiz A, Ribas M: Comparison of dose calculation algorithms in phantoms with lung equivalent heterogeneities under conditions of lateral electronic disequilibrium. Med Phys 2004,31(10):2899-911. 10.1118/1.1788932View ArticlePubMedGoogle Scholar
- Krieger T, Sauer OA: Monte Carlo-versus pencil-beam-/collapsed-cone-dose calculation in a heterogeneous multi-layer phantom. Physics in Medicine and Biology 2005,50(5):859-868. 10.1088/0031-9155/50/5/010View ArticlePubMedGoogle Scholar
- Fraass BA, Smathers J, Deye J: Summary and recommendations of a National Cancer Institute workshop on issues limiting the clinical use of Monte Carlo dose calculation algorithms for megavoltage external beam radiation therapy. Med Phys 2003,30(12):3206-16. 10.1118/1.1626990View ArticlePubMedGoogle Scholar
- Van Esch A, Tillikainen L, Pyykkonen J, Tenhunen M, Helminen H, Siljamäki S, Alakuijala J, Paiusco M, Lori M, Huyskens DP: Testing of the analytical anisotropic algorithm for photon dose calculation. Med Phys 2006,33(11):4130-48. 10.1118/1.2358333View ArticlePubMedGoogle Scholar
- Fippel M, Haryanto F, Dohm O, Nusslin F, Kriesen S: A virtual photon energy fluence model for Monte Carlo dose calculation. Medical Physics 2003,30(3):301-311. 10.1118/1.1543152View ArticlePubMedGoogle Scholar
- Sikora M, Dohm O, Alber M: A virtual photon source model of an Elekta linear accelerator with integrated mini MLC for Monte Carlo based IMRT dose calculation. Phys Med Biol 2007,52(15):4449-63. 10.1088/0031-9155/52/15/006View ArticlePubMedGoogle Scholar
- Laub W, Alber M, Birkner M, Nüsslin F: Monte Carlo dose computation for IMRT optimization. Phys Med Biol 2000,45(7):1741-54. 10.1088/0031-9155/45/7/303View ArticlePubMedGoogle Scholar
- Alber M, Birkner M, Laub W, Nüsslin F: Hyperion: An integrated IMRT planning tool. In 13th International Conference on the Use of Computers in Radiation Therapy, Heidelberg Germany. Edited by: schlegel w, bortfeld t. Springer-Verlag Berlin Heidelberg New York; 2000:46-48.View ArticleGoogle Scholar
- Fippel M: Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm. Med Phys 1999,26(8):1466-75. 10.1118/1.598676View ArticlePubMedGoogle Scholar
- Söhn M, Weinmann M, Alber M: Intensity-modulated radiotherapy optimization in a quasi-periodically deforming patient model. Int J Radiat Oncol Biol Phys 2009,75(3):906-14.View ArticlePubMedGoogle Scholar
- Jeleń U, Söhn M, Alber M: A finite size pencil beam for IMRT dose optimization. Phys Med Biol 2005,50(8):1747-66. 10.1088/0031-9155/50/8/009View ArticlePubMedGoogle Scholar
- Jeleń U, Alber M: A finite size pencil beam algorithm for IMRT dose optimization: density corrections. Phys Med Biol 2007,52(3):617-33. 10.1088/0031-9155/52/3/006View ArticlePubMedGoogle Scholar
- Bakai A, Alber M, Nüsslin F: A revision of the gamma-evaluation concept for the comparison of dose distributions. Phys MedBiol 2003,48(21):3543-53.Google Scholar
- Alber M, Nüsslin F: An objective function for radiation treatment optimization based on local biological measures. Phys Med Biol 1999,44(2):479-93. 10.1088/0031-9155/44/2/014View ArticlePubMedGoogle Scholar
- Alber M: Normal tissue dose-effect models in biological dose optimisation. Z Med Phys 2008,18(2):102-10.View ArticlePubMedGoogle Scholar
- Niemierko A: Reporting and analyzing dose distributions: a concept of equivalent uniform dose. Med Phys 1997, 24: 103-10. 10.1118/1.598063View ArticlePubMedGoogle Scholar
- Alber M, Nüsslin F: Optimization of intensity modulated radiotherapy under constraints for static and dynamic MLC delivery. Phys Med Biol 2001,46(12):3229-39. 10.1088/0031-9155/46/12/311View ArticlePubMedGoogle Scholar
- Martens C, Reynaert N, De Wagter C, Nilsson P, Coghe M, Palmans H, Thierens H, De Neve W: 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. Med Phys 2002,29(7):1528-35. 10.1118/1.1487421View ArticlePubMedGoogle Scholar
- Vanderstraeten B, Reynaert N, Paelinck L, Madani I, De Wagter C, De Gersem W, De Neve W, Thierens H: Accuracy of patient dose calculation for lung IMRT: A comparison of Monte Carlo, convolution/superposition, and pencil beam computations. Med Phys 2006,33(9):3149-58. 10.1118/1.2241992View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.