Simultaneous dose and dose rate optimization (SDDRO) for FLASH proton therapy

Biological optimization for hybrid proton-photon radiotherapy

Objective.Hybrid proton-photon radiotherapy (RT) is a cancer treatment option to broaden access to proton RT. Additionally, with a refined treatment planning method, hybrid RT has the potential to offer superior plan quality compared to proton-only or photon-only RT, particularly in terms of target coverage and sparing organs-at-risk (OARs), when considering robustness to setup and range uncertainties. However, there is a concern regarding the underestimation of the biological effect of protons on OARs, especially those in close proximity to targets. This study seeks to develop a hybrid treatment planning method with biological dose optimization, suitable for clinical implementation on existing proton and photon machines, with each photon or proton treatment fraction delivering a uniform target dose.Approach.The proposed hybrid biological dose optimization method optimized proton and photon plan variables, along with the number of fractions for each modality, minimizing biological dose to the OARs and surrounding normal tissues. To mitigate underestimation of hot biological dose spots, proton biological dose was minimized within a ring structure surrounding the target. Hybrid plans were designed to be deliverable separately and robustly on existing proton and photon machines, with enforced uniform target dose constraints for the proton and photon fraction doses. A probabilistic formulation was utilized for robust optimization of setup and range uncertainties for protons and photons. The nonconvex optimization problem, arising from minimum monitor unit constraint and dose-volume histogram constraints, was solved using an iterative convex relaxation method.Main results.Hybrid planning with biological dose optimization effectively eliminated hot spots of biological dose, particularly in normal tissues surrounding the target, outperforming proton-only planning. It also provided superior overall plan quality and OAR sparing compared to proton-only or photon-only planning strategies.Significance.This study presents a novel hybrid biological treatment planning method capable of generating plans with reduced biological hot spots, superior plan quality to proton-only or photon-only plans, and clinical deliverability on existing proton and photon machines, separately and robustly.

Very high-energy electron therapy as light-particle alternative to transmission proton FLASH therapy - An evaluation of dosimetric performances

PURPOSE

Clinical translation of FLASH-radiotherapy (RT) to deep-seated tumours is still a technological challenge. One proposed solution consists of using ultra-high dose rate transmission proton (TP) beams of about 200-250 MeV to irradiate the tumour with the flat entrance of the proton depth-dose profile. This work evaluates the dosimetric performance of very high-energy electron (VHEE)-based RT (50-250 MeV) as a potential alternative to TP-based RT for the clinical transfer of the FLASH effect.

METHODS

Basic physics characteristics of VHEE and TP beams were compared utilizing Monte Carlo simulations in water. A VHEE-enabled research treatment planning system was used to evaluate the plan quality achievable with VHEE beams of different energies, compared to 250 MeV TP beams for a glioblastoma, an oesophagus, and a prostate cancer case.

RESULTS

Like TP, VHEE above 100 MeV can treat targets with roughly flat (within ± 20 %) depth-dose distributions. The achievable dosimetric target conformity and adjacent organs-at-risk (OAR) sparing is consequently driven for both modalities by their lateral beam penumbrae. Electron beams of 400[500] MeV match the penumbra of 200[250] MeV TP beams and penumbra is increased for lower electron energies. For the investigated patient cases, VHEE plans with energies of 150 MeV and above achieved a dosimetric plan quality comparable to that of 250 MeV TP plans. For the glioblastoma and the oesophagus case, although having a decreased conformity, even 100 MeV VHEE plans provided a similar target coverage and OAR sparing compared to TP.

CONCLUSIONS

VHEE-based FLASH-RT using sufficiently high beam energies may provide a lighter-particle alternative to TP-based FLASH-RT with comparable dosimetric plan quality.

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.

Feasibility study of multiple-energy Bragg peak proton FLASH on a superconducting gantry with large momentum acceptance

BACKGROUND

While the Bragg peak proton beam (BP) is capable of superior target conformity and organs-at-risk sparing than the transmission proton beam (TB), its efficacy in FLASH-RT is hindered by both a slow energy switching process and the beam current. A universal range shifter (URS) can pull back the high-energy proton beam while preserving the beam current. Meanwhile, a superconducting gantry with large momentum acceptance (LMA-SC gantry) enables fast energy switching.

PURPOSE

This study explores the feasibility of multiple-energy BP FLASH-RT on the LMA-SC gantry.

METHOD AND MATERIALS

A simultaneous dose and spot map optimization algorithm was developed for BP FLASH-RT treatment planning to improve the dose delivery efficiency. The URS was designed to be 0-27 cm thick, with 1 cm per step. BP plans using the URS were optimized using single-field optimization (SFO) and multiple-field optimization (MFO) for ten prostate cancer patients and ten lung cancer patients. The plan delivery parameters, dose, and dose rate metrics of BP plans were compared to those of TB plans using the parameters of the LMA-SC gantry.

RESULTS

Compared to TB plans, BP plans significantly reduced MUs by 42.7% (P < 0.001) with SFO and 33.3% (P < 0.001) with MFO for prostate cases. For lung cases, the reduction in MUs was 56.8% (P < 0.001) with SFO and 36.4% (P < 0.001) with MFO. BP plans also outperformed TB plans by reducing mean normal tissue doses. BP-SFO plans achieved a reduction of 56.7% (P < 0.001) for prostate cases and 57.7% (P < 0.001) for lung cases, while BP-MFO plans achieved a reduction of 54.2% (P < 0.001) for the prostate case and 40.0% (P < 0.001) for lung cases. For both TB and BP plans, normal tissues in prostate and lung cases received 100.0% FLASH dose rate coverage (>40 Gy/s).

CONCLUSIONS

By utilizing the URS and the LMA-SC gantry, it is possible to perform multiple-energy BP FLASH-RT, resulting in better normal tissue sparing, as compared to TB plans.

A Novel Inverse Algorithm To Solve the Integrated Optimization of Dose, Dose Rate, and Linear Energy Transfer of Proton FLASH Therapy With Sparse Filters

PURPOSE

The recently proposed Integrated Physical Optimization Intensity Modulated Proton Therapy (IPO-IMPT) framework allows simultaneous optimization of dose, dose rate, and linear energy transfer (LET) for ultra-high dose rate (FLASH) treatment planning. Finding solutions to IPO-IMPT is difficult because of computational intensiveness. Nevertheless, an inverse solution that simultaneously specifies the geometry of a sparse filter and weights of a proton intensity map is desirable for both clinical and preclinical applications. Such solutions can reduce effective biologic dose to organs at risk in patients with cancer as well as reduce the number of animal irradiations needed to derive extra biologic dose models in preclinical studies.

METHODS AND MATERIALS

Unlike the initial forward heuristic, this inverse IPO-IMPT solution includes simultaneous optimization of sparse range compensation, sparse range modulation, and spot intensity. The daunting computational tasks vital to this endeavor were resolved iteratively with a distributed computing framework to enable Simultaneous Intensity and Energy Modulation and Compensation (SIEMAC). SIEMAC was demonstrated on a human patient with central lung cancer and a minipig.

RESULTS

SIEMAC simultaneously improves maps of spot intensities and patient-field-specific sparse range compensators and range modulators. For the patient with lung cancer, at our maximum nozzle current of 300 nA, dose rate coverage above 100 Gy/s increased from 57% to 96% in the lung and from 93% to 100% in the heart, and LET coverage above 4 keV/µm dropped from 68% to 9% in the lung and from 26% to <1% in the heart. For a simple minipig plan, the full-width half-maximum of the dose, dose rate, and LET distributions decreased by 30%, 1.6%, and 57%, respectively, again with similar target dose coverage, thus reducing uncertainty in these quantities for preclinical studies.

CONCLUSIONS

The inverse solution to IPO-IMPT demonstrated the capability to simultaneously modulate subspot proton energy and intensity distributions for clinical and preclinical studies.

Lattice position optimization for LATTICE therapy

BACKGROUND

LATTICE radiation therapy delivers 3D heterogenous dose of high peak-to-valley dose ratio (PVDR) to the tumor target, with peak dose at lattice vertices inside the target and valley dose for the rest of the target. Although the lattice vertex positions can impact PVDR inside the target and sparing of organs-at-risk (OAR), they are fixed as constants and not optimized during treatment planning in current clinical practice.

PURPOSE

This work proposes a new LATTICE plan optimization method that can optimize lattice vertex positions during LATTICE treatment planning, which is the first lattice position optimization study to the best of our knowledge.

METHODS

The new LATTICE treatment planning method optimizes lattice vertex positions as well as other plan variables (e.g., photon fluences or proton spot weights), with optimization objectives for target PVDR and OAR sparing. To satisfy mathematical differentiability, the lattice vertices are approximated in sigmoid functions. For geometric feasibility, proper geometry constraints are enforced onto lattice vertex positions. The lattice position optimization problem is solved by iterative convex relaxation (ICR) method and alternating direction method of multipliers (ADMM), and lattice vertex positions and photon/proton plan variables are jointly updated via the Quasi-Newton method.

RESULTS

Both photon and proton LATTICE RT were considered, and the optimal lattice vertex positions in terms of plan objectives were found by solving all possible combinations on given discrete positions via exhaustive searching based on standard IMRT/IMPT, which served as the ground truth for validating the new LATTICE method. The results show that the new method indeed provided the optimal lattice vertex positions with the smallest optimization objective, the largest target PVDR, and the best OAR sparing.

CONCLUSIONS

A new LATTICE treatment planning method is proposed and validated that can optimize lattice vertex positions as well as other photon or proton plan variables for improving target PVDR and OAR sparing.

Very high-energy electron dose calculation using the Fermi-Eyges theory of multiple scattering and a simplified pencil beam model

BACKGROUND

Very high-energy electrons (VHEE) radiotherapy, in the energy range of 100-200 MeV is currently considered a promising technique for the future of radiation therapy and could benefit from the promises of ultra-high dose rate FLASH therapy. However, to our knowledge, no analytical calculation models have been tested for this type of application and the approximations proposed for multiple scattering with electron beams have not been extensively evaluated at these high energies.

PURPOSE

In this work, we discuss the derivation of a simple and fast algorithm based on the Fermi-Eyges theory of multiple Coulomb scattering for fast dose calculation for VHEE beams (up to 200 MeV). Similar to the Gaussian pencil beam models used for electron or proton beams, this pencil beam kernel is separated into a central and an off-axis term. Monte Carlo simulations are performed to compare the analytical calculations with simulations and to determine the parametrizations used in the model at the highest electron energies.

METHODS

The normalized electron planar fluence distribution is described in water according to the Fermi-Eyges theory of multiple Coulomb scattering and a double Gaussian distribution model. The main quantities used in the model and their calculation (mass angular scattering power, mean energy, range straggling) are discussed and tested for electron energies up to 200 MeV. The TOPAS/Geant4 Monte Carlo (MC) toolkit is used to compare analytical calculations with MC simulations for a theoretical pencil beam irradiation and to find the best parameters describing the range straggling. The model is then tested on a realistic simulation of a pencil beam scanning beamline with treatment field dimensions up to 15 × 15 cm2  and for deep-seated targets.

RESULTS

Radial dose distributions of a pencil beam in water were calculated with the model and compared with the results of a complete Monte Carlo simulation. A good agreement (within 2%/2 mm gamma passing rate superior to 90%, and a mean deviation between calculated and simulated pencil beam radial spread smaller than 0.6 mm) was observed between analytical dose distributions and simulations for energies up to 200 MeV and field sizes up to 15 × 15 cm2 .

CONCLUSIONS

A parameterization of an electron source and an analytical pencil beam model were proposed in this work, thereby allowing a suitable reproduction of the lateral fluence of a VHEE beam and good agreement between calculations and simulated data. Further improvement of the method would require the consideration of a model describing the large-angle scattering of the electrons. The results of this work could support future research into VHEE radiotherapy and might be of interest for use together with VHEE broad beams produced by scanned narrow pencil beams.

A Novel Ultrahigh-Dose-Rate Proton Therapy Technology: Spot-Scanning Proton Arc Therapy + FLASH (SPLASH)

PURPOSE

To take full advantage of FLASH dose rate (40 Gy/s) and high-dose conformity, we introduce a novel optimization and delivery technique, the spot-scanning proton arc therapy (SPArc) + FLASH (SPLASH).

METHODS AND MATERIALS

SPLASH framework was implemented in an open-source proton planning platform (MatRad, Department of Medical Physics in Radiation Oncology, German Cancer Research Center). It optimizes with the clinical dose-volume constraint based on dose distribution and the dose-average dose rate by minimizing the monitor unit constraint on spot weight and accelerator beam current sequentially, enabling the first dynamic arc therapy with voxel-based FLASH dose rate. This new optimization framework minimizes the overall cost function value combined with plan quality and voxel-based dose-rate constraints. Three representative cases (brain, liver, and prostate cancer) were used for testing purposes. Dose-volume histogram, dose-rate-volume histogram, and dose-rate map were compared among intensity modulated proton radiation therapy (IMPT), SPArc, and SPLASH.

RESULTS

SPLASH/SPArc could offer superior plan quality over IMPT in terms of dose conformity. The dose-rate-volume histogram results indicated SPLASH could significantly improve V40 Gy/s in the target and region of interest for all tested cases compared with SPArc and IMPT. The optimal beam current per spot is simultaneously generated, which is within the existing proton machine specifications in the research version (<200 nA).

CONCLUSIONS

SPLASH offers the first voxel-based ultradose-rate and high-dose conformity treatment using proton beam therapy. Such a technique has the potential to fit the needs of a broad range of disease sites and simplify clinical workflow without applying a patient-specific ridge filter, which has never before been demonstrated.

Superiorization of projection algorithms for linearly constrained inverse radiotherapy treatment planning

Objective

We apply the superiorization methodology to the constrained intensity-modulated radiation therapy (IMRT) treatment planning problem. Superiorization combines a feasibility-seeking projection algorithm with objective function reduction: The underlying projection algorithm is perturbed with gradient descent steps to steer the algorithm towards a solution with a lower objective function value compared to one obtained solely through feasibility-seeking.

Approach

Within the open-source inverse planning toolkit matRad, we implement a prototypical algorithmic framework for superiorization using the well-established Agmon, Motzkin, and Schoenberg (AMS) feasibility-seeking projection algorithm and common nonlinear dose optimization objective functions. Based on this prototype, we apply superiorization to intensity-modulated radiation therapy treatment planning and compare it with (i) bare feasibility-seeking (i.e., without any objective function) and (ii) nonlinear constrained optimization using first-order derivatives. For these comparisons, we use the TG119 water phantom, the head-and-neck and the prostate patient of the CORT dataset.

Main results

Bare feasibility-seeking with AMS confirms previous studies, showing it can find solutions that are nearly equivalent to those found by the established piece-wise least-squares optimization approach. The superiorization prototype solved the linearly constrained planning problem with similar dosimetric performance to that of a general-purpose nonlinear constrained optimizer while showing smooth convergence in both constraint proximity and objective function reduction.

Significance

Superiorization is a useful alternative to constrained optimization in radiotherapy inverse treatment planning. Future extensions with other approaches to feasibility-seeking, e.g., with dose-volume constraints and more sophisticated perturbations, may unlock its full potential for high performant inverse treatment planning.

A FLASH model of radiolytic oxygen depletion and reactive oxygen species for differential tumor and normal-tissue response

Objective

FLASH-RT can potentially improve the sparing of normal tissues while preserving the tumoricidal efficiency, owing to the radiation with ultra-high dose rate. However, the FLASH mechanism remains to be solved. A popular FLASH model is based on radiolytic oxygen depletion (ROD), which explains for radiation protection of normal tissues under FLASH-RT. However, ROD does not explain the preservation of tumoricidal efficiency for tumors. This work will develop a ROS+ROD FLASH model that can explain the differential tumor and normal-tissue response.

Approach

The new FLASH model utilizes reactive oxygen species (ROS) in addition to ROD, and takes into account that ROS level decreases during FLASH-RT. Specifically, the differential-equation model takes into account that the basic ROS level is lower during FLASH-RT and the degeneration rates of ROS are different in tumor cells and healthy cells. Based on this ROS+ROD FLASH model, the surviving fractions of tumor and normal cells are respectively compared between conventional radiotherapy (CONV-RT) and FLASH-RT.

Main results

While ROD alone does not distinguish the response of tumors and normal tissues to FLASH-RT, the proposed new FLASH model based on ROD and ROS successfully explained the differential response of tumors and normal tissues to FLASH-RT, i.e., the preserved tumoricidal capability, which cannot be explained by ROD alone, and the extra normal-tissue protection owing to the ultra-high dose rate.

Significance

Since the ROS level decreases slower in tumors than in normal tissues, during FLASH-RT, ROS decreases more in normal tissue, thus can get more protection. By incorporating ROS in addition to ROD, the new FLASH model can not only recover all results by previous FLASH model with ROD alone, but also explain the differential response: preserved lethality of FLASH-RT to tumors and improved protection to normal tissues.

Framework for Quality Assurance of Ultrahigh Dose Rate Clinical Trials Investigating FLASH Effects and Current Technology Gaps

FLASH radiation therapy (FLASH-RT), delivered with ultrahigh dose rate (UHDR), may allow patients to be treated with less normal tissue toxicity for a given tumor dose compared with currently used conventional dose rate. Clinical trials are being carried out and are needed to test whether this improved therapeutic ratio can be achieved clinically. During the clinical trials, quality assurance and credentialing of equipment and participating sites, particularly pertaining to UHDR-specific aspects, will be crucial for the validity of the outcomes of such trials. This report represents an initial framework proposed by the NRG Oncology Center for Innovation in Radiation Oncology FLASH working group on quality assurance of potential UHDR clinical trials and reviews current technology gaps to overcome. An important but separate consideration is the appropriate design of trials to most effectively answer clinical and scientific questions about FLASH. This paper begins with an overview of UHDR RT delivery methods. UHDR beam delivery parameters are then covered, with a focus on electron and proton modalities. The definition and control of safe UHDR beam delivery and current and needed dosimetry technologies are reviewed and discussed. System and site credentialing for large, multi-institution trials are reviewed. Quality assurance is then discussed, and new requirements are presented for treatment system standard analysis, patient positioning, and treatment planning. The tables and figures in this paper are meant to serve as reference points as we move toward FLASH-RT clinical trial performance. Some major questions regarding FLASH-RT are discussed, and next steps in this field are proposed. FLASH-RT has potential but is associated with significant risks and complexities. We need to redefine optimization to focus not only on the dose but also on the dose rate in a manner that is robust and understandable and that can be prescribed, validated, and confirmed in real time. Robust patient safety systems and access to treatment data will be critical as FLASH-RT moves into the clinical trials.

An orthogonal matching pursuit optimization method for solving minimum-monitor-unit problems: Applications to proton IMPT, ARC and FLASH

BACKGROUND

The intensities (i.e., number of protons in monitor unit [MU]) of deliverable proton spots need to be either zero or meet a minimum-MU (MMU) threshold, which is a nonconvex problem. Since the dose rate is proportionally associated with the MMU threshold, higher-dose-rate proton radiation therapy (RT) (e.g., efficient intensity modulated proton therapy (IMPT) and ARC proton therapy, and high-dose-rate-induced FLASH effect needs to solve the MMU problem with larger MMU threshold, which however makes the nonconvex problem more difficult to solve.

PURPOSE

This work will develop a more effective optimization method based on orthogonal matching pursuit (OMP) for solving the MMU problem with large MMU thresholds, compared to state-of-the-art methods, such as alternating direction method of multipliers (ADMM), proximal gradient descent method (PGD), or stochastic coordinate descent method (SCD).

METHODS

The new method consists of two essential components. First, the iterative convex relaxation (ICR) method is used to determine the active sets for dose-volume planning constraints and decouple the MMU constraint from the rest. Second, a modified OMP optimization algorithm is used to handle the MMU constraint: the non-zero spots are greedily selected via OMP to form the solution set to be optimized, and then a convex constrained subproblem is formed and can be conveniently solved to optimize the spot weights restricted to this solution set via OMP. During this iterative process, the new non-zero spots localized via OMP will be adaptively added to or removed from the optimization objective.

RESULTS

The new method via OMP is validated in comparison with ADMM, PGD and SCD for high-dose-rate IMPT, ARC, and FLASH problems of large MMU thresholds, and the results suggest that OMP substantially improved the plan quality from PGD, ADMM and SCD in terms of both target dose conformality (e.g., quantified by max target dose and conformity index) and normal tissue sparing (e.g., mean and max dose). For example, in the brain case, the max target dose for IMPT/ARC/FLASH was 368.0%/358.3%/283.4% respectively for PGD, 154.4%/179.8%/150.0% for ADMM, 134.5%/130.4%/123.0% for SCD, while it was <120% in all scenarios for OMP; compared to PGD/ADMM/SCD, OMP improved the conformity index from 0.42/0.52/0.33 to 0.65 for IMPT and 0.46/0.60/0.61 to 0.83 for ARC.

CONCLUSIONS

A new OMP-based optimization algorithm is developed to solve the MMU problems with large MMU thresholds, and validated using examples of IMPT, ARC, and FLASH with substantially improved plan quality from ADMM, PGD, and SCD.

Feasibility study of hybrid inverse planning with transmission beams and single-energy spread-out Bragg peaks for proton FLASH radiotherapy

BACKGROUND

Ultra-high dose rate (FLASH) proton planning with only transmission beams (TBs) has limitations in normal tissue sparing. The single-energy spread-out Bragg peaks (SESOBPs) of the FLASH dose rate have been demonstrated feasible for proton FLASH planning.

PURPOSE

To investigate the feasibility of combining TBs and SESOBPs for proton FLASH treatment.

METHODS

A hybrid inverse optimization method was developed to combine the TBs and SESOBPs (TB-SESOBP) for FLASH planning. The SESOBPs were generated field-by-field from spreading out the BPs by pre-designed general bar ridge filters (RFs) and placed at the central target by range shifters (RSs) to obtain a uniform dose within the target. The SESOBPs and TBs were fully placed field-by-field allowing automatic spot selection and weighting in the optimization process. A spot reduction strategy was conducted in the optimization process to push up the minimum MU/spot assuring the plan deliverability at beam current of 165 nA. The TB-SESOBP plans were validated in comparison with the TB only (TB-only) plans and the plans with the combination of TBs and BPs (TB-BP plans) regarding 3D dose and dose rate (dose-averaged dose rate) distributions for five lung cases. The FLASH dose rate coverage (V40Gy/s ) was evaluated in the structure volume receiving > 10% of the prescription dose.

RESULTS

Compared to the TB-only plans, the mean spinal cord D1.2cc drastically reduced by 41% (P < 0.05), the mean lung V7Gy and V7.4 Gy moderately reduced by up to 17% (P < 0.05), and the target dose homogeneity slightly increased in the TB-SESOBP plans. Comparable dose homogeneity was achieved in both TB-SESOBP and TB-BP plans. Besides, prominent improvements were achieved in lung sparing for the cases of relatively large targets by the TB-SESOBP plans compared to the TB-BP plans. The targets and the skin were fully covered with the FLASH dose rate in all three plans. For the OARs, V40Gy/s  = 100% was achieved by the TB-only plans while V40Gy/s  > 85% was obtained by the other two plans.

CONCLUSION

We have demonstrated that the hybrid TB-SESOBP planning was feasible to achieve FLASH dose rate for proton therapy. With pre-designed general bar RFs, the hybrid TB-SESOBP planning could be implemented for proton adaptive FLASH radiotherapy. As an alternative FLASH planning approach to TB-only planning, the hybrid TB-SESOBP planning has great potential in dosimetrically improving OAR sparing while maintaining high target dose homogeneity.

Beam angle optimization for proton therapy via group-sparsity based angle generation method

BACKGROUND

In treatment planning, beam angle optimization (BAO) refers to the selection of a subset with a given number of beam angles from all available angles that provides the best plan quality. BAO is a NP-hard combinatorial problem. Although exhaustive search (ES) can exactly solve BAO by exploring all possible combinations, ES is very time-consuming and practically infeasible.

PURPOSE

To the best of our knowledge, (1) no optimization method has been demonstrated that can provide the exact solution to BAO, and (2) no study has validated an optimization method for solving BAO by benchmarking with the optimal BAO solution (e.g., via ES), both of which will be addressed by this work.

METHODS

This work considers BAO for proton therapy, for example, the selection of 2-4 beam angles for IMPT. The optimal BAO solution is obtained via ES and serves as the ground truth. A new BAO algorithm, namely angle generation (AG) method, is proposed, and demonstrated to provide nearly-exact solutions for BAO in reference to the ES solution. AG iteratively optimizes the angular set via group-sparsity (GS) regularization, until the planning objective does not decrease further.

RESULTS

Since GS alone can also solve BAO, AG was validated and compared with GS for 2-angle brain, 3-angle lung, and 4-angle brain cases, in reference to the optimal BAO solutions obtained by ES: the AG solution had the rank (1/276, 1/2024, 4/10 626), while the GS solution had the rank (42/276, 279/2024, 4328/10 626).

CONCLUSIONS

A new BAO algorithm called AG is proposed and shown to provide substantially improved accuracy for BAO from current methods with nearly-exact solutions to BAO, in reference to the ground truth of optimal BAO solution via ES.

A treatment plan optimization method with direct minimization of number of energy jumps for proton arc therapy

Objective. The optimization of energy layer distributions is crucial to proton arc therapy: on one hand, a sufficient number of energy layers is needed to ensure the plan quality; on the other hand, an excess number of energy jumps (NEJ) can substantially slow down the treatment delivery. This work will develop a new treatment plan optimization method with direct minimization of (NEJ), which will be shown to outperform state-of-the-art methods in both plan quality and delivery efficiency.Approach. The proposed method jointly optimizes the plan quality and minimizes the NEJ. To minimize NEJ, (1) the proton spotsxis summed per energy layer to form the energy vectory; (2)yis binarized via sigmoid transform intoy1; (3)y1is multiplied with a predefined energy order vector via dot product intoy2; (4)y2is filtered through the finite-differencing kernel intoy3in order to identify NEJ; (5) only the NEJ ofy3is penalized, whilexis optimized for plan quality. The solution algorithm to this new method is based on iterative convex relaxation.Main results. The new method is validated in comparison with state-of-the-art methods called energy sequencing (ES) method and energy matrix (EM) method. In terms of delivery efficiency, the new method had fewer NEJ, less energy switching time, and generally less total delivery time. In terms of plan quality, the new method had smaller optimization objective values, lower normal tissue dose, and generally better target coverage.Significance. We have developed a new treatment plan optimization method with direct minimization of NEJ, and demonstrated that this new method outperformed state-of-the-art methods (ES and EM) in both plan quality and delivery efficiency.

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.

Treatment planning consideration for very high-energy electron FLASH radiotherapy

PURPOSE

Very high-energy electron (VHEE) can make up the insufficient treatment depth of the low-energy electron while offering an intermediate dosimetric advantage between photon and proton. Combining FLASH with VHEE, a quantitative comparison between different energies was made, with regard to plan quality, dose rate distribution (both in PTV and OAR), and total duration of treatment (beam-on time).

METHODS

In two patient cases (head and lung), we created the treatment plans utilizing the scanning pencil beam via the Monte Carlo simulation and a PTV-based optimization algorithm. Geant4 was used to simulate VHEE pencil beams and sizes of 0.3-5 mm defined by the full width at half maximum (FWHM). Monoenergetic beams with Gaussian distribution in x and y directions (ISOURC = 19) were used as the source of electrons. A large-scale non-linear solver (IPOPT) was used to calculate the optimal spot weights. After optimization, a quantitative comparison between different energies was made regarding treatment plan quality, dose rate distribution (both in PTV and OAR), and total beam duration.

RESULTS

For head (80 MeV, 100 MeV, and 120 MeV) and lung cases (100 MeV, 120 MeV, and 140 MeV), the minimum beam intensity needs to be ∼2.5 × 10¹¹ electrons/s and ∼9.375 × 10¹¹ electrons/s to allow > 90 % volume of PTV reaching the average dose rate (DADR) higher than 40 Gy/s. At this beam intensity (fraction dose: 10 Gy), the overall irradiation time for the head case is 5258.75 ms (80 MeV), 5149.75 ms (100 MeV), and 4976.75 ms (120 MeV), including scanning time 872.75 ms. For lung cases, this number is 1034.25 ms (100 MeV), 981.55 ms (120 MeV), and 928.15 ms (140 MeV), including scanning time 298.75 ms. The plan of higher energy always performs with a higher dose rate (both in PTV and OAR) and thereby costs less delivery time (beam-on time).

CONCLUSION

The study systematically investigated the currently known FLASH parameters for VHEE radiotherapy and successfully established a benchmark reference for its FLASH dose rate performance.

TVL1-IMPT: Optimization of Peak-to-Valley Dose Ratio Via Joint Total-Variation and L1 Dose Regularization for Spatially Fractionated Pencil-Beam-Scanning Proton Therapy

PURPOSE

Proton minibeam radiation therapy (pMBRT) is a novel proton modality of spatially fractionated RT. pMBRT can reduce the radiation damage to normal tissues via biological dose sparing of high peak-to-valley dose ratio (PVDR). This work will develop a new spatially fractionated IMPT treatment planning method for pMBRT that jointly optimizes the plan quality and maximizes the PVDR.

METHODS

The new optimization method simultaneously maximizes the normal-tissue PVDR and optimizes the dose distribution at tumor targets and organs at risk. The PVDR maximization is through the joint total variation (TV) and L1 regularization with respect to the normal-tissue dose. That is, the beam-eye view projects dose slices of several depths for each beam angle; the TV of dose is maximized, corresponding to the PVDR maximization; and the L1 of dose is minimized, corresponding to the minimization of the organs-at-risk dose and maximization of survival fraction (SF).

RESULTS

The new IMPT method with TV and L1 regularization was validated in comparison with the conventional IMPT method for pMBRT in several clinical cases. The results show that TVL1 provided larger PVDR and SF than the conventional IMPT method for biological sparing of normal tissues, with preserved plan quality in terms of physical dose distribution.

CONCLUSIONS

A new spatially fractionated IMPT treatment planning method was developed for pMBRT that can optimize and improve normal-tissue PVDR and SF by incorporating TV and L1 dose regularization with properly chosen regularization parameters into IMPT.

Pencil-beam Delivery Pattern Optimization Increases Dose Rate for Stereotactic FLASH Proton Therapy

PURPOSE

FLASH dose rates >40 Gy/s are readily available in proton therapy (PT) with cyclotron-accelerated beams and pencil-beam scanning (PBS). The PBS delivery pattern will affect the local dose rate, as quantified by the PBS dose rate (PBS-DR), and therefore needs to be accounted for in FLASH-PT with PBS, but it is not yet clear how. Our aim was to optimize patient-specific scan patterns for stereotactic FLASH-PT of early-stage lung cancer and lung metastases, maximizing the volume irradiated with PBS-DR >40 Gy/s of the organs at risk voxels irradiated to >8 Gy (FLASH coverage).

METHODS AND MATERIALS

Plans to 54 Gy/3 fractions with 3 equiangular coplanar 244 MeV proton shoot-through transmission beams for 20 patients were optimized with in-house developed software. Planning target volume-based planning with a 5 mm margin was used. Planning target volume ranged from 4.4 to 84 cc. Scan-pattern optimization was performed with a Genetic Algorithm, run in parallel for 20 independent populations (islands). Mapped crossover, inversion, swap, and shift operators were applied to achieve (local) optimality on each island, with migration between them for global optimality. The cost function was chosen to maximize the FLASH coverage per beam at >8 Gy, >40 Gy/s, and 40 nA beam current. The optimized patterns were evaluated on FLASH coverage, PBS-DR distribution, and population PBS-DR-volume histograms, compared with standard line-by-line scanning. Robustness against beam current variation was investigated.

RESULTS

The optimized patterns have a snowflake-like structure, combined with outward swirling for larger targets. A population median FLASH coverage of 29.0% was obtained for optimized patterns compared with 6.9% for standard patterns, illustrating a significant increase in FLASH coverage for optimized patterns. For beam current variations of 5 nA, FLASH coverage varied between -6.1%-point and 2.2%-point for optimized patterns.

CONCLUSIONS

Significant improvements on the PBS-DR and, hence, on FLASH coverage and potential healthy-tissue sparing are obtained by sequential scan-pattern optimization. The optimizer is flexible and may be further fine-tuned, based on the exact conditions for FLASH.

An iterative convex relaxation method for proton LET optimization

Objective:A constant relative biological effectiveness of 1.1 in current clinical practice of proton radiotherapy (RT) is a crude approximation and may severely underestimate the biological dose from proton RT to normal tissues, especially near the treatment target at the end of Bragg peaks that exhibits high linear energy transfer (LET). LET optimization can account for biological effectiveness of protons during treatment planning, for minimizing biological proton dose and hot spots to normal tissues. However, the LET optimization is usually nonlinear and nonconvex to solve, for which this work will develop an effective optimization method based on iterative convex relaxation (ICR).Approach: In contrast to the generic nonlinear optimization method, such as Quasi-Newton (QN) method, that does not account for specific characteristics of LET optimization, ICR is tailored to LET modeling and optimization in order to effectively and efficiently solve the LET problem. Specifically, nonlinear dose-averaged LET term is iteratively linearized and becomes convex during ICR, while nonconvex dose-volume constraint and minimum-monitor-unit constraint are also handled by ICR, so that the solution for LET optimization is obtained by solving a sequence of convex and linearized convex subproblems. Since the high LET mostly occurs near the target, a 1 cm normal-tissue expansion of clinical target volume (CTV) (excluding CTV), i.e. CTV1cm, is defined to as an auxiliary structure during treatment planning, where LET is minimized.Main results: ICR was validated in comparison with QN for abdomen, lung, and head-and-neck cases. ICR was effective for LET optimization, as ICR substantially reduced the LET and biological dose in CTV1cm the ring, with preserved dose conformality to CTV. Compared to QN, ICR had smaller LET, physical and biological dose in CTV1cm, and higher conformity index values; ICR was also computationally more efficient, which was about 3 times faster than QN.Significance: A LET-specific optimization method based on ICR has been developed for solving proton LET optimization, which has been shown to be more computationally efficient than generic nonlinear optimizer via QN, with better plan quality in terms of LET, biological and physical dose conformality.

Fraction optimization for hybrid proton-photon treatment planning

BACKGROUND

Hybrid proton-photon radiotherapy (RT) can provide better plan quality than proton or photon only RT, in terms of robustness of target coverage and sparing of organs-at-risk (OAR).

PURPOSE

This work develops a hybrid treatment planning method that can optimize the number of proton and photon fractions as well as proton and photon plan variables, so that the hybrid plans can be clinically delivered day-to-day using either proton or photon machine.

METHODS

In the new hybrid treatment planning method, the total dose distribution (sum of proton dose and photon dose) is optimized for robust target coverage and optimal OAR sparing, by jointly optimizing proton spots and photon fluences, while the target dose uniformity is also enforced individually on both proton dose and photon dose, so that the hybrid plans can be separately and robustly delivered on proton or photon machine. To ensure the target dose uniformity for proton and photon plans, the number of proton and photon fractions is explicitly modeled and optimized, so that the target dose for proton and photon dose components is constrained to be a constant fraction of the total prescription dose while the plan quality based on total dose is optimized. The feasibility of hybrid planning using the proposed method is validated with representative clinical cases including abdomen, lung, head-and-neck (HN), and brain.

RESULTS

For all cases, hybrid plans provided better overall plan quality and OAR sparing than proton-only or photon-only plans, better target dose uniformity and robustness than proton-only plans, quantified by treatment planning objectives and dosimetric parameters. Moreover, for HN and brain cases, hybrid plans also had better target coverage than photon-only plans.

CONCLUSIONS

We have developed a new hybrid treatment planning method that optimizes number of proton and photon fractions as well as proton spots and photon fluences, for generating hybrid plans that can be separately and robustly delivered on proton or photon machines. Preliminary results have demonstrated that hybrid plans generated by the new method have better plan quality than proton-only or photon-only plans.

An Integrated Physical Optimization Framework for Proton Stereotactic Body Radiation Therapy FLASH Treatment Planning Allows Dose, Dose Rate, and Linear Energy Transfer Optimization Using Patient-Specific Ridge Filters

PURPOSE

Patient-specific ridge filters provide a passive means to modulate proton energy to obtain a conformal dose. Here we describe a new framework for optimization of filter design and spot maps to meet the unique demands of ultrahigh-dose-rate (FLASH) radiation therapy. We demonstrate an integrated physical optimization Intensity-modulated proton therapy (IMPT) (IPO-IMPT) approach for optimization of dose, dose-averaged dose rate (DADR), and dose-averaged linear energy transfer (LET_d).

METHODS AND MATERIALS

We developed an inverse planning software to design patient-specific ridge filters that spread the Bragg peak from a fixed-energy, 250-MeV beam to a proximal beam-specific planning target volume. The software defines patient-specific ridge filter pin shapes and uses a Monte Carlo calculation engine, based on Geant4, to provide dose and LET influence matrices. Plan optimization, using matRAD, accommodates the IPO-IMPT objective function considering dose, dose rate, and LET simultaneously with minimum monitor unit constraints. The framework enables design of both regularly spaced and sparse-optimized ridge filters, from which some pins are omitted to allow faster delivery and selective LET optimization. To demonstrate the framework, we designed ridge filters for 3 example patients with lung cancer and optimized the plans using IPO-IMPT.

RESULTS

The IPO-IMPT framework selectively spared the organs at risk by reducing LET and increasing dose rate, relative to IMPT planning. Sparse-optimized ridge filters were superior to regularly spaced ridge filters in dose rate. Depending on which parameter is prioritized, volume distributions and histograms for dose, DADR, and LET_d, using evaluation structures specific to heart, lung, and esophagus, show high levels of FLASH dose-rate coverage and/or reduced LET_d, while maintaining dose coverage within the beam specific planning target volume.

CONCLUSIONS

This proof-of-concept study demonstrates the feasibility of using an IPO-IMPT framework to accomplish proton FLASH stereotactic body proton therapy, accounting for dose, DADR, and LET_d simultaneously.

Dose and dose rate objectives in Bragg peak and shoot-through beam orientation optimization for FLASH proton therapy

PURPOSE

The combined use of Bragg peak (BP) and shoot-through (ST) beams has previously been shown to increase the normal tissue volume receiving FLASH dose rates while maintaining dose conformality compared to conventional intensity-modulated proton therapy (IMPT) methods. However, the fixed beam optimization method has not considered the effects of beam orientation on the dose and dose rates. To maximize the proton FLASH effect, here, we incorporate dose rate objectives into our beam orientation optimization framework.

METHODS

From our previously developed group-sparsity dose objectives, we add upper and lower dose rate terms using a surrogate dose-averaged dose rate definition and solve using the fast-iterative shrinking threshold algorithm. We compare the dosimetry for three head-and-neck cases between four techniques: (1) spread-out BP IMPT (BP), (2) dose rate optimization using BP beams only (BP-DR), (3) dose rate optimization using ST beams only (ST-DR), and (4) dose rate optimization using combined BP and ST (BPST-DR), with the goal of sparing organs at risk without loss of tumor coverage and maintaining high dose rate within a 10 mm region of interest (ROI) surrounding the clinical target volume (CTV).

RESULTS

For BP, BP-DR, ST-DR, and BPST-DR, CTV homogeneity index and Dmax were found to be on average 0.886, 0.867, 0.687, and 0.936 and 107%, 109%, 135%, and 101% of prescription, respectively. Although ST-DR plans were not able to meet dosimetric standards, BPST-DR was able to match or improve either maximum or mean dose in the right submandibular gland, left and right parotids, constrictors, larynx, and spinal cord compared to BP plans. Volume of ROIs receiving greater than 40 Gy/s ( V γ 0 ) ${V_{\gamma 0}})$ was 51.0%, 91.4%, 95.5%, and 92.1% on average.

CONCLUSIONS

The dose rate techniques, particularly BPST-DR, were able to significantly increase dose rate without compromising physical dose compared with BP. Our algorithm efficiently selects beams that are optimal for both dose and dose rate.

Trade-off in healthy tissue sparing of FLASH and fractionation in stereotactic proton therapy of lung lesions with transmission beams

PURPOSE AND OBJECTIVE

Besides a dose-rate threshold of 40-100 Gy/s, the FLASH effect may require a dose >3.5-7 Gy. Even in hypofractioned treatments, with all beams delivered in each fraction (ABEF), most healthy tissue is irradiated to a lower fraction dose. This can be circumvented by single-beam-per-fraction (SBPF) delivery, with a loss of healthy tissue sparing by fractionation. We investigated the trade-off between FLASH and loss of fractionation in SBPF stereotactic proton therapy of lung cancer and determined break-even FLASH-enhancement ratios (FERs).

MATERIALS AND METHODS

Treatment plans for 12 patients were generated. GTV delineations were available and a 5 mm GTV-PTV margin was applied. Equiangular arrangements of 3, 5, 7, and 9 244 MeV proton transmission beams were used. To facilitate SBPF, the number of fractions was equal to the number of beams. Iso-effective fractionation schedules with a single field uniform dose prescription were used: D_95%,PTV = 100%D_pres per beam. All plans were evaluated in terms of dose to lung and conformity of dose to target of FLASH-enhanced biologically equivalent dose (EQD2).

RESULTS

Compared to ABEF, SBPF resulted in a median increase of EQD2_mean to healthy lung of 56%, 58%, 55% and 54% in plans with 3, 5, 7 and 9 fractions respectively and of 90%, 108%, 106% and 102% in V_{100% EQD2}, quantifying conformity. This can be compensated for by FERs of at least 1.28, 1.32, 1.30 and 1.23 respectively for EQD2_mean and 1.29, 1.18, 1.28 and 1.15 for V_100%,EQD2.

CONCLUSION

A FLASH effect outweighing the loss of fractionation in SBPF may be achieved in stereotactic lung treatments. The trade-off with fractionation depends on the conditions under which the FLASH effect occurs. Better understanding of the underlying biology and the impact of delivery conditions is needed.

Dose rate and dose robustness for proton transmission FLASH-RT treatment in lung cancer

Purposes

To evaluate the plan quality and robustness of both dose and dose rate of proton pencil beam scanning (PBS) transmission FLASH delivery in lung cancer treatment.

Methods and materials

An in-house FLASH planning platform was used to optimize 10 lung cancer patients previously consecutively treated with proton stereotactic body radiation therapy (SBRT) to receive 3 and 5 transmission beams (Trx-3fds and Trx-5fds, respectively) to 34 Gy in a single fraction. Perturbation scenarios (n=12) for setup and range uncertainties (5 mm and 3.5%) were introduced, and dose-volume histogram and dose-rate-volume histogram bands were generated. Conventional proton SBRT clinical plans were used as a reference. RTOG 0915 dose metrics and 40 Gy/s dose rate coverage (V_40Gy/s) were used to assess the dose and dose rate robustness.

Results

Trx-5fds yields a comparable iCTV D_2% of 105.3%, whereas Trx-3fds resulted in inferior D_2% of 111.9% to the clinical SBRT plans with D_2% of 105.6% (p<0.05). Both Trx-5fds and Trx-3fds plans had slightly worse dose metrics to organs at risk than SBRT plans. Trx-5fds achieved superior dosimetry robustness for iCTV, esophagus, and spinal cord doses than both Trx-3fds and conventional SBRT plans. There was no significant difference in dose rate robustness for V_40Gy/s coverage between Trx-3fds and Trx-5fds. Dose rate distribution has similar distributions to the dose when perturbation exists.

Conclusion

Transmission plans yield overall modestly inferior plan quality compared to the conventional proton SBRT plans but provide improved robustness and the potential for a toxicity-sparing FLASH effect. By using more beams (5- versus 3-field), both dose and dose rate robustness for transmission plans can be achieved.

FLASH radiotherapy treatment planning and models for electron beams

The FLASH effect designates normal tissue sparing at ultra-high dose rate (UHDR, >40 Gy/s) compared to conventional dose rate (∼0.1 Gy/s) irradiation while maintaining tumour control and has the potential to improve the therapeutic ratio of radiotherapy (RT). UHDR high-energy electron (HEE, 4-20 MeV) beams are currently a mainstay for investigating the clinical potential of FLASH RT for superficial tumours. In the future very-high energy electron (VHEE, 50-250 MeV) UHDR beams may be used to treat deep-seated tumours. UHDR HEE treatment planning focused at its initial stage on accurate dosimetric modelling of converted and dedicated UHDR electron RT devices for the clinical transfer of FLASH RT. VHEE treatment planning demonstrated promising dosimetric performance compared to clinical photon RT techniques in silico and was used to evaluate and optimise the design of novel VHEE RT devices. Multiple metrics and models have been proposed for a quantitative description of the FLASH effect in treatment planning, but an improved experimental characterization and understanding of the FLASH effect is needed to allow for an accurate and validated modelling of the effect in treatment planning. The importance of treatment planning for electron FLASH RT will augment as the field moves forward to treat more complex clinical indications and target sites. In this review, TPS developments in HEE and VHEE are presented considering beam models, characteristics, and future FLASH applications.

Advanced pencil beam scanning Bragg peak FLASH-RT delivery technique can enhance lung cancer planning treatment outcomes compared to conventional multiple-energy proton PBS techniques

PURPOSE

To investigate the dosimetric characteristics between an advanced proton pencil beam scanning (PBS) Bragg peak FLASH technique and conventional PBS planning technique in lung tumors. To evaluate the "FLASHness" of single-field in a multiple-field delivery scheme for a hypofractionation regimen and move a step forward to clinical application.

METHODS

Single-energy PBS Bragg peak FLASH treatment plans were optimized based on a novel Bragg peak tracking technique to enable Bragg peaks to stop at the distal edge of the target. Inverse treatment planning using multiple-field optimization (MFO) can achieve sufficient FLASH dose rate and intensity-modulated proton therapy (IMPT)-equivalent dosimetric quality. The dose rate of organs-at-risk (OARs) and the target were calculated under FLASH machine parameters. A group of 10 consecutive lung SBRT patients was optimized to 34 Gy/fraction using a standard treatment of PBS technique with multiple energy layers as references to the Bragg peak plans. The dosimetric quality was compared between Bragg peak FLASH and conventional plans based on RTOG0915 dose metrics. FLASH dose rate ratios (V_40Gy/s) were calculated as a metric of the FLASH-sparing effect.

RESULTS

For higher dose thresholds, the Bragg peak plans achieved greater V_40Gy/s FLASH coverage for all major OARs. The V40Gy/s was close to 100% for all OARs when the dose thresholds were > 5 Gy for full plan and single beam evaluations. The less "FLASHness" region was associated with a low dose distribution, mainly occurring in the PBS field penumbra region. The conventional IMPT treatment plans yielded slightly superior target dose uniformity with a D_2%(%) of 108.02% versus that of Bragg peak 300 MU plans of 111.81% (p < 0.01) and that of Bragg peak 1200 MU plans of 115.95% (p < 0.01). No significant difference in dose metrics was found between Bragg peak and IMPT treatment plans for the spinal cord, esophagus, heart, or lung-GTV (all p > 0.05).

CONCLUSION

Hypofractionated lung Bragg peak plans can maintain comparable plan quality to conventional PBS while achieving sufficient FLASH dose rate coverage for major OARs for each field under the multiple-field delivery scheme. The novel Bragg peak FLASH technique has the potential to enhance lung cancer planning treatment outcomes compared to standard PBS treatment techniques.

Treatment planning considerations for the development of FLASH proton therapy

With increasing focus on the translation of the observed FLASH effect into clinical practice, this paper presents treatment planning considerations for its development using proton therapy. Potential requirements to induce a FLASH effect are discussed along with the properties of existing proton therapy delivery systems and the changes in planning and delivery approaches required to satisfy these prerequisites. For the exploration of treatment planning approaches for FLASH, developments in treatment planning systems are needed. Flexibility in adapting to new information will be important in such an evolving area. Variations in definitions, threshold values and assumptions can make it difficult to compare different published studies and to interpret previous studies in the context of new information. Together with the fact that much is left to be understood about the underlying mechanism behind the FLASH effect, a systematic and comprehensive approach to information storage is encouraged. Collecting and retaining more detailed information on planned and realised dose delivery as well as reporting the assumptions made in planning studies creates the potential for research to be revisited and re-evaluated in the light of future improvements in understanding. Forward thinking at the time of study development can help facilitate retrospective analysis. This, we hope, will increase the available evidence and accelerate the translation of the FLASH effect into clinical benefit.

Use of single-energy proton pencil beam scanning Bragg peak for intensity-modulated proton therapy FLASH treatment planning in liver hypofractionated radiation therapy

PURPOSE

The transmission proton FLASH technique delivers high doses to the normal tissue distal to the target, which is less conformal compared to the Bragg peak technique. To investigate FLASH RT planning using single-energy Bragg peak beams with a similar beam arrangement as clinical intensity-modulated proton therapy (IMPT) in liver stereotactic body radiation therapy (SBRT) and to characterize the plan quality, dose sparing of organs-at-risk (OARs), and FLASH dose rate percentage.

MATERIALS AND METHODS

An in-house platform was developed to enable inverse IMPT-FLASH planning using single-energy Bragg peaks. A universal range shifter and range compensators were utilized to effectively align the Bragg peak to the distal edge of the target. Two different minimum MU settings of 400 and 800 MU/spot (Bragg-400MU and Bragg-800MU) plans were investigated on 10 consecutive hepatocellular carcinoma patients previously treated by IMPT-SBRT to evaluate the FLASH dose and dose rate coverage for OARs. The IMPT-FLASH using single-energy Bragg peaks delivered 50 Gy in 5 fractions with similar or identical beam arrangement to the clinical IMPT-SBRT plans. NRG GI003 dose constraint metrics were used. Three dose rate calculation methods, including average dose rate (ADR), dose threshold dose rate (DTDR), and dose-averaged dose rate (DADR), were all studied.

RESULTS

The novel spot map optimization can fulfill the inverse planning using single-energy Bragg peaks. All the Bragg peak FLASH plans achieved similar results for the liver-GTV D_mean and heart D_0.5cc , compared to SBRT-IMPT. The Bragg-800MU plans resulted in 18.3% higher CTV D_2cc compared with SBRT (p < 0.05), and no significant difference was found between Bragg-400MU and SBRT plans. For the CTV D_max , SBRT plans resulted in 10.3% (p<0.01) less than Bragg-400MU plans and 16.6% (p<0.01) less than Bragg-800MU plans. The Bragg-800MU plans generally achieved higher ADR, DADR, and DTDR dose rates than Bragg-400MU plans, and DADR mostly led to the highest V_40Gy/s compared to other dose rate calculation methods, whereas ADR led to the lowest. The lower dose rate portions in certain OARs are related to the lower dose deposited due to the farther distances from targets, especially in the penumbra of the beams.

CONCLUSION

Single-energy Bragg peak IMPT-FLASH plans eliminate the exit dose in normal tissues, maintaining comparable dose metrics to the conventional IMPT-SBRT plans while achieving a sufficient FLASH dose rate for liver cancers. This study demonstrates the feasibility of and sufficiently high dose rate when applying Bragg peak FLASH treatment for liver cancer hypofractionated FLASH therapy. The advancement of this novel method has the potential to optimize treatment for liver cancer patients. This article is protected by copyright. All rights reserved.

Energy layer optimization via energy matrix regularization for proton spot-scanning arc therapy

PURPOSE

Spot-scanning arc therapy (SPArc) is an emerging proton modality that can potentially offer a combination of advantages in plan quality and delivery efficiency, compared to traditional IMPT of a few beam angles. Unlike IMPT, frequent low-to-high energy layer switching (so called switch-up (SU)) can degrade delivery efficiency for SPArc. However, it is a tradeoff between the minimization of SU times and the optimization of plan quality. This work will consider the energy layer optimization (ELO) problem for SPArc and develop a new ELO method via energy matrix (EM) regularization to improve plan quality and delivery efficiency.

METHODS

The major innovation of EM method for ELO is to design an energy matrix that encourages desirable energy-layer map with minimal SU during SPArc, and then incorporate this energy matrix into the SPArc treatment planning to simultaneously minimize the number of SU and optimize plan quality. The EM method is solved by the fast iterative shrinkage-thresholding algorithm and validated in comparison with a state-of-the-art method, so-called energy sequencing (ES).

RESULTS

EM is validated and compared with ES using representative clinical cases. In terms of delivery efficiency, EM had fewer SU than ES with an average of 35% reduction of SU. In terms of plan quality, compared to ES, EM had smaller optimization objective values and better target dose conformality, and generally lower dose to organs-at-risk and lower integral dose to body. In terms of computational efficiency, EM was substantially more efficient than ES by at least ten-fold.

CONCLUSION

We have developed a new ELO method for SPArc using energy matrix regularization, and shown that this new method EM can improve both delivery efficiency and plan quality, with substantially reduced computational time, compared to ES. This article is protected by copyright. All rights reserved.

Design of static and dynamic ridge filters for FLASH-IMPT: a simulation study

PURPOSE

This paper focused on the design and optimization of ridge filter-based intensity-modulated proton therapy (IMPT), and its potential applications for FLASH. Differing from the standard pencil beam scanning (PBS) mode, no energy/layer switching is required and total treatment time can be shortened.

METHODS

Unique dose influence matrices were generated as a proton beam traverse through slabs of different thicknesses (i.e. modulation by different layers). To establish the references for comparison, conventional IMPT plans (single field) were created using a large-scale non-linear solver. The spot weights from the reference IMPT plans were used as inputs for optimizing the design of ridge filters. Two designs were evaluated: model A (static) and model B (dynamic). The ridge filters designs were first verified (by GEANT4 simulation) in a water phantom and then in a H&N case. Direct comparison was made between the GEANT4 simulation results of two models and their respective references, with regard to plan quality, dose-averaged dose rate (DADR), and total treatment time.

RESULTS

In both the water phantom and the H&N case, two models are able to modulate dose distributions with high conformity, showing no significant difference relative to the reference plans. Dose rate volume histograms (DRVHs) suggest that in order to achieve a dose rate of 40 Gy/s over 90% PTV, the beam intensity needs to be 2.5×10¹¹ protons/s for both models. For a fraction dose of 10 Gy, the total treatment time (including both irradiation time and dead time) can be shortened by a factor of 4.9 (model A) and 6.5 (model B), relative to the reference plans.

CONCLUSION

Two proposed designs (both static and dynamic) can be used for PBS-IMPT requiring no layer switching. They are promising candidates for FLASH-IMPT capable of reducing treatment time and achieving high dose rates, while maintaining dose conformity simultaneously. This article is protected by copyright. All rights reserved.

Ultra‐High dose rate radiation production and delivery systems intended for FLASH

Higher dose rates, a trend for radiotherapy machines, can be beneficial in shortening treatment times for radiosurgery and mitigating the effects of motion. Recently, even higher doses (e.g., 100 times greater) have become targeted because of their potential to generate the FLASH effect (FE). We refer to these physical dose rates as ultra‐high (UHDR). The complete relationship between UHDR and the FE is unknown. But UHDR systems are needed to explore the relationship further and to deliver clinical UHDR treatments, where indicated. Despite the challenging set of unknowns, the authors seek to make reasonable assumptions to probe how existing and developing technology can address the UHDR conditions needed to provide beam generation capable of producing the FE in preclinical and clinical applications. As a preface, this paper discusses the known and unknown relationships between UHDR and the FE. Based on these, different accelerator and ionizing radiation types are then discussed regarding the relevant UHDR needs. The details of UHDR beam production are discussed for existing and potential future systems such as linacs, cyclotrons, synchrotrons, synchrocyclotrons, and laser accelerators. In addition, various UHDR delivery mechanisms are discussed, along with required developments in beam diagnostics and dose control systems.

Simultaneous dose and dose rate optimization (SDDRO) of the FLASH effect for pencil-beam-scanning proton therapy

PURPOSE

Compared to CONV-RT (with conventional dose rate), FLASH-RT (with ultra-high dose rate) can provide biological dose sparing for organs-at-risk (OARs) via the so-called FLASH effect, in addition to physical dose sparing. However, the FLASH effect only occurs, when both dose and dose rate meet certain minimum thresholds. This work will develop a simultaneous dose and dose rate optimization (SDDRO) method accounting for both FLASH dose and dose rate constraints during treatment planning for pencil-beam-scanning proton therapy.

METHODS

SDDRO optimizes the FLASH effect (specific to FLASH-RT) as well as the dose distribution (similar to CONV-RT). The nonlinear dose rate constraint is linearized, and the reformulated optimization problem is efficiently solved via iterative convex relaxation powered by alternating direction method of multipliers. To resolve and quantify the generic tradeoff of FLASH-RT between FLASH and dose optimization, we propose the use of FLASH effective dose based on dose modifying factor (DMF) owing to the FLASH effect.

RESULTS

FLASH-RT via transmission beams (TB) (IMPT-TB or SDDRO) and CONV-RT via Bragg peaks (BP) (IMPT-BP) were evaluated for clinical prostate, lung, head-and-neck (HN), and brain cases. Despite the use of TB, which is generally suboptimal to BP for normal tissue sparing, FLASH-RT via SDDRO considerably reduced FLASH effective dose for high-dose OAR adjacent to the target. For example, in the lung SBRT case, the max esophageal dose constraint 27 Gy was only met by SDDRO (24.8 Gy), compared to IMPT-BP (35.3 Gy) or IMPT-TB (36.6 Gy); in the brain SRS case, the brain constraint V12Gy≤15cc was also only met by SDDRO (13.7cc), compared to IMPT-BP (43.9cc) or IMPT-TB (18.4cc). In addition, SDDRO substantially improved the FLASH coverage from IMPT-TB, e.g., an increase from 37.2% to 67.1% for lung, from 39.1% to 58.3% for prostate, from 65.4% to 82.1% for HN, from 50.8% to 73.3% for the brain.

CONCLUSIONS

Both FLASH dose and dose rate constraints are incorporated into SDDRO for FLASH-RT that jointly optimizes the FLASH effect and physical dose distribution. FLASH effective dose via FLASH DMF is introduced to reconcile the tradeoff between physical dose sparing and FLASH sparing, and quantify the net effective gain from CONV-RT to FLASH-RT.

Development of Ultra-High Dose-Rate (FLASH) Particle Therapy

Research efforts in FLASH radiotherapy have increased at an accelerated pace recently. FLASH radiotherapy involves ultra-high dose rates and has shown to reduce toxicity to normal tissue while maintaining tumor response in pre-clinical studies when compared to conventional dose rate radiotherapy. The goal of this review is to summarize the studies performed to-date with proton, electron, and heavy ion FLASH radiotherapy, with particular emphasis on the physical aspects of each study and the advantages and disadvantages of each modality. Beam delivery parameters, experimental set-up, and the dosimetry tools used are described for each FLASH modality. In addition, modeling efforts and treatment planning for FLASH radiotherapy is discussed along with potential drawbacks when translated into the clinical setting. The final section concludes with further questions that have yet to be answered before safe clinical implementation of FLASH radiotherapy.

Treatment planning for Flash radiotherapy: general aspects and applications to proton beams

The increased radioresistence of healthy tissues when irradiated at very high dose rates (known as the Flash effect) is a radiobiological mechanism that is currently investigated in order to increase the therapeutic ratio of radiotherapy treatments. To maximize the benefits of the clinical application of Flash, a patient-specific balance between different properties of the dose distribution should be found, i.e. Flash needs to be one of the variables considered in treatment planning. We investigated the Flash potential of three proton therapy planning and beam delivery techniques, each on a different anatomical region. Based on a set of beam delivery parameters, on hypotheses on the dose and dose rate thresholds needed for the Flash effect to occur, and on two definitions of Flash dose rate, we generated exemplary illustrations of the capabilities of current proton therapy equipment to generate Flash dose distributions. All techniques investigated could both produce dose distributions comparable with a conventional proton plan and reach the Flash regime, to an extent that was strongly dependent on the dose per fraction and the Flash dose threshold. The beam current, Flash dose rate threshold and dose rate definition typically had a more moderate effect on the amount of Flash dose in normal tissue. A systematic estimation of the impact of Flash on different patient anatomies and treatment protocols is possible only if Flash-specific treatment planning features become readily available. Planning evaluation tools such as a voxel-based dose delivery time structure, and the inclusion in the optimization cost function of parameters directly associated with Flash (e.g. beam current, spot delivery sequence and scanning speed), are needed to generate treatment plans that are taking full advantage of the potential benefits of the Flash effect. This article is protected by copyright. All rights reserved.

A Universal Range Shifter and Range Compensator Can Enable Proton Pencil Beam Scanning Single-Energy Bragg Peak FLASH-RT Treatment Using Current Commercially Available Proton Systems

PURPOSE

Transmission beams have been proposed for ultra-high dose (or FLASH) proton planning, limiting the organ sparing potentials of proton therapy. By pulling back the ranges of the highest energy proton beams and compensating proton ranges to adapt to the target distally, the exit dose of proton beams can be eliminated to better protect organs at risk while still preserving FLASH dose rate delivery.

METHOD AND MATERIALS

An inverse planning tool was developed to optimize intensity modulated proton therapy using a single-energy layer for FLASH radiation therapy planning. The range pull-backs were calculated to stop single-energy proton beams at the distal edge of the target. The spot map and weights of each field were optimized to achieve a sufficient dose rate using proton beam Bragg peaks. A C-shape target in phantom, along with 6 consecutive lung cancer patients previously treated using proton stereotactic body radiation therapy were planned using this novel Bragg Peak method and also transmission technique. Dosimetry characteristics and 3-dimensional dose rate were investigated.

RESULTS

The minimum monitor units (MU) for transmission and Bragg peak plans were 400 MU/spot and 1200 MU/spot, respectively, corresponding to spot peak dose rates of 670 GyRBE (relative biological effectiveness) per second and 1950 GyRBE per second. Bragg peak plans yield a generally comparable target uniformity while significantly reducing dose spillage volume from the low to medium dose level. For all the 6 lung cases delivery of 34 GyRBE in 1 fraction, assessing Radiation Therapy Oncology Group 0915 constraints, the lung V_7GyRBE volume was reduced by up to 32% (P = .001) for Bragg peak plans. The transmission plans tended to generate 2.4% higher FLASH dose rate coverage (V_40GyRBE/s) versus Bragg peak plans over the major organs at risk. However, Bragg peak plans could also reach the FLASH radiation therapy threshold of V_40GyRBE/s using a higher MU/spot and sophisticated dose-rate optimization algorithm.

CONCLUSIONS

This first proof-of-concept study has demonstrated this novel method of combining range pull-back and powerful inverse optimization capable of achieving FLASH dose rate based on currently available machine parameters using a single-energy Bragg peak. Similar target coverage and uniformity can be maintained by Bragg peak FLASH plans while substantially improving the sparing of organs at risk compared with transmission plans.

Minimum-monitor-unit optimization via a stochastic coordinate descent method

Objective. Deliverable proton spots are subject to the minimum monitor-unit (MMU) constraint. The MMU optimization problem with relatively large MMU threshold remains mathematically challenging due to its strong nonconvexity. However, the MMU optimization is fundamental to proton radiotherapy (RT), including efficient IMPT and proton arc delivery (ARC). This work aims to develop a new optimization algorithm that is effective in solving the MMU problem.Approach.Our new algorithm is primarily based on stochastic coordinate decent (SCD) method. It involves three major steps: first to decouple the determination of active sets for dose-volume-histogram (DVH) planning constraints from the MMU problem via iterative convex relaxation method; second to handle the nonconvexity of the MMU constraint via SCD to localize the index set of nonzero spots; third to solve convex subproblems projected to this convex set of nonzero spots via projected gradient descent method.Main results.Our new method SCD is validated and compared with alternating direction method of multipliers (ADMM) for IMPT and ARC. The results suggest SCD had better plan quality than ADMM, e.g. the improvement of conformal index (CI) from 0.56 to 0.69 during IMPT, and from 0.28 to 0.80 during ARC for the lung case. Moreover, SCD successfully handled the nonconvexity from large MMU threshold that ADMM failed to handle, in the sense that (1) the plan quality from ARC was worse than IMPT (e.g. CI was 0.28 with IMPT and 0.56 with ARC for the lung case), when ADMM was used; (2) in contrast, with SCD, ARC achieved better plan quality than IMPT (e.g. CI was 0.69 with IMPT and 0.80 with ARC for the lung case), which is compatible with more optimization degrees of freedom from ARC compared to IMPT.Significance. To the best of our knowledge, our new MMU optimization method via SCD can effectively handle the nonconvexity from large MMU threshold that none of the current methods can solve. Therefore, we have developed a unique MMU optimization algorithm via SCD that can be used for efficient IMPT, proton ARC, and other particle RT applications where large MMU threshold is desirable (e.g. for the delivery of high dose rates or/and a large number of spots).

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.

FLASH Radiotherapy Using Single-Energy Proton PBS Transmission Beams for Hypofractionation Liver Cancer: Dose and Dose Rate Quantification

PURPOSE

This work aims to study the dose and ultra-high-dose rate characteristics of transmission proton pencil beam scanning (PBS) FLASH radiotherapy (RT) for hypofractionation liver cancer based on the parameters of a commercially available proton system operating under FLASH mode.

METHODS AND MATERIALS

An in-house treatment planning software (TPS) was developed to perform intensity-modulated proton therapy (IMPT) FLASH-RT planning. Single-energy transmission proton PBS plans of 4.5 Gy × 15 fractions were optimized for seven consecutive hepatocellular carcinoma patients, using 2 and 5 fields combined with 1) the minimum MU/spot chosen between 100 and 400, and minimum spot time (MST) of 2 ms, and 2) the minimum MU/spot of 100, and MST of 0.5 ms, based upon considerations in target uniformities, OAR dose constraints, and OAR FLASH dose rate coverage. Then, the 3D average dose rate distribution was calculated. The dose metrics for the mean dose of Liver-GTV and other major OARs were characterized to evaluate the dose quality for the different combinations of field numbers and minimum spot times compared to that of conventional IMPT plans. Dose rate quality was evaluated using 40 Gy/s volume coverage (V40Gy/s).

RESULTS

All plans achieved favorable and comparable target uniformities, and target uniformity improved as the number of fields increased. For OARs, no significant dose differences were observed between plans of different field numbers and the same MST. For plans using shorter MST and the same field numbers, better sparing was generally observed in most OARs and was statistically significant for the chest wall. However, the FLASH dose rate coverage V40Gy/s was increased by 20% for 2-field plans compared to 5-field plans in most OARs with 2-ms MST, which was less evident in the 0.5-ms cases. For 2-field plans, dose metrics and V40Gy/s of select OARs have large variations due to the beam angle selection and variable distances to the targets. The transmission plans generally yielded inferior dosimetric quality to the conventional IMPT plans.

CONCLUSION

This is the first attempt to assess liver FLASH treatment planning and demonstrates that it is challenging for hypofractionation with smaller fractional doses (4.5 Gy/fraction). Using fewer fields can allow higher minimum MU/spot, resulting in higher OAR FLASH dose rate coverages while achieving similar plan quality compared to plans with more fields. Shorter MST can result in better plan quality and comparable or even better FLASH dose rate coverage.

An adaptive spot placement method on Cartesian grid for pencil beam scanning proton therapy

Pencil beam scanning proton radiotherapy (RT) offers flexible proton spot placement near treatment targets for delivering tumoricidal radiation dose to tumor targets while sparing organs-at-risk. Currently the spot placement is mostly based on a non-adaptive sampling (NS) strategy on a Cartesian grid. However, the spot density or spacing during NS is a constant for the Cartesian grid that is independent of the geometry of tumor targets, and thus can be suboptimal in terms of plan quality (e.g. target dose conformality) and delivery efficiency (e.g. number of spots). This work develops an adaptive sampling (AS) spot placement method on the Cartesian grid that fully accounts for the geometry of tumor targets. Compared with NS, AS places (1) a relatively fine grid of spots at the boundary of tumor targets to account for the geometry of tumor targets and treatment uncertainties (setup and range uncertainty) for improving dose conformality, and (2) a relatively coarse grid of spots in the interior of tumor targets to reduce the number of spots for improving delivery efficiency and robustness to the minimum-minitor-unit (MMU) constraint. The results demonstrate that (1) AS achieved comparable plan quality with NS for regular MMU and substantially improved plan quality from NS for large MMU, using merely about 10% of spots from NS, where AS was derived from the same Cartesian grid as NS; (2) on the other hand, with similar number of spots, AS had better plan quality than NS consistently for regular and large MMU.

A Novel Proton Pencil Beam Scanning FLASH RT Delivery Method Enables Optimal OAR Sparing and Ultra-High Dose Rate Delivery: A Comprehensive Dosimetry Study for Lung Tumors

PURPOSE

While transmission proton beams have been demonstrated to achieve ultra-high dose rate FLASH therapy delivery, they are unable to spare normal tissues distal to the target. This study aims to compare FLASH treatment planning using single energy Bragg peak proton beams versus transmission proton beams in lung tumors and to evaluate Bragg peak plan optimization, characterize plan quality, and quantify organ-at-risk (OAR) sparing.

MATERIALS AND METHODS

Both Bragg peak and transmission plans were optimized using an in-house platform for 10 consecutive lung patients previously treated with proton stereotactic body radiation therapy (SBRT). To bring the dose rate up to the FLASH-RT threshold, Bragg peak plans with a minimum MU/spot of 1200 and transmission plans with a minimum MU/spot of 400 were developed. Two common prescriptions, 34 Gy in 1 fraction and 54 Gy in 3 fractions, were studied with the same beam arrangement for both Bragg peak and transmission plans (n = 40 plans). RTOG 0915 dosimetry metrics and dose rate metrics based on different dose rate calculations, including average dose rate (ADR), dose-averaged dose rate (DADR), and dose threshold dose rate (DTDR), were investigated. We then evaluated the effect of beam angular optimization on the Bragg peak plans to explore the potential for superior OAR sparing.

RESULTS

Bragg peak plans significantly reduced doses to several OAR dose parameters, including lung V7.4Gy and V7Gy by 32.0% (p < 0.01) and 30.4% (p < 0.01) for 34Gy/fx plans, respectively; and by 40.8% (p < 0.01) and 41.2% (p < 0.01) for 18Gy/fx plans, respectively, compared with transmission plans. Bragg peak plans have ~3% less in DADR and ~10% differences in mean OARs in DTDR and DADR relative to transmission plans due to the larger portion of lower dose regions of Bragg peak plans. With angular optimization, optimized Bragg peak plans can further reduce the lung V7Gy by 20.7% (p < 0.01) and V7.4Gy by 19.7% (p < 0.01) compared with Bragg peak plans without angular optimization while achieving a similar 3D dose rate distribution.

CONCLUSION

The single-energy Bragg peak plans achieve superior dosimetry performances in OARs to transmission plans with comparable dose rate performances for lung cancer FLASH therapy. Beam angle optimization can further improve the OAR dosimetry parameters with similar 3D FLASH dose rate coverage.

Dosimetry, Efficacy, Safety, and Cost-Effectiveness of Proton Therapy for Non-Small Cell Lung Cancer

Non-small cell lung cancer (NSCLC) is the most common malignancy which requires radiotherapy (RT) as an important part of its multimodality treatment. With the advent of the novel irradiation technique, the clinical outcome of NSCLC patients who receive RT has been dramatically improved. The emergence of proton therapy, which allows for a sharper dose of build-up and drop-off compared to photon therapy, has potentially improved clinical outcomes of NSCLC. Dosimetry studies have indicated that proton therapy can significantly reduce the doses for normal organs, especially the lung, heart, and esophagus while maintaining similar robust target volume coverage in both early and advanced NSCLC compared with photon therapy. However, to date, most studies have been single-arm and concluded no significant changes in the efficacy for early-stage NSCLC by proton therapy over stereotactic body radiation therapy (SBRT). The results of proton therapy for advanced NSCLC in these studies were promising, with improved clinical outcomes and reduced toxicities compared with historical photon therapy data. However, these studies were also mainly single-arm and lacked a direct comparison between the two therapies. Currently, there is much emerging evidence focusing on dosimetry, efficacy, safety, and cost-effectiveness of proton therapy for NSCLC that has been published, however, a comprehensive review comparing these therapies is, to date, lacking. Thus, this review focuses on these aspects of proton therapy for NSCLC.

SDDRO-joint: simultaneous dose and dose rate optimization with the joint use of transmission beams and Bragg peaks for FLASH proton therapy

Cancer radiotherapy (RT) with the irradiation at ultra-high dose rates, namely FLASH-RT, can substantially reduce radiation-induced normal tissue toxicities while maintaining tumor response. Currently, clinical FLASH-RT on deep-seated tumors can only be performed with proton beams. One way to achieve ultra-high dose rates at depth is through the use of high-energy transmission beams (TB), where the Bragg peaks (BP) fall outside the body. However, planning with TB alone does not fully leverage the degrees of freedom for dose shaping as traditional intensity modulated proton therapy (IMPT) which uses the BP of multi-energy proton beams at the tumor target. This work will develop a simultaneous dose and dose rate optimization (SDDRO) method with the joint use of TB and BP, namely SDDRO-Joint. Specifically, BP are placed inside tumor targets to improve the target dose conformality and sparse the normal-tissue dose, while TB primarily cover the tumor boundary to achieve ultra-high dose rate coverage of organs-at-risk (OAR) close to tumor targets. The sparing of OAR and other normal tissues via SDDRO-Joint is jointly by TB and BP, i.e. the FLASH sparing by TB and the dose sparing by BP. The results suggest that the addition of BP substantially increased the target dose conformality for SDDRO. Noticeably SDDRO-Joint also provided slightly higher conformal index values than the conventional IMPT method with BP alone.

Ultra-High Dose Rate Transmission Beam Proton Therapy for Conventionally Fractionated Head and Neck Cancer: Treatment Planning and Dose Rate Distributions

Simple Summary

Standard intensity-modulated proton therapy (IMPT) places the Bragg-peak in the target. However, it is also possible to use high energy proton transmission beams (TBs), where the Bragg-peak is placed outside the patient, irradiating with the beam section proximal to the Bragg-peak. TBs use only one energy, increase robustness, are insensitive to density changes and have sharper penumbras. TBs can also be delivered at ultra-high dose-rates (UHDRs, e.g., ≥40 Gy/s), which is one of the requirements for the FLASH-effect. The aim of this work was twofold: (1) comparison of TB-plan quality to IMPT and photon volumetric-modulated arc therapy (VMAT) for conventionally fractionated head-and-neck cancer; (2) analysis of TB-plan UHDR-metrics. We showed that TB-plan quality was comparable to IMPT for contoured organs at risk and better than VMAT. Any potential FLASH-effect would only further improve plan quality. TB plans can also be delivered quickly, which might facilitate higher patient through-put and enhance patient comfort.

Abstract

Transmission beam (TB) proton therapy (PT) uses single, high energy beams with Bragg-peak behind the target, sharp penumbras and simplified planning/delivery. TB facilitates ultra-high dose-rates (UHDRs, e.g., ≥40 Gy/s), which is a requirement for the FLASH-effect. We investigated (1) plan quality for conventionally-fractionated head-and-neck cancer treatment using spot-scanning proton TBs, intensity-modulated PT (IMPT) and photon volumetric-modulated arc therapy (VMAT); (2) UHDR-metrics. VMAT, 3-field IMPT and 10-field TB-plans, delivering 70/54.25 Gy in 35 fractions to boost/elective volumes, were compared (n = 10 patients). To increase spot peak dose-rates (SPDRs), TB-plans were split into three subplans, with varying spot monitor units and different gantry currents. Average TB-plan organs-at-risk (OAR) sparing was comparable to IMPT: mean oral cavity/body dose were 4.1/2.5 Gy higher (9.3/2.0 Gy lower than VMAT); most other OAR mean doses differed by <2 Gy. Average percentage of dose delivered at UHDRs was 46%/12% for split/non-split TB-plans and mean dose-averaged dose-rate 46/21 Gy/s. Average total beam-on irradiation time was 1.9/3.8 s for split/non-split plans and overall time including scanning 8.9/7.6 s. Conventionally-fractionated proton TB-plans achieved comparable OAR-sparing to IMPT and better than VMAT, with total beam-on irradiation times <10s. If a FLASH-effect can be demonstrated at conventional dose/fraction, this would further improve plan quality and TB-protons would be a suitable delivery system.

Translational Research in FLASH Radiotherapy-From Radiobiological Mechanisms to In Vivo Results

FLASH radiotherapy, or the administration of ultra-high dose rate radiotherapy, is a new radiation delivery method that aims to widen the therapeutic window in radiotherapy. Thus far, most in vitro and in vivo results show a real potential of FLASH to offer superior normal tissue sparing compared to conventionally delivered radiation. While there are several postulations behind the differential behaviour among normal and cancer cells under FLASH, the full spectra of radiobiological mechanisms are yet to be clarified. Currently the number of devices delivering FLASH dose rate is few and is mainly limited to experimental and modified linear accelerators. Nevertheless, FLASH research is increasing with new developments in all the main areas: radiobiology, technology and clinical research. This paper presents the current status of FLASH radiotherapy with the aforementioned aspects in mind, but also to highlight the existing challenges and future prospects to overcome them.