Citations

Abstract The anticipated clinical demand for carbon ion radiotherapy (CIRT) in the United States is expected to include a substantial proportion of clinically complex disease sites that require accurate, highly conformal, and adaptive treatment. These clinical needs, together with the unique physical and biological characteristics of carbon ion beams, form the foundation of the design philosophy for the Mayo Clinic Florida (MCF) CIRT program. The program remains under active development, with clinical deployment planned for 2028. Building upon the clinical rationale and anticipated U.S. patient population discussed in Part I of this two‐part review, this Vision Article describes how MCF intends to develop an integrated CIRT treatment system. The program is guided by two overarching principles: fast, accurate, and safe beam delivery and accurate treatment planning and treatment delivery, supported by an integrated patient‐specific quality assurance framework. Key elements include high‐dose‐rate dose‐driven continuous scanning with advanced beam monitoring and irradiation control, delivery‐aware treatment planning, the Monte Carlo‐based BioPhys computational platform integrating physical and biological dose modeling, a hierarchical adaptive workflow, and patient‐specific quality assurance spanning pretreatment, on‐treatment, and post‐treatment verification. This framework provides a technical roadmap for the integrated clinical implementation of CIRT at MCF.

Abstract Purpose Modeling relative biological effectiveness (RBE) is central to carbon ion radiotherapy treatment planning. The modified Microdosimetric Kinetic Model (mMKM) is a clinically established RBE framework that has guided treatment protocols at many existing carbon centers, while the Mayo Clinic Florida Microdosimetric Kinetic Model (MCF MKM) is a recently developed alternative. This work aims to implement the MCF MKM in the open‐source treatment planning system matRad and to quantitatively compare its RBE‐weighted dose predictions with those of the clinically established mMKM using identical physical dose distributions across multiple disease sites. These findings will help assess their dosimetric equivalence and inform protocol development for carbon ion radiotherapy at MCF. Methods Monte Carlo simulations of the MCF carbon beamline were performed to generate physical (IDD, LET, lateral spread) and biological base data for integration of each RBE model into matRad . Treatment plans were generated for six patients, each corresponding to a different disease site, using clinical beam configurations with carbon‐reference dose prescriptions, and plans were optimized using the MCF MKM. To isolate differences attributable solely to the RBE model, the resulting physical dose distributions were held fixed and RBE‐weighted doses were recalculated using the mMKM. Dose volume histogram (DVH) metrics and spatial dose‐difference maps were used to compare target coverage and organ‐at‐risk doses between the two models. Results Across patient cases, RBE‐weighted dose distributions from MCF MKM and mMKM showed strong agreement. Differences in target coverage were small, with CTV D95% differing by less than 1.6% across all disease sites and maximum target dose differences not exceeding 0.88%. Organ‐at‐risk dose deviations were limited, with differences of 3.0% or less across evaluated DVH metrics. Spatial dose‐difference maps showed that the largest discrepancies occurred in regions of steep dose gradients near target to organ‐at‐risk interfaces, while overall dose conformity and plan quality remained comparable between the two models. Conclusions: This study served as the first systematic model comparison of the MCF MKM and mMKM within a treatment planning environment. These findings suggest that the MCF MKM and mMKM produce dosimetrically consistent RBE‐weighted dose predictions under realistic planning conditions using carbon‐reference parameters. Accordingly, the fractionation schemes developed from years of clinical experience with mMKM implementations may serve as a practical foundation for protocol development at MCF.

Introduction Proton dose driven continuous scanning (DDCS) is a form of proton pencil beam scanning (PBS) in which the beam may remain on during transitions between successive spots. Mayo Clinic Florida (MCF) is preparing to commission a synchrotron-based proton therapy system designed for DDCS, incorporating novel irradiation control features. Given the novelty and complexity of this system, a systematic risk-based analysis was performed to inform upcoming commissioning strategies. Methods A prospective risk analysis was conducted using the American Association of Physicists in Medicine (AAPM) Task Group 100 (TG 100) framework. A comprehensive process map of the DDCS irradiation control workflow was developed. Failure modes and effects analysis (FMEA) was performed to identify potential delivery system failure modes and to score each failure mode based on occurrence, severity, and lack of detectability. Targeted mitigation and commissioning strategies were then developed to address the highest-risk items. Results Twenty-nine delivery system failure modes were identified. The highest-risk failure modes were specific to DDCS operation. Targeted commissioning strategies were defined, including scan path measurement and modeling, validation of move and flap dose, dose monitor collection efficiency testing, spot position monitor performance evaluation, and verification of control system intervention behavior. Incorporation of these strategies into the commissioning plan reduced the mean RPN by more than 50%. Discussion This work provides a detailed description of the MCF proton DDCS irradiation control system and uses a TG 100 based approach for identifying system-specific risks. As new systems incorporate technologies for which limited guidance exists, institutions must perform their own risk analyses to determine appropriate quality management strategies.

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.

Abstract Background Cone‐beam computed tomography (CBCT) is essential for image‐guided particle therapy (IGPT), providing daily patient positioning, anatomical monitoring, and treatment verification. However, conventional single‐detector CBCT suffers from poor soft‐tissue contrast, long acquisition times, and motion artifacts, reducing its effectiveness for adaptive radiotherapy workflows. Addressing these challenges requires an advanced CBCT reconstruction approach capable of enhancing image quality while reducing scan time and radiation exposure. Purpose This study introduces a Dual Detector Combined Scanning (DDCS) CBCT reconstruction algorithm to overcome conventional CBCT limitations. By integrating a dual‐source, orthogonal imaging setup, DDCS significantly reduces scanning time and imaging artifacts. A modified Parker‐weighting algorithm further improves image reconstruction accuracy, ensuring high‐fidelity visualization. The goal is to enhance adaptive radiotherapy by improving real‐time patient positioning, dose verification, and motion management. Methods A dual‐detector CBCT system was designed to acquire orthogonal projections simultaneously, reducing scan time and improving data completeness. A modified Parker‐weighting algorithm corrected angular overlap distortions, ensuring better reconstruction accuracy. The system was evaluated using numerical and physical phantom studies to assess spatial resolution, contrast‐to‐noise ratio (CNR), and artifact reduction. Additionally, a dynamic phantom simulated respiratory motion, validating motion robustness for adaptive radiotherapy. Performance was compared against single‐detector CBCT, focusing on image fidelity, noise suppression, and computational efficiency. Results DDCS‐CBCT demonstrated higher CNR, reduced motion artifacts, and improved spatial resolution, leading to more accurate anatomical visualization. The system's efficiency enables faster, more reliable patient setup in IGPT while maintaining a lower imaging dose. Conclusion The proposed DDCS‐CBCT approach significantly improves imaging accuracy, reduces scan time, and enhances real‐time volumetric guidance in IGRT/IGPT. These findings support the clinical feasibility of the DDCS system and its integration into online adaptive radiotherapy workflows.

Abstract Background All Hitachi proton pencil beam scanning facilities currently use discrete spot scanning (DSS). Mayo Clinic Florida (MCF) is installing a Hitachi particle therapy system with advanced technologies, including fast scan speeds, high beam intensity, rapid beam off control (RBOC), a skip spot function, and proton pencil beam scanning using dose driven continuous scanning (DDCS). A potential concern of RBOC is the generation of a shoulder at the end of the normal spot delivery due to a flap spot (FS) with a flap dose (FD), which has been investigated for carbon synchrotron but not for proton delivery. While investigated, for instance, for Hitachi's installation at MCF, this methodology could be applicable for all future high intensity proton deliveries. Purpose No Hitachi proton facility currently uses the proposed RBOC. This study aimed to understand the dosimetric impact of proton FD at MCF by simulating the FS with a Hitachi proton machine in research mode, reflecting the higher proton intensities expected with RBOC at MCF. Method Experiments were conducted to simulate MCF RBOC at Kyoto Prefecture University of Medicine (KPUM) in research mode, reducing delay time (Td) from 1.5 ms to 0.1 ms. 5,000 contiguous spots were delivered on the central axis for proton energies of 70.2, 142.5, and 220.0 MeV; at normal, high dose rate (HDR), and ultra‐high dose rate (uHDR) intensities; and at vertical and horizontal gantry angles for different Td. Measurements were taken using a fast oscilloscope and the nozzle's spot position monitor (SPM) and dose monitor (DM). A model was developed to predict FD dependence on beam intensity and assess the dosimetric impact for prostate and brain treatment plans. Two simulation types were planned: a flap DSS plan with FS at every spot and a flap DDCS plan with FS only at the end of each layer. Result FD was observed for RBOC with Td = 0.1 ms, showing no gantry angle dependence. FD increased with higher delayed dose rate (DDR), that is, beam intensity. The planning study showed dose volume histogram deterioration with increased FD compared to the clinical plan, but it was only significant for uHDR intensities. Deterioration was marginal in flap DSS plans for the HDR intensities planned at MCF, and flap DDCS plans were even less sensitive than flap DSS plans. Conclusion MCF is installing proton DDCS with higher beam intensities, a skip spot function, and fast beam‐off control. The resulting FD had an insignificant impact on dose distribution for two patient plans with both DSS and DDCS at the anticipated MCF intensities. However, significant dependence was observed in the case of uHDR. A method to measure the position and dose of the FS during commissioning is described in addition to recommendations for regular QA and log‐based proton patient‐specific quality assurance.

AbstractBackgroundAdaptive radiotherapy (ART) workflows have been increasingly adopted to achieve dose escalation and tissue sparing under shifting anatomic conditions, but the necessity of recontouring and the associated time burden hinders a real‐time or online ART workflow. In response to this challenge, approaches to auto‐segmentation involving deformable image registration, atlas‐based segmentation, and deep learning‐based segmentation (DLS) have been developed. Despite the particular promise shown by DLS methods, implementing these approaches in a clinical setting remains a challenge, namely due to the difficulty of curating a data set of sufficient size and quality so as to achieve generalizability in a trained model.PurposeTo address this challenge, we have developed an intentional deep overfit learning (IDOL) framework tailored to the auto‐segmentation task. However, certain limitations were identified, particularly the insufficiency of the personalized dataset to effectively overfit the model. In this study, we introduce a personalized hyperspace learning (PHL)‐IDOL segmentation framework capable of generating datasets that induce the model to overfit specific patient characteristics for medical image segmentation.MethodsThe PHL‐IDOL model is trained in two stages. In the first, a conventional, general model is trained with a diverse set of patient data (n = 100 patients) consisting of CT images and clinical contours. Following this, the general model is tuned with a data set consisting of two components: (a) selection of a subset of the patient data (m < n) using the similarity metrics (mean square error (MSE), peak signal‐to‐noise ratio (PSNR), structural similarity index (SSIM), and the universal quality image index (UQI) values); (b) adjust the CT and the clinical contours using a deformed vector generated from the reference patient and the selected patients using (a). After training, the general model, the continual model, the conventional IDOL model, and the proposed PHL‐IDOL model were evaluated using the volumetric dice similarity coefficient (VDSC) and the Hausdorff distance 95% (HD95%) computed for 18 structures in 20 test patients.ResultsImplementing the PHL‐IDOL framework resulted in improved segmentation performance for each patient. The Dice scores increased from 0.810.05 with the general model, 0.83 for the continual model, 0.83 for the conventional IDOL model to an average of 0.87 with the PHL‐IDOL model. Similarly, the Hausdorff distance decreased from 3.06 with the general model, 2.84 for the continual model, 2.79 for the conventional IDOL model and 2.36 for the PHL‐IDOL model. All the standard deviations were decreased by nearly half of the values comparing the general model and the PHL‐IDOL model.ConclusionThe PHL‐IDOL framework applied to the auto‐segmentation task achieves improved performance compared to the general DLS approach, demonstrating the promise of leveraging patient‐specific prior information in a task central to online ART workflows.

AbstractPurposeThe time structures of proton spot delivery in proton pencil beam scanning (PBS) radiation therapy are essential in many clinical applications. This study aims to characterize the time structures of proton PBS delivered by both synchrotron and synchrocyclotron accelerators using a non‐invasive technique based on scattered particle tracking.MethodsA pixelated semiconductor detector, AdvaPIX‐Timepix3, with a temporal resolution of 1.56 ns, was employed to measure time of arrival of secondary particles generated by a proton beam. The detector was placed laterally to the high‐flux area of the beam in order to allow for single particle detection and not interfere with the treatment. The detector recorded counts of radiation events, their deposited energy and the timestamp associated with the single events. Individual recorded events and their temporal characteristics were used to analyze beam time structures, including energy layer switch time, magnet switch time, spot switch time, and the scanning speeds in the x and y directions. All the measurements were repeated 30 times on three dates, reducing statistical uncertainty.ResultsThe uncertainty of the measured energy layer switch times, magnet switch time, and the spot switch time were all within 1% of average values. The scanning speeds uncertainties were within 1.5% and are more precise than previously reported results. The measurements also revealed continuous sub‐milliseconds proton spills at a low dose rate for the synchrotron accelerator and radiofrequency pulses at 7 µs and 1 ms repetition time for the synchrocyclotron accelerator.ConclusionThe AdvaPIX‐Timepix3 detector can be used to directly measure and monitor time structures on microseconds scale of the PBS proton beam delivery. This method yielded results with high precision and is completely independent of the machine log files.

AbstractBackgroundIntensity Modulated Proton Therapy (IMPT) is a sophisticated radiation treatment allowing for precise dose distributions. However, conventional spot selection strategies in IMPT face challenges, particularly with minimum monitor unit (MU) constraints, affecting planning quality and efficiency.PurposeThis study introduces an innovative Two‐Stage Mixed Integer Linear Programming (MILP) method to optimize spot intensity in IMPT with Lower Bound (LB) constraints. This method seeks to improve treatment planning efficiency and precision, overcoming limitations of existing strategies.MethodsOur approach evaluates prevalent IMPT spot selection strategies, identifying their limitations, especially concerning MU constraints. We integrated LB constraints into a MILP framework, using a novel three‐phase strategy for spot pool selection, to enhance performance over traditional heuristic methods and L1 + L∞ strategies. The method's efficacy was tested in eight study cases, using Dose‐Volume Histograms (DVHs), spot selection efficiency, and computation time analysis for benchmarking against established methods.ResultsThe proposed method showed superior performance in DVH quality, adhering to LB constraints while maintaining high‐quality treatment plans. It outperformed existing techniques in spot selection, reducing unnecessary spots and balancing precision with efficiency. Cases studies confirmed the method's effectiveness in producing clinically feasible plans with enhanced dose distributions and reduced hotspots, especially in cases with elevated LB constraints.ConclusionsOur Two‐Stage MILP strategy signifies a significant advancement in IMPT treatment planning. By incorporating LB constraints directly into the optimization process, it achieves superior plan quality and deliverability compared to current methods. This approach is particularly advantageous in clinical settings requiring minimum spot number and high MU LB constraints, offering the potential for improved patient outcomes through more precise and efficient radiation therapy plans.