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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Assessing DFT-D3 Damping Functions Across Widely Used Density Functionals: Can We Do Better?

With the aim of improving the utility of the DFT-D3 empirical dispersion correction, we herein generalize the DFT-D3 damping function by optimizing an additional parameter, an exponent, which controls the rate at which the dispersion tail is activated. This method- DFT-D3(op), shorthand for “optimized power,” where power refers to the newly introduced exponent- is then parametrized for use with ten popular density functional approximations across a small set of noncovalent interactions and isomerization energies; the resulting methods are then evaluated across a large independent test set of 2475 noncovalent binding energies and isomerization energies. We find that the DFT-D3(op) tail represents a substantial improvement over existing damping functions, as it affords significant reductions in errors associated with noncovalent interaction energies and geometries. Finally, the revPBE0-D3(op) and MS2-D3(op) methods in particular stand out, and our extensive testing indicates they are competitive with other modern density functionals.

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

Self-interaction correction in water–ion clusters

We study the importance of self-interaction errors in density functional approximations for various water–ion clusters. We have employed the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in conjunction with LSDA, PBE, and SCAN to describe binding energies of hydrogen-bonded water–ion clusters, i.e., water–hydronium, water– hydroxide, water–halide, as well as non-hydrogen-bonded water–alkali clusters. In the hydrogen-bonded water–ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules, because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water– ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher–level calculations. Here, the non-hydrogen-bonded water–alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies (≤ 50 meV/H 2 O) for hydrogen-bonded water–ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement.

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

G 0 W 0 Ionization Potentials of First-Row Transition Metal Aqua Ions

Here, we report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V–Cu) aqua ions in their most common oxidation states, i.e., V 3+ , Cr 3+ , Cr 2+ , Mn 2+ , Fe 3+ , Fe 2+ , Co 2+ , Ni 2+ , and Cu 2+ . The d-orbital occupancy of these systems spans a broad range from d 2 to d 9 . All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G 0 W 0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G 0 W 0 @PBE, G 0 W 0 @PBE0, and G 0 W 0 @r 2 SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G 0 W 0 @PBE0 and G 0 W 0 @r 2 SCAN approaches, with the r 2 SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.

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

Accuracy of Density Functional Theory Methods for the Calculation of Magnetic Exchange Couplings in Binuclear Iron(III) Complexes

We assess the performance of popular hybrid density functional approximations for the prediction of magnetic exchange couplings (J) in FeIII complexes. To this end, eleven dinuclear oxo-bridged iron(III) complexes with accurately determined J values were chosen to cover cases of weak and strong antiferromagnetic couplings ranging from -6 cm -1 to -132 cm -1 . Six representative exchange-correlation functionals ranging from the hybrids PBE and B3LYP to meta hybrids (M06, M062X), range separated functionals (HSE, LC-!PBE), and LSDA with 30 % of HF exchange were used for this assessment. These functionals have been suggested as good performers for this particular property in transition metal complexes and hence this assessment aims to narrow the interval of con dence for the particular case of Fe complexes. In addition, we examined the performance of an alternative method based on di erential local spin rotations to calculate magnetic exchange coupling parameters in these complexes. Our results show that, although some of these functionals perform on par, the M06 and HSE functionals yields magnetic exchange couplings in closer agreement with experimental data, with a mean absolute percentage errors of 8.6 and 9.2 %, respectively.

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

Spin–Orbit Matrix Elements for a Combined Spin-Flip and IP/EA approach

We introduce a practical approach for computing the Breit–Pauli spin–orbit matrix elements of multiconfigurational systems with both spin and spatial degeneracies based on our recently developed RAS-nSF-IP/EA method (Houck, S. E.; et al. J. Chem. Theory Comput. 2019, 15, 2278). The spin–orbit matrix elements over all the multiplet components are computed using a single one-particle reduced density matrix as a result of the Wigner–Eckart theorem. A mean field spin–orbit approximation was used to account for the two-electron contributions. Basis set dependence as well as the effect of including additional excitations is presented. The effect of correlating the core and semicore orbitals is also examined. Surprisingly accurate results are obtained for spin–orbit coupling constants, despite the fact that the efficient wave function approximations we research neglect the bulk of dynamical correlation.

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

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2018

Thermal expansion in dispersion-bound molecular crystals

In this paper, we explore how anharmonicity, nuclear quantum effects (NQE), many-body dispersion interactions, and Pauli repulsion influence thermal properties of dispersion-bound molecular crystals. Accounting for anharmonicity with ab initio molecular dynamics yields cell parameters accurate to within 2% of experiment for a set of pyridinelike molecular crystals at finite temperatures and pressures. From the experimental thermal expansion curve, we find that pyridine-I has a Debye temperature just above its melting point, indicating sizable NQE across the entire crystalline range of stability. We find that NQE lead to a substantial volume increase in pyridine-I (≈ 40% more than classical thermal expansion at 153 K) and attribute this to intermolecular Pauli repulsion promoted by intramolecular quantum fluctuations. Finally, when predicting delicate properties such as the thermal expansivity, we show that many-body dispersion interactions and more sophisticated density functional approximations improve the accuracy of the theoretical model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Quantum Chemical Calculation of p K a s of Environmentally Relevant Functional Groups: Carboxylic Acids, Amines, and Thiols in Aqueous Solution

Developing accurate quantum chemical approaches for calculating p K a s is of broad interest. Useful accuracy can be obtained by using density functional theory (DFT) in combination with a polarizable continuum solvent model. However, some classes of molecules present problems for this approach, yielding errors greater than 5 p K units. Various methods have been developed to improve the accuracy of the combined strategy. These methods perform well but either do not generalize or introduce additional degrees of freedom, increasing the computational cost. The Solvation Model based on Density (SMD) has emerged as one of the most commonly used continuum solvent models. Nevertheless, for some classes of organic compounds, e.g., thiols, the p K a s calculated with the original SMD model show errors of 6–10 p K units, and we traced these errors to inaccuracies in the solvation free energies of the anions. To improve the accuracy of p K a s calculated with DFT and the SMD model, we developed a scaled solvent-accessible surface approach for constructing the solute–solvent boundary. By using a “direct” approach, in which all quantities are computed in the presence of the continuum solvent, the use of thermodynamic cycles is avoided. Furthermore, no explicit water molecules are required. Three benchmark data sets of experimentally measured p K a values, including 28 carboxylic acids, 10 aliphatic amines, and 45 thiols, were used to assess the optimized SMD model, which we call SMD with a scaled solvent-accessible surface (SMD sSAS ). Of the methods tested, the M06-2X density functional approximation, 6-31+G(d,p) basis set, and SMD sSAS solvent model provided the most accurate p K a s for each set, yielding mean unsigned errors of 0.9, 0.4, and 0.5 p K units, respectively, for carboxylic acids, aliphatic amines, and thiols. As a result, this approach is therefore useful for efficiently calculating the p K a s of environmentally relevant functional groups.

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

Scalar field in Reissner–Nordström spacetime: Bound state and scattering state (with appendix on eliminating oscillation in partial sum approximation of periodic function)

Highlights: • Bound-state and scattering-state solutions of massive scalar fields in R-N spacetime. • Bound-state wave functions and eigenvalues of massive scalar fields in R-N spacetime. • Solving explicit expressions of scattering phase shifts by integral equation methods. • Introducing tortoise coordinates for R-N spacetime. • Eliminating oscillations in the partial sum approximation of periodic functions. In this paper, we solve the massive scalar field in the Reissner–Nordström spacetime. The scalar field in the Reissner–Nordström spacetime has both bound states and scattering states. For bound states, we solve the bound-state wave function and the eigenvalue spectrum. For scattering states, we solve the scattering wave function and give an explicit expression for scattering phase shift by the integral equation method. Especially, we introduce the tortoise coordinate for the Reissner–Nordström spacetime. Moreover, in the calculation of scattering cross sections, we encounter a difficulty in partial sum approximation. If the sum of partial waves cannot be performed exactly, one has to turn to the partial sum approximation which approximates a function by the first several terms of the series. However, in the partial sum approximation there exists an incorrect oscillation which cannot be eliminated by keeping more terms. In the appendix we suggest an approach for eliminating such oscillations in the partial sum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Model selection and parameter estimation in tumor growth models using approximate Bayesian computation-ABC

Cancer is one of the most fatal diseases in the world. Governments and researchers from various areas have continuously concentrated efforts to better understand the disease and propose diagnostic and treatment techniques. The use of mathematical models of tumor growth is of great importance for the development of such techniques. Due to the variety of models nowadays available in the literature, the problems of model selection and parameter estimation come into picture, aiming at suitably predicting the patient’s status of the disease. As the available data on dependent variables of existing models might not justify the use of common likelihood functions, approximate Bayesian computation (ABC) becomes a very attractive tool for model selection and model calibration (parameter estimation) in tumor growth models. In the present study, a Monte Carlo approximate Bayesian computation (ABC) algorithm is applied to select among competing models of tumor growth, with and without chemotherapy treatment. Simulated measurements are used in this work. The results obtained show that the algorithm correctly selects the model and estimates the parameters used to generate the simulated measurements.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Assessing the Limitations of Self-Interaction-Corrected Functionals for Describing the Hydrated Electron

Simulating the hydrated electron using density functional theory is challenging due to the prevalence of self-interaction error in standard functionals. Hybrid functionals like PBEh(40) can reasonably describe the chemistry of an excess electron in water and partially mitigate self-interaction error by incorporating exact Hartree–Fock exchange, but they are computationally expensive making them impractical for large-scale and long-time ab initio molecular dynamics simulations. Explicit self-interaction correction schemes that are applied on an orbital-by-orbital basis offer a potential alternative when the correction is limited to the singly occupied molecular orbital obtained with a generalized gradient approximation functional. Here, we examine whether the Perdew–Zunger self-interaction correction scheme applied to the revPBE functional can provide a computationally efficient and physically sensible alternative to PBEh(40) for the hydrated electron. We find that functionals incorporating a self-interaction correction scheme should be viewed with caution when applied to the hydrated electron and its reactivity. Furthermore, we show that it is critical to consider extensive sampling and diverse chemical environments when validating their performance.

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

An efficient approach for solving Klein–Gordon equation arising in quantum field theory using wavelets

A new approach for the solution of Klein–Gordon equation using Legendre wavelet-based approximation method is presented. The properties of Legendre wavelets are used to reduce the problem to the solution of system of algebraic equations. Usually 2{sup (k−1)}M{sup 2} connection coefficients are required to solve the Klein–Gordon equation by Legendre wavelets method as mentioned in the literature. But, our proposed method based on Legendre wavelets and algebraic polynomials require only 2{sup (k−1)}M connection coefficients instead of 2{sup (k−1)}M{sup 2} . Also the convergence analysis and error estimation for the proposed function approximation through the truncated series have been discussed and approved with the exact solution. Illustrative examples are discussed to demonstrate the validity and applicability of the technique with lesser computational effort.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Direct NeTS sampling of nuclear graphite $S(α, β, T)$ in Serpent

For advanced reactor applications, Neural Thermal Scattering (NeTS) modules were developed to predict the thermal scattering law (TSL or $S(α, β, T)$) of a nuclear graphite neutron moderator. NeTS are multi-layer, feedforward artificial neural networks, which act as universal function approximators designed for TSL datasets. In this case, a 4-layer neural network with 164 neurons per layer is trained using FLASSH evaluated data in PyTorch and serialized as a torchscript dictionary to predict $S(α, β, T)$ on-the-fly. Relative, absolute and maximum percent deviations of NeTS from File 7 data generated using the FLASSH code are on the order of 0.01%, 0.1% and 1%, respectively, with low inference latencies of 0.000172 s per $S(α, β, T)$ at a given temperature. Capturing the full dimensionality of possible inelastic neutron-lattice interactions, NeTS functionality is embedded in the Serpent Monte Carlo code, where $S(α, β, T)_{NeTS}$ sampling is conducted on-the-fly and compared to ACE look-up-tables for predicting TREAT criticality. k-eff differences between sampling algorithms of 6 pcm are observed and are within the order of Monte Carlo uncertainty. Compared to discrete and continuous-energy ACE files (30 MB and 131 MB per temperature), the NeTS format is on the order of 200–300 kB for a continuous-temperature, interpolation-free representation of $S(α, β, T)$ and cross sections. NeTS-in-Serpent runtimes comparable with ACE look-up tables are achieved by scaling NeTS for high performance computing architectures with hybrid OpenMP + MPI parallelization. This work validates a novel, self-contained reactor physics framework for predictive cross sections, and demonstrates a general methodology for embedding modern machine learning libraries within existing neutronic analysis frameworks.

Nuclear Criticality Safety Program (NCSP)↗
DOE OSTI2021

Deep Kronecker neural networks: A general framework for neural networks with adaptive activation functions

Here we propose a new type of neural networks, Kronecker neural networks (KNNs), that form a general framework for neural networks with adaptive activation functions. KNNs employ the Kronecker product, which provides an efficient way of constructing a very wide network while keeping the number of parameters low. Our theoretical analysis reveals that under suitable conditions, KNNs induce a faster decay of the loss than that by the feed-forward networks. This is also empirically verified through a set of computational examples. Furthermore, under certain technical assumptions, we establish global convergence of gradient descent for KNNs. As a specific case, we propose the Rowdy activation function that is designed to get rid of any saturation region by injecting sinusoidal fluctuations, which include trainable parameters. The proposed Rowdy activation function can be employed in any neural network architecture like feed-forward neural networks, Recurrent neural networks, Convolutional neural networks etc. The effectiveness of KNNs with Rowdy activation is demonstrated through various computational experiments including function approximation using feed-forward neural networks, solution inference of partial differential equations using the physics-informed neural networks, and standard deep learning benchmark problems using convolutional and fully-connected neural networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Some Problems in Density Functional Theory

Though calculations based on density functional theory (DFT) are used remarkably widely in chemistry, physics, materials science, and biomolecular research and though the modern form of DFT has been studied for almost 60 years, some mathematical problems remain. From a physical science perspective, it is far from clear whether those problems are of major import. For context, we provide an outline of the basic structure of DFT as it is presented and used conventionally in physical sciences, note some unresolved mathematical difficulties with those conventional demonstrations, then pose several questions regarding both the time-independent and time-dependent forms of DFT that could benefit from attention in applied mathematics. Furthermore, progress on any of these would aid in development of better approximate functionals and in interpretation of DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Defect-induced states, defect-induced phase transition, and excitonic states in bent tungsten disulfide (WS 2 ) nanoribbons: Density functional vs. many body theory

Two-dimensional (2D) transition metal dichalcogenide (TMD) materials have versatile electronic and optical properties. TMD nanoribbons exhibit interesting properties due to reduced dimensionality, quantum confinement, and edge states, which make them suitable for various electronic and optoelectronic applications. In a previous work conducted by our group, we demonstrated that the edge bands evolved with bending can tune the optical properties for various widths of TMD nanoribbons. Defects are commonly present in 2D TMD materials, and can dramatically change the material properties. Here, in this following work, we investigate the interaction between the edge and the defect states in tungsten disulfide (WS 2 ) nanoribbons with lines of W- or S-atom vacancies defects under different bending conditions, using density functional theory (DFT). To gain understanding about the limits of density functional approximations, we compare results on band gaps and energies of defect states with quasiparticle GW. We reveal interesting semiconductor-metal phase transitions, suggesting potential applications in nanoelectronics or molecular electronics. We also calculate the optical absorption of the bent and defective nanoribbons with the many-body GW-BSE (Bethe-Salpeter equation) approach, revealing a tunable optical spectrum and diverse exciton states in the defective WS 2 nanoribbons.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Fermi-Löwdin orbital self-interaction correction using the strongly constrained and appropriately normed meta-GGA functional

Despite the success of density functional approximations (DFAs) in describing the electronic properties of many-electron systems, the most widely used approximations suffer from self-interaction errors (SIE) that limit their predictive power. Herein we describe the effects of removing SIE from the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA) using the Fermi-Lowdin Orbital Self-Interaction Correction (FLOSIC) method. FLOSIC is a size-extensive implementation of the Perdew-Zunger self-interaction correction (PZ-SIC) formalism. We find that FLOSIC-SCAN calculations require careful treatment of numerical details and describe an integration grid that yields reliable accuracy with this approach. We investigate the performance of FLOSIC-SCAN for predicting a wide array of properties and find that it provides better results than FLOSIC-LDA and FLOSIC-PBE in nearly all cases. Furthermore, it gives better predictions than SCAN for orbital energies and dissociation energies where self-interaction effects are known to be important, but total energies and atomization energies are made worse. For these properties, we also investigate the use of the self-consistent FLOSIC-SCAN density in the SCAN functional and find that this DFA@FLOSIC-DFA approach yields improved results compared to pure, self-consistent SCAN calculations. Thus FLOSIC-SCAN provides improved results over the parent SCAN functional in cases where SIEs are dominant, and even when they are not, if the SCAN@FLOSIC-SCAN method is used.

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

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗