Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Amorphous-crystalline transition-driven synthesis of Co single-atom catalysts on MoO 3 for enhanced hydrogen evolution in acidic and alkaline media

Single atom catalysts (SACs) dispersed in metal oxide supports offer not only maximized catalyst utilization but also extended modulation capabilities through interactions between the SA and its support. However, achieving facile preparation and ensuring electrochemical durability of such SACs, particularly for hydrogen evolution reaction (HER) applications, remain formidable challenges. Here, our study addresses these issues by employing an amorphous-to-crystalline phase transition in MoO 3 to synthesize a cobalt (Co) SA catalyst. This method facilitates the production of highly selective SA catalysts at low temperatures and ensures their enhanced stability in HER applications. The Co SA MoO 3 catalyst exhibits superior performance in HER, operating effectively in both acidic and alkaline environments. Significantly, it maintains stable HER activity across these diverse electrolytic conditions. Our Density Functional Theory (DFT) calculations provide insights into the exceptional HER performance of Co SA MoO 3 . These calculations reveal that the strong affinity for hydrogen and water, facilitated by the modulation of the p-band orbitals at specific oxygen sites adjacent to the Co SA, establishes a thermoneutral pathway for HER. This study represents a pivotal advancement, showcasing a highly practical and robust single-atom catalyst, marking a significant stride towards sustainable energy solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Transient kinetic insights into selective propene oxidation over industrial bismuth molybdate catalysts

Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.

42 - ENGINEERING↗
DOE OSTI2026

Elucidating interfacial active sites in ruthenium–boron nitride nanotube catalysts for efficient low-temperature ammonia-to-hydrogen conversion

Tailoring the interaction between metal nanoparticles and catalyst support presents a prominent strategy to enhance both the activity and durability in hydrogen (H 2 ) production catalysts. In this work, ruthenium nanoparticles (NPs) supported on boron nitride nanotubes (Ru/BNNT) are introduced as efficient and thermally robust catalysts for low-temperature ammonia (NH 3 ) decomposition. The unique curvature and ionic nature of BNNTs enable uniform Ru dispersion and metal-support interactions (MSIs), resulting in exceptional H 2 generation efficiency and long-term operational stability. In-situ transmission electron microscopy (TEM) reveals remarkable thermal resistance of Ru/BNNT with minimal nanoparticle sintering, while density functional theory (DFT) calculations uncover a dual-site mechanism in which interfacial Ru atoms promote NH 3 dissociation and adjacent Ru sites facilitate 2H* recombination and H 2 desorption. This cooperative interaction between metal NPs and the BNNT support underpins the outstanding catalytic performance and durability observed. In conclusion, the findings highlight the strategic potential of BNNTs as versatile supports for high-performance and stable catalysts in sustainable H 2 energy conversion and related catalytic processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The role of the metal core in the performance of WO x inverse catalysts

In this study, metal-metal oxide inverse catalysts are valuable in diverse applications but selecting pairs is challenging. Here we investigate the role of metal cores (Pt, Ru, Rh, Pd, Ni) in the formation and dynamics of Brønsted acid sites on WO x overlayer on carbon. By in situ/ex situ characterizations, probe chemistry kinetics, and DFT calculations, we demonstrate that metals with varying work functions promote the WO x dispersion, enhancing the Brønsted acidity than the 3D crystallites on carbon at submonolayer coverages by tuning the oxidation state of WO x . The metal core impacts the water splitting and hydroxylation of WO x . Water splitting on Ru-WO x is thermodynamically favorable compared to Ni-WO x and increases the dehydration rate on the former, leading to an activity enhancement by increasing the H 2 pulsing frequency. This study provides insights into optimizing inverse catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhancing carbon nanotube production from carbon dioxide and ethane using bimetallic catalysts

Converting CO 2 into carbon nanotubes (CNTs) offers a promising way for CO 2 utilization and sequestration, potentially mitigating environmental impacts from anthropogenic emissions. This study reports that bimetallic CoFe catalysts can increase CNT production from the reaction of CO 2 and C 2 H 6 by an order of magnitude compared to their monometallic counterparts. The active sites and CNT morphologies are composition-dependent: Co-rich catalysts (Co/Fe ratio ≥ 5) form stable face-centered cubic (fcc) CoFe alloys, producing cylindrical CNTs; Fe-containing catalysts (Co/Fe ≤ 2) favor body-centered cubic (bcc) CoFe alloy upon reduction, which transforms into carbides, resulting in bamboo-like CNTs. Experimental evidence and DFT calculations reveal that adjacent Fe and Co atoms modulate CO and C x H y adsorption, regulating CNT production pathways through the CO Boudouard reaction and C 2 H 6 decomposition. In conclusion, these results highlight the dual benefits of bimetallic catalysts in enhancing CNT yield and controlling CNT morphology through adjustment of catalyst compositions.

58 GEOSCIENCES↗
DOE OSTI2021

Machine learning elastic constants of multi-component alloys

The present manuscript explores application of machine learning methods for determining elastic constants and other derived mechanical properties of multi-component alloys. Here, a number of machine learning models, including linear regression, neural network and random forest based models, are trained and tested on a dataset of binary alloys generated using density functional theory (DFT) calculations and spanning over a large number of elemental species in the periodic table. Starting with a wide range of simple and easily accessible compositionally-averaged elemental features, a correlation-based feature selection strategy was used to systematically down-select a set of most relevant features towards the prediction of the elasticity tensor components. The true predictive performance and the associated uncertainties of the models were established by testing on unseen data and bootstrapping, respectively. A single and pair-wise feature partial dependence analysis was performed to visualize the average property trends in the multi-dimensional feature space in order to further understand the achieved predictive performance. The utility of the trained model is further demonstrated by obtaining sufficiently accurate yet highly efficient approximations for bulk modulus, Young’s modulus, shear modulus and Poisson’s ratio for alloys beyond the binary space (i.e., two-component alloys) on which the model was originally trained. More importantly, we test and validate the predictive performance of the developed model directly against the experimentally measured elastic constants of technologically relevant multi-component alloys (such as, Ni- and Ti-based alloys). Finally, utility of such a data-enabled route is demonstrated by predicting the possible range of various elastic properties for vast composition space available within the five component Ni-Cr-Fe-Mo-W alloy system in a high-throughput manner.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of phase stability and chemical segregation in the Mo-V alloys using a generalized embedded atom method potential

A new interatomic potential for the Mo-V system is introduced to facilitate the study of phase stability and mechanical properties at lower temperatures. This potential is based on a generalization of the embedded atom method and includes contributions from embedding energy, explicit two- and three-body interactions and nonlocal many-body interaction terms. The parameters of the potential are optimized by using data from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, such as elastic constants, equation of states, phonon dispersion curves, point defect properties and melting temperatures for different compositions. Even though our potential is trained on a small dataset, its accuracy is comparable to available machine learning potentials for Mo and V. Furthermore, our results show that an ordered B2 phase is stable at low temperatures in alloys containing 50% V, but the solid solution phase is stable above 800 K. However, such long-range ordering is not observed in V-rich or Mo-rich alloys. In addition, our results show that V segregates to dislocation cores and grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2020

GPU-acceleration of the ELPA2 distributed eigensolver for dense symmetric and hermitian eigenproblems

The solution of eigenproblems is often a key computational bottleneck that limits the tractable system size of numerical algorithms, among them electronic structure theory in chemistry and in condensed matter physics. Large eigenproblems can easily exceed the capacity of a single compute node, thus must be solved on distributed-memory parallel computers. We here present GPU-oriented optimizations of the ELPA two-stage tridiagonalization eigensolver (ELPA2). On top of cuBLAS-based GPU offloading, we add a CUDA kernel to speed up the back-transformation of eigenvectors, which can be the computationally most expensive part of the two-stage tridiagonalization algorithm. Furthermore, we benchmark the performance of this GPU-accelerated eigensolver on two hybrid CPU–GPU architectures, namely a compute cluster based on Intel Xeon Gold CPUs and NVIDIA Volta GPUs, and the Summit supercomputer based on IBM POWER9 CPUs and NVIDIA Volta GPUs. Consistent with previous benchmarks on CPU-only architectures, the GPU-accelerated two-stage solver exhibits a parallel performance superior to the one-stage counterpart. Finally, we demonstrate the performance of the GPU-accelerated eigensolver developed in this work for routine semi-local KS-DFT calculations comprising thousands of atoms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Radionuclide incorporation in negative thermal expansion α-Zr(WO 4 ) 2 : A density functional theory study

The incorporation of uranium, plutonium and technetium in the negative thermal expansion (NTE) α-Zr (WO 4 ) 2 has been investigated within the framework of density functional theory (DFT). It is found that the vacancy formation energies of the charged vacancies are overall larger than that of its counterpart neutral Frenkel defects and Schottky defects. DFT calculations suggest that U and Pu substitutions for the Zr site are preferred in α-Zr (WO 4 ) 2 . In case of Tc substitution, both Tc(IV) for the Zr site and Tc(VII) for the W site are considered under oxygen-poor and oxygen-rich conditions, while Tc(VII) substitution can be improved significantly by including Y 2 O 3 (charge compensation).

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

Azo–hydrazone tautomerism in organometallic complexes triggered by a -Re(CO) 3 (L) core: A spectroscopic and theoretical study

The spectroscopic properties and tautomeric behavior of five novel Rhenium (I) tricarbonyl complexes bearing an azo ligand are presented. The organic ligand is stable in solution as the hydrazone tautomer. It remains as hydrazone in different medium conditions (solvent, concentration, pH, etc.) without the formation of detectable amounts of the azo tautomer. However, the complexation of this ligand to the strong electron-withdrawing fragment -Re(CO) 3 X (X = Cl – , Br – ) causes tautomerism to appear in the organic moiety. Two well-defined regions in the electronic spectra for both tautomers were observed, allowing the azo/enol-keto/hydrazone equilibrium to be followed. TD-DFT calculations indicate that for the keto/hydrazone form, the main absorption band is attributed to an IL transition. In contrast, the azo/enol species shows a major contribution of the 1 MLCT (dπ(Re) → dπ*(NN)) transition. The tautomeric equilibrium is easily shifted using solvents with different dielectric constant and hydrogen bond donor/acceptor (HBD/HBA) abilities. Also, the tautomerization process is deeply influenced by the electronic properties of the axial ligand. Therefore, this behavior represents a different strategy for the design of novel materials with optical properties. Furthermore, tautomerism is affected by the concentration of complexes and the presence of water in solvents. Finally, the acid-base behavior of these compounds in a 50% v/v ethanol-buffer system showed that the azo/enol tautomer is stabilized in acidic media. At the same time, an increase of pH promotes tautomerization toward the keto/hydrazone, followed by forming the anionic form at pH > 8.

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

The magnetic order in multiferroic DyMnO 3

We report that with flexibility in tuning their electric and magnetic properties, multiferroics can be used in information exchange and storage in ways that are very different from the present electronic materials. Here we use resonant soft X-ray scattering spectroscopy to study the F-type (0,τ,0) and C-type (0, 1–2τ, 0) diffraction peaks from sinusoidal antiferromagnetic spin order in multiferroic DyMnO 3 . By comparing the temperature dependence of ordering wave vectors τ, peak intensities I, and correlation lengths λ measured at Mn L 2 - , O K-, and Dy M 5 -edges, we show that the nearly perfect locking between the ordering wave vectors from Dy 4f states and Mn 3d orbitals manifesting the second harmonic diffraction peak implies the notable orbital involvement in the coupling between Mn and Dy spins. Our DFT calculations further suggest that the lattice response to different antiferromagnetic ground states (A-type versus E-type) is much weaker in TbMnO 3 , in agreement with previous claim that the symmetric exchange interaction can be an important factor for understanding the ferroelectricity in DyMnO 3 than in TbMnO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lithium superoxide-based high rate Li-Air batteries enabled by Di-iridium sulfur bridge active sites

Li-oxygen (Li-O 2 ) batteries can potentially provide much higher energy density than Li-ion batteries; however, the practical application of these batteries is hindered due to several drawbacks such as low current rates and high overpotential for the charging process. Here, in this paper, we report a novel Li-Air battery system that operates under high current rates (up to 1mAcm -2 ) with LiO 2 as the primary discharge product instead of the commonly reported Li 2 O 2 . This LiO 2 based battery at high rates is through a combination of an as-synthesized new one-dimensional (1D) transition metal trichalcogenide mid-entropy alloy of SnIrS 3.6 as a cathode catalyst and an electrolyte blend with a SnI 2 bi-functional additive. It is revealed that SnIrS 3.6 has a microporous structure composed of six- and five-coordinated metal atoms, forming octahedral and triangular bipyramids which has not been observed in other layered chalcogeide materials. DFT calculations reveal that the SnIrS 3.6 structure can result in LiO 2 formation through di-iridium sulfur bridge active sites that results in strong binding of O 2 and LiO 2 preventing disproportionation to Li 2 O 2 and enabling high rates. This finding will open a new perspective in designing advanced LiO 2 -based Li-O 2 batteries for real practices.

25 ENERGY STORAGE↗
DOE OSTI2022

Magnetocaloric response with significant mechanical efficiency in frustrated intermetallic compound $\mathrm{Pr_{2}Co_{0.86}Si_{2.88}}$

We report magnetically frustrated materials are considered as a promising unconventional members of caloric materials. Here, the magnetocaloric properties of the frustrated Pr 2 Co 0.86 Si 2.88 have been investigated and discussed with the aid of density functional theory (DFT) calculations. The material exhibits a magnetic entropy change (-ΔS M ) of 13.1 J/kg-K for ΔH = 70 kOe around low-temperature transition T L ~ 4 K associated with the antiferromagnetic coupling between localized Pr-4f and itinerant Co-3d moments. Despite the absence of any long-range magnetic order, the obtained - ΔS M is one of the highest among the known Pr-based good magnetocaloric materials at cryogenic temperature range. It also exhibits a relative cooling power (RCP) of ~ 201 J/kg and an adiabatic temperature change of 6.3 K around T L at 70 kOe. Moreover, the compound also exhibits a high mechanical efficiency and moderate electrical efficiency, being beneficial for possible technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cubic to hexagonal tuning in Fe 2 Mn(Si 1– x Ge x ) Heusler alloys

Here, the competition between the stability of the cubic and hexagonal full Heusler alloys and the implications concerning their magnetic properties were systematically studied through the detailed structural and magnetic characterization of the Fe 2 Mn(Si 1– x Ge x ) system. This system was specifically chosen as the parent compositions are cubic ( x = 0) and hexagonal ( x = 1). It is found that the formation of hexagonal phases occurs for the x ≥ 0.6 samples, whereas its phase fraction monotonically increases with x until the pure hexagonal Fe 2 MnGe is formed. The change in structure results in high sensitiveness of both the saturation of magnetization ($M_S$) and Curie temperature ($T_C$) with x values, related to a strong magnetocrystalline anisotropy of the hexagonal phase. Both cubic and hexagonal magnetic features were qualitatively reproduced by Density Functional Theory (DFT) calculations. This work provides an experimental and theoretical foundation for further design of Heusler systems with controlled structures and magnetic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigating the non-classical M-H 2 bonding in OsClH 3 (PPh 3 ) 3

Dihydrogen ligation to metal centers results in a wide range of metal-hydrogen interactions. Compounds containing non-classical dihydrogen ligands in which the H-H contacts are significantly longer than the 0.74 A distance found in solid H 2 are classified as Kubas-like compounds, stretched dihydrogen ligands, or compressed dihydrides, depending on the H-H contact length. Compounds containing stretched dihydrogen and compressed dihydride ligands are somewhat rare. Developing a more complete understanding of the metal-hydrogen and hydrogen-hydrogen interactions in these species may provide insights into the mechanisms for the oxidative addition reaction and hydrogen storage. Here, we use diffraction and inelastic neutron scattering (INS) measurements, paired with both molecular and solid-state density functional theory (DFT), to fully characterize the nature of the metal-hydrogen and hydrogen-hydrogen interactions in a previously reported compound, OsClH 3 (PPh 3 ) 3 . In that report, the H-H distance was measured as 1.48(2) A by single-crystal neutron diffraction and was described as a "stretched" dihydrogen. The INS-generated phonon density of states is well described by the DFT calculations, including those normal modes dominated by the vibrational motions of the three Os-ligated hydrogen atoms. Additionally, the resulting calculated electron charge density landscape indicates that there is no remnant bonding character between the Os-ligated hydrogen atoms in the compound. Based on these findings, we update the understanding of the metal-hydrogen bonding in this complex and reclassify it as a compressed dihydride. As such, OsClH 3 (PPh 3 ) 3 represents a rare example of a non-classical compressed dihydride species.

08 HYDROGEN↗
DOE OSTI2020

Oxide-catalyzed self- and cross-condensation of cycloketones. Kinetically relevant steps that determine product distribution

Metal oxides with acid-base properties can be used to catalyze aldol condensation and potentially obtain valuable chemicals from biomass-derived oxygenates. Understanding the reaction mechanisms could enable the development of effective catalysts with controllable selectivity. In this work, we investigate the reaction parameters that affect the kinetics and product distribution of self- and cross-aldol condensation of cyclopentanone (CPO) and cyclohexanone (CHO) over CeO 2 and ZrO 2 catalysts by combining detailed analysis of experiments and DFT calculations. It is well-known that either the deprotonation of an α-C–H that generates an enolate or the C–C coupling that produces the dimers can be rate-limiting step. Indeed, the detailed kinetics analysis shows that the observed data can be fitted by unimolecular or bimolecular Langmuir-Hinshelwood models. Specifically, on the ZrO 2 catalyst, the unimolecular α-C–H activation is the rate limiting step for the CPO self-condensation, while the bimolecular C–C coupling limits the CHO self-condensation. By contrast, on CeO 2 , the bimolecular C–C coupling step seems to be rate-limiting for both CPO and CHO self-condensation reactions. The adsorption parameters derived from kinetics fitting indicate that when the two reactants are co-fed, the surface is preferentially covered by CHO. So, the condensation rate of [CPO]-activated products ([CPO]CPO + [CPO]CHO) decreases in the presence of CHO, due to the competitive adsorption. However, despite a higher surface coverage, CHO is not an effective electrophile because it imposes stronger steric constraints to the formation of the C–C bond than CPO. Therefore, for a given enolate, the product distribution shows that [CPO]CPO > [CPO]CHO and, similarly [CHO]CPO > [CHO]CHO. Likewise, the yields of [CPO]- and [CHO]-activated products in mixed CPO/CHO feeds are similar because the higher surface coverage and stability of the [CHO] enolate is counterbalanced by the inhibition of C–C coupling by steric constraints when CHO is involved. In summary, the overall rank of products observed experimentally can be explained in terms of a combination of the following effects: relative surface coverage of the reactants, steric constraints for the C–C coupling, stability of the enolate, and effectiveness of the electrophile to accept an electron.

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

Prussian blue analogues as platform materials for understanding and developing oxygen evolution reaction electrocatalysts

Transition metal based materials containing Fe have drawn great attention as oxygen evolution reaction (OER) catalysts. The nature of the electrocatalytic active species remains under debate due to the ambiguous physicochemical properties of the catalyst materials, such as the oxidation states and crystal structures. Here, in order to address this issue, transition metal Prussian blue analogues (TM-PBA, Na(TM)(Fe)(CN) 6 , TM = V, Fe, Co, and Ni) with an isomorphous structure are investigated for OER catalysis. Our combined experimental measurements and density functional theory (DFT) calculations reveal that TM-PBAs exhibit volcano-like OER activity with Ni-PBA located near the top of the volcano. Such a volcano-like activity profile can be attributed to the distinctive binding energy difference between *O and *OH on different TM-PBAs surfaces. This research demonstrates that TM-PBAs can be used as platform materials for understanding structure-property-activity relationships in OER catalysts.

25 ENERGY STORAGE↗