Abstract:Purpose Proton dose distribution is sensitive to uncertainties in range estimation and patient positioning. Currently, the proton robustness is managed by worst-case scenario optimization methods, which are computationally inefficient. To overcome these challenges, we develop a novel intensity-modulated proton therapy (IMPT) optimization method that integrates dose fidelity with a sensitivity term that describes dose perturbation as the result of range and positioning uncertainties. Methods In the integrated optimization framework, the optimization cost function is formulated to include two terms: a dose fidelity term and a robustness term penalizing the inner product of the scanning spot sensitivity and intensity. The sensitivity of an IMPT scanning spot to perturbations is defined as the dose distribution variation induced by range and positioning errors. To evaluate the sensitivity, the spatial gradient of the dose distribution of a specific spot is first calculated. The spot sensitivity is then determined by the total absolute value of the directional gradients of all affected voxels. The fast iterative shrinkage-thresholding algorithm is used to solve the optimization problem. This method was tested on three skull base tumor (SBT) patients and three bilateral head-and-neck (H&N) patients. The proposed sensitivity-regularized method (SenR) was implemented on both clinic target volume (CTV) and planning target volume (PTV). They were compared with conventional PTV-based optimization method (Conv) and CTV-based voxel-wise worst-case scenario optimization approach (WC). Results Under the nominal condition without uncertainties, the three methods achieved similar CTV dose coverage, while the CTV-based SenR approach better spared organs at risks (OARs) compared with the WC approach, with an average reduction of [Dmean, Dmax] of [4.72, 3.38] GyRBE for the SBT cases and [2.54, 3.33] GyRBE for the H&N cases. The OAR sparing of the PTV-based SenR method was comparable with the WC method. The WC method, and SenR approaches all improved the plan robustness from the conventional PTV-based method. On average, under range uncertainties, the lowest [D95%, V95%, V100%] of CTV were increased from [93.75%, 88.47%, 47.37%] in the Conv method, to [99.28%, 99.51%, 86.64%] in the WC method, [97.71%, 97.85%, 81.65%] in the SenR-CTV method and [98.77%, 99.30%, 85.12%] in the SenR-PTV method, respectively. Under setup uncertainties, the average lowest [D95%, V95%, V100%] of CTV were increased from [95.35%, 94.92%, 65.12%] in the Conv method, to [99.43%, 99.63%, 87.12%] in the WC method, [96.97%, 97.13%, 77.86%] in the SenR-CTV method, and [98.21%, 98.34%, 83.88%] in the SenR-PTV method, respectively. The runtime of the SenR optimization is eight times shorter than that of the voxel-wise worst-case method. Conclusion We developed a novel computationally efficient robust optimization method for IMPT. The robustness is calculated as the spot sensitivity to both range and shift perturbations. The dose fidelity term is then regularized by the sensitivity term for the flexibility and trade-off between the dosimetry and the robustness. In the stress test, SenR is more resilient to unexpected uncertainties. These advantages in combination with its fast computation time make it a viable candidate for clinical IMPT planning.

Optimization of patient-specific range modulators for conformal FLASH proton therapy

Optimization of pencil beam scanning pattern for FLASH proton therapy

Conformal Dose Modulator for Proton Beam Therapy Part 1: A Simulation Study

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

An Integrated Biological Optimization framework for proton SBRT FLASH treatment planning allows dose, dose rate, and LET optimization using patient-specific ridge filters

A Novel Dose Rate Optimization Method to Maximize Ultra-High Dose Rate Coverage of Critical Organs-at-risk Without Compromising Dosimetry Metrics in Proton Pencil Beam Scanning FLASH Radiotherapy

Simultaneous Dose and Dose Rate Optimization (SDDRO) of the FLASH Effect for Pencil-Beam-scanning Proton Therapy

Simultaneous Dose and Dose Rate Optimization (SDDRO) for FLASH Proton Therapy.

Simultaneous dose and dose rate optimization via dose modifying factor modeling for FLASH effective dose

A Fast 3D Range-Modulator Delivery Approach: Validation of the FLUKA Model on a Varian ProBeam System Including a Robustness Analysis

Concurrent Monte Carlo Transport and Fluence Optimization with Fluence Adjusting Scalable Transport Monte Carlo.

Robust Optimization for Intensity-Modulated Proton Therapy with Soft Spot Sensitivity Regularization

Validation of the Quantum Physics Processes Underlying the Integrated Optimization of Proton FLASH Radiotherapy

Using Modularized Pin Ridge Filter in Proton FLASH Planning for Liver Stereotactic Ablative Body Radiotherapy

A novel fast robust optimization algorithm for intensity-modulated proton therapy with minimum monitor unit constraint

MO-D-351-05: A “slicing” Approach for IMRT Optimization Using Linear Programming

Feasibility and constraints of Bragg peak FLASH proton therapy treatment planning

Plan-Delivery-Time-Constrained Inverse Optimization for High-Throughput PBS Proton RT

Feasibility Study of Modularized Pin Ridge Filter Implementation in Proton FLASH Planning for Liver Stereotactic Ablative Body Radiotherapy

Effectiveness of Robust Optimization in Intensity-Modulated Proton Therapy Planning for Head and Neck Cancers

A Fast Robust Optimizer for Intensity Modulated Proton Therapy Using GPU.