Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
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 OSTI2016

Improvements on the minimax algorithm for the Laplace transformation of orbital energy denominators

We present a robust and non-heuristic algorithm that finds all extremum points of the error distribution function of numerically Laplace-transformed orbital energy denominators. The extremum point search is one of the two key steps for finding the minimax approximation. If pre-tabulation of initial guesses is supposed to be avoided, strategies for a sufficiently robust algorithm have not been discussed so far. We compare our non-heuristic approach with a bracketing and bisection algorithm and demonstrate that 3 times less function evaluations are required altogether when applying it to typical non-relativistic and relativistic quantum chemical systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
DOE OSTI2019

Positivity of the Spectral Densities of Retarded Floquet Green Functions

Periodically driven nonequilibrium many-body systems are interesting because they have quasi-energy spectra, which can be tailored by controlling the external driving fields. We derive the general spectral representation of retarded Green functions in the Floquet regime, thereby generalizing the well-known Lehmann representation from equilibrium many-body physics. Here, the derived spectral Floquet representation allows us to prove the non-negativity of spectral densities and to determine exact spectral sum rules, which can be employed to benchmark the accuracy of approximations to the exact Floquet many-body Green functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Permafrost microbial communities and functional genes are structured by latitudinal and soil geochemical gradients

Permafrost underlies approximately one quarter of Northern Hemisphere terrestrial surfaces and contains 25–50% of the global soil carbon (C) pool. Permafrost soils and the C stocks within are vulnerable to ongoing and future projected climate warming. The biogeography of microbial communities inhabiting permafrost has not been examined beyond a small number of sites focused on local-scale variation. Permafrost is different from other soils. Perennially frozen conditions in permafrost dictate that microbial communities do not turn over quickly, thus possibly providing strong linkages to past environments. Thus, the factors structuring the composition and function of microbial communities may differ from patterns observed in other terrestrial environments. Here, we analyzed 133 permafrost metagenomes from North America, Europe, and Asia. Permafrost biodiversity and taxonomic distribution varied in relation to pH, latitude and soil depth. The distribution of genes differed by latitude, soil depth, age, and pH. Genes that were the most highly variable across all sites were associated with energy metabolism and C-assimilation. Specifically, methanogenesis, fermentation, nitrate reduction, and replenishment of citric acid cycle intermediates. This suggests that adaptations to energy acquisition and substrate availability are among some of the strongest selective pressures shaping permafrost microbial communities. The spatial variation in metabolic potential has primed communities for specific biogeochemical processes as soils thaw due to climate change, which could cause regional- to global- scale variation in C and nitrogen processing and greenhouse gas emissions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The impacts of convex piecewise linear cost formulations on AC optimal power flow

Despite strong connections through shared application areas, research efforts on power market optimization (e.g., unit commitment) and power network optimization (e.g., optimal power flow) remain largely independent. A notable illustration of this is the treatment of power generation cost functions, where nonlinear network optimization has largely used polynomial representations and market optimization has adopted piecewise linear encodings. This work combines state-of-the-art results from both lines of research to understand the best mathematical formulations of the nonlinear AC optimal power flow problem with piecewise linear generation cost functions. An extensive numerical analysis of non-convex models, linear approximations, and convex relaxations across fifty-four realistic test cases illustrates that nonlinear optimization methods are surprisingly sensitive to the mathematical formulation of piecewise linear functions. The results indicate that a poor formulation choice can slow down algorithm performance by a factor of ten, increasing the runtime from seconds to minutes. Furthermore, these results provide valuable insights into the best formulations of nonlinear optimal power flow problems with piecewise linear cost functions, an important step towards building a new generation of energy markets that incorporate the nonlinear AC power flow model.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2017

Kinetic-energy release of fragments from electron-impact dissociation of the molecular hydrogen ion and its isotopologues

We calculate here the kinetic-energy release distributions of fragments produced for electron-impact dissociation of the vibrationally excited molecular hydrogen ion ${{\mathrm{H}}_{2}}^{+}$ and its isotopologues ${{\mathrm{D}}_{2}}^{+}$ and ${{\mathrm{T}}_{2}}^{+}$. Here we apply the adiabatic-nuclei convergent close-coupling method and compare results with several different methods, including the ${\delta}$ approximation. Results are presented for a number of dissociative excitation transitions and dissociative ionization as a function of the initial vibrational state of the molecule. We confirm that the square root approximation is a good approximation for the adiabatic-nuclei kinetic-energy release cross sections of ${{\mathrm{H}}_{2}}^{+}$. Agreement with experiment, where available, is good.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Assessment of electronic structure methods for the determination of the ground spin states of Fe( ii ), Fe( iii ) and Fe( iv ) complexes

Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe 2+ ion, gaseous FeO and 14 iron complexes using Kohn–Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe(II), five are Fe(III) and two are Fe(IV) complexes. These are calculated using 20 exchange–correlation functionals. In particular, we use a local spin density approximation (LSDA) – GVWN5, four generalized gradient approximations (GGAs) – BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) – GAM and N12, two meta-GGAs – M06-L and M11-L, a meta-NGA – MN15-L, five hybrid GGAs – B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs – M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA – MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe 2+ ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. Here, the local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state.

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

Activation of N 2 on Manganese Nitride-Supported Ni 3 and Fe 3 Clusters and Relevance to Ammonia Formation

Dual-site models were constructed to represent manganese nitride (Mn 4 N)-supported Ni 3 and Fe 3 clusters for NH 3 synthesis. Density functional theory calculations produced an energy barrier of approximately 0.55 eV for N–N bond activation at the interfacial nitrogen vacancy sites (N v ); also, the hydrogenation and removal of interfacial N is promoted by earth-abundant Ni and Fe metals. Steady-state microkinetic modeling revealed that the turnover frequencies of NH 3 production follow an order of Fe 3 @Mn 4 N ≈ Ni 3 @Mn 4 N > Mn 4 N > Fe >> Ni. Moreover, we present clear evidence that, before NH 3 formation, NH migrates from N v onto the metallic sites. Using N binding energy (BE N ) and the transition-state energy of N 2 activation (E TS ) as descriptors, we concluded that the beneficial effects owing to interfacial N v sites are the most pronounced when BE N is either too strong or too weak while E TS is high; otherwise, excessive N v sites may hinder catalyst performance.

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

Thermodynamics of Uranium Tri-Iodide from Density-Functional Theory

Density-functional theory (DFT) is employed to investigate the thermodynamic and ground-state properties of bulk uranium tri-iodide, UI3. The theory is fully relativistic and electron correlations, beyond the DFT and generalized gradient approximation, are addressed with orbital polarization. The electronic structure indicates anti-ferromagnetism, in agreement with neutron diffraction, with band gaps and a non-metallic system. Furthermore, the formation energy, atomic volume, crystal structure, and heat capacity are calculated in reasonable agreement with experiments, whereas for the elastic constants experimental data are unavailable for comparison. The thermodynamical properties are modeled within a quasi-harmonic approximation and the heat capacity and Gibbs free energy as functions of temperature agree with available calculation of phase diagram (CALPHAD) thermodynamic assessment of the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗