Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tracking errors”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

GPU Acceleration of Mean Free Path Based Kernel Density Estimators for Monte Carlo Neutronics Simulations

Kernel Density Estimators (KDEs) are a non-parametric density estimation technique that has recently been applied to Monte Carlo radiation transport simulations. Kernel density estimators are an alternative to histogram tallies for obtaining global solutions in Monte Carlo tallies. With KDEs, a single event, either a collision or particle track, can contribute to the score at multiple tally points with the uncertainty at those points being independent of the desired resolution of the solution. Thus, KDEs show potential for obtaining estimates of a global solution with reduced variance when compared to a histogram. Previously, KDEs have been applied to neutronics for one-group reactor physics problems and fixed source shielding applications. However, little work was done to obtain reaction rates using KDEs. This paper introduces a new form of the MFP KDE that is capable of handling general geometries. Furthermore, extending the MFP KDE to 2-D problems in continuous energy introduces inaccuracies to the solution. An ad-hoc solution to these inaccuracies is introduced that produces errors smaller than 4% at material interfaces.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Potential and limitations of satellite laser altimetry for monitoring water surface dynamics: ICESat for US lakes

Elevation measurements from the Ice, Cloud and Land Elevation Satellite (ICESat) have been applied to monitor dynamics of lakes and other surface water bodies. Despite such potential, the true utility of ICESat--more generally, satellite laser altimetry--for tracking surface water dynamics over time has not been adequately assessed, especially in the continental or global contexts. Here, we analyzed ICESat elevation data for the conterminous United States and examined the potential and limitations of satellite laser altimetry in measuring water-level dynamics. Owing to a lack of spatially-explicit ground-based water-level data, we first resorted to high-fidelity land elevation data acquired by airborne lidar to quantify ICESat’s ranging accuracy. We then performed trend and frequency analyses to evaluate how reliably ICESat could capture water-level dynamics over a range of temporal scales, as compared to in-situ gauge measurements. Our analyses showed that ICESat had a vertical ranging error of 0.16 m at the footprint level—a limit on the detectable range of water-level dynamics. The sparsity of data over time was identified as a major factor limiting the use of ICESat for water dynamics studies. Of all the US lakes, only 361 had quality ICESat measurements for more than two flight passes. Even for those lakes with sufficient temporal coverage, ICESat failed to capture the true interannual water-level dynamics in 68% of the cases. Our frequency analysis suggested that even with a repeat cycle of two months, ICESat could capture only 60% of the variations in water-level dynamics for at most 34 % of the US lakes. To capture 60% of the water-level variation for most of the US lakes, a weekly repeat cycle (e.g., less than 5 days) is needed – a requirement difficult to meet in current designs of spaceborne laser altimetry. Overall, our results highlight that current or near-future satellite laser missions, though with high ranging accuracies, are unlikely to fulfill the general needs in remotely monitoring water surface dynamics for lakes or reservoirs.

ICESat↗
DOE OSTI2012

A WASHINGTON PHOTOMETRIC SURVEY OF THE LARGE MAGELLANIC CLOUD FIELD STAR POPULATION

We present photometry for an unprecedented database of some 5.5 million stars distributed throughout the Large Magellanic Cloud main body, from 21 fields covering a total area of 7.6 deg{sup 2}, obtained from Washington CT{sub 1} T{sub 2} CTIO 4 m MOSAIC data. Extensive artificial star tests over the whole mosaic image data set and the observed behavior of the photometric errors with magnitude demonstrate the accuracy of the morphology and clearly delineate the position of the main features in the color-magnitude diagrams (CMDs). The representative T{sub 1}(MS TO) mags are on average {approx}0.5 mag brighter than the T{sub 1} mags for the 100% completeness level of the respective field, allowing us to derive an accurate age estimate. We have analyzed the CMD Hess diagrams and used the peaks in star counts at the main sequence turnoff and red clump (RC) locations to age date the most dominant sub-population (or 'representative' population) in the stellar population mix. The metallicity of this representative population is estimated from the locus of the most populous red giant branch track. We use these results to derive age and metallicity estimates for all of our fields. The analyzed fields span age and metallicity ranges covering most of the galaxy's lifetime and chemical enrichment, i.e., ages and metallicities between {approx}1 and 13 Gyr and {approx}-0.2 and -1.2 dex, respectively. We show that the dispersions associated with the mean ages and metallicities represent in general a satisfactory estimate of the age/metallicity spread ({approx}1-3 Gyr/0.2-0.3 dex), although a few subfields have a slightly larger age/metallicity spread. Finally, we revisit the study of the vertical structure (VS) phenomenon, a striking feature composed of stars that extend from the bottom, bluest end of the RC to {approx}0.45 mag fainter. We confirm that the VS phenomenon is not clearly seen in most of the studied fields and suggest that its occurrence is linked to some other condition(s) in addition to the appropriate age, metallicity, and the necessary red giant star density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Prediction of the Lorentz Force Detuning and pressure sensitivity for a Pillbox cavity

The Lorentz Force Detuning (LFD) and the pressure sensitivity are two critical concerns during the design of a Superconducting Radio Frequency (SRF) cavity resonator. The mechanical deformation of the bare Niobium cavity walls, due to the electromagnetic fields and fluctuation of the external pressure in the Helium bath, can dynamically and statically detune the frequency of the cavity and can cause beam phase errors. The frequency shift can be compensated by additional RF power, that is required to maintain the accelerating gradient, or by sophisticated tuning mechanisms and control-compensation algorithms. Passive stiffening is one of the simplest and most effective tools that can be used during the early design phase, capable of satisfying the Radio Frequency (RF) requisites. This approach requires several multiphysics simulations as well as a deep mechanical and RF knowledge of the phenomena involved. In this paper, is presented a new numerical model for a pillbox cavity that can predict the frequency shifts caused by the LFD and external pressure. This method allows to greatly reduce the computational effort, which is necessary to meet the RF requirements and to keep track of the frequency shifts without using the time consuming multiphysics simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

SU-G-JeP4-02: An Investigation of Respiratory Surrogate Motion Data Requirements for Multiple-Step Ahead Prediction

Purpose: Respiratory-gated radiotherapy and dynamic tracking employ real-time imaging and surrogate motion-monitoring methods with tumor motion prediction in advance of real-time. This study investigated respiratory motion data length on prediction accuracy of tumor motion. Methods: Predictions generated from the algorithm are validated against a one-dimensional surrogate signal of amplitude versus time. Prediction consists of three major components: extracting top-ranked subcomponents from training data matching the last respiratory cycle; calculating weighting factors from best-matched subcomponents; fusing data proceeding best-matched subcomponents with respective weighting factors to form predictions. Predictions for one respiratory cycle (∼3-6seconds) were assessed using 351 patient data from the respiratory management device. Performance was evaluated for correlation coefficient and root mean square error (RMSE) between prediction and final respiratory cycle. Results: Respiratory prediction results fell into two classes, where best predictions for 70 cycles or less performed using relative prediction and greater than 70 cycles are predicted similarly using relative and derivative relative. For 70 respiratory cycles or less, the average correlation between prediction and final respiratory cycle was 0.9999±0.0001, 0.9999±0.0001, 0.9988±0.0003, 0.9985±0.0023, and 0.9981±0.0023 with RMSE values of 0.0091±0.0030, 0.0091±0.0030, 0.0305±0.0051, 0.0299±0.0259, and 0.0299±0.0259 for equal, relative, pattern, derivative equal and derivative relative weighting methods, respectively. Respectively, the total best prediction for each method was 37, 65, 20, 22, and 22. For data with greater than 70 cycles average correlation was 0.9999±0.0001, 0.9999±0.0001, 0.9988±0.0004, 0.9988±0.0020, and 0.9988±0.0020 with RMSE values of 0.0081±0.0031, 0.0082±0.0033, 0.0306±0.0056, 0.0218±0.0222, and 0.0218±0.0222 for equal, relative, pattern, derivative equal and derivative relative weighting methods, respectively. Respectively, the total best prediction for each method was 24, 44, 42, 30, and 45. Conclusion: The prediction algorithms are effective in estimating surrogate motion in advance. These results indicate an advantage in using relative prediction for shorter data and either relative or derivative relative prediction for longer data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

An Unmanned Aerial System (UAS) for concurrent measurements of solar-induced chlorophyll fluorescence and hyperspectral reflectance toward improving crop monitoring

Unmanned aerial system (UAS)-based remote sensing can provide high spatial- and temporal-resolution crop monitoring for precision management. Existing crop monitoring UAS primarily use conventional techniques such as multi-spectral broadband vegetation indices (VIs) to track canopy structure/biomass. Recently developed lightweight hyperspectral sensors can track crop physiology and performance via complementary signals such as solar-induced chlorophyll fluorescence (SIF) and photochemical reflectance index (PRI) that can be derived from hyperspectral reflectance. However, few existing UAS can acquire high-quality SIF and hyperspectral reflectance. Here, we designed a novel UAS that simultaneously captures far-red SIF and hyperspectral reflectance, using a single bifurcated fiber and motorized arm to measure downwelling and upwelling irradiance. The UAS was tested over (1) a heterogeneous corn field with spatially varying crop yield potential and (2) a set of soybean and corn plots under differing nutrient treatments. Our UAS can maintain stability under low-to-moderate winds (2–3 m/s) with <0.4 m geolocation accuracy, <0.1 m altitude drift, and <1.5° error of the fiber optic from nadir position during scanning. Weekly seasonal campaigns reveal that SIF outperforms conventional VIs such as the normalized difference vegetation index (NDVI) for distinguishing plots with different crop yield potential, especially after canopy closure when NDVI tends to saturate. The UAS can capture the diurnal dynamics of SIF and PRI for both corn and soybean. At the seasonal scale, SIF acquired from the UAS correlates well with that from a fixed SIF tower system ( R 2 = 0.81), which measured the same target albeit with different footprints, demonstrating the capability of UAS in characterizing the seasonal progression of crop activity. In conclusion, our newly developed UAS provides high-quality SIF and hyperspectral reflectance, facilitating mechanistic understanding of the physiological control on photosynthesis dynamics. Our findings also imply that this UAS could enable extension of fixed tower-based systems to monitor multiple targets from diurnal to seasonal scales and improve monitoring of heterogeneous fields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Off-line tracking of series parameters in distribution systems using AMI data

Electric distribution systems have historically lacked measurement points, and equipment is often operated to its failure point, resulting in customer outages. The widespread deployment of sensors at the distribution level is enabling observability. This paper presents an off-line parameter value tracking procedure that takes advantage of the increasing number of measurement devices being deployed at the distribution level to estimate changes in series impedance parameter values over time. The tracking of parameter values enables non-diurnal and non-seasonal change to be flagged for investigation. The presented method uses an unbalanced Distribution System State Estimation (DSSE) and a measurement residual-based parameter estimation procedure. Measurement residuals from multiple measurement snapshots are combined in order to increase the effective local redundancy and improve the robustness of the calculations in the presence of measurement noise. Data from devices on the primary distribution system and from customer meters, via an AMI system, form the input data set. Results of simulations on the IEEE 13-Node Test Feeder are presented to illustrate the proposed approach applied to changes in series impedance parameters. A 5% change in series resistance elements can be detected in the presence of 2% measurement error when combining less than 1 day of measurement snapshots into a single estimate.

Williams, Tess L.↗
DOE OSTI2023

Evaluation of Uncertainties in the Anthropogenic SO 2 Emissions in the USA from the OMI Point Source Catalog

Satellite remote sensing is a promising method of monitoring emissions that may be missing in inventories, but the accuracy of these estimates is often not clear. We demonstrate here a comprehensive evaluation of errors in anthropogenic sulfur dioxide (SO 2 ) emission estimates from NASA’s OMI point source catalog for the contiguous US by comparing emissions from the catalog with high-quality emission inventory data over different dimensions including size of individual sources, aggregate vs individual source errors, and potential bias in individual source estimates over time. For sources that are included in the catalog, we find that errors in aggregate (sum of error for all included sources) are relatively low. Errors for individual sources in any given year can be substantial, however, with over- or underestimates in terms of total error ranging from –80 to 110 kt (roughly 10–90th percentile). In this study, we find that these errors are not necessarily random over time and that there can be consistently positive or negative biases for individual sources. We did not find any overall statistical relationship between the degree of isolation of a source and bias, either at a 40 or 70 km scales. For a sub-set of sources where inventory emissions over a radius of 70 km around an OMI detection are larger than twice the emissions within 40 km, the OMI value is consistently overestimated. We find, as expected, that emission sources not included in the catalog are the largest aggregate source of difference between the satellite estimates and inventories, especially in more recent years where source emission magnitudes have been decreasing and note that trends in satellite detections do not necessarily track trends in total emissions. We find that the OMI-based SO 2 emissions are accurate in aggregate, when summed over a number of sources, but must be interpreted more cautiously at the individual source level. Similar analyses would be valuable for other satellite emission estimates; however, in many cases, the appropriate high-quality reference data may need to be generated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

In vivo evaluation of an adaptive resuscitation controller using whole blood and crystalloid infusates for hemorrhagic shock

Introduction Hemorrhage remains the leading cause of preventable death on the battlefield. The most effective means to increase survivability is early hemorrhage control and fluid resuscitation. Unfortunately, fluid resuscitation requires constant adjustments to ensure casualty is properly managed, which is often not feasible in the pre-hospital setting. In this study, we showed how an adaptive closed-loop controller for hemorrhage resuscitation can be used to automate hemodynamic management using a swine hemorrhagic shock injury model. Methods The adaptive resuscitation controller (ARC) was previously developed to track pressure–volume responsiveness in real time and adjust its infusion rate to reach the target mean arterial pressure (MAP). Swine while maintained under a surgical plane of anesthesia and analgesia underwent a splenectomy, followed by two hemorrhage and resuscitation events. For the first resuscitation event, hemorrhage was induced to reduce the MAP to 35 mmHg until arterial lactate reached 4 mmol/L. The ARC system then infused whole blood (WB) to reach the target MAP and maintained the subject using crystalloids for 120 min. For the second resuscitation event, the subjects were hemorrhaged again but resuscitated using only crystalloid infusion to reach the target MAP and 120-min maintenance. Results The ARC was effective at WB resuscitation, reaching the target MAP in 2.0 ± 1.0 min. The median performance error was 1.1% ± 4.6%, and target overshoot was 14.4% ± 7.0% of the target MAP. The ARC maintained all animals throughout the 120 min maintenance period. For the second crystalloid-based resuscitation, ARC required a longer time to reach the target MAP, at an average rise time of 4.3 ± 4.0 min. However, target overshoot was reduced to 8.4% ± 7.3% of the target MAP. Much higher flow rates were required to maintain the target MAP during the second resuscitation event than during the first resuscitation event. Discussion The ARC was able to rapidly reach and maintain the target MAP effectively. However, this sometimes required large volumes of fluid as the ARC’s only goal was to reach the target MAP. Further clinical insight is needed regarding the preferred aggression level to achieve the target MAP. In conclusion, the ARC was successful in its programmed objective of reaching and maintaining the target MAP for extended periods of time in vivo , a critical next step toward improving hemorrhage treatment in the pre-hospital environment.

Snider, Eric J.↗
DOE OSTI2015

SU-C-213-04: Application of Depth Sensing and 3D-Printing Technique for Total Body Irradiation (TBI) Patient Measurement and Treatment Planning

Purpose: To develop and validate an innovative method of using depth sensing cameras and 3D printing techniques for Total Body Irradiation (TBI) treatment planning and compensator fabrication. Methods: A tablet with motion tracking cameras and integrated depth sensing was used to scan a RANDOTM phantom arranged in a TBI treatment booth to detect and store the 3D surface in a point cloud (PC) format. The accuracy of the detected surface was evaluated by comparison to extracted measurements from CT scan images. The thickness, source to surface distance and off-axis distance of the phantom at different body section was measured for TBI treatment planning. A 2D map containing a detailed compensator design was calculated to achieve uniform dose distribution throughout the phantom. The compensator was fabricated using a 3D printer, silicone molding and tungsten powder. In vivo dosimetry measurements were performed using optically stimulated luminescent detectors (OSLDs). Results: The whole scan of the anthropomorphic phantom took approximately 30 seconds. The mean error for thickness measurements at each section of phantom compare to CT was 0.44 ± 0.268 cm. These errors resulted in approximately 2% dose error calculation and 0.4 mm tungsten thickness deviation for the compensator design. The accuracy of 3D compensator printing was within 0.2 mm. In vivo measurements for an end-to-end test showed the overall dose difference was within 3%. Conclusion: Motion cameras and depth sensing techniques proved to be an accurate and efficient tool for TBI patient measurement and treatment planning. 3D printing technique improved the efficiency and accuracy of the compensator production and ensured a more accurate treatment delivery.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Motion as a perturbation: Measurement-guided dose estimates to moving patient voxels during modulated arc deliveries

Purpose: To present a framework for measurement-guided VMAT dose reconstruction to moving patient voxels from a known motion kernel and the static phantom data, and to validate this perturbation-based approach with the proof-of-principle experiments. Methods: As described previously, the VMAT 3D dose to a static patient can be estimated by applying a phantom measurement-guided perturbation to the treatment planning system (TPS)-calculated dose grid. The fraction dose to any voxel in the presence of motion, assuming the motion kernel is known, can be derived in a similar fashion by applying a measurement-guided motion perturbation. The dose to the diodes in a helical phantom is recorded at 50 ms intervals and is transformed into a series of time-resolved high-density volumetric dose grids. A moving voxel is propagated through this 4D dose space and the fraction dose to that voxel in the phantom is accumulated. The ratio of this motion-perturbed, reconstructed dose to the TPS dose in the phantom serves as a perturbation factor, applied to the TPS fraction dose to the similarly situated voxel in the patient. This approach was validated by the ion chamber and film measurements on four phantoms of different shape and structure: homogeneous and inhomogeneous cylinders, a homogeneous cube, and an anthropomorphic thoracic phantom. A 2D motion stage was used to simulate the motion. The stage position was synchronized with the beam start time with the respiratory gating simulator. The motion patterns were designed such that the motion speed was in the upper range of the expected tumor motion (1-1.4 cm/s) and the range exceeded the normally observed limits (up to 5.7 cm). The conformal arc plans for X or Y motion (in the IEC 61217 coordinate system) consisted of manually created narrow (3 cm) rectangular strips moving in-phase (tracking) or phase-shifted by 90 Degree-Sign (crossing) with respect to the phantom motion. The XY motion was tested with the computer-derived VMAT MLC sequences. For all phantoms and plans, time-resolved (10 Hz) ion chamber dose was collected. In addition, coronal (XY) films were exposed in the cube phantom to a VMAT beam with two different starting phases, and compared to the reconstructed motion-perturbed dose planes. Results: For the X or Y motions with the moving strip and geometrical phantoms, the maximum difference between perturbation-reconstructed and ion chamber doses did not exceed 1.9%, and the average for any motion pattern/starting phase did not exceed 1.3%. For the VMAT plans on the cubic and thoracic phantoms, one point exhibited a 3.5% error, while the remaining five were all within 1.1%. Across all the measurements (N = 22), the average disagreement was 0.5 {+-} 1.3% (1 SD). The films exhibited {gamma}(3%/3 mm) passing rates {>=}90%. Conclusions: The dose to an arbitrary moving voxel in a patient can be estimated with acceptable accuracy for a VMAT delivery, by performing a single QA measurement with a cylindrical phantom and applying two consecutive perturbations to the TPS-calculated patient dose. The first one accounts for the differences between the planned and delivered static doses, while the second one corrects for the motion.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TH-C-BRD-07: Minimizing Dose Uncertainty for Spot Scanning Beam Proton Therapy of Moving Tumor with Optimization of Delivery Sequence

Purpose: Intensity modulated proton therapy (IMPT) has been shown to be able to reduce dose to normal tissue compared to intensity modulated photon radio-therapy (IMRT), and has been implemented for selected lung cancer patients. However, respiratory motion-induced dose uncertainty remain one of the major concerns for the radiotherapy of lung cancer, and the utility of IMPT for lung patients was limited because of the proton dose uncertainty induced by motion. Strategies such as repainting and tumor tracking have been proposed and studied but repainting could result in unacceptable long delivery time and tracking is not yet clinically available. We propose a novel delivery strategy for spot scanning proton beam therapy. Method: The effective number of delivery (END) for each spot position in a treatment plan was calculated based on the parameters of the delivery system, including time required for each spot, spot size and energy. The dose uncertainty was then calculated with an analytical formula. The spot delivery sequence was optimized to maximize END and minimize the dose uncertainty. 2D Measurements with a detector array on a 1D moving platform were performed to validate the calculated results. Results: 143 2D measurements on a moving platform were performed for different delivery sequences of a single layer uniform pattern. The measured dose uncertainty is a strong function of the delivery sequence, the worst delivery sequence results in dose error up to 70% while the optimized delivery sequence results in dose error of <5%. END vs. measured dose uncertainty follows the analytical formula. Conclusion: With optimized delivery sequence, it is feasible to minimize the dose uncertainty due to motion in spot scanning proton therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-112: Intensity-Based Pulmonary Image Registration: An Evaluation Study

Purpose: Accurate alignment of thoracic CT images is essential for dose tracking and to safely implement adaptive radiotherapy in lung cancers. At the same time it is challenging given the highly elastic nature of lung tissue deformations. The objective of this study was to assess the performances of three state-of-art intensity-based algorithms in terms of their ability to register thoracic CT images subject to affine, barrel, and sinusoid transformation. Methods: Intensity similarity measures of the evaluated algorithms contained sum-of-squared difference (SSD), local mutual information (LMI), and residual complexity (RC). Five thoracic CT scans obtained from the EMPIRE10 challenge database were included and served as reference images. Each CT dataset was distorted by realistic affine, barrel, and sinusoid transformations. Registration performances of the three algorithms were evaluated for each distortion type in terms of intensity root mean square error (IRMSE) between the reference and registered images in the lung regions. Results: For affine distortions, the three algorithms differed significantly in registration of thoracic images both visually and nominally in terms of IRMSE with a mean of 0.011 for SSD, 0.039 for RC, and 0.026 for LMI (p<0.01; Kruskal-Wallis test). For barrel distortion, the three algorithms showed nominally no significant difference in terms of IRMSE with a mean of 0.026 for SSD, 0.086 for RC, and 0.054 for LMI (p=0.16) . A significant difference was seen for sinusoid distorted thoracic CT data with mean lung IRMSE of 0.039 for SSD, 0.092 for RC, and 0.035 for LMI (p=0.02). Conclusion: Pulmonary deformations might vary to a large extent in nature in a daily clinical setting due to factors ranging from anatomy variations to respiratory motion to image quality. It can be appreciated from the results of the present study that the suitability of application of a particular algorithm for pulmonary image registration is deformation-dependent.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Depth-resolved registration of transesophageal echo to x-ray fluoroscopy using an inverse geometry fluoroscopy system

Purpose: Image registration between standard x-ray fluoroscopy and transesophageal echocardiography (TEE) has recently been proposed. Scanning-beam digital x-ray (SBDX) is an inverse geometry fluoroscopy system designed for cardiac procedures. This study presents a method for 3D registration of SBDX and TEE images based on the tomosynthesis and 3D tracking capabilities of SBDX. Methods: The registration algorithm utilizes the stack of tomosynthetic planes produced by the SBDX system to estimate the physical 3D coordinates of salient key-points on the TEE probe. The key-points are used to arrive at an initial estimate of the probe pose, which is then refined using a 2D/3D registration method adapted for inverse geometry fluoroscopy. A phantom study was conducted to evaluate probe pose estimation accuracy relative to the ground truth, as defined by a set of coregistered fiducial markers. This experiment was conducted with varying probe poses and levels of signal difference-to-noise ratio (SDNR). Additional phantom and in vivo studies were performed to evaluate the correspondence of catheter tip positions in TEE and x-ray images following registration of the two modalities. Results: Target registration error (TRE) was used to characterize both pose estimation and registration accuracy. In the study of pose estimation accuracy, successful pose estimates (3D TRE < 5.0 mm) were obtained in 97% of cases when the SDNR was 5.9 or higher in seven out of eight poses. Under these conditions, 3D TRE was 2.32 ± 1.88 mm, and 2D (projection) TRE was 1.61 ± 1.36 mm. Probe localization error along the source-detector axis was 0.87 ± 1.31 mm. For the in vivo experiments, mean 3D TRE ranged from 2.6 to 4.6 mm and mean 2D TRE ranged from 1.1 to 1.6 mm. Anatomy extracted from the echo images appeared well aligned when projected onto the SBDX images. Conclusions: Full 6 DOF image registration between SBDX and TEE is feasible and accurate to within 5 mm. Future studies will focus on real-time implementation and application-specific analysis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within {approx}30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H{sub 2}O{sup +}+ H{sub 2}O {yields} OH + H{sub 3}O{sup +}. The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM/MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with {approx}40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H{sub 2}O{sup +}{sub (aq)} species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Measurement of the Neutron Electromagnetic Form Factor Ratio at High Momentum Transfer

The inner structure of the nucleon (proton and neutron) remains a topic of great interest in nuclear and particle physics, after many decades of study. For example, understanding the quark-gluon dynamics inside the nucleon would shed light on how 99% of the nucleon mass is created. The neutron electromagnetic form factors, Gn E and Gn M , give important insights into the neutron structure. The Super BigBite Spectrometer (SBS) program at Jefferson Lab (JLab) seeks to extend the form factor measurements for both the proton and the neutron. The neutron electric form actor, Gn E , has been historically difficult to measure due to the short lifetime of the free neutron and the small value of Gn E . The GEn-II experiment is part of the SBS program and seeks to measure Gn E , significantly increasing the high momentum transfer coverage. A newly designed polarized 3He target increased the figure of merit by three times compared to previous measurements. The analysis of this data is especially challenging due to the unprecedented high-rate environment caused by the open nature of the spectrometer with a direct line of sight to the target. This required developing new Gas Electron Multiplier (GEM) particle trackers which can cover large areas demanded by this setup and handle particle rates up to 500 kHz/cm2. Rates this high over a large area is unprecedented in particle tracking systems and came with a number of challenges. Data taken in the SBS program was critical to understanding hardware and software solutions that improved the track reconstruction efficiency to be >97% with a position resolution of 70 ?m. In previous experiments the proton electromagnetic form factors, Gp E and Gp M were measured up to Q2 = 8.5 GeV2 and Q2 = 30 GeV2, respectively, while Gn E has only been measured up to Q2 = 3.4 GeV2. The GEn-II experiment has measured the neutron form factor ratio, Gn E/Gn M, at Q2 values of 2.90, 6.50, and 9.47 GeV2 by scattering a polarized electron beam with a polarized 3He target, used here as an effective polarized neutron target, and measuring the double spin asymmetry of the cross section. Previous Gn E measurements do not extend above Q2 = 3.4 GeV2, and therefore this analysis has extended the world data by almost three times. The background correction is especially difficult at the higher Q2 settings leading to large systematic errors. As very exploratory results from this early analysis of the data, we find for Q2 = 2.90 GeV2, Gn E = 0.0157 ±stat 0.0016 ±sys 0.0011, for Q2 = 6.50 GeV2, Gn E = 0.0067 ±stat 0.0019 ±sys 0.0005, and for Q2 = 9.46 GeV2, Gn E = 0.0046 ±stat 0.0023 ±sys 0.0005. These results are compared to predictions from the Dyson-Schwinger Equations (DSE) model and a Relativistic Constituent Quark Model (RCQM).

Jeffas, Sean↗
DOE OSTI2015

Geometric validation of self-gating k-space-sorted 4D-MRI vs 4D-CT using a respiratory motion phantom

Purpose: MRI is increasingly being used for radiotherapy planning, simulation, and in-treatment-room motion monitoring. To provide more detailed temporal and spatial MR data for these tasks, we have recently developed a novel self-gated (SG) MRI technique with advantage of k-space phase sorting, high isotropic spatial resolution, and high temporal resolution. The current work describes the validation of this 4D-MRI technique using a MRI- and CT-compatible respiratory motion phantom and comparison to 4D-CT. Methods: The 4D-MRI sequence is based on a spoiled gradient echo-based 3D projection reconstruction sequence with self-gating for 4D-MRI at 3 T. Respiratory phase is resolved by using SG k-space lines as the motion surrogate. 4D-MRI images are reconstructed into ten temporal bins with spatial resolution 1.56 × 1.56 × 1.56 mm{sup 3}. A MRI-CT compatible phantom was designed to validate the performance of the 4D-MRI sequence and 4D-CT imaging. A spherical target (diameter 23 mm, volume 6.37 ml) filled with high-concentration gadolinium (Gd) gel is embedded into a plastic box (35 × 40 × 63 mm{sup 3}) and stabilized with low-concentration Gd gel. The phantom, driven by an air pump, is able to produce human-type breathing patterns between 4 and 30 respiratory cycles/min. 4D-CT of the phantom has been acquired in cine mode, and reconstructed into ten phases with slice thickness 1.25 mm. The 4D images sets were imported into a treatment planning software for target contouring. The geometrical accuracy of the 4D MRI and CT images has been quantified using target volume, flattening, and eccentricity. The target motion was measured by tracking the centroids of the spheres in each individual phase. Motion ground-truth was obtained from input signals and real-time video recordings. Results: The dynamic phantom has been operated in four respiratory rate (RR) settings, 6, 10, 15, and 20/min, and was scanned with 4D-MRI and 4D-CT. 4D-CT images have target-stretching, partial-missing, and other motion artifacts in various phases, whereas the 4D-MRI images are visually free of those artifacts. Volume percentage difference for the 6.37 ml target ranged from 5.3% ± 4.3% to 10.3% ± 5.9% for 4D-CT, and 1.47 ± 0.52 to 2.12 ± 1.60 for 4D-MRI. With an increase of respiratory rate, the target volumetric and geometric deviations increase for 4D-CT images while remaining stable for the 4D-MRI images. Target motion amplitude errors at different RRs were measured with a range of 0.66–1.25 mm for 4D-CT and 0.2–0.42 mm for 4D-MRI. The results of Mann–Whitney tests indicated that 4D-MRI significantly outperforms 4D-CT in phase-based target volumetric (p = 0.027) and geometric (p < 0.001) measures. Both modalities achieve equivalent accuracy in measuring motion amplitude (p = 0.828). Conclusions: The k-space self-gated 4D-MRI technique provides a robust method for accurately imaging phase-based target motion and geometry. Compared to 4D-CT, the current 4D-MRI technique demonstrates superior spatiotemporal resolution, and robust resistance to motion artifacts caused by fast target motion and irregular breathing patterns. The technique can be used extensively in abdominal targeting, motion gating, and toward implementing MRI-based adaptive radiotherapy.

60 APPLIED LIFE SCIENCES↗