No abstract

Monte Carlo (MC) simulation method plays an essential role in the refinement and development of positron emission tomography (PET) systems. However, most existing MC simulation packages suffer from long execution time for practical PET simulations. To fully address this issue, we developed and validated gPET, a graphics processing unit (GPU)-based MC simulation tool for PET. gPET was built on the NVidia CUDA platform. The simulation process was modularized into three functional parts and carried out by the GPU parallel threads: (1) source management, including positron decay, transport and annihilation; (2) gamma transport inside the phantom; and (3) signal detection and processing inside the detector. A hybrid of voxelized (for patient phantoms) and parametrized (for detectors) geometries were employed to sufficiently support particle navigations. Multiple inputs and outputs were available. Hence, a user can flexibly examine different aspects of a PET simulation. We evaluated the performance of gPET in three test cases with benchmark work from GATE8.0, in terms of the testing of the functional modules, the physics models used for gamma transport inside the detector, and the geometric configuration of an irregularly shaped PET detector. Both accuracy and efficiency were quantified. In all test cases, the differences between gPET and GATE for the coincidences with respect to the energy and crystal index distributions are below 3.18% and 2.54%, respectively. The speedup factor is 500 for gPET on a single Titan Xp GPU (1.58 GHz) over GATE8.0 on a single core of Intel i7-6850K CPU (3.6 GHz) for all test cases. In summary, gPET is an accurate and efficient MC simulation tool for PET.

Purpose An automated, in‐vivo system to detect patient anatomy changes and machine output was developed using novel analysis of in‐vivo electronic portal imaging device (EPID) images for every fraction of treatment on a Varian Halcyon. In‐vivo approach identifies errors that go undetected by routine quality assurance (QA) to compliment daily machine performance check (MPC), with minimal physicist workload. Methods Images for all fractions treated on a Halcyon were automatically downloaded and analyzed at the end of treatment day. For image analysis, compared to first fraction, the mean difference of high‐dose region of interest is calculated. This metric has shown to predict changes in planning treatment volume (PTV) mean dose. Flags are raised for: (Type‐A) treatment fraction whose mean difference exceeds 10%, to protect against large errors, and (Type‐B) patients with three consecutive fractions with mean exceeding ±3%, to protect against systematic trends. If a threshold is exceeded, a physicist is e‐mailed, a report for flagged patients, for investigation. To track machine output changes, for all patients treated on a day, the average and standard deviations are uploaded to a QA portal, along with the reviewed MPC, ensuring comprehensive QA for the Halcyon. To guide clinical implementation, a retrospective study from November 2017 till December 2020 was conducted, which grouped errors by treatment site. This framework has been used prospectively since January 2021. Results From retrospective data of 1633 patients (35 759 fractions), no Type‐A errors were found and only 45 patients (2.76%) had Type‐B errors. These Type‐B deviations were due to head‐and‐neck weight loss. For 6 months of prospective use (345 patients), 13 patients (3.7%) had Type‐B errors and no Type‐A errors. Conclusions This automated system protects against errors that can occur in vivo to provide a more comprehensive QA. This fully automated framework can be implemented in other centers with a Halcyon, requiring a desktop computer and analysis scripts.

No abstract

Range uncertainty remains a big concern in particle therapy, as it may cause target dose degradation and normal tissue overdosing. Positron emission tomography (PET) and prompt gamma imaging (PGI) are two promising modalities for range verification. However, the relatively long acquisition time of PET and the relatively low yield of PGI pose challenges for real-time range verification. In this paper, we explore using the primary Carbon-11 (C-11) ion beams to enhance the gamma yield compared to the primary C-12 ion beams to improve PET and PGI by using Monte Carlo simulations of water and PMMA phantoms at four incident energies (95, 200, 300, and 430 MeV u−1). Prompt gammas (PGs) and annihilation gammas (AGs) were recorded for post-processing to mimic PGI and PET imaging, respectively. We used both time-of-flight (TOF) and energy selections for PGI, which boosted the ratio of PGs to background neutrons to 2.44, up from 0.87 without the selections. At the lowest incident energy (100 MeVu-1), PG yield from C-11 was 0.82 times of that from C-12, while AG yield from C-11 was 6 ∼ 11 folds higher than from C-12 in PMMA. At higher energies, PG differences between C-11 and C-12 were much smaller, while AG yield from C-11 was 30%∼90% higher than from C-12 using minute-acquisition. With minute-acquisition, the AG depth distribution of C-11 showed a sharp peak coincident with the Bragg peak due to the decay of the primary C-11 ions, but that of C-12 had no such one. The high AG yield and distinct peaks could lead to more precise range verification of C-11 than C-12. These results demonstrate that using C-11 ion beams for potentially combined PGI and PET has great potential to improve online single-spot range verification accuracy and precision.