Citations

Abstract Purpose This study investigates a novel optimization framework that integrates the newly developed MCF MKM relative biological effectiveness (RBE) model into intensity modulated particle therapy (IMPT) planning for carbon ion therapy. Unlike conventional RBE models such as MKM or LEM, the MCF MKM provides detailed microdosimetric information, which improves predictive accuracy of cell survival but elevates both memory and time requirement in the optimization process. Furthermore, to the implicit nature of the MCF MKM, the gradient information, which is essential in inverse planning, is hardly accessible. In response, a randomized gradient‐free optimization method (RGFM) was developed to mitigate the computational challenge and achieve improved RBE dose distributions. Materials and Methods An RBE dose optimization problem was formulated based on the discretized MCF MKM, known as the Abridged Microdosimetry Distribution Model (AMDM). The RGFM was derived to find feasible and high‐quality settings of monitor unit for IMPT plans without the need for explicit gradient calculations. To verify the clinical plausibility of the proposed method, it was tested on two clinical scenarios—prostate and brain cases in comparison with benchmark plans derived solely from physical dose optimization. Computational times and dose‐volume histograms (DVHs) for both physical and RBE doses were evaluated. Results The RGFM algorithm successfully optimized the MCF MKM‐based IMPT plans. Compared to the benchmark, the RBE dose distributions showed improved coverage of the target and reduced hot spots. Such improvements were observed to coincide with trade‐offs in physical dose distributions. The complexity of the MCF MKM extended the computational time, yet high‐quality solutions were attainable in a clinically reasonable timeframe. Multi‐field optimization (MFO) further enhanced target coverage and critical organ sparing over single‐field optimization (SFO). Conclusions This proof‐of‐principle study confirmed the feasibility of using a randomized gradient‐free method to optimize RBE dose based on the MCF MKM in IMPT planning for carbon ion therapy. The approach provides a pathway for clinical adoption of advanced RBE models, potentially leading to more biologically informed and effective carbon ion therapy treatments. Future work involves refining algorithms, incorporating parallel computing, and performing additional case studies to improve computational efficiency and generalizability.

Abstract Background Dose‐driven continuous scanning (DDCS) enhances the efficiency and precision of proton pencil beam delivery by reducing beam pauses inherent in discrete spot scanning (DSS). However, current DDCS optimization studies using traveling salesman problem (TSP) formulations often rely on fixed beam intensity and computationally expensive interpolation for move spot generation, limiting efficiency and methodological robustness. Purpose This study introduces a Break Spot‐Guided (BSG) method, combined with two acceleration strategies—dose rate skipping and bounding—to optimize beam intensity while minimizing beam delivery time (BDT). In addition, a ‐method is proposed for efficient move spot generation, aiming to balance computational efficiency and dose calculation accuracy. Methods The BSG framework simplifies the original mixed‐integer nonlinear programming (MINLP) problem by transforming it into a series of TSP evaluations using break‐spot dose rates. The two acceleration strategies are integrated to reduce computational complexity. The ‐method is introduced to efficiently generate move spots without compromising dose fidelity. The performance of the BSG method, with and without the acceleration strategies, and the ‐method, was evaluated across multiple clinical cases. Results By jointly optimizing the scan path and dose rate, the proposed BSG‐S‐B method offered the most efficient delivery, achieving a 64% reduction in BDT compared to DSS and a 57% reduction compared to Liu's method in the prostate case. These gains were consistently observed across all cases. Utilizing acceleration strategies, BSG‐S‐B significantly reduced computation time from hours to minutes (e.g., from 10 525.7 to 108.5 s in the prostate case, and from 6202.4 to 129.8 s in the head‐and‐neck case). When combined with the ‐method, it reduced the number of move spots by more than 50%, while still ensuring dosimetric quality. Conclusions The proposed framework significantly improves computational efficiency, making scan path optimization with dose rate selection feasible for clinical DDCS delivery. The ‐method effectively reduces computational complexity while maintaining dose calculation accuracy, demonstrating strong potential for clinical adoption in proton therapy treatment planning.