Patient-specific quality assurance of dynamically-collimated proton therapy treatment plans.

BACKGROUND: The dynamic collimation system (DCS) provides energy layer-specific collimation for pencil beam scanning (PBS) proton therapy using two pairs of orthogonal nickel trimmer blades. While excellent measurement-to-calculation agreement has been demonstrated for simple cube-shaped DCS-trimmed dose distributions, no comparison of measurement and dose calculation has been made for patient-specific treatment plans. PURPOSE: To validate a patient-specific quality assurance (PSQA) process for DCS-trimmed PBS treatment plans and evaluate the agreement between measured and calculated dose distributions. METHODS: Three intracranial patient cases were considered. Standard uncollimated PBS and DCS-collimated treatment plans were generated for each patient using the Astroid treatment planning system (TPS). Plans were recalculated in a water phantom and delivered at the Miami Cancer Institute (MCI) using an Ion Beam Applications (IBA) dedicated nozzle system and prototype DCS. Planar dose measurements were acquired at two depths within low-gradient regions of the target volume using an IBA MatriXX ion chamber array. RESULTS: Measured and calculated dose distributions were compared using 2D gamma analysis with 3%/3 mm criteria and low dose threshold of 10% of the maximum dose. Median gamma pass rates across all plans and measurement depths were 99.0% (PBS) and 98.3% (DCS), with a minimum gamma pass rate of 88.5% (PBS) and 91.2% (DCS). CONCLUSIONS: The PSQA process has been validated and experimentally verified for DCS-collimated PBS. Dosimetric agreement between the measured and calculated doses was demonstrated to be similar for DCS-collimated PBS to that achievable with noncollimated PBS.

Simulation and experimental benchmarking of a proton pencil beam scanning nozzle model for development of MR-integrated proton therapy.

BACKGROUND: MR-integrated proton therapy is under development. It consists of the unique challenge of integrating a proton pencil beam scanning (PBS) beam line nozzle with an magnetic resonance imaging (MRI) scanner. The magnetic interaction between these two components is deemed high risk as the MR images can be degraded if there is cross-talk during beam delivery and image acquisition. PURPOSE: To create and benchmark a self-consistent proton PBS nozzle model for empowering the next stages of MR-integrated proton therapy development, namely exploring and de-risking complete integrated prototype system designs including magnetic shielding of the PBS nozzle. MATERIALS AND METHODS: Magnetic field (COMSOL Multiphysics ${\text{Multiphysics}}$ ) and radiation transport (Geant4) models of a proton PBS nozzle located at OncoRay (Dresden, Germany) were developed according to the manufacturers specifications. Geant4 simulations of the PBS process were performed by using magnetic field data generated by the COMSOL Multiphysics ${\text{Multiphysics}}$ simulations. In total 315 spots were simulated which consisted of a 40 × 30 cm 2 $40\times 30\,{\text{cm}}^{2}$ scan pattern with 5 cm spot spacings and for proton energies of 70, 100, 150, 200, and 220 MeV. Analysis of the simulated deflection at the beam isocenter plane was performed to determine the self-consistency of the model. The magnetic fringe field from a sub selection of 24 of the 315 spot simulations were directly compared with high precision magnetometer measurements. These focused on the maximum scanning setting of ± $\pm$  20 cm beam deflection as generated from the second scanning magnet in the PBS for a proton beam energy of 220 MeV. Locations along the beam line central axis (CAX) were measured at beam isocenter and downstream of 22, 47, 72, 97, and 122 cm. Horizontal off-axis positions were measured at 22 cm downstream of isocenter ( ± $\pm$  50, ± $\pm$  100, and ± $\pm$  150 cm from CAX). RESULTS: The proton PBS simulations had good spatial agreement to the theoretical values in all 315 spots examined at the beam line isocenter plane (0-2.9 mm differences or within 1.5 % of the local spot deflection amount). Careful analysis of the experimental measurements were able to isolate the changes in magnetic fields due solely to the scanning magnet contribution, and showed 1.9  ± $\pm$  1.2 µ T $\bf{\mu} {\text{T}}$ -9.4 ± $\pm$  1.2 µ T $\bf{\mu} {\text{T}}$ changes over the range of measurement locations. Direct comparison with the equivalent simulations matched within the measurement apparatus and setup uncertainty in all but one measurement point. CONCLUSIONS: For the first time a robust, accurate and self-consistent model of a proton PBS nozzle assembly has been created and successfully benchmarked for the purposes of advancing MR-integrated proton therapy research. The model will enable confidence in further simulation based work on fully integrated designs including MRI scanners and PBS nozzle magnetic shielding in order to de-risk and realize the full potential of MR-integrated proton therapy.

A novel treatment planning method via scissor beams for uniform-target-dose proton GRID with peak-valley-dose-ratio optimization.

BACKGROUND: Proton spatially fractionated RT (SFRT) can potentially synergize the unique advantages of using proton Bragg peak and SFRT peak-valley dose ratio (PVDR) to reduce the radiation-induced damage for normal tissues. Uniform-target-dose (UTD) proton GRID is a proton SFRT modality that can be clinically desirable and conveniently adopted since its UTD resembles target dose distribution in conventional proton RT (CONV). However, UTD proton GRID is not used clinically, which is likely due to the lack of an effective treatment planning method. PURPOSE: This work will develop a novel treatment planning method using scissor beams (SB) for UTD proton GRID, with the joint optimization of PVDR and dose objectives. METHODS: The SB method for spatial dose modulation in normal tissues with UTD has two steps: (1) a primary beam (PB) is halved with interleaved beamlets, to generate spatial dose modulation in normal tissues; (2) a complementary beam (CB) is added to fill in previously valley-dose positions in the target to generate UTD, while the CB is angled slightly from the PB, to maintain spatial dose modulation in normal tissues. A treatment planning method with PVDR optimization via the joint total variation and L1 (TVL1) regularization is developed to jointly optimize PVDR and dose objectives. The plan optimization solution is obtained using an iterative convex relaxation algorithm. RESULTS: The new methods SB and SB-TVL1 were validated in comparison with CONV. Compared to CONV of relatively homogeneous dose distribution, SB had modulated spatial dose pattern in normal tissues with UTD and comparable plan quality. Compared to SB, SB-TVL1 further maximized PVDR, with comparable dose-volume parameters. CONCLUSIONS: A novel SB method is proposed that can generate modulated spatial dose pattern in normal tissues to achieve UTD proton GRID. A treatment planning method with PVDR optimization capability via TVL1 regularization is developed that can jointly optimize PVDR and dose objectives for proton GRID.

Impact of intensity-modulated proton therapy in reducing radiation-induced lymphopenia in glioma patients.

Background: Current standard management in adult grades 2-4 gliomas includes maximal safe resection followed by adjuvant radiotherapy (RT) and chemotherapy. Radiation-induced lymphopenia (RIL) has been shown to possibly affect treatment outcomes adversely. Proton beam therapy (PBT) may reduce the volume of the normal brain receiving moderate radiation doses, and consequently RIL. Our aim was to evaluate the incidence and severity of RIL during proton beam therapy (PBT). Methods: We identified patients with grades 2-4 glioma treated with PBT at our center between January 2019 and December 2021. We evaluated the incidence and severity of RIL from weekly complete blood count (CBC) data collected during PBT and compared it to the patients who were treated with photon-based RT (XRT) at our center during the same time. Results: The incidence of any degree of lymphopenia (48% in PBT, vs. 81.2% in XRT, P value = .001) and severe lymphopenia (8% in PBT, vs. 24.6% in XRT, P value = .093) were both significantly lesser in patients who received PBT. Severe RIL in patients receiving PBT was seen in only CNS WHO Gr-4 tumors. Mean whole brain V20GyE and V25GyE inversely correlated to nadir ALC and were both significantly lower with PBT. Patients with lymphopenia during PBT showed a trend toward poorer progression-free survival (P = .053) compared to those with maintained lymphocyte counts. Conclusions: Proton therapy seems to have a superior sparing of normal brain to moderate dose radiation than photon-based RT and reduces the incidence of lymphopenia. Glioma patients with lymphopenia possibly have worse outcomes than the ones with maintained lymphocyte counts.

Two-dimensional oxygen-diffusion modelling for FLASH proton therapy with pencil beam scanning-Impact of diffusive tissue properties, dose, dose rate and scan patterns.

Objective.Oxygen depletion is generally believed to play an important role in the FLASH effect-a differential reduction of the radiosensitivity of healthy tissues, relative to that of the tumour under ultra-high dose-rate (UHDR) irradiation conditions. In proton therapy (PT) with pencil-beam scanning (PBS), the deposition of dose, and, hence, the degree of (radiolytic) oxygen depletion varies both spatially and temporally. Therefore, the resulting oxygen concentration and the healthy-tissue sparing effect through radiation-induced hypoxia varies both spatially and temporally as well.Approach.We propose and numerically solve a physical oxygen diffusion model to study these effects and their dependence on tissue parameters and the scan pattern in pencil-beam delivery. Since current clinical FLASH PT (FLASH-PT) is based on 250 MeV shoot-through (transmission) beams, for which dose and dose rate (DR) hardly vary with depth compared to the variation transverse to the beam axis, we focus on the two-dimensional case. We numerically integrate the model to obtain the oxygen concentration in each voxel as a function of time and extract voxel-based and spatially and temporarily integrated metrics for oxygen (FLASH) enhanced dose. Furthermore, we evaluate the impact on oxygen enhancement of standard pencil-beam delivery patterns and patterns that were optimised on dose-rate. Our model can contribute to the identification of tissue properties and pencil-beam delivery parameters that are critical for FLASH-PT and it may be used for the optimisation of FLASH-PT treatment plans and their delivery.Main results.(i) the diffusive properties of oxygen are critical for the steady state concentration and therefore the FLASH effect, even more so in two dimensions when compared to one dimension. (ii) The FLASH effect through oxygen depletion depends primarily on dose and less on other parameters. (iii) At a fixed fraction dose there is a slight dependence on DR. (iv) Scan patterns optimised on DR slightly increase the oxygen induced FLASH effect.Significance.To our best knowledge, this is the first study assessing the impact of scan-pattern optimization (SPO) in FLASH-PT with PBS on a biological FLASH model. While the observed impact of SPO is relatively small, a larger effect is expected for larger target volumes. A better understanding of the FLASH effect and the role of oxygen (depletion) therein is essential for the further development of FLASH-PT with PBS, and SPO.

Fetal dose assessment in a pregnant patient with brain tumor: A comparative study of proton PBS and 3DCRT/VMAT radiation therapy techniques.

PURPOSE: The treatment of brain tumors in pregnant patients poses challenges, as the out-of-field dose exposure to the fetus can potentially be harmful. A pregnant patient with prior radiation treatment was presented with a brain tumor at our clinic. This work reports on our pre-treatment study that compared fetal dose exposure between intensity-modulated proton therapy (IMPT) using pencil beam scanning (PBS) and conventional photon 3D conformal radiation therapy (3DCRT) and volumetric-modulated arc therapy (VMAT), and the subsequent pregnant patient's radiation treatment. MATERIALS AND METHODS: Pre-treatment measurements of clinical plans, 3DCRT, VMAT, and IMPT, were conducted on a phantom. Measurements were performed using a device capable of neutron detections, closely following AAPM guidelines, TG158. For photon measurements, fetus shielding was utilized. On patient treatment days, which was determined to be proton treatment, shielding was used only during daily imaging for patient setup. Additionally, an in vivo measurement was conducted on the patient. RESULTS: Measurements showed that IMPT delivered the lowest fetal dose, considering both photon and neutron out-of-field doses to the fetus, even when shielding was implemented for photon measurements. Additionally, the proton plans demonstrated superior treatment for the mother, a reirradiation case. CONCLUSION: The patient was treated with proton therapy, and the baby was subsequently delivered at full term with no complications. This case study supports previous clinical findings and advocates for the expanded use of proton therapy in this patient population.

Validation and reproducibility of in vivo dosimetry for pencil beam scanned FLASH proton treatment in mice.

PURPOSE: To investigate quality assurance (QA) techniques for in vivo dosimetry and establish its routine uses for proton FLASH small animal experiments with a saturated monitor chamber. METHODS AND MATERIALS: 227 mice were irradiated at FLASH or conventional (CONV) dose rates with a 250 MeV FLASH-capable proton beamline using pencil beam scanning to characterize the proton FLASH effect on abdominal irradiation and examining various endpoints. A 2D strip ionization chamber array (SICA) detector was positioned upstream of collimation and used for in vivo dose monitoring during irradiation. Before each irradiation series, SICA signal was correlated with the isocenter dose at each delivered dose rate. Dose, dose rate, and 2D dose distribution for each mouse were monitored with the SICA detector. RESULTS: Calibration curves between the upstream SICA detector signal and the delivered dose at isocenter had good linearity with minimal R2 values of 0.991 (FLASH) and 0.985 (CONV), and slopes were consistent for each modality. After reassigning mice, standard deviations were less than 1.85 % (FLASH) and 0.83 % (CONV) for all dose levels, with no individual subject dose falling outside a ± 3.6 % range of the designated dose. FLASH fields had a field-averaged dose rate of 79.0 ± 0.8 Gy/s and mean local average dose rate of 160.6 ± 3.0 Gy/s. In vivo dosimetry allowed for the accurate detection of variation between the delivered and the planned dose. CONCLUSION: In vivo dosimetry benefits FLASH experiments through enabling real-time dose and dose rate monitoring allowing mouse cohort regrouping when beam fluctuation causes delivered dose to vary from planned dose.

Proton beam spot size and position measurements using a multi-strip ionization chamber.

PURPOSE: To develop and characterize a large-area multi-strip ionization chamber (MSIC) for efficient measurement of proton beam spot size and position at a synchrotron-based proton therapy facility. METHODS AND MATERIALS: A 420 mm x 320 mm MSIC was designed with 240 vertical strips and 180 horizontal strips at 1.75 mm pitch. The MSIC was characterized by irradiating a grid of proton spots across 17 energies from 73.5 MeV to 235 MeV and comparing to simultaneous measurements made with a reference Gafchromic EBT3 film. Beam profiles, spot sizes, and positions were analyzed. Short term measurement stability and sensitivity were evaluated. RESULTS: Excellent agreement was demonstrated between the MSIC and EBT3 film for both spot size and position measurements. Spot sizes agreed within ± 0.18 mm for all energies tested. Measured beam spot positions agreed within ± 0.17 mm. The detector showed good short term measurement stability and low noise performance. CONCLUSION: The large-area MSIC enables efficient and accurate proton beam spot characterization across the clinical energy range. The results indicate the MSIC is suitable for pencil beam scanning proton therapy commissioning and quality assurance applications requiring fast spot size and position quantification.

Beam Position Projection Algorithms in Proton Pencil Beam Scanning.

Beam position uncertainties along the beam trajectory arise from the accelerator, beamline, and scanning magnets (SMs). They can be monitored in real time, e.g., through strip ionization chambers (ICs), and treatments can be paused if needed. Delivery is more reliable and accurate if the beam position is projected from monitored nozzle parameters to the isocenter, allowing for accurate online corrections to be performed. Beam position projection algorithms are also used in post-delivery log file analyses. In this paper, we investigate the four potential algorithms that can be applied to all pencil beam scanning (PBS) nozzles. For some combinations of nozzle configurations and algorithms, however, the projection uses beam properties determined offline (e.g., through beam tuning or technical commissioning). The best algorithm minimizes either the total uncertainty (i.e., offline and online) or the total offline uncertainty in the projection. Four beam position algorithms are analyzed (A1-A4). Two nozzle lengths are used as examples: a large nozzle (1.5 m length) and a small nozzle (0.4 m length). Three nozzle configurations are considered: IC after SM, IC before SM, and ICs on both sides. Default uncertainties are selected for ion chamber measurements, nozzle entrance beam position and angle, and scanning magnet angle. The results for other uncertainties can be determined by scaling these results or repeating the error propagation. We show the propagation of errors from two locations and the SM angle to the isocenter for all the algorithms. The best choice of algorithm depends on the nozzle length and is A1 and A3 for the large and small nozzles, respectively. If the total offline uncertainty is to be minimized (a better choice if the offline uncertainty is not stable), the best choice of algorithm changes to A1 for the small nozzle for some hardware configurations. Reducing the nozzle length can help to reduce the gantry size and make proton therapy more accessible. This work is important for designing smaller nozzles and, consequently, smaller gantries. This work is also important for log file analyses.

Optimal minimum MU for intensity-modulated proton therapy with pencil-beam scanning proton beams.

PURPOSE: A higher minimum monitor unit (minMU) for pencil-beam scanning proton beams in intensity-modulated proton therapy is preferred for more efficient delivery. However, plan quality may be compromised when the minMU is too large. This study aimed to identify the optimal minMU (OminMU) to improve plan delivery efficiency while maintaining high plan quality. METHODS: We utilized clinical plans including six anatomic sites (brain, head and neck, breast, lung, abdomen, and prostate) from 23 patients previously treated with the Varian ProBeam system. The minMU of each plan was increased from the current clinical minMU of 1.1 to 3-24 MU depending on the daily prescribed dose (DPD). The dosimetric parameters of the plans were evaluated for consistency against a 1.1-minMU plan for target coverage as well as organs-at-risk dose sparing. DPD/minMU was defined as the ratio of DPD to minMU (cGy/MU) to find the OminMU by ensuring that dosimetric parameters did not differ by >1% compared to those of the 1.1-minMU plan. RESULTS: All plans up to 5 minMU showed no significant dose differences compared to the 1.1-minMU plan. Plan qualities remained acceptable when DPD/minMU ≥35 cGy/MU. This suggests that the 35 cGy/MU criterion can be used as the OminMU, which implies that 5 MU is the OminMU for a conventional fraction dose of 180 cGy. Treatment times were decreased by an average of 32% (max 56%, min 7%) and by an average of 1.6 min when the minMU was increased from 1.1 to OminMU. CONCLUSION: A clinical guideline for OminMU has been established. The minMU can be increased by 1 MU for every 35 cGy of DPD without compromising plan quality for most cases analyzed in this study. Significant treatment time reduction of up to 56% was observed when the suggested OminMU is used.

Proton therapy re-irradiation outcomes and genomic landscape of patients with recurrent head and neck cancer.

INTRODUCTION: Re-irradiation (re-RT) for recurrent head and neck cancer (rHNC) is challenging. We describe clinical outcomes and toxicity of proton therapy (PT) for recurrent HNC, and report genomic alterations associated with patterns of failure. MATERIALS & METHODS: We performed a retrospective analysis of rHNC patients treated with PT. Outcomes were estimated using the Kaplan-Meier method. Univariate (UVA) and multivariate analyses (MVA) were performed to assess multiple patient factors. Next-generation sequencing and genomic analyses were performed on available samples. RESULTS: Eighty-nine patients treated with PBS-PT for rHNC with a median follow-up of 12 mo (0-71 mo) were included. The 1- and 2-y local control (LC) rates were 80.8 % (95 % CI: 70.8-90.8) and 66.2 % (95 % CI: 50.7-81.7), and 1- and 2-y distant metastasis-free survival (DMFS) were 41.0 % (95 % CI: 30.0-52.0) and 26.3 % (95 % CI: 15.7-36.9). The median overall survival (OS) was 13 mo (95 % CI: 9.3-16.7). On UVA and MVA, smaller gross tumor volume (GTV) was associated with improved OS (HR 1.002, P = 0.004), DMFS (HR 1.002, P = 0.004), and PFS (HR 1.002, P = 0.014). There were 35 late Gr3 + toxicity events (30.3 %). Patients with higher candidate gene-specific mutation burden (genes with [OR] > 2, P < 0.05) had inferior PFS. TP53, NOTCH4, and ARID1B mutations were associated with inferior DMFS (OR > 2, P < 0.05). CONCLUSIONS: PBS-PT is effective at achieving LC for rHNC with favorable toxicity. Distant metastases are common, and associated with TP53, NOTCH4, and ARID1B mutations. Inclusion of genomic alterations in the clinical decision process may be warranted.

Cardinality-constrained plan-quality and delivery-time optimization method for proton therapy.

BACKGROUND: While minimizing plan delivery time is beneficial for proton therapy in terms of motion management, patient comfort, and treatment throughput, it often poses a tradeoff with optimizing plan quality. A key component of plan delivery time is the energy switching time, which is approximately proportional to the number of energy layers, that is, the cardinality. PURPOSE: This work aims to develop a novel optimization method that can efficiently compute the pareto surface between plan quality and energy layer cardinality, for the planner to navigate through this quality-and-efficiency tradeoff and select the appropriate plan of a balanced tradeoff. METHODS: A new IMPT method CARD is proposed that (1) explicitly incorporates the minimization of energy layer cardinality as an optimization objective, and (2) automatically generates a set of plans sequentially with a descending order in number of energy layers. The energy layer cardinality is penalized through the l1,0-norm regularization with an upper bound, and the upper bound is monotonically decreased to compute a series of treatment plans with gradually decreased energy layer cardinality on the quality-and-efficiency pareto surface. For any given treatment plan, the plan optimality is enforced using dose-volume planning objectives and the plan deliverability is imposed through minimum-monitor-unit (MMU) constraints, with optimization solution algorithm based on iterative convex relaxation. RESULTS: The new method CARD was validated in comparison with the benchmark plan of all energy layers (P0), and a state-of-the-art method called MMSEL, using prostate, head-and-neck (HN), lung, pancreas, liver and brain cases. While labor-intensive and time-consuming manual parameter tuning was needed for MMSEL to generate plans of predefined energy layer cardinality, CARD automatically and efficiently computed all plans with sequentially decreasing predefined energy layer cardinality all at once. With the acceptable plan quality (i.e., no more than 110% of total optimization objective value from P0), CARD achieved the reduction of number of energy layers to 52% (from 77 to 40), 48% (from 135 to 65), 59% (from 85 to 50), 67% (from 52 to 35), 80% (from 50 to 40), and 30% (from 66 to 20), for prostate, HN, lung, pancreas, liver, and brain cases, respectively, compared to P0, with overall better plan quality than MMSEL. Moreover, due to the nonconvexity of the MMU constraint, CARD provided the similar or even smaller optimization objective than P0, at the same time with fewer number of energy layers, that is, 55 versus 77, 85 versus 135, 45 versus 52, and 25 versus 66 for prostate, HN, pancreas, and brain cases, respectively. CONCLUSIONS: We have developed a novel optimization algorithm CARD that can efficiently and automatically compute a series of treatment plans of any given energy layer sequentially, which allows the planner to navigate through the plan-quality and energy-layer-cardinality tradeoff and select the appropriate plan of a balanced tradeoff.

Comparing different boost concepts and beam configurations for proton therapy of pancreatic cancer.

Background and Purpose: Interfractional geometrical and anatomical variations impact the accuracy of proton therapy for pancreatic cancer. This study investigated field-in-field (FIF) and simultaneous integrated boost (SIB) concepts for scanned proton therapy treatment with different beam configurations. Materials and Methods: Robustly optimized treatment plans for fifteen patients were generated using FIF and SIB techniques with two, three, and four beams. The prescribed dose in 20 fractions was 60 Gy(RBE) for the internal gross tumor volume (IGTV) and 46 Gy(RBE) for the internal clinical target volume. Verification computed tomography (vCT) scans was performed on treatment days 1, 7, and 16. Initial treatment plans were recalculated on the rigidly registered vCTs. V100% and D95% for targets and D2cm3 for the stomach and duodenum were evaluated. Robustness evaluations (range uncertainty of 3.5 %) were performed to evaluate the stomach and duodenum dose-volume parameters. Results: For all techniques, IGTV V100% and D95% decreased significantly when recalculating the dose on vCTs (p < 0.001). The median IGTV V100% and D95% over all vCTs ranged from 74.2 % to 90.2 % and 58.8 Gy(RBE) to 59.4 Gy(RBE), respectively. The FIF with two and three beams, and SIB with two beams maintained the highest IGTV V100% and D95%. In robustness evaluations, the ΔD2cm3 of stomach was highest in two beams plans, while the ΔD2cm3 of duodenum was highest in four beams plans, for both concepts. Conclusion: Target coverage decreased when recalculating on CTs at different time for both concepts. The FIF with three beams maintained the highest IGTV coverage while sparing normal organs the most.

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.

Combined optimization of spot positions and weights for better FLASH proton therapy.

Objective.In Intensity Modulated Proton Therapy (IMPT), the weights of individual pencil-beams or spots are optimized to fulfil dosimetric constraints. Theses spots are usually located on a regular lattice and their positions are fixed during optimization. In many cases, the range of spot weights may however be limited, leading sometimes to sub-optimal plan quality. An emblematic use case is the delivery of a plan at ultra-high dose rate (FLASH-RT), for which the spot weights are typically constrained to high values.Approach. To improve further the quality of IMPT FLASH plans, we propose here a novel algorithm to optimize both the spot weights and positions directly based on the objectives defined by the treatment planner.Main results. For all cases considered, optimizing the spot positions lead to an enhanced dosimetric score, while maintaining a high dose rate.Significance. Overall, this approach resulted in a substantial plan quality improvement compared to optimizing only the spot weights, and in a similar execution time.

Esophageal Cancer Outcomes After Definitive Chemotherapy With Intensity Modulated Proton Therapy.

Purpose: The effectiveness of intensity-modulated proton therapy (IMPT) for esophageal cancer treated with definitive concurrent chemoradiation therapy remains inadequately explored. We investigated long-term outcomes and toxicity experienced by patients who received IMPT as part of definitive esophageal cancer treatment. Patients and Methods: We retrospectively identified and analyzed 34 patients with locally advanced esophageal cancer who received IMPT with concurrent chemotherapy as a definitive treatment regimen at The University of Texas MD Anderson Cancer Center from 2011 to 2021. The median IMPT dose was 50.4 GyRBE in 28 fractions; concurrent chemotherapy consisted of fluorouracil and/or taxane and/or platinum. Survival outcomes were determined by the Kaplan-Meier method, and toxicity was scored according to the Common Terminology Criteria for Adverse Events version 4.0. Results: The median age of all patients was 71.5 years. Most patients had stage III (cT3 cM0) adenocarcinoma of the lower esophagus. At a median follow-up time of 39 months, the 5-year overall survival rate was 41.1%; progression-free survival, 34.6%; local regional recurrence-free survival, 78.1%; and distant metastasis-free survival, 65.0%. Common acute chemoradiation therapy-related toxicities included hematologic toxicity, esophagitis (and late-onset), fatigue, weight loss, and nausea (and late-onset); grade 3 toxicity rates were 26.0% for hematologic, 18.0% for esophagitis and 9.0% for nausea. No patient had grade ≥3 wt loss or radiation pneumonitis, and no patients had pulmonary fibrosis or esophageal fistula. No grade ≥4 events were observed except for hematologic toxicity (lymphopenia) in 2 patients. Conclusion: Long-term survival and toxicity were excellent after IMPT for locally advanced esophageal cancer treated definitively with concurrent chemoradiation therapy. When available, IMPT should be offered to such patients to minimize treatment-related cardiopulmonary toxicity without sacrificing outcomes.

Repurposing DailyQA3 for an efficient and spot position sensitive daily quality assurance tool for proton therapy.

INTRODUCTION: Daily quality assurance is an integral part of a radiotherapy workflow to ensure the dose is delivered safely and accurately to the patient. It is performed before the first treatment of the day and needs to be time and cost efficient for a multiple gantries proton center. In this study, we introduced an efficient method to perform QA for output constancy, range verification, spot positioning accuracy and imaging and proton beam isocenter coincidence with DailyQA3. METHODS: A stepped acrylic block of specific dimensions is fabricated and placed on top of the DailyQA3 device. Treatment plans comprising of two different spread-out Bragg peaks and five individual spots of 1.0 MU each are designed to be delivered to the device. A mathematical framework to measure the 2D distance between the detectors and individual spot is introduced and play an important role in realizing the spot positioning and centering QA. Lastly, a 5 months trends of the QA for two gantries are presented. RESULTS: The outputs are monitored by two ion chambers in the DailyQA3 and a tolerance of ± 3 % $ \pm 3\% $ are used. The range of the SOBPs are monitored by the ratio of ion chamber signals and a tolerance of ± 1 mm $ \pm 1\ {\mathrm{mm}}$ is used. Four diodes at ± 10 cm $ \pm 10\ {\mathrm{cm}}$ from the central ion chambers are used for spot positioning QA, while the central ion chamber is used for imaging and proton beam isocenter coincidence QA. Using the framework, we determined the absolute signal threshold corresponding to the offset tolerance between the individual proton spot and the detector. A 1.5 mm $1.5\ {\mathrm{mm}}$ tolerances are used for both the positioning and centering QA. No violation of the tolerances is observed in the 5 months trends for both gantries. CONCLUSION: With the proposed approach, we can perform four QA items in the TG224 within 10 min.

Development and characterization of the first proton minibeam system for single-gantry proton facility.

BACKGROUND: Minibeam represents a preclinical spatially fractionated radiotherapy modality with great translational potential. The advantage lies in its high therapeutic index (compared to GRID and LATTICE) and ability to treat at greater depth (compared to microbeam). Proton minibeam radiotherapy (pMBRT) is a synergy of proton and minibeam. While the single-gantry proton facility has gained popularity due to its affordability and compact design, it often has limited beam time available for research purposes. Conversely, given the current requirement of pMBRT on specific minibeam hardware collimators, necessitates a reproducible and fast setup to minimize pMBRT treatment time and streamline the switching time between pMBRT and conventional treatment for clinically translation. PURPOSE: The contribution of this work is the development and characterization of the first pMBRT system tailored for single-gantry proton facility. The system allows for efficient and reproducible plug-and-play setup, achievable within minutes. METHODS: The single room pMBRT system is constructed based on IBA ProteusONE proton machine. The end of nozzle is attached with beam modifying accessories though an accessory drawer. A small snout is attached to the accessory drawer and used to hold apertures and range shifters. The minibeam aperture consists of two components: a fitting ring and an aperture body. Three minibeam apertures were manufactured. The first-generation apertures underwent qualitatively analysis with film, and the second generation aperture underwent more comprehensive quantitative measurement. The reproducibility of the setup is accessed, and the film measurements are performed to characterize the pMBRT system in cross validation with Monte Carlo (MC) simulations. RESULTS: We presented initial results of large field pMBRT aperture and the film measurements indicates the effect of source-to-isocenter distance = 930 cm in Y proton scanning direction. Consistent with TOPAS MC simulation, the dose uniformity of pMBRT field <2 cm is demonstrated to be better than 2%, rendering its suitability for pre-clinical studies. Subsequently, we developed the second generation of aperture with five slits and characterized the aperture with film dosimetry studies and compared the results to the benchmark MC. Comprehensive film measurements were also performed to evaluate the effect of divergence, air gap and gantry-angle dependency and repeatability and revealing a consistent performance within 5%. Furthermore, the 2D gamma analysis indicated a passing rate exceeding 99% using 3% dose difference and 0.2 mm distance agreement criteria. We also establish the peak valley dose ratio and the depth dose profile measurements, and the results are within 10% from MC simulation. CONCLUSIONS: We have developed the first pMBRT system tailored for a single-gantry proton facility, which has demonstrated accuracy in benchmark with MC simulations, and allows for efficient plug-and-play setup, emphasizing efficiency.

An artificial neural network based approach for predicting the proton beam spot dosimetric characteristics of a pencil beam scanning technique.

Utilising Machine Learning (ML) models to predict dosimetric parameters in pencil beam scanning proton therapy presents a promising and practical approach. The study developed Artificial Neural Network (ANN) models to predict proton beam spot size and relative positional errors using 9000 proton spot data. The irradiation log files as input variables and corresponding scintillation detector measurements as the label values. The ANN models were developed to predict six variables: spot size in thex-axis,y-axis, major axis, minor axis, and relative positional errors in thex-axis andy-axis. All ANN models used a Multi-layer perception (MLP) network using one input layer, three hidden layers, and one output layer. Model performance was validated using various statistical tools. The log file recorded spot size and relative positional errors, which were compared with scintillator-measured data. The Root Mean Squared Error (RMSE) values for the x-spot and y-spot sizes were 0.356 mm and 0.362 mm, respectively. Additionally, the maximum variation for the x-spot relative positional error was 0.910 mm, while for the y-spot, it was 1.610 mm. The ANN models exhibit lower prediction errors. Specifically, the RMSE values for spot size prediction in the x, y, major, and minor axes are 0.053 mm, 0.049 mm, 0.053 mm, and 0.052 mm, respectively. Additionally, the relative spot positional error prediction model for the x and y axes yielded maximum errors of 0.160 mm and 0.170 mm, respectively. The normality of models was validated using the residual histogram and Q-Q plot. The data over fit, and bias were tested using K (k = 5) fold cross-validation, and the maximum RMSE value of the K fold cross-validation among all the six ML models was less than 0.150 mm (R-Square 0.960). All the models showed excellent prediction accuracy. Accurately predicting beam spot size and positional errors enhances efficiency in routine dosimetric checks.

Experimental comparison of cylindrical and plane parallel ionization chambers for reference dosimetry in continuous and pulsed scanned proton beams.

Objective. In this experimental work we compared the determination of absorbed dose to water using four ionization chambers (ICs), a PTW-34045 Advanced Markus, a PTW-34001 Roos, an IBA-PPC05 and a PTW-30012 Farmer, irradiated under the same conditions in one continuous- and in two pulsed-scanned proton beams.Approach. The ICs were positioned at 2 cm depth in a water phantom in four square-field single-energy scanned-proton beams with nominal energies between 80 and 220 MeV and in the middle of 10 × 10 × 10 cm3dose cubes centered at 10 cm or 12.5 cm depth in water. The water-equivalent thickness (WET) of the entrance window and the effective point of measurement was considered when positioning the plane parallel (PP) ICs and the cylindrical ICs, respectively. To reduce uncertainties, all ICs were calibrated at the same primary standards laboratory. We used the beam quality (kQ) correction factors for the ICs under investigation from IAEA TRS-398, the newly calculated Monte Carlo (MC) values and the anticipated IAEA TRS-398 updated recommendations.Main results. Dose differences among the four ICs ranged between 1.5% and 3.7% using both the TRS-398 and the newly recommendedkQvalues. The spread among the chambers is reduced with the newlykQvalues. The largest differences were observed between the rest of the ICs and the IBA-PPC05 IC, obtaining lower dose with the IBA-PPC05.Significance. We provide experimental data comparing different types of chambers in different proton beam qualities. The observed dose differences between the ICs appear to be related to inconsistencies in the determination of thekQvalues. For PP ICs, MC studies account for the physical thickness of the entrance window rather than the WET. The additional energy loss that the wall material invokes is not negligible for the IBA-PPC05 and might partially explain the lowkQvalues determined for this IC. To resolve this inconsistency and to benchmark MC values,kQvalues measured using calorimetry are needed.