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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A Data-Driven Approach for High-Impedance Fault Localization in Distribution Systems: Preprint

Accurate and quick identification of high-impedance faults (HIFs) is critical for the reliable operation of distribution systems. Unlike other faults in power grids, HIFs are very difficult to detect by conventional overcurrent relays due to the low fault current. Although HIFs can be affected by various factors, the voltage-current characteristics can substantially imply how the system responds to the disturbance and thus provides opportunities to effectively localize HIFs. In this work, we propose a data-driven approach for the identification of HIF events. To tackle the nonlinearity of the voltage-current trajectory, first, we formulate optimization problems to approximate the trajectory with piecewise functions. Then we collect the function features of all segments as inputs and use the support vector machine approach to efficiently identify HIFs at different locations. Numerical studies on the IEEE 123-node test feeder demonstrate the validity and accuracy of the proposed approach for real-time HIF identification.

explainable artificial intelligence↗

DOE OSTI2023

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Characterizing Peatland Microtopography Using Gradient and Microform-Based Approaches

Abstract Peatlands represent an important component of the global carbon cycle, storing 180–621 Gt of carbon (C). Small-scale spatial variations in elevation, frequently referred to as microtopography, influence ecological processes associated with the peatland C cycle, including Sphagnum photosynthesis and methane flux. Microtopography can be characterized with measures of topographic variability and by using conceptual classes (microforms) linked to function: most commonly hummocks and hollows. However, the criteria used to define these conceptual classes are often poorly described, if at all, and vary between studies. Such inconsistencies compel development of explicit quantitative methods to classify microforms. Furthermore, gradient-based characterizations that describe spatial variability without the use of microforms are lacking in the literature. Therefore, the objectives of this study were to (1) calculate peatland microtopographical elevation gradients and measures of spatial variability, (2) develop three microform classification methods intended for specific purposes, and (3) evaluate and contrast classification methods. Our results suggest that at spatial scales much larger than microforms, elevation distributions are unimodal and are well approximated with parametric probability density functions. Results from classifications were variable between methods and years and exhibited significant differences in mean hollow areal coverages of a raised ombrotrophic bog. Our results suggest that the conceptualization and classification of microforms can significantly influence microtopographic structural metrics. The three explicit methods for microform classification described here may be used and built upon for future applications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Estimates of Quantum Tunneling Effects for Hydrogen Diffusion in PuO 2

We detail the estimation of activation energies and quantum nuclear vibrational tunneling effects for hydrogen diffusion in PuO 2 based on Density Functional Theory calculations and a quantum double well approximation. We find that results are relatively insensitive to choice of exchange correlation functional. In addition, the representation of spin in the system and use of an extended Hubbard U correction has only a small effect on hydrogen point defect formation energies when the PuO 2 lattice is held fixed at the experimental density. We then compute approximate activation energies for transitions between hydrogen interstitial sites seeded by a semi-empirical quantum model and determine the quantum tunneling enhancement relative to classical kinetic rates. Our model indicates that diffusion rates in H/PuO 2 systems could be enhanced by more than one order of magnitude at ambient conditions and that these effects persist at high temperature. The method we propose here can be used as a fast screening tool for assessing possible quantum nuclear vibrational effects in any number of condensed phase materials and surfaces, where hydrogen hopping tends to follow well defined minimum energy pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Structural relaxation of water during rapid cooling from ambient temperatures

Experiments investigating the properties of deeply supercooled liquid water are needed to develop a comprehensive understanding of water’s anomalous properties. One approach involves transiently heating nanoscale water films into the supercooled region for several nanoseconds at a time and then interrogating the water films after they have quenched to cryogenic temperatures. To relate the results obtained with this approach to other experiments and simulations on supercooled water, it is important to understand how closely the quenched structure tracks the (metastable) equilibrium structure of water as a function of the transient heating temperature. A key step involves quantifying the extent to which water that is transiently heated to ambient temperatures [hyperquenched water (HQW)] subsequently relaxes toward the structure of low-density amorphous (LDA) ice as it cools. We analyzed the infrared reflection–absorption spectra of LDA, HQW, and crystalline ice films to determine their complex indices of refraction. With this information, we estimate that HQW retains ~50%–60% of a structural motif characteristic of water at high temperatures with the balance comprised of a low-temperature motif. This result, along with results from x-ray diffraction experiments on water and amorphous ices, allows one to quantify the fraction of the high-temperature motif at approximately zero pressure as a function of temperature from 150 to 350 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Orbital contributions in the element-resolved valence electronic structure of Bi 2 Se 3

In this work, we studied the bulk band structure of a topological insulator (TI) Bi 2 Se 3 and determined the contributions of the Bi and Se orbital states to the valence bands using standing-wave excited hard x-ray photoemission spectroscopy (SW-HAXPES). This SW technique can provide the element-resolved information and extract individual Bi and Se contributions to the Bi2Se3 valence band. Comparisons with density-functional theory calculations (local density approximation and GW) reveal that the Bi 6s, Bi 6p, and Se 4p states are dominant in the Bi 2 Se 3 HAXPES valence band. Furthermore, these findings pave a way for studying the element-resolved band structure and orbital contributions of this class of TIs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗

DOE OSTI2020

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Accelerating defect predictions in semiconductors using graph neural networks

First-principles computations reliably predict the energetics of point defects in semiconductors but are constrained by the expense of using large supercells and advanced levels of theory. Machine learning models trained on computational data, especially ones that sufficiently encode defect coordination environments, can be used to accelerate defect predictions. Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of group IV, III–V, and II–VI zinc blende semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely crystal graph convolutional neural network, materials graph network, and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98% given the range of values within the dataset, improving significantly on the state-of-the-art. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. The current models are based on the semi-local generalized gradient approximation-Perdew–Burke–Ernzerhof (PBE) functional but are highly promising because of the correlation of computed energetics and defect levels with higher levels of theory and experimental data, the accuracy and necessity of discovering novel metastable and low energy defect structures at the PBE level of theory before advanced methods could be applied, and the ability to train multi-fidelity models in the future with new data from non-local functionals. The DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially optimized defective structures, helping identify electronically active defects in technologically important semiconductors.

Rahman, Md Habibur (ORCID:000000027705984X)↗

DOE OSTI2025

Power Flow Geometry and Approximation

Here, the power flow equations are important in numerous power systems problems of practical interest which consider alternating current power flow (ACPF) physics. Perhaps the most well studied being the alternating current optimal power flow problem (ACOPF), seeking to optimize the operation of an electric power system. Due to their non-linearity, problems which include the power flow equations are typically challenging, particularly in optimization. Interestingly, the set of solutions to the power flow equations forms a smooth manifold. As a result, differential geometry can be used to describe and analyze this set of equations. This approach has proven effective in several engineering applications (e.g., solving ACOPF and analyzing the solution space boundary). Central to the success of this approach is an understanding of the power flow manifold's geometry. In this work, we develop the geometric and topological properties of this manifold using concepts from differential geometry. After demonstrating the convenience of this manifold's representation as a function's graph, computational methods are emphasized: we develop retractions, error bounds for linear approximation, and formulas for evaluating the Riemannian metric (including associated objects such as geodesics and the curvature tensor). Scalar curvature and the second fundamental form play a new role in quantifying the quality of linear approximations, like the popular direct current approximation. All functions are implemented in Julia and available in an online repository. Proofs are included for completeness.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

First-Principles Calculations to Investigate the Structural, Electronic, and Half-Metallic Properties of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} (x = 0, 0.25, 0.5, 0.75, and 1) Quaternary Heusler Alloys

The structural, electronic, and magnetic properties of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} (x = 0, 0.25, 0.5, 0.75, and 1) quaternary Heusler alloys have been investigated using first-principles density functional theory within generalized gradient approximation (GGA) and local spin-density approximation with Hubbard-U corrections (LSDA+U). The calculations have been performed using the full-potential linearized augmented Plane wave plus local orbitals method. The crystal structures of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} were found to be cubic in (F-43m) symmetry for x = 0, 1 and cubic in (P-43m) symmetry for x = 0.25, 0.75, while tetragonal in (P-4m2) symmetry for x = 0.5. The results show that all Ti{sub 2}RhSn1-xSix, Ti{sub 2}RhGe{sub 0.75}Si{sub 0.25}, and Ti{sub 2}RhGe alloys are half-metallic ferrimagnetic. These new materials are good candidates for potential applications in spintronic.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Magneto-Electronic, Thermodynamic, and Thermoelectric Properties of 5f-Electron System BaBkO{sub 3}

In the present work, the structural, electronic, magnetic, mechanical, thermodynamic, and thermoelectric properties of cubic perovskite BaBkO{sub 3} were investigated within the framework of density functional theory using generalized gradient approximation, onsite Coulomb interaction, and modified Becke-Johnson methods. BaBkO{sub 3} is stable in the ferromagnetic phase with a magnetic moment of 7 μ{sub B}. The spin-polarized electronic band structure reveals 100% spin polarization of the material with metallic character in spin-up state and semiconducting in spin-down state. Pugh and Cauchy relation displays that the material is highly ductile. The quasi-harmonic Debye model was used to determine the thermodynamic properties of BaBkO{sub 3} under high pressure and temperature. BaBkO{sub 3} presents high value of Seebeck coefficient and power factor at room temperature with a value of 135 μVK{sup − 1} and 1.5 × 10{sup 12}μWcm{sup − 1} K{sup − 2} s{sup − 1}, respectively. The possession of half metallicity and high Seebeck coefficient makes BaBkO{sub 3} suitable for spintronic and thermoelectric device applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

A Finite-Field Approach for GW Calculations beyond the Random Phase Approximation

In this work, we describe a finite-field approach to compute density response functions, which allows for efficient G0W0 and G0W0Γ0 calculations beyond the random phase approximation. The method is easily applicable to density functional calculations performed with hybrid functionals. Lastly, we present results for the electronic properties of molecules and solids and we discuss a general scheme to overcome slow convergence of quasiparticle energies obtained from G0W0Γ0 calculations, as a function of the basis set used to represent the dielectric matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Pre-equilibrium evolution of conserved charges with initial conditions in the ICCING Monte Carlo event generator

Heavy-ion collisions can be well described through relativistic viscous hydrodynamics, but questions still remain when hydrodynamics is applicable because the initial state may begin very far from equilibrium. Thus, a pre-equilibrium evolution phase is used to bridge the gap between the initial state and hydrodynamics. KøMPøST is one such pre-equilibrium model that propagates the energy-momentum tensor by decomposing it into the background and fluctuations around that background, whose evolution is captured by Green's functions. We extend this formalism to include conserved charges and calculate the corresponding nonequilibrium Green's functions in the relaxation-time approximation. The ICCING algorithm initializes conserved charges in the initial state by sampling $g$ → $q$$\overline{q}$ splitting probabilities and is, thus, perfectly positioned to implement Green's functions for charge propagation. As a result, we show that this method alters the initial-state charge geometries and is applicable in central to mid-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

First-Principles Study of Pressure Dependence of Optical Spectra of MnS

MnS is a semiconductor material that can be used in solar cell coatings as window/buffer material and in many other applications. The present contribution deals with the optical properties of MnS in the zinc-blende structure using ab initio calculations within the generalized gradient approximation. Features such as dielectric function, refractive index, reflectivity, optical absorption, conductivity, and electron loss function spectra have been reported and their pressure dependence has been examined and discussed. The information gathered from the present study can be useful for photovoltaic applications.

14 SOLAR ENERGY↗

DOE OSTI2018

Structural, Elastic, Thermodynamic, Electronic, and Magnetic Investigations of Full-Heusler Compound Ag{sub 2}CeAl: FP-LAPW Method

Structural, elastic, thermodynamic, electronic, and magnetic properties of the full-Heusler compound Ag{sub 2}CeAl were determined using generalized gradient approximation with exchange-correlation functional GGA (PBEsol) with spin-orbit coupling (SOC) correction. The elastic modulus and their pressure dependence are calculated. From the elastic parameter behavior, it is inferred that this compound is elastically stable and ductile in nature. Through the quasi-harmonic Debye model, in which the phononic effect is considered the effect of pressure P (0 to 50) and temperature T (0 to 1000) on the lattice constant, the elastic parameters, bulk modulus B, heat capacity and thermal expansion α, internal energy U, entropy S, Debye temperature 𝜃{sub D}, Helmholtz free energy A, and Gibbs free energy G are investigated. The thermodynamic properties show that the compound Ag{sub 2}CeAl is a heavy fermion material. The density of state (DOS), magnetic momentum, and band structure are computed, to investigate the magnetic and metallic characteristics. The calculated polarization of the compound is 77.34%. The obtained results are the first predictions of the physical properties for the rare-earth-based (Ce) full-Heusler Ag{sub 2}CeAl.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗