Abstract
Patient’s treatment plan verification covers substantial amount of the quality assurance (QA) resources; this is especially true for Intensity-Modulated Proton Therapy (IMPT). The use of Monte Carlo (MC) simulations in supporting QA has been widely discussed, and several methods have been proposed. In this paper, we studied an alternative approach from the one being currently applied clinically at Centro Nazionale di Adroterapia Oncologica (CNAO). We reanalyzed the previously published data (Molinelli et al. (1)), where 9 patient plans were investigated in which the warning QA threshold of 3% mean dose deviation was crossed. The possibility that these differences between measurement and calculated dose were related to dose modeling (Treatment Planning Systems (TPS) vs. MC), limitations on dose delivery system, or detectors mispositioning was originally explored, but other factors, such as the geometric description of the detectors, were not ruled out. For the purpose of this work, we compared ionization chambers’ measurements with different MC simulation results. It was also studied that some physical effects were introduced by this new approach, for example, inter-detector interference and the delta ray thresholds. The simulations accounting for a detailed geometry typically are superior (statistical difference – p-value around 0.01) to most of the MC simulations used at CNAO (only inferior to the shift approach used). No real improvement was observed in reducing the current delta ray threshold used (100 keV), and no significant interference between ion chambers in the phantom were detected (p-value 0.81). In conclusion, it was observed that the detailed geometrical description improves the agreement between measurement and MC calculations in some cases. But in other cases, position uncertainty represents the dominant uncertainty. The inter-chamber disturbance was not detected for the therapeutic protons energies, and the results from the current delta threshold are acceptable for MC simulations in IMPT.
1 Introduction
Delivering an appropriate radiation therapy dose starts by preparing the most suitable treatment plan for each patient. This is done by conforming the delivered dose to the clinical target volume and avoiding critical organs (2) in order to limit the observed side effects on the surrounding tissue in the patient. Proton beams, with their defined range, can play an important part in increasing this conformity (3). Monte Carlo (MC) simulations are one of the proposed three different dose calculation algorithms, alongside ray trace and pencil beam. Although MC has been considered the gold standard between these approaches in respect to its accuracy, pencil beam model is mostly used in the treatment plan system (TPS) due to its compromise between accuracy and computational time.
After finding the best solution for how to deliver the dose, a verification process is needed in order to check if the equipment is able to deliver the planned treatment fields. Several methods have been proposed, such as the ones by PSI (4) and MD Anderson (5), but at Italian National Center for Oncological Hadron Therapy (CNAO), the method developed by GSI and used at HIT (1, 6), is adopted. CNAO is a hospital-based hadrontherapy facility equipped with a custom synchrotron and Dose Delivery System (DDS) to provide actively scanned proton beams with energies of 62–227 MeV/u and carbon with 115–400 MeV/u, corresponding to ranges in water of 3–32 and 3–27 cm for protons and carbon ions, respectively (7).
Individual treatment plan verifications in the experimental environment can be very time and manpower intensive, and it is prone to dose delivery uncertainties and setup errors. Molinelli et al. (1) presented CNAO’s quality assurance results for all the patients treatment plans verification that have been performed in CNAO with proton beams concerning 1 year (September 2011–August 2012). Nine cases have been found where the quality assurance warning threshold was exceeded, which is fixed at 3% mean absolute deviation between measurements and TPS. Originally, the possibility explored was that these differences between measurement and calculated dose were related to dose modeling (TPS vs. MC), limitations on DDS, or detectors mispositioned (shift), but other factors were not ruled out, such as oversimplification of the dose modeling.
FLUKA (8, 9) was the MC code chosen for this work due to its demonstrated capabilities (10, 11) and available powerful graphical interface (12).
In this work, we have evaluated if improvements could be applied to the MC simulations in order to get better agreement with the measured data on these previously described cases (Figure 1). More specifically, we studied the use of more detailed representation of the detectors (13, 14) and the effect of physical processes that these could introduce in the MC simulation results, for example, the required threshold settings for specific scenarios (15).
Figure 1
2 Materials and Methods
2.1 TPS Patient Plan Verification
The TPS used in CNAO is the CE-marked syngo® RT Planning by Siemens AG Healthcare (Erlangen, Germany) version VB10, which is based on TRiP98 (16, 17).
The current CNAO quality assurance procedure (1, 18) specifies that for each patient, plan verification will be performed (6). For this, a water tank with a 3D detector block controlled by a motorized arm (PTW) is used. This enables to measure the deposited dose at different known depths and positions. This detector block provides a support holder for the ionization chambers (IC), in such way that an individual IC do not mask the direct path of the beam to other the IC (PTW pin point IC – Figure 2). IC measurement values are then compared with the ones calculated by the TPS (equation (1)). For data set analysis, the mean deviation is calculated as the difference between measured (dmeasi) and calculated dose (dcalci), normalized to the maximum beam dose (dmax) and averaged over N IC positions i:
Figure 2
The number of points, N, included in the calculation can be equal or lower than 12, depending on the data set. The TPS provides a 3D-averaged dose gradient for each IC position. Points with a calculated gradient higher than 0.04 Gy mm−1 are excluded from the analysis, since they could not be measured sufficiently accurately due to the finite size of the detector sensitive volume and experimental setup uncertainties. For QA measurements in reference conditions, the applied acceptance threshold is 5% for both mean deviation and SD over a data set.
2.2 Monte Carlo Simulations
FLUKA is a multipurpose MC transport code originally designed for high-energy physics but with extensive use in medical applications (10, 11). For the purpose of this paper, the HADROTHErapy suite of physical settings (known as Defaults) was selected. All geometry updates and modifications were completed with the aid of FLAIR (a graphics user interface of FLUKA).
2.2.1 Current MC-Based Plan Verification
A complete detailed description of CNAO facility, including accelerator design and rooms layout, can be found in the literature (1, 7). For the simulation purposes, the geometry description accounted for the different structures, mainly from the monitors of the DDS, present in the beam path. The validation of this DDS description with FLUKA has been described previously (1). In Figure 3, the photo of the end of the nozzle in one of the treatment rooms is shown together with its description in a 3-dimensional model and its description within the FLUKA simulation.
Figure 3
Current MC patient plan verifications, as per TPS, use a simpler approach to geometrically represent the IC when calculating the dose deposition. All detectors’ structures and holders are not included, and the detector dose is sampled from the dose distribution in a water tank. By doing that, the structure and materials of the IC are not taken into consideration for the simulation. MC obtains the deposited dose in the chambers by calculating the average dose to water over several voxels, corresponding to the active volume of the detector, situated in the positions where the chambers is located.
2.2.2 New Detailed Geometry
The previously described approach, with its geometric approximations and simplifications, obtained deviations below 3% for the majority of studied cases. But for these nine cases where the agreement between the TPS calculations and measurements was above this threshold, we decided to investigate the impact of using a detailed geometry in order to account for the dose disturbances, mainly from scattered particles produced in the wall of the IC and detector holder. In the new detailed geometry, all geometry described above is kept with the inclusion of the PTW3D block and IC description (respecting all structures, dimensions, and material compositions). Detailed technical drawings were obtained from the manufactures (PTW Freiburg). Flair geometry editor was instrumental and extremely helpful in dealing with drawing and 3D visualization (Figure 4). As for the original MC approach, the absolute mean dose deviation is calculated applying equation (1).
Figure 4
2.2.3 Delta Rays
Delta rays are defined as electrons that acquire sufficiently high kinetic energies through collisions so as to enable them to carry this energy a significant distance away from the track of the primary particle and produce their own ionization of absorber atoms (19). The FLUKA “HADROTHErapy” option uses per default delta ray production and transport cuts of 100 keV. We have chosen to vary the threshold limits in order to evaluate if the observed variation between measurements and FLUKA simulations was influenced by the delta rays threshold value.
The dose to water was calculated by averaging the dose deposited in the sensitive volume of the IC. And in order to study the effect of the delta rays threshold, all regions surrounding the sensitive volume had their threshold changed. In this work, we studied the effect of using 10 and 1000 keV in comparison to the default of 100 keV.
2.2.4 Organization of This Work
In total, nine fields from different patients’ plans were analyzed, including MC simulations (for both described geometries) and dose deposition matrices from TPS and IC results (from plan verification quality assurance). The analysis of the data and this section are divided as follows:
Geometry effect – Section 3.1: compared data obtained from the MC with and without the complete geometry, with experimental data and TPS calculated dose values.
Influence of the δ-rays thresholds – Section 3.2: compared data obtained from MC simulations of two different fields for different thresholds (10, 100, and 1000 keV). For that, all regions in, or in direct contact with, the active volume had its threshold modified for the purpose of these simulations (Figure 5).
Chamber–chamber effect – Section 3.3: at CNAO, the measurements are made with half of the proposed number of the detectors (24); so data from MC simulations were compared with two different setups, one with all 24 IC versus the same setup with 12 IC.
Figure 5
3 Results
3.1 The Influence of a Detailed Geometry Implementation
In Figure 6, it was compared for each data set results obtained implementing detailed geometry, Section 2.2.4, in relation to the ones obtained by Molinelli et al. (1).
Figure 6
An advantage was noticed when using a more detailed geometry (MC-NGeo-DDS), as it can be seen by the 6 cases where better results were obtained in comparison to the previous MC-DDS (MC simulations based on DDS log files). In order to understand if the difference between these MC simulations and if the measurements are significant, the relative difference between MC and measurements was calculated and analyzed.
A p-value of 0.003 was obtained between MC-DDS (current MC geometry description) and MC-NGeo-DDS (more detailed geometry description) by using a 2-tailed t-test, which describes that the obtained results by the more detailed geometry approach are significantly better in relation to the current MC Geometry description.
3.2 The Influence of Delta Rays Threshold
The effect in the absorbed dose and computational time was analyzed for two patients’ data sets with different δ-ray thresholds. As described previously, in Section 2.2, with FLUKA MC code, the user is able to set different thresholds for both production and transport of different particles. Initially as expected, some differences were noted between the individual measurements for each threshold and comparison to the current default threshold, set at 100 keV (Figure 7). When compared to the measurements individually and as data set (Figure 7), no comparable advantage was noticed by using different thresholds.
Figure 7
When comparing the computational time when the δ-ray threshold is changed, it was observed that by increasing the threshold (from 100 to 1000 keV), the average time to simulate all primaries reduced by 11.45 ± 3.39%, and when the threshold was reduced (from 100 to 10 keV), the average time to simulate all primaries increased by 43.49 ± 22.02%.
3.3 Chamber–Chamber Effect
Another aspect analyzed was the fact that instead of using the full 24 positions available in the ionization chambers holder (see Figure 4), only 12 positions were used, allowing for investigating the influence of chamber–chamber effect. A 2-tailed t-test was performed, as in Section 3.1, and no statistical significant difference (p-value 0.996) was found between both simulations with 12 or 24 chambers. Figure 8 shows the calculated deviations for the different data sets for both MC and TPS.
Figure 8
4 Discussion
4.1 The Influence of a Detailed Geometry
In this work, we evaluated the effect of the IC geometry description in MC simulation for patient plan verification. We compared our geometrical description with the current approach used. Figure 6 showed that by improving the details of the detectors geometry description, on average, we obtained a mean deviation of 1.90% with 0.63% 1 SD for the 9 cases in comparison to the current method (1), which obtained a mean deviation of 2.36% (0.75%).
Another source of uncertainty is the positioning of the phantom. In order to evaluate if the deviations found were introduced not only by the MC simulations but also by the position of these detectors during measurements, we simulate the effect of introducing this uncertainty.
We simulated different detector positions within ±1 mm in all direction for one of the data set (Data Set 2) around the original position, in total, 27 positions miming uncertainties in the detector positioning. A new minimum of 1.94% was found at (1, −1, 1) as dx, dy, and dz, respectively, in comparison to 2.40% as reported by Molinelli et al. (1).
The obtained mean deviation in respect to the applied offset can be seen in Figure 9. It can be seen that the obtained deviation varies with the positioning of the detector block in a systematic manner, where a minimum and a maximum deviation can be obtained by optimizing the detector position. This also shows that for this data set, the importance of a proper positioning of the water phantom.
Figure 9
4.2 Chamber–Chamber Effect
In addition to the benefits of using a more detailed geometry description, additional points needed to be evaluated as possible contribution to errors and uncertainties. The first one was the interference seen by a detector from the interaction of beam to previously positioned detectors. Although this effect had been evaluated for carbon ions (20) in the case of protons which are more susceptible to the broadening, it had not been evaluated. In our study, no significant difference was seen between the measurements with and without the extra detectors.
4.3 The Influence of Delta Rays Threshold
Another possible factor which will influence the MC simulation results with more detailed detector geometry is the choice of delta rays threshold. For this reason, we evaluated threshold in our detailed geometry. We found that current thresholds used by the default, which have been previously analyzed (11), are still sufficient for this detailed geometry, and no improvement was observed by the reduction of these.
5 Conclusion
The use of MC simulations in aiding patient plan verification has been evaluated. In this work, we studied the effect of improving the detectors geometry description in the MC simulations. We showed that even in the most challenging scenarios of very non-uniform fields, a more detailed geometric description of the detectors results in better agreement with the measurements, although at the cost of more computational time (18.8% in average). If taken into considerations that only 9 patient in an entire year period crossed the threshold, this increase of time should not limit the use of a more detailed geometry description. Additionally, we saw that for few cases where the uncertainty of mispositioning was more relevant than the modeling uncertainties, the use of detailed geometry description in the MC simulation was not able to improve agreement with measurements. For these cases, it was only possible to obtain a better agreement after the detector position was shifted.
Statements
Author contributions
All the authors participated in the different stages of the work, from conception, analysis, and writing.
Funding
This research project has been supported by a Marie Curie Early Initial Training Network Fellowship of the European Community’s Seventh Programme under contract number PITN-GA-2010-264552-ENTERVISION. In addition, TVML was also supported by UCL and the Wellcome Trust ISSF Award (097815/Z/11/B).
Acknowledgments
We would like to thank the help and input from Mr. A. Fedynich, Mr. L. Morejon, and Dr. P. Velten. In addition, this work would not be made possible without the valuable help from Dr. Vlachoudis at CERN and the Medical Physics group at CNAO.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
MolinelliSMairaniAMirandolaAVilches FreixasGTessonnierTGiordanengoSet alDosimetric accuracy assessment of a treatment plan verification system for scanned proton beam radiotherapy: one-year experimental results and Monte Carlo analysis of the involved uncertainties. Phys Med Biol (2013) 58(11):3837–47.10.1088/0031-9155/58/11/3837
2
ICRU. Report 78: prescribing, recording, and reporting proton-beam therapy. J ICRU (2007) 7(2):135–139.
3
KraftG. Tumor therapy with heavy charged particles. Prog Part Nucl Phys (2000) 45:S473–544.10.1016/S0146-6410(00)00112-5
4
LomaxABöhringerTBolsiACorayDEmertFGoiteinGet alTreatment planning and verification of proton therapy using spot scanning: initial experiences. Med Phys (2004) 31(11):3150.10.1118/1.1779371
5
ArjomandyBSahooNCiangaruGZhuRSongXGillinM. Verification of patient-specific dose distributions in proton therapy using a commercial two-dimensional ion chamber array. Med Phys (2010) 37(11):5831.10.1118/1.3505011
6
KargerCPJäkelOHartmannGH. A system for three-dimensional dosimetric verification of treatment plans in intensity-modulated radiotherapy with heavy ions. Med Phys (1999) 26(10):2125–32.10.1118/1.598728
7
RossiS. The status of CNAO. Eur Phys J Plus (2011) 126(8):78.10.1140/epjp/i2011-11078-8
8
FerrariASalaPRFassòARanftJ. FLUKA: A Multi-Particle Transport Code (Program Version 2005). Geneva: CERN (2005).
9
BöhlenTTCeruttiFChinMPWFassòAFerrariAOrtegaPGet alThe FLUKA code: developments and challenges for high energy and medical applications. Nucl Data Sheets (2014) 120(C):211–4.10.1016/j.nds.2014.07.049
10
MairaniABöhlenTTSchiaviATessonnierTMolinelliSBronsSet alA Monte Carlo-based treatment planning tool for proton therapy. Phys Med Biol (2013) 58(8):2471–90.10.1088/0031-9155/58/8/2471
11
ParodiKMairaniABronsSHaschBGSommererFNaumannJet alMonte Carlo simulations to support start-up and treatment planning of scanned proton and carbon ion therapy at a synchrotron-based facility. Phys Med Biol (2012) 57(12):3759–84.10.1088/0031-9155/57/12/3759
12
VlachoudisV. FLAIR: a powerful but user friendly graphical interface for FLUKA. International Conference on Mathematics, Computational Methods & Reactor Physics, Saratoga Springs, New York, May 3–7, 2009. LaGrange Park, IL: American Nuclear Society, (2009).
13
GoiteinM. A technique for calculating the influence of thin inhomogeneities on charged particle beams. Med Phys (1978) 5(4):258–64.10.1118/1.594509
14
UrieMMSistersonJMKoehlerAMGoiteinMZoesmanJ. Proton beam penumbra: effects of separation between patient and beam modifying devices. Med Phys (1986) 13(5):734–41.10.1118/1.595974
15
BöhlenTTDosanjhMFerrariAGudowskaIMairaniA. FLUKA simulations of the response of tissue-equivalent proportional counters to ion beams for applications in hadron therapy and space. Phys Med Biol (2011) 56(20):6545–61.10.1088/0031-9155/56/20/002
16
KrämerMScholzM. Treatment planning for heavy-ion radiotherapy: calculation and optimization of biologically effective dose. Phys Med Biol (2000) 45(11):3319–30.10.1088/0031-9155/45/11/313
17
KrämerMJäkelOHabererTKraftGSchardtDWeberU. Treatment planning for heavy-ion radiotherapy: physical beam model and dose optimization. Phys Med Biol (2000) 45(11):3299–317.10.1088/0031-9155/45/11/313
18
CioccaMMairaniAMirandolaAMolinelliSFreixasGVDonettiMet alOC-0057 commissioning and quality assurance of scanned proton beams produced by a synchrotron for particle radiotherapy. Radiother Oncol (2012) 103:S22–3.10.1016/S0167-8140(12)70396-1
19
IAEA. Radiation Oncology Physics: A Handbook for Teachers and Students. 882 Technical report. Vienna, Austria: IAEA (2005).
20
ParcerisaDS. Experimental and Computational Investigation of Water-to-Air Stopping Power Ratio for Ion Chamber Dosimetry in Carbon Ion Radiotherapy. PhD thesis. Germany: University of Heidelberg (2012).
Summary
Keywords
Monte Carlo calculations, treatment plan verification, proton therapy, delta ray effect, dose disturbance
Citation
Lima TVM, Dosanjh M, Ferrari A, Molineli S, Ciocca M and Mairani A (2016) Monte Carlo Calculations Supporting Patient Plan Verification in Proton Therapy. Front. Oncol. 6:62. doi: 10.3389/fonc.2016.00062
Received
29 September 2015
Accepted
04 March 2016
Published
18 March 2016
Volume
6 - 2016
Edited by
Francis A. Cucinotta, University of Nevada Las Vegas, USA
Reviewed by
Dalong Pang, Georgetown University Hospital, USA; John Eley, University of Maryland School of Medicine, USA
Updates
Copyright
© 2016 Lima, Dosanjh, Ferrari, Molineli, Ciocca and Mairani.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Thiago V. M. Lima, thiago.vmlima@ksa.ch
Specialty section: This article was submitted to Radiation Oncology, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.