A mixed integer programming approach to minibeam aperture optimization for multi-collimator proton minibeam radiotherapy

TL;DR

This paper introduces a mixed integer programming method to optimize minibeam aperture selection in proton radiotherapy, significantly improving dose conformity and organ sparing over traditional heuristic approaches.

Contribution

The study presents the first MIP-based optimization framework for minibeam aperture selection in MC-pMBRT, enhancing plan quality and computational efficiency.

Findings

Improved conformity index in head-and-neck case from 0.61 to 0.70.

Enhanced PVDR and reduced doses in OAR.

Demonstrated superiority over heuristic methods in treatment planning.

Abstract

Background: Multi-collimator proton minibeam radiation therapy (MC-pMBRT) has recently emerged as a versatile technique for dose shaping, enabling peak-valley dose patterns in organs-at-risk (OAR) while maintaining a uniform dose distribution in tumor. MC-pMBRT leverages a set of generic multi-slit collimators (MSC) with varying center-to-center distances. However, the current method for minibeam aperture optimization (MAO), i.e., the selection of MSC per beam angle, is manual and heuristic, resulting in computational inefficiencies and no guarantee of optimality. This work introduces a novel mixed integer programming (MIP) approach to MAO for optimizing MC-pMBRT plan quality. Methods: The proposed MIP approach jointly optimizes dose distributions, peak-to-valley dose ratio (PVDR), and selects the optimal set of MSC per beam angle. The optimization problem includes decision variables for MSC selection per beam angle and spot weights. The proposed MIP approach is a two-step process: Step1: the binary variables are optimally determined to select MSC for each beam angle; Step 2: the continuous variables are solved to determine the spot weights. Both steps utilize iterative convex relaxation and the alternating direction method of multipliers to solve the problems. Results: The proposed MIP method for MAO (MIP-MAO) was validated against the conventional heuristic method (CONV) for MC-pMBRT treatment planning. Results indicate that MIP-MAO enhances the conformity index (CI) for the target and improves PVDR for OAR. For instance, in a head-and-neck case, CI improved from 0.61 (CONV) to 0.70 (MIP-MAO); in an abdomen case, CI improved from 0.78 (CONV) to 0.83 (MIP-MAO). Additionally, MIP-MAO reduced mean doses in the body and OAR. Conclusions: A novel MIP approach for MAO in MC-pMBRT is presented, showing demonstrated improvements in plan quality and PVDR compared to the heuristic method.