Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy density functional”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗
DOE OSTI2020

The role of ligands in atomically precise nanocluster-catalyzed CO2 electrochemical reduction

Ligand effects are of major interest in catalytic reactions owing to their potential critical role in determining the reaction activity and selectivity. Herein, we report ligand effects in the CO 2 electrochemical reduction reaction at the atomic level with three unique Au 25 nanoclusters comprising the same kernel but different protecting ligands (–XR, where X = S or Se, and R represents the carbon tail). It is observed that a change in the carbon tail shows no obvious impact on the catalytic selectivity and activity, but the anchoring atom (X = S or Se) strongly affects the electrocatalytic selectivity. Specifically, the S site acts as the active site and sustains CO selectivity, while the Se site shows a higher tendency of hydrogen evolution. Density functional theory (DFT) calculations reveal that the energy penalty associated with the *COOH formation is lower on the S site by 0.26 eV compared to that on the Se site. Additionally, the formation energy of the product (*CO) is lower on the sulfur-based Au nanocluster by 0.43 eV. Furthermore, we attribute these energetic differences to the higher electron density on the sulfur sites of the Au nanocluster, resulting in a modified bonding character of the reaction intermediates that reduce the energetic penalty for the *COOH and *CO formation. Overall, this work demonstrates that S/Se atoms at the metal–ligand interface can play an important role in determining the overall electrocatalytic performance of Au nanoclusters.

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

Interplay between Facets and Defects during the Dissociative and Molecular Adsorption of Water on Metal Oxide Surfaces

Both surface terminations and defects play a central role in determining how water interacts with metal oxides, thereby setting important properties of the interface that govern reactivity such as the type and distribution of hydroxyl groups. However, the interconnections between facets and defects remain poorly understood, limiting the usefulness of conventional notions such as that hydroxylation is controlled by metal cation exposure at the surface. Here, using hematite (α-Fe 2 O 3 ) as a model system, we show how oxygen vacancies overwhelm surface cation-dependent hydroxylation behavior. Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy was used to monitor the adsorption of molecular water and its dissociation to form hydroxyl groups in situ on (001), (012), or (104) facet-engineered hematite nanoparticles. Supported by density functional theory calculations of the respective surface energies and oxygen vacancy formation energies, the findings show how oxygen vacancies are more prone to form on higher energy facets and induce surface hydroxylation at extremely low relative humidity values of 5 x10 - 5 %. Further, when these vacancies are eliminated, the extent of surface hydroxylation across the facets is as expected from the areal density of exposed iron cations at the surface. These findings help answer fundamental questions about the nature of reducible metal oxide-water interfaces in natural and technological settings and lay the groundwork for rational design of improved oxide-based catalysts.

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

Role of ripples in altering the electronic and chemical properties of graphene

Ripples of graphene are known to manipulate electronic and hydrogenation properties of graphitic materials. More detailed work is needed to elucidate the structure–property relationship of these systems. In this work, the density functional theory is used to compute the energy and electronic structure of the graphene models with respect to variable curvatures and hydrogen adsorption sites. The magnitude of finite bandgap opening depends on the orientation of ripples, and the hydrogen adsorption energy depends on the local curvature of graphene. An adsorbed hydrogen alters the local curvature, resulting in relatively weakened adsorption on the neighboring three sites, which gives a rationale to experimentally observed dynamic equilibrium stoichiometry (H:C = 1:4) of hydrogenated graphene. The surface diffusion transition state energy of adsorbed hydrogen is computed, which indicates that the Eley–Rideal surface recombination mechanism may be important to establish the dynamic equilibrium, instead of the commonly assumed Langmuir–Hinshelwood mechanism.

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

Modulating the Electron Affinity of Small Bipyridyl Molecules on Single Gold Nanoparticles for Plasmon-Driven Electron Transfer

Developing controlled platforms for plasmon-driven chemistry is of great importance in catalytic reactions at the nanoscale. We report anion radical formation for five bipyridyl complexes of varying degrees of electron affinity utilizing optically fo-cused intraband (594 nm) and interband (532 nm) pump excitation of single gold nanoparticles. The surface-enhanced Raman scattering (SERS) of anion radicals for the five non-resonant adsorbed molecules 2,2’-bipyridine (22BPY), 4,4’-bipyridine (44BPY), trans-1,2-bis(4-pyridyl)ethylene (BPE), 1,2-bis(4-pyridyl)acetylene (BPA), and 1,2-bis(4-pyridyl)ethane (BPEt) were detected using localized surface-plasmon resonance (LSPR) excitation with 785 nm. The electron affinity of the five bipyridyl complexes were determined using electrochemistry. Molecules with low electron affinity experienced high-er instances of radical anion formation under a plasmon-coupled intraband electron transfer excitation (594 nm) whereas molecules with high electron affinity showed a preference for anion radical formation under direct interband electron trans-fer excitation (532 nm). The lowest unoccupied molecular orbital (LUMO) energy levels for low electron affinity surface-bound molecules (22BPY, BPEt) are on average ~0.43 eV higher for high electron affinity surface-bound molecules (BPA, BPE, 44BPY) as calculated using time-dependent density functional theory, elucidating the importance of plasmon coupling to energy levels that facilitate charge transfer pathways. We also show the ability to ‘activate’ high vs low electron affinity single nanoparticles with the choice of pump excitation wavelength. The findings show the complex interplay between molecular electron affinity, orbital overlap with the density of states of the plasmonic metal, and excitation energetics of the pump laser wavelength. Potential applications of this work include enhanced control over molecular scale catalysis, biosensor design, and solar energy capture.

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

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗
DOE OSTI2023

First-Principles-Based Machine Learning Models for Phase Behavior and Transport Properties of CO 2

In this work, we construct distinct first-principles-based machine-learning models of CO 2 , reproducing the potential energy surface of the PBE-D3, BLYP-D3, SCAN and SCAN-rvv10 approximations of density functional theory. Here we employ the Deep Potential methodology to develop the models and consequently achieve a significant computational efficiency over ab initio molecular dynamics (AIMD) that allows for larger system sizes and time scales to be explored. Although our models are trained only with liquid phase configurations, they are able to simulate a stable interfacial system and predict vapor-liquid equilibrium properties, in good agreement with results from the literature. Because of the computational efficiency of the models, we are also able to obtain transport properties, such as viscosity and diffusion coefficients. We find that the SCAN-based model presents a temperature shift in the position of the critical point, while the SCAN-rvv10-based model shows improvement but still exhibits a temperature shift that remains approximately constant for all properties investigated in this work. We find that the BLYP-D3-based model generally performs better for liquid phase and vapor-liquid equilibrium properties, but the PBE-D3-based model is better suited for predicting transport properties.

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

Nonlocal correlations in iron pnictides and chalcogenides

Deviations of low-energy electronic structurse of iron-based superconductors from density-functional-theory predictions have been parametrized in terms of band- and orbital-dependent mass renormalizations and energy shifts. The former have typically been described in terms of a local self-energy within the framework of dynamical mean field theory, while the latter appears to require nonlocal effects due to interband scattering. By calculating the renormalized band structure in both random phase approximation (RPA) and the two-particle self-consistent approximation (TPSC), we show that correlations in pnictide systems like LaFeAsO and LiFeAs can be described rather well by a nonlocal self-energy. In particular, Fermi pocket shrinkage as seen in experiments occurs due to repulsive interband finite-energy scattering. For the canonical iron chalcogenide system FeSe in its bulk tetragonal phase, the situation is, however, more complex since even including momentum-dependent band renormalizations cannot explain experimental findings. We propose that the nearest-neighbor Coulomb interaction may play an important role in band-structure renormalization in FeSe. Finally, we further compare our evaluations of nonlocal quasiparticle scattering lifetime within RPA and TPSC with experimental data for LiFeAs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational descriptor for electrochemical currents of carbon dioxide reduction on Cu facets

Computation screening is crucial for designing efficient electrochemical catalysts for carbon dioxide (CO 2 R) reduction that produce valuable hydrocarbons and oxygenates. In this work, leveraging density functional theory calculations for the CO adsorption energy ΔE CO on seventeen Cu terminations, we discover a strong linear correlation between ΔE CO and the experimentally measured CO 2 R electrochemical currents (ACS Catal. 2022, 12, 11, 6578–6588). Examining ab initio thermodynamics of early critical intermediates CO*, COH*, and CHO*, we find that CO* → CHO* is the thermodynamically controlling step. Beyond the general CO adsorption energy that only shows a linear trend with CO 2 R activity, we show that the reaction free energy of CO* → CHO* is the descriptor for the overall CO 2 R activity for Cu facets, as it displays a volcano relationship with the experimental current. Importantly, we show that high step and kink density of the Cu terminations not only enhances CO adsorption strength but also modulates the CO* → CHO* pathway, as respectively exemplified in the (941) and (741) facets. In addition, we explain that the high activity of (741) is due to its relatively low hydrogen evolution reaction activity compared with the other Cu surfaces.

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

Factors influencing quantum evaporation of helium from polar semiconductors from first principles

While there is much indirect evidence for the existence of dark matter (DM), to date it has evaded detection. Current efforts focus on DM masses over ~GeV—to push the sensitivity of DM searches to lower masses, new DM targets and detection schemes are needed. In this work, we focus on the latter—a novel detection scheme recently proposed to detect 10–100 meV phonons in polar target materials. Previous work showed that well-motivated models of DM can interact with polar semiconductors to produce an athermal population of phonons. This new sensing scheme proposes that these phonons then facilitate quantum evaporation of 3 He from a van der Waals film deposited on the target material. However, a fundamental understanding of the underlying process is still unclear, with several uncertainties related to the precise rate of evaporation and how it can be controlled. In this work, we use ab initio density functional theory calculations to compare the adsorption energies of helium atoms on a polar target material, sodium iodide, to understand the underlying evaporation physics. We explore the role of surface termination, monolayer coverage, and elemental species on the rate of He evaporation from the target material. Using this, we discuss the optimal target features for He-evaporation experiments and their range of tunability through chemical and physical modifications such as applied field and surface termination.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗
NASA NTRS2018

Li(+)-ligand Binding Energies and the Effect of Ligand Fluorination on the Binding Energies

The Li(+)-ligand binding energies are computed for seven ligands and their perfluoro analogs using Density Functional Theory. The bonding is mostly electrostatic in origin. Thus the size of the binding energy tends to correlate with the ligand dipole moment, however, the charge-induced dipole contribution can be sufficiently large to affect the dipole-binding energy correlation. The perfluoro species are significantly less strongly bound than their parents, because the electron withdrawing power of the fluorine reduces the ligand dipole moment.

Charles W Bauschlicher↗
DOE OSTI2022

Feasibility of Prelithiation in LiFePO 4

Lithium iron phosphate (LiFePO 4 ) is widely applied as the cathode material for the energy storage Li-ion batteries due to its low cost and high cycling stability. However, the low theoretical specific capacity of LiFePO 4 makes its initial capacity loss more concerning. Therefore, lithium compensation by way of prelithiation and applications of sacrificial Li-rich additives in LiFePO 4 is imminent in elevating the energy density and/or prolonging the lifetime of the LiFePO 4 -based Li-ion batteries (LIBs). Prelithiation in LiFePO 4 is herein carried out by electrochemical and chemical methods and its feasibility is proved on the basis of the electrochemical evaluations such as the initial charge capacity and the cycling stability. In addition, the site of the pre-intercalated Li-ions is found via comprehensive physical characterizations and the density functional theory (DFT) calculations. Furthermore, these findings open a new avenue for elevating the energy density and/or prolonging the lifetime of the high-energy-density batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Critical Assessment of the Thermodynamics of Vacancy Formation in Fe 2 O 3 Using Hybrid Density Functional Theory

Fe 2 O 3 hematite is a technologically important material with applications in energy storage as well as being the key phase formed in the rust of iron-based materials. Despite this central importance, there is much that is still unknown regarding the properties of defects in Fe 2 O 3 and consequently oxide growth. In this work, using screened hybrid density functional theory (HSE06), we consider the thermodynamics of vacancies in Fe 2 O 3 , considering the effects of ionic and electronic chemical potentials on both iron and oxygen vacancy formation. We find that, in the oxygen-rich limit, iron vacancies are easier to form, though the difference in formation energy between the two vacancies is only about 1 eV. In contrast, in the Fe-rich limit, oxygen vacancies have an extremely low formation energy, only 0.07 eV at mid-gap, and would spontaneously form as the Fermi level is reduced, while Fe vacancies require over 5 eV to form. Consistent with experiment, this indicates that Fe 2 O 3 is relatively easily reduced but not oxidized. However, the theoretical picture is very different when considering other exchange–correlation functionals (GGA + U or SCAN), emphasizing the critical role of the exchange–correlation functional in describing this system accurately.

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

Use of FLOSIC for understanding anion-solvent interactions

An Achille’s heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3−, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth’s atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi–Löwdin–Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree–Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.

Chemistry↗
DOE OSTI2020

Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

Long range-corrected (LRC) or range-separated hybrid (RSH) functionals where the long-range (LR) limit of electronic interactions is set to the exact exchange have been shown to correct the tendency of traditional density functional theory (DFT) to underestimate the frontier orbital gap. Consequently, the use of such functionals in calculating electronic excited states using linear response based time-dependent DFT (TDDFT) has been successful in correcting the tendency for underestimating the energies of charge transfer states by DFT-based calculations. More recently formulations of functionals that attenuate the LR limit to address condensed-phase effects to polarize the electronic density have been reported. In particular screened RSH (SRSH) combined with polarizable continuum model (PCM) was benchmarked successfully in reproducing the fundamental gap and charge transfer state energies of molecular systems in the condensed phase. Here we use SRSH-PCM to address triplet excited states, and show its success in obtaining correspondence of the low-lying triplet states to the singlet–triplet gap in a similar way that the fundamental orbital gap corresponds to electron removal and addition energies. Importantly, the accuracy of the SRSH-PCM in calculating triplet excitations stands on the polarization consistent framework in addressing the scalar dielectric constant and without affecting the optimal tuning by triplet energies. Finally, the prospect of even further improving the SRSH-PCM accuracy in calculating triplet states can be achieved by optimal tuning on the basis of the spin multiplicity gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗