Quantitative assessment of radionuclide production yields in in-beam and offline PET measurements at different proton irradiation facilities.

Objective. Reliable radionuclide production yield data are a prerequisite for positron-emission-tomography (PET) basedin vivoproton treatment verification. In this context, activation data acquired at two different treatment facilities with different imaging systems were analyzed to provide experimentally determined radionuclide yields in thick targets and were compared with each other to investigate the impact of the respective imaging technique.Approach.Homogeneous thick targets (PMMA, gelatine, and graphite) were irradiated with mono-energetic proton pencil-beams at two distinct energies. Material activation was measured (i)in-beamduring and after beam delivery with a double-head prototype PET camera and (ii)offlineshortly after beam delivery with a commercial full-ring PET/CT scanner. Integral as well as depth-resolvedß+-emitter yields were determined for the dominant positron-emitting radionuclides11C,15O,13N and (in-beamonly)10C.In-beamdata were used to investigate the qualitative impact of different monitoring time schemes on activity depth profiles and their quantitative impact on count rates and total activity.Main results.Production yields measured with thein-beamcamera were comparable to or higher compared to respectiveofflineresults. Depth profiles of radionuclide-specific yields obtained from thedouble-headcamera showed qualitative differences to data acquired with thefull-ringcamera with a more convex profile shape. Considerable impact of the imaging timing scheme on the activity profile was observed for gelatine only with a range variation of up to 3.5 mm. Evaluation of the coincidence rate and the total number of observed events in the considered workflows confirmed a strongly decreasing rate in targets with a large oxygen fraction.Significance. The observed quantitative and qualitative differences between the datasets underline the importance of a thorough system commissioning. Due to the lack of reliable cross-section data, in-house phantom measurements are still considered a gold standard for careful characterization of the system response and to ensure a reliable beam range verification.

In-beam PET at clinical proton beams with pile-up rejection.

Positron emission tomography (PET) is a means of imaging the ß+-activity produced by the radiation field in ion beam therapy and therefore for treatment verification. Prompt γ-rays that are emitted during beam application challenge the detectors and electronics of PET systems, since those are designed for low and medium count rates. Typical PET detectors operated according to a modified Anger principle suffer from multiple events at high rates. Therefore, in-beam PET systems using such detectors rely on a synchronization of beam status and measurement to reject deteriorated data. In this work, a method for pile-up rejection is applied to conventional Anger logic block detectors. It allows for an in-beam data acquisition without further synchronization. Though cyclotrons produce a continuous wave beam, the radiation field shaping technique introduces breaks in the application. Time regimes mimicking synchrotrons as well as cyclotron based ones using double-scattering or pencil beam scanning field shaping at dose rates of 0.5, 1.0 and 2.0Gy/min were investigated. Two types of inhomogeneous phantoms were imaged. The first one simulates cavity structures, the other one mimics a static lung irradiation. It could be shown that, depending on the dose rate and the beam time structure, in-beam measurement including a few seconds decay time only, yield images which revealed all inhomogeneities in the phantoms. This technique can be the basis for the development of an in-beam PET system with traditional detectors and off-the-shelf electronics.

Feasibility Study on Cardiac Arrhythmia Ablation Using High-Energy Heavy Ion Beams.

High-energy ion beams are successfully used in cancer therapy and precisely deliver high doses of ionizing radiation to small deep-seated target volumes. A similar noninvasive treatment modality for cardiac arrhythmias was tested here. This study used high-energy carbon ions for ablation of cardiac tissue in pigs. Doses of 25, 40, and 55 Gy were applied in forced-breath-hold to the atrioventricular junction, left atrial pulmonary vein junction, and freewall left ventricle of intact animals. Procedural success was tracked by (1.) in-beam positron-emission tomography (PET) imaging; (2.) intracardiac voltage mapping with visible lesion on ultrasound; (3.) lesion outcomes in pathohistolgy. High doses (40-55 Gy) caused slowing and interruption of cardiac impulse propagation. Target fibrosis was the main mediator of the ablation effect. In irradiated tissue, apoptosis was present after 3, but not 6 months. Our study shows feasibility to use high-energy ion beams for creation of cardiac lesions that chronically interrupt cardiac conduction.

Compton Camera and Prompt Gamma Ray Timing: Two Methods for In Vivo Range Assessment in Proton Therapy.

Proton beams are promising means for treating tumors. Such charged particles stop at a defined depth, where the ionization density is maximum. As the dose deposit beyond this distal edge is very low, proton therapy minimizes the damage to normal tissue compared to photon therapy. Nevertheless, inherent range uncertainties cast doubts on the irradiation of tumors close to organs at risk and lead to the application of conservative safety margins. This constrains significantly the potential benefits of protons over photons. In this context, several research groups are developing experimental tools for range verification based on the detection of prompt gammas, a nuclear by-product of the proton irradiation. At OncoRay and Helmholtz-Zentrum Dresden-Rossendorf, detector components have been characterized in realistic radiation environments as a step toward a clinical Compton camera. On the one hand, corresponding experimental methods and results obtained during the ENTERVISION training network are reviewed. On the other hand, a novel method based on timing spectroscopy has been proposed as an alternative to collimated imaging systems. The first tests of the timing method at a clinical proton accelerator are summarized, its applicability in a clinical environment for challenging the current safety margins is assessed, and the factors limiting its precision are discussed.

First test of the prompt gamma ray timing method with heterogeneous targets at a clinical proton therapy facility.

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.

Systematic analysis on the achievable accuracy of PT-PET through automated evaluation techniques.

INTRODUCTION: Particle Therapy Positron Emission Tomography (PT-PET) is currently the only clinically applied method for in vivo verification of ion-beam radiotherapy during or close in time to the treatment. Since a direct deduction of the delivered dose from the measured activity is not feasible, images are compared to a reference distribution. The achievable accuracy of two image analysis approaches was investigated by means of reproducible phantom benchmark tests. This is an objective method that excludes patient related factors of influence. MATERIAL AND METHODS: Two types of phantoms were designed to produce well defined deviations in the activity distributions. Pure range differences were simulated using the first phantom type while the other emulated cavity structures. The phantoms were irradiated with (12)C-ions. PT-PET measurements were performed by means of a camera system installed at the beamline. Different measurement time scenarios were investigated, assuming a PET scanner directly at the irradiation site or placed within the treatment room. The images were analyzed by means of the Pearson Correlation Coefficient (PCC) and a range calculation algorithm combined with a dedicated cavity filling detection method. RESULTS: Range differences could be measured with an error of less than 2 mm. The range comparison algorithm yielded slightly better results than the PCC method. The filling of a cavity structure could be safely detected if its inner diameter was at least 5 mm. CONCLUSION: Both approaches evaluate the PT-PET data in an objective way and deliver promising results for in-beam and in-room PET for clinical realistic dose rates.

Range assessment in particle therapy based on prompt γ-ray timing measurements.

Proton and ion beams open up new vistas for the curative treatment of tumors, but adequate technologies for monitoring the compliance of dose delivery with treatment plans in real time are still missing. Range assessment, meaning the monitoring of therapy-particle ranges in tissue during dose delivery (treatment), is a continuous challenge considered a key for tapping the full potential of particle therapies. In this context the paper introduces an unconventional concept of range assessment by prompt-gamma timing (PGT), which is based on an elementary physical effect not considered so far: therapy particles penetrating tissue move very fast, but still need a finite transit time--about 1-2 ns in case of protons with a 5-20 cm range--from entering the patient's body until stopping in the target volume. The transit time increases with the particle range. This causes measurable effects in PGT spectra, usable for range verification. The concept was verified by proton irradiation experiments at the AGOR cyclotron, KVI-CART, University of Groningen. Based on the presented kinematical relations, we describe model calculations that very precisely reproduce the experimental results. As the clinical treatment conditions entail measurement constraints (e.g. limited treatment time), we propose a setup, based on clinical irradiation conditions, capable of determining proton range deviations within a few seconds of irradiation, thus allowing for a fast safety survey. Range variations of 2 mm are expected to be clearly detectable.

Dependence of simulated positron emitter yields in ion beam cancer therapy on modeling nuclear fragmentation.

In ion beam cancer therapy, range verification in patients using positron emission tomography (PET) requires the comparison of measured with simulated positron emitter yields. We found that (1) changes in modeling nuclear interactions strongly affected the positron emitter yields and that (2) Monte Carlo simulations with SHIELD-HIT10Areasonably matched the most abundant PET isotopes (11)C and (15)O. We observed an ion-energy (i.e., depth) dependence of the agreement between SHIELD-HIT10Aand measurement. Improved modeling requires more accurate measurements of cross-section values.

Automated evaluation of setup errors in carbon ion therapy using PET: feasibility study.

PURPOSE: To investigate the possibility of detecting patient mispositioning in carbon-ion therapy with particle therapy positron emission tomography (PET) in an automated image registration based manner. METHODS: Tumors in the head and neck (H&N), pelvic, lung, and brain region were investigated. Biologically optimized carbon ion treatment plans were created with TRiP98. From these treatment plans, the reference ß(+)-activity distributions were calculated using a Monte Carlo simulation. Setup errors were simulated by shifting or rotating the computed tomography (CT). The expected ß(+) activity was calculated for each plan with shifts. Finally, the reference particle therapy PET images were compared to the "shifted" ß(+)-activity distribution simulations using the Pearson's correlation coefficient (PCC). To account for different PET monitoring options the inbeam PET was compared to three different inroom scenarios. Additionally, the dosimetric effects of the CT misalignments were investigated. RESULTS: The automated PCC detection of patient mispositioning was possible in the investigated indications for cranio-caudal shifts of 4 mm and more, except for prostate tumors. In the rather homogeneous pelvic region, the generated ß(+)-activity distribution of the reference and compared PET image were too much alike. Thus, setup errors in this region could not be detected. Regarding lung lesions the detection strongly depended on the exact tumor location: in the center of the lung tumor misalignments could be detected down to 2 mm shifts while resolving shifts of tumors close to the thoracic wall was more challenging. Rotational shifts in the H&N and lung region of +6° and more could be detected using inroom PET and partly using inbeam PET. Comparing inroom PET to inbeam PET no obvious trend was found. However, among the inroom scenarios a longer measurement time was found to be advantageous. CONCLUSIONS: This study scopes the use of various particle therapy PET verification techniques in four indications. The automated detection of patients' setup errors was investigated in a broad accumulation of data sets. The evaluation of introduced setup errors is performed automatically, which is of utmost importance to introduce highly required particle therapy monitoring devices into the clinical routine.

Analysis of metabolic washout of positron emitters produced during carbon ion head and neck radiotherapy.

PURPOSE: Particle Therapy Positron Emission Tomography (PT-PET) is a suitable method for verification of therapeutic dose delivery by measurements of irradiation-induced ß(+)-activity. Due to metabolic processes in living tissue ß(+)-emitters can be removed from the place of generation. This washout is a limiting factor for image quality. The purpose of this study is to investigate whether a washout model obtained by animal experiments is applicable to patient data. METHODS: A model for the washout has been developed by Mizuno et al. [Phys. Med. Biol. 48(15), 2269-2281 (2003)] and Tomitani et al. [Phys. Med. Biol. 48(7), 875-889 (2003)]. It is based upon measurements in a rabbit in living and dead conditions. This model was modified and applied to PET data acquired during the experimental therapy project at GSI Helmholtzzentrum für Schwerionenforschung Darmstadt, Germany. Three components are expected: A fast one with a half life of 2 s, a medium one in the range of 2-3 min, and a slow component of the order of 2-3 h. Ten patients were selected randomly for investigation of the fast component. To analyze the other two components, 12 one-of-a-kind measurements from a single volunteer patient are available. RESULTS: A fast washout on the time scale of a few seconds was not observed in the patient data. The medium processes showed a mean half life of 155.7 ± 4.6 s. This is in the expected range. Fractions of the activity not influenced by the washout were found. CONCLUSIONS: On the time scale of an in-beam or in-room measurement only the medium-time washout processes play a remarkable role. A slow component may be neglected if the measurements do not exceed 20 min after the end of the irradiation. The fast component is not observed due to the low relative blood filled volume in the brain.

Comparison of PHITS, GEANT4, and HIBRAC simulations of depth-dependent yields of ß(+)-emitting nuclei during therapeutic particle irradiation to measured data.

For quality assurance in particle therapy, a non-invasive, in vivo range verification is highly desired. Particle therapy positron-emission-tomography (PT-PET) is the only clinically proven method up to now for this purpose. It makes use of the ß(+)-activity produced during the irradiation by the nuclear fragmentation processes between the therapeutic beam and the irradiated tissue. Since a direct comparison of ß(+)-activity and dose is not feasible, a simulation of the expected ß(+)-activity distribution is required. For this reason it is essential to have a quantitatively reliable code for the simulation of the yields of the ß(+)-emitting nuclei at every position of the beam path. In this paper results of the three-dimensional Monte-Carlo simulation codes PHITS, GEANT4, and the one-dimensional deterministic simulation code HIBRAC are compared to measurements of the yields of the most abundant ß(+)-emitting nuclei for carbon, lithium, helium, and proton beams. In general, PHITS underestimates the yields of positron-emitters. With GEANT4 the overall most accurate results are obtained. HIBRAC and GEANT4 provide comparable results for carbon and proton beams. HIBRAC is considered as a good candidate for the implementation to clinical routine PT-PET.

Using statistical measures for automated comparison of in-beam PET data.

PURPOSE: Positron emission tomography (PET) is considered to be the state of the art technique to monitor particle therapy in vivo. To evaluate the beam delivery the measured PET image is compared to a predicted ß(+)-distribution. Nowadays the range assessment is performed by a group of experts via visual inspection. This procedure is rather time consuming and requires well trained personnel. In this study an approach is presented to support human decisions in an automated and objective way. METHODS: The automated comparison presented uses statistical measures, namely, Pearson's correlation coefficient (PCC), to detect ion beam range deviations. The study is based on 12 in-beam PET patient data sets recorded at GSI and 70 artificial beam range modifications per data set. The range modifications were 0, 4, 6, and 10 mm water equivalent path length (WEPL) in positive and negative beam directions. The reference image to calculate the PCC was both an unmodified simulation of the activity distribution (Test 1) and a measured in-beam PET image (Test 2). Based on the PCCs sensitivity and specificity were calculated. Additionally the difference between modified and unmodified data sets was investigated using the Wilcoxon rank sum test. RESULTS: In Test 1 a sensitivity and specificity over 90% was reached for detecting modifications of ±10 and ±6 mm WEPL. Regarding Test 2 a sensitivity and specificity above 80% was obtained for modifications of ±10 and -6 mm WEPL. The limitation of the method was around 4 mm WEPL. CONCLUSIONS: The results demonstrate that the automated comparison using PCC provides similar results in terms of sensitivity and specificity compared to visual inspections of in-beam PET data. Hence the method presented in this study is a promising and effective approach to improve the efficiency in the clinical workflow in terms of particle therapy monitoring by means of PET.

Implementation and workflow for PET monitoring of therapeutic ion irradiation: a comparison of in-beam, in-room, and off-line techniques.

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.

On the effectiveness of ion range determination from in-beam PET data.

At present, in-beam positron emission tomography (PET) is the only method for in vivo and in situ range verification in ion therapy. At the GSI Helmholtzzentrum für Schwerionenforschung GmbH (GSI) Darmstadt, Germany, a unique in-beam PET installation has been operated from 1997 until the shut down of the carbon ion therapy facility in 2008. Therapeutic irradiation by means of (12)C ion beams of more than 400 patients have been monitored. In this paper a first quantitative study on the accuracy of the in-beam PET method to detect range deviations between planned and applied treatment in clinically relevant situations using simulations based on clinical data is presented. Patient treatment plans were used for performing simulations of positron emitter distributions. For each patient a range difference of + or - 6 mm in water was applied and compared to simulations without any changes. The comparisons were performed manually by six experienced evaluators for data of 81 patients. The number of patients required for the study was calculated using the outcome of a pilot study. The results indicate a sensitivity of (91 + or - 3)% and a specificity of (96 + or - 2)% for detecting an overrange, a reduced range is recognized with a sensitivity of (92 + or - 3)% and a specificity of (96 + or - 2)%. The positive and the negative predictive value of this method are 94% and 87%, respectively. The interobserver coefficient of variation is between 3 and 8%. The in-beam PET method demonstrated a high sensitivity and specificity for the detection of range deviations. As the range is a most indicative factor of deviations in the dose delivery, the promising results shown in this paper confirm the in-beam PET method as an appropriate tool for monitoring ion therapy.

In-beam PET measurements of biological half-lives of 12C irradiation induced beta+-activity.

One of the long-standing problems in carbon-ion therapy is the monitoring of the treatment, i.e. of the delivered dose to a given tissue volume within the patient. Over the last 8 years, in-beam positron emission tomography (PET) has been used at the experimental carbon ion treatment facility at the Gesellschaft fur Schwerionenforschung (GSI) Darmstadt and has become a valuable quality assurance tool. In order to determine and evaluate the correct delivery of the patient dose, a simulation of the positron emitter distribution has been compared to the measurement. One particular effect is the blurring as well as the reduction of the measured activity distribution via washout. The objective of this study is the investigation of tissue dependent effective half-lives from patient data. We find no significant dependence of the effective half-life on the Hounsfield unit but on the local dose. The biological half-life within the high dose region is longer than in the low dose region. Furthermore, the influence of the overall treatment time on the kinetics of the positron emitter is reported. There are indications for a metabolic response of the tissue on the irradiation. Taking into account the biological half-life in the simulation leads to an improvement of the quality of the PET-images in some cases.

Direct time-of-flight for quantitative, real-time in-beam PET: a concept and feasibility study.

We extrapolate the impact of recent detector and scintillator developments, enabling sub-nanosecond coincidence timing resolution (tau), onto in-beam positron emission tomography (in-beam PET) for monitoring charged-hadron radiation therapy. For tau < or = 200 ps full width at half maximum, the information given by the time-of-flight (TOF) difference between the two opposing gamma-rays enables shift-variant, artefact-free in-beam tomographic imaging by means of limited-angle, dual-head detectors. We present the corresponding fast, TOF-based and backprojection-free, 3D reconstruction algorithm that, coupled with a real-time data acquisition and a fast detector encoding scheme, allows the sampled beta+-activity to be visualized in the object during the course of the irradiation. Despite the very low statistics scenario typical of in-beam PET, real-treatment simulations show that in-beam TOF-PET enables high-precision images to be obtained in real-time, either with closed-ring or with fixed, dual-head in-beam TOF-PET systems. The latter greatly alleviates the installation of in-beam PET at radiotherapeutic sites.