Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Detecting macroevolutionary genotype–phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype–phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω C , which measures the error-corrected convergence rate of protein evolution. While ω C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype–phenotype associations, even in lineages that diverged for hundreds of millions of years.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Accelerating error correction in tomographic reconstruction

Abstract Spurred by recent advances in detector technology and X-ray optics, upgrades to scanning-probe-based tomographic imaging have led to an exponential growth in the amount and complexity of experimental data and have created a clear opportunity for tomographic imaging to approach single-atom sensitivity. The improved spatial resolution, however, is highly susceptible to systematic and random experimental errors, such as center of rotation drifts, which may lead to imaging artifacts and prevent reliable data extraction. Here, we present a model-based approach that simultaneously optimizes the reconstructed specimen and sinogram alignment as a single optimization problem for tomographic reconstruction with center of rotation error correction. Our algorithm utilizes an adaptive regularizer that is dynamically adjusted at each alternating iteration step. Furthermore, we describe its implementation in a software package targeting high-throughput workflows for execution on distributed-memory clusters. We demonstrate the performance of our solver on large-scale synthetic problems and show that it is robust to a wide range of noise and experimental drifts with near-ideal throughput.

Ali, Sajid (ORCID:0000000321864636)↗
DOE OSTI2016

Kohn–Sham calculations of NMR shifts for paramagnetic 3d metal complexes: protocols, delocalization error, and the curious amide proton shifts of a high-spin iron(ii) macrocycle complex

A theory for the nuclear chemical shifts of molecules in arbitrary spin states is applied to a set of paramagnetic organometallic complexes of 3d metals. Ligand chemical shifts are calculated and analyzed using Kohn–Sham (KS) density functional theory with and without relativistic corrections. The roles of the KS delocalization error, Gaussian-type versus Slater-type basis sets, relativistic effects (scalar and spin–orbit), and zero field splitting (ZFS) are investigated. A strong functional dependence of the chemical shifts is apparent and correlated with the delocalization error. Here, the functional dependence is between one and two orders of magnitude larger than variations of the NMR shifts due the other influences that are investigated. ZFS effects are negligible in the determination of the NMR chemical shifts of the complexes except at very low temperatures. The DFT calculated shifts agree reasonably well with experiment. A 73 ppm difference in the NMR shifts of the two protons in the amide groups of a high-spin Fe(II) macrocycle complex arises from selective O → Fe dative bonding that only involves the transfer of β spin density, along with orbital delocalization throughout the ligand bonding framework which electronically couples the coordinating oxygen lone pair orbitals directly to the amide trans proton.

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

Delocalization error poisons the density-functional many-body expansion

The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion–water interactions, typified by F − (H 2 O) N with N ≳ 15. We attribute these oscillations to self-interaction error in the density-functional approximation. The effect is minor or negligible in small water clusters, explaining why it has not been noticed previously, but grows to catastrophic proportion in clusters that are only moderately larger. This behavior can be counteracted with hybrid functionals but only if the fraction of exact exchange is ≳50%, whereas modern meta-generalized gradient approximations including ωB97X-V, SCAN, and SCAN0 are insufficient to eliminate divergent behavior. Other mitigation strategies including counterpoise correction, density correction (i.e., exchange–correlation functionals evaluated atop Hartree–Fock densities), and dielectric continuum boundary conditions do little to curtail the problematic oscillations. In contrast, energy-based screening to cull unimportant subsystems can successfully forestall divergent behavior. These results suggest that extreme caution is warranted when the many-body expansion is combined with density functional theory.

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

Informed total-error-minimizing priors: Interpretable cosmological parameter constraints despite complex nuisance effects

While Bayesian inference techniques are standard in cosmological analyses, it is common to interpret resulting parameter constraints with a frequentist intuition. This intuition can fail, for example, when marginalizing high-dimensional parameter spaces onto subsets of parameters, because of what has come to be known as projection effects or prior volume effects. We present the method of informed total-error-minimizing (ITEM) priors to address this problem. An ITEM prior is a prior distribution on a set of nuisance parameters, such as those describing astrophysical or calibration systematics, intended to enforce the validity of a frequentist interpretation of the posterior constraints derived for a set of target parameters (e.g., cosmological parameters). Our method works as follows. For a set of plausible nuisance realizations, we generate target parameter posteriors using several different candidate priors for the nuisance parameters. We reject candidate priors that do not accomplish the minimum requirements of bias (of point estimates) and coverage (of confidence regions among a set of noisy realizations of the data) for the target parameters on one or more of the plausible nuisance realizations. Of the priors that survive this cut, we select the ITEM prior as the one that minimizes the total error of the marginalized posteriors of the target parameters. As a proof of concept, we applied our method to the density split statistics measured in Dark Energy Survey Year 1 data. We demonstrate that the ITEM priors substantially reduce prior volume effects that otherwise arise and that they allow for sharpened yet robust constraints on the parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Quantum error correcting codes and 4-dimensional arithmetic hyperbolic manifolds

Using 4-dimensional arithmetic hyperbolic manifolds, we construct some new homological quantum error correcting codes. They are low density parity check codes with linear rate and distance n{sup ε}. Their rate is evaluated via Euler characteristic arguments and their distance using Z{sub 2}-systolic geometry. This construction answers a question of Zémor [“On Cayley graphs, surface codes, and the limits of homological coding for quantum error correction,” in Proceedings of Second International Workshop on Coding and Cryptology (IWCC), Lecture Notes in Computer Science Vol. 5557 (2009), pp. 259–273], who asked whether homological codes with such parameters could exist at all.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Density-functional errors in ionization potential with increasing system size

This work investigates the effects of molecular size on the accuracy of density-functional ionization potentials for a set of 28 hydrocarbons, including series of alkanes, alkenes, and oligoacenes. As the system size increases, delocalization error introduces a systematic underestimation of the ionization potential, which is rationalized by considering the fractional-charge behavior of the electronic energies. The computation of the ionization potential with many density-functional approximations is not size-extensive due to excessive delocalization of the incipient positive charge. While inclusion of exact exchange reduces the observed errors, system-specific tuning of long-range corrected functionals does not generally improve accuracy. These results emphasize that good performance of a functional for small molecules is not necessarily transferable to larger systems.

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

Three sources of errors in the Ehrenfest treatment of inelastic scattering and possible ways of resolving them

In order to identify the origin of possible errors in the mixed quantum/classical approach to inelastic scattering [A. Semenov and D. Babikov, J. Chem. Phys. 140, 044306 (2014) and A. Semenov, M.-L. Dubernet, and D. Babikov, J. Chem. Phys. 141, 114304 (2014)], a simplified model is considered that consists of one intermolecular degree of freedom and two intramolecular states, coupled by a simple potential. For this system, analytic derivations are carried out to determine (i) the exact quantum mechanical solution of the inelastic scattering problem, (ii) a simplified version of it with all oscillatory terms neglected, and (iii) the Ehrenfest solution in which the translational motion is described by the mean-field trajectory while the internal molecular motion is treated by the time-dependent Schrodinger equation. It is shown that the appropriate choice of velocity for the mean-field trajectory permits to enforce microscopic reversibility and gives results in excellent agreement with full-quantum results. The average velocity method of Billing is rigorously derived as a limiting case (of this more general approach), when reversibility is enforced locally, at the initial moment of time only. It is demonstrated that errors of state-to-state transition probabilities in the Ehrenfest approach occur at lower values of total energy E if the magnitudes of excitation energy ΔE, potential energy difference between the two states ΔV, and coupling of two states V12 are large. Possible ways of applying this concept to rotational transitions in real molecules are explored, using examples from CO + CO inelastic scattering.

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

Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C

With the aim of mitigating the basis set error in density functional theory (DFT) calculations employing local basis sets, we herein develop two empirical corrections for basis set superposition error (BSSE) in the def2-SVPD basis, a basis which - when stripped of BSSE - is capable of providing near-complete-basis DFT results for non-covalent interactions. Specifically, we adapt the existing pairwise geometrical counterpoise (gCP) approach to the def2-SVPD basis, and we develop a beyond-pairwise approach, DFT-C, which we parameterize across a small set of intermolecular interactions. Both gCP and DFT-C are evaluated against the traditional Boys-Bernardi counterpoise correction across a set of 3402 non-covalent binding energies and isomerization energies. We find that the DFT-C method represents a significant improvement over gCP, particularly for non-covalently-interacting molecular clusters. Furthermore, DFT-C is transferable among density functionals and can be combined with existing functionals - such as B97M-V - to recover large-basis results at a fraction of the cost.

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

A neoclassical drift-magnetohydrodynamical fluid model of the interaction of a magnetic island chain with a resonant error-field in a high temperature tokamak plasma

A two-fluid, neoclassical theory of the interaction of a single magnetic island chain with a resonant error-field in a quasi-cylindrical, low-β, tokamak plasma is presented. The plasmas typically found in large hot tokamaks lie in the so-called weak neoclassical flow-damping regime in which the neoclassical ion stress tensor is not the dominant term in the ion parallel equation of motion. Nevertheless, flow-damping in such plasmas dominates ion perpendicular viscosity, and is largely responsible for determining the phase velocity of a freely rotating island chain (which is in the ion diamagnetic direction relative to the local E × B frame at the rational surface). The critical vacuum island width required to lock the island chain is mostly determined by the ion neoclassical poloidal flow damping rate at the rational surface. The stabilizing effect of the average field-line curvature, as well as the destabilizing effect of the perturbed bootstrap current, is the same for a freely rotating, a non-uniformly rotating, and a locked island chain. The destabilizing effect of the error-field averages to zero when the chain is rotating and only manifests itself when the chain locks. The perturbed ion polarization current has a small destabilizing effect on a freely rotating island chain, but a large destabilizing effect on both a non-uniformly rotating and a locked island chain. Finally, this behavior may account for the experimentally observed fact that locked island chains are much more unstable than corresponding freely rotating chains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Subtlety of TiO 2 phase stability: Reliability of the density functional theory predictions and persistence of the self-interaction error

TiO 2 is an important material with broad applications that can exist in different phases with dramatically different properties. Theoretical prediction of their polymorph energetics is therefore critical for the material design and for identifying thermodynamically accessible structures. Determining TiO 2 relative phase stabilities remains challenging for first-principles methods, and density functional theory is the only approach available for studying phase stabilities at finite temperatures with acceptable computational efficiency. Here, we show that density functional theory using the recently developed efficient strongly constrained and appropriately normed (SCAN) [Sun et al., Phys. Rev. Lett. 115, 036402 (2015)] exchange-correlation functional for the first time predicts the phase stability in qualitative agreement with the experimental results at realistic conditions. Further analysis shows that the self-interaction error intrinsic in the density functional persists in the stability prediction. By correcting the self-interaction error through an empirical approach, SCAN predicts the relative stability as well as defect properties in excellent agreement with the experimental results.

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

Rates of adsorption and desorption: Entropic contributions and errors due to mean-field approximations

We have performed exact classical rate calculations to compute adsorption and desorption rate constants with a model representative of a real system. We compute the desorption rate using transition-state theory by taking the dividing-surface far from the surface of the solid. We find that using a mean-field assumption, i.e., applying potential of mean force to transition state theory, could lead to two orders-of-magnitude errors in the rate constant owing to large fluctuations in the desorption barrier. Furthermore, we compute the adsorption rate by including a dynamical factor which reflects the probability of sticking to the solid surface. We find that the sticking probability is highly sensitive to the coverage. Also, we find that the adsorption rate computed from the mean-field assumption is not very different from the exact adsorption rate. We also compute entropic contribution to desorption rates and compare it to that obtained from two limiting models of adsorption—2D ideal gas and 2D ideal lattice gas. We show that at high temperatures (700 K), the entropic contribution to desorption rates computed from the exact calculations is very close to that obtained from the 2D ideal gas model. However, for lower to intermediate temperatures from 200 K to 500 K, the entropic contributions cover a wide range which lies in between the two limiting models and could lead to over two-orders-of-magnitude errors in the rate coefficient.

Activation energies↗
DOE OSTI2019

The effect of self-interaction error on electrostatic dipoles calculated using density functional theory

Spurious electron self-interaction in density functional approximations (DFAs) can lead to inaccurate predictions of charge transfer in heteronuclear molecules that manifest as errors in calculated electrostatic dipoles. Here we show the magnitude of these errors on dipoles computed for a diverse set of 47 molecules taken from the recent benchmark study of Hait and Head-Gordon. We compare results of Perdew-Wang local spin density approximation (PW92), Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and strongly constrained and appropriatedly normed (SCAN) meta-GGA dipole calculations, along with those of their respective self-interaction-corrected (SIC) counterparts, to reference values from accurate wave function-based methods. The SIC calculations were carried out using the Fermi-Löwdin orbital (FLO-SIC) approach. We find that correcting for self-interaction generally increases the degree of charge transfer, thereby increasing the size of calculated dipole moments. The FLO-SIC-PW92 and FLO-SIC-PBE dipoles are in better agreement with reference values than their uncorrected DFA counterparts, particularly for strongly ionic molecules where significant improvement is seen. Applying FLO-SIC to SCAN does not improve dipole values overall. In addition, we show that removing self-interaction improves the description of the dipole for stretched-bond geometries and recovers the physically correct separated atom limit of zero dipole. Finally, we find that the best agreement between the FLO-SIC-DFA and reference dipoles occurs when the molecular geometries are optimized using the FLO-SIC-DFA.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Sources of error and methods to improve accuracy in interface state density analysis using quasi-static capacitance–voltage measurements in wide bandgap semiconductors

Characterizing interface trap states in commercial wide bandgap devices using frequency-based measurements requires unconventionally high probing frequencies to account for both fast and slow traps associated with wide bandgap materials. The C - ψs technique has been suggested as a viable quasi-static method for determining the interface trap state densities in wide bandgap systems, but the results are shown to be susceptible to errors in the analysis procedure. Here this work explores the primary sources of errors present in the C - ψs technique using an analytical model that describes the apparent response for wide bandgap MOS capacitor devices. Measurement noise is shown to greatly impact the linear fitting routine of the 1/C$^{*2}_s$ vs Ψ s plot to calibrate the additive constant in the surface potential/gate voltage relationship, and an inexact knowledge of the oxide capacitance is also shown to impede interface trap state analysis near the band edge. In addition, a slight nonlinearity that is typically present throughout the 1/C$^{*2}_s$ vs Ψ s plot hinders the accurate estimation of interface trap densities, which is demonstrated for a fabricated n-SiC MOS capacitor device. Methods are suggested to improve quasi-static analysis, including a novel method to determine an approximate integration constant without relying on a linear fitting routine.

42 ENGINEERING↗
DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

TYPE Ia SUPERNOVA DISTANCE MODULUS BIAS AND DISPERSION FROM K-CORRECTION ERRORS: A DIRECT MEASUREMENT USING LIGHT CURVE FITS TO OBSERVED SPECTRAL TIME SERIES

We estimate systematic errors due to K-corrections in standard photometric analyses of high-redshift Type Ia supernovae. Errors due to K-correction occur when the spectral template model underlying the light curve fitter poorly represents the actual supernova spectral energy distribution, meaning that the distance modulus cannot be recovered accurately. In order to quantify this effect, synthetic photometry is performed on artificially redshifted spectrophotometric data from 119 low-redshift supernovae from the Nearby Supernova Factory, and the resulting light curves are fit with a conventional light curve fitter. We measure the variation in the standardized magnitude that would be fit for a given supernova if located at a range of redshifts and observed with various filter sets corresponding to current and future supernova surveys. We find significant variation in the measurements of the same supernovae placed at different redshifts regardless of filters used, which causes dispersion greater than ∼0.05 mag for measurements of photometry using the Sloan-like filters and a bias that corresponds to a 0.03 shift in w when applied to an outside data set. To test the result of a shift in supernova population or environment at higher redshifts, we repeat our calculations with the addition of a reweighting of the supernovae as a function of redshift and find that this strongly affects the results and would have repercussions for cosmology. We discuss possible methods to reduce the contribution of the K-correction bias and uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PULSAR TIMING ERRORS FROM ASYNCHRONOUS MULTI-FREQUENCY SAMPLING OF DISPERSION MEASURE VARIATIONS

Free electrons in the interstellar medium cause frequency-dependent delays in pulse arrival times due to both scattering and dispersion. Multi-frequency measurements are used to estimate and remove dispersion delays. In this paper, we focus on the effect of any non-simultaneity of multi-frequency observations on dispersive delay estimation and removal. Interstellar density variations combined with changes in the line of sight from pulsar and observer motions cause dispersion measure (DM) variations with an approximately power-law power spectrum, augmented in some cases by linear trends. We simulate time series, estimate the magnitude and statistical properties of timing errors that result from non-simultaneous observations, and derive prescriptions for data acquisition that are needed in order to achieve a specified timing precision. For nearby, highly stable pulsars, measurements need to be simultaneous to within about one day in order for the timing error from asynchronous DM correction to be less than about 10 ns. We discuss how timing precision improves when increasing the number of dual-frequency observations used in DM estimation for a given epoch. For a Kolmogorov wavenumber spectrum, we find about a factor of two improvement in precision timing when increasing from two to three observations but diminishing returns thereafter.

79 ASTRONOMY AND ASTROPHYSICS↗