Search

Objective. Verification of delivered proton therapy treatments is essential for reaping the many benefits of the modality, with the most widely proposed in vivo verification 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 σ 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.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., (Creative Commons Attribution license.)

In vivo dose delivery verification in proton therapy can be performed by positron emission tomography (PET) of the positron-emitting nuclei produced by the proton beam in the patient. A PET scanner installed in the treatment position of a proton therapy facility that takes data with the beam on will see very short-lived nuclides as well as longer-lived nuclides. The most important short-lived nuclide for proton therapy is 12 N (Dendooven et al 2015 Phys. Med. Biol. 60 8923-47), which has a half-life of 11 ms. The results of a proof-of-principle experiment of beam-on PET imaging of short-lived 12 N nuclei are presented. The Philips Digital Photon Counting Module TEK PET system was used, which is based on LYSO scintillators mounted on digital SiPM photosensors. A 90 MeV proton beam from the cyclotron at KVI-CART was used to investigate the energy and time spectra of PET coincidences during beam-on. Events coinciding with proton bunches, such as prompt gamma rays, were removed from the data via an anti-coincidence filter with the cyclotron RF. The resulting energy spectrum allowed good identification of the 511 keV PET counts during beam-on. A method was developed to subtract the long-lived background from the 12 N image by introducing a beam-off period into the cyclotron beam time structure. We measured 2D images and 1D profiles of the 12 N distribution. A range shift of 5 mm was measured as 6  ±  3 mm using the 12 N profile. A larger, more efficient, PET system with a higher data throughput capability will allow beam-on 12 N PET imaging of single spots in the distal layer of an irradiation with an increased signal-to-background ratio and thus better accuracy. A simulation shows that a large dual panel scanner, which images a single spot directly after it is delivered, can measure a 5 mm range shift with millimeter accuracy: 5.5  ±  1.1 mm for 1  ×  10 8 protons and 5.2  ±  0.5 mm for 5  ×  10 8 protons. This makes fast and accurate feedback on the dose delivery during treatment possible.