1
|
Younkin JE, Morales DH, Shen J, Ding X, Stoker JB, Yu NY, Sio TT, Daniels TB, Bues M, Fatyga M, Schild SE, Liu W. Technical Note: Multiple energy extraction techniques for synchrotron-based proton delivery systems may exacerbate motion interplay effects in lung cancer treatments. Med Phys 2021; 48:4812-4823. [PMID: 34174087 DOI: 10.1002/mp.15056] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Revised: 03/12/2021] [Accepted: 06/09/2021] [Indexed: 11/07/2022] Open
Abstract
PURPOSE The multiple energy extraction (MEE) delivery technique for synchrotron-based proton delivery systems reduces beam delivery time by decelerating the beam multiple times during one accelerator spill, but this might cause additional plan quality degradation due to intrafractional motion. We seek to determine whether MEE causes significantly different plan quality degradation compared to single energy extraction (SEE) for lung cancer treatments due to the interplay effect. METHODS Ten lung cancer patients treated with IMPT at our institution were nonrandomly sampled based on a representative range of tumor motion amplitudes, tumor volumes, and respiratory periods. Dose-volume histogram (DVH) indices from single-fraction SEE and MEE four-dimensional (4D) dynamic dose distributions were compared using the Wilcoxon signed-rank test. Distributions of monitor units (MU) to breathing phases were investigated for features associated with plan quality degradation. SEE and MEE DVH indices were compared in fractionated deliveries of the worst-case patient treatment scenario to evaluate the impact of fractionation. RESULTS There were no clinically significant differences in target mean dose, target dose conformity, or dose to organs-at-risk between SEE and MEE in single-fraction delivery. Three patients had significantly worse dose homogeneity with MEE compared to SEE (single-fraction mean D5% -D95% increased by up to 9.6% of prescription dose), and plots of MU distribution to breathing phases showed synchronization patterns with MEE but not SEE. However, after 30 fractions the patient in the worst-case scenario had clinically acceptable target dose homogeneity and coverage with MEE (mean D5% -D95% increased by 1% compared to SEE). CONCLUSIONS For some patients with breathing periods close to the mean spill duration, MEE resulted in significantly worse single-fraction target dose homogeneity compared to SEE due to the interplay effect. However, this was mitigated by fractionation, and target dose homogeneity and coverage were clinically acceptable after 30 fractions with MEE.
Collapse
Affiliation(s)
- James E Younkin
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | | | - Jiajian Shen
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Xiaoning Ding
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Joshua B Stoker
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Nathan Y Yu
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Terence T Sio
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Thomas B Daniels
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Martin Bues
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Mirek Fatyga
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Steven E Schild
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| | - Wei Liu
- Department of Radiation Oncology, Mayo Clinic Arizona, Phoenix, AZ, USA
| |
Collapse
|
2
|
Alina G, Krieger M, Jud C, Duetschler A, Salomir R, Bieri O, Bauman G, Nguyen D, Weber DC, Lomax AJ, Zhang Y, Cattin PC. Liver-ultrasound based motion modelling to estimate 4D dose distributions for lung tumours in scanned proton therapy. ACTA ACUST UNITED AC 2020; 65:235050. [DOI: 10.1088/1361-6560/abaa26] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
3
|
Smith BR, Hyer DE, Flynn RT, Hill PM, Culberson WS. Trimmer sequencing time minimization during dynamically collimated proton therapy using a colony of cooperating agents. Phys Med Biol 2019; 64:205025. [PMID: 31484170 DOI: 10.1088/1361-6560/ab416d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
The dynamic collimation system (DCS) can be combined with pencil beam scanning proton therapy to deliver highly conformal treatment plans with unique collimation at each energy layer. This energy layer-specific collimation is accomplished through the synchronized motion of four trimmer blades that intercept the proton beam near the target boundary in the beam's eye view. However, the corresponding treatment deliveries come at the cost of additional treatment time since the translational speed of the trimmer is slower than the scanning speed of the proton pencil beam. In an attempt to minimize the additional trimmer sequencing time of each field while still maintaining a high degree of conformity, a novel process utilizing ant colony optimization (ACO) methods was created to determine the most efficient route of trimmer sequencing and beamlet scanning patterns for a collective set of collimated proton beamlets. The ACO process was integrated within an in-house treatment planning system optimizer to determine the beam scanning and DCS trimmer sequencing patterns and compared against an analytical approximation of the trimmer sequencing time should a contour-like scanning approach be assumed instead. Due to the stochastic nature of ACO, parameters where determined so that they could ensure good convergence and an efficient optimization of trimmer sequencing that was faster than an analytical contour-like trimmer sequencing. The optimization process was tested using a set of three intracranial treatment plans which were planned using a custom research treatment planning system and were successfully optimized to reduce the additional trimmer sequencing time to approximately 60 s per treatment field while maintaining a high degree of target conformity. Thus, the novel use of ACO techniques within a treatment planning algorithm has been demonstrated to effectively determine collimation sequencing patterns for a DCS in order to minimize the additional treatment time required for trimmer movement during treatment.
Collapse
Affiliation(s)
- Blake R Smith
- Department of Medical Physics, School of Medicine and Public Health, University of Wisconsin-Madison, Madison, Wisconsin 53705. Author to whom correspondence should be addressed
| | | | | | | | | |
Collapse
|
4
|
Boria AJ, Uh J, Pirlepesov F, Stuckey JC, Axente M, Gargone MA, Hua CH. Interplay Effect of Target Motion and Pencil-Beam Scanning in Proton Therapy for Pediatric Patients. Int J Part Ther 2018; 5:1-10. [PMID: 30800718 PMCID: PMC6383772 DOI: 10.14338/ijpt-17-00030.1] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022] Open
Abstract
Purpose: To investigate the effect of interplay between spot-scanning proton beams and respiration-induced tumor motion on internal target volume coverage for pediatric patients. Materials and Methods: Photon treatments for 10 children with representative tumor motions (1–13 mm superior-inferior) were replanned to simulate single-field uniform dose–optimized proton therapy. Static plans were designed by using average computed tomography (CT) data sets created from 4D CT data to obtain nominal dose distributions. The motion interplay effect was simulated by assigning each spot in the static plan delivery sequence to 1 of 10 respiratory-phase CTs, using the actual patient breathing trace and specifications of a synchrotron-based proton system. Dose distributions for individual phases were deformed onto the space of the average CT and summed to produce the accumulated dose distribution, whose dose-volume histogram was compared with the one from the static plan. Results: Tumor motion had minimal impact on the internal target volume hot spot (D2), which deviated by <3% from the nominal values of the static plans. The cold spot (D98) was also minimally affected, except in 2 patients with diaphragmatic tumor motion exceeding 10 mm. The impact on tumor coverage was more pronounced with respect to the V99 rather than the V95. Decreases of 10% to 49% in the V99 occurred in multiple patients for whom the beam paths traversed the lung-diaphragm interface and were, therefore, more sensitive to respiration-induced changes in the water equivalent path length. Fractionation alone apparently did not mitigate the interplay effect beyond 6 fractions. Conclusion: The interplay effect is not a concern when delivering scanning proton beams to younger pediatric patients with tumors located in the retroperitoneal space and tumor motion of <5 mm. Children and adolescents with diaphragmatic tumor motion exceeding 10 mm require special attention, because significant declines in target coverage and dose homogeneity were seen in simulated treatments of such patients.
Collapse
Affiliation(s)
- Andrew J Boria
- Department of Radiation Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA.,School of Health Sciences, Purdue University, West Lafayette, IN, USA
| | - Jinsoo Uh
- Department of Radiation Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Fakhriddin Pirlepesov
- Department of Radiation Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - James C Stuckey
- Department of Radiation Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA.,Department of Physics, Rhodes College, Memphis, TN, USA
| | - Marian Axente
- Department of Radiation Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Melissa A Gargone
- Department of Radiation Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Chia-Ho Hua
- Department of Radiation Oncology, St. Jude Children's Research Hospital, Memphis, TN, USA
| |
Collapse
|
5
|
Cummings D, Tang S, Ichter W, Wang P, Sturgeon JD, Lee AK, Chang C. Four-dimensional Plan Optimization for the Treatment of Lung Tumors Using Pencil-beam Scanning Proton Radiotherapy. Cureus 2018; 10:e3192. [PMID: 30402360 PMCID: PMC6200439 DOI: 10.7759/cureus.3192] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Purpose This study aimed to evaluate the effectiveness of four-dimensional (4D) robust optimization for proton pencil-beam scanning (PBS) treatment of lung tumors. Patients and methods In seven patients with lung cancer, proton beam therapy was planned using 4D robust optimization over 4D computed tomography (CT) data sets. The gross target volume (GTV) was contoured based on individual breathing phases, and a 5-mm expansion was used to generate the clinical target volume (CTV) for each phase. The 4D optimization was conducted directly on the 4D CT data set. The robust optimization settings included a CT Hounsfield unit (HU) uncertainty of 4% and a setup uncertainty of 5 mm to obtain the CTV. Additional target dose objectives such as those for the internal target volume (ITV) as well as the organ-at-risk (OAR) dose requirements were placed on the average CT. For comparison, three-dimensional (3D) robust optimization was also performed on the average CT. An additional verification 4D CT was performed to verify plan robustness against inter-fractional variations. Results Target coverages were generally higher for 4D optimized plans. The difference was most pronounced for ITV V70Gy when evaluating individual breathing phases. The 4D optimized plans were shown to be able to maintain the ITV coverage at full prescription, while 3D optimized plans could not. More importantly, this difference in ITV V70Gy between the 4D and 3D optimized plans was also consistently observed when evaluating the verification 4D CT, indicating that the 4D optimized plans were more robust against inter-fractional variations. Less difference was seen between the 4D and 3D optimized plans in the lungs criteria: V5Gy and V20Gy. Conclusion The proton PBS treatment plans optimized directly on the 4D CT were shown to be more robust when compared to those optimized on a regular 3D CT. Robust 4D optimization can improve the target coverage for the proton PBS lung treatments.
Collapse
Affiliation(s)
| | - Shikui Tang
- Medical Physics, Texas Center for Proton Therapy, Irving, USA
| | | | - Peng Wang
- Physics, Texas Center for Proton Therapy, Irving, USA
| | | | - Andrew K Lee
- Radiation Oncology, Texas Center for Proton Therapy, Irving, USA
| | - Chang Chang
- Medical Physics, Texas Center for Proton Therapy, Irving, USA
| |
Collapse
|
6
|
Li H, Chang JY. Accounting for, Mitigating, and Choice of Margins for Moving Tumors. Semin Radiat Oncol 2018; 28:194-200. [DOI: 10.1016/j.semradonc.2018.02.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
7
|
Krieger M, Klimpki G, Fattori G, Hrbacek J, Oxley D, Safai S, Weber DC, Lomax AJ, Zhang Y. Experimental validation of a deforming grid 4D dose calculation for PBS proton therapy. ACTA ACUST UNITED AC 2018; 63:055005. [DOI: 10.1088/1361-6560/aaad1e] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
8
|
Birgani MJT, Chegeni N, Zabihzadeh M, Tahmasbi M. Analytical investigation of magnetic field effects on Proton lateral deflection and penetrating depth in the water phantom: A relativistic approach. Electron Physician 2017; 9:5932-5939. [PMID: 29560144 PMCID: PMC5843418 DOI: 10.19082/5932] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2017] [Accepted: 09/12/2017] [Indexed: 11/20/2022] Open
Abstract
Background Integrated proton therapy - MRI systems are capable of delivering high doses to the target tissues near sensitive organs and achieve better therapeutic results; however, the applied magnetic field for imaging, influences the protons path, changes the penetration depth and deflects the particles, laterally, leading to dose distribution variations. Objective To determine the effects of a magnetic field on the range and the lateral deflection of protons, analytically. Methods An analytical survey based on protons energy and range power law relation, without using small angle assumption was done. The penetration depth and lateral deflection of protons with therapeutic energy ranges 60–250 MeV in the presence of uniform magnetic fields of 0–10T intensities, were calculated analytically. Calculations were done for relativistic conditions with Mathematica software version 7.0, and MATLAB 7.0 was applied to plot curves and curve fittings. Results In the presence of a magnetic field, the depth of Bragg peak was decreased and it was shifted laterally. A second order polynomial model with power equation for its coefficients and a power model with quadratic polynomial coefficients predicted the maximum lateral deflection (ymax) and maximum penetration depth (zmax) variations with energy and magnetic field intensity, respectively. Conclusion The applied correction for deflection angle will give more reliable results in initial energy of 250 MeV and 3T magnetic field intensity. For lower energies and magnetic field intensities the differences are negligible, clinically.
Collapse
Affiliation(s)
- Mohammad Javad Tahmasebi Birgani
- Ph.D., Professor, Department of Radiation Therapy, Golestan Hospital, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
| | - Nahid Chegeni
- Ph.D., Assistant Professor, Department of Medical Physics, Faculty of Medicine, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
| | - Mansour Zabihzadeh
- Ph.D., Assistant Professor, Department of Medical Physics, Faculty of Medicine, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
| | - Marziyeh Tahmasbi
- Ph.D. Candidate, Department of Medical Physics, Faculty of Medicine, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
| |
Collapse
|
9
|
Zhou S, Zhu X, Zhang M, Zheng D, Lei Y, Li S, Bennion N, Verma V, Zhen W, Enke C. Estimation of internal organ motion-induced variance in radiation dose in non-gated radiotherapy. Phys Med Biol 2016; 61:8157-8179. [PMID: 27804926 DOI: 10.1088/0031-9155/61/23/8157] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
In the delivery of non-gated radiotherapy (RT), owing to intra-fraction organ motion, a certain degree of RT dose uncertainty is present. Herein, we propose a novel mathematical algorithm to estimate the mean and variance of RT dose that is delivered without gating. These parameters are specific to individual internal organ motion, dependent on individual treatment plans, and relevant to the RT delivery process. This algorithm uses images from a patient's 4D simulation study to model the actual patient internal organ motion during RT delivery. All necessary dose rate calculations are performed in fixed patient internal organ motion states. The analytical and deterministic formulae of mean and variance in dose from non-gated RT were derived directly via statistical averaging of the calculated dose rate over possible random internal organ motion initial phases, and did not require constructing relevant histograms. All results are expressed in dose rate Fourier transform coefficients for computational efficiency. Exact solutions are provided to simplified, yet still clinically relevant, cases. Results from a volumetric-modulated arc therapy (VMAT) patient case are also presented. The results obtained from our mathematical algorithm can aid clinical decisions by providing information regarding both mean and variance of radiation dose to non-gated patients prior to RT delivery.
Collapse
Affiliation(s)
- Sumin Zhou
- Department of Radiation Oncology, University of Nebraska Medical Center, Omaha, NE, USA
| | | | | | | | | | | | | | | | | | | |
Collapse
|
10
|
Dowdell S, Grassberger C, Sharp G, Paganetti H. Fractionated Lung IMPT Treatments. Technol Cancer Res Treat 2016. [DOI: 10.1177/1533034615595761] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Treatment uncertainties in radiotherapy are either systematic or random. This study evaluates the sensitivity of fractionated intensity-modulated proton therapy (IMPT) lung treatments to systematic and random setup uncertainties. Treatments in which single-field homogeneity was restricted to within ±20% (IMPT20%) were compared to full IMPT (IMPTfull) for 10 patients with lung cancer. Four-dimensional Monte Carlo calculations were performed using patient computed tomography geometries with ±5 mm systematic or random setup uncertainties applied over a 35 × 2.5Gy(RBE) treatment course. Fifty fractionated courses were simulated for each patient using both IMPT delivery methods with random setup uncertainties applied each fraction and for 3 energy-dependent spot sizes (big spots, σ≈18-9 mm; intermediate spots, σ≈11-5 mm; and small spots, σ≈4-2 mm). These results were compared to Monte Carlo recalculations of the original treatment plan assuming zero setup uncertainty. Results are presented as the difference in equivalent uniform dose (ΔEUD), V95 (ΔV95), and target dose homogeneity (ΔD1-D99). Over the whole patient cohort, the ΔEUD was 2.0 ± 0.5 (big spots), 1.9 ± 0.7 (intermediate spots), and 1.3 ± 0.4 (small spots) times more sensitive to ±5 mm systematic setup uncertainties in IMPTfull compared to IMPT20%. IMPTfull is 1.9 ± 0.9 (big spots), 2.1 ± 1.1 (intermediate spots), and 1.5 ± 0.6 (small spots) times more sensitive to random setup uncertainties than IMPT20% over a fractionated treatment course. The ΔV95 is at least 1.4 times more sensitive to systematic and random setup uncertainties for IMPTfull for all spot sizes considered. The ΔD1-D99 values coincided within uncertainty limits for both IMPT delivery methods for the 3 spot sizes considered, with higher mean values always observed for IMPTfull. The paired t-test indicated that variations observed between IMPTfull and IMPT20% were significantly different for the majority of scenarios. Significantly larger variations were observed in ΔEUD and ΔV95 in IMPTfull lung treatments in addition to higher mean ΔD1−D99. The steep intra-target dose gradients in IMPTfull make it more susceptible to systematic and random setup uncertainties.
Collapse
Affiliation(s)
- Stephen Dowdell
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA, USA
- Department of Radiation Oncology Medical Physics, Shoalhaven Cancer Care Centre, Illawarra Shoalhaven Cancer & Haematology Network, Nowra, NSW, Australia
| | - Clemens Grassberger
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA, USA
| | - Greg Sharp
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA, USA
| | - Harald Paganetti
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA, USA
| |
Collapse
|
11
|
Dueck J, Knopf AC, Lomax A, Albertini F, Persson GF, Josipovic M, Aznar M, Weber DC, Munck af Rosenschöld P. Robustness of the Voluntary Breath-Hold Approach for the Treatment of Peripheral Lung Tumors Using Hypofractionated Pencil Beam Scanning Proton Therapy. Int J Radiat Oncol Biol Phys 2016; 95:534-541. [DOI: 10.1016/j.ijrobp.2015.11.015] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2015] [Revised: 11/04/2015] [Accepted: 11/09/2015] [Indexed: 12/25/2022]
|
12
|
Rasouli FS, Farhad Masoudi S, Keshazare S, Jette D. Effect of elemental compositions on Monte Carlo dose calculations in proton therapy of eye tumors. Radiat Phys Chem Oxf Engl 1993 2015. [DOI: 10.1016/j.radphyschem.2015.08.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
13
|
Abstract
In this review, image guidance and motion management in radiotherapy for lung cancer is discussed. Motion characteristics of lung tumours and image guidance techniques to obtain motion information are elaborated. Possibilities for management of image guidance and motion in the various steps of the treatment chain are explained, including imaging techniques and beam delivery techniques. Clinical studies using different motion management techniques are reviewed, and finally future directions for image guidance and motion management are outlined.
Collapse
Affiliation(s)
- S S Korreman
- Department of Science, Systems and Models, Roskilde University, Roskilde, Denmark
| |
Collapse
|
14
|
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
|
15
|
Han MC, Yeom YS, Kim CH, Kim S, Sohn JW. New approach based on tetrahedral-mesh geometry for accurate 4D Monte Carlo patient-dose calculation. Phys Med Biol 2015; 60:1601-12. [DOI: 10.1088/0031-9155/60/4/1601] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
16
|
Bert C, Graeff C, Riboldi M, Nill S, Baroni G, Knopf AC. Advances in 4D treatment planning for scanned particle beam therapy - report of dedicated workshops. Technol Cancer Res Treat 2014; 13:485-95. [PMID: 24354749 PMCID: PMC4527425 DOI: 10.7785/tcrtexpress.2013.600274] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2013] [Revised: 09/27/2013] [Accepted: 10/01/2013] [Indexed: 11/25/2022] Open
Abstract
We report on recent progress in the field of mobile tumor treatment with scanned particle beams, as discussed in the latest editions of the 4D treatment planning workshop. The workshop series started in 2009, with about 20 people from 4 research institutes involved, all actively working on particle therapy delivery and development. The first workshop resulted in a summary of recommendations for the treatment of mobile targets, along with a list of requirements to apply these guidelines clinically. The increased interest in the treatment of mobile tumors led to a continuously growing number of attendees: the 2012 edition counted more than 60 participants from 20 institutions and commercial vendors. The focus of research discussions among workshop participants progressively moved from 4D treatment planning to complete 4D treatments, aiming at effective and safe treatment delivery. Current research perspectives on 4D treatments include all critical aspects of time resolved delivery, such as in-room imaging, motion detection, beam application, and quality assurance techniques. This was motivated by the start of first clinical treatments of hepato cellular tumors with a scanned particle beam, relying on gating or abdominal compression for motion mitigation. Up to date research activities emphasize significant efforts in investigating advanced motion mitigation techniques, with a specific interest in the development of dedicated tools for experimental validation. Potential improvements will be made possible in the near future through 4D optimized treatment plans that require upgrades of the currently established therapy control systems for time resolved delivery. But since also these novel optimization techniques rely on the validity of the 4DCT, research focusing on alternative 4D imaging technique, such as MRI based 4DCT generation will continue.
Collapse
Affiliation(s)
- Christoph Bert
- University Clinic Erlangen, Radiation Oncology, Universitatsstrasse 27, 91054 Erlangen, Germany.
| | | | | | | | | | | |
Collapse
|
17
|
Whitaker TJ, Beltran C, Tryggestad E, Bues M, Kruse JJ, Remmes NB, Tasson A, Herman MG. Comparison of two methods for minimizing the effect of delayed charge on the dose delivered with a synchrotron based discrete spot scanning proton beam. Med Phys 2014; 41:081703. [PMID: 25086513 DOI: 10.1118/1.4885961] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Delayed charge is a small amount of charge that is delivered to the patient after the planned irradiation is halted, which may degrade the quality of the treatment by delivering unwarranted dose to the patient. This study compares two methods for minimizing the effect of delayed charge on the dose delivered with a synchrotron based discrete spot scanning proton beam. METHODS The delivery of several treatment plans was simulated by applying a normally distributed value of delayed charge, with a mean of 0.001(SD 0.00025) MU, to each spot. Two correction methods were used to account for the delayed charge. Method one (CM1), which is in active clinical use, accounts for the delayed charge by adjusting the MU of the current spot based on the cumulative MU. Method two (CM2) in addition reduces the planned MU by a predicted value. Every fraction of a treatment was simulated using each method and then recomputed in the treatment planning system. The dose difference between the original plan and the sum of the simulated fractions was evaluated. Both methods were tested in a water phantom with a single beam and simple target geometry. Two separate phantom tests were performed. In one test the dose per fraction was varied from 0.5 to 2 Gy using 25 fractions per plan. In the other test the number fractions were varied from 1 to 25, using 2 Gy per fraction. Three patient plans were used to determine the effect of delayed charge on the delivered dose under realistic clinical conditions. The order of spot delivery using CM1 was investigated by randomly selecting the starting spot for each layer, and by alternating per layer the starting spot from first to last. Only discrete spot scanning was considered in this study. RESULTS Using the phantom setup and varying the dose per fraction, the maximum dose difference for each plan of 25 fractions was 0.37-0.39 Gy and 0.03-0.05 Gy for CM1 and CM2, respectively. While varying the total number of fractions, the maximum dose difference increased at a rate of 0.015 Gy and 0.0018 Gy per fraction for CM1 and CM2, respectively. For CM1, the largest dose difference was found at the location of the first spot in each energy layer, whereas for CM2 the difference in dose was small and showed no dependence on location. For CM1, all of the fields in the patient plans had an area where their excess dose overlapped. No such correlation was found when using CM2. Randomly selecting the starting spot reduces the maximum dose difference from 0.708 to 0.15 Gy. Alternating between first and last spot reduces the maximum dose difference from 0.708 to 0.37 Gy. In the patient plans the excess dose scaled linearly at 0.014 Gy per field per fraction for CM1 and standard delivery order. CONCLUSIONS The predictive model CM2 is superior to a cumulative irradiation model CM1 for minimizing the effects of delayed charge, particularly when considering maximal dose discrepancies and the potential for unplanned hot-spots. This study shows that the dose discrepancy potentially scales at 0.014 Gy per field per fraction for CM1.
Collapse
Affiliation(s)
| | - Chris Beltran
- Mayo Clinic, 200 First Street Southwest, Rochester, Minnesota 55905
| | - Erik Tryggestad
- Mayo Clinic, 200 First Street Southwest, Rochester, Minnesota 55905
| | - Martin Bues
- Mayo Clinic, 5777 East Mayo Boulevard, Phoenix, Arizona 85054
| | - Jon J Kruse
- Mayo Clinic, 200 First Street Southwest, Rochester, Minnesota 55905
| | | | | | - Michael G Herman
- Mayo Clinic, 200 First Street Southwest, Rochester, Minnesota 55905
| |
Collapse
|
18
|
Richter D, Schwarzkopf A, Trautmann J, Krämer M, Durante M, Jäkel O, Bert C. Upgrade and benchmarking of a 4D treatment planning system for scanned ion beam therapy. Med Phys 2013; 40:051722. [PMID: 23635270 DOI: 10.1118/1.4800802] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Upgrade and benchmarking of a research 4D treatment planning system (4DTPS) suitable for realistic patient treatment planning and treatment simulations taking into account specific requirements for scanned ion beam therapy, i.e., modeling of dose heterogeneities due to interplay effects and range changes caused by patient motion and dynamic beam delivery. METHODS The 4DTPS integrates data interfaces to 4D computed tomography (4DCT), deformable image registration and clinically used motion monitoring devices. The authors implemented a novel data model for 4D image segmentation using Boolean mask volume datasets and developed an algorithm propagating a manually contoured reference contour dataset to all 4DCT phases. They further included detailed treatment simulation and dose reconstruction functionality, based on the irregular patient motion and the temporal structure of the beam delivery. The treatment simulation functionality was validated against experimental data from irradiation of moving radiographic films in air, 3D moving ionization chambers in a water phantom, and moving cells in a biological phantom with a scanned carbon ion beam. The performance of the program was compared to results obtained with predecessor programs. RESULTS The measured optical density distributions of the radiographic films were reproduced by the simulations to (-2 ± 12)%. Compared to earlier versions of the 4DTPS, the mean agreement improved by 2%, standard deviations were reduced by 7%. The simulated dose to the moving ionization chambers in water showed an agreement with the measured dose of (-1 ± 4)% for the typical beam configuration. The mean deviation of the simulated from the measured biologically effective dose determined via cell survival was (617 ± 538) mGy relative biological effectiveness corresponding to (10 ± 9)%. CONCLUSIONS The authors developed a research 4DTPS suitable for realistic treatment planning on patient data and capable of simulating dose delivery to a moving patient geometry for scanned ion beams. The accuracy and reliability of treatment simulations improved considerably with respect to earlier versions of the 4DTPS.
Collapse
Affiliation(s)
- D Richter
- GSI Helmholtzzentrum für Schwerionenforschung GmbH, Abt. Biophysik, Planckstraße 1, 64291 Darmstadt, Germany
| | | | | | | | | | | | | |
Collapse
|
19
|
Zhang Y, Knopf A, Tanner C, Boye D, Lomax AJ. Deformable motion reconstruction for scanned proton beam therapy using on-line x-ray imaging. Phys Med Biol 2013; 58:8621-45. [DOI: 10.1088/0031-9155/58/24/8621] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
20
|
Li H, Li Y, Zhang X, Li X, Liu W, Gillin MT, Zhu XR. Dynamically accumulated dose and 4D accumulated dose for moving tumors. Med Phys 2013; 39:7359-67. [PMID: 23231285 DOI: 10.1118/1.4766434] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE The purpose of this work was to investigate the relationship between dynamically accumulated dose (dynamic dose) and 4D accumulated dose (4D dose) for irradiation of moving tumors, and to quantify the dose uncertainty induced by tumor motion. METHODS The authors established that regardless of treatment modality and delivery properties, the dynamic dose will converge to the 4D dose, instead of the 3D static dose, after multiple deliveries. The bounds of dynamic dose, or the maximum estimation error using 4D or static dose, were established for the 4D and static doses, respectively. Numerical simulations were performed (1) to prove the principle that for each phase, after multiple deliveries, the average number of deliveries for any given time converges to the total number of fractions (K) over the number of phases (N); (2) to investigate the dose difference between the 4D and dynamic doses as a function of the number of deliveries for deliveries of a "pulsed beam"; and (3) to investigate the dose difference between 4D dose and dynamic doses as a function of delivery time for deliveries of a "continuous beam." A Poisson model was developed to estimate the mean dose error as a function of number of deliveries or delivered time for both pulsed beam and continuous beam. RESULTS The numerical simulations confirmed that the number of deliveries for each phase converges to K∕N, assuming a random starting phase. Simulations for the pulsed beam and continuous beam also suggested that the dose error is a strong function of the number of deliveries and∕or total deliver time and could be a function of the breathing cycle, depending on the mode of delivery. The Poisson model agrees well with the simulation. CONCLUSIONS Dynamically accumulated dose will converge to the 4D accumulated dose after multiple deliveries, regardless of treatment modality. Bounds of the dynamic dose could be determined using quantities derived from 4D doses, and the mean dose difference between the dynamic dose and 4D dose as a function of number of deliveries and∕or total deliver time was also established.
Collapse
Affiliation(s)
- Heng Li
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA.
| | | | | | | | | | | | | |
Collapse
|
21
|
Aitkenhead AH, Bugg D, Rowbottom CG, Smith E, Mackay RI. Modelling the throughput capacity of a single-accelerator multitreatment room proton therapy centre. Br J Radiol 2013; 85:e1263-72. [PMID: 23175492 DOI: 10.1259/bjr/27428078] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
OBJECTIVE We describe a model for evaluating the throughput capacity of a single-accelerator multitreatment room proton therapy centre with the aims of (1) providing quantitative estimates of the throughput and waiting times and (2) providing insight into the sensitivity of the system to various physical parameters. METHODS A Monte Carlo approach was used to compute various statistics about the modelled centre, including the throughput capacity, fraction times for different groups of patients and beam waiting times. A method of quantifying the saturation level is also demonstrated. RESULTS Benchmarking against the MD Anderson Cancer Center showed good agreement between the modelled (140 ± 4 fractions per day) and reported (133 ± 35 fractions per day) throughputs. A sensitivity analysis of that system studied the impact of beam switch time, the number of treatment rooms, patient set-up times and the potential benefit of having a second accelerator. Finally, scenarios relevant to a potential UK facility were studied, finding that a centre with the same four-room, single-accelerator configuration as the MD Anderson Cancer Center but handling a more complex UK-type caseload would have a throughput reduced by approximately 19%, but still be capable of treating in excess of 100 fractions per 16-h treatment day. CONCLUSIONS The model provides a useful tool to aid in understanding the operating dynamics of a proton therapy facility, and for investigating potential scenarios for prospective centres. ADVANCES IN KNOWLEDGE The model helps to identify which technical specifications should be targeted for future improvements.
Collapse
Affiliation(s)
- A H Aitkenhead
- Christie Medical Physics and Engineering, The Christie NHS Foundation Trust, Manchester, UK.
| | | | | | | | | |
Collapse
|
22
|
Echeverria AE, McCurdy M, Castillo R, Bernard V, Ramos NV, Buckley W, Castillo E, Liu P, Martinez J, Guerrero T. Proton therapy radiation pneumonitis local dose-response in esophagus cancer patients. Radiother Oncol 2012; 106:124-9. [PMID: 23127772 DOI: 10.1016/j.radonc.2012.09.003] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2011] [Revised: 08/03/2012] [Accepted: 09/06/2012] [Indexed: 12/25/2022]
Abstract
PURPOSE This study quantifies pulmonary radiation toxicity in patients who received proton therapy for esophagus cancer. MATERIALS/METHODS We retrospectively studied 100 esophagus cancer patients treated with proton therapy. The linearity of the enhanced FDG uptake vs. proton dose was evaluated using the Akaike Information Criterion (AIC). Pneumonitis symptoms (RP) were assessed using the Common Toxicity Criteria for Adverse Events version 4.0 (CTCAEv4). The interaction of the imaging response with dosimetric parameters and symptoms was evaluated. RESULTS The RP scores were: 0 grade 4/5, 7 grade 3, 20 grade 2, 37 grade 1, and 36 grade 0. Each dosimetric parameter was significantly higher for the symptomatic group. The AIC winning models were 30 linear, 52 linear quadratic, and 18 linear logarithmic. There was no significant difference in the linear coefficient between models. The slope of the FDG vs. proton dose response was 0.022 for the symptomatic and 0.012 for the asymptomatic (p=0.014). Combining dosimetric parameters with the slope did not improve the sensitivity or accuracy in identifying symptomatic cases. CONCLUSIONS The proton radiation dose response on FDG PET/CT imaging exhibited a predominantly linear dose response on modeling. Symptomatic patients had a higher dose response slope.
Collapse
|
23
|
Shin J, Perl J, Schümann J, Paganetti H, Faddegon BA. A modular method to handle multiple time-dependent quantities in Monte Carlo simulations. Phys Med Biol 2012; 57:3295-308. [PMID: 22572201 DOI: 10.1088/0031-9155/57/11/3295] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
A general method for handling time-dependent quantities in Monte Carlo simulations was developed to make such simulations more accessible to the medical community for a wide range of applications in radiotherapy, including fluence and dose calculation. To describe time-dependent changes in the most general way, we developed a grammar of functions that we call 'Time Features'. When a simulation quantity, such as the position of a geometrical object, an angle, a magnetic field, a current, etc, takes its value from a Time Feature, that quantity varies over time. The operation of time-dependent simulation was separated into distinct parts: the Sequence samples time values either sequentially at equal increments or randomly from a uniform distribution (allowing quantities to vary continuously in time), and then each time-dependent quantity is calculated according to its Time Feature. Due to this modular structure, time-dependent simulations, even in the presence of multiple time-dependent quantities, can be efficiently performed in a single simulation with any given time resolution. This approach has been implemented in TOPAS (TOol for PArticle Simulation), designed to make Monte Carlo simulations with Geant4 more accessible to both clinical and research physicists. To demonstrate the method, three clinical situations were simulated: a variable water column used to verify constancy of the Bragg peak of the Crocker Lab eye treatment facility of the University of California, the double-scattering treatment mode of the passive beam scattering system at Massachusetts General Hospital (MGH), where a spinning range modulator wheel accompanied by beam current modulation produces a spread-out Bragg peak, and the scanning mode at MGH, where time-dependent pulse shape, energy distribution and magnetic fields control Bragg peak positions. Results confirm the clinical applicability of the method.
Collapse
Affiliation(s)
- J Shin
- UCSF Helen Diller Family Comprehensive Cancer Center, San Francisco, CA 94143-1708, USA
| | | | | | | | | |
Collapse
|
24
|
Schippers JM. Beam Delivery Systems for Particle Radiation Therapy: Current Status and Recent Developments. ACTA ACUST UNITED AC 2012. [DOI: 10.1142/s1793626809000211] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
An overview is given of different techniques of dose delivery applied in currently operating and planned particle therapy systems. Their advantages and disadvantages will be compared and consequences of the methods for the rest of the instrumentation will be discussed. The interrelationship between beam delivery at the patient and the accelerator system is shown by means of several concrete examples. Apart from a description of several subsystems in a particle therapy facility, design rules for optimizing the reliability of an accelerator and beam delivery system will be discussed, as well as some remarks concerning how to deal with future developments.
Collapse
|
25
|
Abstract
An increasing number of proton therapy facilities are being planned and built at hospital based centers. Most facilities are employing traditional dose delivery methods. A second generation of dose application techniques, based on pencil beam scanning, is slowly being introduced into the commercially available proton therapy systems. New developments in accelerator physics are needed to accommodate and fully exploit these new techniques. At the same time new developments such as the development of small cyclotrons, Dielectric Wall Accelerator (DWA) and laser driven systems, aim for smaller, single room treatment units. In general the benefits of proton therapy could be exploited optimally when achieving a higher level in accuracy, beam energy, beam intensity, safety and system reliability. In this review an overview of the current developments will be given followed by a discussion of upcoming new technologies and needs, like increase of energy, on-line MRI and proton beam splitting for independent uses of treatment rooms.
Collapse
|
26
|
|
27
|
Vadapalli R, Yepes P, Newhauser W, Lichti R. GRID-ENABLED TREATMENT PLANNING FOR PROTON THERAPY USING MONTE CARLO SIMULATIONS. NUCL TECHNOL 2011; 175:16-21. [PMID: 25505349 PMCID: PMC4259269 DOI: 10.13182/nt11-a12263] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Grid computing is an emerging technology that enables computational tasks to be accomplished in a collaborative approach by using a distributed network of computers. The grid approach is especially important for computationally intensive problems that are not tractable with a single computer or even with a small cluster of computers, e.g., radiation transport calculations for cancer therapy. The objective of this work was to extend a Monte Carlo (MC) transport code used for proton radiotherapy to utilize grid computing techniques and demonstrate its promise in reducing runtime from days to minutes. As proof of concept we created the Medical Grid between Texas Tech University and Rice University. Preliminary computational experiments were carried out in the GEANT4 simulation environment for transport of 25 ×106 200 MeV protons in a prostate cancer treatment plan. The simulation speedup was approximately linear; deviations were attributed to the spectrum of parallel runtimes and communication overhead due to Medical Grid computing. The results indicate that ~3 × 105 to 5 × 105 proton events with processor core would result in 65 to 83% efficiency. Extrapolation of our results indicates that about 103 processor cores of the class used here would reduce the MC simulation runtime from 18.3 days to ~1 h.
Collapse
Affiliation(s)
- Ravi Vadapalli
- Texas Tech University, High Performance Computing Center, Box 41167 Lubbock, Texas 79409-1167
| | - Pablo Yepes
- Rice University, Department of Physics and Astronomy, MS 315, 6100 Main Street Houston, Texas 77005
| | - Wayne Newhauser
- The University of Texas MD Anderson Cancer Center, Department of Radiation Physics Unit 1202, 1515 Holcombe Boulevard, Houston, Texas 77030
| | - Roger Lichti
- Texas Tech University, Department of Physics, Box 41051, Lubbock, Texas 79409-1051
| |
Collapse
|
28
|
Munck af Rosenschöld P, Aznar MC, Nygaard DE, Persson GF, Korreman SS, Engelholm SA, Nyström H. A treatment planning study of the potential of geometrical tracking for intensity modulated proton therapy of lung cancer. Acta Oncol 2010; 49:1141-8. [PMID: 20831506 DOI: 10.3109/0284186x.2010.500620] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
BACKGROUND Proton therapy of lung cancer holds the potential for a reduction of the volume of irradiated normal lung tissue. In this work we investigate the robustness of intensity modulated proton therapy (IMPT) plans to motion, and evaluate a geometrical tumour tracking method to compensate for tumour motion. MATERIAL AND METHODS Seven patients with a nine targets with 4DCT scans were selected. IMPT plans were made on the midventilation phase using a 3-field technique. The plans were transferred and calculated on the remaining nine phases of the 4DCT, and the combined dose distribution was summed using deformable image registration (DIR). An additional set of plans were made in which the proton beam was simply geometrically shifted to the centre of the gross tumour volume (GTV), i.e. simulating tracking of the tumour motion but without on-line adjustment of the proton energies. A possible interplay effect between the dynamics of the spot scanning delivery and the tumour motion has not been considered in this work. RESULTS Around 97-100% of the GTV was covered by 95% of the prescribed dose (V95) for a tumour displacement of less than about 1 cm with a static beam. For the remaining three of nine targets with a larger motion the tracking method studied provided a marked improvement over static beam; raising the GTV V95 from 95 to 100%, 82 to 98% and 51 to 97%, respectively. CONCLUSION The possibility of performing DIR and summing the dose on the 4DCT data set was shown to be feasible. The fairly simplistic tracking method suggested here resulted in a marked improvement in GTV coverage for tumours with large intra-fractional motion (>1 cm displacement), indicating that on-line adjustment of the proton energies may be redundant.
Collapse
Affiliation(s)
- Per Munck af Rosenschöld
- Department of Radiation Oncology, Copenhagen University Hospital (Rigshospitalet), Copenhagen, Denmark.
| | | | | | | | | | | | | |
Collapse
|
29
|
Lobo J, Popescu IA. Two new DOSXYZnrc sources for 4D Monte Carlo simulations of continuously variable beam configurations, with applications to RapidArc, VMAT, TomoTherapy and CyberKnife. Phys Med Biol 2010; 55:4431-43. [DOI: 10.1088/0031-9155/55/16/s01] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
30
|
Yepes P, Randeniya S, Taddei PJ, Newhauser WD. A TRACK-REPEATING ALGORITHM FOR FAST MONTE CARLO DOSE CALCULATIONS OF PROTON RADIOTHERAPY. NUCL TECHNOL 2009; 168:736-740. [PMID: 20865140 PMCID: PMC2943387 DOI: 10.13182/nt09-a9298] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Monte Carlo codes are utilized for accurate dose calculations in proton radiation therapy research. While they are superior in accuracy to commonly used analytical dose calculations, they require significantly longer computation times. The aim of this work is to characterize a Monte Carlo track-repeating algorithm to increase computation speed without compromising dosimetric accuracy. The track-repeating approach reduced the CPU time required for a complete dose calculation in voxelized patient anatomy by more than two orders of magnitude, while on average reproducing the results from the traditional Monte Carlo approach within 4% dose difference and within 1-mm distance to agreement.
Collapse
Affiliation(s)
- Pablo Yepes
- Rice University, 6100 Main Street, Houston, Texas 77005
| | | | - Phillip J. Taddei
- University of Texas, M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Unit 94, Houston, Texas 77030
| | - Wayne D. Newhauser
- University of Texas, M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Unit 94, Houston, Texas 77030
| |
Collapse
|
31
|
Fassò A, Ferrari A, Sala PR. Radiation transport calculations and simulations. RADIATION PROTECTION DOSIMETRY 2009; 137:118-133. [PMID: 19778934 DOI: 10.1093/rpd/ncp190] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
This article is an introduction to the Monte Carlo method as used in particle transport. After a description at an elementary level of the mathematical basis of the method, the Boltzmann equation and its physical meaning are presented, followed by Monte Carlo integration and random sampling, and by a general description of the main aspects and components of a typical Monte Carlo particle transport code. In particular, the most common biasing techniques are described, as well as the concepts of estimator and detector. After a discussion of the different types of errors, the issue of quality assurance is briefly considered.
Collapse
Affiliation(s)
- A Fassò
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA.
| | | | | |
Collapse
|
32
|
Proton beam radiotherapy versus three-dimensional conformal stereotactic body radiotherapy in primary peripheral, early-stage non-small-cell lung carcinoma: a comparative dosimetric analysis. Int J Radiat Oncol Biol Phys 2009; 75:950-8. [PMID: 19801106 DOI: 10.1016/j.ijrobp.2009.04.023] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2008] [Revised: 04/05/2009] [Accepted: 04/10/2009] [Indexed: 11/23/2022]
Abstract
PURPOSE Proton radiotherapy (PT) and stereotactic body radiotherapy (SBRT) have the capacity to optimize the therapeutic ratio. We analyzed the dosimetric differences between PT and SBRT in treating primary peripheral early-stage non-small-cell lung cancer. METHODS AND MATERIALS Eight patients were simulated, planned, and treated with SBRT according to accepted techniques. SBRT treatments were retrospectively planned using heterogeneity corrections. PT treatment plans were generated using single-, two-, and three-field passively scattered and actively scanned proton beams. Calculated dose characteristics were compared. RESULTS Comparable planning target volume (PTV) median minimum and maximum doses were observed between PT and SBRT plans. Higher median maximum doses 2 cm from the PTV were observed for PT, but higher median PTV doses were observed for SBRT. The total lung mean and V5 doses were significantly lower with actively scanned PT. The lung V13 and V20 were comparable. The dose to normal tissues was lower with PT except to skin and ribs. Although the maximum doses to skin and ribs were similar or higher with PT, the median doses to these structures were higher with SBRT. Passively scattered plans, compared with actively scanned plans, typically demonstrated higher doses to the PTV, lung, and organs at risk. CONCLUSIONS Single-, two-, or three-field passively or actively scanned proton therapy delivered comparable PTV dose with generally less dose to normal tissues in these hypothetic treatments. Actively scanned beam plans typically had more favorable dose characteristics to the target, lung, and other soft tissues compared with the passively scanned plans. The clinical significance of these findings remains to be determined.
Collapse
|
33
|
Peterson SW, Polf J, Bues M, Ciangaru G, Archambault L, Beddar S, Smith A. Experimental validation of a Monte Carlo proton therapy nozzle model incorporating magnetically steered protons. Phys Med Biol 2009; 54:3217-29. [DOI: 10.1088/0031-9155/54/10/017] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
34
|
Riboldi M, Chen GTY, Baroni G, Paganetti H, Seco J. Design and testing of a simulation framework for dosimetric motion studies integrating an anthropomorphic computational phantom into four-dimensional Monte Carlo. Technol Cancer Res Treat 2009; 7:449-56. [PMID: 19044324 DOI: 10.1177/153303460800700606] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We have designed a simulation framework for motion studies in radiation therapy by integrating the anthropomorphic NCAT phantom into a 4D Monte Carlo dose calculation engine based on DPM. Representing an artifact-free environment, the system can be used to identify class solutions as a function of geometric and dosimetric parameters. A pilot dynamic conformal study for three lesions ( approximately 2.0 cm) in the right lung was performed (70 Gy prescription dose). Tumor motion changed as a function of tumor location, according to the anthropomorphic deformable motion model. Conformal plans were simulated with 0 to 2 cm margin for the aperture, with additional 0.5 cm for beam penumbra. The dosimetric effects of intensity modulated radiotherapy (IMRT) vs. conformal treatments were compared in a static case. Results show that the Monte Carlo simulation framework can model tumor tracking in deformable anatomy with high accuracy, providing absolute doses for IMRT and conformal radiation therapy. A target underdosage of up to 3.67 Gy (lower lung) was highlighted in the composite dose distribution mapped at exhale. Such effects depend on tumor location and treatment margin and are affected by lung deformation and ribcage motion. In summary, the complexity in the irradiation of moving targets has been reduced to a controlled simulation environment, where several treatment options can be accurately modeled and quantified The implemented tools will be utilized for extensive motion study in lung/liver irradiation.
Collapse
Affiliation(s)
- M Riboldi
- Department of Bioengineering, Politecnico di Milano University, P.za Leonardo da Vinci 32, 20133 Milano, Italy.
| | | | | | | | | |
Collapse
|
35
|
|
36
|
Paganetti H, Jiang H, Parodi K, Slopsema R, Engelsman M. Clinical implementation of full Monte Carlo dose calculation in proton beam therapy. Phys Med Biol 2008; 53:4825-53. [PMID: 18701772 DOI: 10.1088/0031-9155/53/17/023] [Citation(s) in RCA: 174] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
37
|
Hollmark M, Gudowska I, Belkić D, Brahme A, Sobolevsky N. An analytical model for light ion pencil beam dose distributions: multiple scattering of primary and secondary ions. Phys Med Biol 2008; 53:3477-91. [DOI: 10.1088/0031-9155/53/13/005] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
38
|
Georg D, Hillbrand M, Stock M, Dieckmann K, Pötter R. Can protons improve SBRT for lung lesions? Dosimetric considerations. Radiother Oncol 2008; 88:368-75. [PMID: 18405986 DOI: 10.1016/j.radonc.2008.03.007] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2007] [Revised: 03/05/2008] [Accepted: 03/06/2008] [Indexed: 10/22/2022]
Abstract
BACKGROUND AND PURPOSE The aim of the present study was to investigate potential dosimetric benefits of proton therapy for hypofractionated stereotactic body radiotherapy (SBRT). MATERIALS AND METHOD Twelve patients undergoing hypofractionated SBRT at the Medical University Vienna were selected. Passively scattered protons (PT) and intensity modulated proton therapy (IMPT) were evaluated against a conformal photon technique (3D-CRT), assuming a fractionation of 3x15Gy, prescribed to the 65% isodose. For all treatment techniques shallow breathing with abdominal compression (SB+AC) was compared with a deep inspiration breath hold technique (DIBH). Treatment planning was done with XiO (CMS, USA). Target conformity, dose-volume histograms (DVH) and various associated dosimetric parameters were considered for the planning target volume (PTV), lung, heart and esophagus. RESULTS For both breathing conditions conformity indices were very similar. They were between 0.75 and 0.78 for IMPT and 3D-CRT and around 0.55 for PT using 2-3 beams. Irrespective of treatment modality, DVHs for the ipsilateral lung were improved with the DIBH technique. For the PT technique, the 2Gy isodose (V2Gy) covered on average 7-9% less lung volume compared to 3D-CRT, for IMPT this reduction was more than 10%. Volumes covered the 4 and 6Gy isodoses were 2-4% smaller for IMPT, but very similar for PT and 3D-CRT. Both proton techniques achieved full sparing of the contralateral lung and superior sparing of the heart. Maximum doses to the heart and esophagus were on average around 3Gy for 3D-CRT and almost 0Gy for both proton techniques. For 3D-CRT average V2Gy values for the heart could be reduced from 64% in shallow breathing to 34% in DIBH. V2Gy for protons was negligible. CONCLUSIONS Only small dosimetric differences were found between photons and protons for SBRT of lung lesions. Whether these small dosimetric benefits translate in reduced side effects or have the potential to improve local control rates remains to be demonstrated in clinical studies.
Collapse
Affiliation(s)
- Dietmar Georg
- Department of Radiotherapy, Medical University Vienna/AKH Vienna, Vienna, Austria.
| | | | | | | | | |
Collapse
|
39
|
Widesott L, Amichetti M, Schwarz M. Proton therapy in lung cancer: clinical outcomes and technical issues. A systematic review. Radiother Oncol 2008; 86:154-64. [PMID: 18241945 DOI: 10.1016/j.radonc.2008.01.003] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2007] [Revised: 01/03/2008] [Accepted: 01/03/2008] [Indexed: 12/25/2022]
Abstract
BACKGROUND AND PURPOSE To determine whether, according to the currently available literature, proton therapy (PT) has a role in the treatment of non-small-cell lung cancer (NSCLC), to assess its safety and efficacy and to evaluate the main technical issues specifically related to this treatment technique. MATERIALS AND METHODS During March 2007, two independent researchers conducted a systematic review of the current data on the treatment of NSCLC with PT. RESULTS In total, 113 reports were retrieved, 17 of which were included in the analysis. There were no prospective trials (randomized or non-randomized). Nine uncontrolled single-arm studies were available from three PT centers, providing clinical outcomes for 214 patients in total. These reports were mainly related to stage I-II tumors, with results comparable to those obtained with surgery, without significant toxicity. In addition, two papers were found that compared photon and proton dose distributions, which showed a potential for dose escalation and/or a sparing of the organ at risk with PT. Finally, six studies analyzed dosimetric and technical issues related with PT, mainly underlining the difficulties in designing dose distributions that are representative of the dose actually delivered during treatment. CONCLUSIONS Although from a physical point of view PT is a good option for the treatment of NSCLC, limited data are available on its application in the clinical practice. Furthermore, the application of PT to lung cancer does present technical challenges. Because of the small number of institutions involved in the treatment of this disease, number of patients, and methodological weaknesses of the trials it is therefore not possible to draw definitive conclusions about the superiority of PT with respect to the photon techniques currently available for the treatment of NSCLC.
Collapse
|
40
|
Zhang J, Xu GX, Shi C, Fuss M. Development of a geometry-based respiratory motion-simulating patient model for radiation treatment dosimetry. J Appl Clin Med Phys 2008. [PMID: 18449164 PMCID: PMC2737526 DOI: 10.1120/jacmp.v9i1.2700] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
Temporal and spatial anatomic changes caused by respiration during radiation treatment delivery can lead to discrepancies between prescribed and actual radiation doses. The present paper documents a study to construct a respiratory‐motion‐simulating, four‐dimensional (4D) anatomic and dosimetry model for the study of the dosimetric effects of organ motion for various radiation treatment plans and delivery strategies. The non‐uniform rational B‐splines (NURBS) method has already been used to reconstruct a three‐dimensional (3D) VIP‐Man (“visible photographic man”) model that can reflect the deformation of organs during respiration by using time‐dependent equations to manipulate surface control points. The EGS4 (Electron Gamma Shower, version 4) Monte Carlo code is then used to apply the 4D model to dose simulation. We simulated two radiation therapy delivery scenarios: gating treatment and 4D image‐guided treatment. For each delivery scenario, we developed one conformal plan and one intensity‐modulated radiation therapy plan. A lesion in the left lung was modeled to investigate the effect of respiratory motion on radiation dose distributions. Based on target dose–volume histograms, the importance of using accurate gating to improve the dose distribution is demonstrated. The results also suggest that, during 4D image‐guided treatment delivery, monitoring of the patient's breathing pattern is critical. This study demonstrates the potential of using a “standard” motion‐simulating patient model for 4D treatment planning and motion management. PACS numbers: 87.53.Bn, 87.53.Kn, 87.53.Tf, 87.53.Wz, 87.57.Gg, 89.80.+h
Collapse
Affiliation(s)
- Juying Zhang
- Nuclear Engineering and Engineering Physics, Rensselaer Polytechnic Institute, Troy, New York, USA
| | | | | | | |
Collapse
|
41
|
Xing L, Siebers J, Keall P. Computational Challenges for Image-Guided Radiation Therapy: Framework and Current Research. Semin Radiat Oncol 2007; 17:245-57. [PMID: 17903702 DOI: 10.1016/j.semradonc.2007.07.004] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
It is arguable that the imaging and delivery hardware necessary for delivering real-time adaptive image-guided radiotherapy is available on high-end linear accelerators. Robust and computationally efficient software is the limiting factor in achieving highly accurate and precise radiotherapy to the constantly changing anatomy of a cancer patient. The limitations are not caused by the availability of algorithms but rather issues of reliability, integration, and calculation time. However, each of the software components is an active area of research and development at academic and commercial centers. This article describes the software solutions in 4 broad areas: deformable image registration, adaptive replanning, real-time image guidance, and dose calculation and accumulation. Given the pace of technological advancement, the integration of these software solutions to develop real-time adaptive image-guided radiotherapy and the associated challenges they bring will be implemented to varying degrees by all major manufacturers over the coming years.
Collapse
Affiliation(s)
- Lei Xing
- Department of Radiation Oncology, Stanford University, Stanford, CA 94305-5304, USA
| | | | | |
Collapse
|
42
|
|
43
|
Hérault J, Iborra N, Serrano B, Chauvel P. Spread-out Bragg peak and monitor units calculation with the Monte Carlo code MCNPX. Med Phys 2007; 34:680-8. [PMID: 17388186 DOI: 10.1118/1.2431473] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
The aim of this work was to study the dosimetric potential of the Monte Carlo code MCNPX applied to the protontherapy field. For series of clinical configurations a comparison between simulated and experimental data was carried out, using the proton beam line of the MEDICYC isochronous cyclotron installed in the Centre Antoine Lacassagne in Nice. The dosimetric quantities tested were depth-dose distributions, output factors, and monitor units. For each parameter, the simulation reproduced accurately the experiment, which attests the quality of the choices made both in the geometrical description and in the physics parameters for beam definition. These encouraging results enable us today to consider a simplification of quality control measurements in the future. Monitor Units calculation is planned to be carried out with preestablished Monte Carlo simulation data. The measurement, which was until now our main patient dose calibration system, will be progressively replaced by computation based on the MCNPX code. This determination of Monitor Units will be controlled by an independent semi-empirical calculation.
Collapse
Affiliation(s)
- J Hérault
- Centre Antoine Lacassagne, Cyclotron Biomedical, 227 avenue de la Lanterne, 06200 Nice, France.
| | | | | | | |
Collapse
|
44
|
|