High-Performance Computing Framework for Predictive Simulation of Healthcare Delivery Innovations
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Cyber-attack accommodation in a cyber-physical system is to ensure system operation, integrity and availability while maintaining a reasonable operational performance under attack. In this paper, we present a novel cyber-attack accommodation algorithm by estimating the true operational states of the system with new boundary & performance constrained resilient estimators while the system is continuously operating and is under attack. Our approach is based on combining data driven machine learning and physics based domain knowledge with traditional resilient estimation. The results were evaluated using a high fidelity model-based simulation environment.
Summary Although β ‐carotene is an essential nutrient for the maintenance of human health, it is insoluble in water and extremely unstable, which limits its application. Here, we prepared a stable emulsification system of whey protein (WP)‐ Lycium barbarum L. leaves flavonoids (LBLF) complex as a carrier to load β ‐carotene. The rheological properties, antioxidant properties and stability of emulsion were studied, and the effects of emulsion on the biological retention of β ‐carotene and stability were investigated. The results indicate that, compared with WP, the WP‐LBLF emulsion has smaller particles, a smaller creaming index, and significantly enhanced chemical stability under light and heat treatment conditions. The thiobarbituric acid value of the WP‐LBLF emulsion (0.031 ± 0.016 mg/kg) was dramatically lower than that of the WP emulsion (0.045 ± 0.012 mg/kg) ( P < 0.05). In the simulated gastrointestinal digestion experiment, the retention rate of β ‐carotene was higher in WP‐LBLF emulsion. Additionally, the WP‐LBLF emulsion loaded with β ‐carotene has good salt tolerance and heat stability. In our study, WP‐LBLF was used as an emulsifier to load β ‐carotene emulsion to improve its bioavailability and stability, which provides a new preparation and application method for a complex functional emulsifier.
The Coriolis force is a significant source of systematic phase errors and dephasing in atom interferometry and is often compensated by counter-rotating the interferometry laser beam against Earth's rotation. We present a novel method for performing Coriolis force compensation for long-baseline atom interferometry, which mitigates atom-beam misalignment due to beam rotation, an effect which is magnified by the long lever arm of the baseline length. The method involves adjustment of the angle of the interferometer beam prior to a magnifying telescope, enabling the beam to pivot around a tunable position along the interferometer baseline. By tuning the initial atom kinematics and adjusting the angle with which the interferometer beam pivots about this point, we can ensure that the atoms align with the center of the beam during the atom optics laser pulses. This approach will be used in the MAGIS-100 atom interferometer and could also be applied to other long-baseline atom interferometers. An additional challenge associated with long baseline interferometry is that since long-baseline atom interferometers are often located outside of typical laboratory environments, facilities constraints may require lasers to be housed in a climate-controlled room a significant distance away from the main experiment. Nonlinear effects in optical fibers restrict the use of fiber-based transport of the high-power interferometry beam from the laser room to the experiment. We present the design of and prototype data from a laser transport system for MAGIS-100 that maintains robustness against alignment drifts despite the absence of a long fiber.
Purpose: To develop an MR accelerator-enabled online planning-todelivery technique for stereotactic palliative radiotherapy treatment of spinal metastases. The technical challenges include; automated stereotactic treatment planning, online MR-based dose calculation and MR guidance during treatment. Methods: Using the CT data of 20 patients previously treated at our institution, a class solution for automated treatment planning for spinal bone metastases was created. For accurate dose simulation right before treatment, we fused geometrically correct online MR data with pretreatment CT data of the target volume (TV). For target tracking during treatment, a dynamic T2-weighted TSE MR sequence was developed. An in house developed GPU based IMRT optimization and dose calculation algorithm was used for fast treatment planning and simulation. An automatically generated treatment plan developed with this treatment planning system was irradiated on a clinical 6 MV linear accelerator and evaluated using a Delta4 dosimeter. Results: The automated treatment planning method yielded clinically viable plans for all patients. The MR-CT fusion based dose calculation accuracy was within 2% as compared to calculations performed with original CT data. The dynamic T2-weighted TSE MR Sequence was able to provide an update of the anatomical location of the TV every 10 seconds. Dose calculation and optimization of the automatically generated treatment plans using only one GPU took on average 8 minutes. The Delta4 measurement of the irradiated plan agreed with the dose calculation with a 3%/3mm gamma pass rate of 86.4%. Conclusions: The development of an MR accelerator-enabled planning-todelivery technique for stereotactic palliative radiotherapy treatment of spinal metastases was presented. Future work will involve developing an intrafraction motion adaptation strategy, MR-only dose calculation, radiotherapy quality-assurance in a magnetic field, and streamlining the entire treatment process on an MR accelerator.
Purpose: Ensuring reproducibility of target localization is critical to accurate stereotactic body radiation treatment (SBRT) for lung and soft tissue metastatic lesions. To characterize interfraction variability in set-up and evaluate PTV margins utilized for SBRT, daily CBCTs were used to calculate delivered target and OAR doses compared to those expected from planning. Methods: CBCT images obtained prior to each fraction of SBRT for a lung and thyroid metastatic lesion were evaluated. The target CTV/ITV and OARs on each of 8 CBCT data sets were contoured. Using MIM fusion software and Pinnacle{sup 3} RTP system, delivered dose distribution was reconstructed on each CBCT, utilizing translational shifts performed prior to treatment. Actual delivered vs. expected doses received by target CTV/ITV and adjacent critical structures were compared to characterize accuracy of pre-treatment translational shifts and PTV margins. Results: The planned CTV/ITV D95% and V100% were 4595cGy and 91.47% for the lung lesion, and 3010cGy and 96.34% for the thyroid lesion. Based on CBCT analysis, actual mean D95% and V100% for lung ITV were 4542±344.4cGy and 91.54±3.45%; actual mean D95% and V100% for thyroid metastasis CTV were 3005±25.98cGy and 95.20±2.522%. For the lung lesion, ipsilateral lung V20, heart V32 (cc) and spinal cord (.03 cc) max were 110.15cc, 3.33cc, and 1680cGy vs. 110.27±14.79cc, 6.74±3.76cc, and 1711±46.56cGy for planned vs. delivered doses, respectively. For the thyroid metastatic lesion, esophagus V18, trachea (.03 cc) max, and spinal cord (.03 cc) max were 0.35cc, 2555cGy, and 850cGy vs. 0.16±0.13cc, 2147±367cGy, and 838±45cGy for planned vs. delivered treatments, respectively. Conclusion: Minimal variability in SBRT target lesion dose delivered based on pre-treatment CBCT-based translational shifts suggests tighter PTV margins may be considered to further decrease dose to surrounding critical structures. Guidelines for optimal target alignment during CBCT-guidance for lung and soft tissue metastatic lesions treated with SBRT are being established.
Purpose: A rapid cycling proton beam has several distinct characteristics superior to a slow extraction synchrotron: The beam energy and energy spread, beam intensity and spot size can be varied spot by spot. The feasibility of using a spot scanning beam from a rapidc-ycling-medical-synchrotron (RCMS) at 10 Hz repetition frequency is investigated in this study for its application in proton therapy. Methods: The versatility of the beam is illustrated by two examples in water phantoms: (1) a cylindrical PTV irradiated by a single field and (2) a spherical PTV irradiated by two parallel opposed fields. A uniform dose distribution is to be delivered to the volumes. Geant4 Monte Carlo code is used to validate the dose distributions in each example. Results: Transverse algorithms are developed to produce uniform distributions in each transverseplane in the two examples with a cylindrical and a spherical PTV respectively. Longitudinally, different proton energies are used in successive transverse planes toproduce the SOBP required to cover the PTVs. In general, uniformity of dosedistribution within 3% is obtained for the cylinder and 3.5% for the sphere. The transversealgorithms requires only few hundred beam spots for each plane The algorithms may beapplied to larger volumes by increasing the intensity spot by spot for the same deliverytime of the same dose. The treatment time can be shorter than 1 minute for any fieldconfiguration and tumor shape. Conclusion: The unique beam characteristics of a spot scanning beam from a RCMS at 10 Hz repetitionfrequency are used to design transverse and longitudinal algorithms to produce uniformdistribution for any arbitrary shape and size of targets. The proposed spot scanning beam ismore versatile than existing spot scanning beams in proton therapy with better beamcontrol and lower neutron dose. This work is supported in part by grants from the US Department of Energy under contract; DE-FG02-12ER41800 and the National Science Foundation NSF PHY-1205431.
Purpose: The purpose of this work was to verify the accuracy of the dose distribution along the field junction in a half beam irradiation technique for breast cancer patients receiving radiation to the breast or chest wall (CW) and the supraclavicular LN region for both free breathing and deep inspiration breath hold (DIBH) technique. Methods: We performed in vivo measurements for nine breast cancer patients receiving radiation to the breast/CW and to the supraclavicular LN region. Six patients were treated to the left breast/CW using DIBH technique and three patients were treated to the right breast/CW in free breath. We used five microMOSFET dosimeters: three located along the field junction, one located 1 cm above the junction and the fifth microMOSFET located 1 cm below the junction. We performed consecutive measurements over several days for each patient and compared the measurements to the TPS calculation (Eclipse, Varian™). Results: The calculated and measured doses along the junction were 0.97±0.08 Gy and 1.02±0.14 Gy, respectively. Above the junction calculated and measured doses were 0.91±0.08 Gy and 0.98±0.09 Gy respectively, and below the junction calculated and measured doses were 1.70±0.15 Gy and 1.61±0.09 Gy, respectively. All differences were not statistically significant. When comparing calculated and measured doses for DIBH patients only, there was still no statistically significant difference between values for all dosimeter locations. Analysis was done using the Mann-Whitney Rank-Sum Test. Conclusion: We found excellent correlation between calculated doses from the TPS and measured skin doses at the junction of several half beam fields. Even for the DIBH technique, where there is more potential for variance due to depth of breath, there is no over or underdose along the field junction. This correlation validates the TPS, as well an accurate, reproducible patient setup.
Purpose: To evaluate the use of standard apertures and range shifters for the treatment of brain metastasis in proton stereotactic radiosurgery. Methods: Five localized brain metastasis patients previously treated using our intracranial proton stereotactic radiosurgery procedure (i.e. with a custom aperture and bolus), were randomly selected from our patient cohort. The custom aperture and bolus treatment plans were used as the standard of care in this case and comparative treatment plans using the standard aperture and range shifter concept were generated. Gantry/table angle and the number of treatment beams were optimized as part of this study to evaluate the ability of the standard aperture/range shifter system to deliver a comparable treatment to the patient. Conformity index, homogeneity index, isodose volumes and integral dose were all evaluated to determine the degree of conformity of the plans created and for comparison to the custom aperture/bolus treatment modality. Results: The generated treatment plans demonstrated that the standard aperture and range shifter combination could be used to produce comparable conformity index and isodose volumes to the custom aperture/bolus case in four out of the five patients studied. In two of the patients a comparative conformity index was achieved by optimizing the angles of the 3 treatment beams, while in two of the cases 1 or 2 additional beams were required. Additionally, this system exhibited efficiency gains of 60-90% over the custom aperture bolus system in reducing the time necessary for treatment planning, device manufacture and QA. Conclusion: This work demonstrated that largely spherical shape of brain metastasis makes this target well suited to an application of standard apertures, while additionally providing efficiency gains in device manufacture and QA for treatment.
Purpose: To investigate the robustness and safety of craniospinal irradiation (CSI) planned with a proton pencil beam scanning (PBS) technique which overcomes the complexity of the planning associated with feathering match lines. Methods: Six CSI patients were planned with gradient-dose matching using PBS technique. Uniform dose coverage to the entire target volumes was achieved with averaged junction lengths of 6.9±0.3 cm. Robustness of the plans was evaluated by shifting the isocenter of each treatment field by ±3 mm in longitudinal direction and compared with the original non-shifted plan with metrics of conformity number (CN) and homogeneity index (HI). An anthropomorphic phantom study using film measurements was also carried out on a plan with 5 cm junction length. Results: For a given junction length, the dose errors were directly proportional to the setup errors. Setup errors of 3 mm from each field caused on average 3.5% lower CN and 2.1% higher HI. Minimal D95% to PTV and D98% to CTV were reduced by 2.2%±1.5% and 2.8%±1.7% respectively. A drop of maximal 6.8%±5.5% on the minimal dose to the cribriform plate was also observed. When the junction length was 5cm or longer, these 3mm setup errors from each field resulted in up to 12% dose errors. Consistent results were reached between film measurements and planned dose profiles in the junction area. Due to near-zero exit doses beyond the target volume, sparing of anterior organs such as heart, liver, lung and kidney were observed. Conclusions: Longitudinal setup errors directly reduce the dosimetric accuracy of the CSI treatment with matched proton fields. The reported technique creates a slow dose gradient in the junction area, which makes the treatment more robust and safe to longitudinal setup errors compared to conventional feathering methods.
Purpose: To evaluate the dosimetric properties of proton dose distributions delivered at extended source-to-skin distances (SSD). Methods: Radiation was delivered with a gantry-mounted proton double scattering system (Proteus 230, IBA). This system has a maximum field diameter of 24 cm at isocenter and a nominal source-to-axis distance of 230 cm. Dose was measured at nominal SSD as well as at -10, +10, +25, +50, +75, and +100 cm for several range and modulation width combinations. Depth dose distributions were measured with a multi-layer ionization chamber (MLIC), and lateral dose distributions with a 2D ionization chamber array as well as with a diode in a water phantom. Results: The maximum field diameter was found to increase from 24.0 cm at nominal SSD to 29.1 cm and 33.3 cm at +50 cm and +100 cm respectively. Field flatness remained below 3% for all SSD. Tilt in the spread-out Bragg peak depth dose distribution increased with SSD up to 0.4 %-per-g/cm2 at +100 cm. The measured range decreased with 1.1x10-3 g/cm2 per centimeter shift in SSD due to proton energy loss in air, very close to the theoretically calculated value of 1.06x10-3 g/cm3. The output and dose rate decrease with the inverse of the SSD squared as expected. Conclusion: Extending the SSD up to 100 cm increases the maximum field diameter from 24.0 cm to 33.3 cm while the dose uniformity remains acceptable.
Purpose: To evaluate the 2D array for radiation response for a period of 2 years. Methods: 45 VMAT plans already treated and quality assured before 2 years, 1.5 years and 1 year were selected for the study. Quality assurances of the plans were done using 2D array combined with Octavius phantom. Verification plans were recalculated without changing any parameters in Eclipse 10.0 TPS using the AAA algorithm. Response of 2D array to the plans treated before 2 years, 1.5 years and 1 year where evaluated. The results were analyzed using the Gamma analysis method with the standard gamma passing criteria of 3mm distance to agreement (DTA) and 3% dose difference (DD). Results: All the plans evaluated passed the gamma analysis with a percentage greater than 95, except for three cases. Higher gamma passing criteria where observed for all the analyzed plans, when analysis done before 2 years, 1.5 years and 1 year. The standard deviation of ± 1.38, ± 1.40,± 0.97 where observed between the plans when verification plans did before 2 years, 1.5 years and 1 year respectively. Same set of plans shows a standard deviation of ± 0.70, ± 1.36, and ± 1.18, when analysis done recently. A significance difference in response of the array when analysis done recently for the verification plans treated and quality assured before 2 years.This indicates a slightly reduced response of 2D array towards radiation response as the array gets older. Conclusion: It is found that 2D array shows a reduced response against radiation detection over a period of years. An onsite calibration of the instrument is recommended before the measurements. A dose correction factor can be applied if necessary to the results if the radiation response and efficiency of the array is very poor.
Purpose: To develop treatment planning workflow for rapid radiotherapy delivered with very-high energy electron (VHEE) scanning beam. Methods: VHEE radiotherapy treatment planning was performed by linking Monte Carlo (MC) dose calculations with inverse optimization in a research version of RayStation. In order to study a number of treatment parameters, a Matlab graphical user interface (GUI) for calculation of VHEE beamlet dose was developed. Through the GUI, EGSnrc MC simulations were run for a number of beam energies, number of beams, beamlet spot and grid sizes, and machine bore sizes. VHEE plans for a pediatric patient with a 4.3 cm{sup 3} brain target optimized with spot-scanning algorithm in RayStation were compared to the clinically delivered 6 MV VMAT plan. Results and Discussion: VHEE beam energy had the largest effect on the quality of dose distributions. For the same target dose, the mean doses to critical organs decreased by 10–15% when planned with 100 MeV compared to 60 MeV. VHEE plans calculated with 36 beams outperformed plans calculated with 13 and 17 beams. While beamlet spacing and bore size had a small effect on VHEE dose distributions, 0.1-3mm beamlet sizes resulted in identical dose distributions. Critical organ doses were by up to 70% lower in the best VHEE plan compared to the clinical 6 MV VMAT plan. Conclusions: We have developed a GUI for MC beamlet generation for treatment planning of VHEE radiotherapy. We have demonstrated that pediatric VHEE plans resulted in significant critical organ dose sparing compared to the clinical VMAT plan.
Partial breast irradiation (PBI) following breast-conserving surgery is emerging as an effective means to achieve local control and reduce irradiated breast volume. Patients are planned on a static CT image; however, treatment is delivered while the patient is free-breathing. Respiratory motion can degrade plan quality by reducing target coverage and/or dose homogeneity. A variety of methods can be used to determine the required margin for respiratory motion in PBI. We derive geometric and dosimetric respiratory 1D margin. We also verify the adequacy of the typical 5 mm respiratory margin in 3D by evaluating plan quality for increasing respiratory amplitudes (2–20 mm). Ten PBI plans were used for dosimetric evaluation. A database of volunteer respiratory data, with similar characteristics to breast cancer patients, was used for this study. We derived a geometric 95%-margin of 3 mm from the population respiratory data. We derived a dosimetric 95%-margin of 2 mm by convolving 1D dose profiles with respiratory probability density functions. The 5 mm respiratory margin is possibly too large when 1D coverage is assessed and could lead to unnecessary normal tissue irradiation. Assessing margins only for coverage may be insufficient; 3D dosimetric assessment revealed degradation in dose homogeneity is the limiting factor, not target coverage. Hotspots increased even for the smallest respiratory amplitudes, while target coverage only degraded at amplitudes greater than 10 mm. The 5 mm respiratory margin is adequate for coverage, but due to plan quality degradation, respiratory management is recommended for patients with respiratory amplitudes greater than 10 mm.
Purpose: To describe and evaluate a novel system for generalized Real-Time Interactive Planning (RTIP) applied to head and neck (H and N) VMAT. Methods: The clinician interactively manipulates dose distributions using DVHs, isodoses, or rate of dose fall-off, which may be subjected to user-defined constraints. Dose is calculated using a fast Achievable Dose Estimate (ADE) algorithm, which simulates the limits of what can be achieved during treatment. After each manipulation contributing fluence elements are modified and the dose distribution updates in effectively real-time. For H and N VMAT planning, structure sets for 11 patients were imported into RTIP. Each dose distribution was interactively modified to minimize OAR dose while constraining target DVHs. The resulting RTIP DVHs were transferred to the Eclipse™ VMAT optimizer, and conventional VMAT optimization was performed. Results: Dose calculation and update times for the ADE algorithm ranged from 2.4 to 22.6 milliseconds, thus facilitating effectively real-time manipulation of dose distributions. For each of the 11 H and N VMAT cases, the RTIP process took ∼2–10 minutes. All RTIP plans exhibited acceptable PTV coverage, mean dose, and max dose. 10 of 11 RTIP plans achieved substantially improved sparing of one or more OARs without compromising dose to targets or other OARs. Importantly, 10 of the 11 RTIP plans required only one or two post-RTIP optimizations. Conclusions: RTIP is a novel system for manipulating and updating achievable dose distributions in real-time. H and N VMAT plans generated using RTIP demonstrate improved OAR sparing and planning efficiency. Disclosures: One author has a commercial interest in the presented materials.
Purpose: To efficiently deliver respiratory-gated radiation during synchrotron-based pulsed heavy-ion radiotherapy, a novel respiratory guidance method combining a personalized audio-visual biofeedback (BFB) system, breath hold (BH), and synchrotron-based gating was designed to help patients synchronize their respiratory patterns with synchrotron pulses and to overcome typical limitations such as low efficiency, residual motion, and discomfort. Methods: In-house software was developed to acquire body surface marker positions and display BFB, gating signals, and real-time beam profiles on a LED screen. Patients were prompted to perform short BHs or short deep breath holds (SDBH) with the aid of BFB following a personalized standard BH/SDBH (stBH/stSDBH) guiding curve or their own representative BH/SDBH (reBH/reSDBH) guiding curve. A practical simulation was performed for a group of 15 volunteers to evaluate the feasibility and effectiveness of this method. Effective dose rates (EDRs), mean absolute errors between the guiding curves and the measured curves, and mean absolute deviations of the measured curves were obtained within 10%–50% duty cycles (DCs) that were synchronized with the synchrotron’s flat-top phase. Results: All maneuvers for an individual volunteer took approximately half an hour, and no one experienced discomfort during the maneuvers. Using the respiratory guidance methods, the magnitude of residual motion was almost ten times less than during nongated irradiation, and increases in the average effective dose rate by factors of 2.39–4.65, 2.39–4.59, 1.73–3.50, and 1.73–3.55 for the stBH, reBH, stSDBH, and reSDBH guiding maneuvers, respectively, were observed in contrast with conventional free breathing-based gated irradiation, depending on the respiratory-gated duty cycle settings. Conclusions: The proposed respiratory guidance method with personalized BFB was confirmed to be feasible in a group of volunteers. Increased effective dose rate and improved overall treatment precision were observed compared to conventional free breathing-based, respiratory-gated irradiation. Because breathing guidance curves could be established based on the respective average respiratory period and amplitude for each patient, it may be easier for patients to cooperate using this technique.
Purpose: The purpose of the current study was to characterize and evaluate intensity-modulated proton grid therapy (IMgPT) using a clinical proton beam. Methods: A TOPAS MC model of a new developmental mode (pre-clinical) of the Hitachi proton therapy system (PROBEAT) was used for simulation and characterization of proton grid therapy. TOPAS simulations of different energy ranges, depths and spot separation distances were performed. LET spectra for various energies and depths were produced with FLUKA MC code for evaluation potential interplay between planning parameters and their effect on the characterization of areas (valley) between spots. IMgPT planning aspects (spot spacing, skin dose, peak-to-valley ratios, beam selection, etc.) were evaluated for different phantom and patient cases. Raysearch software (v4.51) was used to perform the evaluation. Results: Calculated beam peak-to-valley ratios scenarios showed strong energy and depth dependence with ratios to be larger for higher energies and shallower depths. Peak-to-valley ratios for R90 range and for spot spacing of 1cm varied from 30% (E = 221.3 MeV, depth 30.6 cm) to 80% (E = 70.3 MeV, depth 4 cm). LET spectra calculations showed spectral hardening with depth, which might potential increase, spot separation distance and improve peak-to-valley ratios. IMgPT optimization, using constant spot spacing, showed skin dose reduction between peak regions of dose due to the irradiation of less skin. Single beam for bulky shallower tumors might be a potential candidate for proton grid therapy. Conclusions: Proton grid therapy using a clinical beam is a promising technique that reduces skin dose between peak regions of dose and may be suitable for the treatment of shallow tumors. IMgPT may be considered for use when bystander effects in off peak regions would be appropriate.