Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional approximations”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

360⁰ -View of Quantum Theory and Ab Initio Simulation at Extreme Conditions: 2014 Sanibel Symposium

The Sanibel Symposium 2014 was held February 16-21, 2014, at the King and Prince, St. Simons Island, GA. It was successful in bringing condensed-matter physicists and quantum chemists together productively to drive the emergence of those specialties. The Symposium had a significant role in preparing a whole generation of quantum theorists. The 54th Sanibel meeting looked to the future in two ways. We had 360⁰-View sessions to honor the exceptional contributions of Rodney Bartlett (70), Bill Butler (70), Yngve Öhrn (80), Fritz Schaefer (70), and Malcolm Stocks (70). The work of these five has greatly impacted several generations of quantum chemists and condensed matter physicists. The “360⁰” is the sum of their ages. More significantly, it symbolizes a panoramic view of critical developments and accomplishments in theoretical and computational chemistry and physics oriented toward the future. Thus, two of the eight 360⁰-View sessions focused specifically on younger scientists. The 360⁰-View program was the major component of the 2014 Sanibel meeting. Another four sessions included a sub-symposium on ab initio Simulations at Extreme Conditions, with focus on getting past the barriers of present-day Born-Oppenheimer molecular dynamics by advances in finite-temperature density functional theory, orbital-free DFT, and new all-numerical approaches.

360-view↗
DOE OSTI2022

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

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

Fast Emulation of Expensive Simulations using Approximate Gaussian Processes [Slides]

Nuclear Computational Low-Energy Initiative (NUCLEI) collaboration uses Density Functional Theory (DFT) simulations to predict the structure and binding energies of nuclei over a wide range of proton (Z) and neutron (N) numbers. The DFT simulations utilize a particular parameterization of a Skyrme energy density functional called UNEDF1 which depends on 12 free parameters that must be fit to data (M Kortelainen et al 2014). Fitting involves comparing (e.g.) predicted binding energies of nuclei to experimentally measured values. We use only binding energies as observables, but DFT with UNEDF1 will predict structure (shape) observables as well. In this work, assessing the capability of approximate GP emulators to balance emulator accuracy with computational speed to facilitate improved UNEDF1 calibration. Sparse GPs are straightforward to train and accurate. Calibration is not straightforward with MCMC (using MH or HMC/NUTS). We produced reusable software for continuing and building on this work as well as accessing and using Darwin cluster compute resources

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

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

Experimental and first principle studies on electronic structure of BaTiO{sub 3}

We have carried out photoemission experiments to obtain valence band spectra of various crystallographic symmetries of BaTiO{sub 3} system which arise as a function of temperature. We also present results of a detailed first principle study of these symmetries of BaTiO{sub 3} using generalized gradient approximation for the exchange-correlation potential. Here we present theoretical results of density of states obtained from DFT based simulations to compare with the experimental valence band spectra. Further, we also perform calculations using post density functional approaches like GGA + U method as well as non-local hybrid exchange-correlation potentials like PBE0, B3LYP, HSE in order to understand the extent of effect of correlation on band gaps of different available crystallographic symmetries (5 in number) of BaTiO{sub 3}.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Improved statistical models for limited datasets in uncertainty quantification using stochastic collocation

This paper presents a data-driven framework for performing uncertainty quantification (UQ) by choosing a stochastic model that accurately describes the sources of uncertainty in a system. This model is propagated through an appropriate response surface function that approximates the behavior of this system using stochastic collocation. Given a sample of data describing the uncertainty in the inputs, our goal is to estimate a probability density function (PDF) using the kernel moment matching (KMM) method so that this PDF can be used to accurately reproduce statistics like mean and variance of the response surface function. Instead of constraining the PDF to be optimal for a particular response function, we show that we can use the properties of stochastic collocation to make the estimated PDF optimal for a wide variety of response functions. We contrast this method with other traditional procedures that rely on the Maximum Likelihood approach, like kernel density estimation (KDE) and its adaptive modification (AKDE). We argue that this modified KMM method tries to preserve what is known from the given data and is the better approach when the available data is limited in quantity. We test the performance of these methods for both univariate and multivariate density estimation by sampling random datasets from known PDFs and then measuring the accuracy of the estimated PDFs, using the known PDF as a reference. Comparing the output mean and variance estimated with the empirical moments using the raw data sample as well as the actual moments using the known PDF, we show that the KMM method performs better than KDE and AKDE in predicting these moments with greater accuracy. This improvement in accuracy is also demonstrated for the case of UQ in electrostatic and electrothermomechanical microactuators. We show how our framework results in the accurate computation of statistics in micromechanical systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Relating polarizability to volume, ionization energy, electronegativity, hardness, moments of momentum, and other molecular properties

Semiquantitative relationships between the mean static dipole polarizability and other molecular properties such as the volume, ionization energy, electronegativity, hardness, and moments of momentum are explored. The relationships are tested using density functional theory computations on the 1641 neutral, ground-state, organic molecules in the TABS database. The best polarizability approximations have median errors under 5%.

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

Neutron diffusion in a randomly inhomogeneous multiplying medium with random phase approximation

Neutron diffusion in a randomly inhomogeneous multiplying medium is studied. By making use of a random phase assumption we show that the average neutron density approximately satisfies an integral equation in Fourier space, which is solved using Kummer functions. We used multi-dimensional formulation. In the case of one dimension, we obtain the result of Rosenbluth and Tao for the mean total density for large t. In the three-dimensional case, a closed form of solution is derived for the mean total neutron density. Its asymptotic behavior is also investigated for large t.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Hybrid functional band gap calculation of SnO{sub 6} containing perovskites and their derived structures

We have studied the properties of SnO{sub 6} octahedra-containing perovskites and their derived structures using ab initio calculations with different density functionals. In order to predict the correct band gap of the materials, we have used B3LYP hybrid density functional, and the results of B3LYP were compared with those obtained using the local density approximation and generalized gradient approximation data. The calculations have been conducted for the orthorhombic ground state of the SnO{sub 6} containing perovskites. We also have expended the hybrid density functional calculation to the ASnO{sub 3}/A'SnO{sub 3} system with different cation orderings. We propose an empirical relationship between the tolerance factor and the band gap of SnO{sub 6} containing oxide materials based on first principles calculation. - Graphical abstract: (a) Structure of ASnO{sub 3} for orthorhombic ground state. The green ball is A (Ba, Sr, Ca) cation and the small (red) ball on edge is oxygen. SnO{sub 6} octahedrons are plotted as polyhedron. (b) Band gap of ASnO{sub 3} as a function of the tolerance factor for different density functionals. The experimental values of the band gap are marked as green pentagons. (c) ASnO{sub 3}/A'SnO{sub 3} superlattices with two types cation arrangement: [001] layered structure and [111] rocksalt structure, respectively. (d) B3LYP hybrid functional band gaps of ASnO{sub 3}, [001] ordered superlattices, and [111] ordered superlattices of ASnO{sub 3}/A'SnO{sub 3} as a function of the effective tolerance factor. Note the empirical linear relationship between the band gap and effective tolerance factor. - Highlights: • We report the hybrid functional band gap calculation of ASnO{sub 3} and ASnO{sub 3}/A'SnO{sub 3}. • The band gap of ASnO{sub 3} using B3LYP functional reproduces the experimental value. • We propose the linear relationship between the tolerance factor and the band gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Ground-state and dynamical properties of two-dimensional dipolar Fermi liquids

We study the ground-state properties of a two-dimensional spin-polarized fluid of dipolar fermions within the Euler–Lagrange Fermi-hypernetted-chain approximation. Our method is based on the solution of a scattering Schrödinger equation for the “pair amplitude” √(g(r)), where g(r) is the pair distribution function. A key ingredient in our theory is the effective pair potential, which includes a bosonic term from Jastrow–Feenberg correlations and a fermionic contribution from kinetic energy and exchange, which is tailored to reproduce the Hartree–Fock limit at weak coupling. Very good agreement with recent results based on quantum Monte Carlo simulations is achieved over a wide range of coupling constants up to the liquid-to-crystal quantum phase transition. Using the fluctuation–dissipation theorem and a static approximation for the effective inter-particle interactions, we calculate the dynamical density–density response function, and furthermore demonstrate that an undamped zero-sound mode exists for any value of the interaction strength, down to infinitesimally weak couplings. -- Highlights: •We have studied the ground state properties of a strongly correlated two-dimensional fluid of dipolar fermions. •We have calculated the effective inter-particle interaction and the dynamical density–density response function. •We have shown that an undamped zero sound mode exists at any value of the interaction strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Computational modeling of mixed-ligand cobalt complexes with o-quinone derivative of corannulene

The structure and energy characteristics of isomers of electrically neutral adducts of tetracoordinate cobalt complexes with o-quinone derivative of corannulene were studied quantum chemically within the density functional theory (DFT) in the UB3LYP*/6-311++G(d,p) approximation. It was shown that the framework rigidity of the redox-active ligand favors convergence of the energy levels of the high-spin and low-spin isomers of the complexes under study as compared to the o-benzoquinone-containing analogues. At the same time, steric strain in the corannulene moiety causes the stability of the adducts to reduce. A compound characterized by high stabilization energy and the ability to switch its spin states via valence tautomeric rearrangements was found by systematic variation of structural parameters of cobalt bis(chelate).

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

Electronic Properties of Ultra‐Wide Bandgap B x Al 1− x N Computed from First‐Principles Simulations

Abstract Ultra‐wide bandgap (UWBG) materials such as AlN and BN hold great promise for future power electronics due to their exceptional properties. They exhibit large bandgaps, high breakdown fields, high thermal conductivity, and high mechanical strengths. AlN and BN have been extensively researched, however, their alloys, B x Al 1− x N, are much less studied despite their ability to offer tunable properties by adjusting x . In this article, the electronic properties of 17 recently predicted ground states of B x Al 1− x N in the x = 0 − 1 range are predicted using first‐principles density functional theory and many‐body perturbation theory within GW approximation. All the B x Al 1− x N structures are found to be UWBG materials and have bandgaps that vary linearly from that of wurtzite‐phase ( w ) AlN (6.19 eV) to that of w ‐BN (7.47 eV). The bandstructures of B x Al 1− x N show that a direct‐to‐indirect bandgap crossover occurs near x = 0.25. Furthermore, it is found that B x Al 1− x N alloys have much larger dielectric constants than the constituent bulk materials (AlN = 9.3 ɛ 0 or BN = 7.3 ɛ 0 ), with values reaching as high as 12.1 ɛ 0 . These alloys are found to exhibit large dielectric breakdown fields in the range 9–35 MV cm −1 with a linear dependence on x . This work provides the much needed advancement in the understanding of the properties of B x Al 1− x N to aid their application in next‐generation devices.

Milne, Cody L.↗
DOE OSTI2016

Why is MP2-Water "Cooler" and "Denser" than DFT-Water?

To maintain water in the liquid phase at the correct (1 g/cm3) density during first-principles simulations, density-functional theory (DFT) with a dispersionless generalized-gradient-approximation (GGA) functional requires a much higher temperature and pressure than the ambient conditions. Conversely, ab initio second-order many-body perturbation (MP2) calculations of liquid water performed by Del Ben et al. [J. Chem. Phys. Lett. 4, 3753 (2013); J. Chem. Phys. 143, 054506 (2015)] and by us [Willow et al., Sci. Rep. 5, 14358 (2015)] required a lower temperature and a negative pressure than DFT to keep water liquid. Here, we present a unifying explanation of these trends derived from classical water simulations using a polarizable force field with different sets of parameters. We show that the calculated temperature and pressure of the liquid phase are strongly correlated with the polarizability of water and the dispersion interaction, respectively. In DFT/GGA, the polarizability and thus the induced dipole moments and the hydrogen-bond strength are all overestimated. This hinders the rotational motion of molecules and requires a higher temperature for water to be liquid. In MP2 and DFT/GGA, the dispersion interaction is stronger and weaker (or lacking), respectively. This explains why liquid water contracts uniformly and becomes too dense in MP2, whereas the opposite is the case for dispersionless DFT/GGA.

Willow, Soohaeng Y.↗
DOE OSTI2016

Thermodynamic Stability of Low- and High-Index Spinel LiMn 2 O 4 Surface Terminations

Density functional theory calculations are performed within the generalized gradient approximation (GGA+U) to determine stable terminations of both low- and high-index spinel LiMn2O4 (LMO) surfaces. A grand canonical thermodynamic approach is employed, permitting a direct comparison of offstoichiometric surfaces with previously reported stoichiometric surface terminations at various environmental conditions. Within this formalism, we have identified trends in the structure of the low-index surfaces as a function of the Li and O chemical potentials. The results suggest that, under a range of chemical potentials for which bulk LMO is stable, Li/O and Li-rich (111) surface terminations are favored, neither of which adopts an inverse spinel structure in the subsurface region. This thermodynamic analysis is extended to identify stable structures for certain high-index surfaces, including (311), (331), (511), and (531), which constitute simple models for steps or defects that may be present on real LMO particles. The low- and high-index results are combined to determine the relative stability of each surface facet under a range of environmental conditions. The relative surface energies are further employed to predict LMO particle shapes through a Wulff construction approach, which suggests that LMO particles will adopt either an octahedron or a truncated octahedron shape at conditions in which LMO is thermodynamically stable. These results are in agreement with the experimental observations of LMO particle shapes.

DFT↗
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 OSTI2020

Interpreting Tafel behavior of consecutive electrochemical reactions through combined thermodynamic and steady state microkinetic approaches

Assessing the reaction pathway of multi-electron-transfer reactions is an essential yet difficult task for the rational design of electrocatalysts. In this work, we develop a heuristic approach that combines thermodynamic adsorption energetics calculated through density functional theory with microkinetic modeling using the steady state approximation to interpret the potential-dependent Tafel behavior of consecutive electrochemical reactions. In doing so, we introduce a kinetic framework for ab initio calculations that ensures self-consistent adsorption energetics based on kinetically limited adsorbate coverages. The approach is applied to experimental results on CoO x (OH) 2-x single crystal electrocatalyst particles yielding coverage dependent mechanistic information and identification of the rate-limiting step with standard rate constants for the oxygen evolution reaction on the ($11\bar{2}0$) surfaces of the β-Co(OH) 2 , β-CoOOH, and CoO 2 bulk phases. Finally, this generalizable method enables catalyst benchmarking based on determining the active species involved and associated intrinsic reaction rate constants in consecutive multi-electron-transfer reactions.

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

Advances in molecular quantum chemistry contained in the Q-Chem 4 program package

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Finally, some specific chemical examples include calculations on the strongly correlated Cr 2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.

computational modelling↗