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

A review of the clinical introduction of 4D particle therapy research concepts

Background and purpose

Many 4D particle therapy research concepts have been recently translated into clinics, however, remaining substantial differences depend on the indication and institute-related aspects. This work aims to summarise current state-of-the-art 4D particle therapy technology and outline a roadmap for future research and developments.

Material and methods

This review focused on the clinical implementation of 4D approaches for imaging, treatment planning, delivery and evaluation based on the 2021 and 2022 4D Treatment Workshops for Particle Therapy as well as a review of the most recent surveys, guidelines and scientific papers dedicated to this topic.

Results

Available technological capabilities for motion surveillance and compensation determined the course of each 4D particle treatment. 4D motion management, delivery techniques and strategies including imaging were diverse and depended on many factors. These included aspects of motion amplitude, tumour location, as well as accelerator technology driving the necessity of centre-specific dosimetric validation. Novel methodologies for X-ray based image processing and MRI for real-time tumour tracking and motion management were shown to have a large potential for online and offline adaptation schemes compensating for potential anatomical changes over the treatment course. The latest research developments were dominated by particle imaging, artificial intelligence methods and FLASH adding another level of complexity but also opportunities in the context of 4D treatments.

Conclusion

This review showed that the rapid technological advances in radiation oncology together with the available intrafractional motion management and adaptive strategies paved the way towards clinical implementation.

Automatic lung segmentation of magnetic resonance images: A new approach applied to healthy volunteers undergoing enhanced Deep-Inspiration-Breath-Hold for motion-mitigated 4D proton therapy of lung tumors

Background and purpose

Respiratory suppression techniques represent an effective motion mitigation strategy for 4D-irradiation of lung tumors with protons. A magnetic resonance imaging (MRI)-based study applied and analyzed methods for this purpose, including enhanced Deep-Inspiration-Breath-Hold (eDIBH). Twenty-one healthy volunteers (41-58 years) underwent thoracic MR scans in four imaging sessions containing two eDIBH-guided MRIs per session to simulate motion-dependent irradiation conditions. The automated MRI segmentation algorithm presented here was critical in determining the lung volumes (LVs) achieved during eDIBH.

Materials and methods

The study included 168 MRIs acquired under eDIBH conditions. The lung segmentation algorithm consisted of four analysis steps: (i) image preprocessing, (ii) MRI histogram analysis with thresholding, (iii) automatic segmentation, (iv) 3D-clustering. To validate the algorithm, 46 eDIBH-MRIs were manually contoured. Sørensen-Dice similarity coefficients (DSCs) and relative deviations of LVs were determined as similarity measures. Assessment of intrasessional and intersessional LV variations and their differences provided estimates of statistical and systematic errors.

Results

Lung segmentation time for 100 2D-MRI planes was ∼ 10 s. Compared to manual lung contouring, the median DSC was 0.94 with a lower 95 % confidence level (CL) of 0.92. The relative volume deviations yielded a median value of 0.059 and 95 % CLs of -0.013 and 0.13. Artifact-based volume errors, mainly of the trachea, were estimated. Estimated statistical and systematic errors ranged between 6 and 8 %.

Conclusions

The presented analytical algorithm is fast, precise, and readily available. The results are comparable to time-consuming, manual segmentations and other automatic segmentation approaches. Post-processing to remove image artifacts is under development.

A bi-directional beam-line energy ramping for efficient patient treatment with scanned proton therapy

Objective.The treatment of mobile tumours using Pencil Beam Scanning (PBS) has become more prevalent in the last decade. However, to achieve the same beam delivery quality as for static tumours, treatments have to be combined with motion mitigation techniques, not limited but including, breath hold, gating and re-scanning, which typically prolong treatment time. In this article we present a novel method of bi-directional energy modulation and demonstrate our initial experience in improvement of treatment efficiency. Approach.At Paul Scherrer Institute Gantry 2 mobile tumours are treated by combining PBS with gating and volumetric re-scanning (VR), where the target volume is irradiated multiple times. Initial implementation of VR used only descending beam energies, creating a substantial dead time due to the beam-line initialization (ramping) before each re-scan. In 2019 we commissioned an energy meandering strategy that allows us to avoid beam line ramping in-between energy series while maintaining beam delivery quality.Main results.The measured beam parameters difference for both energy sequence are in the order of the typical daily variations: 0.2 mm in beam position and 0.2 mm in range. Using machine log files, we performed point-to-point dose difference calculations between original and new applications where we observed dose differences of less than 2%. After three years of operation employing bi-directional energy modulation, we have analysed the individual beam delivery time for 181 patients and have compared this to simulations of the timing behaviour assuming uni-directional energy sequence application. Depending on treatment complexity, we obtained plan delivery time reductions of up to 55%, with a median time gain of 17% for all types of treatments.Significance. Bi-directional energy modulation can help improving patient treatment efficiency by reducing delivery times especially for complex and specialised irradiations. It could be implemented in many existing facilities without significant additional hardware upgrades.

High-speed low-noise optical respiratory monitoring for spot scanning proton therapy

PURPOSE

To investigate a novel optical markerless respiratory sensor for surface guided spot scanning proton therapy and to measure its main technical characteristics.

METHODS

The main characteristics of the respiratory sensor including sensitivity, linearity, noise, signal-to-noise, and time delay were measured using a dynamic phantom and electrical measuring equipment on a laboratory stand. The respiratory signals of free breathing and deep-inspiration breath-hold patterns were acquired for various distances with a volunteer. A comparative analysis of this sensor with existing commercially available and experimental respiratory monitoring systems was carried out based on several criteria including principle of operation, patient contact, application to proton therapy, distance range, accuracy (noise, signal-to-noise ratio), and time delay (sampling rate).

RESULTS

The sensor provides optical respiratory monitoring of the chest surface over a distance range of 0.4-1.2 m with the RMS noise of 0.03-0.60 mm, SNR of 40-15 dB (for motion with peak-to-peak of 10 mm), and time delay of 1.2 ± 0.2 ms.

CONCLUSIONS

The investigated optical respiratory sensor was found to be appropriate to use in surface guided spot scanning proton therapy. This sensor combined with a fast respiratory signal processing algorithm may provide accurate beam control and a fast response in patients' irregular breathing movements. A careful study of correlation between the respiratory signal and 4DCT data of tumor position will be required before clinical implementation.

Parameter based 4D dose calculations for proton therapy

Background and purpose Retrospective log file-based analysis provides the actual dose delivered based on the patient's breathing and the daily beam-delivery dynamics. To predict the motion sensitivity of the treatment plan on a patient-specific basis before treatment start a prospective tool is required. Such a parameter-based tool has been investigated with the aim to be used in clinical routine. Materials and Methods 4D dose calculations (4DDC) were performed for seven cancer patients with small breathing motion treated with scanned pulsed proton beams. Validation of the parameter-based 4DDC (p-4DDC) method was performed with an anthropomorphic phantom and patient data employing measurements and a log file-based 4DDC tool. The dose volume histogram parameters (Dx%) were investigated for the target and the organs at risk, compared to static and the file-based approach. Results The difference between the measured and the p-4DDC dose was within the deviation of the measurements. The maximum deviation was 0.4Gy. For the planning target volume D98% varied up to 15% compared to the static scenario, while the results from the log file and p-4DDC agreed within 2%. For the liver patients, D33%liver deviated up to 35% compared to static and 10% comparing the two 4DDC tools, while for the pancreas patients the D1%stomach varied up to 45% and 11%, respectively. Conclusion The results showed that p-4DDC could be used prospectively. The next step will be the clinical implementation of the p-4DDC tool, which can support a decision to either adapt the treatment plan or apply motion mitigation strategies.

A survey of practice patterns for real-time intrafractional motion-management in particle therapy

Background and purpose

Organ motion compromises accurate particle therapy delivery. This study reports on the practice patterns for real-time intrafractional motion-management in particle therapy to evaluate current clinical practice and wishes and barriers to implementation.

Materials and methods

An institutional questionnaire was distributed to particle therapy centres worldwide (7/2020-6/2021) asking which type(s) of real-time respiratory motion management (RRMM) methods were used, for which treatment sites, and what were the wishes and barriers to implementation. This was followed by a three-round DELPHI consensus analysis (10/2022) to define recommendations on required actions and future vision. With 70 responses from 17 countries, response rate was 100% for Europe (23/23 centres), 96% for Japan (22/23) and 53% for USA (20/38).

Results

Of the 68 clinically operational centres, 85% used RRMM, with 41% using both rescanning and active methods. Sixty-four percent used active-RRMM for at least one treatment site, mostly with gating guided by an external marker. Forty-eight percent of active-RRMM users wished to expand or change their RRMM technique. The main barriers were technical limitations and limited resources. From the DELPHI analysis, optimisation of rescanning parameters, improvement of motion models, and pre-treatment 4D evaluation were unanimously considered clinically important future focus. 4D dose calculation was identified as the top requirement for future commercial treatment planning software.

Conclusion

 A majority of particle therapy centres have implemented RRMM. Still, further development and clinical integration were desired by most centres. Joint industry, clinical and research efforts are needed to translate innovation into efficient workflows for broad-scale implementation.

The Influence of Motion on the Delivery Accuracy When Comparing Actively Scanned Carbon Ions versus Protons at a Synchrotron-Based Radiotherapy Facility

Simple Summary

The interplay of breathing and beam motion reduces the efficacy of particle irradiation in moving tumours. The effect of motion on protons and carbon ion treatments was investigated dosimetrically and the results were benchmarked against each other by employing an anthropomorphic thorax phantom that was able to simulate tumour, rib, and lung motion. The critical question was whether target coverage and organ-at-risk sparing could be maintained when the application of simple motion mitigation was addressed. Special focus was put on unique synchrotron characteristics, such as pulsed beam delivery and beam intensity variations. It could be demonstrated that the effect of motion was greater for carbon ions than for protons. These findings demonstrated the need for applying motion mitigation techniques depending on the motion amplitude, particle type, and treatment prescription considering complex time correlations.

Abstract

Motion amplitudes, in need of mitigation for moving targets irradiated with pulsed carbon ions and protons, were identified to guide the decision on treatment and motion mitigation strategy. Measurements with PinPoint ionisation chambers positioned in an anthropomorphic breathing phantom were acquired to investigate different tumour motion scenarios, including rib and lung movements. The effect of beam delivery dynamics and spot characteristics was considered. The dose in the tumour centre was deteriorated up to 10% for carbon ions but only up to 5% for protons. Dose deviations in the penumbra increased by a factor of two when comparing carbon ions to protons, ranging from 2 to 30% for an increasing motion amplitude that was strongly dependent on the beam intensity. Layer rescanning was able to diminish the dose distortion caused by tumour motion, but an increase in spot size could reduce it even further to 5% within the target and 10% at the penumbra. An increased need for motion mitigation of carbon ions compared to protons was identified to assure target coverage and sparing of adjacent organs at risk in the penumbra region and outside the target. For the clinical implementation of moving target treatments at a synchrotron-based particle facility complex, time dependencies needed to be considered.

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.

Enhanced Deep-Inspiration Breath Hold Superior to High-Frequency Percussive Ventilation for Respiratory Motion Mitigation: A Physiology-Driven, MRI-Guided Assessment Toward Optimized Lung Cancer Treatment With Proton Therapy

Background: To safely treat lung tumors using particle radiation therapy (PRT), motion-mitigation strategies are of critical importance to ensure precise irradiation. Therefore, we compared applicability, effectiveness, reproducibility, and subjects' acceptance of enhanced deep-inspiration breath hold (eDIBH) with high-frequency percussive ventilation (HFPV) by MRI assessment within 1 month.

Methods: Twenty-one healthy subjects (12 males/9 females; age: 49.5 ± 5.8 years; BMI: 24.7 ± 3.3 kg/m⁻²) performed two 1.5 T MRI scans in four visits at weekly intervals under eDIBH and HFPV conditions, accompanied by daily, home-based breath-hold training and spirometric assessments over a 3-week period. eDIBH consisted of 8-min 100% O₂ breathing (3 min resting ventilation, 5 min controlled hyperventilation) prior to breath hold. HFPV was set at 200–250 pulses min⁻¹ and 0.8–1.2 bar. Subjects' acceptance and preference were evaluated by questionnaire. To quantify inter- and intrafractional changes, a lung distance metric representing lung topography was computed for 10 reference points: a motion-invariant spinal cord and nine lung structure contours (LSCs: apex, carina, diaphragm, and six vessels as tumor surrogates distributed equally across the lung). To parameterize individual LSC localizability, measures of their spatial variabilities were introduced and lung volumes calculated by automated MRI analysis.

Results: eDIBH increased breath-hold duration by > 100% up to 173 ± 73 s at visit 1, and to 217 ± 67 s after 3 weeks of home-based training at visit 4 (p < 0.001). Measures of vital capacity and lung volume remained constant over the 3-week period. Two vessels in the lower lung segment and the diaphragm yielded a two- to threefold improved positional stability with eDIBH, whereby absolute distance variability was significantly smaller for five LSCs; ≥70% of subjects showed significantly better intrafractional lung motion mitigation under reproducible conditions with eDIBH compared with HFPV with smaller ranges most apparent in the anterior-posterior and cranial-caudal directions. Approximately 80% of subjects preferred eDIBH over HFPV, with “less discomfort” named as most frequent reason.

Conclusions: Both, eDIBH, and HFPV were well-tolerated. eDIBH duration was long enough to allow for potential PRT. Variability in lung volume was smaller and position of lung structures more precise with eDIBH. Subjects preferred eDIBH over HFPV. Thus, eDIBH is a very promising tool for lung tumor therapy with PRT, and further investigation of its applicability in patients is warranted.

Technical Note: Investigating interplay effects in pencil beam scanning proton therapy with a 4D XCAT phantom within the RayStation treatment planning system

PURPOSE

Pencil beam scanning (PBS) for moving targets is known to be impacted by interplay effects. Four-dimensional computed tomography (4DCT)-based motion evaluation is crucial for understanding interplay and developing mitigation strategies. Availability of high-quality 4DCTs with variable breathing traces is limited. Purpose of this work is the development of a framework for interplay analysis using 4D-XCAT phantoms in conjunction with time-resolved irradiation patterns in a commercial treatment planning system (TPS). Four-dimensional dynamically accumulated dose distributions (4DDDs) are simulated in an in-silico study for a PBS liver treatment.

METHODS

An XCAT phantom with 50 phases, varying linearly in amplitude each by 1 mm, was combined with the RayStation TPS (7.99.10). Deformable registration was used with time-resolved dose calculation, mapping XCAT phases to motion signals. To illustrate the applicability of the method a two-field liver irradiation plan was used. A variety sin⁴ type motion signals, varying in amplitude (1-20 mm), period (1.6-5.2 s) and phase (0-2π) were applied. Either single variable variations or random combinations were selected. The interplay effect within a clinical target (5 cm diameter) was characterized in terms of homogeneity index (HI5), with and without five paintings. In total 2092 scenarios were analyzed within RayStation.

RESULTS

A framework is presented for interplay research, allowing for flexibility in determining motion management techniques, increasing reproducibility, and enabling comparisons of different methods. A case study showed the interplay effect was correlated with amplitude and strongly affected by the starting phase, leading to large variance. The average of all scenarios (single fraction) resulted in HI5 of 0.31 (±0.11), while introduction of five times layered repainting reduced this to 0.11(±0.03).

CONCLUSION

The developed framework, which uses the XCAT phantom and RayStation, allows detailed analysis of motion in context of PBS with comparable results to clinical cases. Flexibility in defining motion patterns for detailed anatomies in combination with time-resolved dose calculation, facilitates investigation of optimal treatment and motion mitigation strategies.

Clinical practice vs. state-of-the-art research and future visions: Report on the 4D treatment planning workshop for particle therapy - Edition 2018 and 2019

The 4D Treatment Planning Workshop for Particle Therapy, a workshop dedicated to the treatment of moving targets with scanned particle beams, started in 2009 and since then has been organized annually. The mission of the workshop is to create an informal ground for clinical medical physicists, medical physics researchers and medical doctors interested in the development of the 4D technology, protocols and their translation into clinical practice. The 10th and 11th editions of the workshop took place in Sapporo, Japan in 2018 and Krakow, Poland in 2019, respectively. This review report from the Sapporo and Krakow workshops is structured in two parts, according to the workshop programs. The first part comprises clinicians and physicists review of the status of 4D clinical implementations. Corresponding talks were given by speakers from five centers around the world: Maastro Clinic (The Netherlands), University Medical Center Groningen (The Netherlands), MD Anderson Cancer Center (United States), University of Pennsylvania (United States) and The Proton Beam Therapy Center of Hokkaido University Hospital (Japan). The second part is dedicated to novelties in 4D research, i.e. motion modelling, artificial intelligence and new technologies which are currently being investigated in the radiotherapy field.

Mitigation of motion effects in pencil-beam scanning - Impact of repainting on 4D robustly optimized proton treatment plans for hepatocellular carcinoma

Proton fields delivered by the active scanning technique can be interfered with the intrafractional motion. This in-silico study seeks to mitigate the dosimetric impacts of motion artifacts, especially its interplay with the time-modulated dose delivery. Here four-dimensional (4d) robust optimization and dose repainting, which is the multiple application of the same field with reduced fluence, were combined. Two types of repainting were considered: layered and volumetric repainting. The time-resolved dose calculation, which is necessary to quantify the interplay effect, was integrated into the treatment planning system and validated. Nine clinical cases of hepatocellular carcinoma (HCC) showing motion in the range of 0.4-1.5cm were studied. It was found that the repainted delivery of 4D robustly optimized plans reduced the impact of interplay effect as quantified by the homogeneity index within the clinical target volume (CTV) to a tolerable level. Similarly, the fractional over- and underdosage was reduced sufficiently for some HCC cases to achieve the purpose of motion management. This holds true for both investigated types of repainting with small dosimetric advantages of volume repainting over layered repainting. Volume repainting, however, cannot be applied clinically in proton centers with slow energy changes. Thus, it served as a reference in the in-silico evaluation. It is recommended to perform the dynamic dose calculation for individual cases to judge if robust optimization in conjunction with repainting is sufficient to keep the interplay effect within bounds.

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

PURPOSE

The 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.

METHODS

We 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.

RESULTS

Acceptable dosimetric performance and safety characteristics were observed for all treatment plans.

CONCLUSIONS

We 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.

Time-resolved dosimetry for validation of 4D dose calculation in PBS proton therapy

Four-dimensional dose calculation (4D-DC) is crucial for predicting the dosimetric outcome in the presence of intra-fractional organ motion. Time-resolved dosimetry can provide significant insights into 4D pencil beam scanning dose accumulation and is therefore irreplaceable for benchmarking 4D-DC. In this study a novel approach of time-resolved dosimetry using five PinPoint ionization chambers (ICs) embedded in an anthropomorphic dynamic phantom was employed and validated against beam delivery details. Beam intensity variations as well as the beam delivery time structure were well reflected with an accuracy comparable to the temporal resolution of the IC measurements. The 4D dosimetry approach was further applied for benchmarking the 4D-DC implemented in the RayStation 6.99 treatment planning system. Agreement between computed values and measurements was investigated for (i) partial doses based on individual breathing phases, and (ii) temporally distributed cumulative doses. For varied beam delivery and patient-related parameters the average unsigned dose difference for (i) was 0.04 ± 0.03 Gy over all considered IC measurement values, while the prescribed physical dose was 2 Gy. By implementing (ii), a strong effect of the dose gradient on measurement accuracy was observed. The gradient originated from scanned beam energy modulation and target motion transversal to the beam. Excluding measurements in the high gradient the relative dose difference between measurements and 4D-DCs for a given treatment plan at the end of delivery was 3.5% on average and 6.6% at maximum over measurement points inside the target. Overall, the agreement between 4D dose measurements in the moving phantom and retrospective 4D-DC was found to be comparable to the static dose differences for all delivery scenarios. The presented 4D-DC has been proven to be suitable for simulating treatment deliveries with various beam- as well as patient-specific parameters and can therefore be employed for dosimetric validation of different motion mitigation techniques.

Dose assessment for patients with stage I non-small cell lung cancer receiving passive scattering carbon-ion radiotherapy using daily computed tomographic images: A prospective study

BACKGROUND AND PURPOSE

This study aimed to assess dose distributions for stage I non-small cell lung cancer (NSCLC) with passive scattering carbon-ion radiotherapy (C-ion RT) using daily computed tomography (CT) images.

MATERIALS AND METHODS

We enrolled 10 patients with stage I NSCLC and acquired a total of 40 pre-fractional CT image series under the same settings as the planning CT images. These CT images were registered with planning CT images for dose evaluation using both bone matching (BM) and tumor matching (TM). Using deformable image registration, we generated accumulated doses. Moreover, the volumetric dose parameters were compared in terms of tumor coverage and lung exposure and statistical analyses were performed.

RESULTS

Overall, 25% of 40 fractional dose distributions were unacceptable with BM, compared with 2.5% with TM (P < 0.001). Using BM, three patients' accumulated dose distributions were unacceptable; however, all were satisfactory with TM (P < 0.001). No differences were observed in water-equivalent path length (WEL). The required margins in patients with poor dose distribution were 5.9 and 4.4 mm for BM and TM, respectively.

CONCLUSIONS

This study establishes that CT image-based TM is robust compared with conventional BM for both daily and accumulated dose distributions. The effects of changes in WEL seem to be limited. Hence, daily CT alignment is recommended for patients with stage I NSCLC receiving C-ion RT.

Virtual 4DCT from 4DMRI for the management of respiratory motion in carbon ion therapy of abdominal tumors

PURPOSE

To evaluate a method for generating virtual four-dimensional computed tomography (4DCT) from four-dimensional magnetic resonance imaging (4DMRI) data in carbon ion radiotherapy with pencil beam scanning for abdominal tumors.

METHODS

Deformable image registration is used to: (a) register each respiratory phase of the 4DMRI to the end-exhale MRI; (b) register the reference end-exhale CT to the end-exhale MRI volume; (c) generate the virtual 4DCT by warping the registered CT according to the obtained deformation fields. A respiratory-gated carbon ion treatment plan is optimized on the planning 4DCT and the corresponding dose distribution is recalculated on the virtual 4DCT. The method was validated on a digital anthropomorphic phantom and tested on eight patients (18 acquisitions). For the phantom, a ground truth dataset was available to assess the method performances from the geometrical and dosimetric standpoints. For the patients, the virtual 4DCT was compared with the planning 4DCT.

RESULTS

In the phantom, the method exhibits a geometrical accuracy within the voxel size and Dose Volume Histograms deviations up to 3.3% for target V_95% (mean dose difference ≤ 0.2% of the prescription dose, gamma pass rate > 98%). For patients, the virtual and the planning 4DCTs show good agreement at end-exhale (3% median D_95% difference), whereas other respiratory phases exhibit moderate motion variability with consequent dose discrepancies, confirming the need for motion mitigation strategies during treatment.

CONCLUSIONS

The virtual 4DCT approach is feasible to evaluate treatment plan robustness against intra- and interfraction motion in carbon ion therapy delivered at the abdominal site.

Dynamic lung phantom commissioning for 4D dose assessment in proton therapy

Anthropomorphic phantoms mimicking organ and tumor motion of patients are essential for end-to-end testing of motion mitigation techniques in ion beam therapy. In this work a commissioning procedure developed with the in-house designed respiratory phantom ARDOS (Advanced Radiation DOSimetry system) is presented. The phantom was tested and benchmarked for 4D dose verification in proton therapy, which included: characterization of the tissue equivalent materials from computed tomography (CT) imaging, assessment of dose calculation accuracy in critical structures of the phantom, and testing various detectors for proton dosimetry in the ARDOS phantom. To prove the validity of the CT calibration curve, measured relative stopping powers (RSP) of the ARDOS materials were compared with values from CTs: original and overwritten with known material parameters. Override of rib- and soft-tissue phantom components improved RSP accuracy while inhomogeneous lung tissue, represented by the balsa wood, was better modelled by the CT Hounsfield units. Monte Carlo (MC) dose calculations were benchmarked against measurements with a reference Farmer chamber embedded in ARDOS material samples showing less than 3% relative dose difference. Differences between MC calculated dose distributions and those calculated by analytical algorithms for the ARDOS geometry were higher than 20% of the prescribed dose, depending on the position in the phantom. Pinpoint ionization chambers and thermoluminescence dosimeters showed differences of up to 5.5% compared to MC dose calculations for all lung setups in the static phantom. They were also able to detect dose distortions due to motion. EBT3 film dosimetry was shown to be suitable for 2D relative dose characterization, which could provide extended information on dose distributions in the penumbra area. The presented methodology and results can be used for drafting general recommendations for dynamic phantom commissioning, which is an essential step towards end-to-end evaluation of motion mitigation techniques in ion beam therapy.

Treatment planning for proton therapy: what is needed in the next 10 years?

Treatment planning is the process where the prescription of the radiation oncologist is translated into a deliverable treatment. With the complexity of contemporary radiotherapy, treatment planning cannot be performed without a computerized treatment planning system. Proton therapy (PT) enables highly conformal treatment plans with a minimum of dose to tissues outside the target volume, but to obtain the most optimal plan for the treatment, there are a multitude of parameters that need to be addressed. In this review areas of ongoing improvements and research in the field of PT treatment planning are identified and discussed. The main focus is on issues of immediate clinical and practical relevance to the PT community highlighting the needs for the near future but also in a longer perspective. We anticipate that the manual tasks performed by treatment planners in the future will involve a high degree of computational thinking, as many issues can be solved much better by e.g. scripting. More accurate and faster dose calculation algorithms are needed, automation for contouring and planning is required and practical tools to handle the variable biological efficiency in PT is urgently demanded just to mention a few of the expected improvements over the coming 10 years.

Real-time intrafraction motion monitoring in external beam radiotherapy

Radiotherapy (RT) aims to deliver a spatially conformal dose of radiation to tumours while maximizing the dose sparing to healthy tissues. However, the internal patient anatomy is constantly moving due to respiratory, cardiac, gastrointestinal and urinary activity. The long term goal of the RT community to 'see what we treat, as we treat' and to act on this information instantaneously has resulted in rapid technological innovation. Specialized treatment machines, such as robotic or gimbal-steered linear accelerators (linac) with in-room imaging suites, have been developed specifically for real-time treatment adaptation. Additional equipment, such as stereoscopic kilovoltage (kV) imaging, ultrasound transducers and electromagnetic transponders, has been developed for intrafraction motion monitoring on conventional linacs. Magnetic resonance imaging (MRI) has been integrated with cobalt treatment units and more recently with linacs. In addition to hardware innovation, software development has played a substantial role in the development of motion monitoring methods based on respiratory motion surrogates and planar kV or Megavoltage (MV) imaging that is available on standard equipped linacs. In this paper, we review and compare the different intrafraction motion monitoring methods proposed in the literature and demonstrated in real-time on clinical data as well as their possible future developments. We then discuss general considerations on validation and quality assurance for clinical implementation. Besides photon RT, particle therapy is increasingly used to treat moving targets. However, transferring motion monitoring technologies from linacs to particle beam lines presents substantial challenges. Lessons learned from the implementation of real-time intrafraction monitoring for photon RT will be used as a basis to discuss the implementation of these methods for particle RT.

Dosimetric uncertainties as a result of temporal resolution in 4D dose calculations for PBS proton therapy

This work investigates the dosimetric impact on 4D dose distribution estimation for pencil beam scanned (PBS) proton therapy as function of the temporal resolution used for the time resolved dose calculation. For three liver patients (CTV volume: 403/122/264 cc), 10-phase 4DCT-MRI datasets with ~15 mm tumour motion were simulated for seven different motion periods (2-8 s). 4D dose distributions were calculated and compared by considering both coarser and finer temporal resolutions (200-800 ms and 20 ms). Single scanned 4D plans for seven fraction doses (0.7/2/4/6/8/10/12 Gy) were investigated, whose dose delivery timelines were simulated by assuming two types of PBS scanning modes: (1) layer-wise raster scanning with varying dose rate per layer and (2) fixed dose rate, discrete scanning. For both delivery scenarios, dosimetric assessments were performed by comparing corresponding dose distributions derived from the two 4D dose calculation (4DDC) results. Differences were quantified as the difference in D5-D95 of the CTV and by comparing total volume of the CTV receiving point-to-point absolute dose difference more than 5%. Our results show that varying temporal resolution in 4DDC has a direct influence on the final accumulated dose distribution. For all scenarios, patients, fraction doses and motion periods studied, pronounced dose differences can be observed between the two 4DDC results. However, the magnitude of differences varies depending on the selected PBS scanning model and prescribed dose per field. For fixed dose rate delivery, the average duration of the delivery of each spot increases for hypo-fractionated treatments, enhancing the benefit of using a finer temporal resolution for 4DDC. In particular, for fraction doses  >4 Gy and motion periods less than 4 s, warping the dose between discrete 4DCT phases can over predict the interplay effect (D5-D95 in CTV) by 3%-10% compared to the use of a finer temporal resolution, resulting in more than 20% of CTV voxels having absolute dose differences of over 5% between the two 4DDC approaches. These findings emphasize the importance for PBS 4DDC using finer temporal resolutions than provided by conventional 4D dose accumulation techniques. In particular, the observed differences in dosimetric effects using the fine temporal resolution provided by dose warping cannot be neglected for hypo-fractionation and short breathing periods, especially when using constant dose rates for dose delivery.

RIDOS: A new system for online computation of the delivered dose distributions in scanning ion beam therapy

PURPOSE

To describe a new system for scanned ion beam therapy, named RIDOS (Real-time Ion DOse planning and delivery System), which performs real time delivered dose verification integrating the information from a clinical beam monitoring system with a Graphic Processing Unit (GPU) based dose calculation in patient Computed Tomography.

METHODS

A benchmarked dose computation algorithm for scanned ion beams has been parallelized and adapted to run on a GPU architecture. A workstation equipped with a NVIDIA GPU has been interfaced through a National Instruments PXI-crate with the dose delivery system of the Italian National Center of Oncological Hadrontherapy (CNAO) to receive in real-time the measured beam parameters. Data from a patient monitoring system are also collected to associate the respiratory phases with each spot during the delivery of the dose. Using both measured and planned spot properties, RIDOS evaluates during the few seconds of inter-spill time the cumulative delivered and prescribed dose distributions and compares them through a fast γ-index algorithm.

RESULTS

The accuracy of the GPU-based algorithms was assessed against the CPU-based ones and the differences were found below 1‰. The cumulative planned and delivered doses are computed at the end of each spill in about 300 ms, while the dose comparison takes approximatively 400 ms. The whole operation provides the results before the next spill starts.

CONCLUSIONS

RIDOS system is able to provide a fast computation of the delivered dose in the inter-spill time of the CNAO facility and allows to monitor online the dose deposition accuracy all along the treatment.