Exploratory Study of 4D versus 3D Robust Optimization in Intensity Modulated Proton Therapy for Lung Cancer

Robustness of intensity modulated proton treatment of esophageal cancer for anatomical changes and breathing motion

BACKGROUND AND PURPOSE

In this study, we assessed the robustness of intensity modulated proton therapy (IMPT) in esophageal cancer for anatomical variations during treatment.

METHODS

The first sixty esophageal cancer patients, treated clinically with chemoradiotherapy were included. The treatment planning strategy was based on an internal target volume (ITV) approach, where the ITV was created from the clinical target volumes (CTVs) delineated on all phases of a 4DCT. For optimization, a 3 mm isotropic margin was added to the ITV, combined with robust optimization using 5 mm setup and 3 % range uncertainty. Each patient received weekly repeat CTs (reCTs). Robust plan re-evaluation on all reCTs, and a robust dose summation was performed. To assess the factors influencing ITV coverage, a multivariate linear regression analysis was performed. Additionally, clinical adaptations were evaluated.

RESULTS

The target coverage was adequate (ITV V94%>98 % on the robust voxel-wise minimum dose) on most reCTs (91 %), and on the summed dose in 92 % of patients. Significant predictors for ITV coverage in the multivariate analysis were diaphragm baseline shift and water equivalent depth (WED) of the ITV in the beam direction. Underdosage of the ITV mainly occurred in week 1 and 4, leading to treatment adaptation of eight patients, all on the first reCT.

CONCLUSION

Our IMPT treatment of esophageal cancer is robust for anatomical variations. Adaptation appears to be most effective in the first week of treatment. Diaphragm baseline shifts and WED are predictive factors for ITV underdosage, and should be incorporated in an adaptation protocol.

Robustness evaluation of pencil beam scanning proton therapy treatment planning: A systematic review

Compared to conventional radiotherapy using X-rays, proton therapy, in principle, allows better conformity of the dose distribution to target volumes, at the cost of greater sensitivity to physical, anatomical, and positioning uncertainties. Robust planning, both in terms of plan optimization and evaluation, has gained high visibility in publications on the subject and is part of clinical practice in many centers. However, there is currently no consensus on the methods and parameters to be used for robust optimization or robustness evaluation. We propose to overcome this deficiency by following the modified Delphi consensus method. This method first requires a systematic review of the literature. We performed this review using the PubMed and Web Of Science databases, via two different experts. Potential conflicts were resolved by a third expert. We then explored the different methods before focusing on clinical studies that evaluate robustness on a significant number of patients. Many robustness assessment methods are proposed in the literature. Some are more successful than others and their implementation varies between centers. Moreover, they are not all statistically or mathematically equivalent. The most sophisticated and rigorous methods have seen more limited application due to the difficulty of their implementation and their lack of widespread availability.

NRG Oncology and Particle Therapy Co-Operative Group Patterns of Practice Survey and Consensus Recommendations on Pencil-Beam Scanning Proton Stereotactic Body Radiation Therapy and Hypofractionated Radiation Therapy for Thoracic Malignancies

Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally fractionated PBSPT because of concerns of amplified uncertainties at the larger dose per fraction. The NRG Oncology and Particle Therapy Cooperative Group Thoracic Subcommittee surveyed proton centers in the United States to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Among other points, the recommendations highlight the need for volumetric image guidance and multiple computed tomography-based robust optimization and robustness tools to minimize further the effect of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques.

Quantifying dose uncertainties resulting from cardiorespiratory motion in intensity-modulated proton therapy for cardiac stereotactic body radiotherapy

Background

Cardiac stereotactic body radiotherapy (CSBRT) with photons efficaciously and safely treats cardiovascular arrhythmias. Proton therapy, with its unique physical and radiobiological properties, can offer advantages over traditional photon-based therapies in certain clinical scenarios, particularly pediatric tumors and those in anatomically challenging areas. However, dose uncertainties induced by cardiorespiratory motion are unknown.

Objective

This study investigated the effect of cardiorespiratory motion on intensity-modulated proton therapy (IMPT) and the effectiveness of motion-encompassing methods.

Methods

We retrospectively included 12 patients with refractory arrhythmia who underwent CSBRT with four-dimensional computed tomography (4DCT) and 4D cardiac CT (4DcCT). Proton plans were simulated using an IBA accelerator based on the 4D average CT. The prescription was 25 Gy in a single fraction, with all plans normalized to ensure that 95% of the target volume received the prescribed dose. 4D dose reconstruction was performed to generate 4D accumulated and dynamic doses. Furthermore, dose uncertainties due to the interplay effect of the substrate target and organs at risk (OARs) were assessed. The differences between internal organs at risk volume (IRV) and OARreal (manually contoured on average CT) were compared. In 4D dynamic dose, meeting prescription requirements entails V25 and D95 reaching 95% and 25 Gy, respectively.

Results

The 4D dynamic dose significantly differed from the 3D static dose. The mean V25 and D95 were 89.23% and 24.69 Gy, respectively, in 4DCT and 94.35% and 24.99 Gy, respectively, in 4DcCT. Eleven patients in 4DCT and six in 4DcCT failed to meet the prescription requirements. Critical organs showed varying dose increases. All metrics, except for Dmean and D50, significantly changed in 4DCT; in 4DcCT, only D50 remained unchanged with regards to the target dose uncertainties induced by the interplay effect. The interplay effect was only significant for the Dmax values of several OARs. Generally, respiratory motion caused a more pronounced interplay effect than cardiac pulsation. Neither IRV nor OARreal effectively evaluated the dose discrepancies of the OARs.

Conclusions

Complex cardiorespiratory motion can introduce dose uncertainties during IMPT. Motion-encompassing techniques may mitigate but cannot entirely compensate for the dose discrepancies. Individualized 4D dose assessments are recommended to verify the effectiveness and safety of CSBRT.

A machine learning-based approach to predict energy layer for each field in spot-scanning proton arc therapy for lung cancer: A feasibility study

BACKGROUND

Determining the optimal energy layer (EL) for each field, under considering both dose constraints and delivery efficiency, is crucial to promoting the development of proton arc therapy (PAT) technology.

PURPOSE

This study aimed to explore the feasibility and potential clinical benefits of utilizing machine learning (ML) technique to automatically select EL for each field in PAT plans of lung cancer.

METHODS

Proton Bragg peak position (BPP) was employed to characterize EL. The ground truth BPPs for each field were determined using the modified ELO-SPAT framework. Features in geometric, water-equivalent thicknesses (WET) and beamlet were defined and extracted. By analyzing the relationship between the extracted features and ground truth, a polynomial regression model with L2-norm regularization (Ridge regression) was constructed and trained. The performance of the regression model was reported as an error between the predictions and the ground truth. Besides, the predictions were used to make PAT plans (PAT_PRED). These plans were compared with those using the ground truth BPPs (PAT_TRUTH) and the mid-WET of the target volumes (PAT_MID) in terms of relative biological effectiveness-weighted dose (RWD) distributions. One hundred ten patients with lung cancer, a total of 7920 samples, were enrolled retrospectively, with 5940 cases randomly selected as the training set and the remaining 1980 cases as the testing set. Nine patients (648 samples) were collected additionally to evaluate the regression model in terms of plan quality and robustness.

RESULTS

With regard to the prediction errors, the root mean squared errors and mean absolute errors between the ML-predicted and ground truth BPPs for the testing set were 9.165 and 6.572 mm, respectively, indicating differences of approximately two to three ELs. As for plan quality, the PAT_TRUTH and PAT_PRED plans performed similarly in terms of plan robustness, target coverage and organs at risk (OARs) protection, with differences smaller than 0.5 Gy(RBE). This trend was also observed for dose conformity and uniformity. The PAT_MID plans produced the lowest robustness index and lowest doses to OARs, along with the highest heterogeneity index, indicating better protection for OARs, improved plan robustness, but compromised dose homogeneity. Additionally, for relatively small tumor sizes, the PAT_MID plan demonstrated a notably poor dose conformity index.

CONCLUSIONS

Within this cohort under investigation, our study demonstrated the feasibility of using ML technique to predict ELs for each field, offering a fast (within 2 s) and memory-efficient reduced way to select ELs for PAT plan.

Clinical Viability of an Active Spot Scanning Beam Delivery System With a Newly Developed Carbon-Ion Treatment Planning System

Purpose

Although active spot scanning irradiation technique is theoretically superior to passive-scattered broad beam irradiation with respect to normal tissue sparing, corroborations of the clinical benefit of carbon-ion spot scanning have remained scarce. This study aims to investigate the feasibility and clinical implementation of an active spot scanning beam calculation algorithm in a homemade carbon-ion treatment planning system by comparing it with a conventional passive uniform scanning technique.

Methods and Materials

Carbon-ion plans were initially formulated using spot/uniform scanning methods in 22 participants enrolled in a prospective observational clinical trial. Subsequently, 2 additional plans were designed, resulting in 3 carbon-ion plans for each participant: uniform and spot scanning with miniridge filters of 2 mm and 4 mm, respectively.

Results

The findings revealed no significant differences in dose homogeneity; however, significant differences in dose conformity were found between the active and passive scanning plans. For dose drop-off outside the target volume, the average gradient index values were 1.94 (95% CI, 1.79%-2.09%), 1.87 (95% CI, 1.73%-2.01%), and 3.20 (95% CI, 2.80%-3.61%) for the miniridge filters of 2 mm and 4 mm, and uniform scanning plans, respectively. The pretreatment tumor volume was 124.7 cm3 (range, 54.2-234 cm3), and the average shrinkage observed was 38.4% (95% CI, 17.6%-59.4%). Seven participants experienced grade 1 acute toxicity, and 4 experienced grade 2 acute toxicity. However, none of the patients developed grade 3 acute toxicity.

Conclusions

Increasing evidence suggests that potential clinical advantages of spot scanning delivery underlie its technical characteristics. As one among the few institutions currently using carbon-ion radiation therapy, the investigation also provides promising safety and efficacy outcomes from the initial groups of treated participants, thereby contributing to the established clinical evidence supporting the effectiveness and superiority of carbon-ion 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.

Technical Note: Generalizable and Promptable Artificial Intelligence Model to Augment Clinical Delineation in Radiation Oncology

BACKGROUND

Efficient and accurate delineation of organs at risk (OARs) is a critical procedure for treatment planning and dose evaluation. Deep learning-based auto-segmentation of OARs has shown promising results and is increasingly being used in radiation therapy. However, existing deep learning-based auto-segmentation approaches face two challenges in clinical practice: generalizability and human-AI interaction. A generalizable and promptable auto-segmentation model, which segments OARs of multiple disease sites simultaneously and supports on-the-fly human-AI interaction, can significantly enhance the efficiency of radiation therapy treatment planning.

PURPOSE

Meta's segment anything model (SAM) was proposed as a generalizable and promptable model for next-generation natural image segmentation. We further evaluated the performance of SAM in radiotherapy segmentation.

METHODS

Computed tomography (CT) images of clinical cases from four disease sites at our institute were collected: prostate, lung, gastrointestinal, and head & neck. For each case, we selected the OARs important in radiotherapy treatment planning. We then compared both the Dice coefficients and Jaccard indices derived from three distinct methods: manual delineation (ground truth), automatic segmentation using SAM's 'segment anything' mode, and automatic segmentation using SAM's 'box prompt' mode that implements manual interaction via live prompts during segmentation.

RESULTS

Our results indicate that SAM's segment anything mode can achieve clinically acceptable segmentation results in most OARs with Dice scores higher than 0.7. SAM's box prompt mode further improves Dice scores by 0.1∼0.5. Similar results were observed for Jaccard indices. The results show that SAM performs better for prostate and lung, but worse for gastrointestinal and head & neck. When considering the size of organs and the distinctiveness of their boundaries, SAM shows better performance for large organs with distinct boundaries, such as lung and liver, and worse for smaller organs with less distinct boundaries, like parotid and cochlea.

CONCLUSIONS

Our results demonstrate SAM's robust generalizability with consistent accuracy in automatic segmentation for radiotherapy. Furthermore, the advanced box-prompt method enables the users to augment auto-segmentation interactively and dynamically, leading to patient-specific auto-segmentation in radiation therapy. SAM's generalizability across different disease sites and different modalities makes it feasible to develop a generic auto-segmentation model in radiotherapy.

Beam mask and sliding window-facilitated deep learning-based accurate and efficient dose prediction for pencil beam scanning proton therapy

BACKGROUND

Accurate and efficient dose calculation is essential for on-line adaptive planning in proton therapy. Deep learning (DL) has shown promising dose prediction results in photon therapy. However, there is a scarcity of DL-based dose prediction methods specifically designed for proton therapy. Successful dose prediction method for proton therapy should account for more challenging dose prediction problems in pencil beam scanning proton therapy (PBSPT) due to its sensitivity to heterogeneities.

PURPOSE

To develop a DL-based PBSPT dose prediction workflow with high accuracy and balanced complexity to support on-line adaptive proton therapy clinical decision and subsequent replanning.

METHODS

PBSPT plans of 103 prostate cancer patients (93 for training and the other 10 for independent testing) and 83 lung cancer patients (73 for training and the other 10 for independent testing) previously treated at our institution were included in the study, each with computed tomography scans (CTs), structure sets, and plan doses calculated by the in-house developed Monte-Carlo dose engine (considered as the ground truth in the model training and testing). For the ablation study, we designed three experiments corresponding to the following three methods: (1) Experiment 1, the conventional region of interest (ROI) (composed of targets and organs-at-risk [OARs]) method. (2) Experiment 2, the beam mask (generated by raytracing of proton beams) method to improve proton dose prediction. (3) Experiment 3, the sliding window method for the model to focus on local details to further improve proton dose prediction. A fully connected 3D-Unet was adopted as the backbone. Dose volume histogram (DVH) indices, 3D Gamma passing rates with a criterion of 3%/3 mm/10%, and dice coefficients for the structures enclosed by the iso-dose lines between the predicted and the ground truth doses were used as the evaluation metrics. The calculation time for each proton dose prediction was recorded to evaluate the method's efficiency.

RESULTS

Compared to the conventional ROI method, the beam mask method improved the agreement of DVH indices for both targets and OARs and the sliding window method further improved the agreement of the DVH indices (for lung cancer, CTV D98 absolute deviation: 0.74 ± 0.18 vs. 0.57 ± 0.21 vs. 0.54 ± 0.15 Gy[RBE], ROI vs. beam mask vs. sliding window methods, respectively). For the 3D Gamma passing rates in the target, OARs, and BODY (outside target and OARs), the beam mask method improved the passing rates in these regions and the sliding window method further improved them (for prostate cancer, targets: 96.93% ± 0.53% vs. 98.88% ± 0.49% vs. 99.97% ± 0.07%, BODY: 86.88% ± 0.74% vs. 93.21% ± 0.56% vs. 95.17% ± 0.59%). A similar trend was also observed for the dice coefficients. This trend was especially remarkable for relatively low prescription isodose lines (for lung cancer, 10% isodose line dice: 0.871 ± 0.027 vs. 0.911 ± 0.023 vs. 0.927 ± 0.017). The dose predictions for all the testing cases were completed within 0.25 s.

CONCLUSIONS

An accurate and efficient deep learning-augmented proton dose prediction framework has been developed for PBSPT, which can predict accurate dose distributions not only inside but also outside ROI efficiently. The framework can potentially further reduce the initial planning and adaptive replanning workload in PBSPT.

Proton Pencil-Beam Scanning Stereotactic Body Radiation Therapy and Hypofractionated Radiation Therapy for Thoracic Malignancies: Patterns of Practice Survey and Recommendations for Future Development from NRG Oncology and PTCOG

Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally-fractionated PBSPT due to concerns of amplified uncertainties at the larger dose per fraction. NRG Oncology and Particle Therapy Cooperative Group (PTCOG) Thoracic Subcommittee surveyed US proton centers to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Amongst other points, the recommendations highlight the need for volumetric image guidance and multiple CT-based robust optimization and robustness tools to minimize further the impact of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques.

Modelling small block aperture in an in-house developed GPU-accelerated Monte Carlo-based dose engine for pencil beam scanning proton therapy

Purpose. To enhance an in-house graphic-processing-unit accelerated virtual particle (VP)-based Monte Carlo (MC) proton dose engine (VPMC) to model aperture blocks in both dose calculation and optimization for pencil beam scanning proton therapy (PBSPT)-based stereotactic radiosurgery (SRS).Methods and materials. A module to simulate VPs passing through patient-specific aperture blocks was developed and integrated in VPMC based on simulation results of realistic particles (primary protons and their secondaries). To validate the aperture block module, VPMC was first validated by an opensource MC code, MCsquare, in eight water phantom simulations with 3 cm thick brass apertures: four were with aperture openings of 1, 2, 3, and 4 cm without a range shifter, while the other four were with same aperture opening configurations with a range shifter of 45 mm water equivalent thickness. Then, VPMC was benchmarked with MCsquare and RayStation MC for 10 patients with small targets (average volume 8.4 c.c. with range of 0.4-43.3 c.c.). Finally, 3 typical patients were selected for robust optimization with aperture blocks using VPMC.Results. In the water phantoms, 3D gamma passing rate (2%/2 mm/10%) between VPMC and MCsquare was 99.71 ± 0.23%. In the patient geometries, 3D gamma passing rates (3%/2 mm/10%) between VPMC/MCsquare and RayStation MC were 97.79 ± 2.21%/97.78 ± 1.97%, respectively. Meanwhile, the calculation time was drastically decreased from 112.45 ± 114.08 s (MCsquare) to 8.20 ± 6.42 s (VPMC) with the same statistical uncertainties of ~0.5%. The robustly optimized plans met all the dose-volume-constraints (DVCs) for the targets and OARs per our institutional protocols. The mean calculation time for 13 influence matrices in robust optimization by VPMC was 41.6 s and the subsequent on-the-fly 'trial-and-error' optimization procedure took only 71.4 s on average for the selected three patients.Conclusion. VPMC has been successfully enhanced to model aperture blocks in dose calculation and optimization for the PBSPT-based SRS.

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.

Artificial general intelligence for radiation oncology

The emergence of artificial general intelligence (AGI) is transforming radiation oncology. As prominent vanguards of AGI, large language models (LLMs) such as GPT-4 and PaLM 2 can process extensive texts and large vision models (LVMs) such as the Segment Anything Model (SAM) can process extensive imaging data to enhance the efficiency and precision of radiation therapy. This paper explores full-spectrum applications of AGI across radiation oncology including initial consultation, simulation, treatment planning, treatment delivery, treatment verification, and patient follow-up. The fusion of vision data with LLMs also creates powerful multimodal models that elucidate nuanced clinical patterns. Together, AGI promises to catalyze a shift towards data-driven, personalized radiation therapy. However, these models should complement human expertise and care. This paper provides an overview of how AGI can transform radiation oncology to elevate the standard of patient care in radiation oncology, with the key insight being AGI's ability to exploit multimodal clinical data at scale.

Deep-learning based fast and accurate 3D CT deformable image registration in lung cancer

BACKGROUND

Deformable Image Registration (DIR) is an essential technique required in many applications of radiation oncology. However, conventional DIR approaches typically take several minutes to register one pair of 3D CT images and the resulting deformable vector fields (DVFs) are only specific to the pair of images used, making it less appealing for clinical application.

PURPOSE

A deep-learning-based DIR method using CT images is proposed for lung cancer patients to address the common drawbacks of the conventional DIR approaches and in turn can accelerate the speed of related applications, such as contour propagation, dose deformation, adaptive radiotherapy (ART), etc. METHODS: A deep neural network based on VoxelMorph was developed to generate DVFs using CT images collected from 114 lung cancer patients. Two models were trained with the weighted mean absolute error (wMAE) loss and structural similarity index matrix (SSIM) loss (optional) (i.e., the MAE model and the M+S model). In total, 192 pairs of initial CT (iCT) and verification CT (vCT) were included as a training dataset and the other independent 10 pairs of CTs were included as a testing dataset. The vCTs usually were taken 2 weeks after the iCTs. The synthetic CTs (sCTs) were generated by warping the vCTs according to the DVFs generated by the pre-trained model. The image quality of the synthetic CTs was evaluated by measuring the similarity between the iCTs and the sCTs generated by the proposed methods and the conventional DIR approaches, respectively. Per-voxel absolute CT-number-difference volume histogram (CDVH) and MAE were used as the evaluation metrics. The time to generate the sCTs was also recorded and compared quantitatively. Contours were propagated using the derived DVFs and evaluated with SSIM. Forward dose calculations were done on the sCTs and the corresponding iCTs. Dose volume histograms (DVHs) were generated based on dose distributions on both iCTs and sCTs generated by two models, respectively. The clinically relevant DVH indices were derived for comparison. The resulted dose distributions were also compared using 3D Gamma analysis with thresholds of 3 mm/3%/10% and 2 mm/2%/10%, respectively.

RESULTS

The two models (wMAE and M+S) achieved a speed of 263.7±163 / 265.8±190 ms and a MAE of 13.15±3.8 / 17.52±5.8 HU for the testing dataset, respectively. The average SSIM scores of 0.987±0.006 and 0.988±0.004 were achieved by the two proposed models, respectively. For both models, CDVH of a typical patient showed that less than 5% of the voxels had a per-voxel absolute CT-number-difference larger than 55 HU. The dose distribution calculated based on a typical sCT showed differences of ≤2cGy[RBE] for clinical target volume (CTV) D95 and D5 , within ±0.06% for total lung V5 , ≤1.5cGy[RBE] for heart and esophagus Dmean , and ≤6cGy[RBE] for cord Dmax compared to the dose distribution calculated based on the iCT. The good average 3D Gamma passing rates (> 96% for 3 mm/3%/10% and > 94% for 2 mm/2%/10%, respectively) were also observed.

CONCLUSION

A deep neural network-based DIR approach was proposed and has been shown to be reasonably accurate and efficient to register the initial CTs and verification CTs in lung cancer.

Selection of breathing phase number in 4D scanned proton treatment planning optimization for lung tumors

The standard four-dimensional (4D) treatment planning includes all breathing states in the optimization process, which is time-consuming. This work was aimed to optimize the number of intermediate phases needed for 4D proton treatment planning optimization to reduce the computational cost. Five 4D optimization strategies adopting different numbers of intermediate states and one three-dimensional (3D) optimization plan were studied for fifteen lung cancer patients treated with scanned protons, optimizing on all ten phases (4D_10), two extreme phases (4D_2), six phases during the exhalation stage (4D_6EX), six phases during the inhalation stage (4D_6IN), two extreme phases plus an intermediate state (4D_3) and average computed tomography image (3D), respectively. The 4D dose evaluation was conducted on all the ten phases, considering the interplay effect. The resulting doses accumulated on the reference phase were computed and compared. Compared to the 4D optimization plans, the 3D optimization plan performed inferiorly in target coverage, but superiorly in organ at risks (OARs) sparing. For the 4D optimization, all the five 4D plans showed similar performance in OARs protection. However, the 4D_6EX and 4D_6IN strategies out-performed the 4D_2 and 4D_3 plans in dose homogeneity. The computing times of the 4D_2, 4D_3, 4D_6EX and 4D_6IN approaches decreased to 32%, 41%, 66% and 67% of the 4D_10 method, respectively. Thus, our study suggested that the use of all phases during inhalation or exhalation stage might be a feasible approach substituting for the full phase strategy to reduce the calculation load while guaranteeing the plan quality for scanned proton therapy.

Modelling small block aperture in an in-house developed GPU-accelerated Monte Carlo-based dose engine for pencil beam scanning proton therapy

Purpose

To enhance an in-house graphic-processing-unit (GPU) accelerated virtual particle (VP)-based Monte Carlo (MC) proton dose engine (VPMC) to model aperture blocks in both dose calculation and optimization for pencil beam scanning proton therapy (PBSPT)-based stereotactic radiosurgery (SRS).

Methods and Materials

A module to simulate VPs passing through patient-specific aperture blocks was developed and integrated in VPMC based on simulation results of realistic particles (primary protons and their secondaries). To validate the aperture block module, VPMC was first validated by an opensource MC code, MCsquare, in eight water phantom simulations with 3cm thick brass apertures: four were with aperture openings of 1, 2, 3, and 4cm without a range shifter, while the other four were with same aperture opening configurations with a range shifter of 45mm water equivalent thickness. Then, VPMC was benchmarked with MCsquare and RayStation MC for 10 patients with small targets (average volume 8.4 cc with range of 0.4 - 43.3 cc). Finally, 3 typical patients were selected for robust optimization with aperture blocks using VPMC.

Results

In the water phantoms, 3D gamma passing rate (2%/2mm/10%) between VPMC and MCsquare was 99.71±0.23%. In the patient geometries, 3D gamma passing rates (3%/2mm/10%) between VPMC/MCsquare and RayStation MC were 97.79±2.21%/97.78±1.97%, respectively. Meanwhile, the calculation time was drastically decreased from 112.45±114.08 seconds (MCsquare) to 8.20±6.42 seconds (VPMC) with the same statistical uncertainties of ~0.5%. The robustly optimized plans met all the dose-volume-constraints (DVCs) for the targets and OARs per our institutional protocols. The mean calculation time for 13 influence matrices in robust optimization by VPMC was 41.6 seconds and the subsequent on-the-fly "trial-and-error" optimization procedure took only 71.4 seconds on average for the selected three patients.

Conclusion

VPMC has been successfully enhanced to model aperture blocks in dose calculation and optimization for the PBSPT-based SRS.

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.

Beam mask and sliding window-facilitated deep learning-based accurate and efficient dose prediction for pencil beam scanning proton therapy

PURPOSE

To develop a DL-based PBSPT dose prediction workflow with high accuracy and balanced complexity to support on-line adaptive proton therapy clinical decision and subsequent replanning.

METHODS

PBSPT plans of 103 prostate cancer patients and 83 lung cancer patients previously treated at our institution were included in the study, each with CTs, structure sets, and plan doses calculated by the in-house developed Monte-Carlo dose engine. For the ablation study, we designed three experiments corresponding to the following three methods: 1) Experiment 1, the conventional region of interest (ROI) method. 2) Experiment 2, the beam mask (generated by raytracing of proton beams) method to improve proton dose prediction. 3) Experiment 3, the sliding window method for the model to focus on local details to further improve proton dose prediction. A fully connected 3D-Unet was adopted as the backbone. Dose volume histogram (DVH) indices, 3D Gamma passing rates, and dice coefficients for the structures enclosed by the iso-dose lines between the predicted and the ground truth doses were used as the evaluation metrics. The calculation time for each proton dose prediction was recorded to evaluate the method's efficiency.

RESULTS

Compared to the conventional ROI method, the beam mask method improved the agreement of DVH indices for both targets and OARs and the sliding window method further improved the agreement of the DVH indices. For the 3D Gamma passing rates in the target, OARs, and BODY (outside target and OARs), the beam mask method can improve the passing rates in these regions and the sliding window method further improved them. A similar trend was also observed for the dice coefficients. In fact, this trend was especially remarkable for relatively low prescription isodose lines. The dose predictions for all the testing cases were completed within 0.25s.

Deep-Learning-based Fast and Accurate 3D CT Deformable Image Registration in Lung Cancer

PURPOSE

In some proton therapy facilities, patient alignment relies on two 2D orthogonal kV images, taken at fixed, oblique angles, as no 3D on-the-bed imaging is available. The visibility of the tumor in kV images is limited since the patient's 3D anatomy is projected onto a 2D plane, especially when the tumor is behind high-density structures such as bones. This can lead to large patient setup errors. A solution is to reconstruct the 3D CT image from the kV images obtained at the treatment isocenter in the treatment position.

METHODS

An asymmetric autoencoder-like network built with vision-transformer blocks was developed. The data was collected from 1 head and neck patient: 2 orthogonal kV images (1024x1024 voxels), 1 3D CT with padding (512x512x512) acquired from the in-room CT-on-rails before kVs were taken and 2 digitally-reconstructed-radiograph (DRR) images (512x512) based on the CT. We resampled kV images every 8 voxels and DRR and CT every 4 voxels, thus formed a dataset consisting of 262,144 samples, in which the images have a dimension of 128 for each direction. In training, both kV and DRR images were utilized, and the encoder was encouraged to learn the jointed feature map from both kV and DRR images. In testing, only independent kV images were used. The full-size synthetic CT (sCT) was achieved by concatenating the sCTs generated by the model according to their spatial information. The image quality of the synthetic CT (sCT) was evaluated using mean absolute error (MAE) and per-voxel-absolute-CT-number-difference volume histogram (CDVH).

RESULTS

The model achieved a speed of 2.1s and a MAE of <40HU. The CDVH showed that <5% of the voxels had a per-voxel-absolute-CT-number-difference larger than 185 HU.

CONCLUSION

A patient-specific vision-transformer-based network was developed and shown to be accurate and efficient to reconstruct 3D CT images from kV images.

Cardiopulmonary Toxicity Following Intensity-Modulated Proton Therapy (IMPT) Versus Intensity-Modulated Radiation Therapy (IMRT) for Stage III Non-Small Cell Lung Cancer

INTRODUCTION

Intensity-modulated proton therapy (IMPT) has the potential to reduce radiation dose to normal organs when compared to intensity-modulated radiation therapy (IMRT). We hypothesized that IMPT is associated with a reduced rate of cardiopulmonary toxicities in patients with Stage III NSCLC when compared with IMRT.

METHODS

We analyzed 163 consecutively treated patients with biopsy-proven, stage III NSCLC who received IMPT (n = 35, 21%) or IMRT (n = 128, 79%). Patient, tumor, and treatment characteristics were analyzed. Overall survival (OS), freedom-from distant metastasis (FFDM), freedom-from locoregional relapse (FFLR), and cardiopulmonary toxicities (CTCAE v5.0) were calculated using the Kaplan-Meier estimate. Univariate cox regressions were conducted for the final model.

RESULTS

Median follow-up of surviving patients was 25.5 (range, 4.6-58.1) months. Median RT dose was 60 (range, 45-72) Gy [RBE]. OS, FFDM, and FFLR were not different based on RT modality. IMPT provided significant dosimetric pulmonary and cardiac sparing when compared to IMRT. IMPT was associated with a reduced rate of grade more than or equal to 3 pneumonitis (HR 0.25, P = .04) and grade more than or equal to 3 cardiac events (HR 0.33, P = .08). Pre-treatment predicted diffusing capacity for carbon monoxide less than equal to 57% (HR 2.8, P = .04) and forced expiratory volume in the first second less than equal to 61% (HR 3.1, P = .03) were associated with an increased rate of grade more than or equal to 3 pneumonitis.

CONCLUSIONS

IMPT is associated with a reduced risk of clinically significant pneumonitis and cardiac events when compared with IMRT without compromising tumor control in stage III NSCLC. IMPT may provide a safer treatment option, particularly for high-risk patients with poor pretreatment pulmonary function.

Virtual particle Monte Carlo: A new concept to avoid simulating secondary particles in proton therapy dose calculation

BACKGROUND

In proton therapy dose calculation, Monte Carlo (MC) simulations are superior in accuracy but more time consuming, compared to analytical calculations. Graphic processing units (GPUs) are effective in accelerating MC simulations but may suffer thread divergence and racing condition in GPU threads that degrades the computing performance due to the generation of secondary particles during nuclear reactions.

PURPOSE

A novel concept of virtual particle (VP) MC (VPMC) is proposed to avoid simulating secondary particles in GPU-accelerated proton MC dose calculation and take full advantage of the computing power of GPU.

METHODS

Neutrons and gamma rays were ignored as escaping from the human body; doses of electrons, heavy ions, and nuclear fragments were locally deposited; the tracks of deuterons were converted into tracks of protons. These particles, together with primary and secondary protons, are considered to be the realistic particles. Histories of primary and secondary protons were replaced by histories of multiple VPs. Each VP corresponded to one proton (either primary or secondary). A continuous-slowing-down-approximation model, an ionization model, and a large angle scattering event model corresponding to nuclear interactions were developed for VPs by generating probability distribution functions (PDFs) based on simulation results of realistic particles using MCsquare. For efficient calculations, these PDFs were stored in the Compute Unified Device Architecture textures. VPMC was benchmarked with TOPAS and MCsquare in phantoms and with MCsquare in 13 representative patient geometries. Comparisons between the VPMC calculated dose and dose measured in water during patient-specific quality assurance (PSQA) of the selected 13 patients were also carried out. Gamma analysis was used to compare the doses derived from different methods and calculation efficiencies were also compared.

RESULTS

Integrated depth dose and lateral dose profiles in both homogeneous and inhomogeneous phantoms all matched well among VPMC, TOPAS, and MCsquare calculations. The 3D-3D gamma passing rates with a criterion of 2%/2 mm and a threshold of 10% was 98.49% between MCsquare and TOPAS and 98.31% between VPMC and TOPAS in homogeneous phantoms, and 99.18% between MCsquare and TOPAS and 98.49% between VPMC and TOPAS in inhomogeneous phantoms, respectively. In patient geometries, the 3D-3D gamma passing rates with 2%/2 mm/10% between dose distributions from VPMC and MCsquare were 98.56 ± 1.09% in patient geometries. The 2D-3D gamma analysis with 3%/2 mm/10% between the VPMC calculated dose distributions and the 2D measured planar dose distributions during PSQA was 98.91 ± 0.88%. VPMC calculation was highly efficient and took 2.84 ± 2.44 s to finish for the selected 13 patients running on four NVIDIA Ampere GPUs in patient geometries.

CONCLUSION

VPMC was found to achieve high accuracy and efficiency in proton therapy dose calculation.

Proton Beam Therapy for Esophageal Cancer

Simple Summary

Early-stage esophageal cancer is managed surgically, with the addition of radiotherapy for locally advanced disease. Current photon-based radiotherapy results in a high treatment-related complications, due to proximal organ involvement. The anatomic location of the esophagus raises challenges due to the anatomical changes associated with diaphragmatic motion, weight loss, tumor changes, and set-up variability. These propelled the interest in proton beam therapy (PBT), which theoretically offers a reduction in the radiation exposure to healthy neighboring tissues with improvements in the therapeutic ratio. In this review, we present the role of PBT for esophageal cancer, including treatment planning, early clinical comparisons with photon-based techniques, ongoing trials, current challenges, toxicities, and issues of equity and health services.

Abstract

Early-stage esophageal cancer is often primarily managed surgically, with the addition of radiotherapy for locally advanced disease. However, current photon-based radiotherapy regimens and surgery results in a high incidence of treatment-related cardiac and pulmonary complications due to the involvement of proximal organs at risk. In addition, the anatomic location of the esophagus raises challenges for radiotherapy due to the anatomical changes associated with diaphragmatic motion, weight loss, tumor changes, and set-up variability. These challenges propelled the interest in proton beam therapy (PBT), which theoretically offers a reduction in the radiation exposure to healthy neighboring tissues with improvements in the therapeutic ratio. Several dosimetric studies support the potential advantages of PBT for esophageal cancer treatment however, translation of these results to improved clinical outcomes remains unclear with limited clinical data, especially in large populations. Studies on the effect on quality of life are likewise lacking. Here, we review the existing and emerging role of PBT for esophageal cancer, including treatment planning, early clinical comparisons of PBT with photon-based techniques, recently concluded and ongoing clinical trials, challenges and toxicities, effects on quality of life, and global inequities in the treatment of esophageal cancer.

Exploratory study of seed spots analysis to characterize dose and linear energy transfer effect in adverse event initialization of pencil beam scanning proton therapy

BACKGROUND

Both dose and linear-energy-transfer (LET) could play a substantial role in adverse event (AE) initialization of cancer patients treated with pencil-beam-scanning proton therapy (PBS). However, not all the voxels within the AE regions are directly induced from the dose and LET effect. It is important to study the synergistic effect of dose and LET in AE initialization by only including a subset of voxels that are dosimetrically important.

PURPOSE

To perform exploratory investigation of the dose and LET effects upon AE initialization in PBS using seed spots analysis.

METHODS

113 head and neck (H&N) cancer patients receiving curative PBS were included. Among them, 20 patients experienced unanticipated CTCAEv4.0 grade≥3 AEs (AE group) and 93 patients did not (control group). Within the AE group, 13 AE patients were included in the seed spot analysis to derive the descriptive features of AE initialization and the remaining 7 mandible osteoradionecrosis patients and 93 control patients were used to derive the feature-based volume constraint of mandible osteoradionecrosis. The AE regions were contoured and the corresponding dose-LET volume histograms (DLVHs) of AE regions were generated for all patients in the AE group. We selected high LET voxels (the highest 5% of each dose bin) with a range of moderate to high dose (≥∼40 Gy[RBE]) as critical voxels. Critical voxels which were contiguous with each other were grouped into clusters. Each cluster was considered as a potential independent seed spot for AE initialization. Seed spots were displayed in a 2D dose-LET plane based on their mean dose and LET to derive the descriptive features of AE initialization. A volume constraint of mandible osteoradionecrosis was then established based on the extracted features using a receiver operating characteristic curve.

RESULTS

The product of dose and LET (xBD) was found to be a descriptive feature of seed spots leading to AE initialization in this preliminary study. The derived xBD volume constraint for mandible osteoradionecrosis showed good performance with an area-under-curve of 0.87 (sensitivity of 0.714 and specificity of 0.807 in the leave-one-out cross validation) for the very limited patient data included in this study.

CONCLUSION

Our exploratory study showed that both dose and LET were observed to be important in AE initializations. The derived xBD volume constraint could predict mandible osteoradionecrosis reasonably well in the very limited H&N cancer patient data treated with PBS included in this study. This article is protected by copyright. All rights reserved.

Reflections on beam configuration optimization for intensity-modulated proton therapy

Presumably, intensity-modulated proton radiotherapy (IMPT) is the most powerful form of proton radiotherapy. In the current state of the art, IMPT beam configurations (i.e. the number of beams and their directions) are, in general, chosen subjectively based on prior experience and practicality. Beam configuration optimization (BCO) for IMPT could, in theory, significantly enhance IMPT’s therapeutic potential. However, BCO is complex and highly computer resource-intensive. Some algorithms for BCO have been developed for intensity-modulated photon therapy (IMRT). They are rarely used clinically mainly because the large number of beams typically employed in IMRT renders BCO essentially unnecessary. Moreover, in the newer form of IMRT, volumetric modulated arc therapy, there are no individual static beams. BCO is of greater importance for IMPT because it typically employs a very small number of beams (2-4) and, when the number of beams is small, BCO is critical for improving plan quality. However, the unique properties and requirements of protons, particularly in IMPT, make BCO challenging. Protons are more sensitive than photons to anatomic changes, exhibit variable relative biological effectiveness along their paths, and, as recently discovered, may spare the immune system. Such factors must be considered in IMPT BCO, though doing so would make BCO more resource intensive and make it more challenging to extend BCO algorithms developed for IMRT to IMPT. A limited amount of research in IMPT BCO has been conducted; however, considerable additional work is needed for its further development to make it truly effective and computationally practical. This article aims to provide a review of existing BCO algorithms, most of which were developed for IMRT, and addresses important requirements specific to BCO for IMPT optimization that necessitate the modification of existing approaches or the development of new effective and efficient ones.

GPU-accelerated Monte Carlo-based online adaptive proton therapy - a feasibility study

PURPOSE

To develop an online Graphic-Processing-Unit (GPU)-accelerated Monte-Carlo-based adaptive radiation therapy (ART) workflow for pencil beam scanning (PBS) proton therapy to address inter-fraction anatomical changes in patients treated with PBS.

METHODS AND MATERIALS

A four-step workflow was developed using our in-house developed GPU-accelerated Monte-Carlo-based treatment planning system to implement online Monte-Carlo-based ART for PBS. The first step conducts diffeomorphic demon-based deformable image registration (DIR) to propagate contours on the initial planning CT (pCT) to the verification CT (vCT) to form a new structure set. The second step performs forward dose calculation of the initial plan on the vCT with the propagated contours after manual approval (possible modifications involved). The third step triggers a re-optimization of the plan depending on whether the verification dose meets the clinical requirements or not. A robust evaluation will be done for both the verification plan in the second step and the re-opotimized plan in the third step. The fourth step involves a two-stage (before and after delivery) patient specific quality assurance (PSQA) of the re-optimized plan. The before-delivery PSQA is to compare the plan dose to the dose calculated using an independent fast open-source Monte Carlo code, MCsquare. The after-delivery PSQA is to compare the plan dose to the dose re-calculated using the log file (spot MU, spot position, and spot energy) collected during the delivery. Jaccard index (JI), Dice similarity coefficients (DSCs), and Hausdorff distance (HD) were used to assess the quality of the propagated contours in the first step. A commercial plan evaluation software, ClearCheck™, was integrated into the workflow to carry out efficient plan evaluation. 3D Gamma analysis was used during the fourth step to ensure the accuracy of the plan dose from re-optimization. Three patients with three different disease sites were chosen to evaluate the feasibility of the online ART workflow for PBS.

RESULTS

For all three patients, the propagated contours were found to have good volume conformance [JI (lowest-highest: 0.833-0.983) and DSC (0.909-0.992)] but sub-optimal boundary coincidence [HD (2.37-20.76 mm)] for organs at risk (OARs). The verification dose evaluated by ClearCheck™ showed significant degradation of the target coverage due to the inter-fractional anatomical changes. Re-optimization on the vCT resulted in great improvement of the plan quality to a clinically acceptable level. 3D Gamma analyses of PSQA confirmed the accuracy of the plan dose before delivery (mean Gamma index = 98.74% with a threshold of 2%/2 mm/10%), and after delivery based on the log files (mean Gamma index = 99.05% with a threshold of 2%/2 mm/10%). The average time cost for the complete execution of the workflow was around 858 seconds, excluding the time for manual intervention.

CONCLUSION

The proposed online ART workflow for PBS was demonstrated to be efficient and effective by generating a re-optimized plan that significantly improved the plan quality. This article is protected by copyright. All rights reserved.

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

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

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.

Empirical Relative Biological Effectiveness (RBE) for Mandible Osteoradionecrosis (ORN) in Head and Neck Cancer Patients Treated With Pencil-Beam-Scanning Proton Therapy (PBSPT): A Retrospective, Case-Matched Cohort Study

Purpose

To retrospectively investigate empirical relative biological effectiveness (RBE) for mandible osteoradionecrosis (ORN) in head and neck (H&N) cancer patients treated with pencil-beam-scanning proton therapy (PBSPT).

Methods

We included 1,266 H&N cancer patients, of which, 931 patients were treated with volumetric-modulated arc therapy (VMAT) and 335 were treated with PBSPT. Among them, 26 VMAT and 9 PBSPT patients experienced mandible ORN (ORN group), while all others were included in the control group. To minimize the impact of the possible imbalance in clinical factors between VMAT and PBSPT patients in the dosimetric comparison between these two modalities and the resulting RBE quantification, we formed a 1:1 case-matched patient cohort (335 VMAT patients and 335 PBSPT patients including both the ORN and control groups) using the greedy nearest neighbor matching of propensity scores. Mandible dosimetric metrics were extracted from the case-matched patient cohort and statistically tested to evaluate the association with mandibular ORN to derive dose volume constraints (DVCs) for VMAT and PBSPT, respectively. We sought the equivalent constraint doses for VMAT so that the critical volumes of VMAT were equal to those of PBSPT at different physical doses. Empirical RBEs of PBSPT for ORN were obtained by calculating the ratio between the derived equivalent constraint doses and physical doses of PBSPT. Bootstrapping was further used to get the confidence intervals.

Results

Clinical variables of age, gender, tumor stage, prescription dose, chemotherapy, hypertension or diabetes, dental extraction, smoking history, or current smoker were not statistically related to the incidence of ORN in the overall patient cohort. Smoking history was found to be significantly associated with the ORN incidence in PBSPT patients only. V40Gy[RBE], V50Gy[RBE], and V60Gy[RBE] were statistically different (p<0.05) between the ORN and control group for VMAT and PBSPT. Empirical RBEs of 1.58(95%CI: 1.34-1.64), 1.34(95%CI: 1.23-1.40), and 1.24(95%: 1.15-1.26) were obtained for proton dose at 40 Gy[RBE=1.1], 50 Gy[RBE=1.1] and 60 Gy[RBE=1.1], respectively.

Conclusions

Our study suggested that RBEs were larger than 1.1 at moderate doses (between 40 and 60 Gy[RBE=1.1]) with high LET for mandible ORN. RBEs are underestimated in current clinical practice in PBSPT. The derived DVCs can be used for PBSPT plan evaluation and optimization to minimize the incidence rate of mandible ORN.

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.

Per-voxel constraints to minimize hot spots in linear energy transfer-guided robust optimization for base of skull head and neck cancer patients in IMPT

PURPOSE

Due to the employment of quadratic programming using soft constraints to implement dose volume constraints and the "trial-and-error" procedure needed to achieve a clinically acceptable plan, conventional dose volume constraints (upper limit) are not adequately effective in controlling small and isolated hot spots in the dose/linear energy transfer (LET) distribution. Such hot spots can lead to adverse events. In order to mitigate the risk of brain necrosis, one of the most clinically significant adverse events in patients receiving intensity-modulated proton therapy (IMPT) for base of skull (BOS) cancer, we propose per-voxel constraints to minimize hot spots in LET-guided robust optimization.

METHODS AND MATERIALS

Ten BOS cancer patients treated with IMPT were carefully selected by meeting one of the following conditions: (1) diagnosis of brain necrosis during follow-up; and (2) considered high risk for brain necrosis by not meeting dose constraints to the brain. An optimizing structure (BrainOPT) and an evaluating structure (BrainROI) that both contained the aforementioned hot dose regions in the brain were generated for optimization and evaluation, respectively. Two plans were generated for every patient: one using conventional dose-only robust optimization, the other using LET-guided robust optimization. The impact of LET was integrated into the optimization via a term of extra biological dose (xBD). A novel optimization tool of per-voxel constraints to control small and isolated hot spots in either the dose, LET, or combined (dose/LET) distribution was developed and used to minimize dose/LET hot spots of the selected structures. Indices from dose-volume histogram (DVH) and xBD dose-volume histogram (xBDVH) were used in the plan evaluation. A newly developed tool of the dose-LET-volume histogram (DLVH) was also adopted to illustrate the underlying mechanism. Wilcoxon signed-rank test was used for statistical comparison of the DVH and xBDVH indices between the conventional dose-only and the LET-guided robustly optimized plans.

RESULTS

Per-voxel constraints effectively and efficiently minimized dose hot spots in both dose-only and LET-guided robust optimization and LET hot spots in LET-guided robust optimization. Compared to the conventional dose-only robust optimization, the LET-guided robust optimization could generate plans with statistically lower xBD hot spots in BrainROI (VxBD,50 Gy[RBE], p = 0.009; VxBD,60 Gy[RBE], p = 0.025; xBD1cc, p = 0.017; xBD2cc, p = 0.022) with comparable dose coverage, dose hot spots in the target, and dose hot spots in BrainROI. DLVH analysis indicated that LET-guided robust optimization could either reduce LET at the same dose level or redistribute high LET from high dose regions to low dose regions.

CONCLUSION

Per-voxel constraint is a powerful tool to minimize dose/LET hot spots in IMPT. The LET-guided robustly optimized plans outperformed the conventional dose-only robustly optimized plans in terms of xBD hot spots control.

Executive Summary of Clinical and Technical Guidelines for Esophageal Cancer Proton Beam Therapy From the Particle Therapy Co-Operative Group Thoracic and Gastrointestinal Subcommittees

Radiation therapy (RT) is an integral component of potentially curative management of esophageal cancer (EC). However, RT can cause significant acute and late morbidity due to excess radiation exposure to nearby critical organs, especially the heart and lungs. Sparing these organs from both low and high radiation dose has been demonstrated to achieve clinically meaningful reductions in toxicity and may improve long-term survival. Accruing dosimetry and clinical evidence support the consideration of proton beam therapy (PBT) for the management of EC. There are critical treatment planning and delivery uncertainties that should be considered when treating EC with PBT, especially as there may be substantial motion-related interplay effects. The Particle Therapy Co-operative Group Thoracic and Gastrointestinal Subcommittees jointly developed guidelines regarding patient selection, treatment planning, clinical trials, and future directions of PBT for EC.

The role of motion management and position verification in lymphoma radiotherapy

In the last decades, the substantial technical progress in radiation oncology offered the opportunity for more accurate planning and delivery of treatment. At the same time, the evolution of systemic treatment and the advent of modern diagnostic tools allowed for more accurate staging and consequently a safe reduction of radiotherapy (RT) target volumes and RT doses in the treatment of lymphomas. As a result, incidental irradiation of organs at risk was reduced, with a consequent reduction of severe late toxicity in long-term lymphoma survivors. Nevertheless, these innovations warrant that professionals pay attention to concurrently ensure precise planning and dose delivery to the target volume and safe sparing of the organs at risk. In particular, target and organ motion should be carefully managed in order to prevent any compromise of treatment efficacy. Several aspects should be taken into account during the treatment pathway to minimise uncertainties and to apply a valuable motion management strategy, when needed. These include: reliable image registration between diagnostic and planning radiologic exams to facilitate the contouring process, image guidance to limit positioning uncertainties and to ensure the accuracy of dose delivery and management of lung motion through procedures of respiratory gating and breath control. In this review, we will cover the current clinical approaches to minimise these uncertainties in patients treated with modern RT techniques, with a particular focus on mediastinal lymphoma. In addition, since uncertainties have a different impact on the dose deposition of protons compared to conventional x-rays, the role of motion management and position verification in proton beam therapy (PBT) will be discussed in a separate section.

Normal Tissue Injury Induced by Photon and Proton Therapies: Gaps and Opportunities

Despite technological advances in radiation therapy (RT) and cancer treatment, patients still experience adverse effects. Proton therapy (PT) has emerged as a valuable RT modality that can improve treatment outcomes. Normal tissue injury is an important determinant of the outcome; therefore, for this review, we analyzed 2 databases: (1) clinical trials registered with ClinicalTrials.gov and (2) the literature on PT in PubMed, which shows a steady increase in the number of publications. Most studies in PT registered with ClinicalTrials.gov with results available are nonrandomized early phase studies with a relatively small number of patients enrolled. From the larger database of nonrandomized trials, we listed adverse events in specific organs/sites among patients with cancer who are treated with photons and protons to identify critical issues. The present data demonstrate dosimetric advantages of PT with favorable toxicity profiles and form the basis for comparative randomized prospective trials. A comparative analysis of 3 recently completed randomized trials for normal tissue toxicities suggests that for early stage non-small cell lung cancer, no meaningful comparison could be made between stereotactic body RT and stereotactic body PT due to low accrual (NCT01511081). In addition, for locally advanced non-small cell lung cancer, a comparison of intensity modulated RTwith passive scattering PT (now largely replaced by spot-scanned intensity modulated PT), PT did not provide any benefit in normal tissue toxicity or locoregional failure over photon therapy. Finally, for locally advanced esophageal cancer, proton beam therapy provided a lower total toxicity burden but did not improve progression-free survival and quality of life (NCT01512589). The purpose of this review is to inform the limitations of current trials looking at protons and photons, considering that advances in technology, physics, and biology are a continuum, and to advocate for future trials geared toward accurate precision RT that need to be viewed as an iterative process in a defined path toward delivering optimal radiation treatment. A foundational understanding of the radiobiologic differences between protons and photons in tumor and normal tissue responses is fundamental to, and necessary for, determining the suitability of a given type of biologically optimized RT to a patient or cohort.

Assessment of a diaphragm override strategy for robustly optimized proton therapy planning for esophageal cancer patients

Purpose

To ensure target coverage in the treatment of esophageal cancer, a density override to the region of diaphragm motion can be applied in the optimization process. Here, we evaluate the benefit of this approach during robust optimization for intensity modulated proton therapy (IMPT) planning.

Materials and methods

For 10 esophageal cancer patients, two robustly optimized IMPT plans were created either using (WDO) or not using (NDO) a diaphragm density override of 1.05 g/cm³ during plan optimization. The override was applied to the excursion of the diaphragm between exhale and inhale. Initial robustness evaluation was performed for plan acceptance (setup errors of 8 mm, range errors of ±3%), and subsequently, on all weekly repeated 4DCTs (setup errors of 2 mm, range errors of ±3%). Target coverage and hotspots were analyzed on the resulting voxel‐wise minimum (Vw_min) and voxel‐wise maximum (Vw_max) dose distributions.

Results

The nominal dose distributions were similar for both WDO and NDO plans. However, visual inspection of the Vw_max of the WDO plans showed hotspots behind the right diaphragm override region. For one patient, target coverage and hotspots improved by applying the diaphragm override. We found no differences in target coverage in the weekly evaluations between the two approaches.

Conclusion

The diaphragm override approach did not result in a clinical benefit in terms of planning and interfractional robustness. Therefore, we do not see added value in employing this approach as a default option during robust optimization for IMPT planning in esophageal cancer.

Technical Note: Multiple energy extraction techniques for synchrotron-based proton delivery systems may exacerbate motion interplay effects in lung cancer treatments

PURPOSE

The multiple energy extraction (MEE) delivery technique for synchrotron-based proton delivery systems reduces beam delivery time by decelerating the beam multiple times during one accelerator spill, but this might cause additional plan quality degradation due to intrafractional motion. We seek to determine whether MEE causes significantly different plan quality degradation compared to single energy extraction (SEE) for lung cancer treatments due to the interplay effect.

METHODS

Ten lung cancer patients treated with IMPT at our institution were nonrandomly sampled based on a representative range of tumor motion amplitudes, tumor volumes, and respiratory periods. Dose-volume histogram (DVH) indices from single-fraction SEE and MEE four-dimensional (4D) dynamic dose distributions were compared using the Wilcoxon signed-rank test. Distributions of monitor units (MU) to breathing phases were investigated for features associated with plan quality degradation. SEE and MEE DVH indices were compared in fractionated deliveries of the worst-case patient treatment scenario to evaluate the impact of fractionation.

RESULTS

There were no clinically significant differences in target mean dose, target dose conformity, or dose to organs-at-risk between SEE and MEE in single-fraction delivery. Three patients had significantly worse dose homogeneity with MEE compared to SEE (single-fraction mean D_5% -D_95% increased by up to 9.6% of prescription dose), and plots of MU distribution to breathing phases showed synchronization patterns with MEE but not SEE. However, after 30 fractions the patient in the worst-case scenario had clinically acceptable target dose homogeneity and coverage with MEE (mean D_5% -D_95% increased by 1% compared to SEE).

CONCLUSIONS

For some patients with breathing periods close to the mean spill duration, MEE resulted in significantly worse single-fraction target dose homogeneity compared to SEE due to the interplay effect. However, this was mitigated by fractionation, and target dose homogeneity and coverage were clinically acceptable after 30 fractions with MEE.

A Critical Review of LET-Based Intensity-Modulated Proton Therapy Plan Evaluation and Optimization for Head and Neck Cancer Management

In this review article, we review the 3 important aspects of linear-energy-transfer (LET) in intensity-modulated proton therapy (IMPT) for head and neck (H&N) cancer management. Accurate LET calculation methods are essential for LET-guided plan evaluation and optimization, which can be calculated either by analytical methods or by Monte Carlo (MC) simulations. Recently, some new 3D analytical approaches to calculate LET accurately and efficiently have been proposed. On the other hand, several fast MC codes have also been developed to speed up the MC simulation by simplifying nonessential physics models and/or using the graphics processor unit (GPU)–acceleration approach. Some concepts related to LET are also briefly summarized including (1) dose-weighted versus fluence-weighted LET; (2) restricted versus unrestricted LET; and (3) microdosimetry versus macrodosimetry. LET-guided plan evaluation has been clinically done in some proton centers. Recently, more and more studies using patient outcomes as the biological endpoint have shown a positive correlation between high LET and adverse events sites, indicating the importance of LET-guided plan evaluation in proton clinics. Various LET-guided plan optimization methods have been proposed to generate proton plans to achieve biologically optimized IMPT plans. Different optimization frameworks were used, including 2-step optimization, 1-step optimization, and worst-case robust optimization. They either indirectly or directly optimize the LET distribution in patients while trying to maintain the same dose distribution and plan robustness. It is important to consider the impact of uncertainties in LET-guided optimization (ie, LET-guided robust optimization) in IMPT, since IMPT is sensitive to uncertainties including both the dose and LET distributions. We believe that the advancement of the LET-guided plan evaluation and optimization will help us exploit the unique biological characteristics of proton beams to improve the therapeutic ratio of IMPT to treat H&N and other cancers.

Technical Note: 4D robust optimization in small spot intensity-modulated proton therapy (IMPT) for distal esophageal carcinoma

PURPOSE

To compare the dosimetric performances of small-spot three-dimensional (3D) and four-dimensional (4D) robustly optimized intensity-modulated proton (IMPT) plans in the presence of uncertainties and interplay effect simultaneously for distal esophageal carcinoma.

METHOD AND MATERIALS

Thirteen (13) patients were selected and re-planned with small-spot ( σ  ~ 2-6 mm) 3D and 4D robust optimization in IMPT, respectively. The internal clinical target volumes (CTV_high3d , CTV_low3d ) were used in 3D robust optimization. Different CTVs (CTV_high4d , CTV_low4d ) were generated by subtracting an inner margin of the motion amplitudes in three cardinal directions from the internal CTVs and used in 4D robust optimization. All patients were prescribed the same dose to CTVs (50 Gy[RBE] for CTV_high3d /CTV_high4d and 45 Gy[RBE] for CTV_low3d /CTV_low4d ). Dose-volume histogram (DVH) indices were calculated to assess plan quality. Comprehensive plan robustness evaluations that consisted of 300 perturbed scenarios (10 different motion patterns to consider irregular motion (sampled from a Gaussian distribution) and 30 different uncertainties scenarios (sampled from a 4D uniform distribution) combined), were performed to quantify robustness to uncertainties and interplay effect simultaneously. Wilcoxon signed-rank test was used for statistical analysis.

RESULTS

Compared to 3D robustly optimized plans, 4D robustly optimized plans had statistically improved target coverage and better sparing of lungs and heart (heart D_mean , P = 0.001; heart V_30Gy[RBE] , P = 0.001) in the nominal scenario. 4D robustly optimized plans had better robustness in target dose coverage (CTV_high4d V_100% , P = 0.002) and the protection of lungs and heart (heart D_mean , P = 0.001; heart V_30Gy[RBE] , P = 0.001) when uncertainties and interplay effect were considered simultaneously.

CONCLUSIONS

Even with small spots in IMPT, 4D robust optimization outperformed 3D robust optimization in terms of normal tissue protection and robustness to uncertainties and interplay effect simultaneously. Our findings support the use of 4D robust optimization to treat distal esophageal carcinoma with small spots in IMPT.

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.

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.

Feasibility of using megavoltage computed tomography to reduce proton range uncertainty: A simulation study

Purpose

To demonstrate that variation in chemical composition has a negligible effect on the mapping curve from relative electron density (RED) to proton stopping power ratio (SPR), and to establish the theoretical framework of using Megavoltage (MV) computed tomography (CT), instead of kilovoltage (kV) dual energy CT, to accurately estimate proton SPR.

Methods

A simulation study was performed to evaluate the effect of chemical composition variation on kVCT number and proton SPR. The simulation study involved both reference and simulated human tissues. The reference human tissues, together with their physical densities and chemical compositions, came from the ICRP publication 23. The simulated human tissues were created from the reference human tissues assuming that elemental percentage weight followed a Gaussian distribution. For all tissues, kVCT number and proton SPR were obtained through (a) theoretical calculation from tissue’s physical density and chemical composition which served as the ground truth, and (b) estimation from RED using the calibration curves established from the stoichiometric method. Deviations of the estimated values from the calculated values were quantified as errors in using RED to estimate kVCT number and proton SPR.

Results

Given a chemical composition variation of 5% (1σ) of the nominal percentage weights, the total estimation error of using RED to estimate kVCT number was 0.34%, 0.62%, and 0.77% and the total estimation error of using RED to estimate proton SPR was 0.30%, 0.22%, and 0.16% for fat tissues, non‐fat soft tissues and bone tissues, respectively.

Conclusion

Chemical composition had a negligible effect on the method of using RED to determine proton SPR. RED itself is sufficient to accurately determine proton SPR. MVCT number maintains a superb linear relationship with RED because it is highly dominated by Compton scattering. Therefore, MVCT has great potential in reducing the proton range uncertainty.

Development of robustness evaluation strategies for enabling statistically consistent reporting

Robustness evaluation of proton therapy treatment plans is essential for ensuring safe treatment delivery. However, available evaluation procedures feature a limited exploration of the actual robustness of the plan and generally do not provide confidence levels. This study compared established and more sophisticated robustness evaluation procedures, with quantified confidence levels. We have evaluated several robustness evaluation methods for 5 bilateral head-and-neck patients optimized considering spot scanning delivery and with a conventional CTV-to-PTV margin of 4 mm. Method (1) good practice scenario selection (GPSS) (e.g. +/- 4 mm setup error 3% range uncertainty); (2) statistically sound scenario selection (SSSS) either only on or both on and inside isoprobability hypersurface encompassing 90% of the possible errors; (3) statistically sound dosimetric selection (SSDS). In the last method, the 90% best plans were selected according to either target coverage quantified by D 95 (SSDS_D 95) or to an approximation of the final objective function (OF) used during treatment optimization (SSDS_OF). For all methods, we have considered systematic setup and systematic range errors. A mix of systematic and random setup errors were also simulated for SSDS, but keeping the same conventional margin of 4 mm. All robustness evaluations have been performed using the fast Monte Carlo dose engine MCsquare. Both SSSS strategies yielded on average very similar results. SSSS and GPSS yield comparable values for target coverage (within 0.5 Gy). The most noticeable differences were found for the CTV between GPSS, on the one hand, and SSDS_D 95 and SSDS_OF, on the other hand (average worst-case D 98 were 2.8 and 2.0 Gy larger than for GPSS, respectively). Simulating explicitly random errors in SSDS improved almost all DVH metrics. We have observed that the width of DVH-bands and the confidence levels depend on the method chosen to sample the scenarios. Statistically sound estimation of the robustness of the plan in the dosimetric space may provide an improved insight on the actual robustness of the plan for a given confidence level.

Proton beam radiotherapy for patients with early-stage and advanced lung cancer: a narrative review with contemporary clinical recommendations

Although lung cancer rates are decreasing nationally, lung cancer remains the leading cause of cancer related death. Despite advancements in treatment and technology, overall survival (OS) for lung cancer remains poor. Proton beam therapy (PBT) is an advanced radiation therapy (RT) modality for treatment of lung cancer with the potential to achieve dose escalation to tumor while sparing critical structures due to higher target conformality. In early and late-stage non-small cell lung cancer (NSCLC), dosimetric studies demonstrated reduced doses to organs at risk (OARs) such as the lung, spinal cord, and heart, and clinical studies report limited toxicities with PBT, including hypofractionated regimens. In limited-stage SCLC, studies showed that regimens chemo RT including PBT were well tolerated, which may help optimize clinical outcomes. Improved toxicity profiles may be beneficial in post-operative radiotherapy, for which initial dosimetric and clinical data are encouraging. Sparing of OARs may also increase the proportion of patients able to complete reirradiation for recurrent disease. However, there are various challenges of using PBT including a higher financial burden on healthcare and limited data supporting its cost-effectiveness. Further studies are needed to identify subgroups that benefit from PBT based on prognostic factors, and to evaluate PBT combined with immunotherapy, in order to elucidate the benefit that PBT may offer future lung cancer patients.

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.

Evaluation of continuous beam rescanning versus pulsed beam in pencil beam scanned proton therapy for lung tumours

The treatment of moving targets with pencil beam scanned proton therapy (PBS-PT) may rely on rescanning strategies to smooth out motion induced dosimetric disturbances. PBS-PT machines, such as Proteus®Plus (PPlus) and Proteus®One (POne), deliver a continuous or a pulsed beam, respectively. In PPlus, scaled (or no) rescanning can be applied, while POne implies intrinsic 'rescanning' due to its pulsed delivery. We investigated the efficacy of these PBS-PT delivery types for the treatment of lung tumours. In general, clinically acceptable plans were achieved, and PPlus and POne showed similar effectiveness.