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

Innovations in dedicated PET instrumentation: from the operating room to specimen imaging

This review casts a spotlight on intraoperative positron emission tomography (PET) scanners and the distinctive challenges they confront. Specifically, these systems contend with the necessity of partial coverage geometry, essential for ensuring adequate access to the patient. This inherently leans them towards limited-angle PET imaging, bringing along its array of reconstruction and geometrical sensitivity challenges. Compounding this, the need for real-time imaging in navigation systems mandates rapid acquisition and reconstruction times. For these systems, the emphasis is on dependable PET image reconstruction (without significant artefacts) while rapid processing takes precedence over the spatial resolution of the system. In contrast, specimen PET imagers are unburdened by the geometrical sensitivity challenges, thanks to their ability to leverage full coverage PET imaging geometries. For these devices, the focus shifts: high spatial resolution imaging takes precedence over rapid image reconstruction. This review concurrently probes into the technical complexities of both intraoperative and specimen PET imaging, shedding light on their recent designs, inherent challenges, and technological advancements.

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.

Range shift verification in spot scanning proton therapy using gamma electron vertex imaging

BACKGROUND

In proton therapy, a highly steep distal dose penumbra can be utilized for dose conformity, given the Bragg peak characteristic of protons. However, the location of the Bragg peak in patients (i.e., the beam range) is very sensitive to range uncertainty. Even a small shift of beam range can produce a significant variation of delivered dose to tumor and normal tissues, thus degrading treatment quality and threatening patient safety. This range uncertainty issue, therefore, is one of the important aspects to be managed in proton therapy.

PURPOSE

For better management of range uncertainty, range verification has been widely studied, and prompt gamma imaging (PGI) is considered one of the promising methods in that effort. In this context, a PGI system named the gamma electron vertex imaging (GEVI) system was developed and recently upgraded for application to pencil-beam scanning (PBS) proton therapy. Here, we report the first experimental results using the therapeutic spot scanning proton beams.

METHODS

A homogeneous slab phantom and an anthropomorphic phantom were employed. Spherical and cubic planning target volumes (PTVs) were defined. Various range shift scenarios were introduced. Prompt gamma (PG) measurement was synchronized with beam irradiation. The measured PG distributions were aggregated to improve the PG statistics. The range shift was estimated based on the relative change of the centroid in the measured PG distribution. The estimated range shifts were analyzed by range shift mapping, confidence interval (CI) estimation, and statistical hypothesis testing.

RESULTS

The range shift mapping results showed an obvious measured range shift tendency following the true shift values. However, some fluctuations were found for spots that had still-low PG statistics after spot aggregation. The 99% CI distributions showed clearly distributed range shift measurement data. The overall accuracy and precision for all investigated scenarios were 0.36 and 0.20 mm, respectively. The results of one-sample t-tests confirmed that every shift scenario could be observed up to 1 mm of shift. The ANOVA results proved that the measured range shift data could be discriminated from one another, except for 16 (of 138) comparison cases having 1-2 mm shift differences.

CONCLUSIONS

This study demonstrated the feasibility of the GEVI system for measurement of range shift in spot scanning proton therapy. Our experimental results showed that the proton beam can be measured up to 1 mm of range shift with high accuracy and precision. We believe that the GEVI system is one of the most promising PGI systems for in vivo range verification. Further research for application to more various cases and patient treatments is planned.

Impact of and interplay between proton arc therapy and range uncertainties in proton therapy for head-and-neck cancer

Objective. Proton therapy reduces the integral dose to the patient compared to conventional photon treatments. However,in vivoproton range uncertainties remain a considerable hurdle. Range uncertainty reduction benefits depend on clinical practices. During intensity-modulated proton therapy (IMPT), the target is irradiated from only a few directions, but proton arc therapy (PAT), for which the target is irradiated from dozens of angles, may see clinical implementation by the time considerable range uncertainty reductions are achieved. It is therefore crucial to determine the impact of PAT on range uncertainty reduction benefits.Approach. For twenty head-and-neck cancer patients, four different treatment plans were created: an IMPT and a PAT treatment plan assuming current clinical range uncertainties of 3.5% (IMPT3.5%and PAT3.5%), and an IMPT and a PAT treatment plan assuming that range uncertainties can be reduced to 1% (IMPT1%and PAT1%). Plans were evaluated with respect to target coverage and organ-at-risk doses as well as normal tissue complication probabilities (NTCPs) for parotid glands (endpoint: parotid gland flow <25%) and larynx (endpoint: larynx edema).Main results. Implementation of PAT (IMPT3.5%-PAT3.5%) reduced mean NTCPs in the nominal and worst-case scenario by 3.2 percentage points (pp) and 4.2 pp, respectively. Reducing range uncertainties from 3.5% to 1% during use of IMPT (IMPT3.5%-IMPT1%) reduced evaluated NTCPs by 0.9 pp and 2.0 pp. Benefits of range uncertainty reductions subsequently to PAT implementation (PAT3.5%-PAT1%) were 0.2 pp and 1.0 pp, with considerably higher benefits in bilateral compared to unilateral cases.Significance. The mean clinical benefit of implementing PAT was more than twice as high as the benefit of a 3.5%-1% range uncertainty reduction. Range uncertainty reductions are expected to remain beneficial even after PAT implementation, especially in cases with target positions allowing for full leveraging of the higher number of gantry angles during PAT.

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.

The quest for multifunctional and dedicated PET instrumentation with irregular geometries

We focus on reviewing state-of-the-art developments of dedicated PET scanners with irregular geometries and the potential of different aspects of multifunctional PET imaging. First, we discuss advances in non-conventional PET detector geometries. Then, we present innovative designs of organ-specific dedicated PET scanners for breast, brain, prostate, and cardiac imaging. We will also review challenges and possible artifacts by image reconstruction algorithms for PET scanners with irregular geometries, such as non-cylindrical and partial angular coverage geometries and how they can be addressed. Then, we attempt to address some open issues about cost/benefits analysis of dedicated PET scanners, how far are the theoretical conceptual designs from the market/clinic, and strategies to reduce fabrication cost without compromising performance.

Near-field coded-mask technique and its potential for proton therapy monitoring

Objective.Prompt-gamma imaging encompasses several approaches to the online monitoring of the beam range or deposited dose distribution in proton therapy. We test one of the imaging techniques - a coded mask approach - both experimentally and via simulations.Approach.Two imaging setups have been investigated experimentally. Each of them comprised a structured tungsten collimator in the form of a modified uniformly redundant array mask and a LYSO:Ce scintillation detector of fine granularity. The setups differed in detector dimensions and operation mode (1D or 2D imaging). A series of measurements with radioactive sources have been conducted, testing the performance of the setups for near-field gamma imaging. Additionally, Monte Carlo simulations of a larger setup of the same type were conducted, investigating its performance with a realistic gamma source distribution occurring during proton therapy.Main results.The images of point-like sources reconstructed from two small-scale prototypes' data using the maximum-likelihood expectation maximisation algorithm constitute the experimental proof of principle for the near-field coded-mask imaging modality, both in the 1D and the 2D mode. Their precision allowed us to calibrate out certain systematic offsets appearing due to the limited alignment accuracy of setup elements. The simulation of the full-scale setup yielded a mean distal falloff retrieval precision of 0.72 mm in the studies for beam energy range 89.5-107.9 MeV and with 1 × 108protons (a typical number for distal spots). The implemented algorithm of image reconstruction is relatively fast-a typical procedure needs several seconds.Significance.Coded-mask imaging appears a valid option for proton therapy monitoring. The results of simulations let us conclude that the proposed full-scale setup is competitive with the knife-edge-shaped and the multi-parallel slit cameras investigated by other groups.

Tumour status prediction by means of carbon-ion beam irradiation: comparison of washout rates between in-beam PET and DCE-MRI in rats

Objective.Tumour response to radiation therapy appears as changes in tumour vascular condition. There are several methods for analysing tumour blood circulatory changes one of which is dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), but there is no method that can observe the tumour vascular condition and physiological changes at the site of radiation therapy. Positron emission tomography (PET) has been applied for treatment verification in charged particle therapy, which is based on the detection of positron emitters produced through nuclear fragmentation reactions in a patient's body. However, the produced positron emitters are washed out biologically depending on the tumour vascular condition. This means that measuring the biological washout rate may allow evaluation of the tumour radiation response, in a similar manner to DCE-MRI. Therefore, this study compared the washout rates in rats between in-beam PET during12C ion beam irradiation and DCE-MRI.Approach.Different vascular conditions of the tumour model were prepared for six nude rats. The tumour of each nude rat was irradiated by a12C ion beam with simultaneous in-beam PET measurement. In 10-12 h, the DCE-MRI experiment was performed for the same six nude rats. The biological washout rate of the produced positron emitters (k2,1st) and the MRI contrast agent (k2a) were derived using the single tissue compartment model.Main results.A linear correlation was observed betweenk2,1standk2a, and they were inversely related to fractional necrotic volume.Significance.This is the first animal study which confirmed the biological washout rate of in-beam PET correlates closely with tumour vascular condition measured with the MRI contrast agent administrated intravenously.

A simulation study of in-beam visualization system for proton therapy by monitoring scattered protons

Recently, in-beam positron emission tomography (PET) has been actively researched for reducing biological washout effects and dose monitoring during irradiation. However, the positron distribution does not precisely reflect the dose distribution since positron production and ionization are completely different physical processes. Thus, a novel in-beam system was proposed to determine proton dose range by measuring scattered protons with dozens of scintillation detectors surrounding the body surface. While previous studies conducted a preliminary experiment with a simple phantom, we simulated more complex situations in this paper. Especially, we conducted three stepwise simulation studies to demonstrate the feasibility of the proposed method. First, a simple rectangular phantom was reproduced on simulation and irradiated with protons for obtaining current values and Monte Carlo (MC) dose. Next, we trained a deep learning model to estimate 2-dimensional-dose range (2D-DL dose) from measured current values for simulation (A). We simulated plastic scintillators as detectors to measure the scattered protons. Second, a rectangular phantom with an air layer was used, and 3D-DL dose was estimated in simulation (B). Finally, a cylindrical phantom that mimics the human body was used for confirming the estimation quality of the simulation (C). Consequently, the position of the Bragg peak was estimated with an error of 1.0 mm in simulation (A). In addition, the position of the air layer, as well as the verifying peak position with an error of 2.1 mm, was successfully estimated in simulation (B). Although the estimation error of the peak position was 12.6 mm in simulation (C), the quality was successfully further improved to 9.3 mm by incorporating the mass density distribution obtained from the computed tomography (CT). These simulation results demonstrated the potential of the as-proposed verification system. Additionally, the effectiveness of CT utilization for estimating the DL dose was also indicated.

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.

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.

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.

Range uncertainty reductions in proton therapy may lead to the feasibility of novel beam arrangements which improve organ-at-risk sparing

PURPOSE

In proton therapy, dose distributions are currently often conformed to organs at risk (OARs) using the less sharp dose fall-off at the lateral beam edge to reduce the effects of uncertainties in the in vivo proton range. However, range uncertainty reductions may make greater use of the sharper dose fall-off at the distal beam edge feasible, potentially improving OAR sparing. We quantified the benefits of such novel beam arrangements.

METHODS

For each of 10 brain or skull base cases, five treatment plans robust to 2 mm setup and 0%-4% range uncertainty were created for the traditional clinical beam arrangement and a novel beam arrangement making greater use of the distal beam edge to conform the dose distribution to the brainstem. Metrics including the brainstem normal tissue complication probability (NTCP) with the endpoint of necrosis were determined for all plans and all setup and range uncertainty scenarios.

RESULTS

For the traditional beam arrangement, reducing the range uncertainty from the current level of approximately 4% to a potentially achievable level of 1% reduced the brainstem NTCP by up to 0.9 percentage points in the nominal and up to 1.5 percentage points in the worst-case scenario. Switching to the novel beam arrangement at 1% range uncertainty improved these values by a factor of 2, that is, to 1.8 percentage points and 3.2 percentage points, respectively. The novel beam arrangement achieved a lower brainstem NTCP in all cases starting at a range uncertainty of 2%.

CONCLUSION

The benefits of novel beam arrangements may be of the same magnitude or even exceed the direct benefits of range uncertainty reductions. Indirect effects may therefore contribute markedly to the benefits of reducing proton range uncertainties.

Optically stimulated luminescence in state-of-the-art LYSO:Ce scintillators enables high spatial resolution 3D dose imaging

In this contribution, we study the optically stimulated luminescence (OSL) exhibited by commercial \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Lu}_{(2-x)}\hbox {Y}_x\hbox {SiO}_5$$\end{document}Lu(2-x)YxSiO5:Ce crystals. This photon emission mechanism, complementary to scintillation, can trap a fraction of radiation energy deposited in the material and provides sufficient signal to develop a novel post-irradiation 3D dose readout. We characterize the OSL emission through spectrally and temporally resolved measurements and monitor the dose linearity response over a broad range. The measurements show that the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Ce}^{3+}$$\end{document}Ce3+ centers responsible for scintillation also function as recombination centers for the OSL mechanism. The capture to OSL-active traps competes with scintillation originating from the direct non-radiative energy transfer to the luminescent centers. An OSL response on the order of 100 ph/MeV is estimated. We demonstrate the imaging capabilities provided by such an OSL photon yield using a proof-of-concept optical readout method. A 0.1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {mm}^3$$\end{document}mm3 spatial resolution for doses as low as 0.5 Gy is projected using a cubic crystal to image volumetric dose profiles. While OSL degrades the intrinsic scintillating performance by reducing the number of scintillation photons emitted following the passage of ionizing radiation, it can encode highly resolved spatial information of the interaction point of the particle. This feature combines ionizing radiation spectroscopy and 3D reusable dose imaging in a single material.

In vivo production of fluorine-18 in a chicken egg tumor model of breast cancer for proton therapy range verification

Range verification of clinical protontherapy systems via positron-emission tomography (PET) is not a mature technology, suffering from two major issues: insufficient signal from low-energy protons in the Bragg peak area and biological washout of PET emitters. The use of contrast agents including ¹⁸O, ⁶⁸Zn or ⁶³Cu, isotopes with a high cross section for low-energy protons in nuclear reactions producing PET emitters, has been proposed to enhance the PET signal in the last millimeters of the proton path. However, it remains a challenge to achieve sufficient concentrations of these isotopes in the target volume. Here we investigate the possibilities of ¹⁸O-enriched water (18-W), a potential contrast agent that could be incorporated in large proportions in live tissues by replacing regular water. We hypothesize that 18-W could also mitigate the problem of biological washout, as PET (¹⁸F) isotopes created inside live cells would remain trapped in the form of fluoride anions (F-), allowing its signal to be detected even hours after irradiation. To test our hypothesis, we designed an experiment with two main goals: first, prove that 18-W can incorporate enough ¹⁸O into a living organism to produce a detectable signal from ¹⁸F after proton irradiation, and second, determine the amount of activity that remains trapped inside the cells. The experiment was performed on a chicken embryo chorioallantoic membrane tumor model of head and neck cancer. Seven eggs with visible tumors were infused with 18-W and irradiated with 8-MeV protons (range in water: 0.74 mm), equivalent to clinical protons at the end of particle range. The activity produced after irradiation was detected and quantified in a small-animal PET-CT scanner, and further studied by placing ex-vivo tumours in a gamma radiation detector. In the acquired images, specific activity of ¹⁸F (originating from 18-W) could be detected in the tumour area of the alive chicken embryo up to 9 h after irradiation, which confirms that low-energy protons can indeed produce a detectable PET signal if a suitable contrast agent is employed. Moreover, dynamic PET studies in two of the eggs evidenced a minimal effect of biological washout, with 68% retained specific ¹⁸F activity at 8 h after irradiation. Furthermore, ex-vivo analysis of 4 irradiated tumours showed that up to 3% of oxygen atoms in the targets were replaced by ¹⁸O from infused 18-W, and evidenced an entrapment of 59% for specific activity of ¹⁸F after washing, supporting our hypothesis that F- ions remain trapped within the cells. An infusion of 18-W can incorporate ¹⁸O in animal tissues by replacing regular water inside cells, producing a PET signal when irradiated with low-energy protons that could be used for range verification in protontherapy. ¹⁸F produced inside cells remains entrapped and suffers from minimal biological washout, allowing for a sharper localization with longer PET acquisitions. Further studies must evaluate the feasibility of this technique in dosimetric conditions closer to clinical practice, in order to define potential protocols for its use in patients.

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.

Management of Motion and Anatomical Variations in Charged Particle Therapy: Past, Present, and Into the Future

The major aim of radiation therapy is to provide curative or palliative treatment to cancerous malignancies while minimizing damage to healthy tissues. Charged particle radiotherapy utilizing carbon ions or protons is uniquely suited for this task due to its ability to achieve highly conformal dose distributions around the tumor volume. For these treatment modalities, uncertainties in the localization of patient anatomy due to inter- and intra-fractional motion present a heightened risk of undesired dose delivery. A diverse range of mitigation strategies have been developed and clinically implemented in various disease sites to monitor and correct for patient motion, but much work remains. This review provides an overview of current clinical practices for inter and intra-fractional motion management in charged particle therapy, including motion control, current imaging and motion tracking modalities, as well as treatment planning and delivery techniques. We also cover progress to date on emerging technologies including particle-based radiography imaging, novel treatment delivery methods such as tumor tracking and FLASH, and artificial intelligence and discuss their potential impact towards improving or increasing the challenge of motion mitigation in charged particle therapy.

Applications of Machine Learning to Improve the Clinical Viability of Compton Camera Based in vivo Range Verification in Proton Radiotherapy

We studied the application of a deep, fully connected Neural Network (NN) to process prompt gamma (PG) data measured by a Compton camera (CC) during the delivery of clinical proton radiotherapy beams. The network identifies 1) recorded "bad" PG events arising from background noise during the measurement, and 2) the correct ordering of PG interactions in the CC to help improve the fidelity of "good" data used for image reconstruction. PG emission from a tissue-equivalent target during irradiation with a 150 MeV proton beam delivered at clinical dose rates was measured with a prototype CC. Images were reconstructed from both the raw measured data and the measured data that was further processed with a neural network (NN) trained to identify "good" and "bad" PG events and predict the ordering of individual interactions within the good PG events. We determine if NN processing of the CC data could improve the reconstructed PG images to a level in which they could provide clinically useful information about the in vivo range and range shifts of the proton beams delivered at full clinical dose rates. Results showed that a deep, fully connected NN improved the achievable contrast to noise ratio (CNR) in our images by more than a factor of 8x. This allowed the path, range, and lateral width of the clinical proton beam within a tissue equivalent target to easily be identified from the PG images, even at the highest dose rates of a 150 MeV proton beam used for clinical treatments. On average, shifts in the beam range as small as 3 mm could be identified. However, when limited by the amount of PG data measured with our prototype CC during the delivery of a single proton pencil beam (~1 × 10⁹ protons), the uncertainty in the reconstructed PG images limited the identification of range shift to ~5 mm. Substantial improvements in CC images were obtained during clinical beam delivery through NN pre-processing of the measured PG data. We believe this shows the potential of NNs to help improve and push CC-based PG imaging toward eventual clinical application for proton RT treatment delivery 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.

The effects of Compton camera data acquisition and readout timing on PG imaging for proton range verification

The purpose of this study was to determine how the characteristics of the data acquisition (DAQ) electronics of a Compton camera (CC) affect the quality of the recorded prompt gamma (PG) interaction data and the reconstructed images, during clinical proton beam delivery. We used the Monte-Carlo-plus-Detector-Effect (MCDE) model to simulate the delivery of a 150 MeV clinical proton pencil beam to a tissue-equivalent plastic phantom. With the MCDE model we analyzed how the recorded PG interaction data changed as two characteristics of the DAQ electronics of a CC were changed: (1) the number of data readout channels; and (2) the active charge collection, readout, and reset time. As the proton beam dose rate increased, the number of recorded PG single-, double-, and triple-scatter events decreased by a factor of 60× for the current DAQ configuration of the CC. However, as the DAQ readout channels were increased and the readout/reset timing decreased, the number of recorded events decreased by <5× at the highest clinical dose rate. The increased number of readout channels and reduced readout/reset timing also resulted in higher quality recorded data. That is, a higher percentage of the recorded double- and triple-scatters were "true" events (caused by a single incident gamma) and not "false" events (caused by multiple incident gammas). The increase in the number and the quality of recorded data allowed higher quality PG images to be reconstructed even at the highest clinical dose rates.

Hadron Therapy Achievements and Challenges: The CNAO Experience

Protons and carbon ions (hadrons) have useful properties for the treatments of patients affected by oncological pathologies. They are more precise than conventional X-rays and possess radiobiological characteristics suited for treating radio-resistant or inoperable tumours. This paper gives an overview of the status of hadron therapy around the world. It focusses on the Italian National Centre for Oncological Hadron therapy (CNAO), introducing operation procedures, system performance, expansion projects, methodologies and modelling to build individualized treatments. There is growing evidence that supports safety and effectiveness of hadron therapy for a variety of clinical situations. However, there is still a lack of high-level evidence directly comparing hadron therapy with modern conventional radiotherapy techniques. The results give an overview of pre-clinical and clinical research studies and of the treatments of 3700 patients performed at CNAO. The success and development of hadron therapy is strongly associated with the creation of networks among hadron therapy facilities, clinics, universities and research institutions. These networks guarantee the growth of cultural knowledge on hadron therapy, favour the efficient recruitment of patients and present available competences for R&D (Research and Development) programmes.

Radioactivity and Space Range of Ultra-Low-Activity for in vivo Off-line PET Verification of Proton and Carbon Ion Beam-A Phantom Study

Purpose: The radioactivity induced by proton and heavy ion beam belongs to the ultra-low-activity (ULA). Therefore, the radioactivity and space range of commercial off-line positron emission tomography (PET) acquisition based on ULA should be evaluated accurately to guarantee the reliability of clinical verification. The purpose of this study is to quantify the radioactivity and space range of off-line PET acquisition by simulating the ULA triggered by proton and heavy ion beam.

Methods: PET equipment validation phantom and low activity ¹⁸F-FDG were used to simulate the ULA with radioactivity of 11.1–1480 Bq/mL. The radioactivity of ULA was evaluated by comparing the radioactivity in the images with the values calculated from the decay function with a radioactivity error tolerance of 5%. The space range of ULA was evaluated by comparing the width of the R50 analyzed activity distribution curve with the actual width of the container with a space range error tolerance of 4 mm.

Results: When radioactivity of ULA was >148 Bq/mL, the radioactivity error was <5%. When radioactivity of ULA was >30 Bq/mL, the space range error was below 4 mm.

Conclusions: Off-line PET can be used to quantify the radioactivity of proton and heavy ion beam when the ULA exceeds 148 Bq/mL, both in radioactivity and in space range.

Experimental demonstration of accurate Bragg peak localization with ionoacoustic tandem phase detection (iTPD)

Accurate knowledge of the exact stopping location of ions inside the patient would allow full exploitation of their ballistic properties for patient treatment. The localized energy deposition of a pulsed particle beam induces a rapid temperature increase of the irradiated volume and leads to the emission of ionoacoustic (IA) waves. Detecting the time-of-flight (ToF) of the IA wave allows inferring information on the Bragg peak location and can henceforth be used forin-vivorange verification. A challenge for IA is the poor signal-to-noise ratio at clinically relevant doses and viable machines. We present a frequency-based measurement technique, labeled as ionoacoustic tandem phase detection (iTPD) utilizing lock-in amplifiers. The phase shift of the IA signal to a reference signal is measured to derive theToF. Experimental IA measurements with a 3.5 MHz lead zirconate titanate (PZT) transducer and lock-in amplifiers were performed in water using 22 MeV proton bursts. A digital iTPD was performedin-silicoat clinical dose levels on experimental data obtained from a clinical facility and secondly, on simulations emulating a heterogeneous geometry. For the experimental setup using 22 MeV protons, a localization accuracy and precision obtained through iTPD deviates from a time-based reference analysis by less than 15μm. Several methodological aspects were investigated experimentally in systematic manner. Lastly, iTPD was evaluatedin-silicofor clinical beam energies indicating that iTPD is in reach of sub-mm accuracy for fractionated doses < 5 Gy. iTPD can be used to accurately measure theToFof IA signals online via its phase shift in frequency domain. An application of iTPD to the clinical scenario using a single pulsed beam is feasible but requires further development to reach <1 Gy detection capabilities.

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.

A proof-of-concept study of an in-situ partial-ring time-of-flight PET scanner for proton beam verification

Development of a PET system capable of in-situ imaging requires a design that can accommodate the proton treatment beam nozzle. Among the several PET instrumentation approaches developed thus far, the dual-panel PET scanner is often used as it is simpler to develop and integrate within the proton therapy gantry. Partial-angle coverage of these systems can however lead to limited-angle artefacts in the reconstructed PET image. We have previously demonstrated via simulations that time-of-flight (TOF) reconstruction reduces the artifacts accompanying limited-angle data, and permits proton range measurement with 1-2 mm accuracy and precision. In this work we show measured results from a small proof-of-concept dual-panel PET system that uses TOF information to reconstruct PET data acquired after proton irradiation. The PET scanner comprises of two detector modules, each comprised of an array of 4×4×30 mm³ lanthanum bromide scintillator. Measurements are performed with an oxygen-rich gel-water, an adipose tissue equivalent material, and in vitro tissue phantoms. For each phantom measurement, 2 Gy dose was deposited using 54 - 100 MeV proton beams. For each phantom, a Monte Carlo simulation generating the expected distribution of PET isotope from the corresponding proton irradiation was also performed. Proton range was calculated by drawing multiple depth-profiles over a central region encompassing the proton dose deposition. For each profile, proton range was calculated using two techniques (a) 50% pick-off from the distal edge of the profile, and (b) comparing the measured and Monte Carlo profile to minimize the absolute sum of differences over the entire profile. A 10 min PET acquisition acquired with minimal delay post proton-irradiation is compared with a 10 min PET scan acquired after a 20 min delay. Measurements show that PET acquisition with minimal delay is necessary to collect ¹⁵O signal, and maximize ¹¹C signal collection with a short PET acquisition. In comparison with the 50% pick-off technique, the shift technique is more robust and offers better precision in measuring the proton range for the different phantoms. Range measurements from PET images acquired with minimal delay, and the shift technique demonstrate the ability to achieve <1.5 mm accuracy and precision in estimating proton range.

A Comparative Study of Two In Vivo PET Verification Methods in Clinical Cases

Purpose

Positron emission tomography (PET) range verification is an important method that can help improve the confidence in proton therapy for clinical applications. Two kinds of verification methods are implemented and compared based on clinical cases in this study.

Method

The study is conducted on 14 breast cancer patients following proton irradiation treatment. Verification is done by calculating the depth error between the numerically predicted values with the measured PET image along the beam direction. Point-based and segment-based methods are applied and compared. The verification results are presented as depth error means and standard deviations in a region of interest (ROI).

Results

The mean value of the depth error of all 14 cases is within the range of [-3, 3] mm for both point-based and segment-based methods, and only one case result calculated by the point-based method is slightly beyond -3 mm. When comparing the mean depth error from the two methods, the paired t-test result shows that the p-value is 0.541, and the standard deviation of the segment-based method is smaller than that of the point-based method.

Conclusion

In breast cancer case verification application, point-based and segment-based methods show no significant difference in the mean value of results. Both methods can quantify the accuracy of proton radiotherapy to the millimeter level.

Radioactive Beams for Image-Guided Particle Therapy: The BARB Experiment at GSI

Several techniques are under development for image-guidance in particle therapy. Positron (β⁺) emission tomography (PET) is in use since many years, because accelerated ions generate positron-emitting isotopes by nuclear fragmentation in the human body. In heavy ion therapy, a major part of the PET signals is produced by β⁺-emitters generated via projectile fragmentation. A much higher intensity for the PET signal can be obtained using β⁺-radioactive beams directly for treatment. This idea has always been hampered by the low intensity of the secondary beams, produced by fragmentation of the primary, stable beams. With the intensity upgrade of the SIS-18 synchrotron and the isotopic separation with the fragment separator FRS in the FAIR-phase-0 in Darmstadt, it is now possible to reach radioactive ion beams with sufficient intensity to treat a tumor in small animals. This was the motivation of the BARB (Biomedical Applications of Radioactive ion Beams) experiment that is ongoing at GSI in Darmstadt. This paper will present the plans and instruments developed by the BARB collaboration for testing the use of radioactive beams in cancer therapy.

Proton range uncertainty reduction benefits for skull base tumors in terms of normal tissue complication probability (NTCP) and healthy tissue doses

PURPOSE

Proton therapy allows for more conformal dose distributions and lower organ at risk and healthy tissue doses than conventional photon-based radiotherapy, but uncertainties in the proton range currently prevent proton therapy from making full use of these advantages. Numerous developments therefore aim to reduce such range uncertainties. In this work, we quantify the benefits of reductions in range uncertainty for treatments of skull base tumors.

METHODS

The study encompassed 10 skull base patients with clival tumors. For every patient, six treatment plans robust to setup errors of 2 mm and range errors from 0% to 5% were created. The determined metrics included the brainstem and optic chiasm normal tissue complication probability (NTCP) with the endpoints of necrosis and blindness, respectively, as well as the healthy tissue volume receiving at least 70% of the prescription dose.

RESULTS

A range uncertainty reduction from the current level of 4% to a potentially achievable level of 1% reduced the probability of brainstem necrosis by up to 1.3 percentage points in the nominal scenario in which neither setup nor range errors occur and by up to 2.9 percentage points in the worst-case scenario. Such a range uncertainty reduction also reduced the optic chiasm NTCP with the endpoint of blindness by up to 0.9 percentage points in the nominal scenario and by up to 2.2 percentage points in the worst-case scenario. The decrease in the healthy tissue volume receiving at least 70% of the prescription dose ranged from -7.8 to 24.1 cc in the nominal scenario and from -3.4 to 38.4 cc in the worst-case scenario.

CONCLUSION

The benefits quantified as part of this study serve as a guideline of the OAR and healthy tissue dose benefits that range monitoring techniques may be able to achieve. Benefits were observed between all levels of range uncertainty. Even smaller range uncertainty reductions may therefore be beneficial.

Towards real-time PGS range monitoring in proton therapy of prostate cancer

Proton therapy of prostate cancer (PCPT) was linked with increased levels of gastrointestinal toxicity in its early use compared to intensity-modulated radiation therapy (IMRT). The higher radiation dose to the rectum by proton beams is mainly due to anatomical variations. Here, we demonstrate an approach to monitor rectal radiation exposure in PCPT based on prompt gamma spectroscopy (PGS). Endorectal balloons (ERBs) are used to stabilize prostate movement during radiotherapy. These ERBs are usually filled with water. However, other water solutions containing elements with higher atomic numbers, such as silicon, may enable the use of PGS to monitor the radiation exposure of the rectum. Protons hitting silicon atoms emit prompt gamma rays with a specific energy of 1.78 MeV, which can be used to monitor whether the ERB is being hit. In a binary approach, we search the silicon energy peaks for every irradiated prostate region. We demonstrate this technique for both single-spot irradiation and real treatment plans. Real-time feedback based on the ERB being hit column-wise is feasible and would allow clinicians to decide whether to adapt or continue treatment. This technique may be extended to other cancer types and organs at risk, such as the oesophagus.

Prompt-gamma emission in GEANT4 revisited and confronted with experiment

PURPOSE

The field of online monitoring of the beam range is one of the most researched topics in proton therapy over the last decade. The development of detectors that can be used for beam range verification under clinical conditions is a challenging task. One promising possible solution are modalities that record prompt-gamma radiation produced by the interactions of the proton beam with the target tissue. A good understanding of the energy spectra of the prompt gammas and the yields in certain energy regions is crucial for a successful design of a prompt-gamma detector. Monte-Carlo simulations are an important tool in development and testing of detector concepts, thus the proper modelling of the prompt-gamma emission in those simulations are of vital importance. In this paper, we confront a number of GEANT4 simulations of prompt-gamma emission, performed with different versions of the package and different physics lists, with experimental data obtained from a phantom irradiation with proton beams of four different energies in the range 70-230 MeV.

METHODS

The comparison is made on different levels: features of the prompt-gamma energy spectrum, gamma emission depth profiles for discrete transitions and the width of the distal fall-off in those profiles.

RESULTS

The best agreement between the measurements and the simulations is found for the GEANT4 version 10.4.2 and the reference physics list QGSP_BIC_HP.

CONCLUSIONS

Modifications to prompt-gamma emission modelling in higher versions of the software increase the discrepancy between the simulation results and the experimental data.

Optimisation of the event-based TOF filtered back-projection for online imaging in total-body J-PET

We perform a parametric study of the newly developed time-of-flight (TOF) image reconstruction algorithm, proposed for the real-time imaging in total-body Jagiellonian PET (J-PET) scanners. The asymmetric 3D filtering kernel is applied at each most likely position of electron-positron annihilation, estimated from the emissions of back-to-back γ-photons. The optimisation of its parameters is studied using Monte Carlo simulations of a 1-mm spherical source, NEMA IEC and XCAT phantoms inside the ideal J-PET scanner. The combination of high-pass filters which included the TOF filtered back-projection (FBP), resulted in spatial resolution, 1.5 times higher in the axial direction than for the conventional 3D FBP. For realistic 10-minute scans of NEMA IEC and XCAT, which require a trade-off between the noise and spatial resolution, the need for Gaussian TOF kernel components, coupled with median post-filtering, is demonstrated. The best sets of 3D filter parameters were obtained by the Nelder-Mead minimisation of the mean squared error between the resulting and reference images. The approach allows training the reconstruction algorithm for custom scans, using the IEC phantom, when the temporal resolution is below 50 ps. The image quality parameters, estimated for the best outcomes, were systematically better than for the non-TOF FBP.

Dose Calculation Algorithms for External Radiation Therapy: An Overview for Practitioners

Radiation therapy (RT) is a constantly evolving therapeutic technique; improvements are continuously being introduced for both methodological and practical aspects. Among the features that have undergone a huge evolution in recent decades, dose calculation algorithms are still rapidly changing. This process is propelled by the awareness that the agreement between the delivered and calculated doses is of paramount relevance in RT, since it could largely affect clinical outcomes. The aim of this work is to provide an overall picture of the main dose calculation algorithms currently used in RT, summarizing their underlying physical models and mathematical bases, and highlighting their strengths and weaknesses, referring to the most recent studies on algorithm comparisons. This handy guide is meant to provide a clear and concise overview of the topic, which will prove useful in helping clinical medical physicists to perform their responsibilities more effectively and efficiently, increasing patient benefits and improving the overall quality of the management of radiation treatment.

Quality assurance method for monitoring of lateral pencil beam positions in scanned carbon-ion radiotherapy using tracking of secondary ions

PURPOSE

Ion beam radiotherapy offers enhances dose conformity to the tumor volume while better sparing healthy tissue compared to conventional photon radiotherapy. However, the increased dose gradient also makes it more sensitive to uncertainties. While the most important uncertainty source is the patient itself, the beam delivery is also subject to uncertainties. Most of the proton therapy centers used cyclotrons, which deliver typically a stable beam over time, allowing a continuous extraction of the beam. Carbon-ion beam radiotherapy (CIRT) in contrast uses synchrotrons and requires a larger and energy-dependent extrapolation of the nozzle-measured positions to obtain the lateral beam positions in the isocenter, since the nozzle-to-isocenter distance is larger than for cyclotrons. Hence, the control of lateral pencil beam positions at isocenter in CIRT is more sensitive to uncertainties than in proton radiotherapy. Therefore, an independent monitoring of the actual lateral positions close to the isocenter would be very valuable and provide additional information. However, techniques capable to do so are scarce, and they are limited in precision, accuracy and effectivity.

METHODS

The detection of secondary ions (charged nuclear fragments) has previously been exploited for the Bragg peak position of C-ion beams. In our previous work, we investigated for the first time the feasibility of lateral position monitoring of pencil beams in CIRT. However, the reported precision and accuracy were not sufficient for a potential implementation into clinical practice. In this work, it is shown how the performance of the method is improved to the point of clinical relevance. To minimize the observed uncertainties, a mini-tracker based on hybrid silicon pixel detectors was repositioned downstream of an anthropomorphic head phantom. However, the secondary-ion fluence rate in the mini-tracker rises up to 1.5 × 10⁵ ions/s/cm² , causing strong pile-up of secondary-ion signals. To solve this problem, we performed hardware changes, optimized the detector settings, adjusted the setup geometry and developed new algorithms to resolve ambiguities in the track reconstruction. The performance of the method was studied on two treatment plans delivered with a realistic dose of 3 Gy (RBE) and averaged dose rate of 0.27 Gy/s at the Heidelberg Ion-Beam Therapy Center (HIT) in Germany. The measured lateral positions were compared to reference beam positions obtained either from the beam nozzle or from a multi-wire proportional chamber positioned at the room isocenter.

RESULTS

The presented method is capable to simultaneously monitor both lateral pencil beam coordinates over the entire tumor volume during the treatment delivery, using only a 2-cm² mini-tracker. The effectivity (defined as the fraction of analyzed pencil beams) was 100%. The reached precision of (0.6 to 1.5) mm and accuracy of (0.5 to 1.2) mm are in line with the clinically accepted uncertainty for QA measurements of the lateral pencil beam positions.

CONCLUSIONS

It was demonstrated that the performance of the method for a non-invasive lateral position monitoring of pencil beams is sufficient for a potential clinical implementation. The next step is to evaluate the method clinically in a group of patients in a future observational clinical study.

A time-of-flight-based reconstruction for real-time prompt-gamma imaging in proton therapy

We propose a novel prompt-gamma (PG) imaging modality for real-time monitoring in proton therapy: PG time imaging (PGTI). By measuring the time-of-flight (TOF) between a beam monitor and a PG detector, our goal is to reconstruct the PG vertex distribution in 3D. In this paper, a dedicated, non-iterative reconstruction strategy is proposed (PGTI reconstruction). Here, it was resolved under a 1D approximation to measure a proton range shift along the beam direction. In order to show the potential of PGTI in the transverse plane, a second method, based on the calculation of the centre of gravity (COG) of the TIARA pixel detectors' counts was also explored. The feasibility of PGTI was evaluated in two different scenarios. Under the assumption of a 100 ps (rms) time resolution (achievable in single proton regime), MC simulations showed that a millimetric proton range shift is detectable at 2σwith 10⁸incident protons in simplified simulation settings. With the same proton statistics, a potential 2 mm sensitivity (at 2σwith 10⁸incident protons) to beam displacements in the transverse plane was found using the COG method. This level of precision would allow to act in real-time if the treatment does not conform to the treatment plan. A worst case scenario of a 1 ns (rms) TOF resolution was also considered to demonstrate that a degraded timing information can be compensated by increasing the acquisition statistics: in this case, a 2 mm range shift would be detectable at 2σwith 10⁹incident protons. By showing the feasibility of a time-based algorithm for the reconstruction of the PG vertex distribution for a simplified anatomy, this work poses a theoretical basis for the future development of a PG imaging detector based on the measurement of particle TOF.

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.

Proton range verification with MACACO II Compton camera enhanced by a neural network for event selection

The applicability extent of hadron therapy for tumor treatment is currently limited by the lack of reliable online monitoring techniques. An active topic of investigation is the research of monitoring systems based on the detection of secondary radiation produced during treatment. MACACO, a multi-layer Compton camera based on LaBr₃ scintillator crystals and SiPMs, is being developed at IFIC-Valencia for this purpose. This work reports the results obtained from measurements of a 150 MeV proton beam impinging on a PMMA target. A neural network trained on Monte Carlo simulations is used for event selection, increasing the signal to background ratio before image reconstruction. Images of the measured prompt gamma distributions are reconstructed by means of a spectral reconstruction code, through which the 4.439 MeV spectral line is resolved. Images of the emission distribution at this energy are reconstructed, allowing calculation of the distal fall-off and identification of target displacements of 3 mm.

Novel On-line PET Imaging for Intra-Beam Range Verification and Delivery Optimization: A Simulation Feasibility Study

On-line PET image-based method uses an initial particle beam to measure the particle beam range (BR) within the same fraction so that any measured range-shift with respect to the predicted BR can be compensated before the rest therapeutic beam deliveries. However, the method requires to use a low-dose initial beam to minimize the risk of beam overshooting, which leads to low image count and inaccurate BR measurement. In this in-silico study, we evaluated the feasibility of a new on-line PET imaging method that measures BR at the mid-plane of a target volume with part of the high-dose therapy beams to verify BR and guide adaptive treatment re-planning. Simulations included various processes of proton beam radiations to a tumor inside a human brain phantom, positron and PET image generation at the mid-plane with initial beams, activity range measurement, and range-shift compensated beam delivery. The results demonstrated that the new method, under the simulated conditions, can achieve ~1.1 mm mid-plane BR measurement accuracy and closely match the delivered range-shift compensated dose distribution with the planned one. Overall, it is promising that this new method may significantly improve particle therapy accuracy.

Denoising PET images for proton therapy using a residual U-net

The use of proton therapy has the advantage of high dose concentration as it is possible to concentrate the dose on the tumor while suppressing damage to the surrounding normal organs. However, the range uncertainty significantly affects the actual dose distribution in the vicinity of the proton range, limiting the benefit of proton therapy for reducing the dose to normal organs. By measuring the annihilation gamma rays from the produced positron emitters, it is possible to obtain a proton induced positron emission tomography (pPET) image according to the irradiation region of the proton beam. Smoothing with a Gaussian filter is generally used to denoise PET images; however, this approach lowers the spatial resolution. Furthermore, other conventional smoothing processing methods may deteriorate the steep region of the pPET images. In this study, we proposed a denoising method based on a Residual U-Net for pPET images. We conducted the Monte Carlo simulation and irradiation experiment on a human phantom to obtain pPET data. The accuracy of the range estimation and the image similarity were evaluated for pPET images using the Residual U-Net, a Gaussian filter, a median filter, the block-matching and 3D-filtering (BM3D), and a total variation (TV) filter. Usage of the Residual U-Net yielded effective results corresponding to the range estimation; however, the results of peak-signal-to-noise ratio were identical to those for the Gaussian filter, median filter, BM3D, and TV filter. The proposed method can contribute to improving the accuracy of treatment verification and shortening the PET measurement time.

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.

A new detector for the beam energy measurement in proton therapy: a feasibility study

The proof of concept of a new device, capable of determining in a few seconds the energy of clinical proton beams by measuring the time of flight (ToF) of protons, is presented. The prototype consists of two thin ultra fast silicon detector (UFSD) pads, aligned along the beam direction in a telescope configuration and readout by a digitizer. The method developed for extracting the energy at the isocenter from the measured ToF, validated by Monte Carlo simulations, and the procedure used to calibrate the system are also presented and discussed in detail. The prototype was tested at the Centro Nazionale di Adroterapia Oncologica (CNAO, Pavia, Italy), at several beam energies, covering the entire clinical range, and using different distances between the sensors. The measured beam energies were benchmarked against the nominal CNAO energy values, obtained during the commissioning of the centre from the measured ranges in water. Deviations of few hundreds of keV have been achieved for all considered proton beam energies for distances between the two sensors larger than 60 cm, indicating a sensitivity to the corresponding beam range in water smaller than the clinical tolerance of 1 mm. Moreover, few seconds of irradiation were necessary to collect the required statistics. These preliminary results indicate that a telescope of UFSDs could achieve in a short time the accuracy required for the clinical application and therefore encourage further investigations towards the improvement and the optimization of the present prototype.

Carbon-11 and Carbon-12 beam range verifications through prompt gamma and annihilation gamma measurements: Monte Carlo simulations

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

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.

Extension of the ML-EM algorithm for dose estimation using PET in proton therapy: application to an inhomogeneous target

Positron emission tomography (PET) has been used for in vivo treatment verification, mainly for range verification, in proton therapy. Evaluating the direct dose from PET measurements remains challenging; however, it is highly desirable from a clinical perspective. In this study, a method for estimating the dose distribution from the positron emitter distributions was developed using the maximum likelihood expectation maximization algorithm. The 1D spatial relationship between positron emitter distributions and a dose distribution in an inhomogeneous target was inputted into the system matrix based on a filter framework. In contrast, spatial resolution of the PET system and total variation regularization (as prior knowledge for dose distribution) were considered in the 3D image-space. The dose estimation was demonstrated using Monte Carlo simulated PET activity distributions with substantial noise in a head and neck phantom. This mimicked the single field irradiation of the spread-out Bragg peak beams at clinical dose levels. Besides the simple implementation of the algorithm, this strategy achieved a high-speed calculation (30 s for a 3D dose estimation) and accurate dose and range estimations (less than 10% and 2 mm errors at 1-σ values, respectively). The proposed method could be key for using PET for in vivo dose monitoring.