Intensity Modulated Carbon Ion Therapy Research Articles

Validation of robust radiobiological optimization algorithms based on the mixed beam model for intensity-modulated carbon-ion therapy.

Currently, treatment planning systems (TPSs) that can compute the intensities of intensity-modulated carbon-ion therapy (IMCT) using scanned carbon-ion beams are limited. In the present study, the computational efficacy of the newly designed IMCT algorithms was analyzed for the first time based on the mixed beam model with respect to the physical and biological doses; moreover, the validity and effectiveness of the robust radiobiological optimization were verified. A dose calculation engine was independently generated to validate a clinical dose determined in the TPS. A biological assay was performed using the HSGc-C5 cell line to validate the calculated surviving fraction (SF). Both spot control (SC) and voxel-wise worst-case scenario (WC) algorithms were employed for robust radiobiological optimization followed by their application in a Radiation Therapy Oncology Group benchmark phantom under homogeneous and heterogeneous conditions and a clinical case for range and position errors. Importantly, for the first time, both SC and WC algorithms were implemented in the integrated TPS platform that can compute the intensities of IMCT using scanned carbon-ion beams for robust radiobiological optimization. For assessing the robustness, the difference between the maximum and minimum values of a dose-volume histogram index in the examined error scenarios was considered as a robustness index. The relative biological effectiveness (RBE) determined by the independent dose calculation engine exhibited a -0.6% difference compared with the RBE defined by the TPS at the isocenter, whereas the measured and the calculated SF were similar. Regardless of the objects, compared with the conventional IMCT, the robust radiobiological optimization enhanced the sensitivity of the examined error scenarios by up to 19% for the robustness index. The computational efficacy of the novel IMCT algorithms was verified according to the mixed beam model with respect to the physical and biological doses. The robust radiobiological optimizations lowered the impact of range and position uncertainties considerably in the examined scenarios. The robustness of the WC algorithm was more enhanced compared with that of the SC algorithm. Nevertheless, the SC algorithm can be used as an alternative to the WC IMCT algorithm with respect to the computational cost.

Fixed Beamline Optimization for Intensity Modulated Carbon-Ion Therapy.

A major obstacle for the adoption of heavy ion therapy is the cost and technical difficulties to construct and maintain a rotational gantry. Many heavy ion treatment facilities instead choose to construct fixed beamlines as a compromise, which we propose to mitigate with optimized treatment couch angle. We formulate the integrated beam orientation and scanning spot optimization problem as a quadratic cost function with a group sparsity regularization term. The optimization problem is efficiently solved using fast iterative shrinkage-thresholding algorithm (FISTA). To test the method, we created the fixed beamline plans with couch rotation (FBCR) and without couch rotation (FB) for intensity modulated carbon-ion therapy (IMCT) and compared with the ideal scenario where both the couch and gantry have 360 degrees of freedom (GCR). FB, FBCR, and GCR IMCT plans were compared for ten pancreas cases. The FBCR plans show comparable PTV coverage and OAR doses for each pancreas case. In conclusion, the dosimetric limitation of fixed beams in heavy ion radiotherapy may be largely mitigated with integrated beam orientation optimization of the couch rotation.

Ion therapy guideline (Version 2020)

Ion therapy guideline (Version 2020)

Robust treatment planning with 4D intensity modulated carbon ion therapy for multiple targets in stage IV non-small cell lung cancer

Intensity modulated particle therapy (IMPT) with carbon ions can generate highly conformal treatment plans; however, IMPT is limited in robustness against range and positioning uncertainty. This is particularly true for moving targets, even though all motion states of a 4DCT are considered in 4D-IMPT. Here, we expand 4D-IMPT to include robust non-linear RBE-weighted optimization to explore its potential in improving plan robustness and sparing critical organs. In this study, robust 4D-optimization—based on worst-case optimization on 9 scenarios—was compared to conventional 4D-optimization with PTV margins using 4D dose calculation and robustness analysis for 21 uncertainty scenarios. Slice-by-slice rescanning was used for motion mitigation. Both 4D-optimization strategies were tested on a cohort of 8 multi-lesion lung cancer patients with the goal of prioritizing OAR sparing in a hypofractionated treatment plan. Planning objectives were to keep the OAR volume doses below corresponding limits while simultaneously achieve CTV coverage with D95% ≥ 95 %. For the conventional plans, average D95% was at 98.7% which fulfilled the target objective in 83.2% of scenarios. For the robust plans, average D95% was reduced to 97.6% which still fulfilled the target objective in 80.7% of cases, but led to significantly improved overall OAR sparing: Volume doses were below the limits in 96.2% of cases for the conventional and 99.5% for the robust plans. When considering the particularly critical smaller airways only, fulfillment rates could be increased from 76.2% to 96% for the robust plans. This study has shown that plan robustness of 4D-IMPT could be improved by using robust 4D-optimization, offering greater control over uncertainties in the actual delivered dose. In some cases, this required sacrificing target coverage for the benefit of better OAR sparing.

A track repeating algorithm for intensity modulated carbon ion therapy

The fast dose calculator (FDC), a track repeating Monte Carlo (MC) algorithm was initially developed for proton therapy. The validation for proton therapy has been demonstrated in a previous work. In this work we presented the extension of FDC to the calculation of dose distributions for ions, particularly for carbon. Moreover the code algorithm is validated by comparing 3D dose distributions and dose volume histograms (DVH) calculated by FDC with Geant4. A total of 19 patients were employed, including three patients of prostate, five of brain, three of head and neck, four of lung and four of spine. We used a gamma-index technique to analyze dose distributions and we performed a dosimetric analysis for DVHs, a more direct and informative quantity for planning system assessment. The gamma-index passing rates of all patients discussed in this paper are above 90% with the criterion 1%/1 mm, above 98% with the criterion 2%/2 mm and over 99.9% with the criterion 3%/3 mm. The root mean square (RMS) of percent difference of dosimetric indices D02, D05, D50, D95 and D98 are 0.75%, 0.70%, 0.79%, 0.83% and 0.76%. All the differences are within clinically accepted norms.

Different optimizing strategies based on spot scanning carbon ion therapy for non-small cell lung cancer

Objective To evaluate the dose variation of target coverage and organs at risk (OARs) among four planning strategies using spot-scanning carbon-ion radiotherapy for non-small cell lung cancer (NSCLC). Methods Ten NSCLC patients utilizing gating motion control were selected to receive dose calculation over multiple acquired 4DCT images. Four optimizing strategies consisted of intensity-modulated carbon-ion therapy (IMCT-NoAS), IMCT combined with internal gross tumor volume (IGTV) assigned muscle (IMCT-ASM), single beam optimization (SBO)(SBO-NoAS) and SBO combined with IGTV assigned muscle (SBO-ASM). The initial plan was re-calculated after the 4DCT data were reviewed and then compared with the initial plan in the dosimetry. Results For re-calculation plans with two reviewing CTs, all four strategies yielded similar planning target volume (PTV) coverage. Merely IMCT-NoAS strategy presented with relatively significant variations in dose distribution. Dose variation for OARs between initial and re-calculated plans: for all four strategies, V20 of ipsilateral lung was increased by approximately 2.0 Gy (relative biological effective dose, RBE), V30 of heart was increased by approximately 1.0 Gy (RBE) for both IGTV assigned muscle strategies, whereas decreased by approximately 0.2 Gy (RBE) for both IGTV non-assigned muscle strategies. The maximum dose of spinal cord was changed by 2.5 Gy (RBE). Conclusions Carbon-ion radiotherapy is sensitive to the anatomic motion within the tumors along the beam path. When the tumor motion along the head-foot (H-F) direction exceeds 8 mm, SBO-ASM strategy provides better dose coverage of target. Strategies with IGTV assignment may result in dose overshoot to a position deeper than the initial planning dose distribution. Key words: Lung neoplasms/intensity-modulated carbon-ion therapy; Single beam optimization; Spot scanning carbon-ion; Plan evaluation

Analytical probabilistic modeling of RBE-weighted dose for ion therapy

Particle therapy is especially prone to uncertainties. This issue is usually addressed with uncertainty quantification and minimization techniques based on scenario sampling. For proton therapy, however, it was recently shown that it is also possible to use closed-form computations based on analytical probabilistic modeling (APM) for this purpose. APM yields unique features compared to sampling-based approaches, motivating further research in this context.This paper demonstrates the application of APM for intensity-modulated carbon ion therapy to quantify the influence of setup and range uncertainties on the RBE-weighted dose. In particular, we derive analytical forms for the nonlinear computations of the expectation value and variance of the RBE-weighted dose by propagating linearly correlated Gaussian input uncertainties through a pencil beam dose calculation algorithm. Both exact and approximation formulas are presented for the expectation value and variance of the RBE-weighted dose and are subsequently studied in-depth for a one-dimensional carbon ion spread-out Bragg peak. With V and B being the number of voxels and pencil beams, respectively, the proposed approximations induce only a marginal loss of accuracy while lowering the computational complexity from order to for the expectation value and from to for the variance of the RBE-weighted dose. Moreover, we evaluated the approximated calculation of the expectation value and standard deviation of the RBE-weighted dose in combination with a probabilistic effect-based optimization on three patient cases considering carbon ions as radiation modality against sampled references. The resulting global γ-pass rates (2 mm,2%) are 99.15% for the expectation value and 94.95% for the standard deviation of the RBE-weighted dose, respectively. We applied the derived analytical model to carbon ion treatment planning, although the concept is in general applicable to other ion species considering a variable RBE.

Full Monte Carlo–Based Biologic Treatment Plan Optimization System for Intensity Modulated Carbon Ion Therapy on Graphics Processing Unit

One of the major benefits of carbon ion therapy is enhanced biological effectiveness at the Bragg peak region. For intensity modulated carbon ion therapy (IMCT), it is desirable to use Monte Carlo (MC) methods to compute the properties of each pencil beam spot for treatment planning, because of their accuracy in modeling physics processes and estimating biological effects. We previously developed goCMC, a graphics processing unit (GPU)-oriented MC engine for carbon ion therapy. The purpose of the present study was to build a biological treatment plan optimization system using goCMC. The repair-misrepair-fixation model was implemented to compute the spatial distribution of linear-quadratic model parameters for each spot. A treatment plan optimization module was developed to minimize the difference between the prescribed and actual biological effect. We used a gradient-based algorithm to solve the optimization problem. The system was embedded in the Varian Eclipse treatment planning system under a client-server architecture to achieve a user-friendly planning environment. We tested the system with a 1-dimensional homogeneous water case and 3 3-dimensional patient cases. Our system generated treatment plans with biological spread-out Bragg peaks covering the targeted regions and sparing critical structures. Using 4 NVidia GTX 1080 GPUs, the total computation time, including spot simulation, optimization, and final dose calculation, was 0.6hour for the prostate case (8282 spots), 0.2hour for the pancreas case (3795 spots), and 0.3hour for the brain case (6724 spots). The computation time was dominated by MC spot simulation. We built a biological treatment plan optimization system for IMCT that performs simulations using a fast MC engine, goCMC. To the best of our knowledge, this is the first time that full MC-based IMCT inverse planning has been achieved in a clinically viable time frame.

Comparison of photon volumetric modulated arc therapy, intensity-modulated proton therapy, and intensity-modulated carbon ion therapy for delivery of hypo-fractionated thoracic radiotherapy

PurposeThe aim of the present study was to compare the dose distribution generated from photon volumetric modulated arc therapy (VMAT), intensity modulated proton therapy (IMPT), and intensity modulated carbon ion therapy (IMCIT) in the delivery of hypo-fractionated thoracic radiotherapy.Methods and materialsTen selected patients who underwent thoracic particle therapy between 2015 and 2016 were re-planned to receive a relative biological effectiveness (RBE) weighted dose of 60 Gy (i.e., GyE) in 15 fractions delivered with VMAT, IMPT, or IMCIT with the same optimization criteria. Treatment plans were then compared.ResultsThere were no significant differences in target volume dose coverage or dose conformity, except improved D95 was found with IMCIT compared with VMAT (p = 0.01), and IMCIT was significantly better than IMPT in all target volume dose parameters. Particle therapy led to more prominent lung sparing at low doses, and this result was most prominent with IMCIT (p < 0.05). Improved sparing of other thoracic organs at risk (OARs) was observed with particle therapy, and IMCIT further lowered the D1cc and D5cc for major blood vessels, as compared with IMPT (p = 0.01).ConclusionAlthough it was comparable to VMAT, IMCIT led to significantly better tumor target dose coverage and conformity than did IMPT. Particle therapy, compared with VMAT, improved thoracic OAR sparing. IMCIT, compared with IMPT, may further improve normal lung and major blood vessel sparing under limited respiratory motion.

Photon Volumetric Modulated Arc Therapy Versus Intensity Modulated Proton Therapy Versus Intensity Modulated Carbon Ion Therapy in the Delivery of Hypofractionated Thoracic Radiation Therapy

Photon Volumetric Modulated Arc Therapy Versus Intensity Modulated Proton Therapy Versus Intensity Modulated Carbon Ion Therapy in the Delivery of Hypofractionated Thoracic Radiation Therapy

Dosimetric Considerations to Determine the Optimal Technique for Localized Prostate Cancer Among External Photon, Proton, or Carbon-Ion Therapy and High-Dose-Rate or Low-Dose-Rate Brachytherapy

To assess the dosimetric differences among volumetric modulated arc therapy (VMAT), scanned proton therapy (intensity-modulated proton therapy, IMPT), scanned carbon-ion therapy (intensity-modulated carbon-ion therapy, IMIT), and low-dose-rate (LDR) and high-dose-rate (HDR) brachytherapy (BT) treatment of localized prostate cancer. Ten patients were considered for this planning study. For external beam radiation therapy (EBRT), planning target volume was created by adding a margin of 5 mm (lateral/anterior-posterior) and 8 mm (superior-inferior) to the clinical target volume. Bladder wall (BW), rectal wall (RW), femoral heads, urethra, and pelvic tissue were considered as organs at risk. For VMAT and IMPT, 78 Gy(relative biological effectiveness, RBE)/2 Gy were prescribed. The IMIT was based on 66 Gy(RBE)/20 fractions. The clinical target volume planning aims for HDR-BT ((192)Ir) and LDR-BT ((125)I) were D(90%) ≥34 Gy in 8.5 Gy per fraction and D(90%) ≥145 Gy. Both physical and RBE-weighted dose distributions for protons and carbon-ions were converted to dose distributions based on 2-Gy(IsoE) fractions. From these dose distributions various dose and dose-volume parameters were extracted. Rectal wall exposure 30-70 Gy(IsoE) was reduced for IMIT, LDR-BT, and HDR-BT when compared with VMAT and IMPT. The high-dose region of the BW dose-volume histogram above 50 Gy(IsoE) of IMPT resembled the VMAT shape, whereas all other techniques showed a significantly lower high-dose region. For all 3 EBRT techniques similar urethra D(mean) around 74 Gy(IsoE) were obtained. The LDR-BT results were approximately 30 Gy(IsoE) higher, HDR-BT 10 Gy(IsoE) lower. Normal tissue and femoral head sparing was best with BT. Despite the different EBRT prescription and fractionation schemes, the high-dose regions of BW and RW expressed in Gy(IsoE) were on the same order of magnitude. Brachytherapy techniques were clearly superior in terms of BW, RW, and normal tissue sparing, with lowest values for HDR-BT.