Search NASA⌕ Search

SEARCH · Search NASA

Results for “ERROR ANALYSIS”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

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

Rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions for cometary and planetary applications, as predicted by mixed quantum-classical theory

We present new calculations of collision cross sections for state-to-state transitions between the rotational states in an H 2 O + H 2 O system, which are used to generate a new database of collisional rate coefficients for cometary and planetary applications. Calculations were carried out using a mixed quantum-classical theory approach that is implemented in the code MQCT. The large basis set of rotational states used in these calculations permits us to predict thermally averaged cross sections for 441 transitions in para- and ortho-H 2 O in a broad range of temperatures. It is found that all state-to-state transitions in the H 2 O + H 2 O system split into two well-defined groups, one with higher cross-section values and lower energy transfer, which corresponds to the dipole-dipole driven processes. The other group has smaller cross sections and higher energy transfer, driven by higher-order interaction terms. We present a detailed analysis of the theoretical error bars, and we symmetrized the state-to-state transition matrixes to ensure that excitation and quenching processes for each transition satisfy the principle of microscopic reversibility. We also compare our results with other data available from the literature for H 2 O + H 2 O collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Density functional theory of water with the machine-learned DM21 functional

The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn–Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. Furthermore, we find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.

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

Picosecond magneto-optic thermometry measurements of nanoscale thermal transport in AlN thin films

The thermal conductivity Λ of wide bandgap semiconductor thin films, such as AlN, affects the performance of high-frequency devices, power devices, and optoelectronics. However, accurate measurements of Λ in thin films with sub-micrometer thicknesses and Λ > 100 W m -1 K -1 is challenging. Widely used pump/probe metrologies, such as time–domain thermoreflectance (TDTR) and frequency–domain thermoreflectance, lack the spatiotemporal resolution necessary to accurately quantify thermal properties of sub-micrometer thin films with high Λ. In this work, we use a combination of magneto-optic thermometry and TiN interfacial layers to significantly enhance the spatiotemporal resolution of pump/probe thermal transport measurements. We use our approach to measure Λ of 100, 400, and 1000 nm AlN thin films. We coat AlN thin films with a ferromagnetic thin-film transducer with the geometry of (1 nm-Pt/0.4 nm-Co) x3 /(2 nm-TiN). This PtCo/TiN transducer has a fast thermal response time of <50 ps, which allows us to differentiate between the thermal response of the transducer, AlN thin film, and substrate. For the 100, 400, and 1000 nm thick AlN films, we determine Λ to be 200 ± 80, 165 ± 35, and 300 ± 70 W m -1 K -1 , respectively. We conclude with an uncertainty analysis that quantifies the errors associated with pump/probe measurements of thermal conductivity, as a function of transducer type, thin-film thermal conductivity, and thin-film thickness. Time resolved magneto-optic Kerr effect experiments can measure films that are three to five times thinner than is possible with standard pump/probe metrologies, such as TDTR. This advance in metrology will enable better characterization of nanoscale heat transfer in high thermal conductivity material systems like wide bandgap semiconductor heterostructures and devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Moments of nucleon isovector structure functions in 2+1+1-flavor QCD

We present results on the isovector momentum fraction, $\langle x \rangle_{u–d}$ , helicity moment, $\langle x \rangle_{Δ u–Δd}$ , and the transversity moment, $\langle x \rangle_{δu–δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC Collaboration using 2 + 1 + 1 -flavors of dynamical highly improved staggered quarks. The correlation functions are calculated using the Wilson-Clover action, and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI'–MOM scheme. The data have been collected at lattice spacings a ≈ 0.15 , 0.12, 0.09, and 0.06 fm and $M_π$ ≈ 310 , 220, and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{\text{MS}}$ scheme at 2 GeV, are $\langle$ x $\rangle$ u – d = 0.173 ( 14 ) ( 07 ) , $\langle x \rangle_{Δu–Δd}$ = 0.213 ( 15 ) ( 22 ) , and $\langle x \rangle_{δu–δd}$ = 0.208 ( 19 ) ( 24 ) , where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

MIONet: Learning Multiple-Input Operators via Tensor Product

As an emerging paradigm in scientific machine learning, neural operators aim to learn operators, via neural networks, that map between infinite-dimensional function spaces. Several neural operators have been recently developed. However, all the existing neural operators are only designed to learn operators defined on a single Banach space; i.e., the input of the operator is a single function. Here, for the first time, we study the operator regression via neural networks for multiple-input operators defined on the product of Banach spaces. We first prove a universal approximation theorem of continuous multiple-input operators. We also provide a detailed theoretical analysis including the approximation error, which provides guidance for the design of the network architecture. Based on our theory and a low-rank approximation, we propose a novel neural operator, MIONet, to learn multiple-input operators. MIONet consists of several branch nets for encoding the input functions and a trunk net for encoding the domain of the output function. Here, we demonstrate that MIONet can learn solution operators involving systems governed by ordinary and partial differential equations. In our computational examples, we also show that we can endow MIONet with prior knowledge of the underlying system, such as linearity and periodicity, to further improve accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Composition Measurements of the LAW Phase 5 Glasses

This report provides the results from the chemical analyses of glass compositions for the Low-Activity Waste Phase 5 study glasses, a series of simulated nuclear waste glasses fabricated at Pacific Northwest National Laboratory. These data will be used in the development of enhanced property/composition models for waste glass vitrification at Hanford. Chemical analyses were performed on a representative sample of each of the quenched glasses to allow for comparisons with targeted compositions. The relative differences between the targeted and measured concentrations of CaO, Cl - , K 2 O, Li 2 O, and SO 3 for some of the glasses were greater than 10%. Several of these glass samples were rerun to verify original measurements; however, there were no significant changes in the measurements to indicate errors in preparation or analysis of the samples. These results can be used in further characterization of this series of glasses, including the normalization of Product Consistency Test results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Characterization of Recovery Human Action Mechanisms in Nuclear Power Plants

Recovery human action is defined as the action that prevents deviant conditions from producing unwanted effects. Analyzing recovery actions has been a critical part in human reliability analysis (HRA). However, there are a couple of limitations to treat recovery actions only depending on the current HRA methods. Representatively, the existing recovery analysis does not specifically consider recovery actions as are occurred in actual nuclear power plants (NPPs). The overall goal of this study aims to develop a novel recovery analysis method to account for human action recoveries in context of scenarios as well as complement the limitations of existing recovery analysis. In this paper, the recovery analysis in current HRA methods and their challenges are introduced. A strategy to achieve the goal is introduced with a modified recovery definition. Then, how we have researched the approach will be introduced in the paper.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Implementing Spatio-Temporal 3D-Convolution Neural Networks and UAV Time Series Imagery to Better Predict Lodging Damage in Sorghum

Unmanned aerial vehicle (UAV)-based remote sensing is gaining momentum in a variety of agricultural and environmental applications. Very-high-resolution remote sensing image sets collected repeatedly throughout a crop growing season are becoming increasingly common. Analytical methods able to learn from both spatial and time dimensions of the data may allow for an improved estimation of crop traits, as well as the effects of genetics and the environment on these traits. Multispectral and geometric time series imagery was collected by UAV on 11 dates, along with ground-truth data, in a field trial of 866 genetically diverse biomass sorghum accessions. We compared the performance of Convolution Neural Network (CNN) architectures that used image data from single dates (two spatial dimensions, 2D) versus multiple dates (two spatial dimensions + temporal dimension, 3D) to estimate lodging detection and severity. Lodging was detected with 3D-CNN analysis of time series imagery with 0.88 accuracy, 0.92 Precision, and 0.83 Recall. This outperformed the best 2D-CNN on a single date with 0.85 accuracy, 0.84 Precision, and 0.76 Recall. The variation in lodging severity was estimated by the best 3D-CNN analysis with 9.4% mean absolute error (MAE), 11.9% root mean square error (RMSE), and goodness-of-fit (R2) of 0.76. This was a significant improvement over the best 2D-CNN analysis with 11.84% MAE, 14.91% RMSE, and 0.63 R2. The success of the improved 3D-CNN analysis approach depended on the inclusion of “before and after” data, i.e., images collected on dates before and after the lodging event. The integration of geometric and spectral features with 3D-CNN architecture was also key to the improved assessment of lodging severity, which is an important and difficult-to-assess phenomenon in bioenergy feedstocks such as biomass sorghum. This demonstrates that spatio-temporal CNN architectures based on UAV time series imagery have significant potential to enhance plant phenotyping capabilities in crop breeding and Precision agriculture applications.

3D-convolution neural networks↗
DOE OSTI2021

Erratum: “On virial analysis at low aspect ratio” [Phys. Plasmas 23, 072508 (2016)]

This article corrects an error in M.W. Bongard et al., 'On Virial Analysis at Low Aspect Ratio,' Phys. Plasmas 23, 072508 (2016) pertaining to a typographical error in its equation 21. Furthermore, all numerical calculations reported in this paper utilized the correct formulation described here. Accordingly, no changes to the conclusions of the original article are required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analyzing rotational bands in odd-mass nuclei using effective field theory and Bayesian methods

We recently developed an Effective Field Theory (EFT) for rotational bands in odd-mass nuclei. Here we use EFT expressions to perform a Bayesian analysis of data on the rotational energy levels of 99 Tc, 155,157 Gd, 159 Dy, 167,169 Er, 167,169 Tm, 183 W, 235 U and 239 Pu. The error model in our Bayesian analysis includes both experimental and EFT truncation uncertainties. It also accounts for the fact that low-energy constants (LECs) at even and odd orders are expected to have different sizes. We use Markov Chain Monte Carlo (MCMC) sampling to explore the joint posterior of the EFT and error-model parameters and show both the LECs and the breakdown scale can be reliably determined. We extract the LECs up to fourth order in the EFT and find that, provided we correctly account for EFT truncation errors in our likelihood, results for lower-order LECs are stable as we go to higher orders. LEC results are also stable with respect to the addition of higher-energy data. We extract the expansion parameter for all the nuclei listed above and find a clear correlation between the extracted and the expected value of the inverse breakdown scale, W, based on the single-particle and vibrational energy scales. However, the W that actually determines the convergence of the EFT expansion is markedly smaller than would be naively expected based on those scales.

EFT↗
DOE OSTI2021

Antiferromagnetic character of the quantum phase transition in the Hubbard model on the honeycomb lattice

In this work we provide a unified, comprehensive treatment of all operators that contribute to the antiferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We use the Hybrid Monte Carlo algorithm to perform a systematic, carefully controlled analysis in the temporal Trotter error and of the thermodynamic limit. We expect our findings to improve the consistency of Monte Carlo determinations of critical exponents. We perform a data collapse analysis and determine the critical exponent β = 0.898 (37) for the semimetal-Mott insulator transition in the hexagonal Hubbard Model. Our methods are applicable to a wide range of lattice theories of strongly correlated electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

VISILIENCE: An Interactive Visualization Framework for Resilience Analysis using Control-Flow Graph

Soft errors have become one of the major concerns for the error resilience of HPC applications, as those errors can cause HPC applications to generate serious outcomes such as Silent Data Corruptions (SDCs). A large body of approaches has been proposed to analyze the resilience of HPC applications. However, existing studies rarely address the challenges of the analysis result perception. Specifically, resilience analysis techniques often produce a massive volume of unstructured data, making it difficult for programmers to conduct the resilience analysis due to non-intuitive raw data. Furthermore, different analysis models produce diverse results with multiple levels of details, which may create hurdles to compare and explore the resilience of HPC program execution. To this end, we present VISILIENCE, an interactive VISual resILIENCE analysis framework to allow programmers to facilitate the resilience analysis of HPC applications. In particular, VISILIENCE leverages an effective visualization approach Control Flow Graph (CFG) to present a function execution. In addition, three widely-used models for resilience analysis (i.e., Y-Branch, IPAS, and TRIDENT) are seamlessly embedded into the framework for resilience analysis and result comparison. Multiple case studies have been conducted to demonstrate the effectiveness of our proposed framework VISILIENCE.

Jiang, Hailong↗
DOE OSTI2020

Infrared spectroscopic method for uranium isotopic analysis

A high performance infrared (HPIR) system was developed and demonstrated for the infrared absorption analysis of {sup 235}U and {sup 238}U isotopes in uranium hexafluoride gas samples. Sweeping the quantum cascade laser light source over the spectral range and sampling via a high-rate analog-to-digital converter provided 0.0005 cm-1 spectral resolution, which allowed for high-precision measurements of the isotopic peak shift. A data analysis method was developed using principal component analysis to predict the isotope weight % content of {sup 235}U. The HPIR precision, accuracy, and error were evaluated for a wide range of isotopic ratio samples (0.287 – 93.7 weight % {sup 235}U), and the results were compared to the International Target Values (ITVS) set forth by the International Atomic Energy Agency (IAEA) for non-destructive and destructive analytical techniques. The method meets or surpasses the IAEA ITVs for non-destructive analysis of samples with isotopic content of depleted to highly enriched. The results also demonstrated the capability of the HPIR system to correctly predict the {sup 235}U weight % content of a mislabeled sample whose isotopic distribution was validated by mass spectroscopic measurements. The HPIR measurement is nondestructive and, thus, allows for confirmatory analyses of the exact sample at a designated IAEA lab if higher-resolution or a certified analysis is needed. (author)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Distributed Berkeley Efficient Long-Read to Long-Read Aligner and Overlapper (DiBELLA) v1.0.0

We present a parallel algorithm and scalable implementation for genome analysis, specifically the problem of finding overlaps and alignments for data from "third generation" long read sequencers. While long sequences of DNA offer enormous advantages for biological analysis and insight, current long read sequencing instruments have high error rates and therefore require different approaches to analysis than their short read counterparts. Our work focuses on an efficient distributed-memory parallelization of an accurate single-node algorithm for overlapping and aligning long reads. We achieve scalability of this irregular algorithm by addressing the competing issues of increasing parallelism, minimizing communication, constraining the memory footprint, and ensuring good load balance. The resulting application, DiBELLA, is the first distributed memory overlapper and aligner specifically designed for long reads and parallel scalability.

Ellis, Marquita↗
DOE OSTI2023

Applications of the Dulmage–Mendelsohn decomposition for debugging nonlinear optimization problems

Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.

42 ENGINEERING↗
DOE OSTI2023

BeyondPlanck: XII. Cosmological parameter constraints with end-to-end error propagation

We present cosmological parameter constraints estimated using the Bayesian BEYONDPLANCK analysis framework. This method supports seamless end-to-end error propagation from raw time-ordered data onto final cosmological parameters. As a first demonstration of the method, we analyzed time-ordered Planck LFI observations, combined with selected external data (WMAP 33–61 GHz, Planck HFI DR4 353 and 857 GHz, and Haslam 408 MHz) in the form of pixelized maps that are used to break critical astrophysical degeneracies. Overall, all the results are generally in good agreement with previously reported values from Planck 2018 and WMAP, with the largest relative difference for any parameter amounting about 1σ when considering only temperature multipoles between 30 ≤ ℓ ≤ 600. In cases where there are differences, we note that the BEYONDPLANCK results are generally slightly closer to the high-ℓ HFI-dominated Planck 2018 results than previous analyses, suggesting slightly less tension between low and high multipoles. Using low-ℓ polarization information from LFI and WMAP, we find a best-fit value of τ = 0.066 ± 0.013, which is higher than the low value of τ = 0.052 ± 0.008 derived from Planck 2018 and slightly lower than the value of 0.069 ± 0.011 derived from the joint analysis of official LFI and WMAP products. Most importantly, however, we find that the uncertainty derived in the BEYONDPLANCK processing is about 30 % greater than when analyzing the official products, after taking into account the different sky coverage. We argue that this uncertainty is due to a marginalization over a more complete model of instrumental and astrophysical parameters, which results in more reliable and more rigorously defined uncertainties. We find that about 2000 Monte Carlo samples are required to achieve a robust convergence for a low-resolution cosmic microwave background (CMB) covariance matrix with 225 independent modes, and producing these samples takes about eight weeks on a modest computing cluster with 256 cores.

79 ASTRONOMY AND ASTROPHYSICS↗