Impact of proton PBS machine operating parameters on the effectiveness of layer rescanning for interplay effect mitigation in lung SBRT treatment.

BACKGROUND: Rescanning is a common technique used in proton pencil beam scanning to mitigate the interplay effect. Advances in machine operating parameters across different generations of particle therapy systems have led to improvements in beam delivery time (BDT). However, the potential impact of these improvements on the effectiveness of rescanning remains an underexplored area in the existing research. METHODS: We systematically investigated the impact of proton machine operating parameters on the effectiveness of layer rescanning in mitigating interplay effect during lung SBRT treatment, using the CIRS phantom. Focused on the Hitachi synchrotron particle therapy system, we explored machine operating parameters from our institution's current (2015) and upcoming systems (2025A and 2025B). Accumulated dynamic 4D dose were reconstructed to assess the interplay effect and layer rescanning effectiveness. RESULTS: Achieving target coverage and dose homogeneity within 2% deviation required 6, 6, and 20 times layer rescanning for the 2015, 2025A, and 2025B machine parameters, respectively. Beyond this point, further increasing the number of layer rescanning did not further improve the dose distribution. BDTs without rescanning were 50.4, 24.4, and 11.4 s for 2015, 2025A, and 2025B, respectively. However, after incorporating proper number of layer rescanning (six for 2015 and 2025A, 20 for 2025B), BDTs increased to 67.0, 39.6, and 42.3 s for 2015, 2025A, and 2025B machine parameters. Our data also demonstrated the potential problem of false negative and false positive if the randomness of the respiratory phase at which the beam is initiated is not considered in the evaluation of interplay effect. CONCLUSION: The effectiveness of layer rescanning for mitigating interplay effect is affected by machine operating parameters. Therefore, past clinical experiences may not be applicable to modern machines.

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.

A method for quantitative evaluations of scanning-proton dose distributions.

PURPOSE: Patient-Specific Quality Assurance (PSQA) measurement analysis depends on generating metrics representative of calculation and measurement agreement. Considering the heightened capability of discrete spot scanning protons to modulate individual dose voxels, a dose plane comparison approach that maintained all of the capabilities of the well-established γ test, but that also provided a more intuitive error parameterization, was desired. METHODS: Analysis was performed for 300 dose planes compared by searching all calculated points within a fixed radius around each measured pixel to determine the dose deviation. Dose plane agreement is reported as the dose difference minimum (DDM) within an empirically established search radius: ΔDmin(r). This per-pixel metric is aggregated into a histogram binned by dose deviation. Search-radius criteria were based on a weighted-beamlet 3σ spatial deviation from imaging isocenter. Equipment setup error was mitigated during analysis using tracked image registration, ensuring beamlet deviations to be the dominant source of spatial error. The percentage of comparison points with <3% dose difference determined pass rate. RESULTS: The mean beamlet radial deviation was 0.38mm from x-ray isocenter, with a standard deviation of 0.19mm, such that 99.9% of relevant pencil beams were within 1 mm of nominal. The dose-plane comparison data showed no change in passing rate between a 3%/1mm ΔDmin(r) analysis (97.6 +/- 3.6%) and a 3%/2mm γ test (97.7 +/- 3.2%). CONCLUSIONS: PSQA dose-comparison agreements corresponding to a search radius outside of machine performance limits are likely false positives. However, the elliptical shape of the γ test is too dose-restrictive with a spatial-error threshold set at 1 mm. This work introduces a cylindrical search shape, proposed herein as more relevant to plan quality, as part of the new DDM planar-dose comparison algorithm. DDM accepts all pixels within a given dose threshold inside the search radius, and carries forward plan-quality metrics in a straightforward manner for evaluation.

Small-spot intensity-modulated proton therapy and volumetric-modulated arc therapies for patients with locally advanced non-small-cell lung cancer: A dosimetric comparative study.

PURPOSE: To compare dosimetric performance of volumetric-modulated arc therapy (VMAT) and small-spot intensity-modulated proton therapy for stage III non-small-cell lung cancer (NSCLC). METHODS AND MATERIALS: A total of 24 NSCLC patients were retrospectively reviewed; 12 patients received intensity-modulated proton therapy (IMPT) and the remaining 12 received VMAT. Both plans were generated by delivering prescription doses to clinical target volumes (CTV) on averaged 4D-CTs. The dose-volume-histograms (DVH) band method was used to quantify plan robustness. Software was developed to evaluate interplay effects with randomized starting phases of each field per fraction. DVH indices were compared using Wilcoxon rank sum test. RESULTS: Compared with VMAT, IMPT delivered significantly lower cord Dmax , heart Dmean , and lung V5 Gy[ RBE ] with comparable CTV dose homogeneity, and protection of other OARs. In terms of plan robustness, the IMPT plans were statistically better than VMAT plans in heart Dmean , but were statistically worse in CTV dose coverage, cord Dmax , lung Dmean , and V5 Gy[ RBE ] . Other DVH indices were comparable. The IMPT plans still met the standard clinical requirements with interplay effects considered. CONCLUSIONS: Small-spot IMPT improves cord, heart, and lung sparing compared to VMAT and achieves clinically acceptable plan robustness at least for the patients included in this study with motion amplitude less than 11 mm. Our study supports the usage of IMPT to treat some lung cancer patients.

Use of a radial projection to reduce the statistical uncertainty of spot lateral profiles generated by Monte Carlo simulation.

Monte Carlo (MC) simulation has been used to generate commissioning data for the beam modeling of treatment planning system (TPS). We have developed a method called radial projection (RP) for postprocessing of MC-simulation-generated data. We used the RP method to reduce the statistical uncertainty of the lateral profile of proton pencil beams with axial symmetry. The RP method takes advantage of the axial symmetry of dose distribution to use the mean value of multiple independent scores as the representative score. Using the mean as the representative value rather than any individual score results in substantial reduction in statistical uncertainty. Herein, we present the concept and step-by-step implementation of the RP method, as well as show the advantage of the RP method over conventional measurement methods for generating lateral profile. Lateral profiles generated by both methods were compared to demonstrate the uncertainty reduction qualitatively, and standard error comparison was performed to demonstrate the reduction quantitatively. The comparisons showed that statistical uncertainty was reduced substantially by the RP method. Using the RP method to postprocess MC data, the corresponding MC simulation time was reduced by a factor of 10 without quality reduction in the generated result from the MC data. We concluded that the RP method is an effective technique to increase MC simulation efficiency for generating lateral profiles for axially symmetric pencil beams.

Exploratory study of the association of volumetric modulated arc therapy (VMAT) plan robustness with local failure in head and neck cancer.

This work is to show which is more relevant to cause local failures (LFs) due to patient setup uncertainty between the planning target volume (PTV) underdosage and the potential target underdosage subject to patient setup uncertainties in head and neck (H&N) cancer treated with volumetric-modulated arc therapy (VMAT). Thirteen LFs in 10 H&N patients treated by VMAT were analyzed. Measures have been taken to minimize the chances of insufficient target delineation for these patients and the patients were clinically determined to have LF based on the PET/CT scan results by an experienced radiologist and then reviewed by a second experienced radiation oncologist. Two methods were used to identify the possible locations of LF due to underdosage: (a) examining the standard VMAT plan, in which the underdosed volume in the nominal dose distribution (UVN) was generated by subtracting the volumes receiving the prescription doses from PTVs, and (b) plan robustness analysis, in which in addition to the nominal dose distribution, six perturbed dose distributions were created by translating the CT iso-center in three cardinal directions by the PTV margin. The coldest dose distribution was represented by the minimum of the seven doses in each voxel. The underdosed volume in the coldest dose distribution (UVC) was generated by subtracting the volumes receiving the prescription doses in the coldest dose distribution from the volumes receiving the prescription doses in the nominal dose distribution. UVN and UVC were subsequently examined for spatial association with the locations of LF. The association was tested using the binominal distribution and the Fisher's exact test of independence. We found that of 13 LFs, 11 were associated with UVCs (P = 0.011), while three were associated with UVNs (P = 0.99). We concluded that the possible target underdosage due to patient setup uncertainties appeared to be a more relevant factor associated with LF in VMAT for H&N cancer than the compromised PTV coverage at least for the patients included in this study.

Implementation of an improved dose-per-MU model for double-scattered proton beams to address interbeamline modulation width variability.

Because treatment planning systems (TPSs) generally do not provide monitor units (MUs) for double-scattered proton plans, models to predict MUs as a function of the range and the nominal modulation width requested of the beam delivery system, such as the one developed by the MGH group, have been proposed. For a given nominal modulation width, however, the measured modulation width depends on the accuracy of the vendor's calibration process and may differ from this nominal value, and also from one beamline to the next. Although such a difference can be replicated in our TPS, the output dependence on range and modulation width for each beam option or suboption has to be modeled separately for each beamline in order to achieve maximal 3% inaccuracy. As a consequence, the MGH output model may not be directly transferable. This work, therefore, serves to extend the model to more general clinic situations. In this paper, a parameterized linear-quadratic transformation is introduced to convert the nominal modulation width to the measured modulation width for each beam option or suboption on a per-beamline basis. Fit parameters are derived for each beamline from measurements of 60 reference beams spanning the minimum and maximum ranges, and modulation widths from 2 cm to full range per option or suboption. Using the modeled modulation width, we extract the MGH parameters for the output dependence on range and modulation width. Our method has been tested with 1784 patient-specific fields delivered across three different beamlines at our facility. For these fields, all measured outputs fall within 3%, and 64.4% fall within 1%, of our model. Using a parameterized linear-quadratic modulation width, MU calculation models can be established on a per-beamline basis for each double scattering beam option or suboption.

Investigation of dosimetric effect of beam current fluctuations in synchrotron-based proton PBS continuous scanning.

Objective.In proton pencil beam scanning (PBS) continuous delivery, the beam is continuously delivered without interruptions between spots. For synchrotron-based systems, the extracted beam current exhibits a spill structure, and recent publications on beam current measurements have demonstrated significant fluctuations around the nominal values. These fluctuations potentially lead to dose deviations from those calculated assuming a stable beam current. This study investigated the dosimetric implications of such beam current fluctuations during proton PBS continuous scanning.Approach.Using representative clinical proton PBS plans, we performed simulations to mimic a worst-case clinical delivery environment with beam current varies from 50% to 250% of the nominal values. The simulations used the beam delivery parameters optimized for the best beam delivery efficiency of the upcoming particle therapy system at Mayo Clinic Florida. We reconstructed the simulated delivered dose distributions and evaluated the dosimetric impact of beam current fluctuations.Main results.Despite significant beam current fluctuations resulting in deviations at each spot level, the overall dose distributions were nearly identical to those assuming a stable beam current. The 1 mm/1% Gamma passing rate was 100% for all plans. Less than 0.2% root mean square error was observed in the planning target volume dose-volume histogram. Minimal differences were observed in all dosimetric evaluation metrics.Significance.Our findings demonstrate that with our beam delivery system and clinical planning practice, while significant beam current fluctuations may result in large local move monitor unit deviations at each spot level, the overall impact on the dose distribution is minimal.

The Status and Challenges for Prostate Stereotactic Body Radiation Therapy Treatments in United States Proton Therapy Centers: An NRG Oncology Practice Survey.

Purpose: To report the current practice pattern of the proton stereotactic body radiation therapy (SBRT) for prostate treatments. Materials and Methods: A survey was designed to inquire about the practice of proton SBRT treatment for prostate cancer. The survey was distributed to all 30 proton therapy centers in the United States that participate in the National Clinical Trial Network in February, 2023. The survey focused on usage, patient selection criteria, prescriptions, target contours, dose constraints, treatment plan optimization and evaluation methods, patient-specific QA, and image-guided radiation therapy (IGRT) methods. Results: We received responses from 25 centers (83% participation). Only 8 respondent proton centers (32%) reported performing SBRT of the prostate. The remaining 17 centers cited 3 primary reasons for not offering this treatment: no clinical need, lack of volumetric imaging, and/or lack of clinical evidence. Only 1 center cited the reduction in overall reimbursement as a concern for not offering prostate SBRT. Several common practices among the 8 centers offering SBRT for the prostate were noted, such as using Hydrogel spacers, fiducial markers, and magnetic resonance imaging (MRI) for target delineation. Most proton centers (87.5%) utilized pencil beam scanning (PBS) delivery and completed Imaging and Radiation Oncology Core (IROC) phantom credentialing. Treatment planning typically used parallel opposed lateral beams, and consistent parameters for setup and range uncertainties were used for plan optimization and robustness evaluation. Measurements-based patient-specific QA, beam delivery every other day, fiducial contours for IGRT, and total doses of 35 to 40 GyRBE were consistent across all centers. However, there was no consensus on the risk levels for patient selection. Conclusion: Prostate SBRT is used in about 1/3 of proton centers in the US. There was a significant consistency in practices among proton centers treating with proton SBRT. It is possible that the adoption of proton SBRT may become more common if proton SBRT is more commonly offered in clinical trials.

The AAPM/ASTRO 2023 Core Physics Curriculum for Radiation Oncology Residents.

PURPOSE: The American Association of Physicists in Medicine Radiation Oncology Medical Physics Education Subcommittee (ROMPES) has updated the radiation oncology physics core curriculum for medical residents in the radiation oncology specialty. METHODS AND MATERIALS: Thirteen physicists from the United States and Canada involved in radiation oncology resident education were recruited to ROMPES. The group included doctorates and master's of physicists with a range of clinical or academic roles. Radiation oncology physician and resident representatives were also consulted in the development of this curriculum. In addition to modernizing the material to include new technology, the updated curriculum is consistent with the format of the American Board of Radiology Physics Study Guide Working Group to promote concordance between current resident educational guidelines and examination preparation guidelines. RESULTS: The revised core curriculum recommends 56 hours of didactic education like the 2015 curriculum but was restructured to provide resident education that facilitates best clinical practice and scientific advancement in radiation oncology. The reference list, glossary, and practical modules were reviewed and updated to include recent literature and clinical practice examples. CONCLUSIONS: ROMPES has updated the core physics curriculum for radiation oncology residents. In addition to providing a comprehensive curriculum to promote best practice for radiation oncology practitioners, the updated curriculum aligns with recommendations from the American Board of Radiology Physics Study Guide Working Group. New technology has been integrated into the curriculum. The updated curriculum provides a framework to appropriately cover the educational topics for radiation oncology residents in preparation for their subsequent career development.

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.

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.

The status and challenges for prostate SBRT treatments in United States proton therapy centers: An NRG Oncology practice survey.

A survey was designed to inquire about the practice of proton SBRT treatment for prostate cancer. The survey was distributed to all 30 proton therapy centers in the United States that participate in the National Clinical Trial Network in Feb. 2023. The survey focused on usage, patient selection criteria, prescriptions, target contours, dose constraints, treatment plan optimization and evaluation methods, patient-specific QA, and IGRT methods. Results: We received responses from 25 centers (83% participation). Only 8 respondent proton centers (32%) reported performing SBRT of the prostate. The remaining 17 centers cited three primary reasons for not offering this treatment: no clinical need, lack of volumetric imaging, and/or lack of clinical evidence. Only 1 center cited the reduction in overall reimbursement as a concern for not offering prostate SBRT. Several common practices among the 8 centers offering SBRT for the prostate were noted, such as using Hydrogel spacers, fiducial markers, and MRI for target delineation. Most proton centers (87.5%) utilized pencil beam scanning (PBS) delivery and completed Imaging and Radiation Oncology Core (IROC) phantom credentialing. Treatment planning typically used parallel opposed lateral beams, and consistent parameters for setup and range uncertainties were used for plan optimization and robustness evaluation. Measurements-based patient-specific QA, beam delivery every other day, fiducial contours for IGRT, and total doses of 35-40 GyRBE were consistent across all centers. However, there was no consensus on the risk levels for patient selection. Conclusion: Prostate SBRT is used in about 1/3 of proton centers in the US. There was a significant consistency in practices among proton centers treating with proton SBRT. It is possible that the adoption of proton SBRT may become more common if proton SBRT is more commonly offered in clinical trials.

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.

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.

Technical note: Delivery benefit and dosimetric implication of synchrotron-based proton pencil beam scanning using continuous scanning mode.

BACKGROUND: Discrete spot scanning (DSS) is the commonly used method for proton pencil beam scanning (PBS). There is lack of data on the dose-driven continuous scanning (DDCS). PURPOSE: To investigate delivery benefits and dosimetric implications of DDCS versus DSS for PBS systems. METHODS: The irradiation duty factor, beam delivery time (BDT), and dose deviation were simulated for eight treatment plans in prostate, head and neck, liver, and lung, with both conventional fractionation and hypofractionation schemes. DDCS results were compared with those of DSS. RESULTS: The DDCS irradiation duty factor (range, 11%-41%) was appreciably improved compared to DSS delivery (range, 4%-14%), within which, hypofractionation schemes had greater improvement than conventional fractionation. With decreasing stop ratio constraints, the DDCS BDT reduction was greater, but dose deviation also increased. With stop ratio constraints of 2, 1, 0.5, and 0, DDCS BDT reduction reached to 6%, 10%, 12%, and 15%, respectively, and dose deviation reached to 0.6%, 1.7%, 3.0%, and 5.2% root mean square error in PTV DVH, respectively. The 3%/2-mm gamma passing rate was greater than 99% with stop ratio constraints of 2 and 1, and greater than 95% with a stop ratio of 0.5. When the stop ratio constraint was removed, five of the eight treatment plans had a 3%/2-mm gamma passing rate greater than 95%, and the other three plans had a 3%/2-mm gamma passing rate between 90% and 95%. CONCLUSIONS: The irradiation duty factor was considerably improved with DDCS. Smaller stop ratio constraints led to shorter BDTs, but with the cost of larger dose deviations. Our finding suggested that a stop ratio of 1 constraint seems to yield acceptable DDCS dose deviation.

Investigation of scan path optimization in improving proton pencil beam scanning continuous delivery.

Objective. To investigate the impact of scan path optimization on the dose accuracy and beam delivery time (BDT) of proton pencil beam scanning in the dose-driven continuous scanning (DDCS).Approach. A diverse set of six clinical plans, representing various spot patterns and treatment sites, was used to evaluate the effectiveness of scan time optimization and scan length optimization. The DDCS dose discrepancy and BDT with optimized scan paths was compared to the default serpentine scan path.Main results. Both scan time optimization and scan path optimization were able to reduce the DDCS dose discrepancy compared to the default serpentine scan path. All plans, except for the layer repainting lung plan, achieved a 2%/2 mm gamma pass rate of over 99% and less than 1% PTV DVH root mean square error (RMSE) through scan path optimization. In the case of the layer repainting lung plan, when compared to the default serpentine scan path, the 2%/2 mm gamma pass rate showed improvements from 91.3% to 93.1% and 95.8%, while the PTV DVH RMSE decreased from 2.1% to 1.7% and 1.1% for scan time optimization and scan length optimization, respectively. Although scan time optimization resulted in shorter total scan times for all plans compared to the default scan path and scan length optimization tended to have longer total scan times. However, due to the short total scan times and their minimal contribution to the total BDT, the impact of scan path optimization on the total BDT was practically negligible.Significance. Both scan time optimization and scan length optimization proved to be effective in minimizing DDCS dose discrepancy. No definitive winner can be determined between these two optimization approaches. Both scan time and scan length optimization had minimal effect on the total BDT.

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.

Technical note: Comprehensive evaluations of gantry and couch rotation isocentricities for implementing proton stereotactic radiosurgery.

BACKGROUND: Mechanical accuracy should be verified before implementing a proton stereotactic radiosurgery (SRS) program. Linear accelerator (Linac)-based SRS systems often use electronic portal imaging devices (EPIDs) to verify beam isocentricity. Because proton therapy systems do not have EPID, beam isocentricity tests of proton SRS may still rely on films, which are not efficient. PURPOSE: To validate that our proton SRS system meets mechanical precision requirements and to present an efficient method to evaluate the couch and gantry's rotational isocentricity for our proton SRS system. METHODS: A dedicated applicator to hold brass aperture for proton SRS system was designed. The mechanical precision of the system was tested using a metal ball and film for 11 combinations of gantry and couch angles. A more efficient quality assurance (QA) procedure was developed, which used a scintillator device to replace the film. The couch rotational isocentricity tests were performed using orthogonal kV x-rays with the couch rotated isocentrically to five positions (0°, 315°, 270°, 225°, and 180°). At each couch position, the distance between the metal ball in kV images and the imaging isocenter was measured. The gantry isocentricity tests were performed using a cone-shaped scintillator and proton beams at five gantry angles (0°, 45°, 90°, 135°, and 180°), and the isocenter position and the distance of each beam path to the isocenter were obtained. Daily QA procedure was performed for 1 month to test the robustness and reproducibility of the procedure. RESULTS: The gantry and couch rotational isocentricity exhibited sub-mm precision, with most measurements within ±0.5 mm. The 1-month QA results showed that the procedure was robust and highly reproducible to within ±0.2 mm. The gantry isocentricity test using the cone-shaped scintillator was accurate and sensitive to variations of ±0.2 mm. The QA procedure was efficient enough to be completed within 30 min. The 1-month isocentricity position variations were within 0.5 mm, which demonstrating that the overall proton SRS system was stable and precise. CONCLUSION: The proton SRS Winston-Lutz QA procedure using a cone-shaped scintillator was efficient and robust. We were able to verify radiation delivery could be performed with sub-mm mechanical precision.

Study of linear energy transfer effect on rib fracture in breast patients receiving pencil-beamscanning proton therapy.

Purpose: To study the effect of proton linear energy transfer (LET) on rib fracture in breast cancer patients treated with pencil-beam scanning proton therapy (PBS) using a novel tool of dose-LET volume histogram (DLVH). Methods: From a prospective registry of patients treated with post-mastectomy proton therapy to the chest wall and regional lymph nodes for breast cancer between 2015 and 2020, we retrospectively identified rib fracture cases detected after completing treatment. Contemporaneously treated control patients that did not develop rib fracture were matched to patients 2:1 considering prescription dose, boost location, reconstruction status, laterality, chest wall thickness, and treatment year.The DLVH index, V(d, l), defined as volume(V) of the structure with at least dose(d) and LET(l), was calculated. DLVH plots between the fracture and control group were compared. Conditional logistic regression (CLR) model was used to establish the relation of V(d, l) and the observed fracture at each combination of d and l. The p-value derived from CLR model shows the statistical difference between fracture patients and the matched control group. Using the 2D p-value map derived from CLR model, the DLVH features associated with the patient outcomes were extracted. Results: Seven rib fracture patients were identified, and fourteen matched patients were selected for the control group. The median time from the completion of proton therapy to rib fracture diagnosis was 12 months (range 5 to 14 months). Two patients had grade 2 symptomatic rib fracture while the remaining 5 were grade 1 incidentally detected on imaging. The derived p-value map demonstrated larger V(0-36Gy[RBE], 4.0-5.0 keV/µm) in patients experiencing fracture (p<0.1). For example, the p value for V(30 Gy[RBE], 4.0 keV/um) was 0.069. Conclusions: In breast cancer patients receiving PBS, a larger volume of chest wall receiving moderate dose and high LET may result in increased risk of rib fracture.