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

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.

In-beam PET treatment monitoring of carbon therapy patients: Results of a clinical trial at CNAO

PURPOSE

Carbon ion therapy treatments can be monitored non-invasively with in-beam Positron Emission Tomography (PET). At CNAO the INSIDE in-beam PET scanner has been used in a clinical trial (NCT03662373) to monitor cancer treatments with proton and carbon therapy. In this work we present the analysis results of carbon therapy data, acquired during the first phase of the clinical trial, analyzing data of nine patients treated at CNAO for various malignant tumors in the head-and-neck region.

MATERIALS AND METHODS

The patient group contained two patients requiring replanning, and seven patients without replanning, based on established protocols. For each patient the PET images acquired along the course of treatment were compared with a reference, applying two analysis methods: the beam-eye-view (BEV) method and the γ-index analysis. Time trends in several parameters were investigated, as well as the agreement with control CTs, if available.

RESULTS

Regarding the BEV-method, the average sigma value σ was 3.7 mm of range difference distributions for patients without changes (sensitivity of the INSIDE detector). The 3D-information obtained from the BEV analysis was partly in agreement with what was observed in the control CT. The data quality and quantity was insufficient for a definite interpretation of the time trends.

CONCLUSION

We analyzed carbon therapy data acquired with the INSIDE in-beam PET detector using two analysis methods. The data allowed to evaluate sensitivity of the INSIDE detector for carbon therapy and to make several recommendations for the future.

Investigation of intra-fractionated range guided adaptive proton therapy: I. On-line PET imaging and range measurement

Objective.Develop a prototype on-line positron emission tomography (PET) scanner and evaluate its capability of on-line imaging and intra-fractionated proton-induced radioactivity range measurement.Approach.Each detector consists of 32 × 32 array of 2 × 2 × 30 mm3Lutetium-Yttrium Oxyorthosilicate scintillators with single-scintillator-end readout through a 20 × 20 array of 3 × 3 mm2Silicon Photomultipliers. The PET can be configurated with a full-ring of 20 detectors for conventional PET imaging or a partial-ring of 18 detectors for on-line imaging and range measurement. All detector-level readout and processing electronics are attached to the backside of the system gantry and their output signals are transferred to a field-programable-gate-array based system electronics and data acquisition that can be placed 2 m away from the gantry. The PET imaging performance and radioactivity range measurement capability were evaluated by both the offline study that placed a radioactive source with known intensity and distribution within a phantom and the online study that irradiated a phantom with proton beams under different radiation and imaging conditions.Main results.The PET has 32 cm diameter and 6.5 cm axial length field-of-view (FOV), ∼2.3-5.0 mm spatial resolution within FOV, 3% sensitivity at the FOV center, 18%-30% energy resolution, and ∼9 ns coincidence time resolution. The offline study shows the PET can determine the shift of distal falloff edge position of a known radioactivity distribution with the accuracy of 0.3 ± 0.3 mm even without attenuation and scatter corrections, and online study shows the PET can measure the shift of proton-induced positron radioactive range with the accuracy of 0.6 ± 0.3 mm from the data acquired with a short-acquisition (60 s) and low-dose (5 MU) proton radiation to a human head phantom.Significance.This study demonstrated the capability of intra-fractionated PET imaging and radioactivity range measurement and will enable the investigation on the feasibility of intra-fractionated, range-shift compensated adaptive proton therapy.

An in-vivo treatment monitoring system for ion-beam radiotherapy based on 28 Timepix3 detectors

Ion-beam radiotherapy is an advanced cancer treatment modality offering steep dose gradients and a high biological effectiveness. These gradients make the therapy vulnerable to patient-setup and anatomical changes between treatment fractions, which may go unnoticed. Charged fragments from nuclear interactions of the ion beam with the patient tissue may carry information about the treatment quality. Currently, the fragments escape the patient undetected. Inter-fractional in-vivo treatment monitoring based on these charged nuclear fragments could make ion-beam therapy safer and more efficient. We developed an ion-beam monitoring system based on 28 hybrid silicon pixel detectors (Timepix3) to measure the distribution of fragment origins in three dimensions. The system design choices as well as the ion-beam monitoring performance measurements are presented in this manuscript. A spatial resolution of 4 mm along the beam axis was achieved for the measurement of individual fragment origins. Beam-range shifts of1.5 mm were identified in a clinically realistic treatment scenario with an anthropomorphic head phantom. The monitoring system is currently being used in a prospective clinical trial at the Heidelberg Ion Beam Therapy Centre for head-and-neck as well as central nervous system cancer patients.

Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study

GOAL

In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes.

METHODS/MATERIALS

We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans.

RESULTS

Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields.

CONCLUSION

The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.

Measurement of the12C(p,n)12N reaction cross section below 150 MeV

Objective. Proton therapy currently faces challenges from clinical complications on organs-at-risk due to range uncertainties. To address this issue, positron emission tomography (PET) of the proton-induced11C and15O activity has been used to provide feedback on the proton range. However, this approach is not instantaneous due to the relatively long half-lives of these nuclides. An alternative nuclide,12N (half-life 11 ms), shows promise for real-timein vivoproton range verification. Development of12N imaging requires better knowledge of its production reaction cross section.Approach. The12C(p,n)12N reaction cross section was measured by detecting positron activity of graphite targets irradiated with 66.5, 120, and 150 MeV protons. A pulsed beam delivery with 0.7-2 × 108protons per pulse was used. The positron activity was measured during the beam-off periods using a dual-head Siemens Biograph mCT PET scanner. The12N production was determined from activity time histograms.Main results. The cross section was calculated for 11 energies, ranging from 23.5 to 147 MeV, using information on the experimental setup and beam delivery. Through a comprehensive uncertainty propagation analysis, a statistical uncertainty of 2.6%-5.8% and a systematic uncertainty of 3.3%-4.6% were achieved. Additionally, a comparison between measured and simulated scanner sensitivity showed a scaling factor of 1.25 (±3%). Despite this, there was an improvement in the precision of the cross section measurement compared to values reported by the only previous study.Significance. Short-lived12N imaging is promising for real-timein vivoverification of the proton range to reduce clinical complications in proton therapy. The verification procedure requires experimental knowledge of the12N production cross section for proton energies of clinical importance, to be incorporated in a Monte Carlo framework for12N imaging prediction. This study is the first to achieve a precise measurement of the12C(p,n)12N nuclear cross section for such proton energies.

Synthetic CT imaging for PET monitoring in proton therapy: a simulation study

Objective.This study addresses a fundamental limitation of in-beam positron emission tomography (IB-PET) in proton therapy: the lack of direct anatomical representation in the images it produces. We aim to overcome this shortcoming by pioneering the application of deep learning techniques to create synthetic control CT images (sCT) from combining IB-PET and planning CT scan data.Approach.We conducted simulations involving six patients who underwent irradiation with proton beams. Leveraging the architecture of a visual transformer (ViT) neural network, we developed a model to generate sCT images of these patients using the planning CT scans and the inter-fractional simulated PET activity maps during irradiation. To evaluate the model's performance, a comparison was conducted between the sCT images produced by the ViT model and the authentic control CT images-serving as the benchmark.Main results.The structural similarity index was computed at a mean value across all patients of 0.91, while the mean absolute error measured 22 Hounsfield Units (HU). Root mean squared error and peak signal-to-noise ratio values were 56 HU and 30 dB, respectively. The Dice similarity coefficient exhibited a value of 0.98. These values are comparable to or exceed those found in the literature. More than 70% of the synthetic morphological changes were found to be geometrically compatible with the ones reported in the real control CT scan.Significance.Our study presents an innovative approach to surface the hidden anatomical information of IB-PET in proton therapy. Our ViT-based model successfully generates sCT images from inter-fractional PET data and planning CT scans. Our model's performance stands on par with existing models relying on input from cone beam CT or magnetic resonance imaging, which contain more anatomical information than activity maps.

Few-seconds range verification with short-lived positron emitters in carbon ion therapy

In-beam PET (Positron Emission Tomography) is one of the most precise techniques for in-vivo range monitoring in hadron therapy. Our objective was to demonstrate the feasibility of a short irradiation run for range verification before a carbon-ion treatment. To do so a PMMA target was irradiated with a 220 MeV/u carbon-ion beam and annihilation coincidences from short-lived positron emitters were acquired after irradiations lasting 0.6 s. The experiments were performed at the synchrotron-based facility CNAO (Italian National Center of Oncological Hadrontherapy) by using the INSIDE in-beam PET detector. The results show that, with 3·107 carbon ions, the reconstructed positron emitting nuclei distribution is in good agreement with the predictions of a detailed FLUKA Monte Carlo study. Moreover, the radio-nuclei production is sufficiently abundant to determine the average ion beam range with a σ of 1 mm with a 6 s measurement of the activity distribution. Since the data were acquired when the beam was off, the proposed rapid calibration method can be applied to hadron beams extracted from accelerators with very different time structures.

Hybrid PET/Compton-camera imaging: an imager for the next generation

Compton cameras can offer advantages over gamma cameras for some applications, since they are well suited for multitracer imaging and for imaging high-energy radiotracers, such as those employed in radionuclide therapy. While in conventional clinical settings state-of-the-art Compton cameras cannot compete with well-established methods such as PET and SPECT, there are specific scenarios in which they can constitute an advantageous alternative. The combination of PET and Compton imaging can benefit from the improved resolution and sensitivity of current PET technology and, at the same time, overcome PET limitations in the use of multiple radiotracers. Such a system can provide simultaneous assessment of different radiotracers under identical conditions and reduce errors associated with physical factors that can change between acquisitions. Advances are being made both in instrumentation developments combining PET and Compton cameras for multimodal or three-gamma imaging systems, and in image reconstruction, addressing the challenges imposed by the combination of the two modalities or the new techniques. This review article summarizes the advances made in Compton cameras for medical imaging and their combination with PET.

Technological advancements in cancer diagnostics: Improvements and limitations

BACKGROUND

Cancer is characterized by the rampant proliferation, growth, and infiltration of malignantly transformed cancer cells past their normal boundaries into adjacent tissues. It is the leading cause of death worldwide, responsible for approximately 19.3 million new diagnoses and 10 million deaths globally in 2020. In the United States alone, the estimated number of new diagnoses and deaths is 1.9 million and 609 360, respectively. Implementation of currently existing cancer diagnostic techniques such as positron emission tomography (PET), X-ray computed tomography (CT), and magnetic resonance spectroscopy (MRS), and molecular diagnostic techniques, have enabled early detection rates and are instrumental not only for the therapeutic management of cancer patients, but also for early detection of the cancer itself. The effectiveness of these cancer screening programs are heavily dependent on the rate of accurate precursor lesion identification; an increased rate of identification allows for earlier onset treatment, thus decreasing the incidence of invasive cancer in the long-term, and improving the overall prognosis. Although these diagnostic techniques are advantageous due to lack of invasiveness and easier accessibility within the clinical setting, several limitations such as optimal target definition, high signal to background ratio and associated artifacts hinder the accurate diagnosis of specific types of deep-seated tumors, besides associated high cost. In this review we discuss various imaging, molecular, and low-cost diagnostic tools and related technological advancements, to provide a better understanding of cancer diagnostics, unraveling new opportunities for effective management of cancer, particularly in low- and middle-income countries (LMICs).

RECENT FINDINGS

Herein we discuss various technological advancements that are being utilized to construct an assortment of new diagnostic techniques that incorporate hardware, image reconstruction software, imaging devices, biomarkers, and even artificial intelligence algorithms, thereby providing a reliable diagnosis and analysis of the tumor. Also, we provide a brief account of alternative low cost-effective cancer therapy devices (CryoPop®, LumaGEM®, MarginProbe®) and picture archiving and communication systems (PACS), emphasizing the need for multi-disciplinary collaboration among radiologists, pathologists, and other involved specialties for improving cancer diagnostics.

CONCLUSION

Revolutionary technological advancements in cancer imaging and molecular biology techniques are indispensable for the accurate diagnosis and prognosis of cancer.

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.

Depth dose measurements in water for 11C and 10C beams with therapy relevant energies

Owing to the favorable depth-dose distribution and the radiobiological properties of heavy ion radiation, ion beam therapy shows an improved success/toxicity ratio compared to conventional radiotherapy. The sharp dose gradients and very high doses in the Bragg peak region, which represent the larger physical advantage of ion beam therapy, make it also extremely sensitive to range uncertainties. The use of β ⁺-radioactive ion beams would be ideal for simultaneous treatment and accurate online range monitoring through PET imaging. Since all the unfragmented primary ions are potentially contributing to the PET signal, these beams offer an improved image quality while preserving the physical and radiobiological advantages of the stable counterparts. The challenging production of radioactive ion beams and the difficulties in reaching high intensities, have discouraged their clinical application. In this context, the project Biomedical Applications of Radioactive ion Beams (BARB) started at GSI (Helmholtzzentrum für Schwerionenforschung GmbH) with the main goal to assess the technical feasibility and investigate possible advantages of radioactive ion beams on the pre-clinical level. During the first experimental campaign ¹¹C and ¹⁰C beams were produced and isotopically separated with the FRagment Separator (FRS) at GSI. The β ⁺-radioactive ion beams were produced with a beam purity of 99% for all the beam investigated (except one case where it was 94%) and intensities potentially sufficient to treat a small animal tumors within few minutes of irradiation time, ∼ 10⁶ particle per spill for the ¹⁰C and ∼ 10⁷ particle per spill for the ¹¹C beam, respectively. The impact of different ion optical parameters on the depth dose distribution was studied with a precision water column system. In this work, the measured depth dose distributions are presented together with results from Monte Carlo simulations using the FLUKA software.

In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam’s Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

Ultrasound-assisted carbon ion dosimetry and range measurement using injectable polymer-shelled phase-change nanodroplets: in vitro study

Methods allowing for in situ dosimetry and range verification are essential in radiotherapy to reduce the safety margins required to account for uncertainties introduced in the entire treatment workflow. This study suggests a non-invasive dosimetry concept for carbon ion radiotherapy based on phase-change ultrasound contrast agents. Injectable nanodroplets made of a metastable perfluorobutane (PFB) liquid core, stabilized with a crosslinked poly(vinylalcohol) shell, are vaporized at physiological temperature when exposed to carbon ion radiation (C-ions), converting them into echogenic microbubbles. Nanodroplets, embedded in tissue-mimicking phantoms, are exposed at 37 °C to a 312 MeV/u clinical C-ions beam at different doses between 0.1 and 4 Gy. The evaluation of the contrast enhancement from ultrasound imaging of the phantoms, pre- and post-irradiation, reveals a significant radiation-triggered nanodroplets vaporization occurring at the C-ions Bragg peak with sub-millimeter shift reproducibility and dose dependency. The specific response of the nanodroplets to C-ions is further confirmed by varying the phantom position, the beam range, and by performing spread-out Bragg peak irradiation. The nanodroplets' response to C-ions is influenced by their concentration and is dose rate independent. These early findings show the ground-breaking potential of polymer-shelled PFB nanodroplets to enable in vivo carbon ion dosimetry and range verification.

Proton therapy monitoring: spatiotemporal emission reconstruction with prompt gamma timing and implementation with PET detectors

Objective. In this study we introduce spatiotemporal emission reconstruction prompt gamma timing (SER-PGT), a new method to directly reconstruct the prompt photon emission in the space and time domains inside the patient in proton therapy.Approach. SER-PGT is based on the numerical optimisation of a multidimensional likelihood function, followed by a post-processing of the results. The current approach relies on a specific implementation of the maximum-likelihood expectation maximisation algorithm. The robustness of the method is guaranteed by the complete absence of any information about the target composition in the algorithm.Main results. Accurate Monte Carlo simulations indicate a range resolution of about 0.5 cm (standard deviation) when considering 10⁷primary protons impinging on an homogeneous phantom. Preliminary results on an anthropomorphic phantom are also reported.Significance. By showing the feasibility for the reconstruction of the primary particle range using PET detectors, this study provides significant basis for the development of an hybrid in-beam PET and prompt photon device.

Biological Rationale and Clinical Evidence of Carbon Ion Radiation Therapy for Adenoid Cystic Carcinoma: A Narrative Review

Adenoid cystic carcinoma (ACC) is a rare, basaloid, epithelial tumor, arising mostly from salivary glands. Radiation therapy can be employed as a single modality for unresectable tumors, in an adjuvant setting after uncomplete resection, in case of high-risk pathological features, or for recurrent tumors. Due to ACC intrinsic radioresistance, high linear energy transfer (LET) radiotherapy techniques have been evaluated for ACC irradiation: while fast neutron therapy has now been abandoned due to toxicity concerns, charged particle beams such as protons and carbon ions are at present the beams used for hadron therapy. Carbon ion radiation therapy (CIRT) is currently increasingly used for ACC irradiation. The aim of this review is to describe the immunological, molecular and clinicopathological bases that support ACC treatment with CIRT, as well as to expose the current clinical evidence that reveal the advantages of using CIRT for treating ACC.

Intra- and inter-fraction relative range verification in heavy-ion therapy using filtered interaction vertex imaging

Heavy-ion therapy, particularly using scanned (active) beam delivery, provides a precise and highly conformal dose distribution, with maximum dose deposition for each pencil beam at its endpoint (Bragg peak), and low entrance and exit dose. To take full advantage of this precision, robust range verification methods are required; these methods ensure that the Bragg peak is positioned correctly in the patient and the dose is delivered as prescribed. Relative range verification allows intra-fraction monitoring of Bragg peak spacing to ensure full coverage with each fraction, as well as inter-fraction monitoring to ensure all fractions are delivered consistently. To validate the proposed filtered interaction vertex imaging (IVI) method for relative range verification, a¹⁶O beam was used to deliver 12 Bragg peak positions in a 40 mm poly-(methyl methacrylate) phantom. Secondary particles produced in the phantom were monitored using position-sensitive silicon detectors. Events recorded on these detectors, along with a measurement of the treatment beam axis, were used to reconstruct the sites of origin of these secondary particles in the phantom. The distal edge of the depth distribution of these reconstructed points was determined with logistic fits, and the translation in depth required to minimize theχ²statistic between these fits was used to compute the range shift between any two Bragg peak positions. In all cases, the range shift was determined with sub-millimeter precision, to a standard deviation of the mean of 220(10)μm. This result validates filtered IVI as a reliable relative range verification method, which should be capable of monitoring each energy step in each fraction of a scanned heavy-ion treatment plan.

Metallic Nanoparticles: A Useful Prompt Gamma Emitter for Range Monitoring in Proton Therapy?

In clinical practice, dose delivery in proton therapy treatment is affected by uncertainties related to the range of the beam in the patient, which requires medical physicists to introduce safety margins on the penetration depth of the beam. Although this ensures an irradiation of the entire clinical target volume with the prescribed dose, these safety margins also lead to the exposure of nearby healthy tissues and a subsequent risk of side effects. Therefore, non-invasive techniques that allow for margin reduction through online monitoring of prompt gammas emitted along the proton tracks in the patient are currently under development. This study provides the proof-of-concept of metal-based nanoparticles, injected into the tumor, as a prompt gamma enhancer, helping in the beam range verification. It identifies the limitations of this application, suggesting a low feasibility in a realistic clinical scenario but opens some avenues for improvement.

Localization of anatomical changes in patients during proton therapy with in-beam PET monitoring: A voxel-based morphometry approach exploiting Monte Carlo simulations

Purpose

In‐beam positron emission tomography (PET) is one of the modalities that can be used for in vivo noninvasive treatment monitoring in proton therapy. Although PET monitoring has been frequently applied for this purpose, there is still no straightforward method to translate the information obtained from the PET images into easy‐to‐interpret information for clinical personnel. The purpose of this work is to propose a statistical method for analyzing in‐beam PET monitoring images that can be used to locate, quantify, and visualize regions with possible morphological changes occurring over the course of treatment.

Methods

We selected a patient treated for squamous cell carcinoma (SCC) with proton therapy, to perform multiple Monte Carlo (MC) simulations of the expected PET signal at the start of treatment, and to study how the PET signal may change along the treatment course due to morphological changes. We performed voxel‐wise two‐tailed statistical tests of the simulated PET images, resembling the voxel‐based morphometry (VBM) method commonly used in neuroimaging data analysis, to locate regions with significant morphological changes and to quantify the change.

Results

The VBM resembling method has been successfully applied to the simulated in‐beam PET images, despite the fact that such images suffer from image artifacts and limited statistics. Three dimensional probability maps were obtained, that allowed to identify interfractional morphological changes and to visualize them superimposed on the computed tomography (CT) scan. In particular, the characteristic color patterns resulting from the two‐tailed statistical tests lend themselves to trigger alarms in case of morphological changes along the course of treatment.

Conclusions

The statistical method presented in this work is a promising method to apply to PET monitoring data to reveal interfractional morphological changes in patients, occurring over the course of treatment. Based on simulated in‐beam PET treatment monitoring images, we showed that with our method it was possible to correctly identify the regions that changed. Moreover we could quantify the changes, and visualize them superimposed on the CT scan. The proposed method can possibly help clinical personnel in the replanning procedure in adaptive proton therapy treatments.

Nuclear Fragmentation Imaging for Carbon-Ion Radiation Therapy Monitoring: an In Silico Study

Purpose

This article presents an in vivo imaging technique based on nuclear fragmentation of carbon ions in irradiated tissues for potential real-time monitoring of carbon-ion radiation therapy (CIRT) treatment delivery and quality assurance purposes in clinical settings.

Materials and Methods

A proof-of-concept imaging and monitoring system (IMS) was devised to implement the technique. Monte Carlo simulations were performed for a prospective pencil-beam scanning CIRT nozzle. The development IMS benchmark considered a 5×5-cm² pixelated charged-particle detector stack positioned downstream from a target phantom and list-mode data acquisition. The abundance and production origins, that is, vertices, of the detected fragments were studied. Fragment trajectories were approximated by straight lines and a beam back-projection algorithm was built to reconstruct the vertices. The spatial distribution of the vertices was then used to determine plan relevant markers.

Results

The IMS technique was applied for a simulated CIRT case, a primary brain tumor. Four treatment plan monitoring markers were conclusively recovered: a depth dose distribution correlated profile, ion beam range, treatment target boundaries, and the beam spot position. Promising millimeter-scale (3-mm, ≤10% uncertainty) beam range and submillimeter (≤0.6-mm precision for shifts <3 cm) beam spot position verification accuracies were obtained for typical therapeutic energies between 150 and 290 MeV/u.

Conclusions

This work demonstrated a viable online monitoring technique for CIRT treatment delivery. The method's strong advantage is that it requires few signal inputs (position and timing), which can be simultaneously acquired with readily available technology. Future investigations will probe the technique's applicability to motion-sensitive organ sites and patient tissue heterogeneities. In-beam measurements with candidate detector-acquisition systems are ultimately essential to validate the IMS benchmark performance and subsequent deployment in the clinic.

Monitoring Carbon Ion Beams Transverse Position Detecting Charged Secondary Fragments: Results From Patient Treatment Performed at CNAO

Particle therapy in which deep seated tumours are treated using ¹²C ions (Carbon Ions RadioTherapy or CIRT) exploits the high conformity in the dose release, the high relative biological effectiveness and low oxygen enhancement ratio of such projectiles. The advantages of CIRT are driving a rapid increase in the number of centres that are trying to implement such technique. To fully profit from the ballistic precision achievable in delivering the dose to the target volume an online range verification system would be needed, but currently missing. The ¹²C ions beams range could only be monitored by looking at the secondary radiation emitted by the primary beam interaction with the patient tissues and no technical solution capable of the needed precision has been adopted in the clinical centres yet. The detection of charged secondary fragments, mainly protons, emitted by the patient is a promising approach, and is currently being explored in clinical trials at CNAO. Charged particles are easy to detect and can be back-tracked to the emission point with high efficiency in an almost background-free environment. These fragments are the product of projectiles fragmentation, and are hence mainly produced along the beam path inside the patient. This experimental signature can be used to monitor the beam position in the plane orthogonal to its flight direction, providing an online feedback to the beam transverse position monitor chambers used in the clinical centres. This information could be used to cross-check, validate and calibrate, whenever needed, the information provided by the ion chambers already implemented in most clinical centres as beam control detectors. In this paper we study the feasibility of such strategy in the clinical routine, analysing the data collected during the clinical trial performed at the CNAO facility on patients treated using ¹²C ions and monitored using the Dose Profiler (DP) detector developed within the INSIDE project. On the basis of the data collected monitoring three patients, the technique potential and limitations will be discussed.

Feasibility of quasi-prompt PET-based range verification in proton therapy

Compared to photon therapy, proton therapy allows a better conformation of the dose to the tumor volume with reduced radiation dose to co-irradiated tissues. In vivo verification techniques including positron emission tomography (PET) have been proposed as quality assurance tools to mitigate proton range uncertainties. Detection of differences between planned and actual dose delivery on a short timescale provides a fast trigger for corrective actions. Conventional PET-based imaging of ¹⁵O (T_1/2 = 2 min) and ¹¹C (T_1/2 = 20 min) distributions precludes such immediate feedback. We here present a demonstration of near real-time range verification by means of PET imaging of ¹²N (T_1/2 = 11 ms). PMMA and graphite targets were irradiated with a 150 MeV proton pencil beam consisting of a series of pulses of 10 ms beam-on and 90 ms beam-off. Two modules of a modified Siemens Biograph mCT PET scanner (21 × 21 cm² each), installed 25 cm apart, were used to image the beam-induced PET activity during the beam-off periods. The modifications enable the detectors to be switched off during the beam-on periods. ¹²N images were reconstructed using planar tomography. Using a 1D projection of the 2D reconstructed ¹²N image, the activity range was obtained from a fit of the activity profile with a sigmoid function. Range shifts due to modified target configurations were assessed for multiples of the clinically relevant 10⁸ protons per pulse (approximately equal to the highest intensity spots in the pencil beam scanning delivery of a dose of 1 Gy over a cubic 1 l volume). The standard deviation of the activity range, determined from 30 datasets obtained from three irradiations on PMMA and graphite targets, was found to be 2.5 and 2.6 mm (1σ) with 10⁸ protons per pulse and 0.9 and 0.8 mm (1σ) with 10⁹ protons per pulse. Analytical extrapolation of the results from this study shows that using a scanner with a solid angle coverage of 57%, with optimized detector switching and spot delivery times much smaller than the ¹²N half-life, an activity range measurement precision of 2.0 mm (1σ) and 1.3 mm (1σ) within 50 ms into an irradiation with 4 × 10⁷ and 10⁸ protons per pencil beam spot can be potentially realized. Aggregated imaging of neighboring spots or, if possible, increasing the number of protons for a few probe beam spots will enable the realization of higher precision range measurement.

Inter-fractional monitoring of [Formula: see text]C ions treatments: results from a clinical trial at the CNAO facility

The high dose conformity and healthy tissue sparing achievable in Particle Therapy when using C ions calls for safety factors in treatment planning, to prevent the tumor under-dosage related to the possible occurrence of inter-fractional morphological changes during a treatment. This limitation could be overcome by a range monitor, still missing in clinical routine, capable of providing on-line feedback. The Dose Profiler (DP) is a detector developed within the INnovative Solution for In-beam Dosimetry in hadronthErapy (INSIDE) collaboration for the monitoring of carbon ion treatments at the CNAO facility (Centro Nazionale di Adroterapia Oncologica) exploiting the detection of charged secondary fragments that escape from the patient. The DP capability to detect inter-fractional changes is demonstrated by comparing the obtained fragment emission maps in different fractions of the treatments enrolled in the first ever clinical trial of such a monitoring system, performed at CNAO. The case of a CNAO patient that underwent a significant morphological change is presented in detail, focusing on the implications that can be drawn for the achievable inter-fractional monitoring DP sensitivity in real clinical conditions. The results have been cross-checked against a simulation study.

The production of positron emitters with millisecond half-life during helium beam radiotherapy

Therapy with helium ions is currently receiving significantly increasing interest because helium ions have a sharper penumbra than protons and undergo less fragmentation than carbon ions and thus require less complicated dose calculations. For any ion of interest in hadron therapy, the accuracy of dose delivery is limited by range uncertainties. This has led to efforts by several groups to develop in vivo verification techniques, including positron emission tomography (PET), for monitoring of the dose delivery. Beam-on PET monitoring during proton therapy through the detection of short-lived positron emitters such as ¹²N (T _1/2  =  11 ms), an emerging PET technique, provides an attractive option given the achievable range accuracy, minimal susceptibility to biological washout and provision of near prompt feedback. Extension of this approach to helium ions requires information on the production yield of relevant short-lived positron emitters. This study presents the first measurements of the production of short-lived positron emitters in water, graphite, calcium and phosphorus targets irradiated with 59 MeV/u ³He and 50 MeV/u ⁴He beams. For these targets, the most produced short-lived nuclides are ¹³O/¹²N (T _1/2  =  8.6/11 ms) on water, ¹³O/¹²N on graphite, ⁴³Ti/⁴¹Sc/⁴²Sc (T _1/2  =  509-680 ms) on calcium, ²⁸P (T _1/2  =  268 ms) on phosphorus. A translation of the results from elemental targets to PMMA and representative tissues such as adipose tissue, muscle, compact and cortical bone, shows the dominance of ¹³O/¹²N in at least the first 20 s of an irradiation with ⁴He and somewhat longer with ³He. As the production of ¹³O/¹²N in a ³He irradiation is 3-4 times higher than in a ⁴He irradiation, from a statistical point of view, range verification using ¹³O/¹²N PET imaging will be about 2 times more precise for a ³He irradiation compared to a ⁴He irradiation.

Offline imaging of positron emitters induced by therapeutic helium, carbon and oxygen ion beams with a full-ring PET/CT scanner: experiments in reference targets

In vivo verification of light ion therapy based on positron-emission tomography (PET) imaging of irradiation induced patient activation relies on activity predictions from Monte-Carlo (MC) or analytical computational engines for comparison to the measurements. In order to achieve the necessary accuracy, experimental data are indispensable for the validation of the calculation models. For this we irradiated thick reference targets with mono-energetic helium, carbon and oxygen ion beams and measured the resulting material activation offline with a commercial full-ring PET/CT scanner located nearby the treatment room. Acquired PET data were analysed over time to separate the activity contribution of different radionuclides. Determined production yields were compared to published findings obtained from in-beam activation measurements with a limited-angle double-head PET camera. In addition, we investigated the time-dependence of the measured radionuclide-specific contributions and of the distal activity range, as well as the lateral spread of the activity signal as a function of beam penetration depth. We present radionuclide-specific depth-resolved activity distributions and production yields for the radionuclides [Formula: see text], [Formula: see text] and [Formula: see text], dominating irradiation-induced patient activation. We observe systematically lower production yields with a ratio between the dual-head and our full-ring PET measurements of, on average, 1.7 and 1.3 for the oxygen and carbon beam irradiations, and 1.7 (2.1) for the high (low) energy helium beam irradiations. Findings on the temporal development of the activity range confirm the expectation, with the oxygen beam induced signal being the most sensitive scenario. The experimental data reported in this work, acquired with a state-of-the-art full ring PET scanner, provide a comprehensive and consistent basis for the benchmarking of PET signal calculation engines. In particular, they can support a fine-tuning of the underlying physics models used by the respective implementation and therefore improve the accuracy of PET-based therapy verifications at current and future treatment facilities.

Secondary radiation measurements for particle therapy applications: Charged secondaries produced by 16O ion beams in a PMMA target at large angles

Particle therapy is a therapy technique that exploits protons or light ions to irradiate tumor targets with high accuracy. Protons and ¹²C ions are already used for irradiation in clinical routine, while new ions like ⁴He and ¹⁶O are currently being considered. Despite the indisputable physical and biological advantages of such ion beams, the planning of charged particle therapy treatments is challenged by range uncertainties, i.e. the uncertainty on the position of the maximal dose release (Bragg Peak - BP), during the treatment. To ensure correct 'in-treatment' dose deposition, range monitoring techniques, currently missing in light ion treatment techniques, are eagerly needed. The results presented in this manuscript indicate that charged secondary particles, mainly protons, produced by an ¹⁶O beam during target irradiation can be considered as candidates for ¹⁶O beam range monitoring. Hereafter, we report on the first yield measurements of protons, deuterons and tritons produced in the interaction of an ¹⁶O beam impinging on a PMMA target, as a function of detected energy and particle production position. Charged particles were detected at 90° and 60° with respect to incoming beam direction, and homogeneous and heterogeneous PMMA targets were used to probe the sensitivity of the technique to target inhomogeneities. The reported secondary particle yields provide essential information needed to assess the accuracy and resolution achievable in clinical conditions by range monitoring techniques based on secondary charged radiation.

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.