A new approach to integrate GPU-based Monte Carlo simulation into inverse treatment plan optimization for proton therapy

Monte Carlo (MC)-based spot dose calculation is highly desired for inverse treatment planning in proton therapy because of its accuracy. Recent studies on biological optimization have also indicated the use of MC methods to compute relevant quantities of interest, e.g. linear energy transfer. Although GPU-based MC engines have been developed to address inverse optimization problems, their efficiency still needs to be improved. Also, the use of a large number of GPUs in MC calculation is not favorable for clinical applications. The previously proposed adaptive particle sampling (APS) method can improve the efficiency of MC-based inverse optimization by using the computationally expensive MC simulation more effectively. This method is more efficient than the conventional approach that performs spot dose calculation and optimization in two sequential steps. In this paper, we propose a computational library to perform MC-based spot dose calculation on GPU with the APS scheme. The implemented APS method performs a non-uniform sampling of the particles from pencil beam spots during the optimization process, favoring those from the high intensity spots. The library also conducts two computationally intensive matrix-vector operations frequently used when solving an optimization problem. This library design allows a streamlined integration of the MC-based spot dose calculation into an existing proton therapy inverse planning process. We tested the developed library in a typical inverse optimization system with four patient cases. The library achieved the targeted functions by supporting inverse planning in various proton therapy schemes, e.g. single field uniform dose, 3D intensity modulated proton therapy, and distal edge tracking. The efficiency was 41.6  ±  15.3% higher than the use of a GPU-based MC package in a conventional calculation scheme. The total computation time ranged between 2 and 50 min on a single GPU card depending on the problem size.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .

[3]  A. Lomax,et al.  Intensity modulation methods for proton radiotherapy. , 1999, Physics in medicine and biology.

[4]  Uwe Oelfke,et al.  Optimization of radiobiological effects in intensity modulated proton therapy. , 2005, Medical physics.

[5]  E Pedroni,et al.  Monte Carlo dose calculations for spot scanned proton therapy , 2005, Physics in medicine and biology.

[6]  H Edwin Romeijn,et al.  A computational implementation and comparison of several intensity modulated proton therapy treatment planning algorithms. , 2008, Medical physics.

[7]  Katia Parodi,et al.  Clinical implementation of full Monte Carlo dose calculation in proton beam therapy , 2008, Physics in medicine and biology.

[8]  Steve B. Jiang,et al.  GPU-based ultrafast IMRT plan optimization , 2009, Physics in medicine and biology.

[9]  Wayne D Newhauser,et al.  Monte Carlo fast dose calculator for proton radiotherapy: application to a voxelized geometry representing a patient with prostate cancer. , 2009, Physics in medicine and biology.

[10]  Uwe Oelfke,et al.  Application of constant vs. variable relative biological effectiveness in treatment planning of intensity-modulated proton therapy. , 2011, International journal of radiation oncology, biology, physics.

[11]  Harald Paganetti,et al.  Variations in linear energy transfer within clinical proton therapy fields and the potential for biological treatment planning. , 2011, International journal of radiation oncology, biology, physics.

[12]  H Paganetti,et al.  TOPAS: an innovative proton Monte Carlo platform for research and clinical applications. , 2012, Medical physics.

[13]  H. Paganetti.  Range uncertainties in proton therapy and the role of Monte Carlo simulations , 2012, Physics in medicine and biology.

[14]  Steve B. Jiang,et al.  GPU-based fast Monte Carlo dose calculation for proton therapy , 2012, Physics in medicine and biology.

[15]  Wei Liu,et al.  Influence of robust optimization in intensity-modulated proton therapy with different dose delivery techniques. , 2012, Medical physics.

[16]  David Craft,et al.  Linear energy transfer-guided optimization in intensity modulated proton therapy: feasibility study and clinical potential. , 2013, International journal of radiation oncology, biology, physics.

[17]  K Parodi,et al.  A Monte Carlo-based treatment planning tool for proton therapy , 2013, Physics in medicine and biology.

[18]  Bengt K Lind,et al.  A model for the relative biological effectiveness of protons: The tissue specific parameter α/β of photons is a predictor for the sensitivity to LET changes , 2013, Acta oncologica.

[19]  Chris Beltran,et al.  A GPU-accelerated and Monte Carlo-based intensity modulated proton therapy optimization system. , 2014, Medical physics.

[20]  G Cuttone,et al.  A Monte Carlo study for the calculation of the average linear energy transfer (LET) distributions for a clinical proton beam line and a radiobiological carbon ion beam line , 2014, Physics in medicine and biology.

[21]  Jan Seuntjens,et al.  Latent uncertainties of the precalculated track Monte Carlo method. , 2015, Medical physics.

[22]  H. Tseung,et al.  A fast GPU-based Monte Carlo simulation of proton transport with detailed modeling of nonelastic interactions. , 2014, Medical physics.

[23]  Steve B. Jiang,et al.  A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy , 2015, Physics in medicine and biology.