Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2020

Engineering the electronic and strained interface for high activity of PdM core @Pt monolayer electrocatalysts for oxygen reduction reaction

Alloyed nanoparticles with core-shell structures provide a favorable model to modulate interfacial interaction and surface structures at the atomic level, which is important for designing electrocatalysts with high activity and durability. Herein, core-shell structured Pd 3 M@Pt/C nanoparticles with binary PdM alloy cores (M = Fe, Ni, and Co) and a monolayer Pt shell were successfully synthesized with diverse interfaces. Among these, Pd 3 Fe@Pt/C exhibited the best oxygen reduction reaction catalytic performance, roughly 5.4 times more than that of the commercial Pt/C catalyst used as reference. The significantly enhanced activity is attributed to the combined effects of strain engineering, interfacial electron transfer, and improved Pt utilization. Density functional theory simulations and extended X-ray absorption fine structure analysis revealed that engineering the alloy core with moderate lattice mismatch and alloy composition (Pd 3 Fe) optimizes the surface oxygen adsorption energy, thereby rendering excellent electrocatalytic activity. Finally, future researches may use this study as a guide on the construction of highly effective core-shell electrocatalysts for various energy conversions and other applications.

25 ENERGY STORAGE↗
DOE OSTI2026

Identification of In Situ Generated High-Spin Fe(II)-N 4 Active Sites for Acidic Oxygen Reduction Reaction via Operando 57 Fe Mössbauer Spectroscopy

Iron-nitrogen-carbon (FeNC) catalysts are considered among the most promising alternative to Pt catalysts in acidic oxygen reduction reaction (ORR), yet the geometric and electronic structures of the true active site under reaction conditions have not been clearly elucidated. Herein, we synthesized a representative FeNC catalyst by pyrolyzing Fe 3+ absorbed on ZIF-8-derived N-doped carbon at a mild temperature under the H 2 /Ar atmosphere, where a formation mechanism of FeN 4 sites through a Zn-mediated Fe nanoparticle atomization process was proposed. The resulting FeNC-750 catalyst shows high acidic ORR activity with a half-wave potential of 0.838 V and a peak power density of 1 W cm -2 in proton-exchange membrane fuel cell (PEMFC). By using operando 57 Fe Mössbauer spectroscopy on FeNC-750, it was revealed that the pyrrolic N-coordinated high-spin Fe 2+ N 4 sites, which are in situ generated from high-spin Fe 3+ N 4 during ORR, are identified as catalytic active states. Density functional theory calculations further verified that, compared to the pyridinic N-coordinated low-spin Fe 2+ N 4 , the pyrrolic N-coordinated high-spin Fe 2+ exhibits optimized adsorption energy for reaction intermediates, thereby lowering the energy barrier of the rate-determining step (RDS) and facilitating OH* desorption. In conclusion, this work provides both experimental and theoretical evidence of the true active site in the FeNC during acidic ORR, offering significant insights for the rational design of high-performance materials for acidic fuel cells.

FeNC catalyst↗
DOE OSTI2023

Prediction of the structures and heats of formation of MO 2 , MO 3 , and M 2 O 5 for M = V, Nb, Ta, Pa

Structures for the mono-, di-, and tri-bridge isomers of M 2 O 5 as well as those for the MO 2 and MO 3 fragments for M = V, Nb, Ta, and Pa were optimized at the density functional theory (DFT) level. Single point CCSD(T) calculations extrapolated to the complete basis set (CBS) limit at the DFT geometries were used to predict the energetics. Here, the lowest energy dimer isomer was the di-bridge for M = V and Nb and the tri-bridge for M = Ta and Pa. The di-bridge isomers were predicted to be composed of MO 2 + and MO 3 - fragments, whereas the mono- and tri-bridge are two MO 2 + fragments linked by an O 2 - . The heats of formation of M 2 O 5 dimers, as well as MO 2 and MO 3 neutral and ionic species were predicted using the Feller–Peterson–Dixon (FPD) approach. The heats of formation of the MF 5 species were calculated to provide additional benchmarks. Dimerization energies to form the M 2 O 5 dimers are predicted to become more negative going down group 5 and range from -29 to -45 kcal mol -1 . The ionization energies (IEs) for VO 2 and TaO 2 are essentially the same at 8.75 eV whereas the IEs for NbO 2 and PaO 2 are 8.10 and 6.25 eV, respectively. The predicted adiabatic electron affinities (AEAs) range from 3.75 eV to 4.45 eV for the MO 3 species and vertical detachment energies from 4.21 to 4.59 eV for MO 3 - . The calculated M[double bond, length as m-dash]O bond dissociation energies increase from 143 kcal mol -1 for M = V to ~170 kcal mol -1 for M = Nb and Ta to ~200 kcal mol -1 for M = Pa. The M–O bond dissociation energies are all similar ranging from 97 to 107 kcal mol -1 . Natural bond analysis provided insights into the types of chemical bonds in terms of their ionic character. Pa 2 O 5 is predicted to behave like an actinyl species dominated by the interactions of approximately linear PaO 2 + groups.

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

High entropy powering green energy: hydrogen, batteries, electronics, and catalysis

A reformation in energy is underway to replace fossil fuels with renewable sources, driven by the development of new, robust, and multi-functional materials. High-entropy materials (HEMs) have emerged as promising candidates for various green energy applications, having unusual chemistries that give rise to remarkable functionalities. This review examines recent innovations in HEMs, focusing on hydrogen generation/storage, fuel cells, batteries, semiconductors/electronics, and catalysis—where HEMs have demonstrated the ability to outperform state-of-the-art materials. We present new master plots that illustrate the superior performance of HEMs compared to conventional systems for hydrogen generation/storage and heat-to-electricity conversion. We highlight the role of computational methods, such as density functional theory and machine learning, in accelerating the discovery and optimization of HEMs. The review also presents current challenges and proposes future directions for the field. We emphasize the need for continued integration of modeling, data, and experiments to investigate and leverage the underlying mechanisms of the HEMs that are powering progress in sustainable energy.

batteries↗
DOE OSTI2026

Toward Chemical Accuracy for Chemi- and Physisorption with an Efficient Density Functional

Understanding molecular adsorption on surfaces underpins many problems in chemistry and materials science. Accurately and efficiently describing the adsorption has been a challenging task for first-principles methods as the process can involve both short-range chemical bond formations and long-range physical interactions, e.g., van der Waals (vdW) interaction. Density functional theory presents an appealing choice for modeling adsorption reactions, although calculations with many exchange-correlation density functional approximations struggle to accurately describe both chemical and physical molecular adsorptions. Here, we propose an efficient density functional approximation that is accurate for both chemical and physical adsorption by concurrently optimizing its semilocal component and the long-range vdW correction against the prototypical adsorption CO/Pt(111) and Ar 2 binding energy curve. The resulting function opens the door to accurate and efficient modeling of general molecular adsorption.

adsorption↗
DOE OSTI2020

Dopant levels in large nanocrystals using stochastic optimally tuned range-separated hybrid density functional theory

In this work, we apply a stochastic version of an optimally tuned range-separated hybrid functional to provide insight on the electronic properties of P- and B- doped Si nanocrystals of experimentally relevant sizes. We show that we can use the range-separation parameter for undoped systems to calculate accurate results for dopant activation energies. We apply this strategy for tuning functionals to study doped nanocrystals up to 2.5 nm in diameter at the hybrid functional level. In this confinement regime, the P- and B- dopants have large activation energies and have strongly localized states that lie deep within the energy gaps. Structural relaxation plays a greater role for B-substituted dopants and contributes to the increase in activation energy when the B dopant is near the nanocrystal surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

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

MCML: Combining physical constraints with experimental data for a multi-purpose meta-generalized gradient approximation

The predictive power of density functional theory for materials properties can be improved without increasing the overall computational complexity by extending the generalized gradient approximation (GGA) for electronic exchange and correlation to density functionals depending on the electronic kinetic energy density in addition to the charge density and its gradient, resulting in a meta-GGA. Here, we propose an empirical meta-GGA model that is based both on physical constraints and on experimental and quantum chemistry reference data. Overall, the resulting optimized meta-GGA MCML yields improved surface and gas phase reaction energetics without sacrificing the accuracy of bulk property predictions of existing meta-GGA approaches.

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

The optimal co-doping of SrFe 1-x Co x O 3-δ oxygen carriers in redox applications

Although the oxygen carrier SrCoO 3 has higher redox activity than SrFeO 3 , cobalt is both more expensive and scarcer than iron, which would hinder the wide implementation of SrCoO 3 . For these reasons, doping SrFeO 3 with Co is a potential compromise, benefitting the redox properties of SrFeO 3 , while still limiting the overall amount of cobalt being used. In this work, to find the optimal level of Co-doping, density functional theory calculations were performed to investigate the Co-doping effect on the oxygen vacancy formation and oxygen migration in SrFe 1- x Co x O 3- δ ( x = 0, 0.125, 0.25, 0.375, 0.5). Our findings show that the oxygen vacancy formation energies ( E f ) decrease with the increase of Co content resulting from the increased composition of the O-2p band at the Fermi level upon Co doping. In particular, the E f decreases nearly 0.5 eV between the x = 0 and x = 0.25 samples while E f only decreases 0.1 eV further as Co content is increased to x = 0.5. We obtain that x = 0.25 is an optimal cost/benefit ratio for Co doping, which is preserved at both low oxygen vacancy concentrations ( δ = 0.0625 values listed above) and at high concentrations of δ = 0.1875 and 0.375. Kinetically, the oxygen migration barrier has slight change upon Co doping due to the similar size of Co and Fe. Therefore, considering both redox activity and economics in reversible oxygen storage applications, x = 0.25 is suggested as the optimal Co-doping value in SrFe 1- x Co x O 3- δ .

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 OSTI2024

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Comparing the Expense and Accuracy of Methods to Simulate Atomic Vibrations in Rubrene

Atomic vibrations can inform about materials properties from hole transport in organic semiconductors to correlated disorder in metal–organic frameworks. Currently, there are several methods for predicting these vibrations using simulations, but the accuracy–efficiency tradeoffs have not been examined in depth. In this work, rubrene is used as a model system to predict atomic vibrational properties using six different simulation methods: density functional theory, density functional tight binding, density functional tight binding with a Chebyshev polynomial-based correction, a trained machine learning model, a pretrained machine learning model called ANI-1, and a classical forcefield model. The accuracy of each method is evaluated by comparison to the experimental inelastic neutron scattering spectrum. All methods discussed here show some accuracy across a wide energy region, though the Chebyshev-corrected tight-binding method showed the optimal combination of high accuracy with low expense. We then offer broad simulation guidelines to yield efficient, accurate results for inelastic neutron scattering spectrum prediction.

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

Binding of uranyl cations to a Zr-based metal-organic framework by density functional theory

Zr-based metal–organic frameworks (Zr-MOFs) have been widely used as ion adsorbents for the removal or extraction of toxic and/or radionuclide species from aqueous solutions. However, the mechanisms by which uranyl cations (UO 2 2+ ) interact with Zr-MOFs have not been established. In this work, the nature of the bonding of uranyl cations with a Zr-MOF was determined using density functional theory for nineteen structurally distinct candidate complexes. Further, the results showed that in all cases the binding energy was of the order of 1.5 eV, but depended on the specific bonding site. The most stable structure involved coordination of the uranyl cation and two structurally distinct oxygens in the Zr-MOF metal node. It was also found that higher degree of deprotonation in Zr-MOF correlated with higher binding energy between the Zr-MOF and uranyl cations. These insights can aid in the design of Zr-MOFs with optimized features for efficient capture of uranyl cations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

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

Phase Diagrams of Alloys and Their Hydrides via On-Lattice Graph Neural Networks and Limited Training Data

Efficient prediction of sampling-intensive thermodynamic properties is needed to evaluate material performance and permit high-throughput materials modeling for a diverse array of technology applications. To alleviate the prohibitive computational expense of high-throughput configurational sampling with density functional theory (DFT), surrogate modeling strategies like cluster expansion are many orders of magnitude more efficient but can be difficult to construct in systems with high compositional complexity. We therefore employ minimal-complexity graph neural network models that accurately predict and can even extrapolate to out-of-train distribution formation energies of DFT-relaxed structures from an ideal (unrelaxed) crystallographic representation. This enables the large-scale sampling necessary for various thermodynamic property predictions that may otherwise be intractable and can be achieved with small training data sets. Two exemplars, optimizing the thermodynamic stability of low-density high-entropy alloys and modulating the plateau pressure of hydrogen in metal alloys, demonstrate the power of this approach, which can be extended to a variety of materials discovery and modeling problems.

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

Exploring Structure-Sensitive Relations for Small Species Adsorption Using Machine Learning

Accurate prediction of adsorption energies on heterogeneous catalyst surfaces is crucial to predicting reactivity and screening materials. Adsorption linear scaling relations have been developed extensively but often lack accuracy and apply to one adsorbate and a single binding site type at a time. These facts undermine their ability to predict structure sensitivity and optimal catalyst structure. Using machine learning on nearly 300 density functional theory calculations, we demonstrate that generalized coordination number scaling relations hold well for oxygen- and high-valency carbon-binding species but fail for others. Here we reveal that the valency and the electronic coupling of a species with the surface, along with the site type and its coordination environment, are critical for small species adsorption. The model simultaneously predicts the adsorption energy and preferred site and significantly outperforms linear scalings in accuracy. It can expose the structure sensitivity of chemical reactions and enable enhanced catalyst activity via engineering particle shape and facet defects. The generality of our methodology is validated by training the model with transition metal data and transferring it to predict adsorption energies on single-atom alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗