Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function approximation”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Block Krylov Subspace Methods for Functions of Matrices II: Modified Block FOM

We analyze an expansion of the generalized block Krylov subspace framework of [Electron. Trans. Numer. Anal., 47 (2017), pp. 100--126]. This expansion allows the use of low-rank modifications of the matrix projected onto the block Krylov subspace and contains, as special cases, the block GMRES method and the new block Radau--Arnoldi method. Within this general setting, we present results that extend the interpolation property from the nonblock case to a matrix polynomial interpolation property for the block case, and we relate the eigenvalues of the projected matrix to the latent roots of these matrix polynomials. Some error bounds for these modified block FOM methods for solving linear systems are presented. We then show how cospatial residuals can be preserved in the case of families of shifted linear block systems. This result is used to derive computationally practical restarted algorithms for block Krylov approximations that compute the action of a matrix function on a set of several vectors simultaneously. Finally, we prove some error bounds and present numerical results showing that two modifications of FOM, the block harmonic and the block Radau--Arnoldi methods for matrix functions, can significantly improve the convergence behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

First-principles calculation of excited states of diatomic molecules: a benchmark for the Gutzwiller conjugate gradient minimisation method

We recently proposed the Gutzwiller conjugate gradient minimisation (GCGM) method for efficient and accurate calculation of the ground state total energy of molecular and bulk systems. The GCGM method is developed under the framework of Gutzwiller wave function but goes beyond the commonly adopted Gutzwiller approximation to improve the accuracy and flexibility in treating the correlation effects. In this conference proceeding, we benchmark the GCGM method with the calculation of excited state potential energy curves of three diatomic molecules, namely H 2 , N 2 , and O 2 . Our calculations demonstrate the flexibility and reasonable accuracy of the method.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

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

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗
DOE OSTI2020

Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments

The new Bayesian Monte Carlo (MC) evaluation framework described in this abstract has been conceived as an attempt to improve nuclear data evaluations of differential crosssection data by removing the following two approximations conventionally employed for nuclear data evaluations: all probability density functions (PDFs) of all data and model parameters, both prior and posterior, are assumed to be normal (i.e., Gaussian) PDFs, and all uncertainties and covariances are propagated using a linear approximation. With these approximations removed, the Bayesian MC (BMC) framework could be used to account for nonlinear effects and would enable improved evaluations of differential cross sections and IBE data that are presently performed based on the assumptions itemized above. The BMC would also improve upon the uniform sampling of IBE parameters from within ranges defined by their evaluated uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

The quantum Hall effect at a weak magnetic field

Using a weak limit for the hopping integral in one direction in the Hofstadter model, we show that the fermion states in the gaps of the spectrum are determined within the Kitaev chain. The proposed approach allows us to study the behavior of Chern insulators (CI) in different classes of symmetry. We consider the Hofstadter model on the square and honeycomb lattices in the case of rational and irrational magnetic fluxes ϕ, and discuss the behavior of the Hall conductance at a weak magnetic field in a sample of finite size. We show that in the semiclassical limit at the center of the fermion spectrum, the Bloch states of fermions turn into chiral Majorana fermion liquid when the magnetic scale (1/ϕ) is equal to the sample size N. We are talking about the dielectric–metal phase transition, which is determined by the behavior of the Landau levels in 2D fermion systems in a transverse magnetic field. When a magnetic scale, which determines the wave function of fermions, exceeds the size of the sample, a jump in the longitudinal conductance occurs. The wave function describes non-localized states of fermions, the sample becomes a conductor, the system changes from the dielectric state to the metallic one. It is shown, that at 1∕ϕ>N the quantum Hall effect and the Landau levels are not realized, which makes possibility to study the behavior of CI in irrational magnetic fluxes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

On the validity of the guiding-center approximation in the presence of strong magnetic gradients

The motion of a charged particle in a nonuniform straight magnetic field with a constant magnetic-field gradient is solved exactly in terms of elliptic functions. The connection between this problem and the guiding-center approximation is discussed. Here, it is shown that, for this problem, the predictions of higher-order guiding-center theory agree very well with the orbit-averaged particle motion and hold well beyond the standard guiding-center limit $\epsilon$ $\equiv$ p/L << 1, where p is the gyromotion length scale and L is the magnetic-field gradient length scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

The Potential of DAS in Teleseismic Studies: Insights From the Goldstone Experiment

Distributed acoustic sensing (DAS) is a recently developed technique that has demonstrated its utility in the oil and gas industry. Here we demonstrate the potential of DAS in teleseismic studies using the Goldstone OpticaL Fiber Seismic experiment in Goldstone, California. By analyzing teleseismic waveforms from the 10 January 2018 M7.5 Honduras earthquake recorded on ~5,000 DAS channels and the nearby broadband station GSC, we first compute receiver functions for DAS channels using the vertical-component GSC velocity as an approximation for the incident source wavelet. The Moho P -to- s conversions are clearly visible on DAS receiver functions. Here, we then derive meter-scale arrival time measurements along the entire 20-km-long array. We are also able to measure path-averaged Rayleigh wave group velocity and local Rayleigh wave phase velocity. The latter, however, has large uncertainties. Our study suggests that DAS will likely play an important role in many fields of passive seismology in the near future.

58 GEOSCIENCES↗
DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Nuclear kinetic density from ab initio theory

The nuclear kinetic density is one of many fundamental, nonobservable quantities in density functional theory (DFT) dependent on the nonlocal nuclear density. Often, approximations may be made when computing the density that may result in spurious contributions in other DFT quantities. With the ability to compute the nonlocal nuclear density from ab initio wave functions, it is now possible to estimate effects of such spurious contributions. Herein, we derive the kinetic density using ab initio nonlocal scalar one-body nuclear densities computed within the no-core shell model (NCSM) approach, utilizing two- and three-nucleon chiral interactions as the sole input. The ability to compute translationally invariant nonlocal densities allows us to gauge the impact of the spurious center-of-mass (c.m.) contributions in DFT quantities, such as the kinetic density, and provide ab initio insight into refining energy density functionals. The nonlocal nuclear densities are derived from the NCSM one-body densities calculated in second quantization. We present a review of c.m. contaminated and translationally invariant nuclear densities. We then derive an analytic expression for the kinetic density using these nonlocal densities, producing an ab initio kinetic density. The ground-state nonlocal densities of He 4 , 6 , 8 , C 12 , and O 16 are used to compute the kinetic densities of the aforementioned nuclei. The impact of c.m. removal techniques in the density are discussed and compared to a procedure applied in DFT. The results of this work can be extended to other fundamental quantities in DFT. The use of a general nonlocal density allows for the calculation of fundamental quantities taken as input in theories such as DFT. This allows benchmarking c.m. removal procedures and provides a bridge for comparison between ab initio and DFT many-body techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Plasma Heating in Solar Microflares: Statistics and Analysis

In this paper we present the results of an analysis of 481 weak solar flares, from A0.01 class flares to the B GOES class, that were observed during the period of extremely low solar activity from 2009 April to July. For all flares we measured the temperature of the plasma in the isothermal and two-temperature approximations and tried to fit its relationship with the X-ray class using exponential and power-law functions. We found that the whole temperature distribution in the range from A0.01 to X-class cannot be fit by one exponential function. The fitting for weak flares below A1.0 is significantly steeper than that for medium and large flares. The power-law approximation seems to be more reliable: the corresponding functions were found to be in good agreement with experimental data both for microflares and for normal flares. Our study predicts that evidence of plasma heating can be found in flares starting from the A0.0002 X-ray class. Weaker events presumably cannot heat the surrounding plasma. We also estimated emission measures for all flares studied and the thermal energy for 113 events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Calculation of the Physicochemical Characteristics of a New Orthorhombic Form of Diamond

A geometrically optimized structure and some properties of a new diamond polymorph have been assessed using density-functional theory calculations in the local-density approximation. The new phase has a primitive orthorhombic lattice (Pmma). The calculation results indicate that the density, sublimation energy, bulk modulus, Vickers hardness, and band gap of this phase are lower than those of cubic diamond by 5–18%. The orthorhombic diamond polymorph can be obtained by compressing graphite in the [001] direction at a pressure of 59.4 GPa. For this phase transition, we calculated theoretical X-ray diffraction patterns of the new phase and intermediate states, which can be used for their identification in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Revisiting Understanding of Electrochemical CO 2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory

Copper (Cu) electrodes, as the most efficacious of CO 2 reduction reaction (CO 2 RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO 2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO 2 RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO 2 RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO 2 RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at –0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Furthermore, our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2018

Photodissociative Cross-Linking of Non-covalent Peptide-Peptide Ion Complexes in the Gas Phase

We report a gas-phase UV photodissociation study investigating non-covalent interactions between neutral hydrophobic pentapeptides and peptide ions incorporating a diazirine-tagged photoleucine residue. Phenylalanine (Phe) and proline (Pro) were chosen as the conformation-affecting residues that were incorporated into a small library of neutral pentapeptides. Gas-phase ion-molecule complexes of these peptides with photo-labeled pentapeptides were subjected to photodissociation. Selective photocleavage of the diazirine ring at 355 nm formed short-lived carbene intermediates that underwent cross-linking by insertion into H–X bonds of the target peptide. The cross-link positions were established from collision-induced dissociation tandem mass spectra (CID-MS{sup 3}) providing sequence information on the covalent adducts. Effects of the amino acid residue (Pro or Phe) and its position in the target peptide sequence were evaluated. For proline-containing peptides, interactions resulting in covalent cross-links in these complexes became more prominent as proline was moved towards the C-terminus of the target peptide sequence. The photocross-linking yields of phenylalanine-containing peptides depended on the position of both phenylalanine and photoleucine. Density functional theory calculations were used to assign structures of low-energy conformers of the (GLPMG + GLL*LK + H){sup +} complex. Born-Oppenheimer molecular dynamics trajectory calculations were used to capture the thermal motion in the complexes within 100 ps and determine close contacts between the incipient carbene and the H–X bonds in the target peptide. This provided atomic-level resolution of potential cross-links that aided spectra interpretation and was in agreement with experimental data. .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Electronic Structure Engine

INQ is an engine for electronic structure calculations. INQ is implemented as a C++ library of components designed to build a parallel, portable code for electronic structure calculations, including ground state and excited states. INQ is based on density functional theory (DFT). It can calculate ground state properties in DFT, and also excited states using time-dependent DFT (TDDFT) in real-time and linear-response. It implements different approximations for the exchange and correlation part: semi-local functionals like LDAs, GGAs and metaGGAs, and hybrid functionals that are implemented through the ACE approach for fast execution.

ANDRADE-VALENCIA, XAVIERI↗
DOE OSTI2019

Electronic Structure and Quadrupole Interactions in Promising Cathode Materials Na{sub x}M{sub y}(MoO{sub 4}){sub 3}, M = Mn, Fe, Co, and Ni

The electronic structure and the magnetic properties of molybdates Na{sub x}M{sub y}(MoO{sub 4}){sub 3} (M = Mn, Fe, Co, and Ni) which are promising materials for sodium batteries have been studied in the framework of the density functional theory with the GGA and GGA+U approximations for the first time. The calculations show that all the compounds are insulators. An important role of the correlation effects, provided by the on-site Coulomb interactions, was established in the formation of the band gap in these compounds. The quadrupole constants of {sup 23}Na nuclei are calculated in the nonmagnetic and ferromagnetic states within the GGA and GGA+U approaches. It is shown that the quadrupole frequencies for nonequivalent crystallographic positions of sodium are in different frequency ranges, which allows to study the diffusion of sodium in these compounds by the Nuclear Magnetic Resonance method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning

Density functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we, in this study, use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO 3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ~1.5 eV for this class of materials, where ~1 eV comes from downward shifting the valence band and ~0.5 eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗