Citations

AbstractBackgroundBoth dose and linear energy transfer (LET) could play a substantial role in adverse event (AE) initialization of cancer patients treated with pencil‐beam‐scanning (PBS) proton therapy. However, not all the voxels within the AE regions are directly induced from the dose and LET effect. It is important to study the synergistic effect of dose and LET in AE initialization by only including a subset of voxels that are dosimetrically important.PurposeTo perform exploratory investigation of the dose and LET effects upon AE initialization in PBS using seed spots analysis.MethodsA total of 113 head‐and‐neck (H&N) cancer patients receiving curative PBS were included. Among them, 20 patients experienced unanticipated CTCAEv4.0 grade ≥3 AEs (AE group) and 93 patients did not (control group). Within the AE group, 13 AE patients were included in the seed spot analysis to derive the descriptive features of AE initialization and the remaining 7 mandible osteoradionecrosis patients and 93 control patients were used to derive the feature‐based volume constraint of mandible osteoradionecrosis. The AE regions were contoured and the corresponding dose–LET volume histograms of AE regions were generated for all patients in the AE group. We selected high LET voxels (the highest 5% of each dose bin) with a range of moderate‐to‐high dose (≥∼40‐Gy relative biological effectiveness) as critical voxels. Critical voxels that were contiguous with each other were grouped into clusters. Each cluster was considered a potential independent seed spot for AE initialization. Seed spots were displayed in a 2D dose–LET plane based on their mean dose and LET to derive the descriptive features of AE initialization. A volume constraint of mandible osteoradionecrosis was then established based on the extracted features using a receiver operating characteristic curve.ResultsThe product of dose and LET (xBD) was found to be a descriptive feature of seed spots leading to AE initialization in this preliminary study. The derived xBD volume constraint for mandible osteoradionecrosis showed good performance with an area under curve of 0.87 (sensitivity of 0.714 and specificity of 0.807 in the leave‐one‐out cross‐validation) for the very limited patient data included in this study.ConclusionOur exploratory study showed that both dose and LET were observed to be important in AE initializations. The derived xBD volume constraint could predict mandible osteoradionecrosis reasonably well in the very limited H&N cancer patient data treated with PBS included in this study.

AbstractPurposeTo develop an online graphic processing unit (GPU)‐accelerated Monte Carlo‐based adaptive radiation therapy (ART) workflow for pencil beam scanning (PBS) proton therapy to address interfraction anatomical changes in patients treated with PBS.Methods and MaterialsA four‐step workflow was developed using our in‐house developed GPU‐accelerated Monte Carlo‐based treatment planning system to implement online Monte Carlo‐based ART for PBS. The first step conducts diffeomorphic demon‐based deformable image registration (DIR) to propagate contours on the initial planning CT (pCT) to the verification CT (vCT) to form a new structure set. The second step performs forward dose calculation of the initial plan on the vCT with the propagated contours after manual approval (possible modifications involved). The third step triggers a reoptimization of the plan depending on whether the verification dose meets the clinical requirements or not. A robust evaluation will be done for both the verification plan in the second step and the reopotimized plan in the third step. The fourth step involves a two‐stage (before and after delivery) patient‐specific quality assurance (PSQA) of the reoptimized plan. The before‐delivery PSQA is to compare the plan dose to the dose calculated using an independent fast open‐source Monte Carlo code, MCsquare. The after‐delivery PSQA is to compare the plan dose to the dose recalculated using the log file (spot MU, spot position, and spot energy) collected during the delivery. Jaccard index (JI), dice similarity coefficients (DSCs), and Hausdorff distance (HD) were used to assess the quality of the propagated contours in the first step. A commercial plan evaluation software, ClearCheck™, was integrated into the workflow to carry out efficient plan evaluation. 3D Gamma analysis was used during the fourth step to ensure the accuracy of the plan dose from reoptimization. Three patients with three different disease sites were chosen to evaluate the feasibility of the online ART workflow for PBS.ResultsFor all three patients, the propagated contours were found to have good volume conformance [JI (lowest‐highest: 0.833–0.983) and DSC (0.909–0.992)] but suboptimal boundary coincidence [HD (2.37–20.76 mm)] for organs‐at‐risk. The verification dose evaluated by ClearCheck™ showed significant degradation of the target coverage due to the interfractional anatomical changes. Reoptimization on the vCT resulted in great improvement of the plan quality to a clinically acceptable level. 3D Gamma analyses of PSQA confirmed the accuracy of the plan dose before delivery (mean Gamma index = 98.74% with a threshold of 2%/2 mm/10%), and after delivery based on the log files (mean Gamma index = 99.05% with a threshold of 2%/2 mm/10%). The average time cost for the complete execution of the workflow was around 858 s, excluding the time for manual intervention.ConclusionThe proposed online ART workflow for PBS was demonstrated to be efficient and effective by generating a reoptimized plan that significantly improved the plan quality.

Abstract Purpose To evaluate the accuracy of the RayStation Monte Carlo dose engine (RayStation MC) in modeling small‐field block apertures in proton pencil beam scanning. Furthermore, we evaluate the suitability of MCsquare as a second check for RayStation MC. Methods We have enhanced MCsquare to model block apertures. To test the accuracy of both RayStation MC and the newly enhanced MCsquare, we compare the dose predictions of each to in‐water dose measurements obtained using diode detectors and radiochromic film. Nine brass apertures with openings of 1, 2, 3, 4, and 5 cm and either 2 cm or 4 cm thickness were used in the irradiation of a water phantom. Two measurement setups were used, one with a range shifter and 119.7 MeV proton beam energy and the other with no range shifter and 147 MeV proton beam energy. To further test the validity of RayStation MC and MCsquare in modeling block apertures and to evaluate MCsquare as a second check tool, 10 small‐field (average target volume 8.3 cm 3 ) patient treatment plans were calculated by each dose engine followed by a statistical comparison. Results Comparing to the absolute dose measurements in water, RayStation MC differed by 1.2% ± 1.0% while MCsquare differed by −1.8% ± 3.7% in the plateau region of a pristine Bragg peak. Compared to the in‐water film measurements, RayStation MC and MCsquare both performed well with an average 2D‐3D gamma passing rate of 99.4% and 99.7% (3%/3 mm), respectively. A t ‐test comparing the agreement with the film measurements between RayStation MC and MCsquare suggested that the relative spatial dose distributions calculated by MCsquare and RayStation MC were statistically indistinguishable. Directly comparing the dose calculations between MCsquare and RayStation MC over 10 patients resulted in an average 3D‐3D gamma passing rates of 98.5% (3%/3 mm) and 94.1% (2%/2 mm), respectively. Conclusion The validity of RayStation MC algorithm for use with patient‐specific apertures has been expanded to include small apertures. MCsquare has been enhanced to model apertures and was found to be an adequate second check of RayStation MC in this scenario.

AbstractPurposeDue to the employment of quadratic programming using soft constraints to implement dose volume constraints and the “trial‐and‐error” procedure needed to achieve a clinically acceptable plan, conventional dose volume constraints (upper limit) are not adequately effective in controlling small and isolated hot spots in the dose/linear energy transfer (LET) distribution. Such hot spots can lead to adverse events. In order to mitigate the risk of brain necrosis, one of the most clinically significant adverse events in patients receiving intensity‐modulated proton therapy (IMPT) for base of skull (BOS) cancer, we propose per‐voxel constraints to minimize hot spots in LET‐guided robust optimization.Methods and materialsTen BOS cancer patients treated with IMPT were carefully selected by meeting one of the following conditions: (1) diagnosis of brain necrosis during follow‐up; and (2) considered high risk for brain necrosis by not meeting dose constraints to the brain. An optimizing structure (BrainOPT) and an evaluating structure (BrainROI) that both contained the aforementioned hot dose regions in the brain were generated for optimization and evaluation, respectively. Two plans were generated for every patient: one using conventional dose‐only robust optimization, the other using LET‐guided robust optimization. The impact of LET was integrated into the optimization via a term of extra biological dose (xBD). A novel optimization tool of per‐voxel constraints to control small and isolated hot spots in either the dose, LET, or combined (dose/LET) distribution was developed and used to minimize dose/LET hot spots of the selected structures. Indices from dose‐volume histogram (DVH) and xBD dose‐volume histogram (xBDVH) were used in the plan evaluation. A newly developed tool of the dose‐LET‐volume histogram (DLVH) was also adopted to illustrate the underlying mechanism. Wilcoxon signed‐rank test was used for statistical comparison of the DVH and xBDVH indices between the conventional dose‐only and the LET‐guided robustly optimized plans.ResultsPer‐voxel constraints effectively and efficiently minimized dose hot spots in both dose‐only and LET‐guided robust optimization and LET hot spots in LET‐guided robust optimization. Compared to the conventional dose‐only robust optimization, the LET‐guided robust optimization could generate plans with statistically lower xBD hot spots in BrainROI (VxBD,50 Gy[RBE], p = 0.009; VxBD,60 Gy[RBE], p = 0.025; xBD1cc, p = 0.017; xBD2cc, p = 0.022) with comparable dose coverage, dose hot spots in the target, and dose hot spots in BrainROI. DLVH analysis indicated that LET‐guided robust optimization could either reduce LET at the same dose level or redistribute high LET from high dose regions to low dose regions.ConclusionPer‐voxel constraint is a powerful tool to minimize dose/LET hot spots in IMPT. The LET‐guided robustly optimized plans outperformed the conventional dose‐only robustly optimized plans in terms of xBD hot spots control.

AbstractPurposeTo measure diode sensitivity degradation (DSD) induced by cumulative proton dose delivered to a commercial daily quality assurance (QA) device.MethodsAt our institution, six Daily QA 3 (DQA3, Sun Nuclear Corporation, Melbourne, FL, USA) devices have been used for daily proton pencil beam scanning QA in four proton gantry rooms over a span of 4 years. DQA3 diode counts were cross‐calibrated using a homogenous field with a known dose of 1 Gy. The DSD rate (ΔR%/100 Gy) was calculated using linear regression on time‐series plots of diode counts and an estimate of cumulative dose per year based on the cross‐calibration. The effect of DSD on daily QA spot position measurements was quantified by converting DSD to baseline spot position shift.ResultsThe average dose delivered to the four inner DQA3 diodes was 104 ± 5 Gy/year, and the rate of DSD was ‐5.1% ± 1.0/100 Gy with the exception of one DQA3 device that had a significantly higher rate of DSD (–12%/100 Gy). The R2s of the linear fit to time‐series plots were between 0.92 and 0.98. The DSD rates were not constant but decreases with accumulated doses. The four center diodes, which received 40% of the cumulative dose received by inner diodes, had a DSD rate of –7.2% ± 0.9/100 Gy. For our daily QA program, 1 year of DSD was equivalent to a 0.2 mm shift in spot position.ConclusionsThe DSD rate of DQA3 diodes determined by long‐term proton daily QA data was about –5%/100 Gy, which is more than 10 times greater than the reported DSD rate from photon irradiation. DQA3 diodes may be used for daily proton QA programs, provided that they are recalibrated at an appropriate frequency that should be determined specifically for different daily QA programs.

AbstractPurposeTo validate breast tissue expander metallic port (MP) models in a commercial treatment planning system (TPS) in proton pencil beam scanning (PBS) treatments for breast cancer patients with breast tissue expanders.Methods and materialsThree types of MPs taken out of a Mentor CPX4, a Natrelle 133, and a PMT Integra breast tissue expanders and a 650 cc saline filled Mentor CPX4 expander were placed on top of acrylic slabs, and scanned using a Siemens Somatom Definition AS Open RT CT scanner. Structure templates for each of the MPs were designed within Eclipse TPS. The CT numbers for the metallic parts were overridden to reflect measured or calculated relative proton stopping powers (RPSPs). Mock targets were contoured in acrylic to represent postmastectomy chest‐wall radiation therapy (PMRT) targets. Plans with different beam incident angles were optimized using the Eclipse TPS to deliver uniform prescription dose to the target using Hitachi Probeat‐V PBS beams. Eclipse calculated doses and an in‐house Monte Carlo (MC) code calculated doses were compared to the measured Gafchromic EBT3 film doses in acrylic.ResultsTPS/MC and film dose comparison results showed that (1) 3%/2 mm/10% threshold Gamma pass rates were better than 90.8% in the acrylic target region for all plans; (2) comparing TPS and film doses for the individual beam plans in the MP dose shadow areas, the area with dose difference above 5% ([ΔA] 5%) ranged from 1.1 to 5.0 cm2, and the maximum dose difference ([ΔD] 0.01 cm2) ranged from 12.5% to 25.0%; (3) comparing MC and film doses for the individual beam plans in the MP dose shadow areas, the (ΔA) 5% varied from 1.1 to 2.9 cm2 and (ΔD) 0.01 cm2 varied from 8.5% to 24.2%; (4) for a plan composed of three individual beams treating through the Mentor CPX4 expander, the TPS (ΔA) 5% was less than 0.13 cm2, and the (ΔD) 0.01 cm2 was less than 6% in the MP dose shadow areas.ConclusionsIt is feasible to treat patients with tissue expanders using multiple PBS beams using a structure template with CT number overridden to represent the measured/calculated RPSP for MPs for PBS treatment planning. MC dose was more accurate than analytical dose in the areas with high dose gradient caused by the density heterogeneity of the breast tissue expander MPs.

Abstract Purpose The multiple energy extraction (MEE) delivery technique for synchrotron‐based proton delivery systems reduces beam delivery time by decelerating the beam multiple times during one accelerator spill, but this might cause additional plan quality degradation due to intrafractional motion. We seek to determine whether MEE causes significantly different plan quality degradation compared to single energy extraction (SEE) for lung cancer treatments due to the interplay effect. Methods Ten lung cancer patients treated with IMPT at our institution were nonrandomly sampled based on a representative range of tumor motion amplitudes, tumor volumes, and respiratory periods. Dose‐volume histogram (DVH) indices from single‐fraction SEE and MEE four‐dimensional (4D) dynamic dose distributions were compared using the Wilcoxon signed‐rank test. Distributions of monitor units (MU) to breathing phases were investigated for features associated with plan quality degradation. SEE and MEE DVH indices were compared in fractionated deliveries of the worst‐case patient treatment scenario to evaluate the impact of fractionation. Results There were no clinically significant differences in target mean dose, target dose conformity, or dose to organs‐at‐risk between SEE and MEE in single‐fraction delivery. Three patients had significantly worse dose homogeneity with MEE compared to SEE (single‐fraction mean D 5% ‐D 95% increased by up to 9.6% of prescription dose), and plots of MU distribution to breathing phases showed synchronization patterns with MEE but not SEE. However, after 30 fractions the patient in the worst‐case scenario had clinically acceptable target dose homogeneity and coverage with MEE (mean D 5% –D 95% increased by 1% compared to SEE). Conclusions For some patients with breathing periods close to the mean spill duration, MEE resulted in significantly worse single‐fraction target dose homogeneity compared to SEE due to the interplay effect. However, this was mitigated by fractionation, and target dose homogeneity and coverage were clinically acceptable after 30 fractions with MEE.

Purpose To compare the dosimetric performances of small‐spot three‐dimensional (3D) and four‐dimensional (4D) robustly optimized intensity‐modulated proton (IMPT) plans in the presence of uncertainties and interplay effect simultaneously for distal esophageal carcinoma. Method and Materials Thirteen (13) patients were selected and re‐planned with small‐spot ( ~ 2–6 mm) 3D and 4D robust optimization in IMPT, respectively. The internal clinical target volumes (CTV high3d , CTV low3d ) were used in 3D robust optimization. Different CTVs (CTV high4d , CTV low4d ) were generated by subtracting an inner margin of the motion amplitudes in three cardinal directions from the internal CTVs and used in 4D robust optimization. All patients were prescribed the same dose to CTVs (50 Gy[RBE] for CTV high3d /CTV high4d and 45 Gy[RBE] for CTV low3d /CTV low4d ). Dose–volume histogram (DVH) indices were calculated to assess plan quality. Comprehensive plan robustness evaluations that consisted of 300 perturbed scenarios (10 different motion patterns to consider irregular motion (sampled from a Gaussian distribution) and 30 different uncertainties scenarios (sampled from a 4D uniform distribution) combined), were performed to quantify robustness to uncertainties and interplay effect simultaneously. Wilcoxon signed‐rank test was used for statistical analysis. Results Compared to 3D robustly optimized plans, 4D robustly optimized plans had statistically improved target coverage and better sparing of lungs and heart (heart D mean , P = 0.001; heart V 30Gy[RBE] , P = 0.001) in the nominal scenario. 4D robustly optimized plans had better robustness in target dose coverage (CTV high4d V 100% , P = 0.002) and the protection of lungs and heart (heart D mean , P = 0.001; heart V 30Gy[RBE] , P = 0.001) when uncertainties and interplay effect were considered simultaneously. Conclusions Even with small spots in IMPT, 4D robust optimization outperformed 3D robust optimization in terms of normal tissue protection and robustness to uncertainties and interplay effect simultaneously. Our findings support the use of 4D robust optimization to treat distal esophageal carcinoma with small spots in IMPT.

Abstract Purpose To evaluate dosimetric properties of intensity‐modulated proton therapy (IMPT) for simulated treatment planning in patients with atrial fibrillation (AF) targeting left atrial‐pulmonary vein junction (LA‐PVJ), in comparison with volumetric‐modulated arc therapy (VMAT) and helical tomotherapy (TOMO). Methods Ten thoracic 4D‐CT scans with respiratory motion and one with cardiac motion were used for the study. Ten respiratory 4D‐CTs were planned with VMAT, TOMO, and IMPT for simulated AF. Targets at the LA‐PVJ were defined as wide‐area circumferential ablation line. A single fraction of 25 Gy was prescribed to all plans. The interplay effects from cardiac motion were evaluated based on the cardiac 4D‐CT scan. Dose‐volume histograms (DVHs) of the ITV and normal tissues were compared. Statistical analysis was evaluated via one‐way Repeated‐Measures ANOVA and Friedman’s test with Bonferroni’s multiple comparisons test. Results The median volume of ITV was 8.72cc. All plans had adequate target coverage (V 23.75Gy ≥ 99%). Compared with VMAT and TOMO, IMPT resulted in significantly lower dose of most normal tissues. For VMAT, TOMO, and IMPT plans, D mean of the whole heart was 5.52 ± 0.90 Gy, 5.89 ± 0.78 Gy, and 3.01 ± 0.57 Gy ( P < 0.001), mean dose of pericardium was 4.74 ± 0.76 Gy, 4.98 ± 0.62 Gy, and 2.59 ± 0.44 Gy ( P < 0.001), and D 0.03cc of left circumflex artery (LCX) was 13.96 ± 5.45 Gy, 14.34 ± 5.91 Gy, and 8.43 ± 7.24 Gy ( P < 0.001), respectively. However, no significant advantage for one technique over the others was observed when examining the D 0.03cc of esophagus and main bronchi. Conclusions IMPT targeting LA‐PVJ for patients with AF has high potential to reduce dose to surrounding tissues compared to VMAT or TOMO. Motion mitigation techniques are critical for a particle‐therapy approach.