GPU accelerated Monte Carlo scoring of positron emitting isotopes produced during proton therapy for PET verification.

McNamara, Keegan; Schiavi, Angelo; Borys, Damian; Brzezinski, Karol; Gajewski, Jan; Kopeć, Renata; Rucinski, Antoni; Skóra, Tomasz; Makkar, Shubhangi; Hrbacek, Jan; Weber, Damien C; Lomax, Antony J; Winterhalter, Carla (2022). GPU accelerated Monte Carlo scoring of positron emitting isotopes produced during proton therapy for PET verification. Physics in medicine and biology, 67(24) Institute of Physics Publishing IOP 10.1088/1361-6560/aca515

Preview

Text McNamara_2022_Phys._Med._Biol._67_244001.pdf - Published Version Available under License Creative Commons: Attribution (CC-BY). Download (1MB) | Preview

Objective.Verification of delivered proton therapy treatments is essential for reaping the many benefits of the modality, with the most widely proposedin vivoverification technique being the imaging of positron emitting isotopes generated in the patient during treatment using positron emission tomography (PET). The purpose of this work is to reduce the computational resources and time required for simulation of patient activation during proton therapy using the GPU accelerated Monte Carlo code FRED, and to validate the predicted activity against the widely used Monte Carlo code GATE.Approach.We implement a continuous scoring approach for the production of positron emitting isotopes within FRED version 5.59.9. We simulate treatment plans delivered to 95 head and neck patients at Centrum Cyklotronowe Bronowice using this GPU implementation, and verify the accuracy using the Monte Carlo toolkit GATE version 9.0.Main results.We report an average reduction in computational time by a factor of 50 when using a local system with 2 GPUs as opposed to a large compute cluster utilising between 200 to 700 CPU threads, enabling simulation of patient activity within an average of 2.9 min as opposed to 146 min. All simulated plans are in good agreement across the two Monte Carlo codes. The two codes agree within a maximum of 0.95σon a voxel-by-voxel basis for the prediction of 7 different isotopes across 472 simulated fields delivered to 95 patients, with the average deviation over all fields being 6.4 × 10-3σ.Significance.The implementation of activation calculations in the GPU accelerated Monte Carlo code FRED provides fast and reliable simulation of patient activation following proton therapy, allowing for research and development of clinical applications of range verification for this treatment modality using PET to proceed at a rapid pace.

Interest & Impact

Downloads

70 since deposited on 22 Dec 2022
30 in the past 12 months

Citations

5 Citations in Web of Science ®
6 Citations in Scopus

Search

in Google Scholar™

Services

Actions (login required)

Edit item

Actions (login required)

Edit item