Robust treatment planning with 4D intensity modulated carbon ion therapy for multiple targets in stage IV non-small cell lung cancer

Investigation of factors related to treatment planning of x-ray SBRT and scanning carbon-ion radiation therapy for early-stage lung cancer patients

This study aimed to compare the treatment plans of x-ray SBRT and scanning carbon ion radiation therapy (CIRT) for localized lung tumors, and to evaluate the dose dependence of tumor size tumor-to-heart distance. For phantom verification, we used a chest phantom with a spherical simulated tumor. Treatment plans for 3-dimensional conformal radiation therapy (3D-CRT), volumetric modulated arc therapy (VMAT), and CIRT were created. GTVs were created in sizes ranging from 0.5 to 5 cm in diameter, and the dependence of the lung dose on GTV diameter was evaluated for each treatment plan. For patient validation, 30 cases of localized lung tumors were analyzed. 3D-CRT, VMAT, and CIRT treatment plans were developed, and DVH parameters were evaluated for each GTV size and GTV-to-heart distance. In both phantom and patient validations, the OAR doses were the lowest for CIRT. The lung dose increased with increasing GTV diameter for all three treatment plans. CIRT had the smallest ratio of lung dose increase to GTV diameter increase among the three treatment plans. Heart dose in CIRT was independent of GTV size and GTV-to-heart distance Conclusions: The results of the present study suggested that the use of scanning CIRT can reduce the OAR dose while guaranteeing the tumor dose compared to x-ray SBRT. In addition, it was suggested that CIRT can treat patients with large GTV sizes while maintaining low lung and heart dose.

Evaluation of motion mitigation strategies for carbon ion therapy of abdominal tumors based on non-periodic imaging data

Objective.To identify suitable combination strategies for treatment planning and beam delivery in scanned carbon ion therapy of moving tumors.Approach. Carbon ion treatment plans for five abdominal tumors were optimized on four-dimensional (4D) computed tomography (CT) data using the following approaches. 4DITV across all phases and within a gating window, single phase uniform dose, and an innovative 4D tracking internal target volume (ITV) strategy. Delivered single-fraction doses were calculated on time-resolved virtual CT images reconstructed from 2D cine-magnetic resonance imaging series, using a deformable image registration pipeline. Treatment plans were combined with various beam delivery techniques: three-dimensional (no motion mitigation), rescanning, gating, beam tracking, and multi-phase 4D delivery with and without residual tracking (MP4D and MP4DRT) to form in total 11 treatment modalities. Single fraction doses were accumulated to simulate a fractionated treatment.Main results. Breath-sampled treatments using the MP4D and MP4DRT delivery techniques were the only to achieveD95> 95% for hypofractionated treatments, with little dependence on the number of fractions. A combination of MP4DRT with the new 4D tracking ITV approach resulting in conformal dose distributions and demonstrated the greatest robustness against irregular motion and anatomical changes.Significance. This study demonstrates, that real-time adaptive beam delivery strategies can deliver conformal doses within single fractions, thereby enabling hypofractionated treatment schemes that are not feasible with conventional strategies.

Feasibility study of 4D-online monitoring of density gradients induced by lung cancer treatment using carbon ions

Tumor motion is a major challenge for scanned ion-beam therapy. In the case of lung tumors, strong under- and overdosage can be induced due to the high density gradients between the tumor- and bone tissues compared to lung tissues. This work proposes a non-invasive concept for 4D monitoring of high density gradients in carbon ion beam therapy, by detecting charged fragments. The method implements CMOS particle trackers that are used to reconstruct the fragment vertices, which define the emission points of nuclear interactions between the primary carbon ions and the patient tissues. A 3D treatment plan was optimized to deliver 2 Gy to a static spherical target volume. The goodness of the method was assessed by comparing reconstructed vertices measured in two static cases to the ones in a non-compensated moving case with an amplitude of 20 mm. The measurements, performed at the Marburg Ion-Beam Therapy Center (MIT), showed promising results to assess the conformity of the delivered dose. In particular to measure overshoots induced by high density gradients due to motion with 83.0 ± 1.5% and 92.0 ± 1.5% reliability based on the ground truth provided by the time-resolved motor position and depending on the considered volume and the iso-energy layers.

Opportunities and challenges of upright patient positioning in radiotherapy

Objective.Upright positioning has seen a surge in interest as a means to reduce radiotherapy (RT) cost, improve patient comfort, and, in selected cases, benefit treatment quality. In particle therapy (PT) in particular, eliminating the need for a gantry can present massive cost and facility footprint reduction. This review discusses the opportunities of upright RT in perspective of the open challenges.Approach.The clinical, technical, and workflow challenges that come with the upright posture have been extracted from an extensive literature review, and the current state of the art was collected in a synergistic perspective from photon and particle therapy. Considerations on future developments and opportunities are provided.Main results.Modern image guidance is paramount to upright RT, but it is not clear which modalities are essential to acquire in upright posture. Using upright MRI or upright CT, anatomical differences between upright/recumbent postures have been observed for nearly all body sites. Patient alignment similar to recumbent positioning was achieved in small patient/volunteer cohorts with prototype upright positioning systems. Possible clinical advantages, such as reduced breathing motion in upright position, have been reported, but limited cohort sizes prevent resilient conclusions on the treatment impact. Redesign of RT equipment for upright positioning, such as immobilization accessories for various body regions, is necessary, where several innovations were recently presented. Few clinical studies in upright PT have already reported promising outcomes for head&neck patients.Significance.With more evidence for benefits of upright RT emerging, several centers worldwide, particularly in PT, are installing upright positioning devices or have commenced upright treatment. Still, many challenges and open questions remain to be addressed to embed upright positioning firmly in the modern RT landscape. Guidelines, professionals trained in upright patient positioning, and large-scale clinical studies are required to bring upright RT to fruition.

Deep learning-based voxel sampling for particle therapy treatment planning

Objective.Scanned particle therapy often requires complex treatment plans, robust optimization, as well as treatment adaptation. Plan optimization is especially complicated for heavy ions due to the variable relative biological effectiveness. We present a novel deep-learning model to select a subset of voxels in the planning process thus reducing the planning problem size for improved computational efficiency.Approach.Using only a subset of the voxels in target and organs at risk (OARs) we produced high-quality treatment plans, but heuristic selection strategies require manual input. We designed a deep-learning model based onP-Net to obtain an optimal voxel sampling without relying on patient-specific user input. A cohort of 70 head and neck patients that received carbon ion therapy was used for model training (50), validation (10) and testing (10). For training, a total of 12 500 carbon ion plans were optimized, using a highly efficient artificial intelligence (AI) infrastructure implemented into a research treatment planning platform. A custom loss function increased sampling density in underdosed regions, while aiming to reduce the total number of voxels.Main results.On the test dataset, the number of voxels in the optimization could be reduced by 84.8% (median) at <1% median loss in plan quality. When the model was trained to reduce sampling in the target only while keeping all voxels in OARs, a median reduction up to 71.6% was achieved, with 0.5% loss in the plan quality. The optimization time was reduced by a factor of 7.5 for the total AI selection model and a factor of 3.7 for the model with only target selection.Significance.The novel deep-learning voxel sampling technique achieves a significant reduction in computational time with a negligible loss in the plan quality. The reduction in optimization time can be especially useful for future real-time adaptation strategies.

Quasi-real-time range monitoring by in-beam PET: a case for 15O

A fast and reliable range monitoring method is required to take full advantage of the high linear energy transfer provided by therapeutic ion beams like carbon and oxygen while minimizing damage to healthy tissue due to range uncertainties. Quasi-real-time range monitoring using in-beam positron emission tomography (PET) with therapeutic beams of positron-emitters of carbon and oxygen is a promising approach. The number of implanted ions and the time required for an unambiguous range verification are decisive factors for choosing a candidate isotope. An experimental study was performed at the FRS fragment-separator of GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany, to investigate the evolution of positron annihilation activity profiles during the implantation of [Formula: see text]O and [Formula: see text]O ion beams in a PMMA phantom. The positron activity profile was imaged by a dual-panel version of a Siemens Biograph mCT PET scanner. Results from a similar experiment using ion beams of carbon positron-emitters [Formula: see text]C and [Formula: see text]C performed at the same experimental setup were used for comparison. Owing to their shorter half-lives, the number of implanted ions required for a precise positron annihilation activity peak determination is lower for [Formula: see text]C compared to [Formula: see text]C and likewise for [Formula: see text]O compared to [Formula: see text]O, but their lower production cross-sections make it difficult to produce them at therapeutically relevant intensities. With a similar production cross-section and a 10 times shorter half-life than [Formula: see text]C, [Formula: see text]O provides a faster conclusive positron annihilation activity peak position determination for a lower number of implanted ions compared to [Formula: see text]C. A figure of merit formulation was developed for the quantitative comparison of therapy-relevant positron-emitting beams in the context of quasi-real-time beam monitoring. In conclusion, this study demonstrates that among the positron emitters of carbon and oxygen, [Formula: see text]O is the most feasible candidate for quasi-real-time range monitoring by in-beam PET that can be produced at therapeutically relevant intensities. Additionally, this study demonstrated that the in-flight production and separation method can produce beams of therapeutic quality, in terms of purity, energy, and energy spread.

Validation of robust radiobiological optimization algorithms based on the mixed beam model for intensity-modulated carbon-ion therapy

Currently, treatment planning systems (TPSs) that can compute the intensities of intensity-modulated carbon-ion therapy (IMCT) using scanned carbon-ion beams are limited. In the present study, the computational efficacy of the newly designed IMCT algorithms was analyzed for the first time based on the mixed beam model with respect to the physical and biological doses; moreover, the validity and effectiveness of the robust radiobiological optimization were verified. A dose calculation engine was independently generated to validate a clinical dose determined in the TPS. A biological assay was performed using the HSGc-C5 cell line to validate the calculated surviving fraction (SF). Both spot control (SC) and voxel-wise worst-case scenario (WC) algorithms were employed for robust radiobiological optimization followed by their application in a Radiation Therapy Oncology Group benchmark phantom under homogeneous and heterogeneous conditions and a clinical case for range and position errors. Importantly, for the first time, both SC and WC algorithms were implemented in the integrated TPS platform that can compute the intensities of IMCT using scanned carbon-ion beams for robust radiobiological optimization. For assessing the robustness, the difference between the maximum and minimum values of a dose-volume histogram index in the examined error scenarios was considered as a robustness index. The relative biological effectiveness (RBE) determined by the independent dose calculation engine exhibited a -0.6% difference compared with the RBE defined by the TPS at the isocenter, whereas the measured and the calculated SF were similar. Regardless of the objects, compared with the conventional IMCT, the robust radiobiological optimization enhanced the sensitivity of the examined error scenarios by up to 19% for the robustness index. The computational efficacy of the novel IMCT algorithms was verified according to the mixed beam model with respect to the physical and biological doses. The robust radiobiological optimizations lowered the impact of range and position uncertainties considerably in the examined scenarios. The robustness of the WC algorithm was more enhanced compared with that of the SC algorithm. Nevertheless, the SC algorithm can be used as an alternative to the WC IMCT algorithm with respect to the computational cost.

Evaluation of proton and carbon ion beam models in TReatment Planning for Particles 4D (TRiP4D) referring to a commercial treatment planning system

PURPOSE

To investigate the accuracy of the treatment planning system (TPS) TRiP4D in reproducing doses computed by the clinically used TPS SyngoRT.

METHODS

Proton and carbon ion beam models in TRiP4D were converted from SyngoRT. Cubic plans with different depths in a water-tank phantom (WP) and previously treated and experimentally verified patient plans from SyngoRT were recalculated in TRiP4D. The target mean dose deviation (ΔD_mean,T) and global gamma index (2%-2 mm for the absorbed dose and 3%-3mm for the RBE-weighted dose with 10% threshold) were evaluated.

RESULTS

The carbon and proton absorbed dose gamma passing rates (γ-PRs) were ≥99.93% and ΔD_mean,T smaller than -0.22%. On average, the RBE-weighted dose D_mean,T was -1.26% lower for TRiP4D than SyngoRT for cubic plans. In TRiP4D, the faster analytical 'low dose approximation' (Krämer, 2006) was used, while SyngoRT used a stochastic implementation (Krämer, 2000). The average ΔD_{mean, T} could be reduced to -0.59% when applying the same biological effect calculation algorithm. However, the dose recalculation time increased by a factor of 79-477. ΔD_mean,T variation up to -2.27% and -2.79% was observed for carbon absorbed and RBE-weighted doses in patient plans. The γ-PRs were ≥93.92% and ≥91.83% for patient plans, except for one proton beam with a range shifter (γ-PR of 64.19%).

CONCLUSION

The absorbed dose between TRiP4D and SyngoRT were identical for both proton and carbon ion plans in the WP. Compared to SyngoRT, TRiP4D underestimated the target RBE-weighted dose; however more efficient in RBE-weighted dose calculation. Large variation for proton beam with range shifter was observed. TRiP4D will be used to evaluate doses delivered to moving targets. Uncertainties inherent to the 4D-dose reconstruction calculation are expected to be significantly larger than the dose errors reported here. For this reason, the residual differences between TRiP4D and SyngoRT observed in this study are considered acceptable. The study was approved by the Institutional Research Board of Shanghai Proton and Heavy Ion Center (approval number SPHIC-MP-2020-04, RS).

Comparison of 3D and 4D robustly optimized proton treatment plans for non-small cell lung cancer patients with tumour motion amplitudes larger than 5 mm

Background and purpose

There is no consensus about an ideal robust optimization (RO) strategy for proton therapy of targets with large intrafractional motion. We investigated the plan robustness of 3D and different 4D RO strategies.

Materials and methods

For eight non-small cell lung cancer patients with clinical target volume (CTV) motion >5 mm, different RO approaches were investigated: 3DRO considering the average CT (AvgCT) with a target density override, 4DRO considering three/all 4DCT phases, and 4DRO considering the AvgCT and three/all 4DCT phases. Robustness against setup/range errors, interplay effects based on breathing and machine log file data for deliveries with/without rescanning, and interfractional anatomical changes were analyzed for target coverage and OAR sparing.

Results

All nominal plans fulfilled the clinical requirements with individual CTV coverage differences <2pp; 4DRO without AvgCT generated the most conformal dose distributions. Robustness against setup/range errors was best for 4DRO with AvgCT (18% more passed error scenarios than 3DRO). Interplay effects caused fraction-wise median CTV coverage loss of 3pp and missed maximum dose constraints for heart and esophagus in 18% of scenarios. CTV coverage and OAR sparing fulfilled requirements in all cases when accumulating four interplay scenarios. Interfractional changes caused less target misses for RO with AvgCT compared to 4DRO without AvgCT (≤42%/33% vs. ≥56%/44% failed single/accumulated scenarios).

Conclusions

All RO strategies provided acceptable plans with equally low robustness against interplay effects demanding other mitigation than rescanning to ensure fraction-wise target coverage. 4DRO considering three phases and the AvgCT provided best compromise on planning effort and robustness.

Experimental Validation of a Real-Time Adaptive 4D-Optimized Particle Radiotherapy Approach to Treat Irregularly Moving Tumors

PURPOSE

Treatment of locally advanced lung cancer is limited by toxicity and insufficient local control. Particle therapy could enable more conformal treatment than intensity modulated photon therapy but is challenged by irregular tumor motion, associated range changes, and tumor deformations. We propose a new strategy for robust, online adaptive particle therapy, synergizing 4-dimensional optimization with real-time adaptive beam tracking. The strategy was tested and the required motion monitoring precision was determined.

METHODS AND MATERIALS

In multiphase 4-dimensional dose delivery (MP4D), a dedicated quasistatic treatment plan is delivered to each motion phase of periodic 4-dimensional computed tomography (4DCT). In the new extension, "MP4D with residual tracking" (MP4DRT), lateral beam tracking compensates for the displacement of the tumor center-of-mass relative to the current phase in the planning 4DCT. We implemented this method in the dose delivery system of a clinical carbon facility and tested it experimentally for a lung cancer plan based on a periodic subset of a virtual lung 4DCT (planned motion amplitude 20 mm). Treatments were delivered in a quality assurance-like setting to a moving ionization chamber array. We considered variable motion amplitudes and baseline drifts. The required motion monitoring precision was evaluated by adding noise to the motion signal. Log-file-based dose reconstructions were performed in silico on the entire 4DCT phantom data set capable of simulating nonperiodic motion. MP4DRT was compared with MP4D, rescanned beam tracking, and internal target volume plans. Treatment quality was assessed in terms of target coverage (D95), dose homogeneity (D5-D95), conformity number, and dose to heart and lung.

RESULTS

For all considered motion scenarios and metrics, MP4DRT produced the most favorable metrics among the tested motion mitigation strategies and delivered high-quality treatments. The conformity was similar to static treatments. The motion monitoring precision required for D95 >95% was 1.9 mm.

CONCLUSIONS

With clinically feasible motion monitoring, MP4DRT can deliver highly conformal dose distributions to irregularly moving targets.

Precision of the PET activity range during irradiation with10C,11C, and12C beams

Objective. Beams of stable ions have been a well-established tool for radiotherapy for many decades. In the case of ion beam therapy with stable¹²C ions, the positron emitters^10,11C are produced via projectile and target fragmentation, and their decays enable visualization of the beam via positron emission tomography (PET). However, the PET activity peak matches the Bragg peak only roughly and PET counting statistics is low. These issues can be mitigated by using a short-lived positron emitter as a therapeutic beam.Approach.An experiment studying the precision of the measurement of ranges of positron-emitting carbon isotopes by means of PET has been performed at the FRS fragment-separator facility of GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany. The PET scanner used in the experiment is a dual-panel version of a Siemens Biograph mCT PET scanner.Main results.High-quality in-beam PET images and activity distributions have been measured from the in-flight produced positron emitting isotopes¹¹C and¹⁰C implanted into homogeneous PMMA phantoms. Taking advantage of the high statistics obtained in this experiment, we investigated the time evolution of the uncertainty of the range determined by means of PET during the course of irradiation, and show that the uncertainty improves with the inverse square root of the number of PET counts. The uncertainty is thus fully determined by the PET counting statistics. During the delivery of 1.6 × 10⁷ions in 4 spills for a total duration of 19.2 s, the PET activity range uncertainty for¹⁰C,¹¹C and¹²C is 0.04 mm, 0.7 mm and 1.3 mm, respectively. The gain in precision related to the PET counting statistics is thus much larger when going from¹¹C to¹⁰C than when going from¹²C to¹¹C. The much better precision for¹⁰C is due to its much shorter half-life, which, contrary to the case of¹¹C, also enables to include the in-spill data in the image formation.Significance. Our results can be used to estimate the contribution from PET counting statistics to the precision of range determination in a particular carbon therapy situation, taking into account the irradiation scenario, the required dose and the PET scanner characteristics.

Potential benefits of using radioactive ion beams for range margin reduction in carbon ion therapy

Sharp dose gradients and high biological effectiveness make ions such as ¹²C an ideal tool to treat deep-seated tumors, however, at the same time, sensitive to errors in the range prediction. Tumor safety margins mitigate these uncertainties, but during the irradiation they lead to unavoidable damage to the surrounding healthy tissue. To fully exploit the Bragg peak benefits, a large effort is put into establishing precise range verification methods. Despite positron emission tomography being widely in use for this purpose in ¹²C therapy, the low count rates, biological washout, and broad activity distribution still limit its precision. Instead, radioactive beams used directly for treatment would yield an improved signal and a closer match with the dose fall-off, potentially enabling precise in vivo beam range monitoring. We have performed a treatment planning study to estimate the possible impact of the reduced range uncertainties, enabled by radioactive ¹¹C ions treatments, on sparing critical organs in tumor proximity. Compared to ¹²C treatments, (i) annihilation maps for ¹¹C ions can reflect sub- millimeter shifts in dose distributions in the patient, (ii) outcomes of treatment planning with ¹¹C significantly improve and (iii) less severe toxicities for serial and parallel critical organs can be expected.

Time-resolved MRI for off-line treatment robustness evaluation in carbon-ion radiotherapy of pancreatic cancer

PURPOSE

In this study, we investigate the use of magnetic resonance imaging (MRI) for the clinical evaluation of gating treatment robustness in carbon-ion radiotherapy (CIRT) of pancreatic cancer. Indeed, MRI allows radiation-free repeated scans and fast dynamic sequences for time-resolved (TR) imaging (cine-MRI), providing information on inter- and intra-fraction cycle-to-cycle variations of respiratory motion. MRI can therefore support treatment planning and verification, overcoming the limitations of the current clinical standard, that is, four-dimensional computed tomography (4DCT), which describes an "average" breathing cycle neglecting breathing motion variability.

METHODS

We integrated a technique to generate a virtual CT (vCT) from 3D MRI with a method for 3D reconstruction from 2D cine-MRI, to produce TR vCTs for dose recalculations. For eight patients, the method allowed evaluating inter-fraction variations at end-exhale and intra-fraction cycle-to-cycle variability within the gating window in terms of tumor displacement and dose to the target and organs at risk.

RESULTS

The median inter-fraction tumor motion was in the range 3.33-12.16 mm, but the target coverage was robust (-0.4% median D95% variation). Concerning cycle-to-cycle variations, the gating technique was effective in limiting tumor displacement (1.35 mm median gating motion) and corresponding dose variations (-3.9% median D95% variation). The larger exposure of organs at risk (duodenum and stomach) was caused by inter-fraction motion, whereas intra-fraction cycle-to-cycle dose variations were limited.

CONCLUSIONS

This study proposed a method for the generation of TR vCTs from MRI, which enabled an off-line evaluation of gating treatment robustness and suggested its feasibility to support treatment planning of pancreatic tumors in CIRT.

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

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

Biological Rationale and Clinical Evidence of Carbon Ion Radiation Therapy for Adenoid Cystic Carcinoma: A Narrative Review

Adenoid cystic carcinoma (ACC) is a rare, basaloid, epithelial tumor, arising mostly from salivary glands. Radiation therapy can be employed as a single modality for unresectable tumors, in an adjuvant setting after uncomplete resection, in case of high-risk pathological features, or for recurrent tumors. Due to ACC intrinsic radioresistance, high linear energy transfer (LET) radiotherapy techniques have been evaluated for ACC irradiation: while fast neutron therapy has now been abandoned due to toxicity concerns, charged particle beams such as protons and carbon ions are at present the beams used for hadron therapy. Carbon ion radiation therapy (CIRT) is currently increasingly used for ACC irradiation. The aim of this review is to describe the immunological, molecular and clinicopathological bases that support ACC treatment with CIRT, as well as to expose the current clinical evidence that reveal the advantages of using CIRT for treating ACC.

Preliminary tests of dosimetric quality and projected therapeutic outcomes of multi-phase 4D radiotherapy with proton and carbon ion beams

Objective. The purpose of this study was to perform preliminary pre-clinical tests to compare the dosimetric quality of two approaches to treating moving tumors with ion beams: synchronously delivering the beam with the motion of a moving planning target volume (PTV) using the recently developed multi-phase 4D dose delivery (MP4D) approach, and asynchronously delivering the ion beam to a motion-encompassing internal tumor volume (ITV) combined with rescanning.Approach. We created 4D optimized treatment plans with proton and carbon ion beams for two patients who had previously received treatment for non-small cell lung cancer. For each patient, we created several treatment plans, using approaches with and without motion mitigation: MP4D, ITV with rescanning, static deliveries to a stationary PTV, and deliveries to a moving tumor without motion compensation. Two sets of plans were optimized with margins or robust uncertainty scenarios. Each treatment plan was delivered using a recently-developed motion-synchronized dose delivery system (M-DDS); dose distributions in water were compared to measurements using gamma index analysis to confirm the accuracy of the calculations. Reconstructed dose distributions on the patient CT were analyzed to assess the dosimetric quality of the deliveries (conformity, uniformity, tumor coverage, and extent of hotspots).Main results. Gamma index analysis pass rates confirmed the accuracy of dose calculations. Dose coverage was >95% for all static and MP4D treatments. The best conformity and the lowest lung doses were achieved with MP4D deliveries. Robust optimization led to higher lung doses compared to conventional optimization for ITV deliveries, but not for MP4D deliveries.Significance. We compared dosimetric quality for two approaches to treating moving tumors with ion beams. Our findings suggest that the MP4D approach, using an M-DDS, provides conformal motion mitigation, with full target coverage and lower OAR doses.

Extension of RBE-weighted 4D particle dose calculation for non-periodic motion

PURPOSE

Highly conformal scanned Carbon Ion Radiotherapy (CIRT) might permit dose escalation and improved local control in advanced stage thoracic tumors, but is challenged by target motion. Dose calculation algorithms typically assume a periodically repeating, regular motion. To assess the effect of realistic, irregular motion, new algorithms of validated accuracy are needed.

METHODS

We extended an in-house treatment planning system to calculate RBE-weighted dose distributions in CIRT on non-periodic CT image sequences. Dosimetric accuracy was validated experimentally on a moving, time-resolved ionization chamber array. Log-file based dose reconstructions were compared by gamma analysis and correlation to measurements at every intermediate detector frame during delivery. The impact of irregular motion on treatment quality was simulated on a virtual 4DCT thorax phantom. Periodic motion was compared to motion with varying amplitude and period ± baseline drift. Rescanning as a mitigation strategy was assessed on all scenarios.

RESULTS

In experimental validation, average gamma pass rates were 99.89+-0.30% for 3%/3 mm and 88.2+-2.2% for 2%/2 mm criteria. Average correlation for integral dose distributions was 0.990±0.002. Median correlation for single 200 ms frames was 0.947±0.006. In the simulations, irregular motion deteriorated V95 target coverage to 81.2%, 76.6% and 79.0% for regular, irregular motion and irregular motion with base-line drift, respectively. Rescanning restored V95 to >98% for both scenarios without baseline drift, but not with additional baseline drift at 83.7%.

CONCLUSIONS

The validated algorithm permits to study the effects of irregular motion and to develop and adapt appropriate motion mitigation techniques.

Dosimetric Validation of a System to Treat Moving Tumors Using Scanned Ion Beams That Are Synchronized With Anatomical Motion

Purpose

The purpose of this study was to validate the dosimetric performance of scanned ion beam deliveries with motion-synchronization to heterogenous targets.

Methods

A 4D library of treatment plans, comprised of up to 10 3D sub-plans, was created with robust and conventional 4D optimization methods. Each sub-plan corresponded to one phase of periodic target motion. The plan libraries were delivered to a test phantom, comprising plastic slabs, dosimeters, and heterogenous phantoms. This phantom emulated range changes that occur when treating moving tumors. Similar treatment plans, but without motion synchronization, were also delivered to a test phantom with a stationary target and to a moving target; these were used to assess how the target motion degrades the quality of dose distributions and the extent to which motion synchronization can improve dosimetric quality. The accuracy of calculated dose distributions was verified by comparison with corresponding measurements. Comparisons utilized the gamma index analysis method. Plan quality was assessed based on conformity, dose coverage, overdose, and homogeneity values, each extracted from calculated dose distributions.

Results

High pass rates for the gamma index analysis confirmed that the methods used to calculate and reconstruct dose distributions were sufficiently accurate for the purposes of this study. Calculated and reconstructed dose distributions revealed that the motion-synchronized and static deliveries exhibited similar quality in terms of dose coverage, overdose, and homogeneity for all deliveries considered. Motion-synchronization substantially improved conformity in deliveries with moving targets. Importantly, measurements at multiple locations within the target also confirmed that the motion-synchronized delivery system satisfactorily compensated for changes in beam range caused by the phantom motion. Specifically, the overall planning and delivery approach achieved the desired dose distribution by avoiding range undershoots and overshoots caused by tumor motion.

Conclusions

We validated a dose delivery system that synchronizes the movement of the ion beam to that of a moving target in a test phantom. Measured and calculated dose distributions revealed that this system satisfactorily compensated for target motion in the presence of beam range changes due to target motion. The implication of this finding is that the prototype system is suitable for additional preclinical research studies, such as irregular anatomic motion.

A Modular System for Treating Moving Anatomical Targets With Scanned Ion Beams at Multiple Facilities: Pre-Clinical Testing for Quality and Safety of Beam Delivery

Background

Quality management and safety are integral to modern radiotherapy. New radiotherapy technologies require new consensus guidelines on quality and safety. Established analysis strategies, such as the failure modes and effects analysis (FMEA) and incident learning systems have been developed as tools to assess the safety of several types of radiation therapies. An extensive literature documents the widespread application of risk analysis methods to photon radiation therapy. Relatively little attention has been paid to performing risk analyses of nascent radiation therapy systems to treat moving tumors with scanned heavy ion beams. The purpose of this study was to apply a comprehensive safety analysis strategy to a motion-synchronized dose delivery system (M-DDS) for ion therapy.

Methods

We applied a risk analysis method to new treatment planning and treatment delivery processes with scanned heavy ion beams. The processes utilize a prototype, modular dose delivery system, currently undergoing preclinical testing, that provides new capabilities for treating moving anatomy. Each step in the treatment process was listed in a process map, potential errors for each step were identified and scored using the risk probability number in an FMEA, and the possible causes of each error were described in a fault tree analysis. Solutions were identified to mitigate the risk of these errors, including permanent corrective actions, periodic quality assurance (QA) tests, and patient specific QA (PSQA) tests. Each solution was tested experimentally.

Results

The analysis revealed 58 potential errors that could compromise beam delivery quality or safety. Each of the 14 binary (pass-or-fail) tests passed. Each of the nine QA and four PSQA tests were within anticipated clinical specifications. The modular M-DDS was modified accordingly, and was found to function at two centers.

Conclusion

We have applied a comprehensive risk analysis strategy to the M-DDS and shown that it is a clinically viable motion mitigation strategy. The described strategy can be utilized at any ion therapy center that operates with the modular M-DDS. The approach can also be adapted for use at other facilities and can be combined with existing safety analysis systems.

Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors

PURPOSE

To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario.

MATERIALS AND METHODS

The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed.

RESULTS

For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D_95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D_95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D_95% and gastro-intestinal OARs D_2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs.

CONCLUSIONS

The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.

The potential of Gantry beamline large momentum acceptance for real time tumour tracking in pencil beam scanning proton therapy

Tumour tracking is an advanced radiotherapy technique for precise treatment of tumours subject to organ motion. In this work, we addressed crucial aspects of dose delivery for its realisation in pencil beam scanning proton therapy, exploring the momentum acceptance and global achromaticity of a Gantry beamline to perform continuous energy regulation with a standard upstream degrader. This novel approach is validated on simulation data from three geometric phantoms of increasing complexity and one liver cancer patient using 4D dose calculations. Results from a standard high-to-low beamline ramping scheme were compared to alternative energy meandering schemes including combinations with rescanning. Target coverage and dose conformity were generally well recovered with tumour tracking even though for particularly small targets, large variations are reported for the different approaches. Meandering in energy while rescanning has a positive impact on target homogeneity and similarly, hot spots outside the targets are mitigated with a relatively fast convergence rate for most tracking scenarios, halving the volume of hot spots after as little as 3 rescans. This work investigates the yet unexplored potential of having a large momentum acceptance in medical beam line, and provides an alternative to take tumour tracking with particle therapy closer to clinical translation.