For the first time, range verification based on prompt gamma imaging was applied for a clinical proton treatment. With the translation from basic physics experiments into clinical operation, the potential to improve the precision of particle therapy with this technique has increased considerably.

An independent assessment of the dose delivery in ion therapy can be performed using positron emission tomography (PET). For that a distribution of positron emitters which appear as the result of interaction between ions of the therapeutic beam and the irradiated tissue is measured during or after the irradiation. Three concepts for PET monitoring implemented in various therapy facilities are considered in this paper. The in-beam PET concept relies on the PET measurement performed simultaneously to the irradiation by means of a PET scanner which is completely integrated into the irradiation site. The in-room PET concept allows measurement immediately after irradiation by a standalone PET scanner which is installed very close to the irradiation site. In the off-line PET scenario the measurement is performed by means of a standalone PET/CT scanner 10-30 min after the irradiation. These three concepts were evaluated according to image quality criteria, integration costs, and their influence onto the workflow of radiotherapy. In-beam PET showed the best performance. However, the integration costs were estimated as very high for this modality. Moreover, the performance of in-beam PET depends heavily on type and duty cycle of the accelerator. The in-room PET is proposed for planned therapy facilities as a good compromise between the quality of measured data and integration efforts. For facilities which are close to the nuclear medicine departments off-line PET can be suggested under several circumstances.

Ion beam therapy promises enhanced tumour coverage compared to conventional radiotherapy, but particle range uncertainties significantly blunt the achievable precision. Experimental tools for range verification in real-time are not yet available in clinical routine. The prompt gamma ray timing method has been recently proposed as an alternative to collimated imaging systems. The detection times of prompt gamma rays encode essential information about the depth-dose profile thanks to the measurable transit time of ions through matter. In a collaboration between OncoRay, Helmholtz-Zentrum Dresden-Rossendorf and IBA, the first test at a clinical proton accelerator (Westdeutsches Protonentherapiezentrum Essen, Germany) with several detectors and phantoms is performed. The robustness of the method against background and stability of the beam bunch time profile is explored, and the bunch time spread is characterized for different proton energies. For a beam spot with a hundred million protons and a single detector, range differences of 5 mm in defined heterogeneous targets are identified by numerical comparison of the spectrum shape. For higher statistics, range shifts down to 2 mm are detectable. A proton bunch monitor, higher detector throughput and quantitative range retrieval are the upcoming steps towards a clinically applicable prototype. In conclusion, the experimental results highlight the prospects of this straightforward verification method at a clinical pencil beam and settle this novel approach as a promising alternative in the field of in vivo dosimetry.

Proton and ion beam therapies become increasingly relevant in radiation therapy. To fully exploit the potential of this irradiation technique and to achieve maximum target volume conformality, the verification of particle ranges is highly desirable. Many research activities focus on the measurement of the spatial distributions of prompt gamma rays emitted during irradiation. A passively collimating knife-edge slit camera is a promising option to perform such measurements. In former publications, the feasibility of accurate detection of proton range shifts in homogeneous targets could be shown with such a camera. We present slit camera measurements of prompt gamma depth profiles in inhomogeneous targets. From real treatment plans and their underlying CTs, representative beam paths are selected and assembled as one-dimensional inhomogeneous targets built from tissue equivalent materials. These phantoms have been irradiated with monoenergetic proton pencil beams. The accuracy of range deviation estimation as well as the detectability of range shifts is investigated in different scenarios. In most cases, range deviations can be detected within less than 2 mm. In close vicinity to low-density regions, range detection is challenging. In particular, a minimum beam penetration depth of 7 mm beyond a cavity is required for reliable detection of a cavity filling with the present setup. Dedicated data post-processing methods may be capable of overcoming this limitation.

A prompt gamma (PG) slit camera prototype recently demonstrated that Bragg Peak position in a clinical proton scanned beam could be measured with 1-2 mm accuracy by comparing an expected PG detection profile to a measured one. The computation of the expected PG detection profile in the context of a clinical framework is challenging but must be solved before clinical implementation. Obviously, Monte Carlo methods (MC) can simulate the expected PG profile but at prohibitively long calculation times. We implemented a much faster method that is based on analytical processing of precomputed MC data that would allow practical evaluation of this range monitoring approach in clinical conditions. Reference PG emission profiles were generated with MC simulations (PENH) in targets consisting of either (12)C, (14)N, (16)O, (31)P or (40)Ca, with 10% of (1)H. In a given geometry, the local PG emission can then be derived by adding the contribution of each element, according to the local energy of the proton obtained by continuous slowing down approximation and the local composition. The actual incident spot size is taken into account using an optical model fitted to measurements and by super sampling the spot with several rays (up to 113). PG transport in the patient/camera geometries and the detector response are modelled by convolving the PG production profile with a transfer function. The latter is interpolated from a database of transfer functions fitted to MC data (PENELOPE) generated for a photon source in a cylindrical phantom with various radiuses and a camera placed at various positions. As a benchmark, the analytical model was compared to MC and experiments in homogeneous and heterogeneous phantoms. Comparisons with MC were also performed in a thoracic CT. For all cases, the analytical model reproduced the prediction of the position of the Bragg peak computed with MC within 1 mm for the camera in nominal configuration. When compared to measurements, the shape of the profiles was well reproduced and agreement for the estimation of the position of the Bragg peak was within 2.7 mm on average (1.4 mm standard deviation). On a non-optimized MATLAB code, computation time with the analytical model is between 0.3 to 10 s depending on the number of rays simulated per spot. The analytical model can be further used to determine which spots are the best candidates to evaluate the range in clinical conditions and eventually correct for over- and under-shoots depending on the acquired PG profiles.