Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Elucidation of Active Sites in Aldol Condensation of Acetone over Single-Facet Dominant Anatase TiO 2 (101) and (001) Catalysts

Aldol condensations of carbonyl compounds for C–C bond formation are a very important class of reactions in organic synthesis and upgrading of biomass-derived feedstocks. However, the atomic level understanding of reaction mechanisms and structure-activity correlation on widely used transition metal oxide catalysts are limited due to the high degree of struc-tural heterogeneity of catalysts such as commercial TiO 2 powders. Here, we provide a deep understanding of the reaction mechanisms, kinetics, and structure–function relationships for vapor phase acetone aldol condensation through the con-trolled synthesis of two catalysts with high surface areas and clean, dominant facets, coupled with detailed characterization and kinetic studies that are further assisted by density functional theory (DFT) calculations. Temperature-dependent diffuse reflectance infrared Fourier transform spectroscopy showed the existence of abundant acetone bonded to surface hydroxyl groups (acetone-O s H) and acetone bonded to Lewis acid sites (acetone-Ti 5c ) on surface of both {101} and {001} facet dom-inant TiO 2 . Intermolecular C–C coupling of enolate intermediate from acetone-Ti 5c and a vicinal acetone-O s H is a kinetically relevant step, which is consistent with kinetic and isotopic studies and DFT calculations. The {001} facet showed a lower apparent activation energy (or higher reactivity) than the {101} facet. This is likely caused by the weaker Lewis acid and Brønsted base strengths of the {001} facet which favors the reprotonation–desorption of the coupled intermediate, making the C–C coupling step more exothermic on the {001} facet and resulting in an earlier transition state with a lower activation barrier. It is also possible that the {001} facet has a smoother surface configuration and less steric hindrance during intermo-lecular C–C bond formation than the {101} facet.

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

Electrochemical CO 2 Reduction Reaction over Cu Nanoparticles with Tunable Activity and Selectivity Mediated by Functional Groups in Polymeric Binder

Electrochemical carbon dioxide reduction reaction (CO 2 RR) using copper (Cu)-based catalysts has received significant attention mainly because Cu is an element capable of producing hydrocarbons and oxygenates. One possible way to control the CO 2 RR performance at the electrode interface is by modifying catalysts with specific functional groups of different polymeric binders, which are necessary components in the process of electrode fabrication. However, the modification effect of the key functional groups on the CO 2 RR activity and selectivity is poorly understood over Cu-based catalysts. In this work, the role of functional groups (e.g., -COOH and -CF 2 groups) in hydrophilic and hydrophobic polymeric binders on the CO 2 RR of Cu-based catalysts is investigated using a combination of electrochemical measurements, in-situ characterization and density functional theory (DFT) calculations. DFT results reveal that functional groups influence the binding energies of key intermediates involved in both CO 2 RR and the competing hydrogen evolution reaction, consistent with experimental observation of binder-dependent product distributions among formic acid, CO, CH 4 , and H 2 . This study provides a fundamental understanding that the selection of desired polymeric binders is a useful strategy for tuning the CO 2 RR activity and selectivity.

25 ENERGY STORAGE↗
DOE OSTI2022

Di- and Tetrameric Molybdenum Sulfide Clusters Activate and Stabilize Dihydrogen as Hydrides

NaY zeolite-encapsulated dimeric (Mo 2 S 4 ) and tetrameric (Mo 4 S 4 ) molybdenum sulfide clusters stabilize hydrogen as hydride binding to Mo atoms. Density functional theory (DFT) calculations and adsorption measurements suggest that stabilization of hydrogen as sulfhydryl (SH) groups, as typical for layered MoS 2 , is thermodynamically disfavored. Competitive adsorption of H 2 and ethene on Mo was probed by quantifying adsorbed CO on partly hydrogen and/or ethene covered samples with IR spectroscopy. During hydrogenation, experiment and theory suggest that Mo is covered predominately with ethene and sparsely with hydride. DFT calculations further predict that under reaction conditions, each Mo x S y cluster can activate only one H 2 , suggesting that the entire cluster (irrespective of its nuclearity) acts as one active site for hydrogenation. The nearly identical turnover frequencies (24.7 ± 3.3 mol ethane ·h -1 ·mol cluster -1 ), apparent activation energies (31-32 kJ·mol -1 ), and reaction orders (~0.5 in ethene and ~1.0 in H 2 ) show that the active sites in both clusters are catalytically indistinguishable.

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

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗
DOE OSTI2020

Prediction of organic homolytic bond dissociation enthalpies at near chemical accuracy with sub-second computational cost

Bond dissociation enthalpies (BDEs) of organic molecules play a fundamental role in determining chemical reactivity and selectivity. However, BDE computations at sufficiently high levels of quantum mechanical theory require substantial computing resources. In this paper, we develop a machine learning model capable of accurately predicting BDEs for organic molecules in a fraction of a second. We perform automated density functional theory (DFT) calculations at the M06-2X/def2-TZVP level of theory for 42,577 small organic molecules, resulting in 290,664 BDEs. A graph neural network trained on a subset of these results achieves a mean absolute error of 0.58 kcal mol-1 (vs DFT) for BDEs of unseen molecules. We further demonstrate the model on two applications: first, we rapidly and accurately predict major sites of hydrogen abstraction in the metabolism of drug-like molecules, and second, we determine the dominant molecular fragmentation pathways during soot formation.

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

In-situ structure and catalytic mechanism of NiFe and CoFe layered double hydroxides during oxygen evolution

NiFe and CoFe (MFe) layered double hydroxides (LDHs) are among the most active electrocatalysts for the alkaline oxygen evolution reaction (OER). Herein, we combine electrochemical measurements, operando X-ray scattering and absorption spectroscopy, and density functional theory (DFT) calculations to elucidate the catalytically active phase, reaction center and the OER mechanism. We provide the first direct atomic-scale evidence that, under applied anodic potentials, MFe LDHs oxidize from as-prepared α-phases to activated γ-phases. The OER-active γ-phases are characterized by about 8% contraction of the lattice spacing and switching of the intercalated ions. DFT calculations reveal that the OER proceeds via a Mars van Krevelen mechanism. The flexible electronic structure of the surface Fe sites, and their synergy with nearest-neighbor M sites through formation of O-bridged Fe-M reaction centers, stabilize OER intermediates that are unfavorable on pure M-M centers and single Fe sites, fundamentally accounting for the high catalytic activity of MFe LDHs.

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

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

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

Redox dynamics and surface structures of an active palladium catalyst during methane oxidation

Catalysts based on palladium are among the most effective in the complete oxidation of methane. Despite extensive studies and notable advances, the nature of their catalytically active species and conceivable structural dynamics remains only partially understood. Here, we combine operando transmission electron microscopy (TEM) with near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) and density functional theory (DFT) calculations to investigate the active state and catalytic function of Pd nanoparticles (NPs) under methane oxidation conditions. We show that the particle size, phase composition and dynamics respond appreciably to changes in the gas-phase chemical potential. In combination with mass spectrometry (MS) conducted simultaneously with in situ observations, we uncover that the catalytically active state exhibits phase coexistence and oscillatory phase transitions between Pd and PdO. Aided by DFT calculations, we provide a rationale for the observed redox dynamics and demonstrate that the emergence of catalytic activity is related to the dynamic interplay between coexisting phases, with the resulting strained PdO having more favorable energetics for methane oxidation.

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

Analytical ab initio hessian from a deep learning potential for transition state optimization

Identifying transition states—saddle points on the potential energy surface connecting reactant and product minima—is central to predicting kinetic barriers and understanding chemical reaction mechanisms. In this work, we train a fully differentiable equivariant neural network potential, NewtonNet, on thousands of organic reactions and derive the analytical Hessians. By reducing the computational cost by several orders of magnitude relative to the density functional theory (DFT) ab initio source, we can afford to use the learned Hessians at every step for the saddle point optimizations. We show that the full machine learned (ML) Hessian robustly finds the transition states of 240 unseen organic reactions, even when the quality of the initial guess structures are degraded, while reducing the number of optimization steps to convergence by 2–3× compared to the quasi-Newton DFT and ML methods. All data generation, NewtonNet model, and ML transition state finding methods are available in an automated workflow.

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

Insights into the mechanisms of NH3 inhibition on Cu-CHA SCR catalysts

This work elucidates the atomic-scale mechanism behind ammonia (NH3) inhibition during the selective catalytic reduction (SCR) of NO? on Cu-CHA catalysts, a key issue limiting low-temperature emission control. Using SCR kinetic analysis, operando electron paramagnetic resonance (EPR) spectroscopy, and density functional theory (DFT), we demonstrate that NH3 inhibition primarily slows the oxidation half-cycle (OHC), while the reduction half-cycle (RHC) remains unaffected. DFT simulations reveal that excess NH3 substantially increases the diffusion barrier for CuI ions, hindering formation of essential CuII-oxo dimer intermediates and thus suppressing OHC kinetics. Operando EPR studies confirm that this inhibition strongly depends on operating temperature and catalyst Cu loading. Our findings highlight strategies to counteract NH3 inhibition, including optimizing Cu loading, precisely managing NH3:NO feed ratios, and enhancing CuI ion mobility through tailored catalyst design and operational adjustments, thereby advancing the efficiency of emission control technologies in automotive applications.

Deka, Dhruba Jyoti↗
DOE OSTI2022

A flexible and scalable scheme for mixing computed formation energies from different levels of theory

Abstract Computational materials discovery efforts are enabled by large databases of properties derived from high-throughput density functional theory (DFT), which now contain millions of calculations at the generalized gradient approximation (GGA) level of theory. It is now feasible to carry out high-throughput calculations using more accurate methods, such as meta-GGA DFT; however recomputing an entire database with a higher-fidelity method would not effectively leverage the enormous investment of computational resources embodied in existing (GGA) calculations. Instead, we propose here a general procedure by which higher-fidelity, low-coverage calculations (e.g., meta-GGA calculations for selected chemical systems) can be combined with lower-fidelity, high-coverage calculations (e.g., an existing database of GGA calculations) in a robust and scalable manner. We then use legacy PBE(+ U ) GGA calculations and new r 2 SCAN meta-GGA calculations from the Materials Project database to demonstrate that our scheme improves solid and aqueous phase stability predictions, and discuss practical considerations for its implementation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗
DOE OSTI2024

Electronic structure prediction of multi-million atom systems through uncertainty quantification enabled transfer learning

The ground state electron density — obtainable using Kohn-Sham Density Functional Theory (KS-DFT) simulations — contains a wealth of material information, making its prediction via machine learning (ML) models attractive. However, the computational expense of KS-DFT scales cubically with system size which tends to stymie training data generation, making it difficult to develop quantifiably accurate ML models that are applicable across many scales and system configurations. Here, we address this fundamental challenge by employing transfer learning to leverage the multi-scale nature of the training data, while comprehensively sampling system configurations using thermalization. Our ML models are less reliant on heuristics, and being based on Bayesian neural networks, enable uncertainty quantification. We show that our models incur significantly lower data generation costs while allowing confident — and when verifiable, accurate — predictions for a wide variety of bulk systems well beyond training, including systems with defects, different alloy compositions, and at multi-million-atom scales. Moreover, such predictions can be carried out using only modest computational resources.

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

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2026

Developing a complete AI-accelerated workflow for superconductor discovery

The quest to identify new superconducting materials with enhanced properties is hindered by the prohibitive cost of computing electron-phonon spectral functions, severely limiting the materials space that can be explored. Here, we introduce a Bootstrapped Ensemble of Equivariant Graph Neural Networks (BEE-NET), a machine-learning model trained to predict the Eliashberg spectral function and superconducting critical temperature with a mean-absolute-error of 0.87 K relative to DFT-based Allen-Dynes calculations. Intriguingly, BEE-NET achieves a true-negative-rate of 99.4%, enabling highly efficient screening for the rare property of superconductivity. Integrated into a multi-stage, AI-accelerated discovery pipeline that incorporates elemental-substitution strategies and machine-learned interatomic potentials, our workflow reduced over 1.3 million candidate structures to 741 dynamically and thermodynamically stable compounds with DFT-confirmed T c > 5 K. We report the successful synthesis and experimental confirmation of superconductivity in two of these previously unreported compounds. This study establishes a data-driven framework that integrates machine learning, quantum calculations, and experiments to systematically accelerate superconductor discovery.

Gibson, Jason B. [Quantum Formatics, Cambridge, MA↗
DOE OSTI2025

Grain boundary zirconia-modified garnet solid-state electrolyte

Here, we report a method for promoting electrochemical stability in garnet Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 solid-state electrolyte based on a composite two-phase oxide–oxide microstructure. Grain boundary precipitation of the controlled distribution of amorphous zirconium oxide microparticles is achieved through the addition of reactive tantalum carbide. During ambient-atmosphere sintering, the carbide decomposes through an in situ reaction, the ‘extra’ Ta substituting for Zr within the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 lattice. Density functional theory (DFT) calculations identify a thermodynamically favourable reaction path and show how substituting Ta 5+ at Zr 4+ sites affects the crystal structure as well as bulk ionic and electronic conductivities. Quantitative stereology highlights that zirconia also acts as a sintering aid, reducing compact porosity. Cryogenic focused-ion-beam scanning electron microscopy and fractography analysis of cycled solid-state electrolytes illustrates that near-universally observed intergranular Li-metal dendrite propagation is suppressed by the two-phase microstructure, favouring transgranular dendrites instead. Importantly, DFT demonstrates that compared with the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 surface, the zirconium oxide surface per se is less electronically conductive and does not trap excess electrons to reduce Li ions. This is a key reason for the substantial improvement in the electrochemical properties over the single-phase baseline.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the liquid–liquid phase transition of dense hydrogen

Until recently the consensus theory/computation interpretation of the challenging liquid liquid phase transition (LLPT) of high-pressure hydrogen was first order. Cheng et al. developed a machine-learnt potential (MLP) that, in larger molecular dynamics (MD) simulations, gives a continuous transition instead. Here, we show that the MLP does not reproduce our still larger MD density-functional theory (MD-DFT) calculations as it should. Since the MLP is not a faithful surrogate for the MD-DFT, the Ref. 6 prediction of a supercritical atomic liquid is unfounded.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

An automatically curated first-principles database of ferroelectrics

Ferroelectric materials have technological applications in information storage and electronic devices. The ferroelectric polar phase can be controlled with external fields, chemical substitution and size-effects in bulk and ultrathin film form, providing a platform for future technologies and for exploratory research. In this work, we integrate spin-polarized density functional theory (DFT) calculations, crystal structure databases, symmetry tools, workflow software, and a custom analysis toolkit to build a library of known, previously-proposed, and newly-proposed ferroelectric materials. With our automated workflow, we screen over 67,000 candidate materials from the Materials Project database to generate a dataset of 255 ferroelectric candidates, and propose 126 new ferroelectric materials. We benchmark our results against experimental data and previous first-principles results. The data provided includes atomic structures, output files, and DFT values of band gaps, energies, and the spontaneous polarization for each ferroelectric candidate. Finally, we contribute our workflow and analysis code to the open-source python packages atomate and pymatgen so others can conduct analogous symmetry driven searches for ferroelectrics and related phenomena.

36 MATERIALS SCIENCE↗