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

Monte Carlo methods for medical imaging research

In radiation-based medical imaging research, computational modeling methods are used to design and validate imaging systems and post-processing algorithms. Monte Carlo methods are widely used for the computational modeling as they can model the systems accurately and intuitively by sampling interactions between particles and imaging subject with known probability distributions. This article reviews the physics behind Monte Carlo methods, their applications in medical imaging, and available MC codes for medical imaging research. Additionally, potential research areas related to Monte Carlo for medical imaging are discussed.

Evaluating an analytical prediction algorithm of positron emitter distributions in patient data for PET monitoring of carbon ion therapy: A simulation study

In vivo treatment monitoring in ion therapy is one of the key issues for improving the treatment quality assurance procedures. Range verification is one of the most relevant and yet complex task used for in vivo treatment monitoring. In carbon ion therapy, positron emission tomography is the most widely used method. This technique exploits the β+-activity of positron emitters created by nuclear interactions between the incoming beam and the irradiated tissue. Currently, high computational efforts and time-consuming Monte Carlo simulation platforms are typically used to predict positron emitter distributions. In order to avoid time-consuming simulations, an extended filtering approach was suggested to analytically predict positron emitter profiles from depth dose distributions in carbon ion therapy. The purpose of this work is to investigate such an analytical prediction model in patient anatomies of varying complexity, highlighting its potential and the need of further improvements, especially in highly heterogeneous anatomies where many air cavities are present in the beam path. The accuracy of range verification showed a mean relative error of ∼3% and a deviation between the simulation and the prediction below 2mm for the three patient cases analysed: a brain case and two head and neck cases. Additional investigations demonstrated the region of applicability of the method for cases of patient data. The analytical method enables range verification in carbon ion therapy by replacing computing-intensive Monte Carlo simulations and thus minimize the PET monitoring burden on the clinical workflow.

First experimental verification of prompt gamma imaging with carbon ion irradiation

Prompt Gamma Imaging (PGI) is a promising technique for range verification in Particle Therapy. This technique was already tested in clinical environment with a knife-edge-collimator camera for proton treatments but remains relatively unexplored for Carbon Ion Radiation Therapy (CIRT). Previous FLUKA simulations suggested that PG profile shifts could be detected in CIRT with a precision of ∼ 4 mm ([Formula: see text]) for a particle statistic equal to [Formula: see text] C-ions using a 10 × 10 cm2 camera. An experimental campaign was carried out at CNAO (Pavia, Italy) to verify these results, using a knife-edge-collimator camera prototype based on a 5 × 5 cm2 pixelated LYSO crystal. PG profiles were measured irradiating a plastic phantom with a C-ion pencil beam at clinical energies and intensities, also moving the detector to extend the FOV to 13 × 5 cm2. The prototype detected Bragg-peak shifts with ∼ 4 mm precision for a statistic of [Formula: see text] C-ions ([Formula: see text] for the extended FOV), slightly larger than expected. Nevertheless, the detector demonstrated significant potential for verifying the precision in dose delivery following a treatment fraction, which remains fundamental in the clinical environment. For the first time to our knowledge, range verification based on PGI was applied to a C-ion beam at clinical energy and intensities.

Feasibility of the J-PET to monitor the range of therapeutic proton beams

PURPOSE

The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J-PET) scanner for intra-treatment proton beam range monitoring.

METHODS

The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J-PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread-Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J-PET scanner prototype dedicated to the proton beam range assessment.

RESULTS

The investigations indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple-layer dual-head geometry. The results indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for the clinical application, CONCLUSIONS:: We performed simulation studies demonstrating that the feasibility of the J-PET detector for PET-based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre-clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double-layer cylindrical and triple-layer dual-head J-PET geometry configurations seem promising for future clinical application.

Experience and new prospects of PET imaging for ion beam therapy monitoring

Pioneering investigations on the usage of positron-emission-tomography (PET) for the monitoring of ion beam therapy with light (protons, helium) and heavier (stable and radioactive neon, carbon and oxygen) ions started shortly after the first realization of planar and tomographic imaging systems, which were able to visualize the annihilation of positrons resulting from irradiation induced or implanted positron emitting nuclei. And while the first clinical experience was challenged by the utilization of instrumentation directly adapted from nuclear medicine applications, new detectors optimized for this unconventional application of PET imaging are currently entering the phase of (pre)clinical testing for more reliable monitoring of treatment delivery during irradiation. Moreover, recent advances in detector technologies and beam production open several new exciting opportunities which will not only improve the performance of PET imaging under the challenging conditions of in-beam applications in ion beam therapy, but will also likely expand its field of application. In particular, the combination of PET and Compton imaging can enable the most efficient utilization of all possible radiative emissions for both stable and radioactive ion beams, while positronium lifetime imaging may enable probing new features of the underlying tumour and normal tissue environment. Thereby, PET imaging will not only provide means for volumetric reconstruction of the delivered treatment and in-vivo verification of the beam range, but can also shed new insights for biological optimization of the treatment or treatment response assessment.

An inception network for positron emission tomography based dose estimation in carbon ion therapy

Objective. We aim to evaluate a method for estimating 1D physical dose deposition profiles in carbon ion therapy via analysis of dynamic PET images using a deep residual learning convolutional neural network (CNN). The method is validated using Monte Carlo simulations of 12C ion spread-out Bragg peak (SOBP) profiles, and demonstrated with an experimental PET image. Approach. A set of dose deposition and positron annihilation profiles for monoenergetic 12C ion pencil beams in PMMA are first generated using Monte Carlo simulations. From these, a set of random polyenergetic dose and positron annihilation profiles are synthesised and used to train the CNN. Performance is evaluated by generating a second set of simulated 12C ion SOBP profiles (one 116 mm SOBP profile and ten 60 mm SOBP profiles), and using the trained neural network to estimate the dose profile deposited by each beam and the position of the distal edge of the SOBP. Next, the same methods are used to evaluate the network using an experimental PET image, obtained after irradiating a PMMA phantom with a 12C ion beam at QST’s Heavy Ion Medical Accelerator in Chiba facility in Chiba, Japan. The performance of the CNN is compared to that of a recently published iterative technique using the same simulated and experimental 12C SOBP profiles. Main results. The CNN estimated the simulated dose profiles with a mean relative error (MRE) of 0.7% ± 1.0% and the distal edge position with an accuracy of 0.1 mm ± 0.2 mm, and estimate the dose delivered by the experimental 12C ion beam with a MRE of 3.7%, and the distal edge with an accuracy of 1.7 mm. Significance. The CNN was able to produce estimates of the dose distribution with comparable or improved accuracy and computational efficiency compared to the iterative method and other similar PET-based direct dose quantification techniques.

ProTheRaMon - a GATE simulation framework for proton therapy range monitoring using PET imaging

OBJECTIVE

This paper reports on the implementation and shows examples of the use of the ProTheRaMon framework for simulating the delivery of proton therapy treatment plans and range monitoring using positron emission tomography (PET). ProTheRaMon offers complete processing of proton therapy treatment plans, patient CT geometries, and intra-treatment PET imaging, taking into account therapy and imaging coordinate systems and activity decay during the PET imaging protocol specific to a given proton therapy facility. We present the ProTheRaMon framework and illustrate its potential use case and data processing steps for a patient treated at the Cyclotron Centre Bronowice (CCB) proton therapy center in Krakow, Poland.

APPROACH

The ProTheRaMon framework is based on GATE Monte Carlo software, the CASToR reconstruction package and in-house developed Python and bash scripts. The framework consists of five separated simulation and data processing steps, that can be further optimized according to the user's needs and specific settings of a given proton therapy facility and PET scanner design.

MAIN RESULTS

ProTheRaMon is presented using example data from a patient treated at CCB and the J-PET scanner to demonstrate the application of the framework for proton therapy range monitoring. The output of each simulation and data processing stage is described and visualized.

SIGNIFICANCE

We demonstrate that the ProTheRaMon simulation platform is a high-performance tool, capable of running on a computational cluster and suitable for multi-parameter studies, with databases consisting of large number of patients, as well as different PET scanner geometries and settings for range monitoring in a clinical environment. Due to its modular structure, the ProTheRaMon framework can be adjusted for different proton therapy centers and/or different PET detector geometries. It is available to the community via github.

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) based in vivo proton 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-beam during and after beam delivery with a double-head prototype PET camera and (ii) offline shortly 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 radionuclides 11C, 15O, 13N and (in-beam only) 10C. In-beam data 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 the in-beam camera were comparable to or higher compared to respective offline results. Depth profiles of radionuclide-specific yields obtained from the double-head camera showed qualitative differences to data acquired with the full-ring camera 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.

Measurement of biological washout rates depending on tumor vascular status in 15O in-beam rat-PET

Objective. The biological washout of positron emitters should be modeled and corrected in order to achieve quantitative dose range verification in charged particle therapy based on positron emission tomography (PET). This biological washout effect is affected by physiological environmental conditions such as blood perfusion and metabolism, but the correlation to tumour pathology has not been studied yet. Approach. The aim of this study was to investigate the dependence of the biological washout rate on tumour vascular status in rat irradiation. Two types of tumour vascularity conditions, perfused and hypoxic, were modelled with nude rats. The rats were irradiated by a radioactive 15O ion beam and time activity curves were acquired by dynamic in-beam PET measurement. Tumour tissue sections were obtained to observe the histology as well. The biological washout rate was derived using a single-compartment model with two decay components (medium decay, k 2m and slow decay, k 2s ). Main results. All k 2m values in the vascular perfused tumour tissue were higher than the values of the normal tissue. All k 2m values in the hypoxic tumour tissue were much lower than the values of the vascular perfused tumour tissue and slightly lower than the values of the normal tissue. Significance. The dependency of the biological washout on the tumour vasculature conditions was experimentally shown.

Technical Design Report for a Carbon-11 Treatment Facility

Particle therapy relies on the advantageous dose deposition which permits to highly conform the dose to the target and better spare the surrounding healthy tissues and organs at risk with respect to conventional radiotherapy. In the case of treatments with heavier ions (like carbon ions already clinically used), another advantage is the enhanced radiobiological effectiveness due to high linear energy transfer radiation. These particle therapy advantages are unfortunately not thoroughly exploited due to particle range uncertainties. The possibility to monitor the compliance between the ongoing and prescribed dose distribution is a crucial step toward new optimizations in treatment planning and adaptive therapy. The Positron Emission Tomography (PET) is an established quantitative 3D imaging technique for particle treatment verification and, among the isotopes used for PET imaging, the 11C has gained more attention from the scientific and clinical communities for its application as new radioactive projectile for particle therapy. This is an interesting option clinically because of an enhanced imaging potential, without dosimetry drawbacks; technically, because the stable isotope 12C is successfully already in use in clinics. The MEDICIS-Promed network led an initiative to study the possible technical solutions for the implementation of 11C radioisotopes in an accelerator-based particle therapy center. We present here the result of this study, consisting in a Technical Design Report for a 11C Treatment Facility. The clinical usefulness is reviewed based on existing experimental data, complemented by Monte Carlo simulations using the FLUKA code. The technical analysis starts from reviewing the layout and results of the facilities which produced 11C beams in the past, for testing purposes. It then focuses on the elaboration of the feasible upgrades of an existing 12C particle therapy center, to accommodate the production of 11C beams for therapy. The analysis covers the options to produce the 11C atoms in sufficient amounts (as required for therapy), to ionize them as required by the existing accelerator layouts, to accelerate and transport them to the irradiation rooms. The results of the analysis and the identified challenges define the possible implementation scenario and timeline.

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

Can a ToF-PET photon attenuation reconstruction test stopping-power estimations in proton therapy? A phantom study

Objective. The aim of the phantom study was to validate and to improve the computed tomography (CT) images used for the dose computation in proton therapy. It was tested, if the joint reconstruction of activity and attenuation images of time-of-flight PET (ToF-PET) scans could improve the estimation of the proton stopping-power.Approach. The attenuation images, i.e. CT images with 511 keV gamma-rays (γCTs), were jointly reconstructed with activity maps from ToF-PET scans. Theβ⁺activity was produced with FDG and in a separate experiment with proton-induced radioactivation. The phantoms contained slabs of tissue substitutes. The use of theγCTs for the prediction of the beam stopping in proton therapy was based on a linear relationship between theγ-ray attenuation, the electron density, and the stopping-power of fast protons.Main results. The FDG based experiment showed sufficient linearity to detect a bias of bony tissue in the heuristic look-up table, which maps between x-ray CT images and proton stopping-power.γCTs can be used for dose computation, if the electron density of one type of tissue is provided as a scaling factor. A possible limitation is imposed by the spatial resolution, which is inferior by a factor of 2.5 compared to the one of the x-ray CT.γCTs can also be derived from off-line, ToF-PET scans subsequent to the application of a proton field with a hypofractionated dose level.Significance. γCTs are a viable tool to support the estimation of proton stopping with radiotracer-based ToF-PET data from diagnosis or staging. This could be of higher potential relevance in MRI-guided proton therapy.γCTs could form an alternative approach to make use of in-beam or off-line PET scans of proton-inducedβ⁺activity with possible clinical limitations due to the low number of coincidence counts.

Roadmap toward the 10 ps time-of-flight PET challenge

Since the seventies, positron emission tomography (PET) has become an invaluable medical molecular imaging modality with an unprecedented sensitivity at the picomolar level, especially for cancer diagnosis and the monitoring of its response to therapy. More recently, its combination with x-ray computed tomography (CT) or magnetic resonance (MR) has added high precision anatomic information in fused PET/CT and PET/MR images, thus compensating for the modest intrinsic spatial resolution of PET. Nevertheless, a number of medical challenges call for further improvements in PET sensitivity. These concern in particular new treatment opportunities in the context personalized (also called precision) medicine, such as the need to dynamically track a small number of cells in cancer immunotherapy or stem cells for tissue repair procedures. A better signal-to-noise ratio (SNR) in the image would allow detecting smaller size tumours together with a better staging of the patients, thus increasing the chances of putting cancer in complete remission. Moreover, there is an increasing demand for reducing the radioactive doses injected to the patients without impairing image quality. There are three ways to improve PET scanner sensitivity: improving detector efficiency, increasing geometrical acceptance of the imaging device and pushing the timing performance of the detectors. Currently, some pre-localization of the electron-positron annihilation along a line-of-response (LOR) given by the detection of a pair of annihilation photons is provided by the detection of the time difference between the two photons, also known as the time-of-flight (TOF) difference of the photons, whose accuracy is given by the coincidence time resolution (CTR). A CTR of about 10 picoseconds FWHM will ultimately allow to obtain a direct 3D volume representation of the activity distribution of a positron emitting radiopharmaceutical, at the millimetre level, thus introducing a quantum leap in PET imaging and quantification and fostering more frequent use of 11C radiopharmaceuticals. The present roadmap article toward the advent of 10 ps TOF-PET addresses the status and current/future challenges along the development of TOF-PET with the objective to reach this mythic 10 ps frontier that will open the door to real-time volume imaging virtually without tomographic inversion. The medical impact and prospects to achieve this technological revolution from the detection and image reconstruction point-of-views, together with a few perspectives beyond the TOF-PET application are discussed.

Inter-patient variations of radiation-induced normal-tissue changes in Gd-EOB-DTPA-enhanced hepatic MRI scans during fractionated proton therapy

•

Radiation-induced effects visible in Gd-EOB-DTPA enhanced MRI during proton therapy.

•

High inter-patient variation in early MRI signal change during therapy.

•

Correlation of signal change with pretreatment IL-6 concentration.

Background and purpose

Previous MRI studies have shown a substantial decrease in normal-tissue uptake of a hepatobiliary-directed contrast agent 6–9 weeks after liver irradiation. In this prospective clinical study, we investigated whether this effect is detectable during the course of proton therapy.

Material and methods

Gd-EOB-DTPA enhanced MRI was performed twice during hypo-fractionated proton therapy of liver lesions in 9 patients (plus two patients with only one scan available). Dose-correlated signal changes were qualitatively scored based on difference images from the two scans. We evaluated the correlation between the MRI signal change with the planned dose map. The GTV was excluded from all analyses. In addition, were examined timing, irradiated liver volume, changes in liver function parameters as well as circulating biomarkers of inflammation.

Results

Strong, moderate or no dose-correlated signal changes were detected for 2, 3 and 5 patients, respectively. Qualitative scoring was consistent with the quantitative dose to signal change correlation. In an exploratory analysis, the strongest correlation was found between the qualitative scoring and pretreatment IL-6 concentration. For all patients, a clear dose-correlated signal decrease was seen in late follow-up scans.

Conclusion

Radiation-induced effects can be detected with Gd-EOB-DTPA enhanced MRI in a subgroup of patients within a few days after proton irradiation. The reason for the large inter-patient variations is not yet understood and will require validation in larger studies.

Prediction of positron emitter distributions for range monitoring in carbon ion therapy: an analytical approach

Range verification is one of the most relevant tasks in ion beam therapy. In the case of carbon ion therapy, positron emission tomography (PET) is the most widely used method for this purpose, which images the [Formula: see text]-activation following nuclear interactions of the ions with the tissue nuclei. Since the positron emitter activity profile is not directly proportional to the dose distribution, until today only its comparison to a prediction of the PET profile allows for treatment verification. Usually, this prediction is obtained from time-consuming Monte Carlo simulations of high computational effort, which impacts the clinical workflow. To solve this issue in proton therapy, a convolution approach was suggested to predict positron emitter activity profiles from depth dose distributions analytically. In this work, we introduce an approach to predict positron emitter distributions from depth dose profiles in carbon ion therapy. While the distal fall-off position of the positron emitter profile is predicted from a convolution approach similar to the one suggested for protons, additional analytical functions are introduced to describe the characteristics of the positron emitter distribution in tissue. The feasibility of this approach is demonstrated with monoenergetic depth dose profiles and spread out Bragg peaks in homogeneous and heterogeneous phantoms. In all cases, the positron emitter profile is predicted with high precision and the distal fall-off position is reproduced with millimeter accuracy.

Particle beam microstructure reconstruction and coincidence discrimination in PET monitoring for hadron therapy

Positron emission tomography is one of the most mature techniques for monitoring the particles range in hadron therapy, aiming to reduce treatment uncertainties and therefore the extent of safety margins in the treatment plan. In-beam PET monitoring has been already performed using inter-spill and post-irradiation data, i.e. while the particle beam is off or paused. The full beam acquisition procedure is commonly discarded because the particle spills abruptly increase the random coincidence rates and therefore the image noise. This is because random coincidences cannot be separated by annihilation photons originating from radioactive decays and cannot be corrected with standard random coincidence techniques due to the time correlation of the beam-induced background with the ion beam microstructure. The aim of this paper is to provide a new method to recover in-spill data to improve the images obtained with full-beam PET acquisitions. This is done by estimating the temporal microstructure of the beam and thus selecting input PET events that are less likely to be random ones. The PET detector we used was the one developed within the INSIDE project and tested at the CNAO synchrotron-based facility. The data were taken on a PMMA phantom irradiated with 72 MeV proton pencil beams. The obtained results confirm the possibility of improving the acquired PET data without any external signal coming from the synchrotron or ad hoc detectors.

An overview of recent developments in FLUKA PET tools

The new developments of the FLUKA Positron-Emission-Tomography (PET) tools are detailed. FLUKA is a fully integrated Monte Carlo (MC) particle transport code, used for an extended range of applications, including Medical Physics. Recently, it provided the medical community with dedicated simulation tools for clinical applications, including the PET simulation package. PET is a well-established imaging technique in nuclear medicine, and a promising method for clinical in vivo treatment verification in hadrontherapy. The application of clinically established PET scanners to new irradiation environments such as hadrontherapy requires further experimental and theoretical research to which MC simulations could be applied. The FLUKA PET tools, besides featuring PET scanner models in its library, allow the configuration of new PET prototypes via the FLUKA Graphical User Interface (GUI) Flair. Both the beam time structure and scan time can be specified by the user, reproducing PET acquisitions in time, in a particle therapy scenario. Furthermore, different scoring routines allow the analysis of single and coincident events, and identification of parent isotopes generating annihilation events. Two reconstruction codes are currently supported: the Filtered Back-Projection (FBP) and Maximum-Likelihood Expectation Maximization (MLEM), the latter embedded in the tools. Compatibility with other reconstruction frameworks is also possible. The FLUKA PET tools package has been successfully tested for different detectors and scenarios, including conventional functional PET applications and in beam PET, either using radioactive sources, or simulating hadron beam irradiations. The results obtained so far confirm the FLUKA PET tools suitability to perform PET simulations in R&D environment.

Carbon ions beam therapy monitoring with the INSIDE in-beam PET

In vivo range monitoring techniques are necessary in order to fully take advantage of the high dose gradients deliverable in hadrontherapy treatments. Positron emission tomography (PET) scanners can be used to monitor beam-induced activation in tissues and hence measure the range. The INSIDE (Innovative Solutions for In-beam DosimEtry in Hadrontherapy) in-beam PET scanner, installed at the Italian National Center of Oncological Hadrontherapy (CNAO, Pavia, Italy) synchrotron facility, has already been successfully tested in vivo during a proton therapy treatment. We discuss here the system performance evaluation with carbon ion beams, in view of future in vivo tests. The work is focused on the analysis of activity images obtained with therapeutic treatments delivered to polymethyl methacrylate (PMMA) phantoms, as well as on the test of an innovative and robust Monte Carlo simulation technique for the production of reliable prior activity maps. Images are reconstructed using different integration intervals, so as to monitor the activity evolution during and after the treatment. Three procedures to compare activity images are presented, namely Pearson correlation coefficient, Beam's eye view and overall view. Images of repeated irradiations of the same treatments are compared to assess the integration time necessary to provide reproducible images. The range agreement between simulated and experimental images is also evaluated, so as to validate the simulation capability to provide sound prior information. The results indicate that at treatment end, or at most 20 s afterwards, the range measurement is reliable within 1-2 mm, when comparing both different experimental sessions and data with simulations. In conclusion, this work shows that the INSIDE in-beam PET scanner performance is promising towards its in vivo test with carbon ions.

Improving the modelling of irradiation-induced brain activation for in vivo PET verification of proton therapy

BACKGROUND AND PURPOSES

A reliable Monte Carlo prediction of proton-induced brain tissue activation used for comparison to particle therapy positron-emission-tomography (PT-PET) measurements is crucial for in vivo treatment verification. Major limitations of current approaches to overcome include the CT-based patient model and the description of activity washout due to tissue perfusion.

MATERIAL AND METHODS

Two approaches were studied to improve the activity prediction for brain irradiation: (i) a refined patient model using tissue classification based on MR information and (ii) a PT-PET data-driven refinement of washout model parameters. Improvements of the activity predictions compared to post-treatment PT-PET measurements were assessed in terms of activity profile similarity for six patients treated with a single or two almost parallel fields delivered by active proton beam scanning.

RESULTS

The refined patient model yields a generally higher similarity for most of the patients, except in highly pathological areas leading to tissue misclassification. Using washout model parameters deduced from clinical patient data could considerably improve the activity profile similarity for all patients.

CONCLUSIONS

Current methods used to predict proton-induced brain tissue activation can be improved with MR-based tissue classification and data-driven washout parameters, thus providing a more reliable basis for PT-PET verification.

Online proton therapy monitoring: clinical test of a Silicon-photodetector-based in-beam PET

Particle therapy exploits the energy deposition pattern of hadron beams. The narrow Bragg Peak at the end of range is a major advantage but range uncertainties can cause severe damage and require online verification to maximise the effectiveness in clinics. In-beam Positron Emission Tomography (PET) is a non-invasive, promising in-vivo technique, which consists in the measurement of the β+ activity induced by beam-tissue interactions during treatment, and presents the highest correlation of the measured activity distribution with the deposited dose, since it is not much influenced by biological washout. Here we report the first clinical results obtained with a state-of-the-art in-beam PET scanner, with on-the-fly reconstruction of the activity distribution during irradiation. An automated time-resolved quantitative analysis was tested on a lacrimal gland carcinoma case, monitored during two consecutive treatment sessions. The 3D activity map was reconstructed every 10 s, with an average delay between beam delivery and image availability of about 6 s. The correlation coefficient of 3D activity maps for the two sessions (above 0.9 after 120 s) and the range agreement (within 1 mm) prove the suitability of in-beam PET for online range verification during treatment, a crucial step towards adaptive strategies in particle therapy.

Three dimensional reconstruction of therapeutic carbon ion beams in phantoms using single secondary ion tracks

Carbon ion beam radiotherapy enables a very localised dose deposition. However, even small changes in the patient geometry or positioning errors can significantly distort the dose distribution. A live, non-invasive monitoring system of the beam delivery within the patient is therefore highly desirable, and could improve patient treatment. We present a novel three-dimensional method for imaging the beam in the irradiated object, exploiting the measured tracks of single secondary ions emerging under irradiation. The secondary particle tracks are detected with a TimePix stack-a set of parallel pixelated semiconductor detectors. We developed a three-dimensional reconstruction algorithm based on maximum likelihood expectation maximization. We demonstrate the applicability of the new method in the irradiation of a cylindrical PMMA phantom of human head size with a carbon ion pencil beam of [Formula: see text] MeV u^-1. The beam image in the phantom is reconstructed from a set of nine discrete detector positions between [Formula: see text] and [Formula: see text] from the beam axis. Furthermore, we demonstrate the potential to visualize inhomogeneities by irradiating a PMMA phantom with an air gap as well as bone and adipose tissue surrogate inserts. We successfully reconstructed a three-dimensional image of the treatment beam in the phantom from single secondary ion tracks. The beam image corresponds well to the beam direction and energy. In addition, cylindrical inhomogeneities with a diameter of [Formula: see text] cm and density differences down to [Formula: see text] g cm^-3 to the surrounding material are clearly visualized. This novel three-dimensional method to image a therapeutic carbon ion beam in the irradiated object does not interfere with the treatment and requires knowledge only of single secondary ion tracks. Even with detectors with only a small angular coverage, the three-dimensional reconstruction of the fragmentation points presented in this work was found to be feasible.

Visualization of air and metal inhomogeneities in phantoms irradiated by carbon ion beams using prompt secondary ions

PURPOSE

Non-invasive methods for monitoring of the therapeutic ion beam extension in the patient are desired in order to handle deteriorations of the dose distribution related to changes of the patient geometry. In carbon ion radiotherapy, secondary light ions represent one of potential sources of information about the dose distribution in the irradiated target. The capability to detect range-changing inhomogeneities inside of an otherwise homogeneous phantom, based on single track measurements, is addressed in this paper.

METHODS

Air and stainless steel inhomogeneities, with PMMA equivalent thickness of 10mm and 4.8mm respectively, were inserted into a PMMA-phantom at different positions in depth. Irradiations of the phantom with therapeutic carbon ion pencil beams were performed at the Heidelberg Ion Beam Therapy Center. Tracks of single secondary ions escaping the phantom under irradiation were detected with a pixelized semiconductor detector Timepix. The statistical relevance of the found differences between the track distributions with and without inhomogeneities was evaluated.

RESULTS

Measured shifts of the distal edge and changes in the fragmentation probability make the presence of inhomogeneities inserted into the traversed medium detectable for both, 10mm air cavities and 1mm thick stainless steel. Moreover, the method was shown to be sensitive also on their position in the observed body, even when localized behind the Bragg-peak.

CONCLUSIONS

The presented results demonstrate experimentally, that the method using distributions of single secondary ion tracks is sensitive to the changes of homogeneity of the traversed material for the studied geometries of the target.

Feasibility study of using fall-off gradients of early and late PET scans for proton range verification

PURPOSE

While positron emission tomography (PET) allows for the imaging of tissues activated by proton beams in terms of monitoring the therapy administered, most endogenous tissue elements are activated by relatively high-energy protons. Therefore, a relatively large distance off-set exists between the dose fall-off and activity fall-off. However, ¹⁶ O(p,2p,2n)¹³ N has a relatively low energy threshold which peaks around 12 MeV and also a residual proton range that is approximately 1 to 2 mm. In this phantom study, we tested the feasibility of utilizing the ¹³ N production peak as well as the differences in activity fall-off between early and late PET scans for proton range verification. One of the main purposes for this research was developing a proton range verification methodology that would not require Monte Carlo simulations.

METHODS AND MATERIALS

Both monoenergetic and spread-out Bragg peak beams were delivered to two phantoms - a water-like gel and a tissue-like gel where the proton ranges came to be approximately 9.9 and 9.1 cm, respectively. After 1 min of postirradiation delay, the phantoms were scanned for a period of 30 min using an in-room PET. Two separate (Early and Late) PET images were reconstructed using two different postirradiation delays and acquisition times; Early PET: 1 min delay and 3 min acquisition, Late PET: 21 min delay and 10 min acquisition. The depth gradients of the PET signals were then normalized and plotted as functions of depth. The normalized gradient of the early PET images was subtracted from that of the late PET images, to observe the ¹³ N activity distribution in relation to depth. Monte Carlo simulations were also conducted with the same set-up as the measurements stated previously.

RESULTS

The subtracted gradients show peaks at 9.4 and 8.6 cm in water-gel and tissue-gel respectively for both pristine and SOBP beams. These peaks are created in connection with the sudden change of ¹³ N signals with depth and consistently occur 2 mm upstream to where ¹³ N signals were most abundantly created (9.6 and 8.8 cm in water-gel and tissue-gel, respectively). Monte Carlo simulations provided similar results as the measurements.

CONCLUSIONS

The subtracted PET signal gradient peaks and the proton ranges for water-gel and tissue-gel show distance off-sets of 4 to 5 mm. This off-set may potentially be used for proton range verification using only the PET measured data without Monte Carlo simulations. More studies are necessary to overcome various limitations, such as perfusion-driven washout, for the feasibility of this technique in living patients.

Full-beam performances of a PET detector with synchrotron therapeutic proton beams

Treatment quality assessment is a crucial feature for both present and next-generation ion therapy facilities. Several approaches are being explored, based on prompt radiation emission or on PET signals by [Formula: see text]-decaying isotopes generated by beam interactions with the body. In-beam PET monitoring at synchrotron-based ion therapy facilities has already been performed, either based on inter-spill data only, to avoid the influence of the prompt radiation, or including both in-spill and inter-spill data. However, the PET images either suffer of poor statistics (inter-spill) or are more influenced by the background induced by prompt radiation (in-spill). Both those problems are expected to worsen for accelerators with improved duty cycle where the inter-spill interval is reduced to shorten the treatment time. With the aim of assessing the detector performance and developing techniques for background reduction, a test of an in-beam PET detector prototype was performed at the CNAO synchrotron-based ion therapy facility in full-beam acquisition modality. Data taken with proton beams impinging on PMMA phantoms showed the system acquisition capability and the resulting activity distribution, separately reconstructed for the in-spill and the inter-spill data. The coincidence time resolution for in-spill and inter-spill data shows a good agreement, with a slight deterioration during the spill. The data selection technique allows the identification and rejection of most of the background originated during the beam delivery. The activity range difference between two different proton beam energies (68 and 72 MeV) was measured and found to be in sub-millimeter agreement with the expected result. However, a slightly longer (2 mm) absolute profile length is obtained for in-spill data when compared to inter-spill data.

INSIDE in-beam positron emission tomography system for particle range monitoring in hadrontherapy

The quality assurance of particle therapy treatment is a fundamental issue that can be addressed by developing reliable monitoring techniques and indicators of the treatment plan correctness. Among the available imaging techniques, positron emission tomography (PET) has long been investigated and then clinically applied to proton and carbon beams. In 2013, the Innovative Solutions for Dosimetry in Hadrontherapy (INSIDE) collaboration proposed an innovative bimodal imaging concept that combines an in-beam PET scanner with a tracking system for charged particle imaging. This paper presents the general architecture of the INSIDE project but focuses on the in-beam PET scanner that has been designed to reconstruct the particles range with millimetric resolution within a fraction of the dose delivered in a treatment of head and neck tumors. The in-beam PET scanner has been recently installed at the Italian National Center of Oncologic Hadrontherapy (CNAO) in Pavia, Italy, and the commissioning phase has just started. The results of the first beam test with clinical proton beams on phantoms clearly show the capability of the in-beam PET to operate during the irradiation delivery and to reconstruct on-line the beam-induced activity map. The accuracy in the activity distal fall-off determination is millimetric for therapeutic doses.

Initial clinical evaluation of PET-based ion beam therapy monitoring under consideration of organ motion

PURPOSE

Intrafractional organ motion imposes considerable challenges to scanned ion beam therapy and demands for a thorough verification of the applied treatment. At the Heidelberg Ion-Beam Therapy Center (HIT), the scanned ion beam delivery is verified by means of postirradiation positron-emission-tomography (PET) imaging. This work presents a first clinical evaluation of PET-based treatment monitoring in ion beam therapy under consideration of target motion.

METHODS

Three patients with mobile liver lesions underwent scanned carbon ion irradiation at HIT and postirradiation PET/CT (x-ray-computed-tomography) imaging with a commercial scanner. Respiratory motion was recorded during irradiation and subsequent image acquisition. This enabled a time-resolved (4D) calculation of the expected irradiation-induced activity pattern and, for one patient where an additional 4D CT was acquired at the PET/CT scanner after treatment, a motion-compensated PET image reconstruction. For the other patients, PET data were reconstructed statically. To verify the treatment, calculated prediction and reconstructed measurement were compared with a focus on the ion beam range.

RESULTS

Results in the current three patients suggest that for motion amplitudes in the order of 2 mm there is no benefit from incorporating respiratory motion information into PET-based treatment monitoring. For a target motion in the order of 10 mm, motion-related effects become more severe and a time-resolved modeling of the expected activity distribution can lead to an improved data interpretation if a sufficient number of true coincidences is detected. Benefits from motion-compensated PET image reconstruction could not be shown conclusively at the current stage.

CONCLUSIONS

The feasibility of clinical PET-based treatment verification under consideration of organ motion has been shown for the first time. Improvements in noise-robust 4D PET image reconstruction are deemed necessary to enhance the clinical potential.

Nuclear physics in particle therapy: a review

Charged particle therapy has been largely driven and influenced by nuclear physics. The increase in energy deposition density along the ion path in the body allows reducing the dose to normal tissues during radiotherapy compared to photons. Clinical results of particle therapy support the physical rationale for this treatment, but the method remains controversial because of the high cost and of the lack of comparative clinical trials proving the benefit compared to x-rays. Research in applied nuclear physics, including nuclear interactions, dosimetry, image guidance, range verification, novel accelerators and beam delivery technologies, can significantly improve the clinical outcome in particle therapy. Measurements of fragmentation cross-sections, including those for the production of positron-emitting fragments, and attenuation curves are needed for tuning Monte Carlo codes, whose use in clinical environments is rapidly increasing thanks to fast calculation methods. Existing cross sections and codes are indeed not very accurate in the energy and target regions of interest for particle therapy. These measurements are especially urgent for new ions to be used in therapy, such as helium. Furthermore, nuclear physics hardware developments are frequently finding applications in ion therapy due to similar requirements concerning sensors and real-time data processing. In this review we will briefly describe the physics bases, and concentrate on the open issues.

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.

Characterizing proton-activated materials to develop PET-mediated proton range verification markers

Conventional proton beam range verification using positron emission tomography (PET) relies on tissue activation alone and therefore requires particle therapy PET whose installation can represent a large financial burden for many centers. Previously, we showed the feasibility of developing patient implantable markers using high proton cross-section materials ((18)O, Cu, and (68)Zn) for in vivo proton range verification using conventional PET scanners. In this technical note, we characterize those materials to test their usability in more clinically relevant conditions. Two phantoms made of low-density balsa wood (~0.1 g cm(-3)) and beef (~1.0 g cm(-3)) were embedded with Cu or (68)Zn foils of several volumes (10-50 mm(3)). The metal foils were positioned at several depths in the dose fall-off region, which had been determined from our previous study. The phantoms were then irradiated with different proton doses (1-5 Gy). After irradiation, the phantoms with the embedded foils were moved to a diagnostic PET scanner and imaged. The acquired data were reconstructed with 20-40 min of scan time using various delay times (30-150 min) to determine the maximum contrast-to-noise ratio. The resultant PET/computed tomography (CT) fusion images of the activated foils were then examined and the foils' PET signal strength/visibility was scored on a 5 point scale by 13 radiologists experienced in nuclear medicine. For both phantoms, the visibility of activated foils increased in proportion to the foil volume, dose, and PET scan time. A linear model was constructed with visibility scores as the response variable and all other factors (marker material, phantom material, dose, and PET scan time) as covariates. Using the linear model, volumes of foils that provided adequate visibility (score 3) were determined for each dose and PET scan time. The foil volumes that were determined will be used as a guideline in developing practical implantable markers.

The FLUKA Code: An Accurate Simulation Tool for Particle Therapy

Monte Carlo (MC) codes are increasingly spreading in the hadrontherapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code for application to hadrontherapy demands accurate and reliable physical models capable of handling all components of the expected radiation field. This becomes extremely important for correctly performing not only physical but also biologically based dose calculations, especially in cases where ions heavier than protons are involved. In addition, accurate prediction of emerging secondary radiation is of utmost importance in innovative areas of research aiming at in vivo treatment verification. This contribution will address the recent developments of the FLUKA MC code and its practical applications in this field. Refinements of the FLUKA nuclear models in the therapeutic energy interval lead to an improved description of the mixed radiation field as shown in the presented benchmarks against experimental data with both (4)He and (12)C ion beams. Accurate description of ionization energy losses and of particle scattering and interactions lead to the excellent agreement of calculated depth-dose profiles with those measured at leading European hadron therapy centers, both with proton and ion beams. In order to support the application of FLUKA in hospital-based environments, Flair, the FLUKA graphical interface, has been enhanced with the capability of translating CT DICOM images into voxel-based computational phantoms in a fast and well-structured way. The interface is capable of importing also radiotherapy treatment data described in DICOM RT standard. In addition, the interface is equipped with an intuitive PET scanner geometry generator and automatic recording of coincidence events. Clinically, similar cases will be presented both in terms of absorbed dose and biological dose calculations describing the various available features.

Investigating the limits of PET/CT imaging at very low true count rates and high random fractions in ion-beam therapy monitoring

PURPOSE

External beam radiotherapy with protons and heavier ions enables a tighter conformation of the applied dose to arbitrarily shaped tumor volumes with respect to photons, but is more sensitive to uncertainties in the radiotherapeutic treatment chain. Consequently, an independent verification of the applied treatment is highly desirable. For this purpose, the irradiation-induced β(+)-emitter distribution within the patient is detected shortly after irradiation by a commercial full-ring positron emission tomography/x-ray computed tomography (PET/CT) scanner installed next to the treatment rooms at the Heidelberg Ion-Beam Therapy Center (HIT). A major challenge to this approach is posed by the small number of detected coincidences. This contribution aims at characterizing the performance of the used PET/CT device and identifying the best-performing reconstruction algorithm under the particular statistical conditions of PET-based treatment monitoring. Moreover, this study addresses the impact of radiation background from the intrinsically radioactive lutetium-oxyorthosilicate (LSO)-based detectors at low counts.

METHODS

The authors have acquired 30 subsequent PET scans of a cylindrical phantom emulating a patientlike activity pattern and spanning the entire patient counting regime in terms of true coincidences and random fractions (RFs). Accuracy and precision of activity quantification, image noise, and geometrical fidelity of the scanner have been investigated for various reconstruction algorithms and settings in order to identify a practical, well-suited reconstruction scheme for PET-based treatment verification. Truncated listmode data have been utilized for separating the effects of small true count numbers and high RFs on the reconstructed images. A corresponding simulation study enabled extending the results to an even wider range of counting statistics and to additionally investigate the impact of scatter coincidences. Eventually, the recommended reconstruction scheme has been applied to exemplary postirradiation patient data-sets.

RESULTS

Among the investigated reconstruction options, the overall best results in terms of image noise, activity quantification, and accurate geometrical recovery were achieved using the ordered subset expectation maximization reconstruction algorithm with time-of-flight (TOF) and point-spread function (PSF) information. For this algorithm, reasonably accurate (better than 5%) and precise (uncertainty of the mean activity below 10%) imaging can be provided down to 80,000 true coincidences at 96% RF. Image noise and geometrical fidelity are generally improved for fewer iterations. The main limitation for PET-based treatment monitoring has been identified in the small number of true coincidences, rather than the high intrinsic random background. Application of the optimized reconstruction scheme to patient data-sets results in a 25% - 50% reduced image noise at a comparable activity quantification accuracy and an improved geometrical performance with respect to the formerly used reconstruction scheme at HIT, adopted from nuclear medicine applications.

CONCLUSIONS

Under the poor statistical conditions in PET-based treatment monitoring, improved results can be achieved by considering PSF and TOF information during image reconstruction and by applying less iterations than in conventional nuclear medicine imaging. Geometrical fidelity and image noise are mainly limited by the low number of true coincidences, not the high LSO-related random background. The retrieved results might also impact other emerging PET applications at low counting statistics.

Monitoring proton therapy with PET

Protons are being used in radiation therapy because of typically better dose conformity and reduced total energy deposited in the patient as compared with photon techniques. Both aspects are related to the finite range of a proton beam. The finite range also allows advanced dose shaping. These benefits can only be fully utilized if the end of range can be predicted accurately in the patient. The prediction of the range in tissue is associated with considerable uncertainties owing to imaging, patient set-up, beam delivery, interfractional changes in patient anatomy and dose calculation. Consequently, a significant range (of the order of several millimetres) is added to the prescribed range in order to ensure tumour coverage. Thus, reducing range uncertainties would allow a reduction of the treatment volume and reduce dose to potential organs at risk.

Comparing the biological washout of β+-activity induced in mice brain after 12C-ion and proton irradiation

In clinical ion beam therapy, protons as well as heavier ions such as carbon are used for treatment. For protons, β(+)-emitters are only induced by fragmentation reactions in the target (target fragmentation), whereas for heavy ions, they are additionally induced by fragmentations of the projectile (further referred to as autoactivation). An approach utilizing these processes for treatment verification, by comparing measured Positron Emission Tomography (PET) data to predictions from Monte Carlo simulations, has already been clinically implemented. For an accurate simulation, it is important to consider the biological washout of β(+)-emitters due to vital functions. To date, mathematical expressions for washout have mainly been determined by using radioactive beams of (10)C- and (11)C-ions, both β(+)-emitters, to enhance the counting statistics in the irradiated area. Still, the question of how the choice of projectile (autoactivating or non-autoactivating) influences the washout coefficients, has not been addressed. In this context, an experiment was carried out at the Heidelberg Ion Beam Therapy Center with the purpose of directly comparing irradiation-induced biological washout coefficients in mice for protons and (12)C-ions. To this aim, mice were irradiated in the brain region with protons and (12)C-ions and measured after irradiation with a PET/CT scanner (Siemens Biograph mCT). After an appropriate waiting time, the mice were sacrificed, then irradiated and measured again under similar conditions. The resulting data were processed and fitted numerically to deduce the main washout parameters. Despite the very low PET counting statistics, a consistent difference could be identified between (12)C-ion and proton irradiated mice, with the (12)C data being described best by a two component fit with a combined medium and slow washout fraction of 0.50 ± 0.05 and the proton mice data being described best by a one component fit with only one (slow) washout fraction of 0.73 ± 0.06.