Recent developments and comprehensive evaluations of a GPU-based Monte Carlo package for proton therapy.
Monte Carlo (MC) simulation is commonly considered as the most accurate dose calculation method for proton therapy. Aiming at achieving fast MC dose calculations for clinical applications, we have previously developed a graphics-processing unit (GPU)-based MC tool, gPMC. In this paper, we report our recent updates on gPMC in terms of its accuracy, portability, and functionality, as well as comprehensive tests on this tool. The new version, gPMC v2.0, was developed under the OpenCL environment to enable portability across different computational platforms. Physics models of nuclear interactions were refined to improve calculation accuracy. Scoring functions of gPMC were expanded to enable tallying particle fluence, dose deposited by different particle types, and dose-averaged linear energy transfer (LETd). A multiple counter approach was employed to improve efficiency by reducing the frequency of memory writing conflict at scoring. For dose calculation, accuracy improvements over gPMC v1.0 were observed in both water phantom cases and a patient case. For a prostate cancer case planned using high-energy proton beams, dose discrepancies in beam entrance and target region seen in gPMC v1.0 with respect to the gold standard tool for proton Monte Carlo simulations (TOPAS) results were substantially reduced and gamma test passing rate (1%/1 mm) was improved from 82.7%-93.1%. The average relative difference in LETd between gPMC and TOPAS was 1.7%. The average relative differences in the dose deposited by primary, secondary, and other heavier particles were within 2.3%, 0.4%, and 0.2%. Depending on source proton energy and phantom complexity, it took 8-17 s on an AMD Radeon R9 290x GPU to simulate [Formula: see text] source protons, achieving less than [Formula: see text] average statistical uncertainty. As the beam size was reduced from 10  ×  10 cm2to 1  ×  1 cm2, the time on scoring was only increased by 4.8% with eight counters, in contrast to a 40% increase using only one counter. With the OpenCL environment, the portability of gPMC v2.0 was enhanced. It was successfully executed on different CPUs and GPUs and its performance on different devices varied depending on processing power and hardware structure.
Department of Radiation Oncology, University of Texas Southwestern Medical Center, Dallas, TX 75390, USA.
Nan
N
Qin
Pablo
P
Botas
Drosoula
D
Giantsoudi
Jan
J
Schuemann
Zhen
Z
Tian
Steve B
SB
Jiang
Harald
H
Paganetti
Xun
X
Jia
Phys Med Biol. 2009 Jan 7;54(1):N21-8
19075361
Phys Med Biol. 2014 Nov 21;59(22):R419-72
25361443
Phys Med Biol. 2014 Feb 21;59(4):R151-82
24486639
Med Phys. 2012 Nov;39(11):6818-37
23127075
Phys Med Biol. 2014 Aug 7;59(15):4007-31
24990623
Phys Med Biol. 2015 Mar 21;60(6):2257-69
25715661
Phys Med Biol. 2005 Mar 7;50(5):1001-10
15798272
Phys Med Biol. 2015 Jul 7;60(13):5007-18
26061350
Phys Med Biol. 2015 Jul 21;60(14):N283-91
26147442
Med Phys. 2004 Aug;31(8):2263-73
15377093
Phys Med Biol. 2015 Oct 7;60(19):7419-35
26352012
Phys Med Biol. 2015 Nov 7;60(21):8399-416
26459756
Phys Med Biol. 2008 Sep 7;53(17):4825-53
18701772
Med Phys. 2015 Jun;42(6):2967-78
26127050
Phys Med Biol. 2003 May 21;48(10):1277-88
12812446
Phys Med Biol. 2011 Mar 7;56(5):1431-41
21317484
Phys Med Biol. 2012 Dec 7;57(23):7783-97
23128424
Phys Med Biol. 2015 Jul 7;60(13):5053-70
26061666
Med Phys. 2013 Dec;40(12):121719
24320505
Int J Radiat Oncol Biol Phys. 2013 Sep 1;87(1):216-22
23790771
Int J Radiat Oncol Biol Phys. 2016 Aug 1;95(5):1535-1543
27325476
Med Phys. 2014 Dec;41(12):121707
25471954
Int J Radiat Oncol Biol Phys. 2015 Aug 1;92(5):1157-1164
26025779