Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set methods”

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.

At least 973 records · Page 54

Self-interaction corrected density functional calculations of molecular Rydberg states

A method is presented for calculating the wave function and energy of Rydberg excited states of molecules. A good estimate of the Rydberg state orbital is obtained using ground state density functional theory including Perdew-Zunger self-interaction correction and an optimized effective potential. The total energy of the excited molecule is obtained using the Delta Self-Consistent Field method where an electron is removed from the highest occupied orbital and placed in the Rydberg orbital. Results are presented for the first few Rydberg states of NH{sub 3}, H{sub 2}O, H{sub 2}CO, C{sub 2}H{sub 4}, and N(CH{sub 3}){sub 3}. The mean absolute error in the energy of the 33 molecular Rydberg states presented here is 0.18 eV. The orbitals are represented on a real space grid, avoiding the dependence on diffuse atomic basis sets. As in standard density functional theory calculations, the computational effort scales as NM{sup 2} where N is the number of orbitals and M is the number of grid points included in the calculation. Due to the slow scaling of the computational effort with system size and the high level of parallelism in the real space grid approach, the method presented here makes it possible to estimate Rydberg electron binding energy in large molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-G-201-02: Application of RayStretch in Clinical Cases: A Calculation for Heterogeneity Corrections in LDR Permanent I-125 Prostate Brachytherapy

Purpose: Tissue heterogeneities and calcifications have significant impact on the dosimetry of low energy brachytherapy (BT). RayStretch is an analytical algorithm developed in our institution to incorporate heterogeneity corrections in LDR prostate brachytherapy. The aim of this work is to study its application in clinical cases by comparing its predictions with the results obtained with TG-43 and Monte Carlo (MC) simulations. Methods: A clinical implant (71 I-125 seeds, 15 needles) from a real patient was considered. On this patient, different volumes with calcifications were considered. Its properties were evaluated in three ways by i) the Treatment planning system (TPS) (TG-43), ii) a MC study using the Penelope2009 code, and iii) RayStretch. To analyse the performance of RayStretch, calcifications located in the prostate lobules covering 11% of the total prostate volume and larger calcifications located in the lobules and underneath the urethra for a total occupied volume of 30% were considered. Three mass densities (1.05, 1.20, and 1.35 g/cm3) were explored for the calcifications. Therefore, 6 different scenarios ranging from small low density calcifications to large high density ones have been discussed. Results: DVH and D90 results given by RayStretch agree within 1% with the full MC simulations. Although no effort has been done to improve RayStretch numerical performance, its present implementation is able to evaluate a clinical implant in a few seconds to the same level of accuracy as a detailed MC calculation. Conclusion: RayStretch is a robust method for heterogeneity corrections in prostate BT supported on TG-43 data. Its compatibility with commercial TPSs and its high calculation speed makes it feasible for use in clinical settings for improving treatment quality. It will allow in a second phase of this project, its use during intraoperative ultrasound planning. This study was partly supported by a fellowship grant from the Spanish Ministry of Education, by the Generalitat Valenciana under Project PROMETEOII/2013/010, by the Spanish Government under Project No. FIS2013-42156 and by the European Commission within the SeventhFramework Program through ENTERVISION (grant agreement number 264552).

60 APPLIED LIFE SCIENCES↗

Satellite/radiosonde comparison

Collocated measurements of radiosonde and TIROS Operational Vertical Sounder (TOVS) data are compared to determine whether significant latitudinal, seasonal, or instrumental biases exist. Long-wave structure determined by each system is being examined for correlation. Because U.S. radiosonde data are archived uncorrected, an effort to correct the measurements is being considered. These corrections will be applied to 7 years of measurements. New comparisons will then be made using corrected radiosonde and TOVS data. A regression method is used to retrieve temperatures, and radiosonde zonal means are used to update the regression coefficients. Preliminary results of auto-correlation analysis indicate that detection of long-wave structure in each set of data is slightly different with some levels and locations experiencing different wave periods or amplitudes. Scatter diagrams and linear regression for each layer show bias in the mean temperatures. Winter-time measurements show poorer correlations due to noisy measurements.

Schmidlin, F. J.↗

The NASA Ocean Data System at the Jet Propulsion Laboratory

Global change studies require satellite oceanographic observations archived in an active data center. This paper describes the activities of an experimental active oceanographic data archive center specializing in distribution of high-level oceanographic data products produced by research oceanographers from satellite altimeter, scatterometer, and microwave and infrared radiometer measurements. Methods to encourage data sharing practices were developed. Services initiated to help oceanographers conduct research more effectively included rapid delivery of data, reprocessing data sets, subsetting of large data sets, and distribution of multiple data sets on a single medium.

Halpern, David↗

Interfacial Reaction Studies Using ONIOM

In this report, we focus on the calculations of the energetics and chemical kinetics of heterogeneous reactions for Organometallic vapor phase epitaxy (OMVPE). The work described in this report builds upon our own previous thermochemical and chemical kinetics studies. The first of these articles refers to the prediction of thermochemical properties, and the latter one deals with the prediction of rate constants for gaseous homolytic dissociation reactions. The calculations of this investigation are at the microscopic level. The systems chosen consisted of a gallium nitride (GaN) substrate, and molecular nitrogen (N2) and ammonia (NH3) as adsorbants. The energetics for the adsorption and the adsorbant dissociation processes were estimated, and reaction rate constants for the dissociation reactions of free and adsorbed molecules were predicted. The energetics for substrate decomposition was also computed. The ONIOM method, implemented in the Gaussian98 program, was used to perform the calculations. This approach has been selected since it allows dividing the system into two layers that can be treated at different levels of accuracy. The atoms of the substrate were modeled using molecular mechanics6 with universal force fields, whereas the adsorbed molecules were approximated using quantum mechanics, based on density functional theory methods with B3LYP functionals and 6-311G(d,p) basis sets. Calculations for the substrate were performed in slabs of several unit cells in each direction. The N2 and NH3 adsorbates were attached to a central location at the Ga-lined surface.

Cardelino, Beatriz H.↗

The Challenges of Air Monitoring at a Soil Remediation Project in a Suburban Setting - 20448

Air monitoring during soil excavation at a chemically and radiologically impacted site in a suburban setting presents unique challenges. Site activities were conducted in a manner that minimized the release of airborne particulates of concern, principally radionuclides of concern (ROC). Dust and air monitoring were performed to ensure the public was not being exposed to the soil contaminants. The project team collected representative air and dust data while being mindful to not impose on the lives of residents in the area. Due to the close proximity to the public, there was a high level of oversight from the client and regulatory agencies. The challenges involved with air sampling at this soil remediation site located in a suburban setting included the number of contaminants, the ranges of risks and hazards presented by the contaminants and the method of remediation, the concerns of the regulators and the public, and the perception of risk from the public. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Uncertainty quantification methodology for hyperbolic systems with application to blood flow in arteries

We present a Stochastic Finite Volume - ADER (SFV-ADER) methodology for Uncertainty Quantification (UQ) in the general framework of systems of hyperbolic balance laws with uncertainty in parameters. The resulting method is weakly intrusive, meaning that deterministic solvers need minor modifications to include random contributions, and has no theoretical accuracy barrier. An illustration of the second-order version for the viscous Burgers equation with uncertain viscosity coefficient is first given in detail, under different deterministic initial conditions; the attainment of the theoretically expected convergence rate is demonstrated empirically; results and features of our scheme are discussed. We then extend the SFV-ADER method to a non-linear hyperbolic system with source terms that models one-dimensional blood flow in arteries, assuming uncertainties in a set of parameters of the problem. Results are then compared with measurements in a well-defined 1:1 experimental replica of the network of largest arteries in the human systemic circulation. For a 99%-confidence level, severe influence of parameter fluctuations is seen in pressure profiles, while secondary effects on flow rate profiles become visible, particularly in terminal branches. In conclusion, results suggest that the proposed methodology could successfully be applied to other physical problems involving hyperbolic balance laws.

60 APPLIED LIFE SCIENCES↗

CONSTRAINTS ON LEMAITRE-TOLMAN-BONDI MODELS FROM OBSERVATIONAL HUBBLE PARAMETER DATA

We use the observational Hubble parameter data (OHD), both the latest observational data set (Stern et al.; referred to as SJVKS) and the simulated data sets, to constrain Lemaitre-Tolman-Bondi (LTB) void models. The necessity of the consistency check on OHD itself in the LTB cosmology is stressed. Three voids are chosen as test models and are constrained using the Union2 data set of SN Ia as well as OHD. Despite their different parameterization, the results from our test models show some similarities; e.g., the best-fit voids obtained from OHD are all considerably broader than those from SN Ia. Due to the small size of the SJVKS data set, the constraints are not conclusive. The constraining power of the future OHD observations are therefore investigated through a figure of merit (FoM) analysis based on the Monte Carlo simulated data. We found that, in the case that the future OHD become more consistent with SN Ia, the results from the test models are almost unanimous: (1) as many as 32 OHD data points at the SJVKS-like uncertainty level are needed to give a higher FoM than the Union2 data set of SN Ia, (2) precise observation helps reduce this required number, (3) increasing the survey depth does not always increase the FoM. On the other hand, if the future OHD and the Union2 data set keep favoring different voids, in a similar manner as they do at present, the 1{sigma} confidence regions obtained from the two probes should finally separate. We test this conjecture and find that the minimum observational requirement (the size of the data set, the uncertainty level, and the survey depth) for this inconsistency to emerge depends strongly on the void model.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of Spectroscopic Techniques for Determining the Peroxide Value of 19 Classes of Naturally Aged, Plant-Based Edible Oils

The peroxide value of edible oils is a measure of the degree of oxidation, which directly relates to the freshness of the oil sample. Several studies previously reported in the literature have paired various spectroscopic techniques with multivariate analyses to rapidly determine peroxide values using field portable and process instrumentation; those efforts presented “best-case scenarios” with oils from narrowly defined training and test sets. The purpose of this paper is to evaluate the use of near- and mid-infrared absorption and Raman scattering spectroscopies on oil samples from different oil classes, including seasonal and vendor variations, to determine which measurement technique or combination thereof is best for predicting peroxide values. Following peroxide value assays of each oil class using an established titration-based method, global and global-subset calibration models were constructed from spectroscopic data collected on the 19 oil classes used in this study. Spectra from each optical technique were used to create partial least squares regression calibration models to predict the peroxide value of unknown oil samples. A global peroxide value model based on near-infrared (8 mm optical path length) oil measurements produced the lowest RMSEP (4.9), followed by 24 mm optical path length near infrared (5.1), Raman (6.9) and 50 µm optical path length mid-infrared (7.3). However, it was determined that the Raman RMSEP resulted from chance correlations. Global peroxide value models based on low-level fusion of the NIR (8 and 24 mm optical path length) data and all infrared data produced the same RMSEP of 5.1. Global subset models, based on any of the spectroscopies and olive oil training sets from any class (pure, extra light, extra virgin), all failed to extrapolate to the non-olive oils. However, the near-infrared global subset model built on extra virgin olive oil could extrapolate to test samples from other olive oil classes. This work demonstrates the difficulty of developing a truly global method for determining peroxide value of oils.

36 MATERIALS SCIENCE↗

An integrated model-driven method for in-treatment upper airway motion tracking using cine MRI in head and neck radiation therapy

Purpose: For the first time, MRI-guided radiation therapy systems can acquire cine images to dynamically monitor in-treatment internal organ motion. However, the complex head and neck (H&N) structures and low-contrast/resolution of on-board cine MRI images make automatic motion tracking a very challenging task. In this study, the authors proposed an integrated model-driven method to automatically track the in-treatment motion of the H&N upper airway, a complex and highly deformable region wherein internal motion often occurs in an either voluntary or involuntary manner, from cine MRI images for the analysis of H&N motion patterns. Methods: Considering the complex H&N structures and ensuring automatic and robust upper airway motion tracking, the authors firstly built a set of linked statistical shapes (including face, face-jaw, and face-jaw-palate) using principal component analysis from clinically approved contours delineated on a set of training data. The linked statistical shapes integrate explicit landmarks and implicit shape representation. Then, a hierarchical model-fitting algorithm was developed to align the linked shapes on the first image frame of a to-be-tracked cine sequence and to localize the upper airway region. Finally, a multifeature level set contour propagation scheme was performed to identify the upper airway shape change, frame-by-frame, on the entire image sequence. The multifeature fitting energy, including the information of intensity variations, edge saliency, curve geometry, and temporal shape continuity, was minimized to capture the details of moving airway boundaries. Sagittal cine MR image sequences acquired from three H&N cancer patients were utilized to demonstrate the performance of the proposed motion tracking method. Results: The tracking accuracy was validated by comparing the results to the average of two manual delineations in 50 randomly selected cine image frames from each patient. The resulting average dice similarity coefficient (93.28% ± 1.46%) and margin error (0.49 ± 0.12 mm) showed good agreement between the automatic and manual results. The comparison with three other deformable model-based segmentation methods illustrated the superior shape tracking performance of the proposed method. Large interpatient variations of swallowing frequency, swallowing duration, and upper airway cross-sectional area were observed from the testing cine image sequences. Conclusions: The proposed motion tracking method can provide accurate upper airway motion tracking results, and enable automatic and quantitative identification and analysis of in-treatment H&N upper airway motion. By integrating explicit and implicit linked-shape representations within a hierarchical model-fitting process, the proposed tracking method can process complex H&N structures and low-contrast/resolution cine MRI images. Future research will focus on the improvement of method reliability, patient motion pattern analysis for providing more information on patient-specific prediction of structure displacements, and motion effects on dosimetry for better H&N motion management in radiation therapy.

60 APPLIED LIFE SCIENCES↗

WE-E-18A-03: How Accurately Can the Peak Skin Dose in Fluoroscopy Be Determined Using Indirect Dose Metrics?

Purpose: Skin dosimetry is important for fluoroscopically-guided interventions, as peak skin doses (PSD) that Result in skin reactions can be reached during these procedures. The purpose of this study was to assess the accuracy of different indirect dose estimates and to determine if PSD can be calculated within ±50% for embolization procedures. Methods: PSD were measured directly using radiochromic film for 41 consecutive embolization procedures. Indirect dose metrics from procedures were collected, including reference air kerma (RAK). Four different estimates of PSD were calculated and compared along with RAK to the measured PSD. The indirect estimates included a standard method, use of detailed information from the RDSR, and two simplified calculation methods. Indirect dosimetry was compared with direct measurements, including an analysis of uncertainty associated with film dosimetry. Factors affecting the accuracy of the indirect estimates were examined. Results: PSD calculated with the standard calculation method were within ±50% for all 41 procedures. This was also true for a simplified method using a single source-to-patient distance (SPD) for all calculations. RAK was within ±50% for all but one procedure. Cases for which RAK or calculated PSD exhibited large differences from the measured PSD were analyzed, and two causative factors were identified: ‘extreme’ SPD and large contributions to RAK from rotational angiography or runs acquired at large gantry angles. When calculated uncertainty limits [−12.8%, 10%] were applied to directly measured PSD, most indirect PSD estimates remained within ±50% of the measured PSD. Conclusions: Using indirect dose metrics, PSD can be determined within ±50% for embolization procedures, and usually to within ±35%. RAK can be used without modification to set notification limits and substantial radiation dose levels. These results can be extended to similar procedures, including vascular and interventional oncology. Film dosimetry is likely an unnecessary effort for these types of procedures.

60 APPLIED LIFE SCIENCES↗

Can fractal dimensions objectivize gastropod shell morphometrics? A case study from Lake Lugu (SW China)

1. Morphometrics are fundamental for the analysis of size and shape in fossils, particularly because soft parts or DNA are rarely preserved and hard parts such as shells are commonly the only source of information. Geometric morphometrics, that is, landmark analysis, is well established for the description of shape but it exhibits a couple of shortcomings resulting from subjective choices during landmarking (number and position of landmarks) and from difficulties in resolving shape at the level of micro-sculpture. 2. With the aid of high-resolution 3D scanning technology and analyses of fractal dimensions, we test whether such shortcomings of linear and landmark morphometrics can be overcome. As a model group, we selected a clade of modern viviparid gastropods from Lake Lugu, with shells that show a high degree of sculptural variation. Linear and landmark analyses were applied to the same shells in order to establish the fractal dimensions. 3. The genetic diversity of the gastropod clade was assessed. The genetic results suggest that the gastropod clade represents a single species. The results of all morphometric methods applied are in line with the genetic results, which is that no specific morphotype could be delimited. Apart from this overall agreement, landmark and fractal dimension analyses do not correspond to each other but represent data sets with different information. Generally, the fractal dimension values quantify the roughness of the shell surface, the resolution of the 3D scans determining the level. In our approach, we captured the micro-sculpture but not the first-order sculptural elements, which explains that fractal dimension and landmark data are not in phase. 4. We can show that analyzing fractal dimensions of gastropod shells opens a window to more detailed information that can be considered in evolutionary and ecological contexts. We propose that using low-resolution 3D scans may successfully substitute landmark analyses because it overcomes the subjective landmarking. Analyses of 3D scans with higher resolution than used in this study will provide surface roughness information at the mineralogical level. We suggest that fractal dimension analyses of a combination of differently resolved 3D models will significantly improve the quality of shell morphometrics.

54 ENVIRONMENTAL SCIENCES↗

Positivity Limiters for Filtered Spectral Approximations of Linear Kinetic Transport Equations

We analyze the properties and compare the performance of several positivity limiters for spectral approximations with respect to the angular variable of linear transport equations. It is well-known that spectral methods suffer from the occurrence of (unphysical) negative spatial particle concentrations due to the fact that the underlying polynomial approximations are not always positive at the kinetic level. Positivity limiters address this defect by enforcing positivity of the polynomial approximation on a finite set of preselected points. With a proper PDE solver, they ensure positivity of the particle concentration at each step in a time integration scheme. We review several known positivity limiters proposed in other contexts and also introduce a modification for one of them. We give error estimates for the consistency of the positive approximations produced by these limiters and compare the theoretical estimates to numerical results. We then solve two benchmark problems with these limiters, make qualitative and quantitative observations about the solutions, and then compare the efficiency of the different limiters.

97 MATHEMATICS AND COMPUTING↗

Time-Adjusted Internal Target Volume: A Novel Approach Focusing on Heterogeneity of Tumor Motion Based on 4-Dimensional Computed Tomography Imaging for Radiation Therapy Planning of Lung Cancer

Purpose: To consider nonuniform tumor motion within the internal target volume (ITV) by defining time-adjusted ITV (TTV), a volume designed to include heterogeneity of tumor existence on the basis of 4-dimensional computed tomography (4D-CT). Methods and Materials: We evaluated 30 lung cancer patients. Breath-hold CT (BH-CT) and free-breathing 4D-CT scans were acquired for each patient. The tumors were manually delineated using a lung CT window setting (window, 1600 HU; level, −300 HU). Tumor in BH-CT images was defined as gross tumor volume (GTV), and the sum of tumors in 4D-CT images was defined as ITV-4D. The TTV images were generated from the 4D-CT datasets, and the tumor existence probability within ITV-4D was calculated. We calculated the TTV{sub 80} value, which is the percentage of the volume with a tumor existence probability that exceeded 80% on ITV-4D. Several factors that affected the TTV{sub 80} value, such as the ITV-4D/GTV ratio or tumor centroid deviation, were evaluated. Results: Time-adjusted ITV images were acquired for all patients, and tumor respiratory motion heterogeneity was visualized. The median (range) ITV-4D/GTV ratio and median tumor centroid deviation were 1.6 (1.0-4.1) and 6.3 mm (0.1-30.3 mm), respectively. The median TTV{sub 80} value was 43.3% (2.9-98.7%). Strong correlations were observed between the TTV{sub 80} value and the ITV-4D/GTV ratio (R=−0.71) and tumor centroid deviation (R=−0.72). The TTV images revealed the tumor motion pattern features within ITV. Conclusions: The TTV images reflected nonuniform tumor motion, and they revealed the tumor motion pattern features, suggesting that the TTV concept may facilitate various aspects of radiation therapy planning of lung cancer while incorporating respiratory motion in the future.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Differential cross sections for electron impact excitation of the electronic bands of phenol

We report results from a joint theoretical and experimental investigation into electron scattering from the important organic species phenol (C{sub 6}H{sub 5}OH). Specifically, differential cross sections (DCSs) have been measured and calculated for the electron-impact excitation of the electronic states of C{sub 6}H{sub 5}OH. The measurements were carried out at energies in the range 15–40 eV, and for scattered-electron angles between 10{sup ∘} and 90{sup ∘}. The energy resolution of those experiments was typically ∼80 meV. Corresponding Schwinger multichannel method with pseudo-potentials calculations, with and without Born-closure, were also performed for a sub-set of the excited electronic-states that were accessed in the measurements. Those calculations were conducted at the static exchange plus polarisation (SEP)-level using a minimum orbital basis for single configuration interaction (MOBSCI) approach. Agreement between the measured and calculated DCSs was typically fair, although to obtain quantitative accord, the theory would need to incorporate even more channels into the MOBSCI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Going Forward with the Nancy Grace Roman Space Telescope Transient Survey: Validation of Precision Forward-modeling Photometry for Undersampled Imaging

The Nancy Grace Roman Space Telescope (Roman) is an observatory for both wide-field observations and coronagraphy that is scheduled for launch in the mid-2020s. Part of the planned survey is a deep, cadenced field or fields that enable cosmological measurements with type Ia supernovae (SNe Ia). With a pixel scale of 0$_.^"$11, the Wide Field Instrument will be undersampled, presenting a difficulty for precisely subtracting the galaxy light underneath the SNe. We use simulated data to validate the ability of a forward-model code (such codes are frequently also called “scene-modeling” codes) to perform precision supernova photometry for the Roman SN survey. Our simulation includes over 760,000 image cutouts around SNe Ia or host galaxies (~10% of a full-scale survey). To have a realistic 2D distribution of underlying galaxy light, we use the VELA simulated high-resolution images of galaxies. We run each set of cutouts through our forward-modeling code which automatically measures time-dependent SN fluxes. Given our assumed inputs of a perfect model of the instrument point-spread functions and calibration, we find biases at the millimagnitude level from this method in four red filters (Y106, J129, H158, and F184), easily meeting the 0.5% Roman inter-filter calibration requirement for a cutting-edge measurement of cosmological parameters using SNe Ia. Simulated data in the bluer Z087 filter shows larger ~ 2–3 mmg biases, also meeting this requirement, but with more room for improvement. Our forward-model code has been released on Zenodo.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-123: Dosimetric Comparison Between Portrait and Landscape Orientations in Radiochromic Film Dosimetry

Purpose: To compare the dosimetric properties of radiochromic films with different orientation. Methods: A sheet of EBT3 film was cut into eight pieces with the following sizes: 15×15 cm2 (one piece), 5x15 cm{sup 2} (two) and 4×5 cm{sup 2} (five). A set of two EBT3 sheets was used at each dose level. Two sets were used changing the delivered doses (1 and 2 Gy). The 5×15 cm{sup 2} pieces were rotated by 90 degrees in relation to each other, such that one had landscape orientation and the other had portrait orientation. All 5×15 cm2 pieces were irradiated with their long side aligned with the x-axis of the radiation field. The 15×15 cm{sup 2} pieces were irradiated rotated at 90 degrees to each other. Five pieces, (a total of ten from two sheets) were used to obtain a calibration curve. The irradiated films were scanned using an Epson ES-2200 scanner and were analyzed using ImageJ software. In this study, no correction was applied for the nonuniform scanner signal that is evident in the direction of the scanner lamp. Each film piece was scanned both in portrait and landscape orientations. Dosimetric comparisons of the beam profiles were made in terms of the film orientations (portrait and landscape) and scanner bed directions (perpendicular and parallel to the scanner movement). Results: In general, portrait orientation exhibited higher noise than landscape and was adversely affected to a great extent by the nonuniformity in the direction of the scanner lamp. A significant difference in the measured field widths between the perpendicular and parallel directions was found for both orientations. Conclusion: Without correction for the nonuniform scanner signal in the direction of the scanner lamp, a landscape orientation is preferable. A more detailed investigation is planned to evaluate quantitatively the effect of orientation on the dosimetric properties of a film.

60 APPLIED LIFE SCIENCES↗

SU-F-T-186: A Treatment Planning Study of Normal Tissue Sparing with Robustness Optimized IMPT, 4Pi IMRT, and VMAT for Head and Neck Cases

Purpose: We performed a retrospective dosimetric comparison study between the robustness optimized Intensity Modulated Proton Therapy (RO-IMPT), volumetric-modulated arc therapy (VMAT), and the non-coplanar 4? intensity modulated radiation therapy (IMRT). These methods represent the most advanced radiation treatment methods clinically available. We compare their dosimetric performance for head and neck cancer treatments with special focus on the OAR sparing near the tumor volumes. Methods: A total of 11 head and neck cases, which include 10 recurrent cases and one bilateral case, were selected for the study. Different dose levels were prescribed to tumor target depending on disease and location. Three treatment plans were created on commercial TPS systems for a novel noncoplanar 4π method (20 beams), VMAT, and RO-IMPT technique (maximum 4 fields). The maximum patient positioning error was set to 3 mm and the maximum proton range uncertainty was set to 3% for the robustness optimization. Line dose profiles were investigated for OARs close to tumor volumes. Results: All three techniques achieved 98% coverage of the CTV target and most photon plans had less than 110% of the hot spots. The RO-IMPT plans show superior tumor dose homogeneity than 4? and VMAT plans. Although RO-IMPT has greater R50 dose spillage to the surrounding normal tissue than 4π and VMAT, the RO-IMPT plans demonstrate better or comparable OAR (parotid, mandible, carotid, oral cavity, pharynx, and etc.) sparing for structures closely abutting tumor targets. Conclusion: The RO-IMPT’s ability of OAR sparing is benchmarked against the C-arm linac based non-coplanar 4π technique and the standard VMAT method. RO-IMPT consistently shows better or comparable OAR sparing even for tissue structures closely abutting treatment target volume. RO-IMPT further reduces treatment uncertainty associated with proton therapy and delivers robust treatment plans to both unilateral and bilateral head and neck cancer patients with desirable treatment time.

60 APPLIED LIFE SCIENCES↗