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

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

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

GOAL: In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes. METHODS/MATERIALS: We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans. RESULTS: Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields. CONCLUSION: The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.

A Data-Driven Fragmentation Model for Carbon Therapy GPU-Accelerated Monte-Carlo Dose Recalculation.

The advent of Graphics Processing Units (GPU) has prompted the development of Monte Carlo (MC) algorithms that can significantly reduce the simulation time with respect to standard MC algorithms based on Central Processing Unit (CPU) hardware. The possibility to evaluate a complete treatment plan within minutes, instead of hours, paves the way for many clinical applications where the time-factor is important. FRED (Fast paRticle thErapy Dose evaluator) is a software that exploits the GPU power to recalculate and optimise ion beam treatment plans. The main goal when developing the FRED physics model was to balance accuracy, calculation time and GPU execution guidelines. Nowadays, FRED is already used as a quality assurance tool in Maastricht and Krakow proton clinical centers and as a research tool in several clinical and research centers across Europe. Lately the core software has been updated including a model of carbon ions interactions with matter. The implementation is phenomenological and based on carbon fragmentation data currently available. The model has been tested against the MC FLUKA software, commonly used in particle therapy, and a good agreement was found. In this paper, the new FRED data-driven model for carbon ion fragmentation will be presented together with the validation tests against the FLUKA MC software. The results will be discussed in the context of FRED clinical applications to 12C ions treatment planning.

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.

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.

Biological Impact of Target Fragments on Proton Treatment Plans: An Analysis Based on the Current Cross-Section Data and a Full Mixed Field Approach.

Clinical routine in proton therapy currently neglects the radiobiological impact of nuclear target fragments generated by proton beams. This is partially due to the difficult characterization of the irradiation field. The detection of low energetic fragments, secondary protons and fragments, is in fact challenging due to their very short range. However, considering their low residual energy and therefore high LET, the possible contribution of such heavy particles to the overall biological effect could be not negligible. In this context, we performed a systematic analysis aimed at an explicit assessment of the RBE (relative biological effectiveness, i.e., the ratio of photon to proton physical dose needed to achieve the same biological effect) contribution of target fragments in the biological dose calculations of proton fields. The TOPAS Monte Carlo code has been used to characterize the radiation field, i.e., for the scoring of primary protons and fragments in an exemplary water target. TRiP98, in combination with LEM IV RBE tables, was then employed to evaluate the RBE with a mixed field approach accounting for fragments' contributions. The results were compared with that obtained by considering only primary protons for the pristine beam and spread out Bragg peak (SOBP) irradiations, in order to estimate the relative weight of target fragments to the overall RBE. A sensitivity analysis of the secondary particles production cross-sections to the biological dose has been also carried out in this study. Finally, our modeling approach was applied to the analysis of a selection of cell survival and RBE data extracted from published in vitro studies. Our results indicate that, for high energy proton beams, the main contribution to the biological effect due to the secondary particles can be attributed to secondary protons, while the contribution of heavier fragments is mainly due to helium. The impact of target fragments on the biological dose is maximized in the entrance channels and for small α/ß values. When applied to the description of survival data, model predictions including all fragments allowed better agreement to experimental data at high energies, while a minor effect was observed in the peak region. An improved description was also obtained when including the fragments' contribution to describe RBE data. Overall, this analysis indicates that a minor contribution can be expected to the overall RBE resulting from target fragments. However, considering the fragmentation effects can improve the agreement with experimental data for high energy proton beams.

A Novel Approach to Design and Evaluate BNCT Neutron Beams Combining Physical, Radiobiological, and Dosimetric Figures of Merit.

(1) Background:The quality of neutron beams for Boron Neutron Capture Therapy (BNCT) is currently defined by its physical characteristics in air. Recommendations exist to define whether a designed beam is useful for clinical treatment. This work presents a new way to evaluate neutron beams based on their clinical performance and on their safety, employing radiobiological quantities. (2) Methods: The case study is a neutron beam for deep-seated tumors from a 5 MeV proton beam coupled to a beryllium target. Physical Figures of Merit were used to design five beams; however, they did not allow a clear ranking of their quality in terms of therapeutic potential. The latter was then evaluated based on in-phantom dose distributions and on the calculation of the Uncomplicated Tumor Control Probability (UTCP). The safety of the beams was also evaluated calculating the in-patient out-of-beam dosimetry. (3) Results: All the beams ensured a UTCP comparable to the one of a clinical beam in phantom; the safety criterion allowed to choose the best candidate. When this was tested in the treatment planning of a real patient treated in Finland, the UTCP was still comparable to the one of the clinical beam. (4) Conclusions: Even when standard physical recommendations are not met, radiobiological and dosimetric criteria demonstrate to be a valid tool to select an effective and safe beam for patient treatment.

Deep Seated Tumour Treatments With Electrons of High Energy Delivered at FLASH Rates: The Example of Prostate Cancer.

Different therapies are adopted for the treatment of deep seated tumours in combination or as an alternative to surgical removal or chemotherapy: radiotherapy with photons (RT), particle therapy (PT) with protons or even heavier ions like 12C, are now available in clinical centres. In addition to these irradiation modalities, the use of Very High Energy Electron (VHEE) beams (100-200 MeV) has been suggested in the past, but the diffusion of that technique was delayed due to the needed space and budget, with respect to standard photon devices. These disadvantages were not paired by an increased therapeutic efficacy, at least when comparing to proton or carbon ion beams. In this contribution we investigate how recent developments in electron beam therapy could reshape the treatments of deep seated tumours. In this respect we carefully explored the application of VHEE beams to the prostate cancer, a well-known and studied example of deep seated tumour currently treated with high efficacy both using RT and PT. The VHEE Treatment Planning System was obtained by means of an accurate Monte Carlo (MC) simulation of the electrons interactions with the patient body. A simple model of the FLASH effect (healthy tissues sparing at ultra-high dose rates), has been introduced and the results have been compared with conventional RT. The study demonstrates that VHEE beams, even in absence of a significant FLASH effect and with a reduced energy range (70-130 MeV) with respect to implementations already explored in literature, could be a good alternative to standard RT, even in the framework of technological developments that are nowadays affordable.

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 12C 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 12C 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 12C 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.

Design of a BNCT irradiation room based on proton accelerator and beryllium target.

Preliminary studies for the design of an accelerator-based BNCT clinical facility are presented. The Beam Shaping Assembly neutron activation was evaluated experimentally and with Monte Carlo simulations. The activations of patient, air and walls in the room, the absorbed doses by the patient and the in-air dose distributions were evaluated. Based on these calculations, different walls compositions were tested to optimize the environmental conditions. Borated concrete, advantageously reducing the thermal flux in the room, was proven the best choice.

Biomedical Research Programs at Present and Future High-Energy Particle Accelerators.

Biomedical applications at high-energy particle accelerators have always been an important section of the applied nuclear physics research. Several new facilities are now under constructions or undergoing major upgrades. While the main goal of these facilities is often basic research in nuclear physics, they acknowledge the importance of including biomedical research programs and of interacting with other medical accelerator facilities providing patient treatments. To harmonize the programs, avoid duplications, and foster collaboration and synergism, the International Biophysics Collaboration is providing a platform to several accelerator centers with interest in biomedical research. In this paper, we summarize the programs of various facilities in the running, upgrade, or construction phase.

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

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

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

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

Addendum: Measurement of charged particle yields from PMMA irradiated by a 220 MeV/u 12C beam.

In this paper we report the re-analysis of the data published in (Piersanti et al. 2014) documenting the charged secondary particles production induced by the interaction of a 220 MeV/u 12C ion beam impinging on a polymethyl methacrylate (PMMA) target, measured in 2012 at the GSI facility in Darmstadt (Germany). This re-analysis takes into account the inhomogeneous light response of the LYSO crystal in the experimental setup measured in a subsequent experiment (2014) performed in the Heidelberg Ion- Beam Therapy Center. A better description of the detector and re-calculation of the geometrical efficiencies have been implemented as well, based on an improved approach that accounts also for the energy dependence of the emission spectrum. The new analysis has small effect on the total secondary charged flux, but has an impact on the production yield and emission velocity distributions of the different particle species (protons, deuterons and tritons) at different angles with respect to the beam direction (60° and 90°). All these observables indeed depend on the particle identification algorithms and hence on the LYSO detector energy response. The results of the data re-analysis presented here are intended to supersede and replace the results published in (Piersanti et al. 2014).

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

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

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.

The FLUKA Code: An Accurate Simulation Tool for Particle Therapy.

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

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.

FLUKA Monte Carlo simulation for the Leksell Gamma Knife Perfexion radiosurgery system: homogeneous media.

The purpose of this work is to investigate the capability of the FLUKA Monte Carlo (MC) code to simulate the Elekta Leksell Gamma Knife Perfexion (LGK-PFX) and reproduce the Treatment Planning System (TPS) Leksell GammaPlan version 8.2 (LGP) dose calculations for the case of a water equivalent phantom target. Thanks to the collaboration with Elekta Instruments AB, the collimation system geometry, the source positions and all the involved material have been simulated in detail. The relative linear dose distribution along the three coordinate axes, for each collimator size, and the Relative Output Factors (ROF) have been investigated. The simulation has been validated comparing simulated linear dose profiles with measurements performed with EBT radiochromic films. The acceptance criterion between experimental data and FLUKA results is based on the gamma index (GI) method. The FLUKA MC calculation for the ROF provided the values of 0.920 for the 8 mm collimators and 0.800 for the 4 mm collimators. These values are in good agreement with the Elekta reference data of 0.924 and 0.805 respectively. The percentage difference between calculated and reference values for the ROF is under 1% and within the FLUKA uncertainty. Also the simulated relative dose profiles show a good agreement with the LGP calculation expressed by means of the gamma index method. This established accuracy proves that FLUKA is a suitable and powerful tool in order to reproduce successfully the LGP calculations for the homogeneous media.