IOeRT conventional and FLASH treatment planning system implementation exploiting fast GPU Monte Carlo: The case of breast cancer

Partial breast irradiation for the treatment of early-stage breast cancer patients can be performed by means of Intra Operative electron Radiation Therapy (IOeRT). One of the main limitations of this technique is the absence of a treatment planning system (TPS) that could greatly help in ensuring a proper coverage of the target volume during irradiation. An IOeRT TPS has been developed using a fast Monte Carlo (MC) and an ultrasound imaging system to provide the best irradiation strategy (electron beam energy, applicator position and bevel angle) and to facilitate the optimisation of dose prescription and delivery to the target volume while maximising the organs at risk sparing. The study has been performed in silico, exploiting MC simulations of a breast cancer treatment. Ultrasound-based input has been used to compute the absorbed dose maps in different irradiation strategies and a quantitative comparison between the different options was carried out using Dose Volume Histograms. The system was capable of exploring different beam energies and applicator positions in few minutes, identifying the best strategy with an overall computation time that was found to be completely compatible with clinical implementation. The systematic uncertainty related to tissue deformation during treatment delivery with respect to imaging acquisition was taken into account. The potential and feasibility of a GPU based full MC TPS implementation of IOeRT breast cancer treatments has been demonstrated in-silico. This long awaited tool will greatly improve the treatment safety and efficacy, overcoming the limits identified within the clinical trials carried out so far.

New Advantage in Stereotactic Treatment of Lung and Pancreatic Cancer. Performance of Ultra-High Energy Electron (VHEE) Therapy Adjuvanted to the FLASH Effect: Clinical Implications and Treatment Plans Analysis

PURPOSE/OBJECTIVE(S)

Very High-Energy Electron (VHEE) beams delivered at ultra-high dose rates and hence profiting from the FLASH effect, may be a viable alternative to conventional treatment plans for the treatment of deep-seated tumors. Experimental data support the evidence of a considerable normal tissue sparing effect when treatments are delivered with dose rates much larger (100 times or more) compared to the conventional ones. Lung cancer and pancreatic cancer are considered the two biggest cancer killers. We urgently need more research in these areas, more awareness that support improvement in treatment strategies. To evaluate the potential of FLASH VHEE irradiation in these two clinical situations, we investigated the achievable sparing of healthy tissues and critical dose-limiting structures, with the goal of performing a higher dose prescription.

MATERIALS/METHODS

The study on the potential of VHEE for the stereotactic treatment of pancreatic and lung lesions was carried out on two clinical cases treated with Volumetric Modulated Arc Therapy (VMAT) techniques at University Hospital Campus Bio-Medico of Rome. The Planning Target Volume (PTV) was identified and the constraints on the Organs at Risk (OAR) and details on the irradiation approach were defined. The VHEE plan was designed to optimize the dose delivery to best activate the modelled FLASH effect based on the current experimental knowledge. In particular the impact on a dose threshold to activate the effect was studied. The VHEE treatment plan was based on an accurate Monte Carlo simulation of the electrons interactions and the results achievable with different FLASH effect models were studied. The simulation allowed the estimation of dose maps, which were used as input to an optimization algorithm that modified the fluence of each beam to meet treatment prescriptions in terms of dose to PTV. At the end the VHEE DVH plans were compared to VMAT plans.

RESULTS

The results demonstrated that FLASH therapy with VHEE beams of 70-130 MeV, could represent a promising alternative to standard radiotherapy allowing a comparable sparing of the healthy tissues. In the case of pancreatic cancer, the Dose Volume Histograms (DVH) showed how such a technique can be effective in sparing the duodenum. In case of lung cancers, the result showed how pulmonary tissue sparing can lead to a substantial reduction of pulmonary toxicity in comparison with the VMAT technique.

CONCLUSION

In the case of pancreatic cancer and assuming a non-negligible contribution from the FLASH effect, the DVH showed how the duodenum healthy tissue sparing could allow a higher dose to be prescribed at the target while keeping the constraints respected, improving the therapeutic ratio. In the case of lung cancer, the advantages of the technique are additionally increased by the significant benefit that could be related to the treatment delivery time reduction (<1s) and to the corresponding advantage coming from a reduced organ movement that translates in a lower risk of lung toxicity.

GPU-accelerated Monte Carlo simulation of electron and photon interactions for radiotherapy applications

Objective. The Monte Carlo simulation software is a valuable tool in radiation therapy, in particular to achieve the needed accuracy in the dose evaluation for the treatment plans optimisation. The current challenge in this field is the time reduction to open the way to many clinical applications for which the computational time is an issue. In this manuscript we present an innovative GPU-accelerated Monte Carlo software for dose valuation in electron and photon based radiotherapy, developed as an update of the FRED (Fast paRticle thErapy Dose evaluator) software.Approach. The code transports particles through a 3D voxel grid, while scoring their energy deposition along their trajectory. The models of electromagnetic interactions in the energy region between 1 MeV-1 GeV available in literature have been implemented to efficiently run on GPUs, allowing to combine a fast tracking while keeping high accuracy in dose assessment. The FRED software has been bench-marked against state-of-art full MC (FLUKA, GEANT4) in the realm of two different radiotherapy applications: Intra-Operative Radio Therapy and Very High Electron Energy radiotherapy applications.Results. The single pencil beam dose-depth profiles in water as well as the dose map computed on non-homogeneous phantom agree with full-MCs at 2% level, observing a gain in processing time from 200 to 5000.Significance. Such performance allows for computing a plan with electron beams in few minutes with an accuracy of ∼%, demonstrating the FRED potential to be adopted for fast plan re-calculation in photon or electron radiotherapy applications.

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.

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.

Review and performance of the Dose Profiler, a particle therapy treatments online monitor

Particle therapy (PT) can exploit heavy ions (such as He, C or O) to enhance the treatment efficacy, profiting from the increased Relative Biological Effectiveness and Oxygen Enhancement Ratio of these projectiles with respect to proton beams. To maximise the gain in tumor control probability a precise online monitoring of the dose release is needed, avoiding unnecessary large safety margins surroundings the tumor volume accounting for possible patient mispositioning or morphological changes with respect to the initial CT scan. The Dose Profiler (DP) detector, presented in this manuscript, is a scintillating fibres tracker of charged secondary particles (mainly protons) that will be operating during the treatment, allowing for an online range monitoring. Such monitoring technique is particularly promising in the context of heavy ions PT, in which the precision achievable by other techniques based on secondary photons detection is limited by the environmental background during the beam delivery. Developed and built at the SBAI department of "La Sapienza", within the INSIDE collaboration and as part of a Centro Fermi flagship project, the DP is a tracker detector specifically designed and planned for clinical applications inside a PT treatment room. The DP operation in clinical like conditions has been tested with the proton and carbon ions beams of Trento proton-therapy center and of the CNAO facility. In this contribution the detector performances are presented, in the context of the carbon ions monitoring clinical trial that is about to start at the CNAO centre.

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.

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.

A novel radioguided surgery technique exploiting β(-) decays

The background induced by the high penetration power of the radiation is the main limiting factor of the current radio-guided surgery (RGS). To partially mitigate it, a RGS with β⁺-emitting radio-tracers has been suggested in literature. Here we propose the use of β⁻-emitting radio-tracers and β⁻ probes and discuss the advantage of this method with respect to the previously explored ones: the electron low penetration power allows for simple and versatile probes and could extend RGS to tumours for which background originating from nearby healthy tissue makes probes less effective. We developed a β⁻ probe prototype and studied its performances on phantoms. By means of a detailed simulation we have also extrapolated the results to estimate the performances in a realistic case of meningioma, pathology which is going to be our first in-vivo test case. A good sensitivity to residuals down to 0.1 ml can be reached within 1 s with an administered activity smaller than those for PET-scans thus making the radiation exposure to medical personnel negligible.

Charged particle's flux measurement from PMMA irradiated by 80 MeV/u carbon ion beam

Hadrontherapy is an emerging technique in cancer therapy that uses beams of charged particles. To meet the improved capability of hadrontherapy in matching the dose release with the cancer position, new dose-monitoring techniques need to be developed and introduced into clinical use. The measurement of the fluxes of the secondary particles produced by the hadron beam is of fundamental importance in the design of any dose-monitoring device and is eagerly needed to tune Monte Carlo simulations. We report the measurements carried out with charged secondary particles produced from the interaction of a 80 MeV/u fully stripped carbon ion beam at the INFN Laboratori Nazionali del Sud, Catania, with a poly-methyl methacrylate target. Charged secondary particles, produced at 90° with respect to the beam axis, have been tracked with a drift chamber, while their energy and time of flight have been measured by means of a LYSO scintillator. Secondary protons have been identified exploiting the energy and time-of-flight information, and their emission region has been reconstructed backtracking from the drift chamber to the target. Moreover, a position scan of the target indicates that the reconstructed emission region follows the movement of the expected Bragg peak position. Exploiting the reconstruction of the emission region, an accuracy on the Bragg peak determination in the submillimeter range has been obtained. The measured differential production rate for protons produced with E(Prod)(kin) > 83 MeV and emitted at 90° with respect to the beam line is dN(P)/(dN(C)dΩ) (E(Prod)(kin) > 83 MeV, θ = 90°) = (2.69 ± 0.08(stat) ± 0.12(sys)) × 10⁻⁴ sr⁻¹.