1
|
A plan verification platform for online adaptive proton therapy using deep learning-based Monte–Carlo denoising. Phys Med 2022; 103:18-25. [DOI: 10.1016/j.ejmp.2022.09.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/12/2022] [Revised: 09/27/2022] [Accepted: 09/28/2022] [Indexed: 11/21/2022] Open
|
2
|
Fuchs H, Elia A, Resch AF, Kuess P, Lühr A, Vidal M, Grevillot L, Georg D. Computer-assisted beam modeling for particle therapy. Med Phys 2020; 48:841-851. [PMID: 33283910 PMCID: PMC7986420 DOI: 10.1002/mp.14647] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Revised: 05/17/2020] [Accepted: 06/22/2020] [Indexed: 11/15/2022] Open
Abstract
Purpose To develop a computer‐driven and thus less user‐dependent method, allowing for a simple and straightforward generation of a Monte Carlo (MC) beam model of a scanned proton and carbon ion beam delivery system. Methods In a first step, experimental measurements were performed for proton and carbon ion energies in the available energy ranges. Data included depth dose profiles measured in water and spot sizes in air at various isocenter distances. Using an automated regularization‐based optimization process (AUTO‐BEAM), GATE/Geant4 beam models of the respective beam lines were generated. These were obtained sequentially by using least square weighting functions with and without regularization, to iteratively tune the beam parameters energy, energy spread, beam sigma, divergence, and emittance until a user‐defined agreement was reached. Based on the parameter tuning for a set of energies, a beam model was semi‐automatically generated. The resulting beam models were validated for all centers comparing to independent measurements of laterally integrated depth dose curves and spot sizes in air. For one representative center, three‐dimensional dose cubes were measured and compared to simulations. The method was applied on one research as well as four different clinical beam lines for proton and carbon ions of three different particle therapy centers using synchrotron or cyclotron accelerator systems: (a) MedAustron ion therapy center, (b) University Proton Therapy Dresden, and (c) Center Antoine Lacassagne Nice. Results Particle beam ranges in the MC beam models agreed on average within 0.2 mm compared to measurements for all energies and beam lines. Spot sizes in air (full‐width at half maximum) at all positions differed by less than 0.4% from the measurements. Dose calculation with the beam model for the clinical beam line at MedAustron agreed better than 1.7% in absolute dose for a representative clinical case treated with protons. For protons, beam model generation, including geometry creation, data conversion, and validation, was possible within three working days. The number of iterations required for the optimization process to converge, was found to be similar for all beam line geometries and particle types. Conclusion The presented method was demonstrated to work independently of the beam optics behavior of the different beam lines, particle types, and geometries. Furthermore, it is suitable for non‐expert users and requires only limited user interaction. Beam model validation for different beam lines based on different beam delivery systems, showed good agreement.
Collapse
Affiliation(s)
- Hermann Fuchs
- Division of Medical Radiation Physics, Department of Radiation Oncology, Medical University of Vienna, Währinger Gürtel 18-20, Vienna, 1090, Austria
| | - Alessio Elia
- MedAustron Ion Therapy Center, Wiener Neustadt, 2700, Austria
| | - Andreas F Resch
- Division of Medical Radiation Physics, Department of Radiation Oncology, Medical University of Vienna, Währinger Gürtel 18-20, Vienna, 1090, Austria
| | - Peter Kuess
- Division of Medical Radiation Physics, Department of Radiation Oncology, Medical University of Vienna, Währinger Gürtel 18-20, Vienna, 1090, Austria
| | - Armin Lühr
- OncoRay - National Center for Radiation Research in Oncology, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, Helmholtz-Zentrum Dresden-Rossendorf, Dresden, 01309, Germany.,Helmholtz-Zentrum Dresden-Rossendorf, Institute of Radiooncology - OncoRay, Dresden, 01309, Germany.,Department of Medical Physics, Faculty of Physics, TU Dortmund University, Dortmund, 44227, Germany
| | - Marie Vidal
- Center Antoine Lacassagne, Nice, 06189, France
| | - Loïc Grevillot
- MedAustron Ion Therapy Center, Wiener Neustadt, 2700, Austria
| | - Dietmar Georg
- Division of Medical Radiation Physics, Department of Radiation Oncology, Medical University of Vienna, Währinger Gürtel 18-20, Vienna, 1090, Austria
| |
Collapse
|
3
|
Guardiola C, De Marzi L, Prezado Y. Verification of a Monte Carlo dose calculation engine in proton minibeam radiotherapy in a passive scattering beamline for preclinical trials. Br J Radiol 2019; 93:20190578. [PMID: 31868523 DOI: 10.1259/bjr.20190578] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
OBJECTIVES Proton minibeam radiation therapy (pMBRT) is a novel therapeutic strategy that combines the benefits of proton therapy with the remarkable normal tissue preservation observed with the use of submillimetric spatially fractionated beams. This promising technique has been implemented at the Institut Curie-Proton therapy centre (ICPO) using a first prototype of a multislit collimator. The purpose of this work was to develop a Monte Carlo-based dose calculation engine to reliably guide preclinical studies at ICPO. METHODS The whole "Y1"-passive beamline at the ICPO, including pMBRT implementation, was modelled using the Monte Carlo GATE v. 7.0 code. A clinically relevant proton energy (100 MeV) was used as starting point. Minibeam generation by means of the brass collimator used in the first experiments was modelled. A virtual source was modelled at the exit of the beamline nozzle and outcomes were compared with dosimetric measurements performed with EBT3 gafchromic films and a diamond detector in water. Dose distributions were recorded in a water phantom and in rat CT images (7-week-old male Fischer rats). RESULTS The dose calculation engine was benchmarked against experimental data and was then used to assess dose distributions in CT images of a rat, resulting from different irradiation configurations used in several experiments. It reduced computational time by an order of magnitude. This allows us to speed up simulations for in vivo trials, where we obtained peak-to-valley dose ratios of 1.20 ± 0.05 and 6.1 ± 0.2 for proton minibeam irradiations targeting the tumour and crossing the rat head. Tumour eradication was observed in the 67 and 22% of the animals treated respectively. CONCLUSION A Monte Carlo dose calculation engine for pMBRT implementation with mechanical collimation has been developed. This tool can be used to guide and interpret the results of in vivo trials. ADVANCES IN KNOWLEDGE This is the first Monte Carlo dose engine for pMBRT that is being used to guide preclinical trials in a clinical proton therapy centre.
Collapse
Affiliation(s)
- Consuelo Guardiola
- Centre National de la Recherche Scientifique (CNRS); Universités Paris 11 and Paris 7, Laboratoire d'Imagerie et Modélisation en Neurobiologie et Cancérologie (IMNC), Orsay Cedex, 91405, France
| | - Ludovic De Marzi
- Institut Curie, PSL Research University, Centre de protonthérapie d'Orsay, Campus universitaire, bâtiment 101, Orsay 91898, France.,Institut Curie, University Paris Saclay, PSL Research University, Inserm U 1021-CNRS UMR 3347, Orsay, France
| | - Yolanda Prezado
- Centre National de la Recherche Scientifique (CNRS); Universités Paris 11 and Paris 7, Laboratoire d'Imagerie et Modélisation en Neurobiologie et Cancérologie (IMNC), Orsay Cedex, 91405, France
| |
Collapse
|
4
|
Experimental characterisation of a proton kernel model for pencil beam scanning techniques. Phys Med 2019; 64:195-203. [DOI: 10.1016/j.ejmp.2019.07.013] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/05/2018] [Revised: 06/20/2019] [Accepted: 07/17/2019] [Indexed: 11/24/2022] Open
|
5
|
Experimental Set-up for FLASH Proton Irradiation of Small Animals Using a Clinical System. Int J Radiat Oncol Biol Phys 2018; 102:619-626. [DOI: 10.1016/j.ijrobp.2018.06.403] [Citation(s) in RCA: 128] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2017] [Revised: 06/15/2018] [Accepted: 06/30/2018] [Indexed: 11/18/2022]
|
6
|
De Marzi L, Patriarca A, Nauraye C, Hierso E, Dendale R, Guardiola C, Prezado Y. Implementation of planar proton minibeam radiation therapy using a pencil beam scanning system: A proof of concept study. Med Phys 2018; 45:5305-5316. [PMID: 30311639 DOI: 10.1002/mp.13209] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2018] [Revised: 07/26/2018] [Accepted: 09/02/2018] [Indexed: 11/06/2022] Open
Abstract
PURPOSE Proton minibeam radiation therapy (pMBRT) is an innovative approach that combines the advantages of minibeam radiation therapy with the more precise ballistics of protons to further reduce the side effects of radiation. One of the main challenges of this approach is the generation of very narrow proton pencil beams with an adequate dose-rate to treat patients within a reasonable treatment time (several minutes) in existing clinical facilities. The aim of this study was to demonstrate the feasibility of implementing pMBRT by combining the pencil beam scanning (PBS) technique with the use of multislit collimators. This proof of concept study of pMBRT with a clinical system is intended to guide upcoming biological experiments. METHODS Monte Carlo simulations (TOPAS v3.1.p2) were used to design a suitable multislit collimator to implement planar pMBRT for conventional pencil beam scanning settings. Dose distributions (depth-dose curves, lateral profiles, Peak-to-Valley Dose Ratio (PVDR) and dose-rates) for different proton beam energies were assessed by means of Monte Carlo simulations and experimental measurements in a water tank using commercial ionization chambers and a new p-type silicon diode, the IBA RAZOR. An analytical intensity-modulated dose calculation algorithm designed to optimize the weight of individual Bragg peaks composing the field was also developed and validated. RESULTS Proton minibeams were then obtained using a brass multislit collimator with five slits measuring 2 cm × 400 μm in width with a center-to-center distance of 4 mm. The measured and calculated dose distributions (depth-dose curves and lateral profiles) showed a good agreement. Spread-out Bragg peaks (SOBP) and homogeneous dose distributions around the target were obtained by means of intensity modulation of Bragg peaks, while maintaining spatial fractionation at shallow depths. Mean dose-rates of 0.12 and 0.09 Gy/s were obtained for one iso-energy layer and a SOBP conditions in the presence of multislit collimator. CONCLUSIONS This study demonstrates the feasibility of implementing pMBRT on a PBS system. It also confirms the reliability of RAZOR detector for pMBRT dosimetry. This newly developed experimental methodology will support the design of future preclinical research with pMBRT.
Collapse
Affiliation(s)
- Ludovic De Marzi
- Institut Curie, PSL Research University, Centre de protonthérapie d'Orsay, Campus universitaire, bâtiment 101, Orsay, 91898, France
| | - Annalisa Patriarca
- Institut Curie, PSL Research University, Centre de protonthérapie d'Orsay, Campus universitaire, bâtiment 101, Orsay, 91898, France
| | - Catherine Nauraye
- Institut Curie, PSL Research University, Centre de protonthérapie d'Orsay, Campus universitaire, bâtiment 101, Orsay, 91898, France
| | - Eric Hierso
- Institut Curie, PSL Research University, Centre de protonthérapie d'Orsay, Campus universitaire, bâtiment 101, Orsay, 91898, France
| | - Rémi Dendale
- Institut Curie, PSL Research University, Centre de protonthérapie d'Orsay, Campus universitaire, bâtiment 101, Orsay, 91898, France
| | - Consuelo Guardiola
- IMNC-UMR 8165, CNRS, Paris 7 and Paris 11 Universities, 15 rue Georges Clemenceau, Orsay Cedex, 91405, France
| | - Yolanda Prezado
- IMNC-UMR 8165, CNRS, Paris 7 and Paris 11 Universities, 15 rue Georges Clemenceau, Orsay Cedex, 91405, France
| |
Collapse
|
7
|
Kim DH, Cho S, Jo K, Shin E, Hong CS, Han Y, Suh TS, Lim DH, Choi DH. Proton range verification in inhomogeneous tissue: Treatment planning system vs. measurement vs. Monte Carlo simulation. PLoS One 2018; 13:e0193904. [PMID: 29505589 PMCID: PMC5837130 DOI: 10.1371/journal.pone.0193904] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2017] [Accepted: 02/05/2018] [Indexed: 11/18/2022] Open
Abstract
In particle radiotherapy, range uncertainty is an important issue that needs to be overcome. Because high-dose conformality can be achieved using a particle beam, a small uncertainty can affect tumor control or cause normal-tissue complications. From this perspective, the treatment planning system (TPS) must be accurate. However, there is a well-known inaccuracy regarding dose computation in heterogeneous media. This means that verifying the uncertainty level is one of the prerequisites for TPS commissioning. We evaluated the range accuracy of the dose computation algorithm implemented in a commercial TPS, and Monte Carlo (MC) simulation against measurement using a CT calibration phantom. A treatment plan was produced for eight different materials plugged into a phantom, and two-dimensional doses were measured using a chamber array. The measurement setup and beam delivery were simulated by MC code. For an infinite solid water phantom, the gamma passing rate between the measurement and TPS was 97.7%, and that between the measurement and MC was 96.5%. However, gamma passing rates between the measurement and TPS were 49.4% for the lung and 67.8% for bone, and between the measurement and MC were 85.6% for the lung and 100.0% for bone tissue. For adipose, breast, brain, liver, and bone mineral, the gamma passing rates computed by TPS were 91.7%, 90.6%, 81.7%, 85.6%, and 85.6%, respectively. The gamma passing rates for MC for adipose, breast, brain, liver, and bone mineral were 100.0%, 97.2%, 95.0%, 98.9%, and 97.8%, respectively. In conclusion, the described procedure successfully evaluated the allowable range uncertainty for TPS commissioning. The TPS dose calculation is inefficient in heterogeneous media with large differences in density, such as lung or bone tissue. Therefore, the limitations of TPS in heterogeneous media should be understood and applied in clinical practice.
Collapse
Affiliation(s)
- Dae-Hyun Kim
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
- Department of Biomedical Engineering and Research Institute of Biomedical Engineering, College of Medicine, The Catholic University of Korea, Seoul, Republic of Korea
| | - Sungkoo Cho
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
| | - Kwanghyun Jo
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
| | - EunHyuk Shin
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
| | - Chae-Seon Hong
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
| | - Youngyih Han
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
- * E-mail: (YH); (TS)
| | - Tae-Suk Suh
- Department of Biomedical Engineering and Research Institute of Biomedical Engineering, College of Medicine, The Catholic University of Korea, Seoul, Republic of Korea
- * E-mail: (YH); (TS)
| | - Do Hoon Lim
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
| | - Doo Ho Choi
- Department of Radiation Oncology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, Republic of Korea
| |
Collapse
|
8
|
Moignier C, Tromson D, de Marzi L, Marsolat F, García Hernández JC, Agelou M, Pomorski M, Woo R, Bourbotte JM, Moignau F, Lazaro D, Mazal A. Development of a synthetic single crystal diamond dosimeter for dose measurement of clinical proton beams. Phys Med Biol 2017; 62:5417-5439. [PMID: 28604370 DOI: 10.1088/1361-6560/aa70cf] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The scope of this work was to develop a synthetic single crystal diamond dosimeter (SCDD-Pro) for accurate relative dose measurements of clinical proton beams in water. Monte Carlo simulations were carried out based on the MCNPX code in order to investigate and reduce the dose curve perturbation caused by the SCDD-Pro. In particular, various diamond thicknesses were simulated to evaluate the influence of the active volume thickness (e AV) as well as the influence of the addition of a front silver resin (250 µm in thickness in front of the diamond crystal) on depth-dose curves. The simulations indicated that the diamond crystal alone, with a small e AV of just 5 µm, already affects the dose at Bragg peak position (Bragg peak dose) by more than 2% with respect to the Bragg peak dose deposited in water. The optimal design that resulted from the Monte Carlo simulations consists of a diamond crystal of 1 mm in width and 150 µm in thickness with the front silver resin, enclosed by a water-equivalent packaging. This design leads to a deviation between the Bragg peak dose from the full detector modeling and the Bragg peak dose deposited in water of less than 1.2%. Based on those optimizations, an SCDD-Pro prototype was built and evaluated in broad passive scattering proton beams. The experimental evaluation led to probed SCDD-Pro repeatability, dose rate dependence and linearity, that were better than 0.2%, 0.4% (in the 1.0-5.5 Gy min-1 range) and 0.4% (for dose higher than 0.05 Gy), respectively. The depth-dose curves in the 90-160 MeV energy range, measured with the SCDD-Pro without applying any correction, were in good agreement with those measured using a commercial IBA PPC05 plane-parallel ionization chamber, differing by less than 1.6%. The experimental results confirmed that this SCDD-Pro is suitable for measurements with standard electrometers and that the depth-dose curve perturbation is negligible, with no energy dependence and no significant dose rate dependence.
Collapse
Affiliation(s)
- Cyril Moignier
- CEA, LIST, System Modelling and Simulation Laboratory, Gif-sur-Yvette, France. Institut Curie, Centre de Protonthérapie d'Orsay, Orsay, France
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
9
|
Lin L, Huang S, Kang M, Hiltunen P, Vanderstraeten R, Lindberg J, Siljamaki S, Wareing T, Davis I, Barnett A, McGhee J, Simone CB, Solberg TD, McDonough JE, Ainsley C. A benchmarking method to evaluate the accuracy of a commercial proton monte carlo pencil beam scanning treatment planning system. J Appl Clin Med Phys 2017; 18:44-49. [PMID: 28300385 PMCID: PMC5689961 DOI: 10.1002/acm2.12043] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2016] [Accepted: 12/16/2016] [Indexed: 11/12/2022] Open
Abstract
AcurosPT is a Monte Carlo algorithm in the Eclipse 13.7 treatment planning system, which is designed to provide rapid and accurate dose calculations for proton therapy. Computational run-time in minimized by simplifying or eliminating less significant physics processes. In this article, the accuracy of AcurosPT was benchmarked against both measurement and an independent MC calculation, TOPAS. Such a method can be applied to any new MC calculation for the detection of potential inaccuracies. To validate multiple Coulomb scattering (MCS) which affects primary beam broadening, single spot profiles in a Solidwater® phantom were compared for beams of five selected proton energies between AcurosPT, measurement and TOPAS. The spot Gaussian sigma in AcurosPT was found to increase faster with depth than both measurement and TOPAS, suggesting that the MCS algorithm in AcurosPT overestimates the scattering effect. To validate AcurosPT modeling of the halo component beyond primary beam broadening, field size factors (FSF) were compared for multi-spot profiles measured in a water phantom. The FSF for small field sizes were found to disagree with measurement, with the disagreement increasing with depth. Conversely, TOPAS simulations of the same FSF consistently agreed with measurement to within 1.5%. The disagreement in absolute dose between AcurosPT and measurement was smaller than 2% at the mid-range depth of multi-energy beams. While AcurosPT calculates acceptable dose distributions for typical clinical beams, users are cautioned of potentially larger errors at distal depths due to overestimated MCS and halo implementation.
Collapse
Affiliation(s)
- Liyong Lin
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Sheng Huang
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Minglei Kang
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | | | | | | | | | - Todd Wareing
- Varian Medical Systems, Palo Alto, California, USA
| | - Ian Davis
- Varian Medical Systems, Palo Alto, California, USA
| | | | - John McGhee
- Varian Medical Systems, Palo Alto, California, USA
| | - Charles B Simone
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Timothy D Solberg
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - James E McDonough
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Christopher Ainsley
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| |
Collapse
|
10
|
Marsolat F, De Marzi L, Patriarca A, Nauraye C, Moignier C, Pomorski M, Moignau F, Heinrich S, Tromson D, Mazal A. Dosimetric characteristics of four PTW microDiamond detectors in high-energy proton beams. Phys Med Biol 2016; 61:6413-29. [DOI: 10.1088/0031-9155/61/17/6413] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
11
|
Verburg JM, Grassberger C, Dowdell S, Schuemann J, Seco J, Paganetti H. Automated Monte Carlo Simulation of Proton Therapy Treatment Plans. Technol Cancer Res Treat 2016; 15:NP35-NP46. [DOI: 10.1177/1533034615614139] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2015] [Revised: 09/11/2015] [Accepted: 10/01/2015] [Indexed: 11/15/2022] Open
Abstract
Simulations of clinical proton radiotherapy treatment plans using general purpose Monte Carlo codes have been proven to be a valuable tool for basic research and clinical studies. They have been used to benchmark dose calculation methods, to study radiobiological effects, and to develop new technologies such as in vivo range verification methods. Advancements in the availability of computational power have made it feasible to perform such simulations on large sets of patient data, resulting in a need for automated and consistent simulations. A framework called MCAUTO was developed for this purpose. Both passive scattering and pencil beam scanning delivery are supported. The code handles the data exchange between the treatment planning system and the Monte Carlo system, which requires not only transfer of plan and imaging information but also translation of institutional procedures, such as output factor definitions. Simulations are performed on a high-performance computing infrastructure. The simulation methods were designed to use the full capabilities of Monte Carlo physics models, while also ensuring consistency in the approximations that are common to both pencil beam and Monte Carlo dose calculations. Although some methods need to be tailored to institutional planning systems and procedures, the described procedures show a general road map that can be easily translated to other systems.
Collapse
Affiliation(s)
- Joost Mathijs Verburg
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Clemens Grassberger
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Stephen Dowdell
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Jan Schuemann
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Joao Seco
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| | - Harald Paganetti
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
| |
Collapse
|
12
|
Tessonnier T, Marcelos T, Mairani A, Brons S, Parodi K. Phase Space Generation for Proton and Carbon Ion Beams for External Users' Applications at the Heidelberg Ion Therapy Center. Front Oncol 2016; 5:297. [PMID: 26793617 PMCID: PMC4707204 DOI: 10.3389/fonc.2015.00297] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2015] [Accepted: 12/14/2015] [Indexed: 12/04/2022] Open
Abstract
In the field of radiation therapy, accurate and robust dose calculation is required. For this purpose, precise modeling of the irradiation system and reliable computational platforms are needed. At the Heidelberg Ion Therapy Center (HIT), the beamline has been already modeled in the FLUKA Monte Carlo (MC) code. However, this model was kept confidential for disclosure reasons and was not available for any external team. The main goal of this study was to create efficiently phase space (PS) files for proton and carbon ion beams, for all energies and foci available at HIT. PSs are representing the characteristics of each particle recorded (charge, mass, energy, coordinates, direction cosines, generation) at a certain position along the beam path. In order to achieve this goal, keeping a reasonable data size but maintaining the requested accuracy for the calculation, we developed a new approach of beam PS generation with the MC code FLUKA. The generated PSs were obtained using an infinitely narrow beam and recording the desired quantities after the last element of the beamline, with a discrimination of primaries or secondaries. In this way, a unique PS can be used for each energy to accommodate the different foci by combining the narrow-beam scenario with a random sampling of its theoretical Gaussian beam in vacuum. PS can also reproduce the different patterns from the delivery system, when properly combined with the beam scanning information. MC simulations using PS have been compared to simulations, including the full beamline geometry and have been found in very good agreement for several cases (depth dose distributions, lateral dose profiles), with relative dose differences below 0.5%. This approach has also been compared with measured data of ion beams with different energies and foci, resulting in a very satisfactory agreement. Hence, the proposed approach was able to fulfill the different requirements and has demonstrated its capability for application to clinical treatment fields. It also offers a powerful tool to perform investigations on the contribution of primary and secondary particles produced in the beamline. These PSs are already made available to external teams upon request, to support interpretation of their measurements.
Collapse
Affiliation(s)
- Thomas Tessonnier
- Department of Radiation Oncology, Heidelberg University Clinic, Heidelberg, Germany; Department of Medical Physics, Ludwig Maximilians University, Munich, Germany
| | - Tiago Marcelos
- Department of Medical Physics, Ludwig Maximilians University , Munich , Germany
| | - Andrea Mairani
- Heidelberg Ion Beam Therapy Center, Heidelberg, Germany; Centro Nazionale di Adroterapia Oncologica, Pavia, Italy
| | - Stephan Brons
- Heidelberg Ion Beam Therapy Center , Heidelberg , Germany
| | - Katia Parodi
- Department of Medical Physics, Ludwig Maximilians University, Munich, Germany; Heidelberg Ion Beam Therapy Center, Heidelberg, Germany
| |
Collapse
|
13
|
Titt U, Sell M, Unkelbach J, Bangert M, Mirkovic D, Oelfke U, Mohan R. Degradation of proton depth dose distributions attributable to microstructures in lung-equivalent material. Med Phys 2015; 42:6425-32. [PMID: 26520732 PMCID: PMC4608968 DOI: 10.1118/1.4932625] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Revised: 08/28/2015] [Accepted: 09/24/2015] [Indexed: 11/07/2022] Open
Abstract
PURPOSE The purpose of the work reported here was to investigate the influence of sub-millimeter size heterogeneities on the degradation of the distal edges of proton beams and to validate Monte Carlo (MC) methods' ability to correctly predict such degradation. METHODS A custom-designed high-resolution plastic phantom approximating highly heterogeneous, lung-like structures was employed in measurements and in Monte Carlo simulations to evaluate the degradation of proton Bragg curves penetrating heterogeneous media. RESULTS Significant differences in distal falloff widths and in peak dose values were observed in the measured and the Monte Carlo simulated curves compared to pristine proton Bragg curves. Furthermore, differences between simulations of beams penetrating CT images of the phantom did not agree well with the corresponding experimental differences. The distal falloff widths in CT image-based geometries were underestimated by up to 0.2 cm in water (corresponding to 0.8-1.4 cm in lung tissue), and the peak dose values of pristine proton beams were overestimated by as much as ˜35% compared to measured curves or depth-dose curves simulated on the basis of true geometry. The authors demonstrate that these discrepancies were caused by the limited spatial resolution of CT images that served as a basis for dose calculations and lead to underestimation of the impact of the fine structure of tissue heterogeneities. A convolution model was successfully applied to mitigate the underestimation. CONCLUSIONS The results of this study justify further development of models to better represent heterogeneity effects in soft-tissue geometries, such as lung, and to correct systematic underestimation of the degradation of the distal edge of proton doses.
Collapse
Affiliation(s)
- Uwe Titt
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030
| | - Martin Sell
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030 and Department of Medical Physics, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg 69120, Germany
| | - Jan Unkelbach
- Department of Radiation Oncology, Massachusetts General Hospital, 55 Fruit Street, Boston, Massachusetts 02114
| | - Mark Bangert
- Department of Medical Physics, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg 69120, Germany
| | - Dragan Mirkovic
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030
| | - Uwe Oelfke
- Department of Medical Physics, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg 69120, Germany and Department of Physics, The Institute of Cancer Research, 123 Old Brompton Road, London SW7 3RP, United Kingdom
| | - Radhe Mohan
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, Texas 77030
| |
Collapse
|
14
|
Fracchiolla F, Lorentini S, Widesott L, Schwarz M. Characterization and validation of a Monte Carlo code for independent dose calculation in proton therapy treatments with pencil beam scanning. Phys Med Biol 2015; 60:8601-19. [PMID: 26501569 DOI: 10.1088/0031-9155/60/21/8601] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We propose a method of creating and validating a Monte Carlo (MC) model of a proton Pencil Beam Scanning (PBS) machine using only commissioning measurements and avoiding the nozzle modeling. Measurements with a scintillating screen coupled with a CCD camera, ionization chamber and a Faraday Cup were used to model the beam in TOPAS without using any machine parameter information but the virtual source distance from the isocenter. Then the model was validated on simple Spread Out Bragg Peaks (SOBP) delivered in water phantom and with six realistic clinical plans (many involving 3 or more fields) on an anthropomorphic phantom. In particular the behavior of the moveable Range Shifter (RS) feature was investigated and its modeling has been proposed. The gamma analysis (3%,3 mm) was used to compare MC, TPS (XiO-ELEKTA) and measured 2D dose distributions (using radiochromic film). The MC modeling proposed here shows good results in the validation phase, both for simple irradiation geometry (SOBP in water) and for modulated treatment fields (on anthropomorphic phantoms). In particular head lesions were investigated and both MC and TPS data were compared with measurements. Treatment plans with no RS always showed a very good agreement with both of them (γ-Passing Rate (PR) > 95%). Treatment plans in which the RS was needed were also tested and validated. For these treatment plans MC results showed better agreement with measurements (γ-PR > 93%) than the one coming from TPS (γ-PR < 88%). This work shows how to simplify the MC modeling of a PBS machine for proton therapy treatments without accounting for any hardware components and proposes a more reliable RS modeling than the one implemented in our TPS. The validation process has shown how this code is a valid candidate for a completely independent treatment plan dose calculation algorithm. This makes the code an important future tool for the patient specific QA verification process.
Collapse
Affiliation(s)
- F Fracchiolla
- Azienda Provinciale per i Servizi Sanitari (APSS) Protontherapy Department, Trento, Italy. Post Graduate School of Medical Physics 'Sapienza' University of Rome, 00185 Roma, Italy
| | | | | | | |
Collapse
|
15
|
Abstract
The physics of proton therapy has advanced considerably since it was proposed in 1946. Today analytical equations and numerical simulation methods are available to predict and characterize many aspects of proton therapy. This article reviews the basic aspects of the physics of proton therapy, including proton interaction mechanisms, proton transport calculations, the determination of dose from therapeutic and stray radiations, and shielding design. The article discusses underlying processes as well as selected practical experimental and theoretical methods. We conclude by briefly speculating on possible future areas of research of relevance to the physics of proton therapy.
Collapse
Affiliation(s)
- Wayne D Newhauser
- Medical Physics Program, Department of Physics and Astronomy, Louisiana State University, 202 Nicholson Hall, Baton Rouge, LA, 70803, USA
- Mary Bird Perkins Cancer Center, 4950 Essen Lane, Baton Rouge, LA, 70809, USA
| | - Rui Zhang
- Medical Physics Program, Department of Physics and Astronomy, Louisiana State University, 202 Nicholson Hall, Baton Rouge, LA, 70803, USA
- Mary Bird Perkins Cancer Center, 4950 Essen Lane, Baton Rouge, LA, 70809, USA
| |
Collapse
|
16
|
Inter-Institutional Comparison of Personalized Risk Assessments for Second Malignant Neoplasms for a 13-Year-Old Girl Receiving Proton versus Photon Craniospinal Irradiation. Cancers (Basel) 2015; 7:407-26. [PMID: 25763928 PMCID: PMC4381265 DOI: 10.3390/cancers7010407] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2014] [Revised: 02/18/2015] [Accepted: 02/19/2015] [Indexed: 01/14/2023] Open
Abstract
Children receiving radiotherapy face the probability of a subsequent malignant neoplasm (SMN). In some cases, the predicted SMN risk can be reduced by proton therapy. The purpose of this study was to apply the most comprehensive dose assessment methods to estimate the reduction in SMN risk after proton therapy vs. photon therapy for a 13-year-old girl requiring craniospinal irradiation (CSI). We reconstructed the equivalent dose throughout the patient’s body from therapeutic and stray radiation and applied SMN incidence and mortality risk models for each modality. Excluding skin cancer, the risk of incidence after proton CSI was a third of that of photon CSI. The predicted absolute SMN risks were high. For photon CSI, the SMN incidence rates greater than 10% were for thyroid, non-melanoma skin, lung, colon, stomach, and other solid cancers, and for proton CSI they were non-melanoma skin, lung, and other solid cancers. In each setting, lung cancer accounted for half the risk of mortality. In conclusion, the predicted SMN risk for a 13-year-old girl undergoing proton CSI was reduced vs. photon CSI. This study demonstrates the feasibility of inter-institutional whole-body dose and risk assessments and also serves as a model for including risk estimation in personalized cancer care.
Collapse
|
17
|
Hyer DE, Hill PM, Wang D, Smith BR, Flynn RT. Effects of spot size and spot spacing on lateral penumbra reduction when using a dynamic collimation system for spot scanning proton therapy. Phys Med Biol 2014; 59:N187-96. [PMID: 25330783 DOI: 10.1088/0031-9155/59/22/n187] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The purpose of this work was to investigate the reduction in lateral dose penumbra that can be achieved when using a dynamic collimation system (DCS) for spot scanning proton therapy as a function of two beam parameters: spot size and spot spacing. This is an important investigation as both values impact the achievable dose distribution and a wide range of values currently exist depending on delivery hardware. Treatment plans were created both with and without the DCS for in-air spot sizes (σair) of 3, 5, 7, and 9 mm as well as spot spacing intervals of 2, 4, 6 and 8 mm. Compared to un-collimated treatment plans, the plans created with the DCS yielded a reduction in the mean dose to normal tissue surrounding the target of 26.2-40.6% for spot sizes of 3-9 mm, respectively. Increasing the spot spacing resulted in a decrease in the time penalty associated with using the DCS that was approximately proportional to the reduction in the number of rows in the raster delivery pattern. We conclude that dose distributions achievable when using the DCS are comparable to those only attainable with much smaller initial spot sizes, suggesting that the goal of improving high dose conformity may be achieved by either utilizing a DCS or by improving beam line optics.
Collapse
Affiliation(s)
- Daniel E Hyer
- Department of Radiation Oncology, University of Iowa, 200 Hawkins Drive, Iowa City, Iowa 52242, USA
| | | | | | | | | |
Collapse
|
18
|
Carnicer A, Angellier G, Thariat J, Sauerwein W, Caujolle JP, Hérault J. Quantification of dose perturbations induced by external and internal accessories in ocular proton therapy and evaluation of their dosimetric impact. Med Phys 2014; 40:061708. [PMID: 23718587 DOI: 10.1118/1.4807090] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Proton scattering on beam shaping devices and protons slowing down on media with different densities within the treatment volume may produce dose perturbations and range variations that are not predicted by treatment planning systems. The aim of this work was to assess the dosimetric impact of elements present in ocular proton therapy treatments that may disturb the prescribed treatment plan. Both distal beam shaping devices and intraocular elements were considered. METHODS A wedge filter, tantalum fiducial marker, hemispherical compensator, two intraocular endotamponades (densities 0.97 and 1.92 g cm(-3)) and an intraocular eye lens (IOL) were considered in the study. For these elements, longitudinal dose distributions were measured and∕or calculated in water in beam alignment for a clinical spread-out Bragg peak. Under the same conditions, the unperturbed dose distributions were similarly measured and∕or calculated in the absence of the element. The dosimetric impact was assessed by comparison of unperturbed and perturbed dose distributions. Measurements and calculations were carried out with two methods. Measurements are based on EBT3 films with dedicated software, which makes use of a calibration curve and correction for the quenching effect. Calculations are based on the Monte Carlo (MC) code MCNPX and reproduce the experimental conditions. Both dose maps are obtained with a resolution of 300 dpi. RESULTS The degree of disturbance of distal beam shaping devices is low for the wedge filter (2% overdose ripple all along the central axis) and moderate for the hemispherical compensator (two bands of variable overdose of up to 10% downstream the compensator lateral edges and -5% underdose on the plateau at off-axis distance of 5 cm). Tantalum clips produce important dose shadows (-20% behind the clip parallel to the beam and range reduction of 1.1 mm) and bands of overdose (15%). The presence of endotamponades modifies the dose distribution very significantly (-5% underdose on the plateau and 3 mm range prolongation for the tamponade with density 0.97 g cm(-3) and -15% underdose on plateau and 8 mm range reduction for that with density 1.92 g cm(-3)). No dose perturbations were found for the IOL. The high performance of EBT3 film and MC tools used was confirmed and good agreement was found between them (percentage of pixels passing the gamma test >87%). CONCLUSIONS The degree of disturbance by external beam shaping devices remains low in ocular proton therapy and can be reduced by bringing accessories closer to the eye. Tantalum fiducial markers must be located in such a way that dose perturbation is not projected on the tumor. The treatment of patients with intraocular endotamponades must be carefully managed. It is fundamental that radiation oncologists and medical physicists are informed about the presence of such substances prior to the treatment.
Collapse
Affiliation(s)
- A Carnicer
- Centre Antoine Lacassagne, Cyclotron Biomédical, 227 Avenue de la Lanterne, 06200 Nice, France
| | | | | | | | | | | |
Collapse
|
19
|
Carnicer A, Angellier G, Gérard A, Garnier N, Dubois C, Amblard R, Hérault J. Development and validation of radiochromic film dosimetry and Monte Carlo simulation tools for acquisition of absolute, high-spatial resolution longitudinal dose distributions in ocular proton therapy. RADIAT MEAS 2013. [DOI: 10.1016/j.radmeas.2013.06.011] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
20
|
Newhauser WD, Rechner L, Mirkovic D, Yepes P, Koch NC, Titt U, Fontenot JD, Zhang R. Benchmark measurements and simulations of dose perturbations due to metallic spheres in proton beams. RADIAT MEAS 2013; 58:37-44. [PMID: 25147474 PMCID: PMC4136527 DOI: 10.1016/j.radmeas.2013.08.001] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
Monte Carlo simulations are increasingly used for dose calculations in proton therapy due to its inherent accuracy. However, dosimetric deviations have been found using Monte Carlo code when high density materials are present in the proton beam line. The purpose of this work was to quantify the magnitude of dose perturbation caused by metal objects. We did this by comparing measurements and Monte Carlo predictions of dose perturbations caused by the presence of small metal spheres in several clinical proton therapy beams as functions of proton beam range, spread-out Bragg peak width and drift space. Monte Carlo codes MCNPX, GEANT4 and Fast Dose Calculator (FDC) were used. Generally good agreement was found between measurements and Monte Carlo predictions, with the average difference within 5% and maximum difference within 17%. The modification of multiple Coulomb scattering model in MCNPX code yielded improvement in accuracy and provided the best overall agreement with measurements. Our results confirmed that Monte Carlo codes are well suited for predicting multiple Coulomb scattering in proton therapy beams when short drift spaces are involved.
Collapse
Affiliation(s)
- Wayne D. Newhauser
- The University of Texas M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Unit 94, Houston, Texas 77030, USA
- The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner, Houston, TX 77030, USA
- Department of Medical Physics, Mary Bird Perkins Cancer Center, Baton Rouge, LA, 70809, USA
| | - Laura Rechner
- The University of Texas M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Unit 94, Houston, Texas 77030, USA
- The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner, Houston, TX 77030, USA
| | - Dragan Mirkovic
- The University of Texas M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Unit 94, Houston, Texas 77030, USA
- The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner, Houston, TX 77030, USA
| | - Pablo Yepes
- Department of Physics and Astronomy, MS 315, Rice University, 6100 Main Street, Houston, TX 77005, USA
| | - Nicholas C. Koch
- Department of Radiation Oncology, Medical University of South Carolina, 169 Ashley Avenue, Charleston, SC, 29425, USA
| | - Uwe Titt
- The University of Texas M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Unit 94, Houston, Texas 77030, USA
- The University of Texas Graduate School of Biomedical Sciences at Houston, 6767 Bertner, Houston, TX 77030, USA
| | - Jonas D. Fontenot
- Department of Medical Physics, Mary Bird Perkins Cancer Center, Baton Rouge, LA, 70809, USA
| | - Rui Zhang
- Department of Medical Physics, Mary Bird Perkins Cancer Center, Baton Rouge, LA, 70809, USA
| |
Collapse
|
21
|
Perl J, Shin J, Schumann J, Faddegon B, Paganetti H. TOPAS: an innovative proton Monte Carlo platform for research and clinical applications. Med Phys 2013; 39:6818-37. [PMID: 23127075 DOI: 10.1118/1.4758060] [Citation(s) in RCA: 573] [Impact Index Per Article: 52.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE While Monte Carlo particle transport has proven useful in many areas (treatment head design, dose calculation, shielding design, and imaging studies) and has been particularly important for proton therapy (due to the conformal dose distributions and a finite beam range in the patient), the available general purpose Monte Carlo codes in proton therapy have been overly complex for most clinical medical physicists. The learning process has large costs not only in time but also in reliability. To address this issue, we developed an innovative proton Monte Carlo platform and tested the tool in a variety of proton therapy applications. METHODS Our approach was to take one of the already-established general purpose Monte Carlo codes and wrap and extend it to create a specialized user-friendly tool for proton therapy. The resulting tool, TOol for PArticle Simulation (TOPAS), should make Monte Carlo simulation more readily available for research and clinical physicists. TOPAS can model a passive scattering or scanning beam treatment head, model a patient geometry based on computed tomography (CT) images, score dose, fluence, etc., save and restart a phase space, provides advanced graphics, and is fully four-dimensional (4D) to handle variations in beam delivery and patient geometry during treatment. A custom-designed TOPAS parameter control system was placed at the heart of the code to meet requirements for ease of use, reliability, and repeatability without sacrificing flexibility. RESULTS We built and tested the TOPAS code. We have shown that the TOPAS parameter system provides easy yet flexible control over all key simulation areas such as geometry setup, particle source setup, scoring setup, etc. Through design consistency, we have insured that user experience gained in configuring one component, scorer or filter applies equally well to configuring any other component, scorer or filter. We have incorporated key lessons from safety management, proactively removing possible sources of user error such as line-ordering mistakes. We have modeled proton therapy treatment examples including the UCSF eye treatment head, the MGH stereotactic alignment in radiosurgery treatment head and the MGH gantry treatment heads in passive scattering and scanning modes, and we have demonstrated dose calculation based on patient-specific CT data. Initial validation results show agreement with measured data and demonstrate the capabilities of TOPAS in simulating beam delivery in 3D and 4D. CONCLUSIONS We have demonstrated TOPAS accuracy and usability in a variety of proton therapy setups. As we are preparing to make this tool freely available for researchers in medical physics, we anticipate widespread use of this tool in the growing proton therapy community.
Collapse
Affiliation(s)
- J Perl
- SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA.
| | | | | | | | | |
Collapse
|
22
|
Yamashita T, Akagi T, Aso T, Kimura A, Sasaki T. Effect of inhomogeneity in a patient's body on the accuracy of the pencil beam algorithm in comparison to Monte Carlo. Phys Med Biol 2012; 57:7673-88. [PMID: 23123683 DOI: 10.1088/0031-9155/57/22/7673] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The pencil beam algorithm (PBA) is reasonably accurate and fast. It is, therefore, the primary method used in routine clinical treatment planning for proton radiotherapy; still, it needs to be validated for use in highly inhomogeneous regions. In our investigation of the effect of patient inhomogeneity, PBA was compared with Monte Carlo (MC). A software framework was developed for the MC simulation of radiotherapy based on Geant4. Anatomical sites selected for the comparison were the head/neck, liver, lung and pelvis region. The dose distributions calculated by the two methods in selected examples were compared, as well as a dose volume histogram (DVH) derived from the dose distributions. The comparison of the off-center ratio (OCR) at the iso-center showed good agreement between the PBA and MC, while discrepancies were seen around the distal fall-off regions. While MC showed a fine structure on the OCR in the distal fall-off region, the PBA showed smoother distribution. The fine structures in MC calculation appeared downstream of very low-density regions. Comparison of DVHs showed that most of the target volumes were similarly covered, while some OARs located around the distal region received a higher dose when calculated by MC than the PBA.
Collapse
Affiliation(s)
- T Yamashita
- Hyogo Ion Beam Medical Center, Shingu-cho, Tatsuno, Hyogo, Japan.
| | | | | | | | | |
Collapse
|
23
|
Jia Y, Beltran C, Indelicato DJ, Flampouri S, Li Z, Merchant TE. Proton therapy dose distribution comparison between Monte Carlo and a treatment planning system for pediatric patients with ependymomaa). Med Phys 2012; 39:4742-7. [DOI: 10.1118/1.4736413] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
24
|
Abstract
The main advantages of proton therapy are the reduced total energy deposited in the patient as compared to photon techniques and the finite range of the proton beam. The latter adds an additional degree of freedom to treatment planning. The range in tissue is associated with considerable uncertainties caused by imaging, patient setup, beam delivery and dose calculation. Reducing the uncertainties would allow a reduction of the treatment volume and thus allow a better utilization of the advantages of protons. This paper summarizes the role of Monte Carlo simulations when aiming at a reduction of range uncertainties in proton therapy. Differences in dose calculation when comparing Monte Carlo with analytical algorithms are analyzed as well as range uncertainties due to material constants and CT conversion. Range uncertainties due to biological effects and the role of Monte Carlo for in vivo range verification are discussed. Furthermore, the current range uncertainty recipes used at several proton therapy facilities are revisited. We conclude that a significant impact of Monte Carlo dose calculation can be expected in complex geometries where local range uncertainties due to multiple Coulomb scattering will reduce the accuracy of analytical algorithms. In these cases Monte Carlo techniques might reduce the range uncertainty by several mm.
Collapse
Affiliation(s)
- Harald Paganetti
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA.
| |
Collapse
|
25
|
Paganetti H. Range uncertainties in proton therapy and the role of Monte Carlo simulations. Phys Med Biol 2012. [PMID: 22571913 DOI: 10.1088/0031‐9155/57/11/r99] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The main advantages of proton therapy are the reduced total energy deposited in the patient as compared to photon techniques and the finite range of the proton beam. The latter adds an additional degree of freedom to treatment planning. The range in tissue is associated with considerable uncertainties caused by imaging, patient setup, beam delivery and dose calculation. Reducing the uncertainties would allow a reduction of the treatment volume and thus allow a better utilization of the advantages of protons. This paper summarizes the role of Monte Carlo simulations when aiming at a reduction of range uncertainties in proton therapy. Differences in dose calculation when comparing Monte Carlo with analytical algorithms are analyzed as well as range uncertainties due to material constants and CT conversion. Range uncertainties due to biological effects and the role of Monte Carlo for in vivo range verification are discussed. Furthermore, the current range uncertainty recipes used at several proton therapy facilities are revisited. We conclude that a significant impact of Monte Carlo dose calculation can be expected in complex geometries where local range uncertainties due to multiple Coulomb scattering will reduce the accuracy of analytical algorithms. In these cases Monte Carlo techniques might reduce the range uncertainty by several mm.
Collapse
Affiliation(s)
- Harald Paganetti
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA 02114, USA.
| |
Collapse
|
26
|
Smeets J, Roellinghoff F, Prieels D, Stichelbaut F, Benilov A, Busca P, Fiorini C, Peloso R, Basilavecchia M, Frizzi T, Dehaes JC, Dubus A. Prompt gamma imaging with a slit camera for real-time range control in proton therapy. Phys Med Biol 2012; 57:3371-405. [DOI: 10.1088/0031-9155/57/11/3371] [Citation(s) in RCA: 238] [Impact Index Per Article: 19.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
27
|
Seravalli E, Robert C, Bauer J, Stichelbaut F, Kurz C, Smeets J, Van Ngoc Ty C, Schaart DR, Buvat I, Parodi K, Verhaegen F. Monte Carlo calculations of positron emitter yields in proton radiotherapy. Phys Med Biol 2012; 57:1659-73. [PMID: 22398196 DOI: 10.1088/0031-9155/57/6/1659] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
28
|
Garcia-Molina R, Abril I, Heredia-Avalos S, Kyriakou I, Emfietzoglou D. A combined molecular dynamics and Monte Carlo simulation of the spatial distribution of energy deposition by proton beams in liquid water. Phys Med Biol 2011; 56:6475-93. [DOI: 10.1088/0031-9155/56/19/019] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
29
|
Grevillot L, Bertrand D, Dessy F, Freud N, Sarrut D. A Monte Carlo pencil beam scanning model for proton treatment plan simulation using GATE/GEANT4. Phys Med Biol 2011; 56:5203-19. [PMID: 21791731 DOI: 10.1088/0031-9155/56/16/008] [Citation(s) in RCA: 123] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
This work proposes a generic method for modeling scanned ion beam delivery systems, without simulation of the treatment nozzle and based exclusively on beam data library (BDL) measurements required for treatment planning systems (TPS). To this aim, new tools dedicated to treatment plan simulation were implemented in the Gate Monte Carlo platform. The method was applied to a dedicated nozzle from IBA for proton pencil beam scanning delivery. Optical and energy parameters of the system were modeled using a set of proton depth-dose profiles and spot sizes measured at 27 therapeutic energies. For further validation of the beam model, specific 2D and 3D plans were produced and then measured with appropriate dosimetric tools. Dose contributions from secondary particles produced by nuclear interactions were also investigated using field size factor experiments. Pristine Bragg peaks were reproduced with 0.7 mm range and 0.2 mm spot size accuracy. A 32 cm range spread-out Bragg peak with 10 cm modulation was reproduced with 0.8 mm range accuracy and a maximum point-to-point dose difference of less than 2%. A 2D test pattern consisting of a combination of homogeneous and high-gradient dose regions passed a 2%/2 mm gamma index comparison for 97% of the points. In conclusion, the generic modeling method proposed for scanned ion beam delivery systems was applicable to an IBA proton therapy system. The key advantage of the method is that it only requires BDL measurements of the system. The validation tests performed so far demonstrated that the beam model achieves clinical performance, paving the way for further studies toward TPS benchmarking. The method involves new sources that are available in the new Gate release V6.1 and could be further applied to other particle therapy systems delivering protons or other types of ions like carbon.
Collapse
Affiliation(s)
- L Grevillot
- Université de Lyon, CREATIS, CNRS UMR5220, Inserm U1044, INSA-Lyon, Université Lyon 1, Centre Léon Bérard, Lyon, France.
| | | | | | | | | |
Collapse
|
30
|
Titt U, Mirkovic D, Sawakuchi GO, Perles LA, Newhauser WD, Taddei PJ, Mohan R. Adjustment of the lateral and longitudinal size of scanned proton beam spots using a pre-absorber to optimize penumbrae and delivery efficiency. Phys Med Biol 2010; 55:7097-106. [PMID: 21076194 DOI: 10.1088/0031-9155/55/23/s10] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
In scanned-beam proton therapy, the beam spot properties, such as the lateral and longitudinal size and the minimum achievable range, are influenced by beam optics, scattering media and drift spaces in the treatment unit. Currently available spot scanning systems offer few options for adjusting these properties. We investigated a method for adjusting the lateral and longitudinal spot size that utilizes downstream plastic pre-absorbers located near a water phantom. The spot size adjustment was characterized using Monte Carlo simulations of a modified commercial scanned-beam treatment head. Our results revealed that the pre-absorbers can be used to reduce the lateral full width at half maximum (FWHM) of dose spots in water by up to 14 mm, and to increase the longitudinal extent from about 1 mm to 5 mm at residual ranges of 4 cm and less. A large factor in manipulating the lateral spot sizes is the drift space between the pre-absorber and the water phantom. Increasing the drift space from 0 cm to 15 cm leads to an increase in the lateral FWHM from 2.15 cm to 2.87 cm, at a water-equivalent depth of 1 cm. These findings suggest that this spot adjustment method may improve the quality of spot-scanned proton treatments.
Collapse
Affiliation(s)
- Uwe Titt
- Department of Radiation Physics, The University of Texas, M D Anderson Cancer Center, Houston, TX 77030, USA.
| | | | | | | | | | | | | |
Collapse
|
31
|
Sawakuchi GO, Mirkovic D, Perles LA, Sahoo N, Zhu XR, Ciangaru G, Suzuki K, Gillin MT, Mohan R, Titt U. An MCNPX Monte Carlo model of a discrete spot scanning proton beam therapy nozzle. Med Phys 2010; 37:4960-70. [PMID: 20964215 DOI: 10.1118/1.3476458] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE The purposes of this study were to validate a discrete spot scanning proton beam nozzle using the Monte Carlo (MC) code MCNPX and use the MC validated model to investigate the effects of a low-dose envelope, which surrounds the beam's central axis, on measurements of integral depth dose (IDD) profiles. METHODS An accurate model of the discrete spot scanning beam nozzle from The University of Texas M. D. Anderson Cancer Center (Houston, Texas) was developed on the basis of blueprints provided by the manufacturer of the nozzle. The authors performed simulations of single proton pencil beams of various energies using the standard multiple Coulomb scattering (MCS) algorithm within the MCNPX source code and a new MCS algorithm, which was implemented in the MCNPX source code. The MC models were validated by comparing calculated in-air and in-water lateral profiles and percentage depth dose profiles for single pencil beams with their corresponding measured values. The models were then further tested by comparing the calculated and measured three-dimensional (3-D) dose distributions. Finally, an IDD profile was calculated with different scoring radii to determine the limitations on the use of commercially available plane-parallel ionization chambers to measure IDD. RESULTS The distance to agreement, defined as the distance between the nearest positions of two equivalent distributions with the same value of dose, between measured and simulated ranges was within 0.13 cm for both MCS algorithms. For low and intermediate pencil beam energies, the MC simulations using the standard MCS algorithm were in better agreement with measurements. Conversely, the new MCS algorithm produced better results for high-energy single pencil beams. The IDD profile calculated with cylindrical tallies with an area equivalent to the area of the largest commercially available ionization chamber showed up to 7.8% underestimation of the integral dose in certain depths of the IDD profile. CONCLUSIONS The authors conclude that a combination of MCS algorithms is required to accurately reproduce experimental data of single pencil beams and 3-D dose distributions for the scanning beam nozzle. In addition, the MC simulations showed that because of the low-dose envelope, ionization chambers with radii as large as 4.08 cm are insufficient to accurately measure IDD profiles for a 221.8 MeV pencil beam in the scanning beam nozzle.
Collapse
Affiliation(s)
- Gabriel O Sawakuchi
- Department of Physics, Carleton University, 1125 Colonel By Drive, Ottawa, Ontario K1S 5B6, Canada.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
32
|
Sawakuchi GO, Titt U, Mirkovic D, Ciangaru G, Zhu XR, Sahoo N, Gillin MT, Mohan R. Monte Carlo investigation of the low-dose envelope from scanned proton pencil beams. Phys Med Biol 2010; 55:711-21. [PMID: 20071752 DOI: 10.1088/0031-9155/55/3/011] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Scanned proton pencil beams carry a low-dose envelope that extends several centimeters from the individual beam's central axis. Thus, the total delivered dose depends on the size of the target volume and the corresponding number and intensity of beams necessary to cover the target volume uniformly. This dependence must be considered in dose calculation algorithms used by treatment planning systems. In this work, we investigated the sources of particles contributing to the low-dose envelope using the Monte Carlo technique. We used a validated model of our institution's scanning beam line to determine the contributions to the low-dose envelope from secondary particles created in a water phantom and particles scattered in beam line components. Our results suggested that, for high-energy beams, secondary particles produced by nuclear interactions in the water phantom are the major contributors to the low-dose envelope. For low-energy beams, the low-dose envelope is dominated by particles undergoing multiple Coulomb scattering in the beam line components and water phantom. Clearly, in the latter situation, the low-dose envelope depends directly on beam line design features. Finally, we investigated the dosimetric consequences of the low-dose envelope. Our results showed that if not modeled properly the low-dose envelope may cause clinically relevant dose disturbance in the target volume. This work suggested that this low-dose envelope is beam line specific for low-energy beams, should be thoroughly experimentally characterized and validated during commissioning of the treatment planning system, and therefore is of great concern for accurate delivery of proton scanning beam doses.
Collapse
Affiliation(s)
- Gabriel O Sawakuchi
- Department of Radiation Physics, The University of Texas M D Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, TX 77030, USA.
| | | | | | | | | | | | | | | |
Collapse
|