A machine learning-based framework for delivery error prediction in proton pencil beam scanning using irradiation log-files

PurposeThis study aims to investigate the use of machine learning models for delivery error prediction in proton pencil beam scanning (PBS) delivery.MethodsA dataset of planned and delivered PBS spot parameters was generated from a set of 20 prostate patient treatments. Planned spot parameters (spot position, MU and energy) were extracted from the treatment planning system (TPS) for each beam. Delivered spot parameters were extracted from irradiation log-files for each beam delivery following treatment. The dataset was used as a training dataset for three machine learning models which were trained to predict delivered spot parameters based on planned parameters. K-fold cross validation was employed for hyper-parameter tuning and model selection where the mean absolute error (MAE) was used as the model evaluation metric. The model with lowest MAE was then selected to generate a predicted dose distribution for a test prostate patient within a commercial TPS.ResultsAnalysis of the spot position delivery error between planned and delivered values resulted in standard deviations of 0.39 mm and 0.44 mm for x and y spot positions respectively. Prediction error standard deviation values of spot positions using the selected model were 0.22 mm and 0.11 mm for x and y spot positions respectively. Finally, a three-way comparison of dose distributions and DVH values for select OARs indicates that the random-forest-predicted dose distribution within the test prostate patient was in closer agreement to the delivered dose distribution than the planned distribution.ConclusionsPBS delivery error can be accurately predicted using machine learning techniques.     

Ultrasound tracking for intra-fractional motion compensation in radiation therapy

Modern techniques as ion beam therapy or 4D imaging require precise target position information. However, target motion particularly in the abdomen due to respiration or patient movement is still a challenge and demands methods that detect and compensate this motion. Ultrasound represents a non-invasive, dose-free and model-independent alternative to fluoroscopy, respiration belt or optical tracking of the patient surface. Thus, ultrasound based motion tracking was integrated into irradiation with actively scanned heavy ions. In a first in vitro experiment, the ultrasound tracking system was used to compensate diverse sinusoidal target motions in two dimensions. A time delay of ∼200 ms between target motion and reported position data was compensated by a prediction algorithm (artificial neural network). The irradiated films proved feasibility of the proposed method. Furthermore, a practicable and reliable calibration workflow was developed to enable the transformation of ultrasound tracking data to the coordinates of the treatment delivery or imaging system – even if the ultrasound probe moves due to respiration. A first proof of principle experiment was performed during time-resolved positron emission tomography (4DPET) to test the calibration workflow and to show the accuracy of an ultrasound based motion tracking in vitro. The results showed that optical ultrasound tracking can reach acceptable accuracies and encourage further research.     

Beam characterisation studies of the 62 MeV proton therapy beamline at the Clatterbridge Cancer Centre

The Clatterbridge Cancer Centre (CCC) in the United Kingdom is the world’s first hospital proton beam therapy facility, providing treatment for ocular cancers since 1989. A 62 MeV beam of protons is produced by a Scanditronix cyclotron and transported through a passive delivery system. In addition to the long history of clinical use, the facility supports a wide programme of experimental work and as such, an accurate and reliable simulation model of the treatment beamline is highly valuable. However, as the facility has seen several changes to the accelerator and beamline over the years, a comprehensive study of the CCC beam dynamics is needed to firstly examine the beam optics. An extensive analysis was required to overcome facility related constraints to determine fundamental beamline parameters and define an optical lattice written with the Methodical Accelerator Design (MAD-X) and the particle tracking Beam Delivery Simulation (BDSIM) code. An optimised case is presented and simulated results of the optical functions, beam distribution, losses and the transverse rms beam sizes along the beamline are discussed. Corresponding optical and beam information was used in TOPAS to simulate transverse beam profiles and compared to EBT3 film measurements. We provide an overview of the magnetic components, beam transport, cyclotron, beam and treatment related parameters necessary for the development of a present day optical model of the facility. This work represents the first comprehensive study of the CCC facility to date, as a basis to determine input beam parameters to accurately simulate and completely characterise the beamline.     

Geometric and dosimetric accuracy and imaging dose of the real-time tumour tracking system of a gimbal mounted linac

PurposeTo suggest a comprehensive testing scheme to evaluate the geometric and dosimetric accuracy and the imaging dose of the VERO dynamic tumour tracking (DTT) for its clinical implementation.MethodsGeometric accuracy was evaluated for gantry 0° and 90° in terms of prediction (E_P), mechanical (E_M) and tracking (E_T) errors for sinusoidal patterns with 10 and 20 mm amplitudes, 2–6 s periods and phase shift up to 1 s and for 3 patient patterns. The automatic 4D model update was investigated simulating changes in the breathing pattern during treatment.Dosimetric accuracy was evaluated with gafchromic films irradiated in static and moving phantom with and without DTT. The entrance skin dose (ESD) was assessed using a solid state detector and gafchromic films.ResultsThe RMS of E_P, E_M, and E_T were up to 0.8, 0.5 and 0.9 mm for all non phased-shifted motion patterns while for the phased-shifted ones, E_P and E_T increased to 2.2 and 2.6 mm. Up to 4 updates are necessary to restore a good correlation model, according to type of change.For 100 kVp and 1 mA s X-ray beam, the ESD per portal due to 20 s fluoroscopy was 16.6 mGy, while treatment verification at a frequency of 1 Hz contributed with 4.2 mGy/min.ConclusionsThe proposed testing scheme highlighted that the VERO DTT system tracks a moving target with high accuracy. The automatic update of the 4D model is a powerful tool to guarantee the accuracy of tracking without increasing the imaging dose.     

Commissioning of the 4-D treatment delivery system for organ motion management in synchrotron-based scanning ion beams

PurposeThe aim of this work was the commissioning of delivery procedures for the treatment of moving targets in scanning pencil beam hadrontherapy.MethodsEBT3 films fixed to the Anzai Respiratory Phantom were exposed to carbon ion scanned homogeneous fields (E = 332 MeV/u). To evaluate the interplay effect, field size and flatness for 3 different scenarios were compared to static condition: gated irradiation or repainting alone and combination of both. Respiratory signal was provided by Anzai pressure sensor or optical tracking system (OTS). End-exhale phase and 1 s gating window were chosen (2.5 mm residual motion). Dose measurements were performed using a PinPoint ionization chamber inserted into the Brainlab ET Gating Phantom. A sub-set of tests was also performed using proton beams.ResultsThe combination of gating technique and repainting (N = 5) showed excellent results (6.1% vs 4.3% flatness, identical field size and dose deviation within 1.3%). Treatment delivery time was acceptable. Dose homogeneity for gated irradiation alone was poor. Both Anzai sensor and OTS appeared suitable for providing respiratory signal. Comparisons between protons and carbon ions showed that larger beam spot sizes represent more favorable condition for minimizing motion effect.ConclusionResults of measurements performed on different phantoms showed that the combination of gating and layered repainting is suitable to treat moving targets using scanning ion beams. Abdominal compression using thermoplastic masks, together with multi-field planning approach and multi-fractionation, have also been assessed as additional strategies to mitigate the effect of patient respiration in the clinical practice.     

Evaluation of scanning resolution,detector choice and detector orientation to be used for accurate and time-efficient commissioning of a 6MV clinical linear accelerator

The present study is aimed at exploring different scanning parameters, detectors and their orientations for time-efficient and accurate commissioning of a 6 MV clinical linear accelerator (LINAC). Beam profiles and percentage depth dose (PDD) curves were measured with a PTW dosimetry diode, a PTW Semiflex and a PinPoint ion chamber in different orientations. To acquire beam data, equidistant (step size of 0.5 mm, 1 mm, 2 mm and 3 mm) and fanline (step size of 2–0.5 mm, 2–1 mm, 3–0.5 mm and 3-1 mm) scanning modes were employed and data measurement time was recorded. Scan time per measurement point was also varied (0.2 s, 0.5 s and 1.0 s) to investigate its effect on the accuracy and acquisition time of beam data. Accuracy of the measured data was analyzed on the basis of the variation between measured data and data modeled by a treatment planning system. Beam profiles (particularly in penumbra region) were found to be sensitive to variation in scanning resolution and showed an improved accuracy with decrease in step size, while PDD curves were affected negligibly. The accuracy of beam data obtained with the PTW dosimetry diode and the PinPoint ion chamber was higher than those obtained with the PTW Semiflex ion chamber for small fields (2?×?2 cm² and 3?×?3 cm²). However, the response of the PTW diode and the PinPoint ion chamber was significantly indifferent in these fields. Furthermore, axial orientation of the PTW Semiflex ion chamber improved accuracy of profiles and PDDs as compared to radial orientation, while such a difference was not significant for the PinPoint ion chamber. It is concluded that a scan time of 0.2 s/point with a fanline scanning resolution of 2–1 mm for beam profiles and 3 mm for PDDs are most favorable in terms of accuracy and time efficiency. For small fields (2?×?2 cm² and 3?×?3 cm²), a PinPoint ion chamber in radial orientation or a dosimetry diode in axial orientation are recommended for both beam profiles and PDDs. If a PinPoint ion chamber and a PTW dosimetry diode are not available, a Semiflex ion chamber in axial orientation may be used for small fields.

     

3D Orbital Tracking in a Modified Two-photon Microscope: An Application to the Tracking of Intracellular Vesicles

The objective of this video protocol is to discuss how to perform and analyze a three-dimensional fluorescent orbital particle tracking experiment using a modified two-photon microscope¹. As opposed to conventional approaches (raster scan or wide field based on a stack of frames), the 3D orbital tracking allows to localize and follow with a high spatial (10 nm accuracy) and temporal resolution (50 Hz frequency response) the 3D displacement of a moving fluorescent particle on length-scales of hundreds of microns². The method is based on a feedback algorithm that controls the hardware of a two-photon laser scanning microscope in order to perform a circular orbit around the object to be tracked: the feedback mechanism will maintain the fluorescent object in the center by controlling the displacement of the scanning beam^3-5. To demonstrate the advantages of this technique, we followed a fast moving organelle, the lysosome, within a living cell^6,7. Cells were plated according to standard protocols, and stained using a commercially lysosome dye. We discuss briefly the hardware configuration and in more detail the control software, to perform a 3D orbital tracking experiment inside living cells. We discuss in detail the parameters required in order to control the scanning microscope and enable the motion of the beam in a closed orbit around the particle. We conclude by demonstrating how this method can be effectively used to track the fast motion of a labeled lysosome along microtubules in 3D within a live cell. Lysosomes can move with speeds in the range of 0.4-0.5 µm/sec, typically displaying a directed motion along the microtubule network⁸.     

Experimental benchmarking of RayStation proton dose calculation algorithms inside and outside the target region in heterogeneous phantom geometries

PurposeThe aim of the presented study was to complement existing literature on benchmarking proton dose by comparing dose calculations with experimental measurements in heterogeneous phantom. Points of interest inside and outside the target were considered to quantify the magnitude of calculation uncertainties in current and previous proton therapy practice that might especially have an impact on the dose in organs at risk (OARs).MethodsThe RayStation treatment planning system (RaySearch Laboratories), offering two dose calculation algorithms for pencil beam scanning in proton therapy, i.e., Pencil Beam (PB) and Monte Carlo (MC), was utilized. Treatment plans for a target located behind the interface of the heterogeneous tissues were generated. Dose measurements within and behind the target were performed in a water phantom with embedded slabs of various tissue equivalent materials and 24 PinPoint ionization chambers (PTW). In total 12 test configurations encompassing two different target depths, oblique beam incidence of 30 degrees and range shifter, were considered.ResultsPB and MC calculated doses agreed equally well with the measurements for all test geometries within the target, including the range shifter (mean dose differences ± 3%). Outside the target, the maximum dose difference of 9% (19%) was observed for MC (PB) for the oblique beam incidence and inserted range shifter.ConclusionThe accuracy of MC dose algorithm was superior compared to the PB algorithm, especially outside the target volumes. MC based dose calculation should therefore be preferred in treatment scenarios with heterogeneities, especially to reduce clinically relevant uncertainties for OARs.     

First-in-human imaging using a MR-compatible e4D ultrasound probe for motion management of radiotherapy

PurposeRespiration-induced tumor or organ positional changes can impact the accuracy of external beam radiotherapy. Motion management strategies are used to account for these changes during treatment. The authors report on the development, testing, and first-in-human evaluation of an electronic 4D (e4D) MR-compatible ultrasound probe that was designed for hands-free operation in a MR and linear accelerator (LINAC) environment.MethodsUltrasound components were evaluated for MR compatibility. Electromagnetic interference (EMI) shielding was used to enclose the entire probe and a factory-fabricated cable shielded with copper braids was integrated into the probe. A series of simultaneous ultrasound and MR scans were acquired and analyzed in five healthy volunteers.ResultsThe ultrasound probe led to minor susceptibility artifacts in the MR images immediately proximal to the ultrasound probe at a depth of <10 mm. Ultrasound and MR-based motion traces that were derived by tracking the salient motion of endogenous target structures in the superior-inferior (SI) direction demonstrated good concordance (Pearson correlation coefficients of 0.95–0.98) between the ultrasound and MRI datasets.ConclusionWe have demonstrated that our hands-free, e4D probe can acquire ultrasound images during a MR acquisition at frame rates of approximately 4 frames per second (fps) without impacting either the MR or ultrasound image quality. This use of this technology for interventional procedures (e.g. biopsies and drug delivery) and motion compensation during imaging are also being explored.     

Clinical implementations of 4D pencil beam scanned particle therapy: Report on the 4D treatment planning workshop 2016 and 2017

In 2016 and 2017, the 8th and 9th 4D treatment planning workshop took place in Groningen (the Netherlands) and Vienna (Austria), respectively. This annual workshop brings together international experts to discuss research, advances in clinical implementation as well as problems and challenges in 4D treatment planning, mainly in spot scanned proton therapy. In the last two years several aspects like treatment planning, beam delivery, Monte Carlo simulations, motion modeling and monitoring, QA phantoms as well as 4D imaging were thoroughly discussed.This report provides an overview of discussed topics, recent findings and literature review from the last two years. Its main focus is to highlight translation of 4D research into clinical practice and to discuss remaining challenges and pitfalls that still need to be addressed and to be overcome.     

Scan path optimization with/without clustering for active beam delivery in charged particle therapy

PurposeTo compare different algorithms to optimize the scanning path in charged particle therapy with quasidiscrete scanning. We implemented a Hybrid Genetic Algorithm with Heuristics (HyGA) and combined it with clustering techniques. The performance was compared to Simulated Annealing (SA) and to commercially available treatment planning system (TPS).MethodsPerformance and clinical implications were assessed using data from 10 patients treated at CNAO (Centro Nazionale di Adroterapia Oncologica). Clinical treatments are performed relying on beam deflection, avoiding irradiation for transitions between adjacent spots larger than 2 cm. A clustering method was implemented with HyGA (HyGA_Cl), which assumes beam deflection during transition between clusters. Clinical performance was determined as the total number of particles delivered during spot transitions and the number of particles wasted due to beam deflection. Results were compared to scan paths obtained with CNAO TPS.ResultsSA and HyGA produced on average shorter paths compared to the currently available TPS. This did not result in a reduction of transit particles, due to the concomitant effect of beam deflection out of the extraction line. HyGA_Cl achieved 2% average reduction in transit particles when compared to CNAO TPS. As a drawback, wasted particles increased, due to more frequent use of beam deflection. Both the SA and HyGA algorithms reduced the number of wasted particles.ConclusionSA and HyGA proved to be the most cost-effective methods in reducing wasted particles, with benefits in terms of shorter scan paths. A decrease in transit particles delivered with beam deflection can be achieved using HyGA_Cl.     

Perturbation analysis of 4D dose distribution for scanned carbon-ion beam radiotherapy

PurposeTo evaluate the patients’ set-up error-induced perturbation effects on 4D dose distributions (4DDD) of range-adapted internal target volume-based (raITV) treatment plan using lung and liver 4DCT data sets.MethodsWe enrolled 20 patients with lung and liver cancer treated with respiratory-gated carbon-ion beam scanning therapy. PTVs were generated by adding a 2 mm range-adapted set-up margin on the raITVs. Set-up errors were simulated by shifting the beam isocenter in three translational directions of ±2 mm, ±4 mm, and ±6 mm. 4DDDs were calculated for both nominal and isocenter-shifted situations. Dose metrics of CTV dose coverage (D95) and normal tissue sparing were evaluated. Statistical significance with p < 0.01 was considered by Wilcoxon signed rank test.ResultsThe CTV dose coverage was more sensitive to set-up errors for lung cases than for liver cases, and more serious in superior-inferior direction. The sufficient CTV-D95 > 98% could be achieved with set-up errors less than ±2 mm in all shift directions both for lung and liver cases. With the increase of set-up error, the CTV dose coverage decreased gradually. The clinical criterial of CTV-D95 > 95% could not be fulfilled with set-up error reached to ±4 mm for lung cases, and ±6 mm for liver cases. OAR doses did not have a significant difference with each set-up error for both lung and liver cases.ConclusionsThe range-adapted set-up margin successfully prevented dose degradation of 4DDDs in the presence of the same magnitude of set-up error for raITV-based carbon-ion beam scanning therapy.     

Motion management with variable cycle-based respiratory guidance method for carbon-ion pencil beam scanning treatment

PurposeA novel variable cycle-based respiratory guidance method was proposed to synchronize the patterns between patients’ breathing and the magnetic excitation of synchrotron under the mode of full-energy depth scanning beam delivery, in order to improve the treatment precision and efficiency for carbon ion therapy.MethodsAudio-visual biofeedback system with variable cycle-based respiratory guidance method was developed. We enrolled 6 healthy volunteers and a simulation study of the fixed cycle-based and variable cycle-based respiratory guidance with three treatment fractions was performed. A total of 72 breathing curves were collected for 4D dose calculations with three 4DCT datasets of lung tumor cases. Target dose coverage (D95: the percent dose covering 95% of the target), dose homogeneity (D5-D95), and treatment time were analyzed. The Wilcoxon signed-rank test was used for statistical difference analysis, and p < 0.05 was considered significant.ResultsWith the variable cycle-based respiratory guidance method, the breath hold phase of breathing curve could be synchronized with the synchrotron flat-top phase over time. The dose homogeneity was improved by factors of 1.94–2.92 compared to the fixed cycle-based respiratory guidance maneuvers alone or in combination with gating technique. Moreover, the treatment efficiency increased by 11–23%, depending on the duty cycle settings of the gating window.ConclusionsThe proposed variable cycle-based respiratory guidance method could improve both the treatment efficiency and precision under the mode of the full-energy depth scanning beam delivery for synchrotron-based carbon ion therapy.     

A patient immobilization device for prone breast radiotherapy: Dosimetric effects and inclusion in the treatment planning system

PurposeTo assess the dosimetric impact of a patient positioning device for prone breast radiotherapy and assess the accuracy of a treatment planning system (TPS) in predicting this impact.MethodsBeam attenuation and build-up dose perturbations, quantified by ionization chamber and radiochromic film dosimetry, were evaluated for 3 components of the patient positioning device: the carbon fiber baseplate, the support cushions and the support wedge for the contralateral breast. Dose calculations were performed using the XVMC dose engine implemented in the Monaco TPS. All components were included during planning CT acquisition.ResultsBeam attenuation amounted to 7.57% (6 MV) and 5.33% (15 MV) for beams obliquely intersecting the couchtop–baseplate combination. Beams traversing large sections of the support wedge were attenuated by 12.28% (6 MV) and 9.37% (15 MV). For the support cushion foam, beam attenuation remained limited to 0.11% (6 MV) and 0.08% (15 MV) per centimeter thickness. A substantial loss of dose build-up was detected when irradiating through any of the investigated components. TPS dose calculations accurately predicted beam attenuation by the baseplate and support wedge. A manual density overwrite was needed to model attenuation by the support cushion foam. TPS dose calculations in build-up regions differed considerably from measurements for both open beams and beams traversing the device components.ConclusionsIrradiating through the components of the positioning device resulted in a considerable degradation of skin sparing. Inclusion of the device components in the treatment planning CT allowed to accurately model the most important attenuation effect, but failed to accurately predict build-up doses.     

Retrospective assessment of a single fiducial marker tracking regimen with robotic stereotactic body radiation therapy for liver tumours

AimTo investigate tumour motion tracking uncertainties in the CyberKnife Synchrony system with single fiducial marker in liver tumours.BackgroundIn the fiducial-based CyberKnife real-time tumour motion tracking system, multiple fiducial markers are generally used to enable translation and rotation corrections during tracking. However, sometimes a single fiducial marker is employed when rotation corrections are not estimated during treatment.Materials and methodsData were analysed for 32 patients with liver tumours where one fiducial marker was implanted. Four-dimensional computed tomography (CT) scans were performed to determine the internal target volume (ITV). Before the first treatment fraction, the CT scans were repeated and the marker migration was determined. Log files generated by the Synchrony system were obtained after each treatment and the correlation model errors were calculated. Intra-fractional spine rotations were examined on the spine alignment images before and after each treatment.ResultsThe mean (standard deviation) ITV margin was 4.1 (2.3) mm, which correlated weakly with the distance between the fiducial marker and the tumour. The mean migration distance of the marker was 1.5 (0.7) mm. The overall mean correlation model error was 1.03 (0.37) mm in the radial direction. The overall mean spine rotations were 0.27° (0.31), 0.25° (0.22), and 0.23° (0.26) for roll, pitch, and yaw, respectively. The treatment time was moderately associated with the correlation model errors and weakly related to spine rotation in the roll and yaw planes.ConclusionsMore caution and an additional safety margins are required when tracking a single fiducial marker.     

A model of time estimation and error feedback in predictive timing behavior

Two key features of sensorimotor prediction are preprogramming and adjusting of performance based on previous experience. Oculomotor tracking of alternating visual targets provides a simple paradigm to study this behavior in the motor system; subjects make predictive eye movements (saccades) at fast target pacing rates (>0.5 Hz). In addition, the initiation errors (latencies) during predictive tracking are correlated over a small temporal window (correlation window) suggesting that tracking performance within this time range is used in the feedback process of the timing behavior. In this paper, we propose a closed-loop model of this predictive timing. In this model, the timing between movements is based on an internal estimation of stimulus timing (an internal clock), which is represented by a (noisy) signal integrated to a threshold. The threshold of the integrate-to-fire mechanism is determined by the timing between movements made within the correlation window of previous performance and adjusted by feedback of recent and projected initiation error. The correlation window size increases with repeated tracking and was estimated by two independent experiments. We apply the model to several experimental paradigms and show that it produces data specific to predictive tracking: a gradual shift from reaction to prediction on initial tracking, phase transition and hysteresis as pacing frequency changes, scalar property, continuation of predictive tracking despite perturbations, and intertrial correlations of a specific form. These results suggest that the process underlying repetitive predictive motor timing is adjusted by the performance and the corresponding errors accrued over a limited time range and that this range increases with continued confidence in previous performance.     

Four-dimensional carbon-ion pencil beam treatment planning comparison between robust optimization and range-adapted internal target volume for respiratory-gated liver and lung treatment

We investigated the dose differences between robust optimization-based treatment planning (4DRO) and range-adapted internal target volume (rITV). We used 4DCT dataset of 20 lung cancer and 20 liver cancer patients, respectively, who had been treated with respiratory-gated carbon-ion pencil beam scanning therapy. 4DRO and rITV plans were created with the same clinical target volume (CTV) and organs at risk (OAR) contours. Four-dimensional dose distribution was calculated using deformable image registration. Dose metrics (e.g. D95, V20) were analyzed. Statistical significance was assessed by the Wilcoxon signed-rank test. For the lung cases, the mean CTV-D95 value for the rITV plan (=98.5%) was same as that for the 4DRO plan (=98.5%, P = 0.106), while the mean D95 value for the CTV + setup margin contour for the rITV plan (=98.2%) was higher than that for the 4DRO plan (95.2%, P < 0.001). For the liver cases, the mean CTV-D95 value for the rITV plan (=98.1%) was slightly lower than that for the 4DRO plan (=98.5%, P < 0.01), while the mean D95 value for the CTV + setup margin contour for the rITV plan (=98.0%) was higher than that for the 4DRO plan (94.1%, P < 0.001). For the doses to the organs at risk (OARs), the ipsilateral lung-V20/liver-V20 values for the rITV plan (=10.1%/19.7%) was significantly higher than that for the 4DRO plan (=8.6%/17.6, P < 0.001). Although the target coverage for 4DRO plan may be worse than that for rITV plan in the presence of the setup error, the 4DRO plan can improve OAR dose while preserving acceptable target dose coverage.     

A GATE/Geant4 beam model for the MedAustron non-isocentric proton treatment plans quality assurance

PurposeTo present a reference Monte Carlo (MC) beam model developed in GATE/Geant4 for the MedAustron fixed beam line. The proposed model includes an absolute dose calibration in Dose-Area-Product (DAP) and it has been validated within clinical tolerances for non-isocentric treatments as routinely performed at MedAustron.Material and MethodsThe proton beam model was parametrized at the nozzle entrance considering optic and energy properties of the pencil beam. The calibration in terms of absorbed dose to water was performed exploiting the relationship between number of particles and DAP by mean of a recent formalism. Typical longitudinal dose distribution parameters (range, distal penumbra and modulation) and transverse dose distribution parameters (spot sizes, field sizes and lateral penumbra) were evaluated. The model was validated in water, considering regular-shaped dose distribution as well as clinical plans delivered in non-isocentric conditions.ResultsSimulated parameters agree with measurements within the clinical requirements at different air gaps. The agreement of distal and longitudinal dose distribution parameters is mostly better than 1 mm. The dose difference in reference conditions and for 3D dose delivery in water is within 0.5% and 1.2%, respectively. Clinical plans were reproduced within 3%.ConclusionA full nozzle beam model for active scanning proton pencil beam is described using GATE/Geant4. Absolute dose calibration based on DAP formalism was implemented. The beam model is fully validated in water over a wide range of clinical scenarios and will be inserted as a reference tool for research and for independent dose calculation in the clinical routine.     

A modular dose delivery system for treating moving targets with scanned ion beams: Performance and safety characteristics,and preliminary tests

PurposeThe purpose of this study was to develop a modular dose-delivery system (DDS) for scanned-ion radiotherapy that mitigates against organ motion artifacts by synchronizing the motion of the beam with that of the moving anatomy.MethodsWe integrated a new motion synchronization system and an existing DDS into two centers. The modular approach to integration utilized an adaptive layer of software and hardware interfaces. The method of synchronization comprised three major tasks, namely, the creation of 3D treatment plans (each representing one phase of respiratory motion and together comprising a 4D plan), monitoring anatomic motion during treatment, and synchronization of the beam to anatomic motion. The synchronization was accomplished in real time by repeatedly selecting and delivering a 3D plan, i.e., the one that most closely corresponded to the current anatomic state, until all plans were delivered. The performance characteristics of the motion mitigation system were tested by delivering 4D treatment plans to a moving phantom and comparing planned and measured dose distributions. Dosimetric performance was considered acceptable when the gamma-index pass rate was >90%, homogeneity-index value was >95%, and conformity-index value was >60%. Selected safety characteristics were tested by introducing errors during treatment and testing DDS response.ResultsAcceptable dosimetric performance and safety characteristics were observed for all treatment plans.ConclusionsWe demonstrated, for the first time, that a modular prototype system, synchronizing scanned ion beams with moving targets can deliver conformal, motion-compensated dose distributions. The prototype system was implemented and characterized at GSI and CNAO.     

Real-time tumor tracking with an artificial neural networks-based method: A feasibility study

The purpose of this study was to develop and assess the performance of a tumor tracking method designed for application in radiation therapy. This motion compensation strategy is currently applied clinically only in conventional photon radiotherapy but not in particle therapy, as greater accuracy in dose delivery is required.We proposed a tracking method that exploits artificial neural networks to estimate the internal tumor trajectory as a function of external surrogate signals. The developed algorithm was tested by means of a retrospective clinical data analysis in 20 patients, who were treated with state of the art infra-red motion tracking for photon radiotherapy, which is used as a benchmark. Integration into a hardware platform for motion tracking in particle therapy was performed and then tested on a moving phantom, specifically developed for this purpose.Clinical data show that a median tracking error reduction up to 0.7 mm can be achieved with respect to state of the art technologies. The phantom study demonstrates that a real-time tumor position estimation is feasible when the external signals are acquired at 60 Hz.The results of this work show that neural networks can be considered a valuable tool for the implementation of high accuracy real-time tumor tracking methodologies.