1
|
Yanheng P, Zhen W, Yisheng H, Shenshen G, Rui Q, Hui Z, Junli L. CPU-GPU-coupled acceleration method for point flux calculation in Monte Carlo particle transport. RADIATION PROTECTION DOSIMETRY 2024; 200:525-537. [PMID: 38411255 DOI: 10.1093/rpd/ncae032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Revised: 11/22/2023] [Accepted: 02/10/2024] [Indexed: 02/28/2024]
Abstract
In Monte Carlo particle transport simulations, point flux tallying is a variance reduction technique that performs well with small detectors and finds broad application in source-detector problems and local point dose calculations. However, its use in large-scale point flux tallying computation adds substantial computational time. To address this issue, we propose a CPU-GPU-coupled acceleration method, which separates the complex logic and computationally intensive parts of particle transport calculation and assigns them to the CPU and GPU, respectively. This proposed method greatly enhances the efficiency of large-scale point flux tallies, providing significant convenience for subsequent dose calculations and other related steps. We validated our method by comparing the performance of a pure CPU program with our CPU-GPU accelerated program using the NUREG/CR-6115 PWR benchmark problem. The results indicate identical outcomes for photon point flux estimation, with the accelerated program being ~50 times faster.
Collapse
Affiliation(s)
- Pu Yanheng
- Department of Engineering Physics, Tsinghua University, Beijing 100084, China
- Key Laboratory of Particle & Radiation Imaging, Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Wu Zhen
- Department of Engineering Physics, Tsinghua University, Beijing 100084, China
- Key Laboratory of Particle & Radiation Imaging, Ministry of Education, Tsinghua University, Beijing 100084, China
- Nuctech Company Limited, Beijing 100084, China
| | - Hao Yisheng
- Department of Engineering Physics, Tsinghua University, Beijing 100084, China
- Key Laboratory of Particle & Radiation Imaging, Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Gao Shenshen
- Department of Engineering Physics, Tsinghua University, Beijing 100084, China
- Key Laboratory of Particle & Radiation Imaging, Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Qiu Rui
- Department of Engineering Physics, Tsinghua University, Beijing 100084, China
- Key Laboratory of Particle & Radiation Imaging, Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Zhang Hui
- Department of Engineering Physics, Tsinghua University, Beijing 100084, China
- Key Laboratory of Particle & Radiation Imaging, Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Li Junli
- Department of Engineering Physics, Tsinghua University, Beijing 100084, China
- Key Laboratory of Particle & Radiation Imaging, Ministry of Education, Tsinghua University, Beijing 100084, China
| |
Collapse
|
2
|
Kollitz E, Roew M, Han H, Pinto M, Kamp F, Kim CH, Schwarz M, Belka C, Newhauser W, Parodi K, Dedes G. Applications of a patient-specific whole-body CT-mesh hybrid computational phantom in second cancer risk prediction. Phys Med Biol 2022; 67. [DOI: 10.1088/1361-6560/ac8851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 08/09/2022] [Indexed: 11/11/2022]
Abstract
Abstract
Objective. CT-mesh hybrid phantoms (or ‘hybrid(s)’) made from integrated patient CT data and mesh-type reference computational phantoms (MRCPs) can be beneficial for patient-specific whole-body dose evaluation, but this benefit has yet to be evaluated for second cancer risk prediction. The purpose of this study is to compare the hybrid’s ability to predict risk throughout the body with a patient-scaled MRCP against ground truth whole-body CTs (WBCTs). Approach. Head and neck active scanning proton treatment plans were created for and simulated on seven hybrids and the corresponding scaled MRCPs and WBCTs. Equivalent dose throughout the body was calculated and input into five second cancer risk models for both excess absolute and excess relative risk (EAR and ERR). The hybrid phantom was evaluated by comparing equivalent dose and risk predictions against the WBCT. Main results. The hybrid most frequently provides whole-body second cancer risk predictions which are closer to the ground truth when compared to a scaled MRCP alone. The performance of the hybrid relative to the scaled MRCP was consistent across ERR, EAR, and all risk models. For all in-field organs, where the hybrid shares the WBCT anatomy, the hybrid was better than or equal to the scaled MRCP for both equivalent dose and risk prediction. For out-of-field organs across all patients, the hybrid’s equivalent dose prediction was superior than the scaled MRCP in 48% of all comparisons, equivalent for 34%, and inferior for 18%. For risk assessment in the same organs, the hybrid’s prediction was superior than the scaled MRCP in 51.8% of all comparisons, equivalent in 28.6%, and inferior in 19.6%. Significance. Whole-body risk predictions from the CT-mesh hybrid have shown to be more accurate than those from a reference phantom alone. These hybrids could aid in risk-optimized treatment planning and individual risk assessment to minimize second cancer incidence.
Collapse
|
3
|
Gordon KB, Saburov VO, Koryakin SN, Gulidov IA, Fatkhudinov TK, Arutyunyan IV, Kaprin AD, Solov’ev AN. Calculation of the Biological Efficiency of the Proton Component from 14.8 MeV Neutron Irradiation in Computational Biology with Help of Video Cards. Bull Exp Biol Med 2022; 173:281-285. [DOI: 10.1007/s10517-022-05534-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Indexed: 10/17/2022]
|
4
|
De Simoni M, Battistoni G, De Gregorio A, De Maria P, Fischetti M, Franciosini G, Marafini M, Patera V, Sarti A, Toppi M, Traini G, Trigilio A, Schiavi A. A Data-Driven Fragmentation Model for Carbon Therapy GPU-Accelerated Monte-Carlo Dose Recalculation. Front Oncol 2022; 12:780784. [PMID: 35402249 PMCID: PMC8990885 DOI: 10.3389/fonc.2022.780784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Accepted: 02/04/2022] [Indexed: 11/13/2022] Open
Abstract
The advent of Graphics Processing Units (GPU) has prompted the development of Monte Carlo (MC) algorithms that can significantly reduce the simulation time with respect to standard MC algorithms based on Central Processing Unit (CPU) hardware. The possibility to evaluate a complete treatment plan within minutes, instead of hours, paves the way for many clinical applications where the time-factor is important. FRED (Fast paRticle thErapy Dose evaluator) is a software that exploits the GPU power to recalculate and optimise ion beam treatment plans. The main goal when developing the FRED physics model was to balance accuracy, calculation time and GPU execution guidelines. Nowadays, FRED is already used as a quality assurance tool in Maastricht and Krakow proton clinical centers and as a research tool in several clinical and research centers across Europe. Lately the core software has been updated including a model of carbon ions interactions with matter. The implementation is phenomenological and based on carbon fragmentation data currently available. The model has been tested against the MC FLUKA software, commonly used in particle therapy, and a good agreement was found. In this paper, the new FRED data-driven model for carbon ion fragmentation will be presented together with the validation tests against the FLUKA MC software. The results will be discussed in the context of FRED clinical applications to 12C ions treatment planning.
Collapse
Affiliation(s)
- Micol De Simoni
- Department of Physics, University of Rome "Sapienza", Rome, Italy.,INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy
| | | | - Angelica De Gregorio
- Department of Physics, University of Rome "Sapienza", Rome, Italy.,INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy
| | - Patrizia De Maria
- Department of Medico-Surgical Sciences and Biotechnologies, Post-Graduate School in Medical Physics, Rome, Italy
| | - Marta Fischetti
- INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy.,Department of Scienze di Base e Applicate per l'Ingegneria (SBAI), University of Rome "Sapienza", Rome, Italy
| | - Gaia Franciosini
- Department of Physics, University of Rome "Sapienza", Rome, Italy.,INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy
| | - Michela Marafini
- INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy.,Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi, Rome, Italy
| | - Vincenzo Patera
- INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy.,Department of Scienze di Base e Applicate per l'Ingegneria (SBAI), University of Rome "Sapienza", Rome, Italy
| | - Alessio Sarti
- INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy.,Department of Scienze di Base e Applicate per l'Ingegneria (SBAI), University of Rome "Sapienza", Rome, Italy
| | - Marco Toppi
- Department of Scienze di Base e Applicate per l'Ingegneria (SBAI), University of Rome "Sapienza", Rome, Italy.,INFN Laboratori Nazionali di Frascati, Frascati, Italy
| | - Giacomo Traini
- INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy
| | - Antonio Trigilio
- Department of Physics, University of Rome "Sapienza", Rome, Italy.,INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy
| | - Angelo Schiavi
- INFN (Istituto Nazionale di Fisica Nucleare) section of Roma 1, Rome, Italy.,Department of Scienze di Base e Applicate per l'Ingegneria (SBAI), University of Rome "Sapienza", Rome, Italy
| |
Collapse
|
5
|
Ramesh P, Liu H, Gu W, Sheng K. Fixed Beamline Optimization for Intensity Modulated Carbon-Ion Therapy. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2022; 6:288-293. [PMID: 36092271 PMCID: PMC9457306 DOI: 10.1109/trpms.2021.3092296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
A major obstacle for the adoption of heavy ion therapy is the cost and technical difficulties to construct and maintain a rotational gantry. Many heavy ion treatment facilities instead choose to construct fixed beamlines as a compromise, which we propose to mitigate with optimized treatment couch angle. We formulate the integrated beam orientation and scanning spot optimization problem as a quadratic cost function with a group sparsity regularization term. The optimization problem is efficiently solved using fast iterative shrinkage-thresholding algorithm (FISTA). To test the method, we created the fixed beamline plans with couch rotation (FBCR) and without couch rotation (FB) for intensity modulated carbon-ion therapy (IMCT) and compared with the ideal scenario where both the couch and gantry have 360 degrees of freedom (GCR). FB, FBCR, and GCR IMCT plans were compared for ten pancreas cases. The FBCR plans show comparable PTV coverage and OAR doses for each pancreas case. In conclusion, the dosimetric limitation of fixed beams in heavy ion radiotherapy may be largely mitigated with integrated beam orientation optimization of the couch rotation.
Collapse
Affiliation(s)
- Pavitra Ramesh
- Physics and Biology in Medicine interdepartmental program, University of California Los Angeles, Los Angeles, CA 90025 USA
| | - Hengjie Liu
- Physics and Biology in Medicine interdepartmental program, University of California Los Angeles, Los Angeles, CA 90025 USA
| | - Wenbo Gu
- Department of Radiation Oncology, University of Pennsylvania, Philadelphia, PA 19104 USA
| | - Ke Sheng
- Physics and Biology in Medicine interdepartmental program, University of California Los Angeles, Los Angeles, CA 90025 USA
| |
Collapse
|
6
|
Wu C, Nguyen D, Xing Y, Montero AB, Schuemann J, Shang H, Pu Y, Jiang S. Improving Proton Dose Calculation Accuracy by Using Deep Learning. MACHINE LEARNING: SCIENCE AND TECHNOLOGY 2021; 2:015017. [PMID: 35965743 PMCID: PMC9374098 DOI: 10.1088/2632-2153/abb6d5] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2020] [Revised: 08/18/2020] [Accepted: 09/09/2020] [Indexed: 12/28/2022] Open
Abstract
Introduction Pencil beam (PB) dose calculation is fast but inaccurate due to the approximations when dealing with inhomogeneities. Monte Carlo (MC) dose calculation is the most accurate method but it is time consuming. The aim of this study was to develop a deep learning model that can boost the accuracy of PB dose calculation to the level of MC dose by converting PB dose to MC dose for different tumor sites. Methods The proposed model uses the PB dose and CT image as inputs to generate the MC dose. We used 290 patients (90 head and neck, 93 liver, 75 prostate and 32 lung) to train, validate, and test the model. For each tumor site, we performed four numerical experiments to explore various combinations of training datasets. Results Training the model on data from all tumor sites together and using the dose distribution of each individual beam as input yielded the best performance for all four tumor sites. The average gamma passing rate (1mm/1%) between the converted and the MC dose was 92.8%, 92.7%, 89.7% and 99.6% for head and neck, liver, lung, and prostate test patients, respectively. The average dose conversion time for a single field was less than 4 seconds. The trained model can be adapted to new datasets through transfer learning. Conclusions Our deep learning-based approach can quickly boost the accuracy of PB dose to that of MC dose. The developed model can be added to the clinical workflow of proton treatment planning to improve dose calculation accuracy.
Collapse
Affiliation(s)
- Chao Wu
- Medical Artificial Intelligence and Automation (MAIA) Laboratory, Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, United States of America
- Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai, People’s Republic of China
- University of Chinese Academy of Sciences, Beijing, People’s Republic of China
| | - Dan Nguyen
- Medical Artificial Intelligence and Automation (MAIA) Laboratory, Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, United States of America
| | - Yixun Xing
- Medical Artificial Intelligence and Automation (MAIA) Laboratory, Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, United States of America
| | - Ana Barragan Montero
- University of Chinese Academy of Sciences, Beijing, People’s Republic of China
- Molecular Imaging Radiation Oncology (MIRO) Laboratory, UCLouvain, Brussels, Belgium
| | - Jan Schuemann
- Department of Radiation Oncology, Massachusetts General Hospital, Boston, MA, United States of America
| | - Haijiao Shang
- Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai, People’s Republic of China
- University of Chinese Academy of Sciences, Beijing, People’s Republic of China
| | - Yuehu Pu
- Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai, People’s Republic of China
| | - Steve Jiang
- Medical Artificial Intelligence and Automation (MAIA) Laboratory, Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, United States of America
| |
Collapse
|
7
|
Ma J, Wan Chan Tseung HS, Courneyea L, Beltran C, Herman MG, Remmes NB. Robust radiobiological optimization of ion beam therapy utilizing Monte Carlo and microdosimetric kinetic model. ACTA ACUST UNITED AC 2020; 65:155020. [DOI: 10.1088/1361-6560/aba08b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
8
|
Shrestha D, Tsai MY, Qin N, Zhang Y, Jia X, Wang J. Dosimetric evaluation of 4D-CBCT reconstructed by Simultaneous Motion Estimation and Image Reconstruction (SMEIR) for carbon ion therapy of lung cancer. Med Phys 2019; 46:4087-4094. [PMID: 31299097 DOI: 10.1002/mp.13706] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2019] [Revised: 06/06/2019] [Accepted: 06/06/2019] [Indexed: 01/03/2023] Open
Abstract
PURPOSE Motion management is critical for the efficacy of carbon ion therapy for moving targets such as lung tumors. We evaluated the feasibility of using four-dimensional cone beam computed tomography (4D-CBCT) reconstructed by Simultaneous Motion Estimation and Image Reconstruction (SMEIR) for dose calculation and accumulation in carbon ion treatment of lung cancer. METHODS Motion-compensated 4D-CBCT images were reconstructed with the SMEIR algorithm to capture the most updated anatomy and motion with an updated interphase motion model on the treatment day. Projections of all CBCT phases were simulated from the planning 4D-CT by the ray tracing technique. Treatment planning and dose calculation were performed with a GPU-based Monte Carlo dose calculation software for carbon ion therapy. The treatment plan was optimized on the average computed tomography (CT) to obtain optimal intensity of the carbon ions. From the optimized plan, dose distributions on individual phases of 4D-CT and 4D-CBCT were calculated by the Monte Carlo-based dose engine. Dose accumulation was performed on 4D-CBCT images using deformable vector fields (DVF) generated by SMEIR. The accumulated planning target volume (PTV) dose based on 4D-CBCT was then compared to the accumulated dose calculated on 4D-CT, where the DVFs between different phases were obtained by the demons deformable registration algorithm. RESULTS Dose value histograms (DVH) as well as absolute deviations of the maximum dose ( Δ D max ), mean dose ( Δ D mean ), and dose coverage metrics ( Δ V 95 % and Δ V 100 % ) for PTV were quantitatively evaluated for the two sets of plans. Good agreement was found between the 4D-CT and 4D-CBCT-based PTV-DVH curves. The average values of Δ D max , Δ D mean , Δ V 95 % , and Δ V 100 % calculated between the 4D-CT and SMEIR-4D-CBCT-based plans were 1.91 % , 3.55 % , 2.12%, and 1.15 % , respectively, for the PTVs of ten patient case studies. CONCLUSIONS Based on these results, SMEIR-reconstructed 4D-CBCTs can potentially be used for motion estimation, dose evaluation, and adaptive treatment planning in lung cancer carbon ion therapy.
Collapse
Affiliation(s)
- Deepak Shrestha
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, 75390, USA
| | - Min-Yu Tsai
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, 75390, USA
| | - Nan Qin
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, 75390, USA
| | - You Zhang
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, 75390, USA
| | - Xun Jia
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, 75390, USA
| | - Jing Wang
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, 75390, USA
| |
Collapse
|
9
|
Durante M, Paganetti H, Pompos A, Kry SF, Wu X, Grosshans DR. Report of a National Cancer Institute special panel: Characterization of the physical parameters of particle beams for biological research. Med Phys 2018; 46:e37-e52. [PMID: 30506898 DOI: 10.1002/mp.13324] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2017] [Revised: 10/28/2018] [Accepted: 11/05/2018] [Indexed: 12/16/2022] Open
Abstract
PURPOSE To define the physical parameters needed to characterize a particle beam in order to allow intercomparison of different experiments performed using different ions at the same facility and using the same ion at different facilities. METHODS At the request of the National Cancer Institute (NCI), a special panel was convened to review the current status of the field and to provide suggested metrics for reporting the physical parameters of particle beams to be used for biological research. A set of physical parameters and measurements that should be performed by facilities and understood and reported by researchers supported by NCI to perform pre-clinical radiobiology and medical physics of heavy ions were generated. RESULTS Standard measures such as radiation delivery technique, beam modifiers used, nominal energy, field size, physical dose and dose rate should all be reported. However, more advanced physical measurements, including detailed characterization of beam quality by microdosimetric spectrum and fragmentation spectra, should also be established and reported. Details regarding how such data should be incorporated into Monte Carlo simulations and the proper reporting of simulation details are also discussed. CONCLUSIONS In order to allow for a clear relation of physical parameters to biological effects, facilities and researchers should establish and report detailed physical characteristics of the irradiation beams utilized including both standard and advanced measures. Biological researchers are encouraged to actively engage facility staff and physicists in the design and conduct of experiments. Modeling individual experimental setups will allow for the reporting of the uncertainties in the measurement or calculation of physical parameters which should be routinely reported.
Collapse
Affiliation(s)
- Marco Durante
- Biophysics Department, GSI Helmholtzzentrum für Schwerionenforschung and Technische Universität Darmstadt, Institute of Condensed Matter Physics, Planckstraße 1, 64291, Darmstadt, Germany
| | - Harald Paganetti
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA, 02114, USA
| | - Arnold Pompos
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX, 75390, USA
| | - Stephen F Kry
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA
| | - Xiaodong Wu
- Department of Medical Physics, Shanghai Proton and Heavy Ion Center, Shanghai, China
| | - David R Grosshans
- Departments of Radiation and Experimental Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, 77054, USA
| |
Collapse
|
10
|
Wedenberg M, Beltran C, Mairani A, Alber M. Advanced Treatment Planning. Med Phys 2018; 45:e1011-e1023. [PMID: 30421811 DOI: 10.1002/mp.12943] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2017] [Revised: 03/22/2018] [Accepted: 04/22/2018] [Indexed: 12/15/2022] Open
Abstract
Treatment planning for protons and heavier ions is adapting technologies originally developed for photon dose optimization, but also has to meet its particular challenges. Since the quality of the applied dose is more sensitive to geometric uncertainties, treatment plan robust optimization has a much more prominent role in particle therapy. This has led to specific planning tools, approaches, and research into new formulations of the robust optimization problems. Tools for solution space navigation and automatic planning are also being adapted to particle therapy. These challenges become even greater when detailed models of relative biological effectiveness (RBE) are included into dose optimization, as is required for heavier ions.
Collapse
Affiliation(s)
| | - Chris Beltran
- Division of Medical Physics, Department of Radiation Oncology, Mayo Clinic, Rochester, MN, USA
| | - Andrea Mairani
- Heidelberg Ion Therapy Center (HIT), Heidelberg University Hospital, Heidelberg, Germany.,National Center for Radiation Research in Oncology (NCRO), Heidelberg Institute for Radiation Oncology (HIRO), Heidelberg, Germany.,The National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy
| | - Markus Alber
- The National Centre for Oncological Hadrontherapy (CNAO), Pavia, Italy.,Section for Medical Physics, Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
| |
Collapse
|
11
|
FRoG-A New Calculation Engine for Clinical Investigations with Proton and Carbon Ion Beams at CNAO. Cancers (Basel) 2018; 10:cancers10110395. [PMID: 30360576 PMCID: PMC6266031 DOI: 10.3390/cancers10110395] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2018] [Revised: 10/15/2018] [Accepted: 10/16/2018] [Indexed: 11/17/2022] Open
Abstract
A fast and accurate dose calculation engine for hadrontherapy is critical for both routine clinical and advanced research applications. FRoG is a graphics processing unit (GPU)-based forward calculation tool developed at CNAO (Centro Nazionale di Adroterapia Oncologica) and at HIT (Heidelberg Ion Beam Therapy Center) for fast and accurate calculation of both physical and biological dose. FRoG calculation engine adopts a triple Gaussian parameterization for the description of the lateral dose distribution. FRoG provides dose, dose-averaged linear energy transfer, and biological dose-maps, -profiles, and -volume-histograms. For the benchmark of the FRoG calculation engine, using the clinical settings available at CNAO, spread-out Bragg peaks (SOBPs) and patient cases for both proton and carbon ion beams have been calculated and compared against FLUKA Monte Carlo (MC) predictions. In addition, FRoG patient-specific quality assurance (QA) has been performed for twenty-five proton and carbon ion fields. As a result, for protons, biological dose values, using a relative biological effectiveness (RBE) of 1.1, agree on average with MC within ~1% for both SOBPs and patient plans. For carbon ions, RBE-weighted dose (DRBE) agreement against FLUKA is within ~2.5% for the studied SOBPs and patient plans. Both MKM (Microdosimetric Kinetic Model) and LEM (Local Effect Model) DRBE are implemented and tested in FRoG to support the NIRS (National Institute of Radiological Sciences)-based to LEM-based biological dose conversion. FRoG matched the measured QA dosimetric data within ~2.0% for both particle species. The typical calculation times for patients ranged from roughly 1 to 4 min for proton beams and 3 to 6 min for carbon ions on a NVIDIA® GeForce® GTX 1080 Ti. This works demonstrates FRoG's potential to bolster clinical activity with proton and carbon ion beams at CNAO.
Collapse
|
12
|
Shrestha D, Qin N, Zhang Y, Kalantari F, Niu S, Jia X, Pompos A, Jiang S, Wang J. Iterative reconstruction with boundary detection for carbon ion computed tomography. Phys Med Biol 2018; 63:055002. [PMID: 29384493 DOI: 10.1088/1361-6560/aaac0f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
In heavy ion radiation therapy, improving the accuracy in range prediction of the ions inside the patient's body has become essential. Accurate localization of the Bragg peak provides greater conformity of the tumor while sparing healthy tissues. We investigated the use of carbon ions directly for computed tomography (carbon CT) to create the relative stopping power map of a patient's body. The Geant4 toolkit was used to perform a Monte Carlo simulation of the carbon ion trajectories, to study their lateral and angular deflections and the most likely paths, using a water phantom. Geant4 was used to create carbonCT projections of a contrast and spatial resolution phantom, with a cone beam of 430 MeV/u carbon ions. The contrast phantom consisted of cranial bone, lung material, and PMMA inserts while the spatial resolution phantom contained bone and lung material inserts with line pair (lp) densities ranging from 1.67 lp cm-1 through 5 lp cm-1. First, the positions of each carbon ion on the rear and front trackers were used for an approximate reconstruction of the phantom. The phantom boundary was extracted from this approximate reconstruction, by using the position as well as angle information from the four tracking detectors, resulting in the entry and exit locations of the individual ions on the phantom surface. Subsequent reconstruction was performed by the iterative algebraic reconstruction technique coupled with total variation minimization (ART-TV) assuming straight line trajectories for the ions inside the phantom. The influence of number of projections was studied with reconstruction from five different sets of projections: 15, 30, 45, 60 and 90. Additionally, the effect of number of ions on the image quality was investigated by reducing the number of ions/projection while keeping the total number of projections at 60. An estimation of carbon ion range using the carbonCT image resulted in improved range prediction compared to the range calculated using a calibration curve.
Collapse
Affiliation(s)
- Deepak Shrestha
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX 75235, United States of America
| | | | | | | | | | | | | | | | | |
Collapse
|
13
|
Qin N, Shen C, Tsai MY, Pinto M, Tian Z, Dedes G, Pompos A, Jiang SB, Parodi K, Jia X. Full Monte Carlo-Based Biologic Treatment Plan Optimization System for Intensity Modulated Carbon Ion Therapy on Graphics Processing Unit. Int J Radiat Oncol Biol Phys 2018; 100:235-243. [PMID: 29079118 DOI: 10.1016/j.ijrobp.2017.09.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2017] [Revised: 08/29/2017] [Accepted: 09/01/2017] [Indexed: 01/29/2023]
Abstract
PURPOSE One of the major benefits of carbon ion therapy is enhanced biological effectiveness at the Bragg peak region. For intensity modulated carbon ion therapy (IMCT), it is desirable to use Monte Carlo (MC) methods to compute the properties of each pencil beam spot for treatment planning, because of their accuracy in modeling physics processes and estimating biological effects. We previously developed goCMC, a graphics processing unit (GPU)-oriented MC engine for carbon ion therapy. The purpose of the present study was to build a biological treatment plan optimization system using goCMC. METHODS AND MATERIALS The repair-misrepair-fixation model was implemented to compute the spatial distribution of linear-quadratic model parameters for each spot. A treatment plan optimization module was developed to minimize the difference between the prescribed and actual biological effect. We used a gradient-based algorithm to solve the optimization problem. The system was embedded in the Varian Eclipse treatment planning system under a client-server architecture to achieve a user-friendly planning environment. We tested the system with a 1-dimensional homogeneous water case and 3 3-dimensional patient cases. RESULTS Our system generated treatment plans with biological spread-out Bragg peaks covering the targeted regions and sparing critical structures. Using 4 NVidia GTX 1080 GPUs, the total computation time, including spot simulation, optimization, and final dose calculation, was 0.6 hour for the prostate case (8282 spots), 0.2 hour for the pancreas case (3795 spots), and 0.3 hour for the brain case (6724 spots). The computation time was dominated by MC spot simulation. CONCLUSIONS We built a biological treatment plan optimization system for IMCT that performs simulations using a fast MC engine, goCMC. To the best of our knowledge, this is the first time that full MC-based IMCT inverse planning has been achieved in a clinically viable time frame.
Collapse
Affiliation(s)
- Nan Qin
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, Texas
| | - Chenyang Shen
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, Texas
| | - Min-Yu Tsai
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, Texas; Department of Computer Science and Information Engineering, National Taiwan University, Taipei, Taiwan
| | - Marco Pinto
- Department of Experimental Physics-Medical Physics, Ludwig-Maximilian University of Munich, Munich, Germany
| | - Zhen Tian
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, Texas
| | - Georgios Dedes
- Department of Experimental Physics-Medical Physics, Ludwig-Maximilian University of Munich, Munich, Germany
| | - Arnold Pompos
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, Texas
| | - Steve B Jiang
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, Texas
| | - Katia Parodi
- Department of Experimental Physics-Medical Physics, Ludwig-Maximilian University of Munich, Munich, Germany
| | - Xun Jia
- Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, Texas.
| |
Collapse
|
14
|
Giordanengo S, Manganaro L, Vignati A. Review of technologies and procedures of clinical dosimetry for scanned ion beam radiotherapy. Phys Med 2017; 43:79-99. [DOI: 10.1016/j.ejmp.2017.10.013] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/30/2017] [Revised: 09/23/2017] [Accepted: 10/18/2017] [Indexed: 12/17/2022] Open
|