Plan Robustness Research Articles

Comparison of Treatment Plan Robustness Between Dynamic Trajectory Radiotherapy and VMAT

Dynamic trajectory radiotherapy (DTRT) was recently developed in our institute for a non-coplanar arc delivery technique utilizing additionally to VMAT table and collimator rotations. Recent studies on the dosimetric performance suggest a substantial improvement in organs at risks (OARs) sparing when compared with conventional C-arm VMAT delivery. In this work, the treatment plan robustness for DTRT and VMAT is compared with respect to setup errors. Robustness was assessed by generating DVH bands as well as so-called robustness maps (rMaps) by means of a false/true evaluation of certain dosimetric user-defined limits considered acceptable for selected structures with respect to a range of setup errors considered. The robustness was assessed for a clinically motivated head and neck, a lung and a prostate case. For each case systematic translational (ranging from -6 to 6 mm) and rotational (ranging from -3° to 3°) setup error scenarios were defined and for each scenario Monte Carlo dose calculations for the DTRT and VMAT treatment plans were performed. The resulting dose distributions served as input for the robustness assessment determining DVH bands and rMaps. The acceptance conditions for the OARs were defined such that the mean doses for parallel OARs and near maximum doses for serial OARs do not deteriorate more than 2 Gy. For the CTV V95% a deterioration of 1% was allowed in order to pass the robustness criteria. For the isotropic CTV to PTV margins 4 mm, 5 mm and 10 mm were applied for the head and neck, the prostate and the lung case, respectively, except for the dorsal direction of the prostate case, where 3 mm was used. For the head and neck and the prostate case the percentage of setup error scenarios passing the acceptance criteria were 29% and 63% (VMAT) and 23% and 58% (DTRT) for translations and 60% and 100% (VMAT) and 62% and 100% (DTRT) for rotations, respectively. Depending on the direction of the translation setup error, the VMAT plan is slightly more robust. However, for the lung case the VMAT treatment plan is more robust leading to passing rates of 48% (VMAT) and 20% (DTRT) for translations and 83% (VMAT) and 55% (DTRT) for rotations, respectively. In all cases, the robustness determined by the rMaps is dominated by the acceptance criteria for the OARs and not by the ones for the CTV. The results of this study indicate that DTRT treatment plans using additional dynamic table and collimator rotations leading to complex trajectories result in similar systematic setup error robustness as VMAT treatment plans for the head and neck and the prostate case. For the lung case investigated, the VMAT treatment plans show an improved robustness.

Evaluation of OAR dose sparing and plan robustness of beam-specific PTV in lung cancer IMRT treatment

PurposeMargins are employed in radiotherapy treatment planning to mitigate the dosimetric effects of geometric uncertainties for the clinical target volume (CTV). Here, we proposed a margin concept that takes into consideration the beam direction, thereby generating a beam-specific planning target volume (BSPTV) on a beam entrance view. The total merged BSPTV was considered a target for optimization. We investigated the impact of this novel approach for lung intensity-modulated radiotherapy (IMRT) treatment, and compared the treatment plans generated using BSPTV with general PTV.Methods and materialsWe generated the BSPTV by expanding the CTV perpendicularly to the incident beam direction using the 2D version of van Herk’s margin concept. The BSPTV and general PTV margin were analyzed using digital phantom simulation. Fifteen lung cancer patients were used in the planning study. First, all patient targets were performed with the CTV projection area analysis to select the suitable beam angles. Then, BSPTV was generated according to the selected beam angles. IMRT plans were optimized with the general PTV and BSPTV as the target volumes, respectively. The dosimetry metrics were calculated and evaluated between these two plans. The plan robustness of both plans for setup uncertainties was evaluated using worst-case analysis.ResultsBoth general PTV and BSPTV plans satisfied the CTV coverage. In addition, the BSPTV plans improved the sparing of high doses to target-surrounding lung tissues compared to the general PTV plans. Both Dmean of Ring PTV and Ring BSPTV were significantly lower in BSPTV plans (38.89 Gy and 39.43 Gy) compared to the general PTV plans (40.27 Gy and 40.68 Gy). The V20, V5, and mean lung dose of the affected lung were significant lower in BSPTV plans (16.20%, 28.75% and 8.93 Gy) compared to general PTV plans (16.69%, 29.22% and 9.18 Gy). In uncertainty scenarios, about 80% of target coverage was achieved for both general PTV and BSPTV plans.ConclusionsThe results suggested that plan robustness can be guaranteed in both the BSPTV and general PTV plans. However, the BSPTV plan spared normal tissues, such as the lungs, significantly better compared to the general PTV plans.

Beam angle comparison for distal esophageal carcinoma patients treated with intensity-modulated proton therapy.

PurposeTo compare the dosimetric performances of intensity‐modulated proton therapy (IMPT) plans generated with two different beam angle configurations (the Right–Left oblique posterior beams and the Superior–Inferior oblique posterior beams) for the treatment of distal esophageal carcinoma in the presence of uncertainties and interplay effect.Methods and MaterialsTwenty patients’ IMPT plans were retrospectively selected, with 10 patients treated with the R‐L oblique posterior beams (Group R‐L) and the other 10 patients treated with the S‐I oblique posterior beams (Group S‐I). Patients in both groups were matched by their clinical target volumes (CTVs—high and low dose levels) and respiratory motion amplitudes. Dose‐volume‐histogram (DVH) indices were used to assess plan quality. DVH bandwidth was calculated to evaluate plan robustness. Interplay effect was quantified using four‐dimensional (4D) dynamic dose calculation with random respiratory starting phase of each fraction. Normal tissue complication probability (NTCP) for heart, liver, and lung was calculated, respectively, to estimate the clinical outcomes. Wilcoxon signed‐rank test was used for statistical comparison between the two groups.ResultsCompared with plans in Group R‐L, plans in Group S‐I resulted in significantly lower liver Dmean and lung V30Gy[RBE] with slightly higher but clinically acceptable spinal cord Dmax. Similar plan robustness was observed between the two groups. When interplay effect was considered, plans in Group S‐I performed statistically better for heart Dmean and V30Gy[RBE], lung Dmean and V5Gy[RBE], and liver Dmean, with slightly increased but clinically acceptable spinal cord Dmax. NTCP for liver was significantly better in Group S‐I.ConclusionsIMPT plans in Group S‐I have better sparing of liver, heart, and lungs at the slight cost of spinal cord maximum dose protection, and are more interplay‐effect resilient compared to IMPT plans in Group R‐L. Our study supports the routine use of the S‐I oblique posterior beams for the treatments of distal esophageal carcinoma.

Intensity-modulated proton therapy (IMPT) interplay effect evaluation of asymmetric breathing with simultaneous uncertainty considerations in patients with non-small cell lung cancer.

Intensity-modulated proton therapy (IMPT) is sensitive to uncertainties from patient setup and proton beam range, as well as interplay effect. In addition, respiratory motion may vary from cycle to cycle, and also from day to day. These uncertainties can severely degrade the original plan quality and potentially affect patient's outcome. In this work, we developed a new tool to comprehensively consider the impact of all these uncertainties and provide plan robustness evaluation under them. We developed a comprehensive plan robustness evaluation tool that considered both uncertainties from patient setup and proton beam range, as well as respiratory motion simultaneously. To mimic patients' respiratory motion, the time spent in each phase was randomly sampled based on patient-specific breathing pattern parameters as acquired during the four-dimensional (4D)-computed tomography (CT) simulation. Spots were then assigned to one specific phase according to the temporal relationship between spot delivery sequence and patients' respiratory motion. Dose in each phase was calculated by summing contributions from all the spots delivered in that phase. The final 4D dynamic dose was obtained by deforming all doses in each phase to the maximum exhalation phase. Three hundred (300) scenarios (10 different breathing patterns with 30 different setup and range uncertainty scenario combinations) were calculated for each plan. The dose-volume histograms (DVHs) band method was used to assess plan robustness. Benchmarking the tool as an application's example, we compared plan robustness under both three-dimensional (3D) and 4D robustly optimized IMPT plans for 10 nonrandomly selected patients with non-small cell lung cancer. The developed comprehensive plan robustness tool had been successfully applied to compare the plan robustness between 3D and 4D robustly optimized IMPT plans for 10 lung cancer patients. In the presence of interplay effect with uncertainties considered simultaneously, 4D robustly optimized plans provided significantly better CTV coverage (D95% , P=0.002), CTV homogeneity (D5% -D95% , P=0.002) with less target hot spots (D5% , P=0.002), and target coverage robustness (CTV D95% bandwidth, P=0.004) compared to 3D robustly optimized plans. Superior dose sparing of normal lung (lung Dmean , P=0.020) favoring 4D plans and comparable normal tissue sparing including esophagus, heart, and spinal cord for both 3D and 4D plans were observed. The calculation time for all patients included in this study was 11.4±2.6min. A comprehensive plan robustness evaluation tool was successfully developed and benchmarked for plan robustness evaluation in the presence of interplay effect, setup and range uncertainties. The very high efficiency of this tool marks its clinical adaptation, highly practical and versatile nature, including possible real-time intra-fractional interplay effect evaluation as a potential application for future use.

Feasibility study: spot-scanning proton arc therapy (SPArc) for left-sided whole breast radiotherapy.

BackgroundThis study investigated the feasibility and potential clinical benefit of utilizing a new proton treatment technique: Spot-scanning proton arc (SPArc) therapy for left-sided whole breast radiotherapy (WBRT) to further reduce radiation dose to healthy tissue and mitigate the probability of normal tissue complications compared to conventional intensity modulated proton therapy (IMPT).MethodsEight patients diagnosed with left-sided breast cancer and treated with breast-preserving surgery followed by whole breast irradiation without regional nodal irradiation were included in this retrospective planning. Two proton treatment plans were generated for each patient: vertical intensity-modulated proton therapy used for clinical treatment (vIMPT, gantry angle 10°–30°) and SPArc for comparison purpose. Both SPArc and vIMPT plans were optimized using the robust optimization of ± 3.5% range and 5 mm setup uncertainties. Root-mean-square deviation dose (RMSD) volume histograms were used for plan robustness evaluation. All dosimetric results were evaluated based on dose-volume histograms (DVH), and the interplay effect was evaluated based on the accumulation of single-fraction 4D dynamic dose on CT50. The treatment beam delivery time was simulated based on a gantry rotation with energy-layer-switching-time (ELST) from 0.2 to 5 s.ResultsThe average D1 to the heart and LAD were reduced to 53.63 cGy and 82.25 cGy compared with vIMPT 110.38 cGy (p = 0.001) and 170.38 cGy (p = 0.001), respectively. The average V5Gy and V20Gy of ipsilateral lung was reduced to 16.77% and 3.07% compared to vIMPT 25.56% (p = 0.001) and 4.68% (p = 0.003). Skin3mm mean and maximum dose were reduced to 3999.38 cGy and 4395.63 cGy compared to vIMPT 4104.25 cGy (p = 0.039) and 4411.63 cGy (p = 0.043), respectively. A significant relative risk reduction (RNTCP = NTCPSPArc/NTCPvIMPT) for organs at risk (OARs) was obtained with SPArc ranging from 0.61 to 0.86 depending on the clinical endpoint. The RMSD volume histogram (RVH) analysis shows SPArc provided better plan robustness in OARs sparing, including the heart, LAD, ipsilateral lung, and skin. The average estimated treatment beam delivery times were comparable to vIMPT plans when the ELST is about 0.5 s.ConclusionSPArc technique can further reduce dose delivered to OARs and the probability of normal tissue complications in patients treated for left-sided WBRT.

A study to investigate the influence of cardiac motion on the robustness of pencil beam scanning proton plans in oesophageal cancer.

While proton therapy offers an excellent dose conformity and sparing of organs at risk, this can be compromised by uncertainties, e.g. organ motion. This study aimed to investigate the influence of cardiac motion on the contoured oesophagus using electrocardiogram-triggered imaging and to assess the impact of this motion on the robustness of proton therapy plans in oesophageal cancer patients. Limited cardiac-induced motion of the oesophagus was observed with a negligible impact on the robustness of proton therapy plans. Therefore, our data suggest that cardiac motion may be safely ignored in the robust optimisation strategy for proton planning in oesophageal cancer.

Clinical implementation of pencil beam scanning proton therapy for liver cancer with forced deep expiration breath hold

PurposeTo present our technique for liver cancer treatments with proton therapy in pencil beam scanning mode and to evaluate the impact of uncertainties on plan quality. Materials and MethodsSeventeen patients affected by liver cancer were included in this study. Patients were imaged and treated in forced breath-hold using the Active Breathing Coordinator system and monitored with an optical tracking system. Three simulation CTs were acquired to estimate the anatomical variability between breath-holds and generate an internal target volume (ITV). The treatment plans were optimized with a Single Field Optimization technique aimed at minimizing the use of range shifter. Plan robustness was tested simulating systematic range and setup uncertainties, as well as the interplay effect between breath-holds. The appropriateness of margin was further verified based on the actual positioning data acquired during treatment. ResultsThe dose distributions of the nominal plans achieved a satisfactory target coverage in 11 out of 17 patients, while in the remaining 6 D95 to the PTV was affected by the constraint on mean liver dose. The constraints for all other organs at risk were always within tolerances.The interplay effect had a limited impact on the dose distributions: the worst case scenario showed a D95 reduction in the ITV < 3.9 GyRBE and no OAR with D1 > 105% of the prescription dose.The robustness analysis showed that for 13 out of 17 patients the ITV coverage in terms of D95 was better than D95 of the PTV in the nominal plan. For the remaining 4 patients, the maximum difference between ITV D95 and PTV D95 was ≤0.7% even for the largest simulated setup error and it was deemed clinically acceptable. Hot spots in the OARs were always lower than 105% of the prescription dose.Positioning images confirmed that the breath hold technique and the PTV margin were adequate to compensate for inter- and intra-breath-hold variations in liver position. ConclusionWe designed and clinically applied a technique for the treatment of liver cancer with proton pencil beam scanning in forced deep expiration breath-hold. The initial data on plan robustness and patient positioning suggest that the choices in terms of planning technique and treatment margins are able to reach the desired balance between target coverage and organ at risk sparing.

Automatic feathering algorithm for VMAT craniospinal irradiation: A comprehensive comparison with other VMAT planning strategies

In craniospinal irradiation, field matching is very sensitive to intrafraction positional uncertainties in cranio-caudal direction, which could lead to severe overdoses/underdoses inside the planning target volume. During the last decade, significant efforts were made to develop volumetric-modulated arc therapy strategies, which were less sensitive to setup uncertainties. In this study, a treatment planning system-integrated method, named automatic feathering (AF) algorithm, was compared against other volumetric-modulated arc therapy strategies. Three patients were retrospectively included. Five different planning techniques were compared, including overlap (O), staggered overlap (SO), gradient optimization (GO), overlap with AF algorithm turned on (O-AF), and staggered overlap with AF algorithm turned on (SO-AF). Three overlapping lengths were considered (5 cm, 7.5 cm, and 10 cm). The middle isocenter was shifted of ±1 mm, ±3 mm, and ±5 mm to simulate setup uncertainties. Plan robustness against simulated uncertainties was evaluated by calculating near maximum and near minimum dose differences between shifted and nonshifted plans (ΔD2%, ΔD98%). Dose differences among combinations of techniques and junction lengths were tested using Wilcoxon signed-rank test. Higher ΔD2% and ΔD98% were obtained using the overlap technique (ΔD2% = 15.4%, ΔD98% = 15.0%). O-AF and SO-AF provided comparable plan robustness to GO technique. Their performance improved significantly for grater overlapping length. For 10-cm overlap and 5-mm shift, GO, O-AF, and SO-AF yielded to the better plan robustness (5.7% < ΔD2% < 6.0%, 6.1% < ΔD98% < 7.6%). SO provided an intermediate plan robustness (9.8% < ΔD2% < 10.8%, 8.9% < ΔD98% < 10.3%). The addition of AF to the overlap technique significantly improves plan robustness especially if larger overlapping lengths are used. Using the AF algorithm, plans become as robust as plans optimized with more sophisticated and time-consuming approaches (like GO).

Pro-con of proton: Dosimetric advantages of intensity-modulation over passive scatter for thoracic malignancies.

The use of passively scattered proton therapy (PSPT) or intensity modulated proton therapy (IMPT) opens the potential for dose escalation or critical structure sparing in thoracic malignancies. While the latter offers greater dose conformality, dose distributions are subjected to greater uncertainties, especially due to interplay effects. Exploration in this area is warranted to determine if there is any dosimetric advantages in using IMPT for thoracic malignancies. This review aims to both compare organs-at-risk sparing and plan robustness between PSPT and IMPT and examine the mitigation strategies for the reduction of interplay effects currently available. Early evidence suggests that IMPT is dosimetrically superior to PSPT in thoracic malignancies. Randomised control trials are required before any clinical benefit of IMPT can be confirmed.

The flexible and real-time commute trip sharing problems

The Commute Trip Sharing Problem (CTSP) was introduced to remove parking pressure on cities, as well as corporate and university campuses. Its goal is to reduce the number of vehicles being used for daily commuting activities. Given a set of inbound and outbound requests, which consists of origin and destination pairs and their departure and return times, the CTSP assigns riders and a driver, as well as inbound and outbound routes, to each vehicle in order to satisfy time-window, capacity, and ride-duration constraints. The CTSP guarantees a ride back for each rider, which is a critical aspect of such a ride-sharing system. This paper generalizes the CTSP to account for uncertainties about the return trip. Each rider is assumed to have a return time specified by a distribution (learned from historical data) and, each day, a percentage of riders will want to preprone or postpone their return trip to accommodate some schedule changes. The paper proposes two generalizations of the CTSP: the Flexible CTSP (FCTSP) and the Real-Time CTSP (RT-CTSP). In the FCTSP, riders must confirm their final return times by a fixed deadline. In the RT-CTSP, riders confirm their new return times in real time with some prior notice. The paper proposes a two-step approach to address the FCTSP and the RT-CTSP. The first step uses a scenario-based stochastic program to choose the drivers and the morning routes in order to maximize the robustness of the driver assignment. The second step reoptimizes the plan at the fixed deadline or in real time once the return times are confirmed. Experiments on a real-world dataset of commute trips demonstrate the effectiveness of the algorithm in generating robust plans and reveal a trade-off between vehicle reduction and plan robustness as the robust plans tend to be conservative. A method is then proposed to evaluate this trade-off using the per-unit price ratio of vehicle increase to uncovered riders.

FRoG: An independent dose and LETd prediction tool for proton therapy at ProBeam® facilities.

Particle therapy is becoming increasingly available world-wide for precise tumor targeting, its favorable depth dose deposition, and increased biological damage to tumor tissue compared to conventional photon therapy. As demand increases for improved robustness and conformality, next-generation secondary dose calculation engines are needed to verify treatment plans independently and provide estimates for clinical decision-making factors, such as dose-averaged linear energy transfer (LETd ) and relative biological effectiveness (RBE). FRoG (Fast dose Recalculation on GPU) has been installed and commissioned at the Danish Centre for Particle Therapy (DCPT). FRoG was developed for synchrotron-based facilities and has previously demonstrated good agreement with gold-standard Monte Carlo simulations and measurements. In this work, additions and modifications to FRoG's pencil beam algorithm to support the ion beam delivery with cyclotron-based technology as used at the DCPT, range shifter (RS) implementation, and robustness analysis methods are presented. FRoG dose predictions are compared to measurements and predictions of the clinical treatment planning system (TPS) Eclipse (Varian Medical Systems, Palo Alto, United States of America, CA, v.13.7.16) in both homogenous and heterogeneous scenarios using a solid-water/water and a half-head anthropomorphic phantom, respectively. Additional capabilities of FRoG are explored by performing a plan robustness analysis, analyzing dose and LETd for ten patients. Mid-target measurements in spread-out Bragg Peaks (SOBP) were on average within -0.19%±0.30% and ≤0.5% of FRoG predictions for irradiations withoutand with RS, respectively. Average 3%/2mm 3D γ-analysis passing rates were 99.1% for ~200 patient plan QA comparisons. Measurement with an anthropomorphic head-phantom yielded a γ-passing rate >98%. Overall, maximum target differences in D02% of <2% between the TPS and FRoG were observed for patient plans. The robustness analysis study accounting for range, delivery, and positioning uncertainties revealed small differences in target dose and a maximum LETd VH02% (LETd received by 2% of the volume having dose larger than 1% of maximum dose) values below 10.1keV/µm to the brain stem. We demonstrate that auxiliary dose calculation systems like FRoG can yield excellent agreement to measurements comparable to clinical beam models. Through this work, application of FRoG as a secondary engine at third party cyclotron-based particle treatment facilities is now established for dose verification as well as providing further insight on LETd and variable RBE distributions for protons, currently absent from the standard clinical TPS.

Retrospective Planning Study of Patients with Superior Sulcus Tumours Comparing Pencil Beam Scanning Protons to Volumetric-Modulated Arc Therapy

AimsTwenty per cent of patients with non-small cell lung cancer present with stage III locally advanced disease. Precision radiotherapy with pencil beam scanning (PBS) protons may improve outcomes. However, stage III is a heterogeneous group and accounting for complex tumour motion is challenging. As yet, it remains unclear as to whom will benefit. In our retrospective planning study, we explored if patients with superior sulcus tumours (SSTs) are a select cohort who might benefit from this treatment. Materials and methodsPatients with SSTs treated with radical radiotherapy using four-dimensional planning computed tomography between 2010 and 2015 were identified. Tumour motion was assessed and excluded if greater than 5 mm. Photon volumetric-modulated arc therapy (VMAT) and PBS proton single-field optimisation plans, with and without inhomogeneity corrections, were generated retrospectively. Robustness analysis was assessed for VMAT and PBS plans involving: (i) 5 mm geometric uncertainty, with an additional 3.5% range uncertainty for proton plans; (ii) verification plans at maximal inhalation and exhalation. Comparative dosimetric and robustness analyses were carried out. ResultsTen patients were suitable. The mean clinical target volume D95 was 98.1% ± 0.4 (97.5–98.8) and 98.4% ± 0.2 (98.1–98.9) for PBS and VMAT plans, respectively. All normal tissue tolerances were achieved. The same four PBS and VMAT plans failed robustness assessment. Inhomogeneity corrections minimally impacted proton plan robustness and made it worse in one case. The most important factor affecting target coverage and robustness was the clinical target volume entering the spinal canal. Proton plans significantly reduced the mean lung dose (by 21.9%), lung V5, V10, V20 (by 47.9%, 36.4%, 12.1%, respectively), mean heart dose (by 21.4%) and thoracic vertebra dose (by 29.2%) (P < 0.05). ConclusionsIn this planning study, robust PBS plans were achievable in carefully selected patients. Considerable dose reductions to the lung, heart and thoracic vertebra were possible without compromising target coverage. Sparing these lymphopenia-related organs may be particularly important in this era of immunotherapy.

Planning constant‐stress accelerated life tests with multiple stresses based on D‐optimal design

AbstractWith the rapid technological advances, products are becoming more reliable. Then, multistress accelerated life testing (MALT) has been adopted in engineering to obtain failure information in a limited time. In order to make the testing procedure more efficient, it is necessary to better design the test plan. However, to date, relevant research on planning of MALT is limited. Multiple stresses will lead to plenty of stress‐level combinations that require too much cost and time to implement. Besides, there may be interactions among multiple stresses, which need more experiments for parameter estimation. To solve these problems, we propose a novel planning method for constant‐stress MALT under lognormal distribution using D‐optimal design, which can reduce required test points efficiently and deal with second‐order effects in models. In ALT, the log‐linear model is often used to describe the life‐stress relationship. Hence, D‐optimal design is adopted in this paper to select test points from the whole test region. Then, optimal unit allocation plans are formulated under V/D‐optimality criterion, respectively, where type I and type II censoring are both discussed. A real case of light‐emitting device (LED) is presented to compare the proposed approach with other two existing methods. The results show that the proposed method performs better than other two existing methods both in prediction accuracy and estimation precision. Moreover, a sensitivity analysis reveals the robustness of the optimal plans determined by the proposed method.

Flexible two‐stage robust model for moving the transmission and reactive power infrastructures expansion planning towards power system integration of renewables

This study proposes a two-stage robust optimisation framework to model uncertainties in construction cost, load and renewable power generations (RPGs) in the coordinated transmission and reactive power infrastructure expansion planning problem. The proposed planning model, based on the network AC model and its voltage constraints, is represented as a mixed-integer linear programing problem with the aim of minimising: (i) the construction costs of lines, (ii) installation cost of reactive power resources and (iii) active load shedding as RPGs curtailment. Using the suggested robust optimisation approach, the best expansion plan is obtained even in the worst case of the defined uncertainty sets and the plan robustness level is controlled via parameter definition for the individual uncertainty sets. The Benders decomposition algorithm is used to solve the two-stage problem where the plan security level is also taken into consideration using the N−1 contingency criterion. The standard Garver 6-bus system and the IEEE 24-bus reliability test system are used to validate the proposed robust framework.

Impact of Gaussian uncertainty assumptions on probabilistic optimization in particle therapy

Range and setup uncertainties in charged particle therapy may induce a discrepancy between the planned and the delivered dose. Countermeasures based on probabilistic (stochastic) optimization usually assume a Gaussian probability density to model the underlying range and setup error. While this standard assumption is generally taken for granted, this study explicitly investigates the dosimetric consequences if the actual range and setup errors obey a different probability density function (PDF) over the course of treatment to the one used during the probabilistic treatment plan optimization. Discrete random sampling was performed for conventionally and probabilistically optimized proton and carbon ion treatment plans utilizing various PDFs that modeled the setup and range error. This method allowed us to assess the treatment plan robustness against different PDFs of conventional and probabilistic plans, which both explicitly assume Gaussian uncertainties. The induced uncertainty in dose was quantified by estimating the expectation value and standard deviation of the RBE-weighted dose for each PDF on the basis of 2500-5000 random dose samples. Probabilistic dose metrics and standard deviation volume histograms were computed to quantify treatment plan robustness of both optimization approaches. It was shown that the classical planning target volume margin extension concept did not compensate the influence of range and setup errors and consequently resulted in a non-negligible average standard deviation in dose of 7.3% throughout the clinical target volume (CTV). In contrast, probabilistic optimization on normally distributed errors yielded treatment plans that not only entailed a lower standard deviation against normally distributed errors accounted for during optimization, but also lower standard deviations for other symmetric PDFs. It was shown that the impact of an incorrect probability distribution assumption is of lower importance after probabilistic optimization as the average uncertainty in the CTV drops to 3.9%. Probabilistic optimization is an effective tool to create robust particle treatment plans. Normally distributed range and setup error assumptions for probabilistic optimization are a reasonable first approximation and yield treatment plans that are also robust against other PDFs.

A class of distribution-free one-sided Cucconi schemes for joint surveillance of location and scale parameters and their application in monitoring cab services

Simultaneous distribution-free surveillance (SDFS) of location and scale parameters of a process based on a single plotting statistic has emerged as a prevalent research issue in the past seven years. Most of the SDFS schemes are, however, appropriate for monitoring two-sided shifts in location-scale. In monitoring service quality or time to an event, among others, decreasing process median or variability are usually considered as a process improvement and not something to worry about in practice. We only desire to have a quick signal when there is an increase in the median or scatter as they refer to a deterioration in the process quality. On the contrary, in network traffic surveillance, a decrease in process variability could be a matter of concern as it may result in network congestion where a decreasing shift may be of more interest. In such instances, specific one-sided shifts are of more importance, and a one-sided SDFS scheme would be more appropriate than traditional two-sided schemes. In this article, we propose several Shewhart-Cucconi type SDFS schemes for jointly monitoring one-sided deviations in location or scale or both parameters of a process. We discuss the implementation steps and establish the in-control robustness of the new surveillance plans. We compare proposed schemes with existing one-sided Shewhart-Lepage type schemes using run-length properties. Numerical results substantiate that the proposed schemes perform quite well in many cases. We explain the application of newly designed schemes in surveillance of waiting time during a taxi ride in some cities.

Weekly robustness evaluation of intensity-modulated proton therapy for oesophageal cancer

Background and purposeIntensity-modulated proton therapy (IMPT) is expected to result in clinical benefits by lowering radiation dose to organs-at-risk (OARs). However, there are concerns about plan robustness due to motion. To address this uncertainty we evaluated the robustness of IMPT compared to the widely clinically used volumetric modulated arc therapy (VMAT) on weekly repeated computed tomographies (CT). Materials and methods19 patients with oesophageal cancer were evaluated. IMPT and VMAT plans were created on a planning 4-Dimensional CT (p4DCT) and evaluated on weekly repeated 4DCTs (r4DCT). In case of inadequate target coverage or unacceptable high dose to normal tissue, re-planning was performed. Dose distributions of the r4DCTs were warped to p4DCT, resulting in an estimated actual given dose (EAGD). ResultsCompared to VMAT, IMPT resulted in significantly lowered dose to heart, lungs, spleen, liver and kidneys. For IMPT, target coverage was adequate (after max 1 replanning) in 17/19 cases. In two cases target coverage remained insufficient. However, in one of these patients the summed dose was insufficient (due to tumor shrinkage) while weekly coverage was adequate. For the other patient the target coverage was also insufficient by VMAT, due to large anatomical changes during treatment. For VMAT, adequate target coverage was achieved in 18/19 cases without re-planning. However, for reasons of high OAR dose re-planning was required in two cases. ConclusionIMPT reduces the dose to OARs significantly, while achieving adequate target coverage in the majority of patients. Re-planning was necessary for both IMPT and VMAT due to anatomical changes.

Investigation of inter-fraction target motion variations in the context of pencil beam scanned proton therapy in non-small cell lung cancer patients.

PurposeFor locally advanced‐stage non‐small cell lung cancer (NSCLC), inter‐fraction target motion variations during the whole time span of a fractionated treatment course are assessed in a large and representative patient cohort. The primary objective is to develop a suitable motion monitoring strategy for pencil beam scanning proton therapy (PBS‐PT) treatments of NSCLC patients during free breathing.MethodsWeekly 4D computed tomography (4DCT; 41 patients) and daily 4D cone beam computed tomography (4DCBCT; 10 of 41 patients) scans were analyzed for a fully fractionated treatment course. Gross tumor volumes (GTVs) were contoured and the 3D displacement vectors of the centroid positions were compared for all scans. Furthermore, motion amplitude variations in different lung segments were statistically analyzed. The dosimetric impact of target motion variations and target motion assessment was investigated in exemplary patient cases.ResultsThe median observed centroid motion was 3.4 mm (range: 0.2–12.4 mm) with an average variation of 2.2 mm (range: 0.1–8.8 mm). Ten of 32 patients (31.3%) with an initial motion <5 mm increased beyond a 5‐mm motion amplitude during the treatment course. Motion observed in the 4DCBCT scans deviated on average 1.5 mm (range: 0.0–6.0 mm) from the motion observed in the 4DCTs. Larger motion variations for one example patient compromised treatment plan robustness while no dosimetric influence was seen due to motion assessment biases in another example case.ConclusionsTarget motion variations were investigated during the course of radiotherapy for NSCLC patients. Patients with initial GTV motion amplitudes of < 2 mm can be assumed to be stable in motion during the treatment course. For treatments of NSCLC patients who exhibit motion amplitudes of > 2 mm, 4DCBCT should be considered for motion monitoring due to substantial motion variations observed.

Alternative method to generate three-dimensional (3D) ultrasonographic images in veterinary medicine

Simultaneous distribution-free surveillance (SDFS) of location and scale parameters of a process based on a single plotting statistic has emerged as a prevalent research issue in the past seven years. Most of the SDFS schemes are, however, appropriate for monitoring two-sided shifts in location-scale. In monitoring service quality or time to an event, among others, decreasing process median or variability are usually considered as a process improvement and not something to worry about in practice. We only desire to have a quick signal when there is an increase in the median or scatter as they refer to a deterioration in the process quality. On the contrary, in network traffic surveillance, a decrease in process variability could be a matter of concern as it may result in network congestion where a decreasing shift may be of more interest. In such instances, specific one-sided shifts are of more importance, and a one-sided SDFS scheme would be more appropriate than traditional two-sided schemes. In this article, we propose several Shewhart-Cucconi type SDFS schemes for jointly monitoring one-sided deviations in location or scale or both parameters of a process. We discuss the implementation steps and establish the in-control robustness of the new surveillance plans. We compare proposed schemes with existing one-sided Shewhart-Lepage type schemes using run-length properties. Numerical results substantiate that the proposed schemes perform quite well in many cases. We explain the application of newly designed schemes in surveillance of waiting time during a taxi ride in some cities.

3D printing for dosimetric optimization and quality assurance in small animal irradiations using megavoltage X-rays

PurposeNew therapeutic options in radiotherapy (RT) are often explored in preclinical in-vivo studies using small animals. We report here on the feasibility of modern megavoltage (MV) linear accelerator (LINAC)-based RT for small animals using easy-to-use consumer 3D printing technology for dosimetric optimization and quality assurance (QA). MethodsIn this study we aimed to deliver 5×2Gy to the half-brain of a rat using a 4MV direct hemi-field X-ray beam. To avoid the beam's build-up in the target and optimize dosimetry, a 1cm thick, customized, 3D-printed bolus was used. A 1:1 scale copy of the rat was 3D printed based on the CT dataset as an end-to-end QA tool. The plan robustness to HU changes was verified. Thermoluminescent dosimeters (TLDs), for both MV irradiations and for kV imaging doses, and a gafchromic film were placed within the phantom for dose delivery verifications. The phantom was designed using a standard treatment planning software, and was irradiated at the LINAC with the target aligned using kV on-board imaging. ResultsThe plan was robust (dose difference<1% for HU modification from 0 to 250). Film dosimetry showed a good concordance between planned and measured dose, with the steep dose gradient at the edge of the hemi-field properly aligned to spare the contralateral half-brain. In the treated region, the mean TLDs percentage dose differences (±2 SD) were 1.3% (±3.8%) and 0.9% (±1.7%) beneath the bolus. The mean (±2 SD) out-of-field dose measurements was 0.05Gy (±0.02Gy) for an expected dose of 0.04Gy. Imaging doses (2mGy) still spared the contralateral-brain. ConclusionsUse of consumer 3D-printers enables dosimetry optimization and QA assessment for small animals MV RT in preclinical studies using standard LINACS.