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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Raman Digital Twin of Monolayer Janus Transition Metal Dichalcogenides

Monolayer transition metal dichalcogenides (TMDs) are a key class of two-dimensional (2D) materials with broad technological potential. Their Janus counterparts exhibit unique properties due to broken out-of-plane symmetry and further enrich the functionalities of TMDs. However, experimental synthesis and identification of Janus TMDs remain challenging. It is thus highly desirable to have a rapid, simple, and in situ characterization technique to monitor, in real time, the conversion process from the parent to Janus structure. Raman spectroscopy stands out for such a task as it is a powerful, nondestructive, and very commonly used tool to characterize 2D materials both in situ and ex situ. To realize the full potential of Raman spectroscopy on rapid characterization of Janus TMDs, we present a computational “Raman digital twin” library for various monolayer Janus TMDs in both 2H and Td phases. We focus on group-6 TMDs: MoS 2 , WS 2 , MoSe 2 , WSe 2 , MoTe 2 , WTe 2 and their Janus variants: MoSSe, MoSTe, MoSeTe, WSSe, WSTe, and WSeTe. Using first-principles density functional theory (DFT), we calculate their vibrational properties and predict distinct Raman fingerprints. These phonon and Raman signatures reflect each material’s structural symmetry and atomic composition, enabling clear identification via Raman spectroscopy. Our theoretical work supports experimental efforts by providing benchmarks for material identification, structural analysis, and quality control. In conclusion, the computational library expedites the discovery and development of Janus 2D materials, facilitating tighter integration between theoretical predictions and experimental validation.

Chalcogenides↗
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 OSTI2022

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Machine-learning enabled thermodynamic model for the design of new rare-earth compounds

We employ a descriptor based machine-learning approach to assess the effect of chemical alloying on formation-enthalpy of rare-earth intermetallics. Application of machine-learning approaches in rare-earth intermetallic design have been sparse due to limited availability of reliable datasets. In this work, we developed an ‘in-house’ rare-earth database with more than 600 + compounds, each entry was populated with formation enthalpy and related atomic features using high-throughput density-functional theory (DFT). The SISSO (sure independence screening and sparsifying operator) based machine-learning method with meaningful atomic features was used for training and testing the formation enthalpies of rare earth compounds. The complex lattice function coupled with the machine-learning model was used to explore the effect of transition metal alloying on the energy stability of Ce based cubic Laves phases (MgCu 2 type). The SISSO predictions show good agreement with high-fidelity DFT calculations and X-ray powder diffraction measurements. Our study provides quantitative guidance for compositional considerations within a machine-learning model and discovering new metastable materials. The electronic-structure of Ce-Fe-Cu based compound was also analyzed to get an in-depth understanding of the electronic origin of phase stability. The interpretable analytical models in combination with density-functional theory and experiments provide a fast and reliable design guide for discovering technologically useful materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the structure sensitivity of and CO coverage effects on formic acid decomposition on Pd surfaces

Using density functional theory calculations, the Pd-catalyzed vapor-phase formic acid decomposition was studied, with a focus on the structure sensitivity and CO coverage effects. A comprehensive reaction network was developed on both the (111) and (100) facets of Pd, at CO coverages of 0 and 5/9 monolayer (ML). Pd(100) was determined to be more reactive than Pd(111) at both CO coverages. The introduction of 5/9 ML CO decreased the activity of both facets significantly, due to destabilization of the surface intermediates and transition states on the CO-decorated surfaces. Three reaction pathways were explored on the clean surfaces: the formate (HCOO) pathway, the carboxyl (COOH) pathway leading to the formation of CO 2 , and the COOH pathway leading to the formation of CO (COOH→CO). Based on the DFT-derived energetics alone, it appears that all three pathways contribute to the reaction on clean Pd, whereas the presence of 5/9 ML of CO inhibits the HCOO pathway on both facets and favors the COOH→CO pathway on the (111) facet, but the COOH→CO 2 one on the (100) facet. Moreover, at high CO coverages, alternative spectator CO-assisted adsorbate decomposition pathways were discovered, which could potentially play a role in formic acid decomposition on Pd catalysts under realistic reaction conditions.

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

Theoretical and Experimental Insights into CO 2 Capture and Methanation over Amine-Grafted Ru-Based Catalysts

Carbon capture and storage (CCS) technologies, along with CO 2 capture and conversion methods, have emerged as crucial research areas to address rising CO 2 emissions. In this study, we seek to understand the mechanistic role of amines in enabling lower-energy pathways for CO 2 conversion. Our research focuses on the development and analysis of dual-functional materials (DFMs) engineered for the reactive capture and conversion (RCC) of CO 2 into methane, utilizing Ru catalysts grafted with amine groups. We employ Density Functional Theory (DFT) calculations using methylamine as a model amine to investigate the impact of amine groups on CO 2 methanation on a Ru(0001) surface, both in the presence and absence of amine groups. The amine ligand alters the carbon coordination environment, promoting direct C–O dissociation and potentially destabilizing the CO* adsorbate, thereby reducing the risk of CO poisoning. Additionally, we observe a preference for hydrogenation, although it becomes more energetically uphill in the amine-bound scenario. Our experiments, however, report similar CO 2 conversion and CH 4 production rates over the synthesized catalysts “Ru/TiO 2 ” and the amine (N-(2-aminoethyl)-3-aminoproplytrimethoxysilane (“diaminosilane”)) deposited catalyst “Diamine−Ru/TiO 2 ”. By constructing comparative reaction-free energy diagrams and performing microkinetic modeling (MKM) simulations, we link our theoretical findings with experimentally observed CO 2 uptake, conversion, and methane production rates. A microkinetic model was employed to investigate the anomaly, showing reduced amine–carbon complex coverage and increased CO 2 coverage at all temperatures. The MKM simulations consistently confirmed these trends. In conclusion, this comprehensive approach offers key insights into the role of the amine-CO 2 bond in methanation, highlighting a pathway toward lower-energy, more efficient CO 2 capture and conversion processes.

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

Machine-learning Kohn–Sham potential from dynamics in time-dependent Kohn–Sham systems

Abstract The construction of a better exchange-correlation potential in time-dependent density functional theory (TDDFT) can improve the accuracy of TDDFT calculations and provide more accurate predictions of the properties of many-electron systems. Here, we propose a machine learning method to develop the energy functional and the Kohn–Sham potential of a time-dependent Kohn–Sham (TDKS) system is proposed. The method is based on the dynamics of the Kohn–Sham system and does not require any data on the exact Kohn–Sham potential for training the model. We demonstrate the results of our method with a 1D harmonic oscillator example and a 1D two-electron example. We show that the machine-learned Kohn–Sham potential matches the exact Kohn–Sham potential in the absence of memory effect. Our method can still capture the dynamics of the Kohn–Sham system in the presence of memory effects. The machine learning method developed in this article provides insight into making better approximations of the energy functional and the Kohn–Sham potential in the TDKS system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Atomistic Mechanism Underlying the Si ( 111 ) - ( 7 × 7 ) Surface Reconstruction Revealed by Artificial Neural-Network Potential

The 7×7 reconstruction of the Si(111) surface represents arguably the most fascinating surface reconstruction so far observed in nature. Yet, the atomistic mechanism underpinning its formation remains unclear after it was discovered sixty years ago. Experimentally, it is observed post priori so that analysis of its formation mechanism can only be carried out in analogy with archaeology. Theoretically, density-functional theory (DFT) correctly predicts the Si(111)-(7×7) ground state but is impractical to simulate its formation process; while empirical potentials failed to produce it as the ground state. Developing an artificial neural-network potential of DFT quality, we carried out accurate large-scale simulations to unravel the formation of the Si(111)-(7×7) surface. We reveal a possible step-mediated atom-pop rate-limiting process that triggers massive nonconserved atomic rearrangements, most remarkably, a critical process of collective vacancy diffusion that mediates a sequence of selective dimer, corner-hole, stacking-fault, and dimer-line pattern formation, to fulfill the 7×7 reconstruction. Our findings may not only solve the long-standing mystery of this famous surface reconstruction but they also illustrate the power of machine learning in studying complex structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2022

From ground to excited electronic state dynamics of electron and ion irradiated graphene nanomaterials

Microscopes utilizing convergent electron and ion beams are emerging as powerful tools for both imaging and manipulating two-dimensional materials with atomic resolution, allowing the ultimate limits of nanofabrication to be realized. In this chapter, we detail the use of time-dependent electronic structure theory to determine the excited state properties and reactivity of functionalized graphene nanostructures. A time-dependent density functional theory treatment of electronic excitations of materials is presented, with specific emphasis on predicting the position-dependent electronic response of two-dimensional nanomaterials to electron beams. The method is demonstrated in a study highlighting the important role that electronic excitation can play in opening reaction pathways relevant to atomically precise defect manipulation in graphene. Finally, we provide some perspective on future development directions for methods of simulating nonequilibrium electronic and vibrational dynamics induced by electron/ion beams.

Lingerfelt, David↗
DOE OSTI2022

Artificial intelligence inferred microstructural properties from voltage–capacity curves

Abstract The quantification of microstructural properties to optimize battery design and performance, to maintain product quality, or to track the degradation of LIBs remains expensive and slow when performed through currently used characterization approaches. In this paper, a convolution neural network-based deep learning approach (CNN) is reported to infer electrode microstructural properties from the inexpensive, easy to measure cell voltage versus capacity data. The developed framework combines two CNN models to balance the bias and variance of the overall predictions. As an example application, the method was demonstrated against porous electrode theory-generated voltage versus capacity plots. For the graphite|LiMn $$_2$$ 2 O $$_4$$ 4 chemistry, each voltage curve was parameterized as a function of the cathode microstructure tortuosity and area density, delivering CNN predictions of Bruggeman’s exponent and shape factor with 0.97 $$R^2$$ R 2 score within 2 s each, enabling to distinguish between different types of particle morphologies, anisotropies, and particle alignments. The developed neural network model can readily accelerate the processing-properties-performance and degradation characteristics of the existing and emerging LIB chemistries.

25 ENERGY STORAGE↗
DOE OSTI2023

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sensing at the Nanoscale Using Nitrogen-Vacancy Centers in Diamond: A Model for a Quantum Pressure Sensor

The sensing of stress under harsh environmental conditions with high resolution has critical importance for a range of applications including earth’s subsurface scanning, geological CO2 storage monitoring, and mineral and resource recovery. Using a first-principles density functional theory (DFT) approach combined with the theoretical modelling of the low-energy Hamiltonian, here, we investigate a novel approach to detect unprecedented levels of pressure by taking advantage of the solid-state electronic spin of nitrogen-vacancy (NV) centers in diamond. We computationally explore the effect of strain on the defect band edges and band gaps by varying the lattice parameters of a diamond supercell hosting a single NV center. A low-energy Hamiltonian is developed that includes the effect of stress on the energy level of a ±1 spin manifold at the ground state. By quantifying the energy level shift and split, we predict pressure sensing of up to 0.3 MPa/Hz using the experimentally measured spin dephasing time. We show the superiority of the quantum sensing approach over traditional optical sensing techniques by discussing our results from DFT and theoretical modelling for the frequency shift per unit pressure. Importantly, we propose a quantum manometer that could be useful to measure earth’s subsurface vibrations as well as for pressure detection and monitoring in high-temperature superconductivity studies and in material sciences. Our results open avenues for the development of a sensing technology with high sensitivity and resolution under extreme pressure limits that potentially has a wider applicability than the existing pressure sensing technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerated Aging of a Polymer-Supported Molecular Amine Direct Air Capture Sorbent: Experiments and Theory

Amine-functionalized polymer sorbents are the preferred materials for CO2 extraction from air (DAC) due to their “high” CO2 adsorption capacity and selectivity at ultra-dilute conditions. The lifetime of adsorbents is a key property of equal importance as adsorption capacity, selectivity and kinetics with direct impact on the economics of commercial-scale operations. Therefore, a detailed investigation into the contributing factors (e.g. temperature, humidity) as well as the mechanism of amine degradation is vital for the development of amines with long-term stability. Herein, in situ attenuated total reflection (ATR) Fourier transform infrared spectroscopy (FTIR) is used to monitor the real time changes to amine-based sorbents under systematically varied conditions. In addition, IR spectra obtained from density functional theory molecular dynamics are used to obtain mechanistic insights pertaining to the degradation of amine groups in DAC sorbents.

Muldoon, Patrick↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

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

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2025

Computational and Experimental Investigation of Chiral and Achiral Two‐Dimensional Organic Lead Bromide Perovskites: Octahedral Distortions and Electronic and Optical Properties

A computational investigation is presented, in conjunction with synthesis and experimental characterization, into the structural, electronic, and optical properties of layered two-dimensional organic lead bromide perovskites. Materials based on the chiral (R/S)-4-fluoro-α-methylbenzylammonium (R/S-FMBA), which have been shown to lead to bright room-temperature circularly polarized luminescence, are contrasted with the similar achiral 4-fluorobenzylammonium (FBA). Using density functional theory (DFT) with van der Waals (vdW) corrections, relaxed structures (compared with X-ray diffraction, XRD) and optical absorption spectra (compared with experiments) are studied, as well as band structure and orbital character of transitions. A Python code is developed and provided to calculate octahedral distortions and compare DFT and XRD results, finding that vdW corrections are important for accuracy and that DFT overestimates octahedral tilt angles. (FMBA) 2 PbBr 4 shows among the largest tilt angle differences (often termed Δ β ) reported, 14°–15°, indicating strong inversion symmetry-breaking, which enables its chiral emission. A large resulting Dresselhaus spin-splitting effect is found. The lowest-energy optical transitions involve the perovskite only and are polarized within the layer. This work furthers understanding of structure-property relations with applications to optoelectronics and spintronics.

UV/vis spectroscopy↗
DOE OSTI2024

Molecular Fe-N 4 Moieties Coupled with Atomic Co-N 4 Sites Toward Improved Oxygen Reduction Performance

Research on high-efficiency and cost-efficient catalysts for oxygen reduction reaction (ORR) is still a vital but challenging issue for commercializing metal–air batteries. Herein, a single-molecule/atom hybrid catalyst is developed to boost the ORR, in which iron phthalocyanine molecules containing molecular Fe-N 4 moieties couple with atomic Co-N 4 sites on the surface of polyhedral carbon. Density functional theory calculations reveal that face-to-face laminated construction of Fe-N 4 and Co-N 4 in the hybrid catalyst can effectively modulate the electronic structure of active iron atoms and reduce the energy barrier of the rate-determining step for ORR. As a result, this hybrid catalyst demonstrates excellent ORR performance, featuring a half-wave potential of 0.904 V, a peak power density of 238.3 mW cm -2 for zinc–air battery, and outstanding electrocatalytic stability. Here, this work offers a distinctive and robust molecular/atomic engineering approach to creating efficient electrocatalysts, advancing the fields of metal–air batteries.

36 MATERIALS SCIENCE↗