Furutani, Keith M

Abstract Objective . To investigate the impact of scan path optimization on the dose accuracy and beam delivery time (BDT) of proton pencil beam scanning in the dose-driven continuous scanning (DDCS). Approach . A diverse set of six clinical plans, representing various spot patterns and treatment sites, was used to evaluate the effectiveness of scan time optimization and scan length optimization. The DDCS dose discrepancy and BDT with optimized scan paths was compared to the default serpentine scan path. Main results . Both scan time optimization and scan path optimization were able to reduce the DDCS dose discrepancy compared to the default serpentine scan path. All plans, except for the layer repainting lung plan, achieved a 2%/2 mm gamma pass rate of over 99% and less than 1% PTV DVH root mean square error (RMSE) through scan path optimization. In the case of the layer repainting lung plan, when compared to the default serpentine scan path, the 2%/2 mm gamma pass rate showed improvements from 91.3% to 93.1% and 95.8%, while the PTV DVH RMSE decreased from 2.1% to 1.7% and 1.1% for scan time optimization and scan length optimization, respectively. Although scan time optimization resulted in shorter total scan times for all plans compared to the default scan path and scan length optimization tended to have longer total scan times. However, due to the short total scan times and their minimal contribution to the total BDT, the impact of scan path optimization on the total BDT was practically negligible. Significance . Both scan time optimization and scan length optimization proved to be effective in minimizing DDCS dose discrepancy. No definitive winner can be determined between these two optimization approaches. Both scan time and scan length optimization had minimal effect on the total BDT.

Pencil beam scanning delivered with continuous scanning has several advantages over conventional discrete spot scanning. Such advantages include improved beam delivery efficiency and reduced beam delivery time. However, a move dose is delivered between consecutive spots with continuous scanning, and current treatment planning systems do not take this into account. Therefore, continuous scanning and discrete spot plans have an inherent dose discrepancy. Using the operating parameters of the state-of-the-art particle therapy system, we conducted a proof-of-concept study in which we systematically generated 28 plans for cubic targets with different combinations of plan parameters and simulated the dose discrepancies between continuous scanning and a planned one. A nomograph to guide the selection of plan parameters was developed to reduce the dose discrepancy. The effectiveness of the nomograph was evaluated with two clinical cases (one prostate and one liver). Plans with parameters guided by the nomograph decreased dose discrepancy than those used standard plan parameters. Specifically, the 2%/2 mm gamma passing rate increased from 96.3% to 100% for the prostate case and from 97.8% to 99.7% for the liver case. The CTV DVH root mean square error decreased from 2.2% to 0.2% for the prostate case and from 1.8% to 0.9% for the liver case. The decreased dose discrepancy may allow the relaxing of the delivery constraint for some cases, leading to greater benefits in continuous scanning. Further investigation is warranted.

Abstract Objective . UNet-based deep-learning (DL) architectures are promising dose engines for traditional linear accelerator (Linac) models. Current UNet-based engines, however, were designed differently with various strategies, making it challenging to fairly compare the results from different studies. The objective of this study is to thoroughly evaluate the performance of UNet-based models on magnetic-resonance (MR)-Linac-based intensity-modulated radiation therapy (IMRT) dose calculations. Approach . The UNet-based models, including the standard-UNet, cascaded-UNet, dense-dilated-UNet, residual-UNet, HD-UNet, and attention-aware-UNet, were implemented. The model input is patient CT and IMRT field dose in water, and the output is patient dose calculated by DL model. The reference dose was calculated by the Monaco Monte Carlo module. Twenty training and ten test cases of prostate patients were included. The accuracy of the DL-calculated doses was measured using gamma analysis, and the calculation efficiency was evaluated by inference time. Results . All the studied models effectively corrected low-accuracy doses in water to high-accuracy patient doses in a magnetic field. The gamma passing rates between reference and DL-calculated doses were over 86% (1%/1 mm), 98% (2%/2 mm), and 99% (3%/3 mm) for all the models. The inference times ranged from 0.03 (graphics processing unit) to 7.5 (central processing unit) seconds. Each model demonstrated different strengths in calculation accuracy and efficiency; Res-UNet achieved the highest accuracy, HD-UNet offered high accuracy with the fewest parameters but the longest inference, dense-dilated-UNet was consistently accurate regardless of model levels, standard-UNet had the shortest inference but relatively lower accuracy, and the others showed average performance. Therefore, the best-performing model would depend on the specific clinical needs and available computational resources. Significance . The feasibility of using common UNet-based models for MR-Linac-based dose calculations has been explored in this study. By using the same model input type, patient training data, and computing environment, a fair assessment of the models’ performance was present.

Abstract Background In order to compute the relative biological effectiveness (RBE) of ion radiation therapy with the Mayo Clinic Florida microdosimetric kinetic model (MCF MKM), it is necessary to process entire microdosimetric distributions. Therefore, a posteriori RBE recalculations (i.e., for a different cell line or another biological endpoint) would require whole spectral information. It is currently not practical to compute and store all this data for each clinical voxel. Purpose To develop a methodology that allows to store a limited amount of physical information without losing accuracy in the RBE calculations nor the possibility of a posteriori RBE recalculations. Methods Computer simulations for four monoenergetic 12 C ion beams and a 12 C ion spread‐out Bragg peak (SOBP) were performed to assess lineal energy distributions as a function of the depth within a water phantom. These distributions were used in combination with the MCF MKM to compute the in vitro clonogenic survival RBE for human salivary gland tumor cells (HSG cell line) and human skin fibroblasts (NB1RGB cell line). The RBE values were also calculated with a new abridged microdosimetric distribution methodology (AMDM) and compared with the reference RBE calculations using the entire distributions. Results The maximum relative deviation between the RBE values computed using the entire distributions and the AMDM was 0.61% (monoenergetic beams) and 0.49% (SOBP) for the HSG cell line, while 0.45% (monoenergetic beams) and 0.26% (SOBP) for the NB1RGB cell line. Conclusion The excellent agreement between the RBE values computed using the entire lineal energy distributions and the AMDM represents a milestone for the clinical implementation of the MCF MKM.

Abstract The Mayo Clinic Florida Integrated Oncology Building will be the home of the first spot-scanning only carbon/proton hybrid therapy system by Hitachi, Ltd. It will provide proton beams up to kinetic energies of 230 MeV and carbon beams up to 430 MeV n −1 for clinical deployment. To provide adequate radiation protection, the Geant4 (v10.6) Monte Carlo toolkit was utilized to quantify the ambient dose equivalent at a 10 mm depth ( H *(10)) for photons and neutrons. To perform accurate calculations of the ambient dose equivalent, three-dimensional computer-aided design files of the entire planned facility were imported into Geant4, as well as certain particle system components such as the bending magnets, fast Faraday cup, and gantry. Particle fluence was scored using 60 cm diameter spheres, which were strategically placed throughout areas of interests. Analytical calculations were performed as first-pass design checks. Major shielding slabs were optimized using Geant4 simulations iteratively, with more than 20 alternative designs evaluated within Geant4. The 430 MeV n −1 carbon beams played the most significant role in concrete thickness Requirements. The primary wall thickness for the carbon fixed beam room is 4 meters. The presence of the proton gantry structure in the simulation caused the ambient dose equivalent to increase by around 67% at the maze entrance, but a decrease in the high energy beam transport corridor. All shielding primary and secondary goals for clinical operations were met per state regulation and national guidelines.

The computation of the relative biological effectiveness (RBE) is a fundamental step in the planning of cancer radiotherapy treatments with accelerated ions. Numerical parameters derived analyzing the dose response of the chosen cell line after irradiation to photons (i.e., α and β, namely the linear and quadratic terms of the linear-quadratic model of cell survival) are generally used as input to biophysical models to predict the ion RBE. The α/β ratio for the photon exposure is generally regarded as an indicator of cell radiosensitivity. However, previous studies suggest that α/β might not be a sufficient parameter to model the RBE of relatively high linear energy transfer (LET) radiation such as carbon ions. For a fixed α/β, the effect of the absolute values of α and β on the computed RBE is underexplored. Furthermore, since α and β are anticorrelated during the fit of the photon-exposed in vitro survival data, different linear-quadratic fits could produce different sets of α and β, thus affecting the RBE calculations. This article reports the combined effect of the α/β ratio and the absolute values α and β on the RBE computed with the Mayo Clinic Florida microdosimetric kinetic model (MCF MKM) for 12C ions of different LET. Furthermore, we introduce a theory-based strategy to potentially mitigate the anticorrelation between α and β during the fit of the photon dose-response biological data.

The relative biological effectiveness (RBE) calculations used during the planning of ion therapy treatments are generally based on the microdosimetric kinetic model (MKM) and the local effect model (LEM). The Mayo Clinic Florida MKM (MCF MKM) was recently developed to overcome the limitations of previous MKMs in reproducing the biological data and to eliminate the need for ion-exposed in vitro data as input for the model calculations. Since we are considering to implement the MCF MKM in clinic, this article presents (a) an extensive benchmark of the MCF MKM predictions against corresponding in vitro clonogenic survival data for 4 rodent and 10 cell lines exposed to ions from 1H to 238U, and (b) a systematic comparison with published results of the latest version of the LEM (LEM IV). Additionally, we introduce a novel approach to derive an approximate value of the MCF MKM model parameters by knowing only the animal species and the mean number of chromosomes. The overall good agreement between MCF MKM predictions and in vitro data suggests the MCF MKM can be reliably used for the RBE calculations. In most cases, a reasonable agreement was found between the MCF MKM and the LEM IV.

Abstract Objective . To investigate similarities and differences in the formalism, processing, and the results of relative biological effectiveness (RBE) calculations with a biological weighting function (BWF), the microdosimetric kinetic model (MKM) and subsequent modifications (non-Poisson MKM, modified MKM). This includes: (a) the extension of the V79-RBE 10% BWF to model the RBE for other clonogenic survival levels; (b) a novel implementation of MKMs as weighting functions; (c) a benchmark against Chinese Hamster lung fibroblast (V79) in vitro data; (d) a study on the effect of pre- or post- processing the average biophysical quantities used for the RBE calculations; (e) a possible modification of the modified MKM parameters to improve the model accuracy at high linear energy transfer (LET). Methodology . Lineal energy spectra were simulated for two spherical targets (diameter = 0.464 or 1.0 μ m) using PHITS for 1 H, 4 He, 12 C, 20 Ne, 40 Ar, 56 Fe and 132 Xe ions. The results of the in silico calculations were compared with published in vitro data. Main results . All models appear to underestimate the RBE α of hydrogen ions. All MKMs generally overestimate the RBE 50% , RBE 10% and RBE 1% for ions with an LET greater than ∼200 keV μ m −1 . This overestimation is greater for small surviving fractions and is likely due to the assumption of a radiation-independent quadratic term of clonogenic survival ( ß ). The overall RBE trends seem to be best described by the novel ‘post-processing average’ implementation of the non-Poisson MKM. In case of calculations with the non-Poisson MKM, pre- or post- processing the average biophysical quantities affects the computed RBE values significantly. Significance . This study presents a systematic analysis of the formalism and results of widely used microdosimetric models of clonogenic survival for ions relevant for cancer particle therapy and space radiation protection. Points for improvements were highlighted and will contribute to the development of upgraded biophysical models.

Stereotactic radiosurgery (SRS) involves the delivery of a highly conformal ablative dose of radiation to both benign and malignant targets. This has traditionally been accomplished in a single fraction; however, fractionated approaches involving five or fewer treatments have been delivered for larger lesions, as well as lesions in close proximity to radiosensitive structures. The clinical utilization of SRS has overwhelmingly involved photon-based sources via dedicated radiosurgery platforms (e.g., Gamma Knife® and Cyberknife®) or specialized linear accelerators. While photon-based methods have been shown to be highly effective, advancements are sought for improved dose precision, treatment duration, and radiobiologic effect, among others, particularly in the setting of repeat irradiation. Particle-based techniques (e.g., protons and carbon ions) may improve many of these shortcomings. Specifically, the presence of a Bragg Peak with particle therapy at target depth allows for marked minimization of distal dose delivery, thus mitigating the risk of toxicity to organs at risk. Carbon ions also exhibit a higher linear energy transfer than photons and protons, allowing for greater relative biological effectiveness. While the data are limited, utilization of proton radiosurgery in the setting of brain metastases has been shown to demonstrate 1-year local control rates >90%, which are comparable to that of photon-based radiosurgery. Prospective studies are needed to further validate the safety and efficacy of this treatment modality. We aim to provide a comprehensive overview of clinical evidence in the use of particle therapy-based radiosurgery.