Interaction vertex imaging (IVI) for carbon ion therapy monitoring: a feasibility study

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.

Uncertainty-aware spot rejection rate as quality metric for proton therapy using a digital tracking calorimeter

Objective.Proton therapy is highly sensitive to range uncertainties due to the nature of the dose deposition of charged particles. To ensure treatment quality, range verification methods can be used to verify that the individual spots in a pencil beam scanning treatment fraction match the treatment plan. This study introduces a novel metric for proton therapy quality control based on uncertainties in range verification of individual spots.Approach.We employ uncertainty-aware deep neural networks to predict the Bragg peak depth in an anthropomorphic phantom based on secondary charged particle detection in a silicon pixel telescope designed for proton computed tomography. The subsequently predicted Bragg peak positions, along with their uncertainties, are compared to the treatment plan, rejecting spots which are predicted to be outside the 95% confidence interval. The such-produced spot rejection rate presents a metric for the quality of the treatment fraction.Main results.The introduced spot rejection rate metric is shown to be well-defined for range predictors with well-calibrated uncertainties. Using this method, treatment errors in the form of lateral shifts can be detected down to 1 mm after around 1400 treated spots with spot intensities of 1 × 107protons. The range verification model used in this metric predicts the Bragg peak depth to a mean absolute error of 1.107 ± 0.015 mm.Significance.Uncertainty-aware machine learning has potential applications in proton therapy quality control. This work presents the foundation for future developments in this area.

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 × 109 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.

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.

Investigation of Suitable Detection Angles for Carbon-Ion Radiotherapy Monitoring in Depth by Means of Secondary-Ion Tracking

The dose conformity of carbon-ion beam radiotherapy, which allows the reduction of the dose deposition in healthy tissue and the escalation of the dose to the tumor, is associated with a high sensitivity to anatomical changes during and between treatment irradiations. Thus, the monitoring of inter-fractional anatomical changes is crucial to ensure the dose conformity, to potentially reduce the size of the safety margins around the tumor and ultimately to reduce the irradiation of healthy tissue. To do so, monitoring methods of carbon-ion radiotherapy in depth using secondary-ion tracking are being investigated. In this work, the detection and localization of a small air cavity of 2 mm thickness were investigated at different detection angles of the mini-tracker relative to the beam axis. The experiments were conducted with a PMMA head phantom at the Heidelberg Ion-Beam Therapy Center (HIT) in Germany. In a clinic-like irradiation of a single field of 3 Gy (RBE), secondary-ion emission profiles were measured by a 2 cm2 mini-tracker composed of two silicon pixel detectors. Two positions of the cavity in the head phantom were studied: in front and in the middle of the tumor volume. The significance of the cavity detection was found to be increased at smaller detection angles, while the accuracy of the cavity localization was improved at larger detection angles. Detection angles of 20° - 30° were found to be a good compromise for accessing both, the detectability and the position of the air cavity along the depth in the head of a patient.

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.

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.

Measurement of PET isotope production cross sections for protons and carbon ions on carbon and oxygen targets for applications in particle therapy range verification

Measured cross sections for the production of the PET isotopes [Formula: see text], [Formula: see text] and [Formula: see text] from carbon and oxygen targets induced by protons (40-220 [Formula: see text]) and carbon ions (65-430 [Formula: see text]) are presented. These data were obtained via activation measurements of irradiated graphite and beryllium oxide targets using a set of three scintillators coupled by a coincidence logic. The measured cross sections are relevant for the PET particle range verification method where accurate predictions of the [Formula: see text] emitter distribution produced by therapeutic beams in the patient tissue are required. The presented dataset is useful for validation and optimization of the nuclear reaction models within Monte Carlo transport codes. For protons the agreement of a radiation transport calculation using the measured cross sections with a thick target PET measurement is demonstrated.

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.

A Monte Carlo feasibility study for neutron based real-time range verification in proton therapy

Uncertainties in the proton range in tissue during proton therapy limit the precision in treatment delivery. These uncertainties result in expanded treatment margins, thereby increasing radiation dose to healthy tissue. Real-time range verification techniques aim to reduce these uncertainties in order to take full advantage of the finite range of the primary protons. In this paper, we propose a novel concept for real-time range verification based on detection of secondary neutrons produced in nuclear interactions during proton therapy. The proposed detector concept is simple; consisting of a hydrogen-rich converter material followed by two charged particle tracking detectors, mimicking a proton recoil telescopic arrangement. Neutrons incident on the converter material are converted into protons through elastic and inelastic (n,p) interactions. The protons are subsequently detected in the tracking detectors. The information on the direction and position of these protons is then utilized in a new reconstruction algorithm to estimate the depth distribution of neutron production by the proton beam, which in turn is correlated with the primary proton range. In this paper, we present the results of a Monte Carlo feasibility study and show that the proposed concept could be used for real-time range verification with millimetric precision in proton therapy.

Mid-range probing-towards range-guided particle therapy

Particle therapy can achieve excellent dose localization but is sensitive to range uncertainty. Therefore, online in vivo range verification before treatment is critical for treatment safety and quality assurance. We introduce a novel range-probing technique that uses mid-range treatment spots selected from the treatment plan as probing beams to be delivered before other treatment spots in pencil beam scanning. The probing spot signal can be acquired by an in-beam positron emission tomography (PET) scanner, and the reconstructed spot positions are compared with pre-calculated positions to measure the range shift. Mid-range probing ensures that the Bragg peaks stay inside the tumor even with significant range variation from the plan. Single-layered spots enable easier spot detection than multi-layered spots without cross-layered spot smearing. With therapeutic dose, the probing beam offers higher positron activities and range detectability than the low-dose imaging beam by up to two orders of magnitude, without exposing patients to extra radiation. Higher positron activities allow sufficient signal statistics in shorter acquisition time, therefore reducing metabolic washout of positron emitters. Thus, range shifts from the plan can be measured easily. We also describe two online range-compensated plan modification methods. We apply correction, if the range shift is above a certain tolerance. We studied feasibility using simulated particle treatment plans with online anatomical changes. For illustration, we demonstrate range shift measurement using simulated probing dose. The proposed range probing and correction effectively handled range shifts in the simulated cases. Both range-compensated adaptation and optimization accounted for online changes so that the delivered dose matched the planned dose. With a dedicated online in-beam PET scanner and phantom and clinical studies, which are currently being developed, this novel strategy may open up a range-guided particle therapy¹ paradigm.

Secondary radiation measurements for particle therapy applications: charged particles produced by 4He and 12C ion beams in a PMMA target at large angle

Proton and carbon ion beams are used in the clinical practice for external radiotherapy treatments achieving, for selected indications, promising and superior clinical results with respect to x-ray based radiotherapy. Other ions, like [Formula: see text] have recently been considered as projectiles in particle therapy centres and might represent a good compromise between the linear energy transfer and the radiobiological effectiveness of [Formula: see text] ion and proton beams, allowing improved tumour control probability and minimising normal tissue complication probability. All the currently used p, [Formula: see text] and [Formula: see text] ion beams allow achieving sharp dose gradients on the boundary of the target volume, however the accurate dose delivery is sensitive to the patient positioning and to anatomical variations with respect to photon therapy. This requires beam range and/or dose release measurement during patient irradiation and therefore the development of dedicated monitoring techniques. All the proposed methods make use of the secondary radiation created by the beam interaction with the patient and, in particular, in the case of [Formula: see text] ion beams are also able to exploit the significant charged radiation component. Measurements performed to characterise the charged secondary radiation created by [Formula: see text] and [Formula: see text] particle therapy beams are reported. Charged secondary yields, energy spectra and emission profiles produced in a poly-methyl methacrylate (PMMA) target by [Formula: see text] and [Formula: see text] beams of different therapeutic energies were measured at 60° and 90° with respect to the primary beam direction. The secondary yield of protons produced along the primary beam path in a PMMA target was obtained. The energy spectra of charged secondaries were obtained from time-of-flight information, whereas the emission profiles were reconstructed exploiting tracking detector information. The obtained measurements are in agreement with results reported in the literature and suggests the feasibility of range monitoring based on charged secondary particle detection: the implications for particle therapy monitoring applications are also discussed.

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

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

MONDO: a neutron tracker for particle therapy secondary emission characterisation

Tumour control is performed in particle therapy using particles and ions, whose high irradiation precision enhances the effectiveness of the treatment, while sparing the healthy tissue surrounding the target volume. Dose range monitoring devices using photons and charged particles produced by the beam interacting with the patient's body have already been proposed, but no attempt has been made yet to exploit the detection of the abundant neutron component. Since neutrons can release a significant dose far away from the tumour region, precise measurements of their flux, production energy and angle distributions are eagerly sought in order to improve the treatment planning system (TPS) software. It will thus be possible to predict not only the normal tissue toxicity in the target region, but also the risk of late complications in the whole body. The aforementioned issues underline the importance of an experimental effort devoted to the precise characterisation of neutron production, aimed at the measurement of their abundance, emission point and production energy. The technical challenges posed by a neutron detector aimed at high detection efficiency and good backtracking precision are addressed within the MONDO (monitor for neutron dose in hadrontherapy) project, whose main goal is to develop a tracking detector that can target fast and ultrafast neutrons. A full reconstruction of two consecutive elastic scattering interactions undergone by the neutrons inside the detector material will be used to measure their energy and direction. The preliminary results of an MC simulation performed using the FLUKA software are presented here, together with the DSiPM (digital SiPM) readout implementation. New detector readout implementations specifically tailored to the MONDO tracker are also discussed, and the neutron detection efficiency attainable with the proposed neutron tracking strategy are reported.

Secondary radiation measurements for particle therapy applications: prompt photons produced by 4He, 12C and 16O ion beams in a PMMA target

Charged particle beams are used in particle therapy (PT) to treat oncological patients due to their selective dose deposition in tissues with respect to the photons and electrons used in conventional radiotherapy. Heavy (Z  >  1) PT beams can additionally be exploited for their high biological effectiveness in killing cancer cells. Nowadays, protons and carbon ions are used in PT clinical routines. Recently, interest in the potential application of helium and oxygen beams has been growing. With respect to protons, such beams are characterized by their reduced multiple scattering inside the body, increased linear energy transfer, relative biological effectiveness and oxygen enhancement ratio. The precision of PT demands online dose monitoring techniques, crucial to improving the quality assurance of any treatment: possible patient mis-positioning and biological tissue changes with respect to the planning CT scan could negatively affect the outcome of the therapy. The beam range confined in the irradiated target can be monitored thanks to the neutral or charged secondary radiation emitted by the interactions of hadron beams with matter. Among these secondary products, prompt photons are produced by nuclear de-excitation processes, and at present, different dose monitoring and beam range verification techniques based on prompt-γ detection are being proposed. It is hence of importance to perform γ yield measurement in therapeutic-like conditions. In this paper we report on the yields of prompt photons produced by the interaction of helium, carbon and oxygen ion beams with a poly-methyl methacrylate (PMMA) beam stopping target. The measurements were performed at the Heidelberg Ion-Beam Therapy Center (HIT) with beams of different energies. An LYSO scintillator, placed at [Formula: see text] and [Formula: see text] with respect to the beam direction, was used as the photon detector. The obtained γ yields for the carbon ion beams are compared with results from the literature, while no other results from helium and oxygen beams have been published yet. A discussion on the expected resolution of a slit camera detector is presented, demonstrating the feasibility of a prompt-γ-based monitoring technique for PT treatments using helium, carbon and oxygen ion beams.

Design of a new tracking device for on-line beam range monitor in carbon therapy

Charged particle therapy is a technique for cancer treatment that exploits hadron beams, mostly protons and carbon ions. A critical issue is the monitoring of the beam range so to check the correct dose deposition to the tumor and surrounding tissues. The design of a new tracking device for beam range real-time monitoring in pencil beam carbon ion therapy is presented. The proposed device tracks secondary charged particles produced by beam interactions in the patient tissue and exploits the correlation of the charged particle emission profile with the spatial dose deposition and the Bragg peak position. The detector, currently under construction, uses the information provided by 12 layers of scintillating fibers followed by a plastic scintillator and a pixelated Lutetium Fine Silicate (LFS) crystal calorimeter. An algorithm to account and correct for emission profile distortion due to charged secondaries absorption inside the patient tissue is also proposed. Finally detector reconstruction efficiency for charged particle emission profile is evaluated using a Monte Carlo simulation considering a quasi-realistic case of a non-homogenous phantom.

Nuclear physics in particle therapy: a review

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

Monitoring of Hadrontherapy Treatments by Means of Charged Particle Detection

The interaction of the incoming beam radiation with the patient body in hadrontherapy treatments produces secondary charged and neutral particles, whose detection can be used for monitoring purposes and to perform an on-line check of beam particle range. In the context of ion-therapy with active scanning, charged particles are potentially attractive since they can be easily tracked with a high efficiency, in presence of a relatively low background contamination. In order to verify the possibility of exploiting this approach for in-beam monitoring in ion-therapy, and to guide the design of specific detectors, both simulations and experimental tests are being performed with ion beams impinging on simple homogeneous tissue-like targets (PMMA). From these studies, a resolution of the order of few millimeters on the single track has been proven to be sufficient to exploit charged particle tracking for monitoring purposes, preserving the precision achievable on longitudinal shape. The results obtained so far show that the measurement of charged particles can be successfully implemented in a technology capable of monitoring both the dose profile and the position of the Bragg peak inside the target and finally lead to the design of a novel profile detector. Crucial aspects to be considered are the detector positioning, to be optimized in order to maximize the available statistics, and the capability of accounting for the multiple scattering interactions undergone by the charged fragments along their exit path from the patient body. The experimental results collected up to now are also valuable for the validation of Monte Carlo simulation software tools and their implementation in Treatment Planning Software packages.

Monte Carlo Simulations of Particle Interactions with Tissue in Carbon Ion Therapy

Monte Carlo simulations are increasingly considered the most accurate tool for calculating particle interactions with tissue. This contribution reviews the basics of Monte Carlo methods and their emerging role for application to several areas of macroscopic simulation in the worldwide rapidly growing field of carbon ion therapy, spanning from dosimetric calculations to imaging of secondary radiation.

Assessment of Geant4 Prompt-Gamma Emission Yields in the Context of Proton Therapy Monitoring

Monte Carlo tools have been long used to assist the research and development of solutions for proton therapy monitoring. The present work focuses on the prompt-gamma emission yields by comparing experimental data with the outcomes of the current version of Geant4 using all applicable proton inelastic models. For the case in study and using the binary cascade model, it was found that Geant4 overestimates the prompt-gamma emission yields by 40.2 ± 0.3%, even though it predicts the prompt-gamma profile length of the experimental profile accurately. In addition, the default implementations of all proton inelastic models show an overestimation in the number of prompt gammas emitted. Finally, a set of built-in options and physically sound Geant4 source code changes have been tested in order to try to improve the discrepancy observed. A satisfactory agreement was found when using the QMD model with a wave packet width equal to 1.3 fm².

Time-resolved imaging of prompt-gamma rays for proton range verification using a knife-edge slit camera based on digital photon counters

Proton range monitoring may facilitate online adaptive proton therapy and improve treatment outcomes. Imaging of proton-induced prompt gamma (PG) rays using a knife-edge slit collimator is currently under investigation as a potential tool for real-time proton range monitoring. A major challenge in collimated PG imaging is the suppression of neutron-induced background counts. In this work, we present an initial performance test of two knife-edge slit camera prototypes based on arrays of digital photon counters (DPCs). PG profiles emitted from a PMMA target upon irradiation with a 160 MeV proton pencil beams (about 6.5 × 10(9) protons delivered in total) were measured using detector modules equipped with four DPC arrays coupled to BGO or LYSO : Ce crystal matrices. The knife-edge slit collimator and detector module were placed at 15 cm and 30 cm from the beam axis, respectively, in all cases. The use of LYSO : Ce enabled time-of-flight (TOF) rejection of background events, by synchronizing the DPC readout electronics with the 106 MHz radiofrequency signal of the cyclotron. The signal-to-background (S/B) ratio of 1.6 obtained with a 1.5 ns TOF window and a 3 MeV-7 MeV energy window was about 3 times higher than that obtained with the same detector module without TOF discrimination and 2 times higher than the S/B ratio obtained with the BGO module. Even 1 mm shifts of the Bragg peak position translated into clear and consistent shifts of the PG profile if TOF discrimination was applied, for a total number of protons as low as about 6.5 × 10(8) and a detector surface of 6.6 cm × 6.6 cm.

Range Verification Methods in Particle Therapy: Underlying Physics and Monte Carlo Modeling

Hadron therapy allows for highly conformal dose distributions and better sparing of organs-at-risk, thanks to the characteristic dose deposition as function of depth. However, the quality of hadron therapy treatments is closely connected with the ability to predict and achieve a given beam range in the patient. Currently, uncertainties in particle range lead to the employment of safety margins, at the expense of treatment quality. Much research in particle therapy is therefore aimed at developing methods to verify the particle range in patients. Non-invasive in vivo monitoring of the particle range can be performed by detecting secondary radiation, emitted from the patient as a result of nuclear interactions of charged hadrons with tissue, including β (+) emitters, prompt photons, and charged fragments. The correctness of the dose delivery can be verified by comparing measured and pre-calculated distributions of the secondary particles. The reliability of Monte Carlo (MC) predictions is a key issue. Correctly modeling the production of secondaries is a non-trivial task, because it involves nuclear physics interactions at energies, where no rigorous theories exist to describe them. The goal of this review is to provide a comprehensive overview of various aspects in modeling the physics processes for range verification with secondary particles produced in proton, carbon, and heavier ion irradiation. We discuss electromagnetic and nuclear interactions of charged hadrons in matter, which is followed by a summary of some widely used MC codes in hadron therapy. Then, we describe selected examples of how these codes have been validated and used in three range verification techniques: PET, prompt gamma, and charged particle detection. We include research studies and clinically applied methods. For each of the techniques, we point out advantages and disadvantages, as well as clinical challenges still to be addressed, focusing on MC simulation aspects.

A review of the use and potential of the GATE Monte Carlo simulation code for radiation therapy and dosimetry applications

In this paper, the authors' review the applicability of the open-source GATE Monte Carlo simulation platform based on the GEANT4 toolkit for radiation therapy and dosimetry applications. The many applications of GATE for state-of-the-art radiotherapy simulations are described including external beam radiotherapy, brachytherapy, intraoperative radiotherapy, hadrontherapy, molecular radiotherapy, and in vivo dose monitoring. Investigations that have been performed using GEANT4 only are also mentioned to illustrate the potential of GATE. The very practical feature of GATE making it easy to model both a treatment and an imaging acquisition within the same framework is emphasized. The computational times associated with several applications are provided to illustrate the practical feasibility of the simulations using current computing facilities.

Automation and uncertainty analysis of a method for in-vivo range verification in particle therapy

We introduce the automation of the range difference calculation deduced from particle-irradiation induced β(+)-activity distributions with the so-called most-likely-shift approach, and evaluate its reliability via the monitoring of algorithm- and patient-specific uncertainty factors. The calculation of the range deviation is based on the minimization of the absolute profile differences in the distal part of two activity depth profiles shifted against each other. Depending on the workflow of positron emission tomography (PET)-based range verification, the two profiles under evaluation can correspond to measured and simulated distributions, or only measured data from different treatment sessions. In comparison to previous work, the proposed approach includes an automated identification of the distal region of interest for each pair of PET depth profiles and under consideration of the planned dose distribution, resulting in the optimal shift distance. Moreover, it introduces an estimate of uncertainty associated to the identified shift, which is then used as weighting factor to 'red flag' problematic large range differences. Furthermore, additional patient-specific uncertainty factors are calculated using available computed tomography (CT) data to support the range analysis. The performance of the new method for in-vivo treatment verification in the clinical routine is investigated with in-room PET images for proton therapy as well as with offline PET images for proton and carbon ion therapy. The comparison between measured PET activity distributions and predictions obtained by Monte Carlo simulations or measurements from previous treatment fractions is performed. For this purpose, a total of 15 patient datasets were analyzed, which were acquired at Massachusetts General Hospital and Heidelberg Ion-Beam Therapy Center with in-room PET and offline PET/CT scanners, respectively. Calculated range differences between the compared activity distributions are reported in a 2D map in beam-eye-view. In comparison to previously proposed approaches, the new most-likely-shift method shows more robust results for assessing in-vivo the range from strongly varying PET distributions caused by differing patient geometry, ion beam species, beam delivery techniques, PET imaging concepts and counting statistics. The additional visualization of the uncertainties and the dedicated weighting strategy contribute to the understanding of the reliability of observed range differences and the complexity in the prediction of activity distributions. The proposed method promises to offer a feasible technique for clinical routine of PET-based range verification.

Particle Detector Applications in Medicine

Selected particle detectors are described which find an application in medicine and have been the topic of presentations at the 2013 Vienna Conference of Instrumentation (VCI).