Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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

DFT + DMFT study of dopant effects in the heavy-fermion compound CeCoIn 5

Here we study the dopant-induced inhomogeneity effect on the electronic properties of heavy fermion CeCoIn 5 using a combined approach of density functional theory (DFT) and dynamical mean-field theory (DMFT). The inhomogeneity of the hybridization between Ce-4f and conduction electrons is introduced to impose the inequivalent Ce atoms with respect to the dopant. From the DFT to the DFT + DMFT results, we demonstrate a variation of the hybridization strength depending on the hole or electron doping. A drastic asymmetric mass renormalization could be reproduced in the DFT + DMFT calculation. Finally, the calculated Kondo temperature reflects the different development of the heavy quasiparticle states, depending on the dopant.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nonprecious Single Atom Catalyst for Methane Pyrolysis

The development of a suitable catalytic system for methane pyrolysis reactions requires a detailed investigation of the activation energy of C-H bonds on catalysts, as well as their stability against sintering and coke formation. In this work, both single-metal Ni atoms and small clusters of Ni atoms deposited on titanium nitride (TiN) plasmonic nanoparticles were characterized for the C-H bond activation of a methane pyrolysis reaction using ab initio spin-polarized density functional theory (DFT) calculations. The present work shows the complete reaction pathway, including energy barriers for C-H bond activation and dehydrogenated fragments, during the methane pyrolysis reaction on catalytic systems. Interestingly, the C-H bond activation barriers were low for both Ni single-atom and Ni-clusters, showing the energy barriers of ~1.10 eV and ~0.88 eV, respectively. Additionally, single-atom Ni-TiN showed weaker binding to adsorbates, and a net endothermic reaction pathway indicated that the single-atom Ni-TiN was expected to resist coke formation on its surface. However, these Ni single-atom catalysts can sinter, aggregate into a small cluster, and form a coke layer from the highly exothermic reaction pathway that the cluster takes despite the facile reaction pathway.

08 HYDROGEN↗
DOE OSTI2022

Quantifying Errors in Effective Cluster Interactions of Lattice Gas Cluster Expansions

The promise of lattice gas (LG) cluster expansions (CEs) is that they can describe a given system property to any level of accuracy since the orthogonal “cluster basis functions” span the complete space of available configurations. Such an approach can be constructed to an arbitrarily large surface of a finite number of distinct adsorption sites. Unfortunately, this is only true for the case of an ideal, fixed lattice decorated with components at precise lattice points (the lattice “sites”) with no distortions or relaxations subsequently allowed. Since most systems, and surfaces specifically, do not conform to such an ideal set of constraints, errors in LG CEs must be expected or CE convergence severely hampered. Beyond this, numerical errors in the provided data can complicate the proper construction of a truly predictive and/or physically significant CE. In this work, we show here how reliance on typical statistical tools like confidence intervals cannot be expected to provide an accurate representation of the uncertainty of effective cluster interactions (ECIs) in the CE due to the nature of the target ab initio data and the nature of CEs themselves. We develop a method for estimating these errors that does not rely on statistical assumptions about the model or data. We then use these ECI errors to quantify fundamental consequences on the uncertainty of ECIs in CEs built from O/Fe(100) data whose surface and adsorbates have been allowed to relax in a typical manner and from O/Fe(100) data whose surface and adsorbates are fixed in ideal lattice positions. We also quantify the effect of using a different density functional theory exchange–correlation functional, using these ECI errors to assess the significance in any deviations. In both cases, our method is shown to have remarkable utility in the quantification of errors in the ECIs of CEs. While we stick to the lattice gas convention in this work, the method is equally applicable to the Ising convention or, in principle, any linear model of sufficient complexity.

08 HYDROGEN↗
DOE OSTI2024

A Semi‐Automated, High‐Throughput Approach for the Synthesis and Identification of Highly Photo‐Cytotoxic Iridium Complexes

Abstract The discovery of new compounds with pharmacological properties is usually a lengthy, laborious and expensive process. Thus, there is increasing interest in developing workflows that allow for the rapid synthesis and evaluation of libraries of compounds with the aim of identifying leads for further drug development. Herein, we apply combinatorial synthesis to build a library of 90 iridium(III) complexes (81 of which are new) over two synthesise‐and‐test cycles, with the aim of identifying potential agents for photodynamic therapy. We demonstrate the power of this approach by identifying highly active complexes that are well‐tolerated in the dark but display very low nM phototoxicity against cancer cells. To build a detailed structure–activity relationship for this class of compounds we have used density functional theory (DFT) calculations to determine some key electronic parameters and study correlations with the experimental data. Finally, we present an optimised semi‐automated synthesise‐and‐test protocol to obtain multiplex data within 72 hours.

Kench, Timothy↗
DOE OSTI2024

A Semi‐Automated, High‐Throughput Approach for the Synthesis and Identification of Highly Photo‐Cytotoxic Iridium Complexes

Abstract The discovery of new compounds with pharmacological properties is usually a lengthy, laborious and expensive process. Thus, there is increasing interest in developing workflows that allow for the rapid synthesis and evaluation of libraries of compounds with the aim of identifying leads for further drug development. Herein, we apply combinatorial synthesis to build a library of 90 iridium(III) complexes (81 of which are new) over two synthesise‐and‐test cycles, with the aim of identifying potential agents for photodynamic therapy. We demonstrate the power of this approach by identifying highly active complexes that are well‐tolerated in the dark but display very low nM phototoxicity against cancer cells. To build a detailed structure–activity relationship for this class of compounds we have used density functional theory (DFT) calculations to determine some key electronic parameters and study correlations with the experimental data. Finally, we present an optimised semi‐automated synthesise‐and‐test protocol to obtain multiplex data within 72 hours.

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

Effect of dilute Rh on oxygen dissociation, spillover, and the oxidation of Cu across many orders of magnitude pressure

Knowledge of how trace amounts of more reactive metals influence the oxidation rate and mechanism of Cu surfaces is essential for developing strategies to optimize the performance of Cu-based catalysts. We find that the addition of 1% Rh to Cu(111) increases the initial O 2 dissociation rate by approximately 9-fold. CO poisoning experiments reveal that single Rh atoms activate O 2 and facilitate the spillover of atomic oxygen to Cu sites. Scanning tunneling microscopy (STM) and in situ X-ray photoelectron spectroscopy (XPS) support this mechanism, showing enhanced surface oxygen near Rh atoms. Here, a density functional theory (DFT)-based model demonstrates that Rh binds the O 2 precursor 0.15 eV more strongly than Cu(111) and lowers the O 2 dissociation barrier by 0.02 eV. Both single-crystal and nanoparticle experiments show that at low oxygen pressures, Rh enhances Cu oxidation, whereas at higher pressures, it inhibits deeper oxidation, as evidenced by in situ ultraviolet-visible (UV-vis) spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pyrite (001) Interface Chemistry is Controlled by a Sulfoxy Termination

Pyrite (FeS 2 ) is the most common sulfide mineral on Earth, forming through inorganic reactions in the crust and oceanic hydrothermal systems and via microbially driven processes in anaerobic sediments. The pyrite–water interface is the site of a wide range of adsorption and reaction processes in Earth systems including oxidation that dramatically affects the geochemistry of surface waters and influences global carbon and oxygen cycles. Mechanistic geochemical models of pyrite interfacial reactivity, however, are limited by the lack of experimentally derived atomistic structures of the reduced and reacting surfaces. Here, in this work, we reveal the atomic-scale structure of the pyrite (001)-water interface that forms at very low oxygen partial pressures, relevant to suboxic environments in Earth. The interface structure and surface speciation were obtained using the crystal truncation rod method supported by ambient-pressure photoelectron spectroscopy and density functional theory calculations. The surface is dominantly composed of disulfide groups bound to a single oxygen atom, forming a sulfoxy group that has no known molecular or bulk mineral analog. This surface is interpreted as the first step in the oxidative dissolution of pyrite. The sulfoxy group is readily protonated through surface acid–base reactions that alter the structure of interfacial water and the free energy of interfacial reactions. Surface iron sites are not oxidized. Surprisingly, this interface can likely develop in equilibrium with bulk pyrite in some reducing and acidic solutions. This termination is therefore likely representative of pyrite surfaces under a vast range of experimental, industrial and Earth conditions.

oxidation↗
DOE OSTI2024

Exceptional Electrical Detection of Trace NO 2 via Mixed Metal MOF-on-MOF Film-Based Sensors

The tunability of metal–organic frameworks (MOFs) makes them exceptional materials for the development of highly selective, low-power sensors for toxic gas detection. Herein, we demonstrate enhanced detection of NO 2 gas by a MOF-based electrical impedance sensor made using a unique mixed metal MOF-on-MOF synthesis. For this work, a combined experimental and computational study was performed using the exemplar Ni x Mg 1–x -MOF-74 to understand the fundamental structure–property relationships behind metal mixing and MOF film synthesis methods on sensor performance. Density functional theory results indicated that the presence of Ni in Mg-MOF-74 increased framework stability and increased the electron density of states at lower energies near the HOMO, as well as enhanced the NO 2 –Mg adsorption interaction. Impedance data of the Ni x Mg 1–x -MOF-74 films with larger Ni contents showed greater impedance change after exposure to 1 ppm of NO 2 gas. Furthermore, when synthesized through either a drop-cast or direct solvothermal film growth approach, the monometallic Ni-based sensors had the best performance. However, the mixed metal Ni x Mg 1–x -MOF-74 sensors synthesized through a MOF-on-MOF approach resulted in the highest impedance change, outperforming all monometallic Ni-based sensors. In particular, the mixed metal Ni-on-Mg-MOF-74 film was the best-performing sensor with an impedance change of 309 upon trace NO 2 exposure. Change in impedance response after NO 2 exposure was improved by 52% compared to the best monometallic Ni-on-Ni-MOF-74 sensor. Structural analysis of the Ni-on-Mg film showed that the first Mg-MOF-74 layer acts as a structural template controlling the structural features of the final film after metal exchange with Ni. This led to improved film quality, evidenced by the greater crystallinity and larger MOF grain sizes, and resulted in enhanced sensor performance which was not achievable through other metal mixing methods. Altogether, this study identifies structure–property relationships and synthetic templating methods that inform MOF-based sensor design, allowing for improved detection of toxic compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integrated design of aluminum-enriched high-entropy refractory B2 alloys with synergy of high strength and ductility

Refractory high-entropy alloys (RHEAs) are promising high-temperature structural materials. Their large compositional space poses great design challenges for phase control and high strength-ductility synergy. The present research pioneers using integrated high-throughput machine learning with Monte Carlo simulations supplemented by ab initio calculations to effectively navigate phase selection and mechanical property predictions, developing single-phase ordered B2 aluminum-enriched RHEAs (Al-RHEAs) demonstrating high strength and ductility. These Al-RHEAs achieve remarkable mechanical properties, including compressive yield strengths up to 1.7 gigapascals, fracture strains exceeding 50%, and notable high-temperature strength retention. They also demonstrate a tensile yield strength of 1.0 gigapascals with a ductility of 9%, albeit with B2 ordering. Furthermore, we identify valence electron count domains for alloy ductility and brittleness with the explanation from density functional theory and provide crucial insights into elemental influence on atomic ordering and mechanical performance. The work sets forth a strategic blueprint for high-throughput alloy design and reveals fundamental principles governing the mechanical properties of advanced structural alloys.

Science & Technology - Other Topics↗
DOE OSTI2020

Highly-Cyclable Room-Temperature Phosphorene Polymer Electrolyte Composites for Li Metal Batteries

Despite significant interest toward solid-state electrolytes owing to their superior safety in comparison to liquid-based electrolytes, sluggish ion diffusion and high interfacial resistance limit their application in durable and high-power density batteries. Here, a novel quasi-solid Li + ion conductive nanocomposite polymer electrolyte containing black phosphorous (BP) nanosheets is reported. The developed electrolyte is successfully cycled against Li metal (over 550 h cycling) at 1 mA cm -2 at room temperature. The cycling overpotential is dropped by 75% in comparison to BP-free polymer composite electrolyte indicating lower interfacial resistance at the electrode/electrolyte interfaces. Molecular dynamics simulations reveal that the coordination number of Li + ions around (trifluoromethanesulfonyl)imide (TFSI - ) pairs and ethylene-oxide chains decreases at the Li metal/electrolyte interface, which facilitates the Li + transport through the polymer host. Here, density functional theory calculations confirm that the adsorption of the LiTFSI molecules at the BP surface leads to the weakening of N and Li atomic bonding and enhances the dissociation of Li + ions. This work offers a new potential mechanism to tune the bulk and interfacial ionic conductivity of solid-state electrolytes that may lead to a new generation of lithium polymer batteries with high ionic conduction kinetics and stable long-life cycling.

2D materials↗
DOE OSTI2020

Passive NO x adsorption on Pd/H-ZSM-5: Experiments and modeling

A combined experimental and modeling study of the passive NO x adsorber (PNA) is presented that advances the understanding and prediction of the effects of various operating parameters and material properties of H-ZSM-5 and Pd/H-ZSM-5. Experiments reveal the NO x uptake process to be kinetically limited. Two microkinetic schemes are developed to explain the NO x uptake and release and NO 2 generation data. Both schemes are incorporated into a 1+1 D monolith model to predict and validate NO x uptake and temperature programmed desorption (TPD) data. Scheme I involves Z – [PdOH] + , Z – Pd 2+ Z – and Z – Pd + with Z – Pd + as the strongest NO binding site. Scheme II involves reduction of PdO 2 to PdO and nitrate formation. Model tuning utilizes a combination of infrared spectroscopy measurements and density functional theory (DFT) estimates of energy barriers. Lastly, the model is validated at different conditions including feed parameters and Pd loading.

42 ENGINEERING↗
DOE OSTI2025

Tailoring MoS 2 for Small-Molecule Electroreduction: The Role of Metal Doping and Heterostructures

The electrification of chemical transformations central to sustainable fuel production and waste valorization, such as overall water splitting (OWS), hydrogen evolution reaction (HER), and electrochemical reduction of CO 2 (CO 2 R), presents a powerful opportunity to advance carbon-neutral energy technologies. Transition metal dichalcogenides (TMDs), particularly MoS 2 , have emerged as promising electrocatalyst candidates, owing to their abundance, tunable active sites, and defect-rich structures. This review highlights recent progress in leveraging metal doping and heterostructure engineering of MoS 2 to enhance the electrocatalytic activity and selectivity. By compiling insights from experimental studies and density functional theory (DFT) predictions, we examine how defect creation, electronic structure modification, and interface design contribute to improved charge transport and catalytic efficiency. Particular emphasis is placed on rational design principles, synthetic strategies, and operando characterization methods that provide a pathway to understanding and optimizing MoS 2 -based materials. We also discuss the challenges of stability, mechanistic ambiguity, and scaling while outlining opportunities to bridge theory and experiment. Collectively, this review underscores how defect and heterostructure engineering of MoS 2 can accelerate the development of efficient, sustainable electrocatalysts for both fuel generation and waste-to-value generation.

CO2 reduction↗
DOE OSTI2026

Mechanistic insights into CO 2 capture and electrochemical conversion in nonaqueous Na–CO 2 batteries

Developing efficient energy storage systems that capture and convert CO 2 is critical for mitigating carbon emissions. Here, we report a Na–CO 2 battery with ruthenium dioxide (RuO 2 ) cathode catalysts and propane-1,3-diamine (PDA) as an electrolyte additive to enhance CO 2 capture and conversion efficiency. The integration of CO 2 adsorption and electrochemical reduction facilitates activation of the inert CO 2 molecule and circumvents gas–solid–liquid ternary-phase reactions at the interface. We employed density functional theory (DFT) calculations to systematically unravel the reaction mechanisms and energetics governing CO 2 reduction, both with and without PDA. Our results reveal an energetically favorable pathway toward the formation of Na 2 CO 3 and C as final discharge products, rather than sodium oxalate (Na 2 C 2 O 4 ). The CO 2 –amine adduct facilitates charge transfer from PDA to CO 2 , which results in activation of CO 2 . The kinetics of CO 2 conversion and regeneration of PDA were found to be significantly enhanced on the RuO 2 surface compared to the bulk electrolyte. More importantly, pre-activation of CO 2 via the amine–CO 2 adduct lowers the total overpotential to 2.44 V, compared to 3.13 V without PDA. This study provides fundamental insights into CO 2 electroreduction in Na–CO 2 batteries and underscores the promise of electrolyte engineering for sustainable CO 2 utilization and high-performance energy storage.

25 ENERGY STORAGE↗
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

Revealing the Formation Energy–Exfoliation Energy–Structure Correlation of MAB Phases Using Machine Learning and DFT

MAB phases are popular as high-temperature materials with high damage tolerance and excellent electrical conductivity that are used to exfoliate two-dimensional (2D) transition metal borides (MBenes); a promising material for developing next- generation nanodevices. In this report, we explore the correlation between the formation energy, exfoliation energy and structural factors of MAB phases using DFT and machine learning. In order to do so, we developed three different machine learning models based on the support vector machine, deep neural network, and random forest regressor to study the stability of the MAB phases by calculating their formation energies. Furthermore, our support vector machine and deep neural network models are capable of predicting the formation energies with mean absolute errors less than 0.1 eV/atom. We demonstrated that the stability of a MAB phase for a given transition metal decreases when the A element changes from Al to Tl. The density functional theory (DFT) re- vealed that M-A and M-B bond strength strongly correlates with the stability of MAB phases. In addition, the exfoliation possibility of 2D MBenes becomes higher when the A element changes from Al to Tl due to weakening of M-A and M-B bonds.

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

Collisional stopping power of ions in warm dense matter

Here, a model for the collisional stopping of ions on free electrons in warm dense matter is developed and explored. It is based on plasma kinetic theory, but with modifications to address the warm dense matter regime. Specifically, it uses the Boltzmann-Uehling-Uhlenbeck kinetic equation to incorporate effects of Fermi degeneracy of electrons. The cross section is computed from quantum scattering of electrons and ions occurring via the potential of mean force derived from an average atom model, which incorporates effects of strong Coulomb correlations. Predictions from this model show comparable accuracy to results from time-dependent density functional theory calculations for deuterium near solid density and a temperature of several electronvolts, at a fraction of the computational cost. Further, the model captures the transition of a plasma from the classical limit to the degenerate limit, including qualitative behaviors of solid state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Carbon-carbon interactions in warm dense titanium carbide

Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.

Ab initio molecular dynamics↗
DOE OSTI2021

Post-Newtonian kinetic theory

A kinetic theory for relativistic gases in the presence of gravitational fields is developed in the second post-Newtonian approximation. The corresponding Boltzmann equation is determined from the evolution of the one-particle distribution function with respect to the proper time along the world line of the particle. From the knowledge of the equilibrium Maxwell–Jüttner distribution function in the second post-Newtonian approximation the components of the particle four-flow and energy–momentum tensor are obtained. The Eulerian hydrodynamic equations for the mass density, mass-energy density and momentum density in the second post-Newtonian approximation are determined from the Boltzmann equation. It is shown that the combination of the hydrodynamic equations of mass and mass-energy densities leads to the hydrodynamic equation for the internal energy density in the first post-Newtonian approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗