Direct probabilistic IMPT treatment planning with setup and range errors for neuro-oncological patients

To show clinical feasibility of a previously proposed probabilistic planning approach that can precisely optimize for clinical goals with patient-specific acceptance probabilities on a neuro-oncological patient group, we compared probabilistic plans with (automated) robust plans for one patient (group A) that could achieve sufficient clinical target coverage and for four patients (group B) where target coverage had to be compromised due to organ-at-risk (OAR) dose constraints.

The probabilistic approach is percentile-based and uses the fact that a (dose) percentile can be approximated as a linear combination of its expected value and standard deviation.

The optimization has a nested structure: the inner optimization optimizes the beam weights for a given percentile estimate, while an outer loop iteratively updates and improves the accuracy of the percentile estimate.

For every outer iteration, the optimization is warm-started from the previous iteration.

Percentiles are efficiently calculated by sampling a polynomial chaos expansion of the dose-influence matrix.

The patient in group A achieved cumulative OAR dose reductions (of OAR-related DVH-metrics) of 19 GyRBE, for identical target coverage.

Target coverage improved for all patients in group B (the 10th percentile of $D_{99.8\%}$ increased up to 0.93 GyRBE), at the same time reaching cumulative OAR dose reductions (of OAR-related DVH-metrics) up to 33 GyRBE.

Probabilistic plans were optimized in 44h to 141h.

For two representative patients, eliminating warm-starting (i.e., the outer loop) from the approach reduced total optimization times to below 10h (which took originally 80h and 141h).

Compared to robust optimization methods, the probabilistic approach achieves improved trade-offs between probabilistic target coverage and OAR sparing, potentially leading to better treatments.