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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Additive Lithography–Organic Monolayer Patterning Coupled with an Area-Selective Deposition

The combination of area-selective deposition (ASD) with a patternable organic monolayer provides a versatile additive lithography platform, enabling the generation of a variety of nanoscale feature geometries. Stearate hydroxamic acid self-assembled monolayers (SAMs) were patterned with extreme ultraviolet (λ = 13.5 nm) or electron beam irradiation and developed with ASD to achieve line space patterns as small as 50 nm. Density functional theory was employed to aid in the synthesis of hydroxamic acid derivatives with optimized packing density to enhance the imaging contrast and improve dose sensitivity. Near-edge X-ray absorption fine structure spectroscopy and infrared spectroscopy reveal that the imaging mechanism is based on improved deposition inhibition provided by the cross-linking of the SAM to produce a more effective barrier during a subsequent deposition step. With patterned substrates composed of coplanar copper lines and silicon spacers, hydroxamic acids selectively formed monolayers on the metal portions and could undergo a pattern-wise exposure followed by ASD in the first combination of a patternable monolayer with ASD. This material system presents an additional capability compared to traditional ASD approaches that generally reflect a starting patterned surface. Finally, this bottoms-up additive approach to lithography may be a viable alternative to subtractive nanoscale feature generation.

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

Constructing a square-like copper cluster to boost C–C coupling for CO 2 electroreduction to ethylene

The CO 2 electroreduction reaction (CO 2 ER) to ethylene (C 2 H 4 ) offers the dual promise of lowering CO 2 emission while storing energy from renewable electricity, for which the development of highly efficient electrocatalysts is of great significance. Herein, by means of density functional theory (DFT) computations, we designed an electrocatalyst for CO 2 -to-C 2 H 4 conversion by anchoring a Cu 5 cluster supported on a MoS 2 monolayer with an S monovacancy (Cu 5 @MoS 2 ). Our results revealed that one Cu atom of the Cu 5 cluster was embedded into the framework of the defective MoS 2 monolayer, while the other four Cu atoms form a square-like island over the substrate surface. Interestingly, the C–C coupling between two *CO species can easily occur on the unique square-like active site with a low kinetic barrier of 0.56 eV to form the key *C 2 O 2 intermediate, which can then be hydrogenated to the C 2 H 4 product with a very low limiting potential (–0.32 eV). Significantly, alkaline conditions (pH = 13) are beneficial to further promote C 2 H 4 synthesis. Finally, our work may offer a new avenue to precisely modulate the structures of Cu clusters for converting CO 2 into high-value target products.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Thermal conductivity of α -U with point defects

Here, we develop a theoretical model for thermal conductivity of α-U that combines density functional theory calculations and the coupled electron–phonon Boltzmann transport equation. The model incorporates both electron and phonon contributions to thermal conductivity and achieves good agreement with experimental data over a wide temperature range. The dominant scattering mechanism governing thermal transport in α-U at different temperatures is examined. By including phonon–defect and electron–defect scatterings in the model, we study the effect of point defects including U-vacancy, U-interstitial, and Zr-substitution on the thermal conductivity of α-U. The degradation of anisotropic thermal conductivity due to point defects as a function of defect concentration, defect type, and temperature is reported. This model provides insights into the impact of defects on both phonon and electron thermal transport. It will promote the fundamental understanding of thermal transport in α-U and provide a ground for investigation of coupled electron–phonon transport in metallic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2024

Intricate short-range order in GeSn alloys revealed by atomistic simulations with highly accurate and efficient machine-learning potentials

GeSn alloys hold promise for silicon-compatible integrated applications in electronics, photonics, and topological quantum devices. However, understanding their intricate structures using density functional theory (DFT) calculations is hindered by spatiotemporal constraints. To overcome this limitation, we develop highly accurate and efficient machine-learning interatomic potentials based on a neuroevolution potential approach with farthest point sampling on a comprehensive DFT data set. The application of the developed machine-learning potential in large-scale atomistic simulations bridges the spatiotemporal gap between modeling and advanced characterizations, and facilitates the discovery of structural intricacies in GeSn alloys. Through extensive statistical sampling, we identify a type of short-range order (SRO) that is distinguished by both its structural signature and electronic band gap from the SRO structure previously predicted. Modeling based on a large simulation cell reveals the coexistence of nano SRO domains with various degrees of ordering, demonstrating a complex spatial heterogeneity of SRO structure. Finally, our study not only reinforces the significance of fine-level structural information in alloys, but it also constitutes an effective framework for exploring SRO in a broad range of complex alloys based on highly accurate and effective machine-learning potentials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Charge Transfer Excitation with Variation After Response Quantum Monte Carlo (Final Technical Report)

The electronically excited states of molecules play critical roles both in spectroscopy and in technological application, ranging from UV/Vis absorption experiments to the production of solar fuels. However, certain types of excitations, such as charge transfer excitations, remain difficult to treat with existing methods such as time dependent density functional theory or linear response coupled cluster. Under this Early Career Research Award, we developed a number of different theoretical methods to address this issue, which began with Variational Monte Carlo methods and branched into multiple types of excited-state-specific quantum chemistry.

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

Developing new and understanding old approximations in TDDFT

When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange–correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn–Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange–correlation potential from a coupling-constant path integration.

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

Development of local hybrid density functionals to treat self-interaction error and many-electron effects

The goal of this project was to improve density functional theory (DFT) for systems where conventional semilocal approximations are limited by self-interaction error (SIE) and by near-degeneracy or strong many-electron effects. While DFT remains the only broadly practical first-principles framework for large-scale materials simulations, its predictive accuracy is often challenged in situations involving stretched bonds, charge transfer, transition-metal chemistry, magnetic couplings, band gaps, and correlated electronic states. This project addressed these limitations by developing physically grounded and numerically robust exchange–correlation functionals that retain the efficiency of modern DFT while extending its predictive scope.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

Recent Developments in DFTB+, a Software Package for Efficient Atomistic Quantum Mechanical Simulations

DFTB+ is a flexible, open-source software package developed by its community, designed for fast and efficient atomistic quantum mechanical simulations. It employs various methods that approximate density functional theory (DFT), such as density functional-based tight binding (DFTB) and the extended tight binding (xTB) approach allowing simulations of large systems over extended time scales with reasonable accuracy, while being significantly faster than traditional ab initio methods. In recent years, several new extensions of the DFTB method have been developed and implemented in the DFTB+ program package in order to improve the accuracy and generality of the available simulation results. In this paper, we review those enhancements, show several use case examples and discuss the strengths and limitations of its features.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2025

The total neutron cross section of liquid and solid ammonia

Ammonia is a material of interest for future neutron moderators at high-power sources due to its high hydrogen density, low melting point, and resistance to polymerization in an intense radiation field. Its performance in such applications cannot currently be calculated due to the absence of suitable computer models for the interaction of neutrons with ammonia under relevant conditions. In an effort to develop suitable scattering kernels for computer simulations of moderator performance, we have conducted a series of Density Functional Theory and Molecular Dynamics calculations of the molecular-level thermal properties of ammonia at various temperatures within both the solid and liquid phases. In this paper, we compare computer calculations for the energy-dependent total neutron cross section of ammonia, based on these models, to experimental measurements of those cross sections at temperatures of 221 K, 180 K, and 35 K. The experimental data were collected over an energy range from 0.1 meV to 10 eV using time-of-flight techniques at the Low Energy Neutron Source (LENS) facility at Indiana University. This comparison provides a first validation in the development of thermal scattering libraries for Monte Carlo source design simulations based on liquid and solid ammonia. In conclusion, we also provide some insights into where additional development of tools for creating such models may be needed.

Ammonia↗
DOE OSTI2022

Insights into the Hydrogen Evolution Reaction on 2D Transition-Metal Dichalcogenides

Understanding hydrogen evolution reaction (HER) behaviors over two-dimensional transition-metal dichalcogenides (2D-TMDs) is critical for the development of nonprecious HER electrocatalysts with better activity. Here, in this work, by combining density functional theory calculations with microkinetic modeling, we thoroughly investigated the HER mechanism on 2D-TMDs. We find an important dependence of simulated cell size on the calculated hydrogen adsorption energy and the activation barrier for MoS 2 . Distinct from previous “H migration” mechanisms proposed for the Heyrovsky reaction, the rate-determining step for MoS 2 , we propose that the Mo site only serves as the stabilized transition state rather than H adsorption. In comparison to transition-metal electrocatalysts, we find that the activation barrier of the Heyrovsky reaction on 2D-TMDs scales with the hydrogen adsorption energy exactly as for transition metals except that all activation energies are displaced upward by ca. 0.4 eV. This higher Heyrovsky activation barrier is responsible for the substantially lower activity of 2D-TMDs. We further show that this higher activation barrier stems from the more positively charged adsorbed hydrogen on the chalcogenides interacting repulsively with the incoming proton. Based on these insights, we discuss potential strategies for the design of nonprecious HER catalysts with activity comparable to Pt.

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

High-Performance Iridium–Molybdenum Oxide Electrocatalysts for Water Oxidation in Acid: Bayesian Optimization Discovery and Experimental Testing

Ir oxides are costly and scarce catalysts for oxygen evolution reaction (OER) in acid. There has been extensive interest in developing alternatives that are either Ir-free or require smaller amounts of Ir to drive the reactions at acceptable rates. One design strategy is to identify Ir-based mixed oxides that achieve similar performance while requiring smaller amounts of Ir. The obstacle to this strategy has been a very large phase space of the Ir-based mixed metal oxides, in terms of the metals combined with Ir and the different crystallographic structures of the mixed oxides, which prevents a thorough exploration of possible materials. In this work, we developed a workflow that uses machine-learning-aided Bayesian optimization in combination with density functional theory to make the exploration of this phase space plausible. This screening identified Mo as a promising dopant for forming acid-tolerant Ir-based oxides for the OER. We synthesized and characterized the Ir–Mo mixed oxides in the form of thin-film electrocatalysts with a known surface area. We show that these mixed oxides exhibited overpotentials ~30 mV lower than a pure Ir control while maintaining 24% lower Ir dissolution rates than the Ir control. Furthermore, these findings suggest that Mo is a promising dopant and highlight the promise of machine learning to guide the experimental exploration and optimization of catalytic materials.

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

Uncovering the True Active Sites in Ni–N–C Catalysts for CO 2 Electroreduction

Understanding and designing active sites in single-atom catalysts (SACs) requires going beyond static models to capture their dynamic evolution under realistic electrochemical conditions. Here, in this work, we develop an integrated theoretical framework that accounts for operational conditions, by combining grand canonical density functional theory (GC-DFT) with machine-learning-accelerated sampling, to uncover structure–activity–stability relationships in Ni–N–C SACs for the CO 2 reduction reaction (CO 2 RR). A library of NiN x C 4–x (x = 0–4) motifs─representing coordination defects likely formed during high-temperature synthesis─was systematically evaluated. Under working conditions, these sites were found to undergo hydrogenation, and NiN 3 C 1_ H 1 was identified as the most probable active site. At reducing potentials, hydrogen adsorbs spontaneously at C–Ni bridge sites rather than Ni top sites, while subsurface hydrogen facilitates bent CO 2 adsorption crucial for activation. High CO 2 RR selectivity toward CO arises from site separation: Ni centers drive CO2RR, while the hydrogen evolution reaction (HER) occurs at the C–Ni bridge or N sites and from thermodynamic suppression of HER at moderate hydrogen coverage. At more negative potentials, a shift in the CO 2 RR rate-determining process (RDP) and Ni out-of-surface displacement induced by coadsorption of H and H 2 O jointly reduce activity and selectivity. Thus, both the high CO2RR selectivity of Ni–N–C catalysts and its reversal with more negative potentials can be rationalized by accounting for hydrogenated surfaces. This highlights the necessity of modeling realistic; in situ conditions. This framework provides generalizable insights into the dynamic behavior of active sites in SACs, offering guidance for the rational design of active and robust catalysts for a wide range of electrochemical reactions.

25 ENERGY STORAGE↗
DOE OSTI2023

Kramers nodal lines and Weyl fermions in SmAlSi

Kramers nodal lines (KNLs) have recently been proposed theoretically as a special type of Weyl line degeneracy connecting time-reversal invariant momenta. KNLs are robust to spin orbit coupling and are inherent to all non-centrosymmetric achiral crystal structures, leading to unusual spin, magneto-electric, and optical properties. However, their existence in in real quantum materials has not been experimentally established. Here we gather the experimental evidence pointing at the presence of KNLs in SmAlSi, a non-centrosymmetric metal that develops incommensurate spin density wave order at low temperature. Using angle-resolved photoemission spectroscopy, density functional theory calculations, and magneto-transport methods, we provide evidence suggesting the presence of KNLs, together with observing Weyl fermions under the broken inversion symmetry in the paramagnetic phase of SmAlSi. We discuss the nesting possibilities regarding the emergent magnetic orders in SmAlSi. Our results provide a solid basis of experimental observations for exploring correlated topology in SmAlSi.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular dynamics study on the interfacial properties of mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders on graphene and graphite surfaces

This study investigates the behavior of two different mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders, polyvinylpyrrolidone:polyvinylidene difluoride (PVP:PVDF) and polyvinylpyrrolidone:polyacrylic acid (PVP:PAA), at graphene and graphite interfaces using classical molecular dynamics simulations. The aim is to identify the best performing monomer binder blend and carbon-based material for the design of battery-optimized energy devices. The PVP:PAA monomer binder blend and graphite are found to have the best interaction energies, densification upon adsorption, and more ordered structure. The adsorption of both monomer binder blends is strongly guided by the higher affinity of PVP and PAA monomeric molecules for the surfaces compared to PVDF. The structure of adsorbed layers of PVP:PVDF monomer binder blend on graphene and graphite develops more quickly than PVP:PAA, indicating faster kinetics. This study complements a previous density functional theory study recently reported by our group and contributes to a better understanding of the nanoscopic features of relevant interfacial regions involving mixtures of monomers of PVP-based battery binders and different carbon-based materials. In conclusion, the effect of a blend of commonly used monomer binders on carbon-based materials is essential for obtaining tightly bound anode and cathode active materials in lithium-ion batteries, which is crucial for designing battery-optimized energy devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗