The 200-MeV proton therapy project at the Paul Scherrer Institute: conceptual design and practical realization

Dosimetric study of synchrotron rapid beam off control and skip spot function for high beam intensity proton therapy

BACKGROUND

All Hitachi proton pencil beam scanning facilities currently use discrete spot scanning (DSS). Mayo Clinic Florida (MCF) is installing a Hitachi particle therapy system with advanced technologies, including fast scan speeds, high beam intensity, rapid beam off control (RBOC), a skip spot function, and proton pencil beam scanning using dose driven continuous scanning (DDCS). A potential concern of RBOC is the generation of a shoulder at the end of the normal spot delivery due to a flap spot (FS) with a flap dose (FD), which has been investigated for carbon synchrotron but not for proton delivery. While investigated, for instance, for Hitachi's installation at MCF, this methodology could be applicable for all future high intensity proton deliveries.

PURPOSE

No Hitachi proton facility currently uses the proposed RBOC. This study aimed to understand the dosimetric impact of proton FD at MCF by simulating the FS with a Hitachi proton machine in research mode, reflecting the higher proton intensities expected with RBOC at MCF.

METHOD

Experiments were conducted to simulate MCF RBOC at Kyoto Prefecture University of Medicine (KPUM) in research mode, reducing delay time (Td) from 1.5 ms to 0.1 ms. 5,000 contiguous spots were delivered on the central axis for proton energies of 70.2, 142.5, and 220.0 MeV; at normal, high dose rate (HDR), and ultra-high dose rate (uHDR) intensities; and at vertical and horizontal gantry angles for different Td. Measurements were taken using a fast oscilloscope and the nozzle's spot position monitor (SPM) and dose monitor (DM). A model was developed to predict FD dependence on beam intensity and assess the dosimetric impact for prostate and brain treatment plans. Two simulation types were planned: a flap DSS plan with FS at every spot and a flap DDCS plan with FS only at the end of each layer.

RESULT

FD was observed for RBOC with Td = 0.1 ms, showing no gantry angle dependence. FD increased with higher delayed dose rate (DDR), that is, beam intensity. The planning study showed dose volume histogram deterioration with increased FD compared to the clinical plan, but it was only significant for uHDR intensities. Deterioration was marginal in flap DSS plans for the HDR intensities planned at MCF, and flap DDCS plans were even less sensitive than flap DSS plans.

CONCLUSION

MCF is installing proton DDCS with higher beam intensities, a skip spot function, and fast beam-off control. The resulting FD had an insignificant impact on dose distribution for two patient plans with both DSS and DDCS at the anticipated MCF intensities. However, significant dependence was observed in the case of uHDR. A method to measure the position and dose of the FS during commissioning is described in addition to recommendations for regular QA and log-based proton patient-specific quality assurance.

Sectored single-energy volumetric-modulated proton arc therapy (VPAT): A preliminary multi-disease-site concept study

PURPOSE

To explore the feasibility of a novel intensity-modulated proton arc technique that uses a single-energy beam from the cyclotron. The beam energy is externally modulated at each gantry angle by a tertiary energy modulator (EM). We hypothesize that irradiating in an arc without requiring an energy change from the cyclotron will achieve a faster delivery (main advantage of our technique) while keeping clinically desirable dosimetric results.

METHODS

In a retrospective cohort of four patients with female pelvis, prostate, lung, and brain cancers, we investigated our volumetric-modulated proton arc therapy (VPAT) technique. Arcs were simulated by sectors of 1°-spaced static beams. Keeping the energy requested from the cyclotron the same for each entire arc was supported by a predesigned EM placed in front of the nozzle. As a feasibility measure, EM thicknesses were calculated. Delivery times and doses to targets and organs at risk (OARs) were compared to those of the clinical plans.

RESULTS

VPAT plans were comparable to their clinical counterparts in achieving target dose conformity, being robust to uncertainties, and meeting clinical dose-volume constraints. Cyclotron energies for the four cases were within 159-220 MeV, and energy modulation range was 69-100 MeV, equivalent to 13-19 cm of water-equivalent thickness (WET). Plan delivery times were reduced from > 5 min in our clinical practice to < 3.5 min in VPAT.

CONCLUSION

For the evaluated plans, the novel VPAT approach achieved shorter delivery times without sacrificing robustness, OAR sparing or target coverage. VPAT's EMs had WETs implementable in a clinical setup.

Craniospinal irradiation using pencil beam Scanning: The PSI experience

INTRODUCTION

We present the dosimetric evaluation of craniospinal irradiation (CSI) treatments delivered with protons at Paul Scherrer Institute (PSI), with special focus on local recurrences and late toxicity outcome.

METHODS

This study included 71 children, adolescents and young adults (c-AYA), who received or intended to receive (3 patients, pts) CSI using PBS-PT at PSI between 2004 and January 2021. The most frequent primary tumours were: medulloblastoma (42 pts), ependymoma (8 pts) and germ cell tumors (6 pts). The patients were treated prone on Gantry1 (G1; 22 pts) up to 2017, and afterwards supine on Gantry2 (G2; 49 pts). Accuracy of prone vs. supine setup was evaluated. Nine patients received CSI for local failure (LF) after a first course of local fractionated radiation therapy (RT). For 59/71 patients (excluding three patients not receiving PBS-PT CSI and nine preirradiated) CSI plans were compared considering gantry and planning technique. Detailed analysis of the full treatment (CSI and boost series) was performed for 8 patients presenting with LFs (4 of them presented also distal failure) and for selected patients presenting with late toxicity (G2 to G4) or asymptomatic radiation-induced radiological findings.

RESULTS

Supine positioning resulted in lower systematic and random errors as compared to prone (0.25 mm and 0.4 mm systematic errors respectively for supine and prone; random errors in PA direction reduced from 1.8 mm for prone to 1.4 mm for supine).

CONCLUSIONS

LFs were not correlated with potential dose inaccuracies or lack of robustness and no correlation of toxicities to enhanced LET have been observed.

Efficient proton arc optimization and delivery through energy layer pre-selection and post-filtering

BACKGROUND

Proton arc therapy (PAT) has emerged as a promising approach for improving dose distribution, but also enabling simpler and faster treatment delivery in comparison to conventional proton treatments. However, the delivery speed achievable in proton arc relies on dedicated algorithms, which currently do not generate plans with a clear speed-up and sometimes even result in increased delivery time.

PURPOSE

This study aims to address the challenge of minimizing delivery time through a hybrid method combining a fast geometry-based energy layer (EL) pre-selection with a dose-based EL filtering, and comparing its performance to a baseline approach without filtering.

METHODS

Three methods of EL filtering were developed: unrestricted, switch-up (SU), and switch-up gap (SU gap) filtering. The unrestricted method filters the lowest weighted EL while the SU gap filtering removes the EL around a new SU to minimize the gantry rotation braking. The SU filtering removes the lowest weighted group of EL that includes a SU. These filters were combined with the RayStation dynamic proton arc optimization framework energy layer selection and spot assignment (ELSA). Four bilateral oropharyngeal and four lung cancer patients' data were used for evaluation. Objective function values, target coverage robustness, organ-at-risk doses and normal tissue complication probability evaluations, as well as comparisons to intensity-modulated proton therapy (IMPT) plans, were used to assess plan quality.

RESULTS

The SU gap filtering algorithm performed best in five out of the eight cases, maintaining plan quality within tolerance while reducing beam delivery time, in particular for the oropharyngeal cohort. It achieved up to approximately 22% and 15% reduction in delivery time for oropharyngeal and lung treatment sites, respectively. The unrestricted filtering algorithm followed closely. In contrast, the SU filtering showed limited improvement, suppressing one or two SU without substantial delivery time shortening. Robust target coverage was kept within 1% of variation compared to the PAT baseline plan while organs-at-risk doses slightly decreased or kept about the same for all patients.

CONCLUSIONS

This study provides insights to accelerate PAT delivery without compromising plan quality. These advancements could enhance treatment efficiency and patient throughput.

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.

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.

Streamlined pin-ridge-filter design for single-energy proton FLASH planning

BACKGROUND

FLASH radiotherapy (FLASH-RT) with ultra-high dose rate has yielded promising results in reducing normal tissue toxicity while maintaining tumor control. Planning with single-energy proton beams modulated by ridge filters (RFs) has been demonstrated feasible for FLASH-RT.

PURPOSE

This study explored the feasibility of a streamlined pin-shaped RF (pin-RF) design, characterized by coarse resolution and sparsely distributed ridge pins, for single-energy proton FLASH planning.

METHODS

An inverse planning framework integrated within a treatment planning system was established to design streamlined pin RFs for single-energy FLASH planning. The framework involves generating a multi-energy proton beam plan using intensity-modulated proton therapy (IMPT) planning based on downstream energy modulation strategy (IMPT-DS), followed by a nested pencil-beam-direction-based (PBD-based) spot reduction process to iteratively reduce the total number of PBDs and energy layers along each PBD for the IMPT-DS plan. The IMPT-DS plan is then translated into the pin-RFs and the single-energy beam configurations for IMPT planning with pin-RFs (IMPT-RF). This framework was validated on three lung cases, quantifying the FLASH dose of the IMPT-RF plan using the FLASH effectiveness model. The FLASH dose was then compared to the reference dose of a conventional IMPT plan to measure the clinical benefit of the FLASH planning technique.

RESULTS

The IMPT-RF plans closely matched the corresponding IMPT-DS plans in high dose conformity (conformity index of <1.2), with minimal changes in V7Gy and V7.4 Gy for the lung (<3%) and small increases in maximum doses (Dmax) for other normal structures (<3.4 Gy). Comparing the FLASH doses to the doses of corresponding IMPT-RF plans, drastic reductions of up to nearly 33% were observed in Dmax for the normal structures situated in the high-to-moderate-dose regions, while negligible changes were found in Dmax for normal structures in low-dose regions. Positive clinical benefits were seen in comparing the FLASH doses to the reference doses, with notable reductions of 21.4%-33.0% in Dmax for healthy tissues in the high-dose regions. However, in the moderate-to-low-dose regions, only marginal positive or even negative clinical benefit for normal tissues were observed, such as increased lung V7Gy and V7.4 Gy (up to 17.6%).

CONCLUSIONS

A streamlined pin-RF design was developed and its effectiveness for single-energy proton FLASH planning was validated, revealing positive clinical benefits for the normal tissues in the high dose regions. The coarsened design of the pin-RF demonstrates potential advantages, including cost efficiency and ease of adjustability, making it a promising option for efficient production.

Complications after proton radiotherapy in children, focusing on severe late complications. A complete Swedish cohort 2008-2019

BACKGROUND

Proton radiotherapy (RT) is an attractive tool to deliver local therapy with minimal dose to uninvolved tissue, however, not suitable for all patients. The aim was to explore complications, especially severe late complications (grades 3-4), following proton RT delivered to a complete Swedish cohort of paediatric patients aged <18 years treated 2008-2019.

MATERIAL AND METHODS

Data was downloaded from a national registry. Complications with a possible causation with RT are reported. Proton treatments until July 2015 was performed with a fixed horizontal 172 MeV beam (The Svedberg Laboratory (TSL), Uppsala) in a sitting position and thereafter with gantry-based pencil-beam scanning technique (Skandion Clinic, Uppsala) in a supine position.

RESULTS

219 courses of proton RT (77 at TSL and 142 at Skandion) were delivered to 212 patients (mean age 9.2 years) with various tumour types (CNS tumours 58%, sarcomas 26%, germ cell tumours 7%). Twenty-five patients had severe acute complications (skin, mucous membrane, pharynx/oesophagus, larynx, upper gastrointestinal canal, lower gastrointestinal canal, eyes, ears). Fifteen patients had severe late complications; with increased proportion over time: 4% at 1-year follow-up (FU), 5% at 3-year, 11% at 5-year. Organs affected were skin (1 patient), subcutaneous tissue (4), salivary glands (1), upper GI (1), bone (7), joints (2), CNS (2), PNS (1), eyes (1) and ears (5). Twenty-one of the 28 patients with 10-year FU had at least one late complication grades 1-4 and fourteen of them had more than one (2-5 each).

CONCLUSION

The most important result of our study is the relatively low proportion of severe late complications, comparable with other proton studies on various tumours. Furthermore, the numbers of late complications are lower than our own data set on a mixed population of photon and proton treated paediatric patients, assuring the safety of using proton therapy also in the clinical practice.

Selecting Optimal Proton Pencil Beam Scanning Plan Parameters to Reduce Dose Discrepancy between Discrete Spot Plan and Continuous Scanning: A Proof-of-Concept Study

Pencil beam scanning delivered with continuous scanning has several advantages over conventional discrete spot scanning. Such advantages include improved beam delivery efficiency and reduced beam delivery time. However, a move dose is delivered between consecutive spots with continuous scanning, and current treatment planning systems do not take this into account. Therefore, continuous scanning and discrete spot plans have an inherent dose discrepancy. Using the operating parameters of the state-of-the-art particle therapy system, we conducted a proof-of-concept study in which we systematically generated 28 plans for cubic targets with different combinations of plan parameters and simulated the dose discrepancies between continuous scanning and a planned one. A nomograph to guide the selection of plan parameters was developed to reduce the dose discrepancy. The effectiveness of the nomograph was evaluated with two clinical cases (one prostate and one liver). Plans with parameters guided by the nomograph decreased dose discrepancy than those used standard plan parameters. Specifically, the 2%/2 mm gamma passing rate increased from 96.3% to 100% for the prostate case and from 97.8% to 99.7% for the liver case. The CTV DVH root mean square error decreased from 2.2% to 0.2% for the prostate case and from 1.8% to 0.9% for the liver case. The decreased dose discrepancy may allow the relaxing of the delivery constraint for some cases, leading to greater benefits in continuous scanning. Further investigation is warranted.

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

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

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.

Long Term Outcome and Quality of Life of Intracranial Meningioma Patients Treated with Pencil Beam Scanning Proton Therapy

The aim of this study was to assess the clinical outcome, including QoL, of patients with intracranial meningiomas WHO grade 1-3 who were treated with Pencil Beam Scanning Proton Therapy (PBS PT) between 1997 and 2022. Two hundred patients (median age 50.4 years, 70% WHO grade 1) were analyzed. Acute and late side effects were classified according to CTCAE version 5.0. Time to event data were calculated. QoL was assessed descriptively by the EORTC-QLQ-C30 and BN20 questionnaires. With a median follow-up of 65 months (range: 3.8-260.8 months) the 5 year OS was 95.7% and 81.8% for WHO grade 1 and grade 2/3, respectively (p < 0.001). Twenty (10%) local failures were observed. Failures occurred significantly (p < 0.001) more frequent in WHO grade 2 or 3 meningioma (WHO grade 1: n = 7, WHO grade 2/3: n = 13), in patients with multiple meningiomas (p = 0.005), in male patients (p = 0.005), and when PT was initiated not as upfront therapy (p = 0.011). There were no high-grade toxicities in the majority (n = 176; 88%) of patients. QoL was assessed for 83 (41.5%) patients and for those patients PT did not impacted QoL negatively during the follow-up. In summary, we observed very few local recurrences of meningiomas after PBS PT, a stable QoL, and a low rate of high-grade toxicity.

General Purpose Transistor Characterized as Dosimetry Sensor of Proton Beams

A commercial pMOS transistor (MOSFET), 3N163 from Vishay (USA), has been characterized as a low-energy proton beam dosimeter. The top of the samples' housing has been removed to guarantee that protons reached the sensitive area, that is, the silicon die. Irradiations took place at the National Accelerator Centre (Seville, Spain). During irradiations, the transistors were biased to improve the sensitivity, and the silicon temperature was monitored activating the parasitic diode of the MOSFET. Bias voltages of 0, 1, 5, and 10 V were applied to four sets of three transistors, obtaining an averaged sensitivity that was linearly dependent on this voltage. In addition, the short-fading effect was studied, and the uncertainty of this effect was obtained. The bias voltage that provided an acceptable sensitivity, (11.4 ± 0.9) mV/Gy, minimizing the uncertainty due to the fading effect (-0.09 ± 0.11) Gy was 1 V for a total absorbed dose of 40 Gy. Therefore, this off-the-shelf electronic device presents promising characteristics as a dosimeter sensor for proton beams.

Transformative Technology for FLASH Radiation Therapy

The general concept of radiation therapy used in conventional cancer treatment is to increase the therapeutic index by creating a physical dose differential between tumors and normal tissues through precision dose targeting, image guidance, and radiation beams that deliver a radiation dose with high conformality, e.g., protons and ions. However, the treatment and cure are still limited by normal tissue radiation toxicity, with the corresponding side effects. A fundamentally different paradigm for increasing the therapeutic index of radiation therapy has emerged recently, supported by preclinical research, and based on the FLASH radiation effect. FLASH radiation therapy (FLASH-RT) is an ultra-high-dose-rate delivery of a therapeutic radiation dose within a fraction of a second. Experimental studies have shown that normal tissues seem to be universally spared at these high dose rates, whereas tumors are not. While dose delivery conditions to achieve a FLASH effect are not yet fully characterized, it is currently estimated that doses delivered in less than 200 ms produce normal-tissue-sparing effects, yet effectively kill tumor cells. Despite a great opportunity, there are many technical challenges for the accelerator community to create the required dose rates with novel compact accelerators to ensure the safe delivery of FLASH radiation beams.

Inter- and intrafractional 4D dose accumulation for evaluating ΔNTCP robustness in lung cancer

BACKGROUND AND PURPOSE

Model-based selection of proton therapy patients relies on a predefined reduction in normal tissue complication probability (NTCP) with respect to photon therapy. The decision is necessarily made based on the treatment plan, but NTCP can be affected when the delivered treatment deviates from the plan due to delivery inaccuracies. Especially for proton therapy of lung cancer, this can be important because of tissue density changes and, with pencil beam scanning, the interplay effect between the proton beam and breathing motion.

MATERIALS AND METHODS

In this work, we verified whether the expected benefit of proton therapy is retained despite delivery inaccuracies by reconstructing the delivered treatment using log-file based dose reconstruction and inter- and intrafractional accumulation. Additionally, the importance of two uncertain parameters for treatment reconstruction, namely deformable image registration (DIR) algorithm and α/β ratio, was assessed.

RESULTS

The expected benefit or proton therapy was confirmed in 97% of all studied cases, despite regular differences up to 2 percent point (p.p.) NTCP between the delivered and planned treatments. The choice of DIR algorithm affected NTCP up to 1.6 p.p., an order of magnitude higher than the effect of α/β ratio.

CONCLUSION

For the patient population and treatment technique employed, the predicted clinical benefit for patients selected for proton therapy was confirmed for 97.0% percent of all cases, although the NTCP based proton selection was subject to 2 p.p. variations due to delivery inaccuracies.

Biological effects of cancer stem cells irradiated by charged particle: a systematic review of in vitro studies

PURPOSE

The existence of cancer stem cells (CSCs) is closely related to tumor recurrence, metastasis, and resistance to chemoradiotherapy. In addition, given the unique physical and biological advantages of charged particle, we hypothesized that charged particle irradiation would produce strong killing effects on CSCs. The purpose of our systematic review is to evaluate the biological effects of CSCs irradiated by charged particle, including proliferation, invasion, migration, and changes in the molecular level.

METHODS

We searched PubMed, EMBASE, and Web of Science until 17 march 2022 according to the key words. Included studies have to be vitro studies of CSCs irradiated by charged particle. Outcomes included one or more of radiation sensitivity, proliferation, metastasis, invasion, and molecular level changes, like DNA damage after been irradiated.

RESULTS

Eighteen studies were included in the final analysis. The 18 articles include 12-carbon ion irradiation, 4-proton irradiation, 1 α-particle irradiation, 1-carbon ion combine proton irradiation.

CONCLUSION

Through the extraction and analysis of data, we came to this conclusion: CSCs have obvious radio-resistance compared with non-CSCs, and charged particle irradiation or in combination with drugs could overcome this resistance, specifically manifested in inhibiting CSCs' proliferation, invasion, migration, and causing more and harder to repair DNA double-stranded breaks (DSB) of CSCs.

Dosimetric impact of systematic spot position errors in spot scanning proton therapy of head and neck tumor

Purpose

The spot position is an important beam parameter in the quality assurance of scanning proton therapy. In this study, we investigated dosimetric impact of systematic 15 spot position errors (SSPE) in spot scanning proton therapy using three types of optimization methods of head and neck tumor.

Materials and Methods

The planning simulation was performed with ± 2 mm model SSPE in the X and Y directions. Treatment plans were created using intensity-modulated proton therapy (IMPT) and single-field uniform dose (SFUD). IMPT plans were created by two optimization methods: with worst-case optimization (WCO-IMPT) and without (IMPT). For clinical target volume (CTV), D95%, D50%, and D2cc were used for analysis. For organs at risk (OAR), Dmean was used to analyze the brain, cochlea, and parotid, and Dmax was used to analyze brainsetem, chiasm, optic nerve, and cord.

Results

For CTV, the variation (1 standard deviation) of D95% was ± 0.88%, 0.97% and 0.97% to WCO-IMPT, IMPT, and SFUD plan. The variation of D50% and D2cc of CTV showed <0.5% variation in all plans. The dose variation due to SSPE was larger in OAR, and worst-case optimization reduced the dose variation, especially in Dmax. The analysis results showed that SSPE has little impact on SFUD.

Conclusions

We clarified the impact of SSPE on dose distribution for three optimization methods. SFUD was shown to be a robust treatment plan for OARs, and the WCO can be used to increase robustness to SSPE in IMPT.

The effect of voxelization in Monte Carlo simulation to validate Bragg peak characteristics for a pencil proton beam

Background

The purpose of this research was to show how the Bragg peak (BP) characteristics were affected by changing the voxel size in longitudinal and transverse directions in Monte Carlo (MC) simulations by using Geant4 and to calculate BP characteristics accurately by considering the voxel size effect for 68 MeV and 235.81 MeV.

Materials and methods

Different interpolation techniques were applied to simulation data to find the closest results to the experimental data.

Results

When the x-size of the voxel was increased 2 times at low energy, the maximum dose increase in the entrance and plateau regions were 17.8% and 17%, respectively, while BP curve shifted to the shallower region, resulting in a 0.5 mm reduction in the curable tumor width (W_80pd). At high energy, the maximum dose increase at the entrance and plateau regions were 0.4% and 0.6%, respectively, while it was observed that W_80pd did not change. When the y-z sizes of the voxel were increased 2 times at low energy, the maximum dose reduction at the entrance and plateau regions was 3.4%, but no change was observed in W_80pd. At high energy, when the y-z sizes of the voxel were increased 2.2 times, the maximum dose reduction at the entrance and plateau regions were 8.9% and 9.1%, respectively, while W_80pd increased by 0.5 mm. When linear, cubic spline, and Akima interpolations were applied to the simulation data, it was found that the results closest to the experimental data were obtained for Akima interpolations for both energies.

Conclusion

it has been shown that the voxel size effect for the longitudinal direction was more effective at low energy than at high energy. However, the voxel size effect for the transverse direction was more effective for high energy.

Ultra-high dose rate dosimetry for pre-clinical experiments with mm-small proton fields

PURPOSE

To characterize an experimental setup for ultra-high dose rate (UHDR) proton irradiations, and to address the challenges of dosimetry in millimetre-small pencil proton beams.

METHODS

At the PSI Gantry 1, high-energy transmission pencil beams can be delivered to biological samples and detectors up to a maximum local dose rate of ∼9000 Gy/s. In the presented setup, a Faraday cup is used to measure the delivered number of protons up to ultra-high dose rates. The response of transmission ion-chambers, as well as of different field detectors, was characterized over a wide range of dose rates using the Faraday cup as reference.

RESULTS

The reproducibility of the delivered proton charge was better than 1 % in the proposed experimental setup. EBT3 films, Al₂O₃:C optically stimulated luminescence detectors and a PTW microDiamond were used to validate the predicted dose. Transmission ionization chambers showed significant volume ion-recombination (>30 % in the tested conditions) which can be parametrized as a function of the maximum proton current density. Over the considered range, EBT3 films, inorganic scintillator-based screens and the PTW microDiamond were demonstrated to be dose rate independent within ±3 %, ±1.8 % and ±1 %, respectively.

CONCLUSIONS

Faraday cups are versatile dosimetry instruments that can be used for dose estimation, field detector characterization and on-line dose verification for pre-clinical experiments in UHDR proton pencil beams. Among the tested detectors, the commercial PTW microDiamond was found to be a suitable option to measure real time the dosimetric properties of narrow pencil proton beams for dose rates up to 2.2 kGy/s.

Validation of pencil beam scanning proton therapy with multi-leaf collimator calculated by a commercial Monte Carlo dose engine

This study aimed to evaluate the clinical beam commissioning results and lateral penumbra characteristics of our new pencil beam scanning (PBS) proton therapy using a multi-leaf collimator (MLC) calculated by use of a commercial Monte Carlo dose engine. Eighteen collimated uniform dose plans for cubic targets were optimized by the RayStation 9A treatment planning system (TPS), varying scan area, modulation widths, measurement depths, and collimator angles. To test the patient-specific measurements, we also created and verified five clinically realistic PBS plans with the MLC, such as the liver, prostate, base-of-skull, C-shape, and head-and-neck. The verification measurements consist of the depth dose (DD), lateral profile (LP), and absolute dose (AD). We compared the LPs and ADs between the calculation and measurements. For the cubic plans, the gamma index pass rates (γ-passing) were on average 96.5% ± 4.0% at 3%/3 mm for the DD and 95.2% ± 7.6% at 2%/2 mm for the LP. In several LP measurements less than 75 mm depths, the γ-passing deteriorated (increased the measured doses) by less than 90% with the scattering such as the MLC edge and range shifter. The deteriorated γ-passing was satisfied by more than 90% at 2%/2 mm using uncollimated beams instead of collimated beams except for three planes. The AD differences and the lateral penumbra width (80%-20% distance) were within ±1.9% and ± 1.1 mm, respectively. For the clinical plan measurements, the γ-passing of LP at 2%/2 mm and the AD differences were 97.7% ± 4.2% on average and within ±1.8%, respectively. The measurements were in good agreement with the calculations of both the cubic and clinical plans inserted in the MLC except for LPs less than 75 mm regions of some cubic and clinical plans. The calculation errors in collimated beams can be mitigated by substituting uncollimated beams.

Universal and dynamic ridge filter for pencil beam scanning particle therapy: a novel concept for ultra-fast treatment delivery

Objective. In pencil beam scanning particle therapy, a short treatment delivery time is paramount for the efficient treatment of moving targets with motion mitigation techniques (such as breath-hold, rescanning, and gating). Energy and spot position change time are limiting factors in reducing treatment time. In this study, we designed a universal and dynamic energy modulator (ridge filter, RF) to broaden the Bragg peak, to reduce the number of energies and spots required to cover the target volume, thus lowering the treatment time. Approach. Our RF unit comprises two identical RFs placed just before the isocenter. Both RFs move relative to each other, changing the Bragg peak’s characteristics dynamically. We simulated different Bragg peak shapes with the RF in Monte Carlo simulation code (TOPAS) and validated them experimentally. We then delivered single-field plans with 1 Gy/fraction to different geometrical targets in water, to measure the dose delivery time using the RF and compare it with the clinical settings. Main results. Aligning the RFs in different positions produces different broadening in the Bragg peak; we achieved a maximum broadening of 2.5 cm. With RF we reduced the number of energies in a field by more than 60%, and the dose delivery time by 50%, for all geometrical targets investigated, without compromising the dose distribution transverse and distal fall-off. Significance. Our novel universal and dynamic RF allows for the adaptation of the Bragg peak broadening for a spot and/or energy layer based on the requirement of dose shaping in the target volume. It significantly reduces the number of energy layers and spots to cover the target volume, and thus the treatment time. This RF design is ideal for ultra-fast treatment delivery within a single breath-hold (5–10 s), efficient delivery of motion mitigation techniques, and small animal irradiation with ultra-high dose rates (FLASH).

Investigation of beam delivery time for synchrotron-based proton pencil beam scanning system with novel scanning mode

Objective.To investigate synchrotron-based proton pencil beam scanning (PBS) beam delivery time (BDT) using novel continuous scanning mode.Approach.A BDT calculation model was developed for the Hitachi particle therapy system. The model was validated against the measured BDT of 36 representative clinical proton PBS plans with discrete spot scanning (DSS) in the current Hitachi proton therapy system. BDTs were calculated with the next generation using Mayo Clinic Florida system operating parameters for conventional DSS, and novel dose driven continuous scanning (DDCS). BDTs of DDCS with and without Break Spots were investigated.Main results.For DDCS without Break Spots, the use of Stop Ratio to control the transit dose largely reduced the beam intensity and consequently, severely prolonged the BDT. DDCS with Break Spots was able to maintain a sufficiently high beam intensity while controlling transit dose. In DDCS with Break Spots, tradeoffs were made between beam intensity and number of Break Spots. Therefore, BDT decreased with increased beam intensity but reached a plateau for beam intensity larger than 10 MU s^-1. Averaging over all clinical plans, BDT was reduced by 10% for DDCS with Break Spots compared to DSS.Significance.DDCS with Break Spots reduced BDT. DDCS has the potential to further reduce BDT under the ideal scenario which requests both stable beam intensity extraction and accurately modelling the transit dose. Further investigation is warranted.

Ultra-fast pencil beam scanning proton therapy for locally advanced non-small-cell lung cancers: field delivery within a single breath-hold

PURPOSE

The use of motion mitigation techniques such as breath-hold can reduce the dosimetric uncertainty of lung cancer proton therapy. We studied the feasibility of pencil beam scanning (PBS) proton therapy field delivery within a single breath-hold at PSI's Gantry 2.

METHODS

In PBS proton therapy, the delivery time for a field is determined by the beam-on time and the dead time between proton spots (the time required to change the energy and/or lateral position). We studied ways to reduce beam-on and lateral scanning time, without sacrificing dosimetric plan quality, aiming at a single field delivery time of 15 seconds at maximum. We tested this approach on 10 lung cases with varying target volumes. To reduce the beam-on time, we increased the beam current at the isocenter by developing new beam optics for PSI's PROSCAN beamline and Gantry 2. To reduce the dead time between the spots, we used spot-reduced plan optimization.

RESULTS

We found that it is possible to achieve conventional fractionated (2 Gy(RBE)/fraction) and hypofractionated (6 Gy(RBE)/fraction) field delivery times within a single breath-hold (<15 sec) for a variety non-small-cell lung cancer cases.

CONCLUSION

In summary, the combination of spot reduction and improved beam line transmission is a promising approach for the treatment of mobile tumours within clinically achievable breath-hold durations.

The dose accumulation and the impact of deformable image registration on dose reporting parameters in a moving patient undergoing proton radiotherapy

INTRODUCTION

Potential changes in patient anatomy during proton radiotherapy may lead to a deviation of the delivered dose. A dose estimate can be computed through a deformable image registration (DIR) driven dose accumulation. The present study evaluates the accumulated dose uncertainties in a patient subject to an inadvertent breathing associated motion.

MATERIALS AND METHODS

A virtual lung tumour was inserted into a pair of single participant landmark annotated computed tomography images depicting opposite breathing phases, with the deep inspiration breath-hold the planning reference and the exhale the off-reference geometry. A novel Monte Carlo N-Particle, Version 6 (MCNP6) dose engine was developed, validated and used in treatment plan optimization. Three DIR methods were compared and used to transfer the exhale simulated dose to the reference geometry. Dose conformity and homogeneity measures from International Committee on Radioactivity Units and Measurements (ICRU) reports 78 and 83 were evaluated on simulated dose distributions registered with different DIR algorithms.

RESULTS

The MCNP6 dose engine handled patient-like geometries in reasonable dose calculation times. All registration methods were able to align image associated landmarks to distances, comparable to voxel sizes. A moderate deterioration of ICRU measures was encountered in comparing doses in on and off-reference anatomy. There were statistically significant DIR driven differences in ICRU measures, particularly a 10% difference in the relative D98% for planning tumour volume and in the 3 mm/3% gamma passing rate.

CONCLUSIONS

T he dose accumulation over two anatomies resulted in a DIR driven uncertainty, important in reporting the associated ICRU measures for quality assurance.

Increase of the transmission and emittance acceptance through a cyclotron‐based proton therapy gantry

Purpose

In proton therapy, the gantry, as the final part of the beamline, has a major effect on beam intensity and beam size at the isocenter. Most of the conventional beam optics of cyclotron‐based proton gantries have been designed with an imaging factor between 1 and 2 from the coupling point (CP) at the gantry entrance to the isocenter (patient location) meaning that to achieve a clinically desirable (small) beam size at isocenter, a small beam size is also required at the CP. Here we will show that such imaging factors are limiting the emittance which can be transported through the gantry. We, therefore, propose the use of large beam size and low divergence beam at the CP along with an imaging factor of 0.5 (2:1) in a new design of gantry beam optics to achieve substantial improvements in transmission and thus increase beam intensity at the isocenter.

Methods

The beam optics of our gantry have been re‐designed to transport higher emittance without the need of any mechanical modifications to the gantry beamline. The beam optics has been designed using TRANSPORT, with the resulting transmissions being calculated using Monte Carlo simulations (BDSIM code). Finally, the new beam optics have been tested with measurements performed on our Gantry 2 at PSI.

Results

With the new beam optics, we could maximize transmission through the gantry for a fixed emittance value. Additionally, we could transport almost four times higher emittance through the gantry compared to conventional optics, whilst achieving good transmissions through the gantry (>50%) with no increased losses in the gantry. As such, the overall transmission (cyclotron to isocenter) can be increased by almost a factor of 6 for low energies. Additionally, the point‐to‐point imaging inherent to the optics allows adjustment of the beam size at the isocenter by simply changing the beam size at the CP.

Conclusion

We have developed a new gantry beam optics which, by selecting a large beam size and low divergence at the gantry entrance and using an imaging factor of 0.5 (2:1), increases the emittance acceptance of the gantry, leading to a substantial increase in beam intensity at low energies. We expect that this approach could easily be adapted for most types of existing gantries.

Dosimetric Comparison of Various Spot Placement Techniques in Proton Pencil Beam Scanning

Purpose

To present quantitative dosimetric evaluations of five proton pencil beam spot placement techniques.

Materials and Methods

The spot placement techniques that were investigated include two grid-based (rectilinear grid and hexagonal grid, both commonly available in commercial planning systems) and three boundary-contoured (concentric contours, hybrid, and optimized) techniques. Treatment plans were created for two different target volumes, one spherical and one conical. An optimal set of planning parameters was defined for all treatment plans and the impact of spot placement techniques on the plan quality was evaluated in terms of lateral/distal dose falloff, normal tissue sparing, conformity and homogeneity of dose distributions, as well as total number of spots used.

Results

The results of this work highlight that for grid-based spot placement techniques, the dose conformity is dependent on target cross-sectional shape perpendicular to beam direction, which changes for each energy layer. This variable conformity problem is mitigated by using boundary contoured spot placement techniques. However, in the case of concentric contours, the conformity is improved but at the cost of decreased homogeneity inside the target. Hybrid and optimized spot placement techniques, which use contoured spots at the boundary and gridlike interior spot patterns, provide more uniform dose distributions inside the target volume while maintaining the improved dose conformity. The optimized spot placement technique improved target coverage, homogeneity of dose, and minimal number of spots. The dependence of these results on spot size is also presented for both target shapes.

Conclusion

This work illustrates that boundary-contoured spot placement techniques offer marked improvement in dosimetry metrics when compared to commercially available grid-based techniques for a range of proton scanned beam spot sizes.

Beam properties within the momentum acceptance of a clinical gantry beamline for proton therapy

PURPOSE

Energy changes in Pencil Beam Scanning (PBS) proton therapy can be a limiting factor in delivery time, hence limiting patient throughput and the effectiveness of motion mitigation techniques requiring fast irradiation. In this study, we investigate the feasibility of performing fast and continuous energy modulation within the momentum acceptance of a clinical beamline for proton therapy.

METHODS

The alternative use of a local beam degrader at the gantry coupling point has been compared with a more common upstream regulation. Focusing on clinically relevant parameters, a complete beam properties characterization has been carried out. In particular, the acquired empirical data allowed to model and parametrize the errors in range and beam current to deliver clinical treatment plans.

RESULTS

For both options, the local and the upstream degrader, depth-dose curves measured in water for off-momentum beams were only marginally distorted (γ(1%,1mm) > 90%) and the errors in the spot position were within the clinical tolerance, even though increasing at the boundaries of the investigated scan range. The impact on the beam size was limited for the upstream degrader while dedicated strategies could be required to tackle the beam broadening through the local degrader. Range correction models were investigated for the upstream regulation. The impaired beam transport required a dedicated strategy for fine range control and compensation of beam intensity losses. Our current parametrization based on empirical data allowed energy modulation within acceptance with range errors (median 0.05 mm) and transmission (median -14%) compatible with clinical operation and remarkably low average 27 ms dead time for small energy changes. The technique, tested for the delivery of a skull glioma treatment, resulted in high gamma pass rates at 1%,1mm compared to conventional deliveries in experimental measurements with about 45% reduction of the energy switching time when regulation could be performed within acceptance.

CONCLUSIONS

Fast energy modulation within beamline acceptance has potential for clinical applications and, when realized with an upstream degrader, does not require modification in the beamline hardware, therefore being potentially applicable in any running facility. Centers with slow energy switching time can particularly profit from such a technique for reducing dead time during treatment delivery. This article is protected by copyright. All rights reserved.

A quantitative FLASH effectiveness model to reveal potentials and pitfalls of high dose rate proton therapy

Purpose

In ultrahigh dose rate radiotherapy, the FLASH effect can lead to substantially reduced healthy tissue damage without affecting tumor control. Although many studies show promising results, the underlying biological mechanisms and the relevant delivery parameters are still largely unknown. It is unclear, particularly for scanned proton therapy, how treatment plans could be optimized to maximally exploit this protective FLASH effect.

Materials and Methods

To investigate the potential of pencil beam scanned proton therapy for FLASH treatments, we present a phenomenological model, which is purely based on experimentally observed phenomena such as potential dose rate and dose thresholds, and which estimates the biologically effective dose during FLASH radiotherapy based on several parameters. We applied this model to a wide variety of patient geometries and proton treatment planning scenarios, including transmission and Bragg peak plans as well as single‐ and multifield plans. Moreover, we performed a sensitivity analysis to estimate the importance of each model parameter.

Results

Our results showed an increased plan‐specific FLASH effect for transmission compared with Bragg peak plans (19.7% vs. 4.0%) and for single‐field compared with multifield plans (14.7% vs. 3.7%), typically at the cost of increased integral dose compared to the clinical reference plan. Similar FLASH magnitudes were found across the different treatment sites, whereas the clinical benefits with respect to the clinical reference plan varied strongly. The sensitivity analysis revealed that the threshold dose as well as the dose per fraction strongly impacted the FLASH effect, whereas the persistence time only marginally affected FLASH. An intermediate dependence of the FLASH effect on the dose rate threshold was found.

Conclusions

Our model provided a quantitative measure of the FLASH effect for various delivery and patient scenarios, supporting previous assumptions about potentially promising planning approaches for FLASH proton therapy. Positive clinical benefits compared to clinical plans were achieved using hypofractionated, single‐field transmission plans. The dose threshold was found to be an important factor, which may require more investigation.

Building a precise machine-specific time structure of the spot and energy delivery model for a cyclotron-based proton therapy system

OBJECTIVE

We proposed an experimental approach to build a precise machine-specific beam delivery time (BDT) prediction and delivery sequence model for standard, volumetric, and layer repainting delivery based on a cyclotron accelerator system. Approach Test fields and clinical treatment plans' log files were used to experimentally derive three main beam delivery parameters that impacted BDT: energy layer switching time (ELST), spot switching time (SSWT), and spot drill time (SDT). This derived machine-specific model includes standard, volumetric, and layer repainting delivery sequences. A total of 103 clinical treatment fields were used to validate the model.

MAIN RESULTS

The study found that ELST is not stochastic in this specific machine. Instead, it is actually the data transmission time or energy selection time, whichever takes longer. The validation showed that the accuracy of each component of the BDT matches well between machine log files and the model's prediction. The average total BDT was about (-0.74±3.33)% difference compared to the actual treatment log files, which is improved from the current commercial proton therapy system's prediction (67.22%±26.19%).

SIGNIFICANCE

An accurate BDT prediction and delivery sequence model was established for an cyclotron-based proton therapy system IBA ProteusPLUS®. Most institutions could adopt this method to build a machine-specific model for their own proton system.

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

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

Purpose

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

Methods

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

Results

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

Conclusions

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

Clinical Outcome of Sacral Chordoma Patients Treated with Pencil Beam Scanning Proton Therapy

AIMS

Sacral chordomas are locally aggressive, radio-resistant tumours. Proton therapy has the potential to deliver high radiation doses, which may improve the therapeutic ratio when compared with conventional radiotherapy. We assessed tumour control and radiation-induced toxicity in a cohort of sacral chordoma patients treated with definitive or postoperative pencil beam scanning proton therapy.

METHODS AND MATERIALS

Sixty patients with histologically proven sacral chordoma treated between November 1997 and October 2018 at the Paul Scherrer Institute with postoperative (n = 50) or definitive proton therapy (n = 10) were retrospectively analysed. Only 10 (17%) patients received combined photon radiotherapy and proton therapy. Survival rates were calculated using the Kaplan-Meier actuarial method. The Log-rank test was used to compare different functions for local control, freedom from distant recurrence and overall survival. Acute and late toxicity were assessed according to the Common Terminology Criteria for Adverse Events (CTCAE) v5.0.

RESULTS

The median follow-up was 48 months (range 4-186). Local recurrence occurred in 20 (33%) patients. The 4-year local control, freedom from distant recurrence and overall survival rates were 77%, 89% and 85%, respectively. On univariate analysis, subtotal resection/biopsy (P = 0.02), tumour extension restricted to bone (P = 0.01) and gross tumour volume >130 ml (P = 0.04) were significant predictors for local recurrence. On multivariate analysis, tumour extension restricted to bone (P = 0.004) and gross total resection (P = 0.02) remained independent favourable prognostic factors for local recurrence. Twenty-four (40%), 28 (47%) and eight (11%) patients experienced acute grade 1, 2 and 3 toxicities, respectively. The 4-year late toxicity-free survival was 91%. Two patients developed secondary malignancies to the bladder 3-7 years after proton therapy.

CONCLUSIONS

Our data indicate that pencil beam scanning proton therapy for sacral chordomas is both safe and effective. Gross total resection, tumour volume <130 ml and tumour restricted to the bone are favourable prognostic factors for local tumour control.

Dose Calculation Algorithms for External Radiation Therapy: An Overview for Practitioners

Radiation therapy (RT) is a constantly evolving therapeutic technique; improvements are continuously being introduced for both methodological and practical aspects. Among the features that have undergone a huge evolution in recent decades, dose calculation algorithms are still rapidly changing. This process is propelled by the awareness that the agreement between the delivered and calculated doses is of paramount relevance in RT, since it could largely affect clinical outcomes. The aim of this work is to provide an overall picture of the main dose calculation algorithms currently used in RT, summarizing their underlying physical models and mathematical bases, and highlighting their strengths and weaknesses, referring to the most recent studies on algorithm comparisons. This handy guide is meant to provide a clear and concise overview of the topic, which will prove useful in helping clinical medical physicists to perform their responsibilities more effectively and efficiently, increasing patient benefits and improving the overall quality of the management of radiation treatment.

Innovations and the Use of Collimators in the Delivery of Pencil Beam Scanning Proton Therapy

Purpose

The development of collimating technologies has become a recent focus in pencil beam scanning (PBS) proton therapy to improve the target conformity and healthy tissue sparing through field-specific or energy-layer–specific collimation. Given the growing popularity of collimators for low-energy treatments, the purpose of this work was to summarize the recent literature that has focused on the efficacy of collimators for PBS and highlight the development of clinical and preclinical collimators.

Materials and Methods

The collimators presented in this work were organized into 3 categories: per-field apertures, multileaf collimators (MLCs), and sliding-bar collimators. For each case, the system design and planning methodologies are summarized and intercompared from their existing literature. Energy-specific collimation is still a new paradigm in PBS and the 2 specific collimators tailored toward PBS are presented including the dynamic collimation system (DCS) and the Mevion Adaptive Aperture.

Results

Collimation during PBS can improve the target conformity and associated healthy tissue and critical structure avoidance. Between energy-specific collimators and static apertures, static apertures have the poorest dose conformity owing to collimating only the largest projection of a target in the beam's eye view but still provide an improvement over uncollimated treatments. While an external collimator increases secondary neutron production, the benefit of collimating the primary beam appears to outweigh the risk. The greatest benefit has been observed for low- energy treatment sites.

Conclusion

The consensus from current literature supports the use of external collimators in PBS under certain conditions, namely low-energy treatments or where the nominal spot size is large. While many recent studies paint a supportive picture, it is also important to understand the limitations of collimation in PBS that are specific to each collimator type. The emergence and paradigm of energy-specific collimation holds many promises for PBS proton therapy.

Commissioning of a clinical pencil beam scanning proton therapy unit for ultra-high dose rates (FLASH)

PURPOSE

The purpose of this work was to provide a flexible platform for FLASH research with protons by adapting a former clinical pencil beam scanning gantry to irradiations with ultra-high dose rates.

METHODS

PSI Gantry 1 treated patients until December 2018. We optimized the beamline parameters to transport the 250 MeV beam extracted from the PSI COMET accelerator to the treatment room, maximizing the transmission of beam intensity to the sample. We characterized a dose monitor on the gantry to ensure good control of the dose, delivered in spot-scanning mode. We characterized the beam for different dose rates and field sizes for transmission irradiations. We explored scanning possibilities in order to enable conformal irradiations or transmission irradiations of large targets (with transverse scanning).

RESULTS

We achieved a transmission of 86% from the cyclotron to the treatment room. We reached a peak dose rate of 9000 Gy/s at 3 mm water equivalent depth, along the central axis of a single pencil beam. Field sizes of up to 5 × 5 mm² were achieved for single-spot FLASH irradiations. Fast transverse scanning allowed to cover a field of 16 × 1.2 cm² . With the use of a nozzle-mounted range shifter, we are able to span depths in water ranging from 19.6 to 37.9 cm. Various dose levels were delivered with precision within less than 1%.

CONCLUSIONS

We have realized a proton FLASH irradiation setup able to investigate continuously a wide dose rate spectrum, from 1 to 9000 Gy/s in single-spot irradiation as well as in the pencil beam scanning mode. As such, we have developed a versatile test bench for FLASH research.

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.

Plan quality effects of maximum monitor unit constraints in pencil beam scanning proton therapy for central nervous system and skull base tumors

PURPOSE/OBJECTIVE(S)

With reports of CNS toxicity in patients treated with proton therapy at doses lower than would be expected based on photon data, it has been proposed that heavy monitor unit (MU) weighting of pencil beam scanning (PBS) proton therapy spots may potentially increase the risk of toxicity. We evaluated the impact of maximum MU weighting per spot (maxMU/spot) restrictions on PBS plan quality, prior to implementing clinic-wide maxMU/spot restrictions.

MATERIALS/METHODS

PBS plans of 11 patients, of which 3 plans included boosts, for a total of 14 PBS sample cases were included. Per sample case, a single dosimetrist created 4 test plans, gradually reducing the maxMU/spot in the plan. Test Plan 1, unrestricted in maxMU/spot, was the reference for all restricted plan comparisons (comparison sets 2 vs. 1; 3 vs. 1; and 4 vs. 1). The impact of MU/spot restrictions on plan quality metrics were analyzed with Wilcoxon signed rank test analyses. Treatment delivery time was modeled for a representative case.

RESULTS

A total of 14 PBS sample cases, 7 (50%) single-field optimized, 7 (50%) multi-field optimized, 9 (64%) delivering > 3500 cGy, 9 (64%) with 3 beams, and 7 (50%) without a range shifter were included. There were no differences in plan quality metrics of target coverage (V95% and V100% prescription), conformality and gradient indices, hot spot volume (V105% prescription), and dose to normal brain (V10%/30%/50%/70%/90%/100% prescription) with reductions of allowable maxMU/spot across all comparison sets (p > 0.05). Max MU/spot restrictions did not increase treatment delivery time when analyzed for a representative case.

CONCLUSION

MaxMU/spot restrictions within the thresholds evaluated in this study did not degrade overall plan quality metrics. Future studies should evaluate spot weighting with linear energy transfer/relative biologic effectiveness-informed planning to determine if spot weighting manipulation impacts clinical outcomes and mitigates toxicity.

FLASH Irradiation with Proton Beams: Beam Characteristics and Their Implications for Beam Diagnostics

FLASH irradiations use dose-rates orders of magnitude higher than commonly used in patient treatments. Such irradiations have shown interesting normal tissue sparing in cell and animal experiments, and, as such, their potential application to clinical practice is being investigated. Clinical accelerators used in proton therapy facilities can potentially provide FLASH beams; therefore, the topic is of high interest in this field. However, a clear FLASH effect has so far been observed in presence of high dose rates (>40 Gy/s), high delivered dose (tens of Gy), and very short irradiation times (<300 ms). Fulfilling these requirements poses a serious challenge to the beam diagnostics system of clinical facilities. We will review the status and proposed solutions, from the point of view of the beam definitions for FLASH and their implications for beam diagnostics. We will devote particular attention to the topics of beam monitoring and control, as well as absolute dose measurements, since finding viable solutions in these two aspects will be of utmost importance to guarantee that the technique can be adopted quickly and safely in clinical practice.

A proton beam energy modulator for rapid proton therapy

We present the design for a rapid proton energy modulator with radiofrequency accelerator cavities, which can deliver the proton radiation dose to varied depth in human tissues much faster than traditional mechanical beam energy degraders. The proton energy modulator is designed as a multi-cell 1-m long accelerator working at 2.856 GHz. Each individual accelerator cavity is powered by a 400 kW compact klystron to provide an accelerating/decelerating gradient of 30 MV/m. The high gradient is enabled by the individual power coupling regime, which provides a high shunt impedance. Beam dynamics simulations were performed, showing that the energy modulator can provide ±30 MeV of beam energy change for a 150 MeV, 7 mm long (full length) proton bunch, and the total energy spread of 3 MeV is satisfactory to clinical needs. A prototype experiment of a single cell has been built and tested, and the low-power microwave measurement results agree very well with simulations. The energy modulator is optimized for the 150 MeV cyclotron proton beam, while this approach can work with different beam energies.

On-line range verification for proton beam therapy using spherical ionoacoustic waves with resonant frequency

In contrast to conventional X-ray therapy, proton beam therapy (PBT) can confine radiation doses to tumours because of the presence of the Bragg peak. However, the precision of the treatment is currently limited by the uncertainty in the beam range. Recently, a unique range verification methodology has been proposed based on simulation studies that exploit spherical ionoacoustic waves with resonant frequency (SPIREs). SPIREs are emitted from spherical gold markers in tumours initially introduced for accurate patient positioning when the proton beam is injected. These waves have a remarkable property: their amplitude is linearly correlated with the residual beam range at the marker position. Here, we present proof-of-principle experiments using short-pulsed proton beams at the clinical dose to demonstrate the feasibility of using SPIREs for beam-range verification with submillimetre accuracy. These results should substantially contribute to reducing the range uncertainty in future PBT applications.

Quantitative analysis of treatments using real-time image gated spot-scanning with synchrotron-based proton beam therapy system log data

A synchrotron-based real-time image gated spot-scanning proton beam therapy (RGPT) system with inserted fiducial markers can irradiate a moving tumor with high accuracy. As gated treatments increase the beam delivery time, this study aimed to investigate the frequency of intra-field adjustments corresponding to the baseline shift or drift and the beam delivery efficiency of a synchrotron-based RGPT system. Data from 118 patients corresponding to 127 treatment plans and 2810 sessions between October 2016 and March 2019 were collected. We quantitatively analyzed the proton beam delivery time, the difference between the ideal beam delivery time based on a simulated synchrotron magnetic excitation pattern and the actual treatment beam delivery time, frequency corresponding to the baseline shift or drift, and the gating efficiency of the synchrotron-based RGPT system according to the proton beam delivery machine log data. The mean actual beam delivery time was 7.1 min, and the simulated beam delivery time in an ideal environment with the same treatment plan was 2.9 min. The average difference between the actual and simulated beam delivery time per session was 4.3 min. The average frequency of intra-field adjustments corresponding to baseline shift or drift and beam delivery efficiency were 21.7% and 61.8%, respectively. Based on our clinical experience with a synchrotron-based RGPT system, we determined the frequency corresponding to baseline shift or drift and the beam delivery efficiency using the beam delivery machine log data. To maintain treatment accuracy within ± 2.0 mm, intra-field adjustments corresponding to baseline shift or drift were required in approximately 20% of cases. Further improvements in beam delivery efficiency may be realized by shortening the beam delivery time.

Characterization of a Low-Cost Plastic Fiber Array Detector for Proton Beam Dosimetry

The Pencil Beam Scanning (PBS) technique in proton therapy uses fast magnets to scan the tumor volume rapidly. Changing the proton energy allows changing to layers in the third dimension, hence scanning the same volume several times. The PBS approach permits adapting the speed and/or current to modulate the delivered dose. We built a simple prototype that measures the dose distribution in a single step. The active detection material consists of a single layer of scintillating fibers (i.e., 1D) with an active length of 100 mm, a width of 18.25 mm, and an insignificant space (20 μm) between them. A commercial CMOS-based camera detects the scintillation light. Short exposure times allow running the camera at high frame rates, thus, monitoring the beam motion. A simple image processing method extracts the dose information from each fiber of the array. The prototype would allow scaling the concept to multiple layers read out by the same camera, such that the costs do not scale with the dimensions of the fiber array. Presented here are the characteristics of the prototype, studied under two modalities: spatial resolution, linearity, and energy dependence, characterized at the Center for Proton Therapy (Paul Scherrer Institute); the dose rate response, measured at an electron accelerator (Swiss Federal Institute of Metrology).

A robust optimization method for weighted-layer-stacking proton beam therapy

The layer-stacking method can provide three-dimensional conformal dose distributions to the target based on a passive scattering method using mini-spread-out Bragg peak (SOBP). The purpose of this work is to demonstrate the effectiveness of a new weight optimization algorithm that can enhance the robustness of dose distributions against layer depth variation in layer-stacking proton beam therapy. In the robustness algorithm, the upper limit of the layer's weight was adapted to the conventional algorithm and varied for 620 weight set evaluations. The optimal weight set was selected by using an analytical objective function based on Gaussian function with σ = 3 mm for WED variation. Then, we evaluated the stabilities of the one-dimensional depth dose distribution against WED variation generated by Gaussian samples. Three-dimensional dose distributions in the water phantom were also evaluated using the Monte-Carlo dose calculation. The variation of dose as well as dose volume histograms for the spherical target and the organ at risk (OAR) were evaluated. The robustness algorithm reduced the change of the dose distribution due to the WED variation by a factor of almost 3/4 compared to those with the conventional procedure. The rate of 91.8% in total samples was maintained within 5% change of the maximum dose, compared with the rate of 64.9% in the conventional algorithm. In the MC calculation, the high dose-volume in the OAR was reduced around the lateral penumbra and distal falloff region by the robustness algorithm. The stability of depth dose distributions was enhanced under the WED variation, compared to the conventional algorithm. This robust algorithm in layer-stacking proton therapy may be useful for treatment in which the sharpness of the distal falloff along the depth distribution needs to be maintained to spare the organ at risk and keep the dose coverage for the target tumor.

The potential of Gantry beamline large momentum acceptance for real time tumour tracking in pencil beam scanning proton therapy

Tumour tracking is an advanced radiotherapy technique for precise treatment of tumours subject to organ motion. In this work, we addressed crucial aspects of dose delivery for its realisation in pencil beam scanning proton therapy, exploring the momentum acceptance and global achromaticity of a Gantry beamline to perform continuous energy regulation with a standard upstream degrader. This novel approach is validated on simulation data from three geometric phantoms of increasing complexity and one liver cancer patient using 4D dose calculations. Results from a standard high-to-low beamline ramping scheme were compared to alternative energy meandering schemes including combinations with rescanning. Target coverage and dose conformity were generally well recovered with tumour tracking even though for particularly small targets, large variations are reported for the different approaches. Meandering in energy while rescanning has a positive impact on target homogeneity and similarly, hot spots outside the targets are mitigated with a relatively fast convergence rate for most tracking scenarios, halving the volume of hot spots after as little as 3 rescans. This work investigates the yet unexplored potential of having a large momentum acceptance in medical beam line, and provides an alternative to take tumour tracking with particle therapy closer to clinical translation.