Motion mitigation for lung cancer patients treated with active scanning proton therapy

Cyclotron and linear accelerator generated scanning proton beams for lung cancer SBRT: Interplay effects and mitigations

BACKGROUND

Pencil beam scanning (PBS) proton therapy for moving targets is known to be impacted by interplay effects between the scanning beam and organ motion. While respiratory motion in the thoracic region is the major cause for organ motion, interplay effects depend on the delivery characteristics of proton accelerators.

PURPOSE

To evaluate the impact of different types of PBS proton accelerators and spot sizes on interplay effects, mitigations, and plan quality for Stereotactic Body Radiation Therapy (SBRT) treatment of non-small cell lung cancer (NSCLC).

METHODS

Twenty NSCLC patients treated with photon SBRT were selected to represent varying tumor volumes and respiratory motion amplitudes (median: 0.6 cm with abdominal compression) for this retrospective study. For each patient, plans were created using: (1) cyclotron-generated proton beams (CPB) with spot sizes of σ = 2.7-7.0 mm; (2) linear accelerator proton beams (LPB) (σ = 2.9-5.5 mm); and (3) linear accelerator proton minibeams (LPMB) (σ = 0.9-3.9 mm). The energy switching time is one second for CPB, and 0.005 s for LPMB and LPB. Plans were robustly optimized on the gross tumor volume (GTV) using each individual phase of four-dimensional computed tomography (4DCT) scans. Initially, single-field optimization (SFO) plans were evaluated; if the plan quality did not meet the dosimetric requirement, multi-field optimization (MFO) was used. MFO plans were created for all patients for comparisons. For each patient, all plans were normalized to have the same dose received by 99% of the GTV. Interplay effects were evaluated by computing the dose on 10 breathing phases, based on the spot distribution. Volumetric repainting (VR) was performed 2-6 times for each plan. We compared volume receiving 100% of the prescribed dose (V100%RX) of the GTV, and normal lung V20Gy.

RESULTS

Twelve of 20 plans can be optimized sufficiently with SFO. SFO plans were less sensitive to the interplay effect compared to MFO plans in terms of target coverage for both LPB and LPMB. The following comparisons showed results utilizing the MFO technique. In the interplay evaluation without repainting, the mean V100%RX of the GTV were 99.42 ± 0.6%, 97.52 ± 3.9%, and 94.49 ± 7.3% for CPB, LPB, and LPMB plans, respectively. Following VR (2 × for CPB; 3 × for LPB; 5 × for LPMB), V100%RX of the GTV were improved (on average) by 0.13%, 1.84%, and 4.63%, respectively, achieving the acceptance criteria of V100%RX > 95%. Because of fast energy switch in linear accelerator proton machines, the delivery time for VR plans was the lowest for LPB plans, while delivery time for LPMB was on average 1 min longer than CPB plans. The advantage of small spot machines was better sparing in normal lung V20Gy, even when VR was applied.

CONCLUSION

In the absence of repainting, proton machines with large spot sizes generated more robust plans against interplay effects. The number of VR increased with decreasing spot sizes to achieve the acceptance criteria. VR improved the plan robustness against interplay effects for modalities with small spot sizes and fast energy changes, preserving the low dose sparing aspect of the LPMB, even when motion is included.

Comparing interplay effects in scanned proton therapy of lung cancer: Free breathing with various layer and volume rescanning versus respiratory gating with different gate widths

PURPOSE

We investigated interplay effects and treatment time (TT) in scanned proton therapy for lung cancer patients. We compared free-breathing (FB) approaches with multiple rescanning strategies and respiratory-gating (RG) methods with various gating widths to identify the superior irradiation technique.

METHODS

Plans were created with 4/1, 2/2, and 1/4 layered/volume rescans of FB (L4V1, L2V2, and L1V4), and 50%, 30%, and 10% gating widths of the total respiratory curves (G50, G30, and G10) of the RG plans with L4V1. We calculated 4-dimensional dynamic doses assuming a constant sinusoidal curve for six irradiation methods. The reconstructed doses per fraction were compared with planned doses in terms of dose differences in 99% clinical-target-volume (CTV) (ΔD99%), near-maximum dose differences (ΔD2%) at organs-at-risk (OARs), and TT.

RESULTS

The mean/minimum CTV ΔD99% values for FB were -1.0%/-4.9%, -0.8%/-4.3%, and -0.1%/-1.0% for L4V1, L2V2, and L1V4, respectively. Those for RG were -0.3%/-1.7%, -0.1%/-1.0%, and 0.0%/-0.5% for G50, G30, and G10, respectively. The CTV ΔD99% of the RGs with less than 50% gate width and the FBs of L1V4 were within the desired tolerance (±3.0%), and the OARs ΔD2% for RG were lower than those for FB. The mean TTs were 90, 326, 824, 158, 203, and 422 s for L4V1, L2V2, L1V4, G50, G30, and G10, respectively.

CONCLUSIONS

FB (L4V1) is the most efficient treatment, but not necessarily the optimal choice due to interplay effects. To satisfy both TT extensions and interplay, RG with a gate width as large as possible within safety limits is desirable.

Experimental validation of a 4D dynamic dose calculation model for proton pencil beam scanning without spot time stamp considering free-breathing motion

PURPOSE

We developed a 4-dimensional dynamic dose (4DDD) calculation model for proton pencil beam scanning (PBS). This model incorporates the spill start time for all energies and uses the remaining irradiated spot time model instead of irradiated spot time logs. This study aimed to validate the calculation accuracy of a log file-based 4DDD model by comparing it with dose measurements performed under free-breathing conditions, thereby serving as an alternative approach to the conventional log file-based system.

METHODS

Three cubic verification plans were created using a heterogeneous block phantom; these plans were created using 10 phase 4D-CT datasets of the phantom. The CIRS dynamic platform was used to simulate motion with amplitudes of 2.5, 3.75, and 5.0 mm. These plans consisted of eight- and two-layered rescanning techniques. The lateral profiles were measured using a 2D ionization chamber array (2D-array) and EBT3 Gafchromic films at four starting phases, including three sinusoidal curves (periods of 3, 4, and 6 s) and a representative patient curve during actual treatment. 4DDDs were calculated using in-house scripting that assigned a time stamp to each spot and performed dose accumulation using deformable image registration. Furthermore, to evaluate the impact of parameter selection on our 4DDD model calculations, simulations were performed assuming a ±10% change in irradiation time stamp (0.8 ± 0.08 s) and spot scan speed. We evaluated the 2D gamma index and the absolute point doses between the calculated values and the measurements.

RESULTS

The 2D-array measurements revealed that the gamma scores for the static plans (no motion) and 4DDD plans exceeded 97.5% and 93.9% at 3%/3 mm, respectively. The average gamma score of the 4DDD plans was at least 96.1%. When using EBT3 films, the gamma scores of the 4DDD model exceeded 92.4% and 98.7% at 2%/2 mm and 3%/3 mm, respectively. Regarding the 4DDD point dose differences, more than 95% of the dose regions exhibited discrepancies within ±5.0% for 97.7% of the total points across all plans. The spot time assignment accuracy of our 4DDD model was acceptable even with ±10% sensitivity. However, the accuracy of the scan speed, when varied within ±10% sensitivity, was not acceptable (minimum gamma scores of 82.6% and maximum dose difference of 12.9%).

CONCLUSIONS

Our 4DDD calculations under free-breathing conditions using amplitudes of less than 5.0 mm were in good agreement with the measurements regardless of the starting phases, breathing curve patterns (between 3 and 6 s periods), and varying numbers of layered rescanning. The proposed system allows us to evaluate actual irradiated doses in various breathing periods, amplitudes, and starting phases, even on PBS machines without the ability to record spot logs.

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.

Emerging technologies for cancer therapy using accelerated particles

Cancer therapy with accelerated charged particles is one of the most valuable biomedical applications of nuclear physics. The technology has vastly evolved in the past 50 years, the number of clinical centers is exponentially growing, and recent clinical results support the physics and radiobiology rationale that particles should be less toxic and more effective than conventional X-rays for many cancer patients. Charged particles are also the most mature technology for clinical translation of ultra-high dose rate (FLASH) radiotherapy. However, the fraction of patients treated with accelerated particles is still very small and the therapy is only applied to a few solid cancer indications. The growth of particle therapy strongly depends on technological innovations aiming to make the therapy cheaper, more conformal and faster. The most promising solutions to reach these goals are superconductive magnets to build compact accelerators; gantryless beam delivery; online image-guidance and adaptive therapy with the support of machine learning algorithms; and high-intensity accelerators coupled to online imaging. Large international collaborations are needed to hasten the clinical translation of the research results.

Limitations of phase-sorting based pencil beam scanned 4D proton dose calculations under irregular motion

Objective.4D dose calculation (4DDC) for pencil beam scanned (PBS) proton therapy is typically based on phase-sorting of individual pencil beams onto phases of a single breathing cycle 4DCT. Understanding the dosimetric limitations and uncertainties of this approach is essential, especially for the realistic treatment scenario with irregular free breathing motion.Approach.For three liver and three lung cancer patient CTs, the deformable multi-cycle motion from 4DMRIs was used to generate six synthetic 4DCT(MRI)s, providing irregular motion (11/15 cycles for liver/lung; tumor amplitudes ∼4-18 mm). 4DDCs for two-field plans were performed, with the temporal resolution of the pencil beam delivery (4-200 ms) or with 8 phases per breathing cycle (500-1000 ms). For the phase-sorting approach, the tumor center motion was used to determine the phase assignment of each spot. The dose was calculated either using the full free breathing motion or individually repeating each single cycle. Additionally, the use of an irregular surrogate signal prior to 4DDC on a repeated cycle was simulated. The CTV volume with absolute dose differences >5% (V_dosediff>5%) and differences in CTVV_95%andD_5%-D_95%compared to the free breathing scenario were evaluated.Main results.Compared to 4DDC considering the full free breathing motion with finer spot-wise temporal resolution, 4DDC based on a repeated single 4DCT resulted inV_dosediff>5%of on average 34%, which resulted in an overestimation ofV_95%up to 24%. However, surrogate based phase-sorting prior to 4DDC on a single cycle 4DCT, reduced the averageV_dosediff>5%to 16% (overestimationV_95%up to 19%). The 4DDC results were greatly influenced by the choice of reference cycle (V_dosediff>5%up to 55%) and differences due to temporal resolution were much smaller (V_dosediff>5%up to 10%).Significance.It is important to properly consider motion irregularity in 4D dosimetric evaluations of PBS proton treatments, as 4DDC based on a single 4DCT can lead to an underestimation of motion effects.

The impact of organ motion and the appliance of mitigation strategies on the effectiveness of hypoxia-guided proton therapy for non-small cell lung cancer

BACKGROUND AND PURPOSE

To investigate the impact of organ motion on hypoxia-guided proton therapy treatments for non-small cell lung cancer (NSCLC) patients.

MATERIALS AND METHODS

Hypoxia PET and 4D imaging data of six NSCLC patients were used to simulate hypoxia-guided proton therapy with different motion mitigation strategies including rescanning, breath-hold, respiratory gating and tumour tracking. Motion-induced dose degradation was estimated for treatment plans with dose painting of hypoxic tumour sub-volumes at escalated dose levels. Tumour control probability (TCP) and dosimetry indices were assessed to weigh the clinical benefit of dose escalation and motion mitigation. In addition, the difference in normal tissue complication probability (NTCP) between escalated proton and photon VMAT treatments has been assessed.

RESULTS

Motion-induced dose degradation was found for target coverage (CTV V_95% up to -4%) and quality of the dose-escalation-by-contour (Q_RMS up to 6%) as a function of motion amplitude and amount of dose escalation. The TCP benefit coming from dose escalation (+4-13%) outweighs the motion-induced losses (<2%). Significant average NTCP reductions of dose-escalated proton plans were found for lungs (-14%), oesophagus (-10%) and heart (-16%) compared to conventional VMAT plans. The best plan dosimetry was obtained with breath hold and respiratory gating with rescanning.

CONCLUSION

NSCLC affected by hypoxia appears to be a prime target for proton therapy which, by dose-escalation, allows to mitigate hypoxia-induced radio-resistance despite the sensitivity to organ motion. Furthermore, substantial reduction in normal tissue toxicity can be expected compared to conventional VMAT. Accessibility and standardization of hypoxia imaging and clinical trials are necessary to confirm these findings in a clinical setting.

Design of a novel 5-camera surface guidance system with multiple imaging isocenters

Purpose/objective(s)

Surface‐guided radiation therapy (SGRT) can track the patient surface noninvasively to complement radiographic image‐guided radiation therapy with a standard 3‐camera system and a single radiation/image isocenter. Here we report the commissioning of a novel SGRT system that monitors three imaging isocenters locations in a proton half‐gantry room with a unique 5‐camera configuration.

Materials/methods

The proton half‐gantry room has three image isocenters, designated ISO‐0, ISO‐1, and ISO‐2, to cover various anatomical sites via a robotic ceiling‐mounted cone‐beam CT. Although ISO‐0 and ISO‐1 are used to image the cranium, head and neck, and thoracic regions, ISO‐2 is often used to image body and extremity sites and contiguous craniospinal target volumes. The five‐camera system was calibrated to the radiographic isocenter by using a stereotactic radiosurgery cube phantom for each image isocenter.

Results

The performance of this 5‐camera system was evaluated for 6 degrees of freedom in three categories: (1) absolute setup accuracy relative to the radiographic kV image isocenter based on the DICOM reference; (2) relative shift accuracy based on a reference surface capture; and (3) isocenter tracking accuracy from one isocenter to another based on a reference surface capture. The evaluation revealed maximum deviations of 0.8, 0.2, and 0.6 mm in translation and 0.2°, 0.1°, and 0.1° in rotation for the first, second, and third categories, respectively. Comparing the dosimetry and latency with static and gated irradiation revealed a 0.1% dose difference and positional differences of 0.8 mm in X and 0.9 mm in Y with less than 50 ms temporal accuracy.

Conclusion

The unique 5‐camera system configuration provides SGRT at the treatment isocenter (ISO‐0) and also imaging isocenter locations (ISO‐0, ISO‐1, and ISO‐2) to ensure correct patient positioning before and after radiographic imaging, especially during transitions from the offset imaging isocenters to the treatment isocenter.

Experimental validation of 4D log file-based proton dose reconstruction for interplay assessment considering amplitude-sorted 4DCTs

PURPOSE

The unpredictable interplay between dynamic proton therapy delivery and target motion in the thorax can lead to severe dose distortions. A fraction-wise four-dimensional (4D) dose reconstruction workflow allows for the assessment of the applied dose after patient treatment while considering the actual beam delivery sequence extracted from machine log files, the recorded breathing pattern and the geometric information from a 4D computed tomography scan (4DCT). Such an algorithm capable of accounting for amplitude-sorted 4DCTs was implemented and its accuracy as well as its sensitivity to input parameter variations was experimentally evaluated.

METHODS

An anthropomorphic thorax phantom with a movable insert containing a target surrogate and a radiochromic film was irradiated with a monoenergetic field for various 1D target motion forms (sin, sin4) and peak-to-peak amplitudes (5/10/15/20/30 mm). The measured characteristic film dose distributions were compared to the respective sections in the 4D reconstructed doses using a 2D γ-analysis (3mm, 3%); γ-pass rates were derived for different dose grid resolutions (1mm/3mm) and deformable image registrations (DIR, automatic/manual) applied during the 4D dose reconstruction process. In an additional analysis, the sensitivity of reconstructed dose distributions against potential asynchronous timing of the motion and machine log files was investigated for both a monoenergetic field and more realistic 4D robustly optimized fields by artificially introduced offsets of ± 1/5/25/50/250 ms. The resulting dose distributions with asynchronized log files were compared to those with synchronized log files by means of a 3D γ-analysis (1mm, 1%) and the evaluation of absolute dose differences.

RESULTS

The induced characteristic interplay patterns on the films were well reproduced by the 4D dose reconstruction with 2D γ-pass rates ≥95% for almost all cases with motion magnitudes ≤15 mm. In general, the 2D γ-pass rates showed a significant decrease for larger motion amplitudes and increase when using a finer dose grid resolution but were not affected by the choice of motion form (sin, sin4). There was also a trend, though not statistically significant, towards the manually defined DIR for better quality of the reconstructed dose distributions in the area imaged by the film. The 4D dose reconstruction results for the monoenergetic as well as the 4D robustly optimized fields were robust against small asynchronies between motion and machine log files of up to 5 ms, which is in the order of potential network latencies.

CONCLUSIONS

We have implemented a 4D log file-based proton dose reconstruction that accounts for amplitude-sorted 4DCTs. Its accuracy was proven to be clinically acceptable for target motion magnitudes of up to 15 mm. Particular attention should be paid to the synchronization of the log file generating systems as the reconstructed dose distribution may vary with log file asynchronies larger than those caused by realistic network delays. This article is protected by copyright. All rights reserved.

Report of AAPM Task Group 290: Respiratory motion management for particle therapy

Dose uncertainty induced by respiratory motion remains a major concern for treating thoracic and abdominal lesions using particle beams. This Task Group report reviews the impact of tumor motion and dosimetric considerations in particle radiotherapy, current motion‐management techniques, and limitations for different particle‐beam delivery modes (i.e., passive scattering, uniform scanning, and pencil‐beam scanning). Furthermore, the report provides guidance and risk analysis for quality assurance of the motion‐management procedures to ensure consistency and accuracy, and discusses future development and emerging motion‐management strategies. This report supplements previously published AAPM report TG76, and considers aspects of motion management that are crucial to the accurate and safe delivery of particle‐beam therapy. To that end, this report produces general recommendations for commissioning and facility‐specific dosimetric characterization, motion assessment, treatment planning, active and passive motion‐management techniques, image guidance and related decision‐making, monitoring throughout therapy, and recommendations for vendors. Key among these recommendations are that: (1) facilities should perform thorough planning studies (using retrospective data) and develop standard operating procedures that address all aspects of therapy for any treatment site involving respiratory motion; (2) a risk‐based methodology should be adopted for quality management and ongoing process improvement.

Small spot size versus large spot size: Effect on plan quality for lung cancer in pencil beam scanning proton therapy

Purpose

The purpose of the current study was to evaluate the impact of spot size on the interplay effect, plan robustness, and dose to the organs at risk for lung cancer plans in pencil beam scanning (PBS) proton therapy

Methods

The current retrospective study included 13 lung cancer patients. For each patient, small spot (∼3 mm) plans and large spot (∼8 mm) plans were generated. The Monte Carlo algorithm was used for both robust plan optimization and final dose calculations. Each plan was normalized, such that 99% of the clinical target volume (CTV) received 99% of the prescription dose. Interplay effect was evaluated for treatment delivery starting in two different breathing phases (T0 and T50). Plan robustness was investigated for 12 perturbed scenarios, which combined the isocenter shift and range uncertainty. The nominal and worst‐case scenario (WCS) results were recorded for each treatment plan. Equivalent uniform dose (EUD) and normal tissue complication probability (NTCP) were evaluated for the total lung, heart, and esophagus.

Results

In comparison to large spot plans, the WCS values of small spot plans at CTV D_95%, D_96%, D_97%, D_98%, and D_99% were higher with the average differences of 2.2% (range, 0.3%–3.7%), 2.3% (range, 0.5%–4.0%), 2.6% (range, 0.6%–4.4%), 2.7% (range, 0.9%–5.2%), and 2.7% (range, 0.3%–6.0%), respectively. The nominal and WCS mean dose and EUD for the esophagus, heart, and total lung were higher in large spot plans. The difference in NTCP between large spot and small spot plans was up to 1.9% for the total lung, up to 0.3% for the heart, and up to 32.8% for the esophagus. For robustness acceptance criteria of CTV D_95% ≥ 98% of the prescription dose, seven small spot plans had all 12 perturbed scenarios meeting the criteria, whereas, for 13 large spot plans, there were ≥2 scenarios failing to meet the criteria. Interplay results showed that, on average, the target coverage in large spot plans was higher by 1.5% and 0.4% in non‐volumetric and volumetric repainting plans, respectively.

Conclusion

For robustly optimized PBS lung cancer plans in our study, a small spot machine resulted in a more robust CTV against the setup and range errors when compared to a large spot machine. In the absence of volumetric repainting, large spot PBS lung plans were more robust against the interplay effect. The use of a volumetric repainting technique in both small and large spot PBS lung plans led to comparable interplay target coverage.

Future Developments in Charged Particle Therapy: Improving Beam Delivery for Efficiency and Efficacy

The physical and clinical benefits of charged particle therapy (CPT) are well recognized. However, the availability of CPT and complete exploitation of dosimetric advantages are still limited by high facility costs and technological challenges. There are extensive ongoing efforts to improve upon these, which will lead to greater accessibility, superior delivery, and therefore better treatment outcomes. Yet, the issue of cost remains a primary hurdle as utility of CPT is largely driven by the affordability, complexity and performance of current technology. Modern delivery techniques are necessary but limited by extended treatment times. Several of these aspects can be addressed by developments in the beam delivery system (BDS) which determines the overall shaping and timing capabilities enabling high quality treatments. The energy layer switching time (ELST) is a limiting constraint of the BDS and a determinant of the beam delivery time (BDT), along with the accelerator and other factors. This review evaluates the delivery process in detail, presenting the limitations and developments for the BDS and related accelerator technology, toward decreasing the BDT. As extended BDT impacts motion and has dosimetric implications for treatment, we discuss avenues to minimize the ELST and overview the clinical benefits and feasibility of a large energy acceptance BDS. These developments support the possibility of advanced modalities and faster delivery for a greater range of treatment indications which could also further reduce costs. Further work to realize methodologies such as volumetric rescanning, FLASH, arc, multi-ion and online image guided therapies are discussed. In this review we examine how increased treatment efficiency and efficacy could be achieved with improvements in beam delivery and how this could lead to faster and higher quality treatments for the future of CPT.

A new emittance selection system to maximize beam transmission for low-energy beams in cyclotron-based proton therapy facilities with gantry

Purpose

In proton therapy, the potential of using high‐dose rates in the cancer treatment is being explored. High‐dose rates could improve efficiency and throughput in standard clinical practice, allow efficient utilization of motion mitigation techniques for moving targets, and potentially enhance normal tissue sparing due to the so‐called FLASH effect. However, high‐dose rates are difficult to reach when lower energy beams are applied in cyclotron‐based proton therapy facilities, because they result in large beam sizes and divergences downstream of the degrader, incurring large losses from the cyclotron to the patient position (isocenter). In current facilities, the emittance after the degrader is reduced using circular collimators; however, this does not provide an optimal matching to the acceptance of the following beamline, causing a low transmission for these energies. We, therefore, propose to use a collimation system, asymmetric in both beam size and divergence, resulting in symmetric emittance in both beam transverse planes as required for a gantry system. This new emittance selection, together with a new optics design for the following beamline and gantry, allows a better matching to the beamline acceptance and an improvement of the transmission.

Methods

We implemented a custom method to design the collimator sizes and shape required to select high emittance, to be transported by the following beamline using new beam optics (designed with TRANSPORT) to maximize acceptance matching. For predicting the transmission in the new configuration (new collimators + optics), we used Monte Carlo simulations implemented in BDSIM, implementing a model of PSI Gantry 2 which we benchmarked against measurements taken in the current clinical scenario (circular collimators + clinical optics).

Results

From the BDSIM simulations, we found that the new collimator system and matching beam optics results in an overall transmission from the cyclotron to the isocenter for a 70 MeV beam of 0.72%. This is an improvement of almost a factor of 6 over the current clinical performance (0.13% transmission). The new optics satisfies clinical beam requirements at the isocenter.

Conclusions

We developed a new emittance collimation system for PSI's PROSCAN beamline which, by carefully selecting beam size and divergence asymmetrically, increases the beam transmission for low‐energy beams in current state‐of‐the‐art cyclotron‐based proton therapy gantries. With these improvements, we could predict almost 1% transmission for low‐energy beams at PSI's Gantry 2. Such a system could easily be implemented in facilities interested in increasing dose rates for efficient motion mitigation and FLASH experiments alike.

Diaphragm-Based Position Verification to Improve Daily Target Dose Coverage in Proton and Photon Radiation Therapy Treatment of Distal Esophageal Cancer

PURPOSE

In modern conformal radiation therapy of distal esophageal cancer, target coverage can be affected by variations in the diaphragm position. We investigated if daily position verification (PV) extended by a diaphragm position correction would optimize target dose coverage for esophageal cancer treatment.

METHODS AND MATERIALS

For 15 esophageal cancer patients, intensity modulated proton therapy (IMPT) and volumetric modulated arc therapy (VMAT) plans were computed. Displacements of the target volume were correlated with diaphragm displacements using repeated 4-dimensional computed tomography images to determine the correction needed to account for diaphragm variations. Afterwards, target coverage was evaluated for 3 PV approaches based on: (1) bony anatomy (PV_B), (2) bony anatomy corrected for the diaphragm position (PV_BD) and (3) target volume (PV_T).

RESULTS

The cranial-caudal mean target displacement was congruent with almost half of the diaphragm displacement (y = 0.459x), which was used for the diaphragm correction in PV_BD. Target dose coverage using PV_B was adequate for most patients with diaphragm displacements up till 10 mm (≥94% of the dose in 98% of the volume [D_98%]). For larger displacements, the target coverage was better maintained by PV_T and PV_BD. Overall, PV_BD accounted best for target displacements, especially in combination with tissue density variations (D_98%: IMPT 94% ± 5%, VMAT 96% ± 5%). Diaphragm displacements of more than 10 mm were observed in 22% of the cases.

CONCLUSIONS

PV_B was sufficient to achieve adequate target dose coverage in case of small deviations in diaphragm position. However, large deviations of the diaphragm were best mitigated by PV_BD. To detect the cases where target dose coverage could be compromised due to diaphragm position variations, we recommend monitoring of the diaphragm position before treatment through online imaging.

High-dose hypofractionated pencil beam scanning carbon ion radiotherapy for lung tumors: Dosimetric impact of different spot sizes and robustness to interfractional uncertainties

PURPOSE

The robustness against setup and motion uncertainties of gated four-dimensional restricted robust optimization (4DRRO) was investigated for hypofractionated carbon ion radiotherapy (CIRT) of lung tumors.

METHODS

CIRT plans of 9 patients were optimized using 4DRRO strategy with 3 mm setup errors, 3% density errors and 3 breathing phases related to the gate window. The prescription was 60 Gy(RBE) in 4 fractions. Standard spots (SS) were compared to big spots (BS). Plans were recalculated on multiple 4DCTs acquired within 3 weeks from treatment simulation and rigidly registered with planning images using bone matching. Warped dose distributions were generated using deformable image registration and accumulated on the planning 4DCTs. Target coverage (D98%, D95% and V95%) and dose to lung were evaluated in the recalculated and accumulated dose distributions.

RESULTS

Comparable target coverage was obtained with both spot sizes (p = 0.53 for D95%). The mean lung dose increased of 0.6 Gy(RBE) with BS (p = 0.0078), still respecting the dose constraint of a 4-fraction stereotactic treatment for the risk of radiation pneumonitis. Statistically significant differences were found in the recalculated and accumulated D95% (p = 0.048 and p = 0.024), with BS showing to be more robust. Using BS, the average degradations of the D98%, D95% and V95% in the accumulated doses were -2.7%, -1.6% and -1.5%.

CONCLUSIONS

Gated 4DRRO was highly robust against setup and motion uncertainties. BS increased the dose to healthy tissues but were more robust than SS. The selected optimization settings guaranteed adequate target coverage during the simulated treatment course with acceptable risk of toxicity.

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.

Physics of Particle Beam and Hypofractionated Beam Delivery in NSCLC

The dosimetric advantages of particle therapy lead to significantly reduced integral dose to normal tissues, making it an attractive treatment option for body sites such as the thorax. With reduced normal tissue dose comes the potential for dose escalation, toxicity reduction, or hypofractionation. While proton and heavy ion therapy have been used extensively for NSCLC, there are challenges in planning and delivery compared with X-ray-based radiation therapy. Particularly, range uncertainties compounded by breathing motion have to be considered. This article summarizes the current state of particle therapy for NSCLC with a specific focus on the impact of dosimetric uncertainties in planning and delivery.

Dosimetric influence of deformable image registration uncertainties on propagated structures for online daily adaptive proton therapy of lung cancer patients

PURPOSE

A major burden of introducing an online daily adaptive proton therapy (DAPT) workflow is the time and resources needed to correct the daily propagated contours. In this study, we evaluated the dosimetric impact of neglecting the online correction of the propagated contours in a DAPT workflow.

MATERIAL AND METHODS

For five NSCLC patients with nine repeated deep-inspiration breath-hold CTs, proton therapy plans were optimised on the planning CT to deliver 60 Gy-RBE in 30 fractions. All repeated CTs were registered with six different clinically used deformable image registration (DIR) algorithms to the corresponding planning CT. Structures were propagated rigidly and with each DIR algorithm and reference structures were contoured on each repeated CT. DAPT plans were optimised with the uncorrected, propagated structures (propagated DAPT doses) and on the reference structures (ideal DAPT doses), non-adapted doses were recalculated on all repeated CTs.

RESULTS

Due to anatomical changes occurring during the therapy, the clinical target volume (CTV) coverage of the non-adapted doses reduces on average by 9.7% (V95) compared to an ideal DAPT doses. For the propagated DAPT doses, the CTV coverage was always restored (average differences in the CTV V95 < 1% compared to the ideal DAPT doses). Hotspots were always reduced with any DAPT approach.

CONCLUSION

For the patients presented here, a benefit of online DAPT was shown, even if the daily optimisation is based on propagated structures with some residual uncertainties. However, a careful (offline) structure review is necessary and corrections can be included in an offline adaption.

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 volumetric repainting to mitigate interplay effect on 4D robustly optimized lung cancer plans in pencil beam scanning proton therapy

Purpose

The interplay effect between dynamic pencil proton beams and motion of the lung tumor presents a challenge in treating lung cancer patients in pencil beam scanning (PBS) proton therapy. The main purpose of the current study was to investigate the interplay effect on the volumetric repainting lung plans with beam delivery in alternating order (“down” and “up” directions), and explore the number of volumetric repaintings needed to achieve acceptable lung cancer PBS proton plan.

Method

The current retrospective study included ten lung cancer patients. The total dose prescription to the clinical target volume (CTV) was 70 Gy(RBE) with a fractional dose of 2 Gy(RBE). All treatment plans were robustly optimized on all ten phases in the 4DCT data set. The Monte Carlo algorithm was used for the 4D robust optimization, as well as for the final dose calculation. The interplay effect was evaluated for both the nominal (i.e., without repainting) as well as volumetric repainting plans. The interplay evaluation was carried out for each of the ten different phases as the starting phases. Several dosimetric metrics were included to evaluate the worst‐case scenario (WCS) and bandwidth based on the results obtained from treatment delivery starting in ten different breathing phases.

Results

The number of repaintings needed to meet the criteria 1 (CR1) of target coverage (D_95% ≥ 98% and D_99% ≥ 97%) ranged from 2 to 10. The number of repaintings needed to meet the CR1 of maximum dose (ΔD_1% < 1.5%) ranged from 2 to 7. Similarly, the number of repaintings needed to meet CR1 of homogeneity index (ΔHI < 0.03) ranged from 3 to 10. For the target coverage region, the number of repaintings needed to meet CR1 of bandwidth (<100 cGy) ranged from 3 to 10, whereas for the high‐dose region, the number of repaintings needed to meet CR1 of bandwidth (<100 cGy) ranged from 1 to 7. Based on the overall plan evaluation criteria proposed in the current study, acceptable plans were achieved for nine patients, whereas one patient had acceptable plan with a minor deviation.

Conclusion

The number of repaintings required to mitigate the interplay effect in PBS lung cancer (tumor motion < 15 mm) was found to be highly patient dependent. For the volumetric repainting with an alternating order, a patient‐specific interplay evaluation strategy must be adopted. Determining the optimal number of repaintings based on the bandwidth and WCS approach could mitigate the interplay effect in PBS lung cancer treatment.

Pro-con of proton: Dosimetric advantages of intensity-modulation over passive scatter for thoracic malignancies

•

Intensity Modulated Proton Therapy (IMPT) results in significant reduction of dose to organ at risk.

•

Improving plan robustness mitigates interplay effects.

•

Blanket use of small spots on a group of patients may severely worsen interplay in selected patients.

•

Hypofractionated regimes have fewer interplay effects in both fractional and overall simulations.

•

Randomised control trials are required before any clinical benefit of IMPT can be confirmed.

The use of passively scattered proton therapy (PSPT) or intensity modulated proton therapy (IMPT) opens the potential for dose escalation or critical structure sparing in thoracic malignancies. While the latter offers greater dose conformality, dose distributions are subjected to greater uncertainties, especially due to interplay effects. Exploration in this area is warranted to determine if there is any dosimetric advantages in using IMPT for thoracic malignancies. This review aims to both compare organs-at-risk sparing and plan robustness between PSPT and IMPT and examine the mitigation strategies for the reduction of interplay effects currently available. Early evidence suggests that IMPT is dosimetrically superior to PSPT in thoracic malignancies. Randomised control trials are required before any clinical benefit of IMPT can be confirmed.

Investigation of inter-fraction target motion variations in the context of pencil beam scanned proton therapy in non-small cell lung cancer patients

Purpose

For locally advanced‐stage non‐small cell lung cancer (NSCLC), inter‐fraction target motion variations during the whole time span of a fractionated treatment course are assessed in a large and representative patient cohort. The primary objective is to develop a suitable motion monitoring strategy for pencil beam scanning proton therapy (PBS‐PT) treatments of NSCLC patients during free breathing.

Methods

Weekly 4D computed tomography (4DCT; 41 patients) and daily 4D cone beam computed tomography (4DCBCT; 10 of 41 patients) scans were analyzed for a fully fractionated treatment course. Gross tumor volumes (GTVs) were contoured and the 3D displacement vectors of the centroid positions were compared for all scans. Furthermore, motion amplitude variations in different lung segments were statistically analyzed. The dosimetric impact of target motion variations and target motion assessment was investigated in exemplary patient cases.

Results

The median observed centroid motion was 3.4 mm (range: 0.2–12.4 mm) with an average variation of 2.2 mm (range: 0.1–8.8 mm). Ten of 32 patients (31.3%) with an initial motion <5 mm increased beyond a 5‐mm motion amplitude during the treatment course. Motion observed in the 4DCBCT scans deviated on average 1.5 mm (range: 0.0–6.0 mm) from the motion observed in the 4DCTs. Larger motion variations for one example patient compromised treatment plan robustness while no dosimetric influence was seen due to motion assessment biases in another example case.

Conclusions

Target motion variations were investigated during the course of radiotherapy for NSCLC patients. Patients with initial GTV motion amplitudes of < 2 mm can be assumed to be stable in motion during the treatment course. For treatments of NSCLC patients who exhibit motion amplitudes of > 2 mm, 4DCBCT should be considered for motion monitoring due to substantial motion variations observed.

Consensus Statement on Proton Therapy in Mesothelioma

PURPOSE

Radiation therapy for mesothelioma remains challenging, as normal tissue toxicity limits the amount of radiation that can be safely delivered to the pleural surfaces, especially radiation dose to the contralateral lung. The physical properties of proton therapy result in better sparing of normal tissues when treating the pleura, both in the postpneumonectomy setting and the lung-intact setting. Compared with photon radiation, there are dramatic reductions in dose to the contralateral lung, heart, liver, kidneys, and stomach. However, the tissue heterogeneity in the thorax, organ motion, and potential for changing anatomy during the treatment course all present challenges to optimal irradiation with protons.

METHODS

The clinical data underlying proton therapy in mesothelioma are reviewed here, including indications, advantages, and limitations.

RESULTS

The Particle Therapy Cooperative Group Thoracic Subcommittee task group provides specific guidelines for the use of proton therapy for mesothelioma.

CONCLUSIONS

This consensus report can be used to guide clinical practice, insurance approval, and future research.

Scanning Beam Proton Therapy versus Photon IMRT for Stage III Lung Cancer: Comparison of Dosimetry, Toxicity, and Outcomes

Purpose

There are limited clinical data on scanning-beam proton therapy (SPT) in treating locally advanced lung cancer, as most published studies have used passive-scatter technology. There is increasing interest in whether the dosimetric advantages of SPT compared with photon therapy can translate into superior clinical outcomes. We present our experience of SPT and photon intensity modulated radiation therapy (IMRT) with clinical dosimetry and outcomes in patients with stage III lung cancer.

Methods and Materials

Patients with stage III lung cancer treated at our center between 2013 and May 2018 were identified in compliance with our institutional review board (64 patients = 34 SPT + 30 IMRT). Most proton patients were treated with pencil beam scanning (28 of 34), and 6 of 34 were treated with uniform scanning. Fisher exact test, χ² test, and Mann-Whitney test were used to compare groups. All tests were 2-sided.

Results

Patient characteristics were similar between the IMRT and SPT patients, except for worse lung function in the IMRT group. Mean dose to lung, heart, and esophagus was lower in the SPT group, with most benefit in the low-dose region (lungs, 9.7 Gy vs 15.7 Gy for SPT vs IMRT, respectively [P = .004]; heart, 7 Gy vs 14 Gy [P = .001]; esophagus, 28.2 Gy vs 30.9 Gy [P = .023]). Esophagitis and dermatitis grades were not different between the 2 groups. Grade 2+ pneumonitis was 21% in the SPT group and 40% in the IMRT group (P = .107). Changes in blood counts were not different between the 2 groups. Overall survival and progression-free survival were not different between SPT and IMRT (median overall survival, 41.6 vs 30.7 months, respectively [P = .52]; median progression-free survival, 19.5 vs 14.6 months [P = .50]).

Conclusions

We report our experience with SPT and IMRT in stage III lung cancer. Our cohort of patients treated with SPT had lower doses to normal organs (lungs, heart, esophagus) than our IMRT cohort. There was no statistically significant difference in toxicity rates or survival, although there may have been a trend toward lower rates of pneumonitis.

Current status of proton therapy techniques for lung cancer

Proton beams have been used for cancer treatment for more than 28 years, and several technological advancements have been made to achieve improved clinical outcomes by delivering more accurate and conformal doses to the target cancer cells while minimizing the dose to normal tissues. The state-of-the-art intensity modulated proton therapy is now prevailing as a major treatment technique in proton facilities worldwide, but still faces many challenges in being applied to the lung. Thus, in this article, the current status of proton therapy technique is reviewed and issues regarding the relevant uncertainty in proton therapy in the lung are summarized.

Comparative photon and proton dosimetry for patients with mediastinal lymphoma in the era of Monte Carlo treatment planning and variable relative biological effectiveness

Background

Existing pencil beam analytical (PBA) algorithms for proton therapy treatment planning are not ideal for sites with heterogeneous tissue density and do not account for the spatial variations in proton relative biological effectiveness (vRBE). Using a commercially available Monte Carlo (MC) treatment planning system, we compared various dosimetric endpoints between proton PBA, proton MC, and photon treatment plans among patients with mediastinal lymphoma.

Methods

Eight mediastinal lymphoma patients with both free breathing (FB) and deep inspiration breath hold (DIBH) CT simulation scans were analyzed. The original PBA plans were re-calculated with MC. New proton plans that used MC for both optimization and dose calculation with equivalent CTV/ITV coverage were also created. A vRBE model, which uses a published model for DNA double strand break (DSB) induction, was applied on MC plans to study the potential impact of vRBE on cardiac doses. Comparative photon plans were generated on the DIBH scan.

Results

Re-calculation of FB PBA plans with MC demonstrated significant under coverage of the ITV V99 and V95. Target coverage was recovered by re-optimizing the PT plan with MC with minimal change to OAR doses. Compared to photons with DIBH, MC-optimized FB and DIBH proton plans had significantly lower dose to the mean lung, lung V5, breast tissue, and spinal cord for similar target coverage. Even with application of vRBE in the proton plans, the putative increase in RBE at the end of range did not decrease the dosimetric advantages of proton therapy in cardiac substructures.

Conclusions

MC should be used for PT treatment planning of mediastinal lymphoma to ensure adequate coverage of target volumes. Our preliminary data suggests that MC-optimized PT plans have better sparing of the lung and breast tissue compared to photons. Also, the potential for end of range RBE effects are unlikely to be large enough to offset the dosimetric advantages of proton therapy in cardiac substructures for mediastinal targets, although these dosimetric findings require validation with late toxicity data.

Proton pencil beam scanning treatment with feedback based voluntary moderate breath hold

Introduction The aim of this article is to introduce a novel protocol for proton pencil beam scanning treatment with moderate deep inspiration breath hold (mDIBH) and report on our clinical implementation results. Methods Three computed tomography (CT) scannings to build the patient's anatomy model were performed during the patient's voluntary mDIBH. All 3 CT scans were used in the optimization during the treatment planning process. Both orthogonal kV imaging and cone-beam computed tomography (CBCT) were implemented for patient alignment with BH prior to the treatment. The BH CBCT images were analyzed for BH reproducibility and the virtual total dose (VTD) retrospectively. To find the VTD, a series of deformable image registrations (DIR) were performed between CBCT and pCT. The effect of the variation of lung density on the dose distribution was also analyzed in the study. Results The values of the mean, standard deviation, maximum, and minimum of the tumor location difference between the CBCT and pCT were 1.9, 1.6, 4.7, and 0.0 mm, respectively. The percentage difference in D99% of CTVs between VTD and the nominal plan was within 1.5%. Conclusions The feedback-based voluntary moderate BH proton PBS treatment was successfully performed in our clinic. This study shows that there is a potential to implement the BH treatment widely in proton centers.

An analytical model for the upper bound estimation of respiratory motion-induced dose uncertainty in spot-scanning proton beam therapy

PURPOSE

We developed an analytical model of a spot-scanning beam delivery system to estimate the upper bound of respiratory motion-induced dose uncertainty for a given treatment plan.

METHODS

The effective delivery time for each spot position in the treatment plan was calculated on the basis of the parameters of the delivery system. The upper bound of the dose uncertainty was then calculated as a function of the effective delivery time. Two-dimensional (2D) measurements with a detector array on a one-dimensional moving platform were obtained to validate the model.

RESULTS

We performed 351 two-dimensional measurements on a moving platform for different delivery sequences of a single-layer uniform pattern and patient treatment field. The measured dose uncertainty was a strong function of the effective delivery time: The shortest effective delivery time resulted in a maximum absolute dose error of >90%, while the longest ones resulted in a maximum absolute dose error of 4.9% for a single layer and 9.7% for a patient field with heterogeneity. The relationship of the effective delivery time and the measured dose uncertainty followed the analytical formula.

CONCLUSIONS

With our analytical model, the upper bound of the dose uncertainty due to motion can be estimated in spot-scanning proton therapy without four-dimensional dynamic dose calculation.

Impact of planned dose reporting methods on Gamma pass rates for IROC lung and liver motion phantoms treated with pencil beam scanning protons

PURPOSE

The purpose of this study is to evaluate the impact of two methods of reporting planned dose distributions on the Gamma analysis pass rates for comparison with measured 2D film dose and simulated delivered 3D dose for proton pencil beam scanning treatment of the Imaging and Radiation Oncology Core (IROC) proton lung and liver mobile phantoms.

METHODS AND MATERIALS

Four-dimensional (4D) computed-tomography (CT) image sets were acquired for IROC proton lung and liver mobile phantoms, which include dosimetry inserts that contains targets, thermoluminescent dosimeters and EBT2 films for plan dose verification. 4DCT measured fixed motion magnitudes were 1.3 and 1.0 cm for the lung and liver phantoms, respectively. To study the effects of motion magnitude on the Gamma analysis pass rate, three motion magnitudes for each phantom were simulated by creating virtual 4DCT image sets with motion magnitudes scaled from the scanned phantom motion by 50, 100, and 200%. The internal target volumes were contoured on the maximum intensity projection CTs of the 4DCTs for the lung phantom and on the minimum intensity projection CTs of the 4DCTs for the liver phantom. Treatment plans were optimized on the average intensity projection (AVE) CTs of the 4DCTs using the RayStation treatment planning system. Plan doses were calculated on the AVE CTs, which was defined as the planned AVE dose (method one). Plan doses were also calculated on all 10 phase CTs of the 4DCTs and were registered using target alignment to and equal-weight-summed on the 50% phase (T50) CT, which was defined as the planned 4D dose (method two). The planned AVE doses and 4D doses for phantom treatment were reported to IROC, and the 2D-2D Gamma analysis pass rates for measured film dose relative to the planned AVE and 4D doses were compared. To evaluate motion interplay effects, simulated delivered doses were calculated for each plan by sorting spots into corresponding respiratory phases using spot delivery time recorded in the log files by the beam delivery system to calculate each phase dose and accumulate dose to the T50 CTs. Ten random beam starting phases were used for each beam to obtain the range of the simulated delivered dose distributions. 3D-3D Gamma analyses were performed to compare the planned 4D/AVE doses with simulated delivered doses.

RESULTS

The planned 4D dose matched better with the measured 2D film dose and simulated delivered 3D dose than the planned AVE dose. Using planned 4D dose as institution reported planned dose to IROC improved IROC film dose 2D-2D Gamma analysis pass rate from 92 to 96% on average for three films for the lung phantom (7% 5 mm), and from 92 to 94% in the sagittal plane for the liver phantom (7% 4 mm), respectively, compared with using the planned AVE dose. The 3D-3D Gamma analysis (3% 3 mm) pass rate showed that the simulated delivered doses for lung and liver phantoms using 10 random beam starting phases for each delivered beam matched the planned 4D dose significantly better than the planned AVE dose for phantom motions larger than 1 cm (p ≤ 0.04).

CONCLUSIONS

It is recommended to use the planned 4D dose as the institution reported planned dose to IROC to compare with the measured film dose for proton mobile phantoms to improve film Gamma analysis pass rate in the IROC credentialing process.

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.

Clinical implementation of respiratory-gated spot-scanning proton therapy: An efficiency analysis of active motion management

Purpose

The aim of this work is to describe the clinical implementation of respiratory‐gated spot‐scanning proton therapy (SSPT) for the treatment of thoracic and abdominal moving targets. The experience of our institution is summarized, from initial acceptance and commissioning tests to the development of standard clinical operating procedures for simulation, motion assessment, motion mitigation, treatment planning, and gated SSPT treatment delivery.

Materials and methods

A custom respiratory gating interface incorporating the Real‐Time Position Management System (RPM, Varian Medical Systems, Inc., Palo Alto, CA, USA) was developed in‐house for our synchrotron‐based delivery system. To assess gating performance, a motion phantom and radiochromic films were used to compare gated vs nongated delivery. Site‐specific treatment planning protocols and conservative motion cutoffs were developed, allowing for free‐breathing (FB), breath‐holding (BH), or phase‐gating (Ph‐G). Room usage efficiency of BH and Ph‐G treatments was retrospectively evaluated using beam delivery data retrieved from our record and verify system and DICOM files from patient‐specific quality assurance (QA) procedures.

Results

More than 70 patients were treated using active motion management between the launch of our motion mitigation program in October 2015 and the end date of data collection of this study in January 2018. During acceptance procedures, we found that overall system latency is clinically‐suitable for Ph‐G. Regarding room usage efficiency, the average number of energy layers delivered per minute was <10 for Ph‐G, 10‐15 for BH and ≥15 for FB, making Ph‐G the slowest treatment modality. When comparing to continuous delivery measured during pretreatment QA procedures, the median values of BH treatment time were extended from 6.6 to 9.3 min (+48%). Ph‐G treatments were extended from 7.3 to 13.0 min (+82%).

Conclusions

Active motion management has been crucial to the overall success of our SSPT program. Nevertheless, our conservative approach has come with an efficiency cost that is more noticeable in Ph‐G treatments and should be considered in decision‐making.

Reproducibility of the lung anatomy under active breathing coordinator control: Dosimetric consequences for scanned proton treatments

PURPOSE

The treatment of moving targets with scanned proton beams is challenging. For motion mitigation, an Active Breathing Coordinator (ABC) can be used to assist breath-holding. The delivery of pencil beam scanning fields often exceeds feasible breath-hold durations, requiring high breath-hold reproducibility. We evaluated the robustness of scanned proton therapy against anatomical uncertainties when treating nonsmall-cell lung cancer (NSCLC) patients during ABC controlled breath-hold.

METHODS

Four subsequent MRIs of five healthy volunteers (3 male, 2 female, age: 25-58, BMI: 19-29) were acquired under ABC controlled breath-hold during two simulated treatment fractions, providing both intrafractional and interfractional information about breath-hold reproducibility. Deformation vector fields between these MRIs were used to deform CTs of five NSCLC patients. Per patient, four or five cases with different tumor locations were modeled, simulating a total of 23 NSCLC patients. Robustly optimized (3 and 5 mm setup uncertainty respectively and 3% density perturbation) intensity-modulated proton plans (IMPT) were created and split into subplans of 20 s duration (assumed breath-hold duration). A fully fractionated treatment was recalculated on the deformed CTs. For each treatment fraction the deformed CTs representing multiple breath-hold geometries were alternated to simulate repeated ABC breath-holding during irradiation. Also a worst-case scenario was simulated by recalculating the complete treatment plan on the deformed CT scan showing the largest deviation with the first deformed CT scan, introducing a systematic error. Both the fractionated breath-hold scenario and worst-case scenario were dosimetrically evaluated.

RESULTS

Looking at the deformation vector fields between the MRIs of the volunteers, up to 8 mm median intra- and interfraction displacements (without outliers) were found for all lung segments. The dosimetric evaluation showed a median difference in D98% between the planned and breath-hold scenarios of -0.1 Gy (range: -4.1 Gy to 2.0 Gy). D98% target coverage was more than 57.0 Gy for 22/23 cases. The D1 cc of the CTV increased for 21/23 simulations, with a median difference of 0.9 Gy (range: -0.3 to 4.6 Gy). For 14/23 simulations the increment was beyond the allowed maximum dose of 63.0 Gy, though remained under 66.0 Gy (110% of the prescribed dose of 60.0 Gy). Organs at risk doses differed little compared to the planned doses (difference in mean doses <0.9 Gy for the heart and lungs, <1.4% difference in V35 [%] and V20 [%] to the esophagus and lung).

CONCLUSIONS

When treating under ABC controlled breath-hold, robustly optimized IMPT plans show limited dosimetric consequences due to anatomical variations between repeated ABC breath-holds for most cases. Thus, the combination of robustly optimized IMPT plans and the delivery under ABC controlled breath-hold presents a safe approach for PBS lung treatments.

Four-dimensional Plan Optimization for the Treatment of Lung Tumors Using Pencil-beam Scanning Proton Radiotherapy

Purpose 

This study aimed to evaluate the effectiveness of four-dimensional (4D) robust optimization for proton pencil-beam scanning (PBS) treatment of lung tumors.

Patients and methods 

In seven patients with lung cancer, proton beam therapy was planned using 4D robust optimization over 4D computed tomography (CT) data sets. The gross target volume (GTV) was contoured based on individual breathing phases, and a 5-mm expansion was used to generate the clinical target volume (CTV) for each phase. The 4D optimization was conducted directly on the 4D CT data set. The robust optimization settings included a CT Hounsfield unit (HU) uncertainty of 4% and a setup uncertainty of 5 mm to obtain the CTV. Additional target dose objectives such as those for the internal target volume (ITV) as well as the organ-at-risk (OAR) dose requirements were placed on the average CT. For comparison, three-dimensional (3D) robust optimization was also performed on the average CT. An additional verification 4D CT was performed to verify plan robustness against inter-fractional variations.

Results 

Target coverages were generally higher for 4D optimized plans. The difference was most pronounced for ITV V70Gy when evaluating individual breathing phases. The 4D optimized plans were shown to be able to maintain the ITV coverage at full prescription, while 3D optimized plans could not. More importantly, this difference in ITV V70Gy between the 4D and 3D optimized plans was also consistently observed when evaluating the verification 4D CT, indicating that the 4D optimized plans were more robust against inter-fractional variations. Less difference was seen between the 4D and 3D optimized plans in the lungs criteria: V5Gy and V20Gy.

Conclusion 

The proton PBS treatment plans optimized directly on the 4D CT were shown to be more robust when compared to those optimized on a regular 3D CT. Robust 4D optimization can improve the target coverage for the proton PBS lung treatments.

4D robust optimization including uncertainties in time structures can reduce the interplay effect in proton pencil beam scanning radiation therapy

PURPOSE

Interplay effects in proton radiotherapy can create large distortions in the dose distribution and severely degrade the plan quality. Standard methods to mitigate these effects include abdominal compression, gating, and rescanning. We propose a new method to include the time structures of the delivery and organ motion in the framework of four-dimensional (4D) robust optimization to generate plans that are robust against interplay effects.

METHODS

The method considers multiple scenarios reflecting the uncertainties in the delivery and in the organ motion. In each scenario, the pencil beam scanning spots are distributed to different phases of the breathing cycle according to each individual spot time stamp, and a partial beam dose is calculated for each phase. The partial beam doses are accumulated on a reference phase through deformable image registrations. Minimax optimization is performed to take all scenarios into account simultaneously. For simplicity, the uncertainties in this proof of concept study are limited to variations in the breathing pattern. The method is evaluated for three different nonsmall cell lung cancer patients and compared to plans using conventional 4D robust optimization both with and without rescanning. We assess the ability of the method to mitigate distortions from the interplay effect over multiple evaluation scenarios using 4D dose calculations. This interplay evaluation is performed in an experimentally validated framework, which is independent of the optimization in the plan generation step.

RESULTS

For the three studied patients, 4D optimization including time structures is efficient, especially for large tumor motions, where rescanning of conventional 4D robustly optimized plans is not sufficient to mitigate the interplay effect. The most efficient approach of the new method is achieved when it is combined with rescanning. For the patient with the largest motion, the mean V95% is 99.2% and mean V107% is 3.65% for the best rescanned 4D plan optimized with time structure. This can be compared to conventional 4D optimized plans with mean V95% of 92.7% and mean V107% of 13.1%.

CONCLUSIONS

The current study shows the potential of reducing interplay effects in proton pencil beam scanning radiotherapy by incorporating organ motion and delivery characteristics in a 4D robust optimization.

Clinical examination of proton pencil beam scanning on a moving anthropomorphic lung phantom

The objective of this study was to examine the use of proton pencil beam scanning for the treatment of moving lung tumors. A single-field uniform dose proton pencil beam scanning (PBS) plan was generated for the standard thorax phantom designed by the Imaging and Radiation Oncology Core (IROC) Houston QA Center. Robust optimization, including range and setup uncertainties as well as volumetric repainting, was used for the plan. Patient-specific quality assurance (QA) measurements were performed using both a water tank and a custom heterogeneous QA phantom. A custom moving phantom was used to find the optimal number of volumetric repainting. Both analytical and Monte Carlo (MC) algorithms were used for dose calculation and their accuracies were compared with actual measurements. A single ionization chamber, a 2-dimensional ionization chamber array, thermoluminescent dosimeters (TLDs), and films were used for dose measurements. The optimal number of volumetric repainting was found to be 4 times in our system. The mean dose overestimations on a moving target by analytical and MC algorithms based on a time-averaged computed tomography (CT) image of the phantom were found to be 4.8% and 2.4%, respectively. The mean gamma indexes for analytical and MC algorithms were 91% and 96%, respectively. The MC dose algorithm calculation was found to have a better agreement with measurements compared with the analytical algorithm. When treating moving lung tumors using proton PBS, the techniques of robust optimization, volumetric repainting, and MC dose calculation were found effective. Extra care needs to be taken when an analytical dose calculation algorithm is used.

Effectiveness of different rescanning techniques for scanned proton radiotherapy in lung cancer patients

Non-small cell lung cancer (NSCLC) is a tumour type thought to be well-suited for proton radiotherapy. However, the lung region poses many problems related to organ motion and can for actively scanned beams induce severe interplay effects. In this study we investigate four mitigating rescanning techniques: (1) volumetric rescanning, (2) layered rescanning, (3) breath-sampled (BS) layered rescanning, and (4) continuous breath-sampled (CBS) layered rescanning. The breath-sampled methods will spread the layer rescans over a full breathing cycle, resulting in an improved averaging effect at the expense of longer treatment times. In CBS, we aim at further improving the averaging by delivering as many rescans as possible within one breathing cycle. The interplay effect was evaluated for 4D robustly optimized treatment plans (with and without rescanning) for seven NSCLC patients in the treatment planning system RayStation. The optimization and final dose calculation used a Monte Carlo dose engine to account for the density heterogeneities in the lung region. A realistic treatment delivery time structure given from the IBA ScanAlgo simulation tool served as basis for the interplay evaluation. Both slow (2.0 s) and fast (0.1 s) energy switching times were simulated. For all seven studied patients, rescanning improves the dose conformity to the target. The general trend is that the breath-sampled techniques are superior to layered and volumetric rescanning with respect to both target coverage and variability in dose to OARs. The spacing between rescans in our breath-sampled techniques is set at planning, based on the average breathing cycle length obtained in conjunction with CT acquisition. For moderately varied breathing cycle lengths between planning and delivery (up to 15%), the breath-sampled techniques still mitigate the interplay effect well. This shows the potential for smooth implementation at the clinic without additional motion monitoring equipment.

Considerations when treating lung cancer with passive scatter or active scanning proton therapy

Lung cancer, due to its poor clinical outcomes and significant toxicity associated with standard photon-based radiation, is a disease site that has the potential to greatly benefit from accurate treatment with proton radiation therapy. The potential of proton therapy is the ability to increase the radiation dose to the tumor while simultaneously decreasing the radiation dose to surrounding healthy tissues. For lung cancer treatment, this could mean significant sparing of the uninvolved healthy lung, which is difficult to achieve with external photon beam therapy, or decreasing the heart dose. In treating lung cancer with proton therapy, some additional considerations need to be made compared to treating patients with external photon beam radiation therapy. These include accounting for the finite range of protons in the patient, understanding temporal effects, potential dose discrepancies and choosing an appropriate treatment planning system for the task. One final consideration is differences between the different available proton therapy delivery systems-passive scattered proton therapy (PSPT) and active scanning proton therapy.

Impact of Spot Size and Spacing on the Quality of Robustly Optimized Intensity Modulated Proton Therapy Plans for Lung Cancer

PURPOSE

To investigate how spot size and spacing affect plan quality, robustness, and interplay effects of robustly optimized intensity modulated proton therapy (IMPT) for lung cancer.

METHODS AND MATERIALS

Two robustly optimized IMPT plans were created for 10 lung cancer patients: first by a large-spot machine with in-air energy-dependent large spot size at isocenter (σ: 6-15 mm) and spacing (1.3 σ), and second by a small-spot machine with in-air energy-dependent small spot size (σ: 2-6 mm) and spacing (5 mm). Both plans were generated by optimizing radiation dose to internal target volume on averaged 4-dimensional computed tomography scans using an in-house-developed IMPT planning system. The dose-volume histograms band method was used to evaluate plan robustness. Dose evaluation software was developed to model time-dependent spot delivery to incorporate interplay effects with randomized starting phases for each field per fraction. Patient anatomy voxels were mapped phase-to-phase via deformable image registration, and doses were scored using in-house-developed software. Dose-volume histogram indices, including internal target volume dose coverage, homogeneity, and organs at risk (OARs) sparing, were compared using the Wilcoxon signed-rank test.

RESULTS

Compared with the large-spot machine, the small-spot machine resulted in significantly lower heart and esophagus mean doses, with comparable target dose coverage, homogeneity, and protection of other OARs. Plan robustness was comparable for targets and most OARs. With interplay effects considered, significantly lower heart and esophagus mean doses with comparable target dose coverage and homogeneity were observed using smaller spots.

CONCLUSIONS

Robust optimization with a small spot-machine significantly improves heart and esophagus sparing, with comparable plan robustness and interplay effects compared with robust optimization with a large-spot machine. A small-spot machine uses a larger number of spots to cover the same tumors compared with a large-spot machine, which gives the planning system more freedom to compensate for the higher sensitivity to uncertainties and interplay effects for lung cancer treatments.

Feasibility of proton pencil beam scanning treatment of free-breathing lung cancer patients

BACKGROUND

The interplay effect might degrade the dose of pencil beam scanning proton therapy to a degree that free-breathing treatment might be impossible without further motion mitigation techniques, which complicate and prolong the treatment. We assessed whether treatment of free-breathing patients without motion mitigation is feasible.

MATERIAL AND METHODS

For 40 lung cancer patients, 4DCT datasets and individual breathing patterns were used to simulate 4D dynamic dose distributions of 3D treatment plans over 33 fractions delivered with an IBA universal nozzle. Evaluation was done by assessing under- and overdosage in the target structure using the parameters V90, V95, V98, D98, D2, V107 and V110. The impact of using beam-specific target volumes and the impact of changes in motion and patient anatomy in control 4DCTs were assessed.

RESULTS

Almost half of the patients had tumour motion amplitudes of less than 5 mm. Under- and overdosage was significantly smaller for patients with tumour motion below 5 mm compared to patients with larger motion (2% vs. 13% average absolute reduction of V95, 2% vs. 8% average increase in V107, p < .01). Simulating a 33-fraction treatment, the dose degradation was reduced but persisted for patients with tumour motion above 5 mm (average ΔV95 of <1% vs. 3%, p < .01). Beam-specific target volumes reduced the dose degradation in a fractionated treatment, but were more relevant for large motion. Repeated 4DCT revealed that changes in tumour motion during treatment might result in unexpected large dose degradations.

CONCLUSION

Tumour motion amplitude is an indicator of dose degradation caused by the interplay effect. Fractionation reduces the dose degradation allowing the unmitigated treatment of patients with small tumour motions of less than 5 mm. The beam-specific target approach improves the dose coverage. The tumour motion and position needs to be assessed during treatment for all patients, to quickly react to possible changes, which might require treatment adaptation.

Density overwrites of internal tumor volumes in intensity modulated proton therapy plans for mobile lung tumors

The purpose of this study was to investigate internal tumor volume density overwrite strategies to minimize intensity modulated proton therapy (IMPT) plan degradation of mobile lung tumors. Four planning paradigms were compared for nine lung cancer patients. Internal gross tumor volume (IGTV) and internal clinical target volume (ICTV) structures were defined encompassing their respective volumes in every 4DCT phase. The paradigms use different planning CT (pCT) created from the average intensity projection (AIP) of the 4DCT, overwriting the density within the IGTV to account for movement. The density overwrites were: (a) constant filling with 100 HU (C100) or (b) 50 HU (C50), (c) maximum intensity projection (MIP) across phases, and (d) water equivalent path length (WEPL) consideration from beam's-eye-view. Plans were created optimizing dose-influence matrices calculated with fast GPU Monte Carlo (MC) simulations in each pCT. Plans were evaluated with MC on the 4DCTs using a model of the beam delivery time structure. Dose accumulation was performed using deformable image registration. Interplay effect was addressed applying 10 times rescanning. Significantly less DVH metrics degradation occurred when using MIP and WEPL approaches. Target coverage ([Formula: see text] Gy(RBE)) was fulfilled in most cases with MIP and WEPL ([Formula: see text] Gy (RBE)), keeping dose heterogeneity low ([Formula: see text] Gy(RBE)). The mean lung dose was kept lowest by the WEPL strategy, as well as the maximum dose to organs at risk (OARs). The impact on dose levels in the heart, spinal cord and esophagus were patient specific. Overall, the WEPL strategy gives the best performance and should be preferred when using a 3D static geometry for lung cancer IMPT treatment planning. Newly available fast MC methods make it possible to handle long simulations based on 4D data sets to perform studies with high accuracy and efficiency, even prior to individual treatment planning.

Multiple-CT optimization of intensity-modulated proton therapy - Is it possible to eliminate adaptive planning?

BACKGROUND AND PURPOSE

We hypothesized that a plan's robustness to anatomical changes can be improved by optimizing with multiple CT scans of a patient. The purpose of this study was to determine whether an intensity modulated proton therapy (IMPT) plan could be developed to meet dose criteria on both planning and adaptive CT plans.

MATERIAL AND METHODS

Eight lung cancer patients who underwent adaptive IMPT were retrospectively selected. Each patient had two CTs: a primary planning CT (PCT) and an adaptive planning CT (ACT), and IMPT plans associated with the scans. PCT and ACT were then used in combination to optimize one plan (MCT plan). The doses to the target and organs at risk from the PCT plan, ACT plan, P-ACT plan (PCT plan calculated on ACT data), and MCT plans calculated on both CTs were compared.

RESULTS

The MCT plan maintained the D_95% on both CTs (mean, 65.99 Gy on PCT and 66.02 Gy on ACT). Target dose coverage on ACT was significantly better with the MCT plan than with the P-ACT plan (p = 0.01). MCT plans had slightly higher lung V₂₀ (0.6%, p = 0.02) than did PCT plans. The various plans showed no statistically significant difference in heart and spinal cord dose.

CONCLUSIONS

The results of this study indicate that an IMPT plan can meet the dose criteria on both PCT and ACT, and that MCT optimization can improve the plan's robustness to anatomical change.

Efficient Interplay Effect Mitigation for Proton Pencil Beam Scanning by Spot-Adapted Layered Repainting Evenly Spread out Over the Full Breathing Cycle

PURPOSE

To develop and implement a practical repainting method for efficient interplay effect mitigation in proton pencil beam scanning (PBS).

METHODS AND MATERIALS

A new flexible repainting scheme with spot-adapted numbers of repainting evenly spread out over the whole breathing cycle (assumed to be 4 seconds) was developed. Twelve fields from 5 thoracic and upper abdominal PBS plans were delivered 3 times using the new repainting scheme to an ion chamber array on a motion stage. One time was static and 2 used 4-second, 3-cm peak-to-peak sinusoidal motion with delivery started at maximum inhalation and maximum exhalation. For comparison, all dose measurements were repeated with no repainting and with 8 repaintings. For each motion experiment, the 3%/3-mm gamma pass rate was calculated using the motion-convolved static dose as the reference. Simulations were first validated with the experiments and then used to extend the study to 0- to 5-cm motion magnitude, 2- to 6-second motion periods, patient-measured liver tumor motion, and 1- to 6-fraction treatments. The effect of the proposed method was evaluated for the 5 clinical cases using 4-dimensional (4D) dose reconstruction in the planning 4D computed tomography scan. The target homogeneity index, HI = (D₂ - D₉₈)/D_mean, of a single-fraction delivery is reported, where D₂ and D₉₈ is the dose delivered to 2% and 98% of the target, respectively, and D_mean is the mean dose.

RESULTS

The gamma pass rates were 59.6% ± 9.7% with no repainting, 76.5% ± 10.8% with 8 repaintings, and 92.4% ± 3.8% with the new repainting scheme. Simulations reproduced the experimental gamma pass rates with a 1.3% root-mean-square error and demonstrated largely improved gamma pass rates with the new repainting scheme for all investigated motion scenarios. One- and two-fraction deliveries with the new repainting scheme had gamma pass rates similar to those of 3-4 and 6-fraction deliveries with 8 repaintings. The mean HI for the 5 clinical cases was 14.2% with no repainting, 13.7% with 8 repaintings, 12.0% with the new repainting scheme, and 11.6% for the 4D dose without interplay effects.

CONCLUSIONS

A novel repainting strategy for efficient interplay effect mitigation was proposed, implemented, and shown to outperform conventional repainting in experiments, simulations, and dose reconstructions. This strategy could allow for safe and more optimal clinical delivery of thoracic and abdominal proton PBS and better facilitate hypofractionated and stereotactic treatments.