Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT simulations”

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

Softening of dd excitation in the resonant inelastic x-ray scattering spectra as a signature of Hund's coupling in nickelates

We investigate the effects of Hund's coupling on the resonant x-ray absorption spectra of the recently discovered family of layered nickelate superconductors. We contrast two scenarios depending on the relative strength of the ratio of the effective Hund's coupling ( J H ) to the crystal fields ( Δ ) in these systems. We carry out the cluster and DFT + DMFT simulations of the RIXS signal at the Ni L edge for different values of Hund's coupling. We find the latter dominates for the parent compound while the former becomes important for sufficiently large doping. Our results are consistent with the observations of a softening of a RIXS peak as a function of doping by Rossi [], only when the Hund coupling is sizable. To interpret the results, we separate the theoretical RIXS signal into spin conserving and nonspin conserving channels and conclude that the infinite layer nickelates are in a regime where Δ and J compete effectively and suggest further experimental tests of the theory. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Thermally Stable Co@C3N4 Single-Atom Catalysts for CO Oxidation: Atomic-Level Insights into Structure and Activity

Single Co atoms supported on C3N4 (Co@C3N4) have demonstrated high activity and selectivity in photocatalysis. However, the investigation of structure–function relationships and reaction mechanisms under photocatalytic conditions is very challenging due to the complex conditions of light absorption, charge transfer, and catalysis. In this study, we employed thermal CO oxidation as a prototypical probe reaction to benchmark the intrinsic catalytic performance and track the active-site evolution of Co@C3N4. Single Co atoms were identified and shown to be the catalytically active sites for CO oxidation based on control experiments and isotope-labeling experiments. The Co sites remained atomically dispersed before, during, and after the reaction with temperatures up to 400 °C, as established by in situ X-ray absorption fine structure (XAFS) combined with density functional theory (DFT), FDMNES simulations, and dynamic-time-warping (DTW)-assisted X-ray absorption near edge structure (XANES) matching. Together with theoretical calculations, the integrated analysis reveals a stable coordination environment under reaction conditions, which correlates with sustained activity, establishing Co@C3N4 single-atom catalysts as thermally stable CO oxidation catalysts. Beyond these findings, the current study provides a workflow for unambiguously assigning active sites in Co@C3N4 for thermal CO oxidation. This workflow will aid the understanding of their behavior in photocatalysis in the future, where light-driven dynamics obscure direct structure–function links. Notably, this study provides fundamental insights for the rational design of robust single-atom catalysts and a foundation for the broader application of Co@C3N4 catalysts in oxidation reactions.

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

Degradation Mechanism of Phosphate‐Based Li‐NASICON Conductors in Alkaline Environment

NASICON-type Li conductors (Li-NASICON) have traditionally been regarded as promising candidates for solid-state Li-air battery applications because of their stability in water and ambient air. However, the presence of water in the cathode of a Li-air battery can induce a highly alkaline environment by modifying the discharge product from Li 2 O 2 to LiOH which can potentially degrade cathode and separator materials. This study investigates the alkaline stability of common Li-NASICON chemistries through a systematic experimental study of LiTi x Ge 2-x (PO 4 ) 3 (LTGP) with varying x = 0–2.0. Density functional theory calculations are combined to gain a mechanistic understanding of the alkaline instability. It is demonstrated that the instability of LTGP in an alkaline environment is mainly driven by the dissolution of PO 4 3– groups, which subsequently precipitate as Li 3 PO 4 . The introduction of Ti facilitates the formation of a Ti-rich compound on the surface that eventually passivates the material, but only after significant bulk degradation. Consequently, phosphate-based Li-NASICON materials exhibit limited alkaline stability, raising concerns about their viability in humid Li-air batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Synthesis of Gold Hydride at High Pressure and High Temperature

Gold is an unreactive metal and its chemical interactions with hydrogen have only recently been explored. Here, in this study, we report the formation of gold hydride above 40 GPa and 2200 K in X-ray free electron laser heated diamond anvil cells using various hydrocarbons as hydrogen sources. Above 40 GPa, a hexagonal phase emerges close to the gold melting point, corresponding to a hydride with stoichiometry Au 2 H x , with x increasing from 0 to near 1 with pressure from 40 to 80 GPa. This is a high-temperature phase which reverts to face centered cubic gold on cooling to 295 K. Accompanying DFT-MD simulations are in excellent agreement with experiment and reveal the structure to consist of an hexagonal close packed gold lattice with atomic hydrogen disordered in the interstices. The hydrogen is superionic and exhibits high diffusivity through the crystalline gold lattice. Our results present the first solid-state binary compound of gold and hydrogen.

Ab initio calculations↗
DOE OSTI2025

Microstructure‐Dependent Sodium Storage Mechanisms in Hard Carbon Anodes

Sustainable energy storage is essential to support the transition to renewables and meet the increasing demand for energy. Sodium‐ion batteries (NIBs) are attractive for grid‐scale energy storage due to the abundance and low cost of sodium, sustainability of other battery components, and electrochemical performance. Hard carbon (HC) is a leading anode material for NIBs, but its complex microstructure complicates the understanding of sodium storage mechanisms. Using X‐ray total scattering and density functional theory calculations, this study clarifies how HC's microstructural variations influence sodium storage across the slope (high potential) and plateau (low potential) regions of the potential capacity curve. In the slope region, sodium initially adsorbs at high‐binding energy defect sites and subsequently intercalates between graphene layers, adsorbing at low‐binding energy defect sites, correlating with different slopes observed during initial sodiation. Initial irreversibility arises from sodium trapping at surface defects and solid electrolyte interface formation. In the plateau region, sodium simultaneously intercalates and fills pores, influenced by pore size, interlayer spacing, and defect concentration. HCs with larger pore sizes form larger sodium clusters. In conclusion, the proposed mechanism underscores the role of microstructure engineering in enhancing HC performance and advancing NIBs for grid‐scale energy storage.

36 MATERIALS SCIENCE↗
DOE OSTI2023

First-principles study of interfaces in Al/SiC metal-matrix composite system

In this work, first-principles calculations were performed on the interfaces between Al and SiC, which is a widely used strengthening agent in aluminum metal-matrix-composites (Al/MMC). C-terminated interfaces have much larger work of adhesion than Si-terminated interfaces, indicating that the former has much stronger interfacial bonding. The electron localization function shows that the chemical bonding between Al and C has a strong covalent character, while the bonding between Al and Si is largely metallic. As a result of the vastly different chemical bonding, the work of adhesion for C-terminated interfaces increases with the number of dangling bonds at the interface, while the opposite trend was observed for Si-terminated interfaces. Additionally, the interface energy for Si-terminated interfaces is comparable to that for C-terminated interfaces, suggesting both types of terminations can coexist in the Al/SiC system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2021

Net Electronic Charge as an Effective Electronic Descriptor for Oxygen Release and Transport Properties of SrFeO 3 -based Oxygen Sorbents

Perovskite oxides, as oxygen sorbents, exhibit excellent potential in thermochemical redox applications such as chemical looping air separation (CLAS), resulting from their excellent redox properties and high tunability. The structural and compositional flexibility of perovskites, while making these mixed oxides highly versatile, also poses challenges for their design and optimization. Moreover, the elevated operating temperature and dynamic oxygen partial pressure swings during redox processes add additional complications for oxygen sorbent design. From a fundamental standpoint, finding simple yet effective descriptors for the effect of cation substitutions on oxygen sorbents’ redox and oxygen transport properties are of significant importance for rational optimization of perovskite-based oxygen sorbents. In the present study, a series of A/B-site-doped SrFeO 3-δ perovskites are investigated using density functional theory (DFT) + U simulations. From these simulations, an effective electronic structural descriptor, the net electronic charge (Δe) of oxygen anions, is extracted based on its excellent correlations with vacancy formation energy, experimentally derived oxygen capacity, and oxygen vacancy migration barrier. Overall, these findings provide an effective tool to correlate the oxygen sorbents’ redox performances with a simple yet physically meaningful descriptor for rational design and optimization of perovskites in the context of CLAS.

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

Tackling the Capacity Fading Issue of Li–S Battery by a Functional Additive—Hexafluorobenzene

Due to the strong electron-withdrawing fluorine groups, the aromatic ring of hexafluorobenzene (HFB) becomes more electron-deficient and strongly interacts with the lone-pair electrons of the medium. When used as an electrolyte additive for Li-S cell, HFB preferably coordinates with low-order lithium sulfides Li 2 S and Li 2 S 2 and subsequently reacts and converts both into active aromatic sulfides-bis(pentafluorophenyl) disulfides and bis-(pentafluorophenyl) polysulfides in situ-thus facilitating the utilization of the electronically insulating Li 2 S and Li 2 S 2 and improving the cycling stability of the Li-S cell. Finally, this research opens a new avenue to tackle the inherent issues associated with the Li-S chemistry.

25 ENERGY STORAGE↗
DOE OSTI2020

Probing Heterogeneous Charge Distributions at the α-Al 2 O 3 (0001)/H 2 O Interface

Unlike metal or semiconductor electrodes, the surface charge resulting from the protonation or deprotonation of insulating mineral oxides is highly localized and heterogeneous in nature. In this work the Stark active C≡N stretch of potassium thiocyanate is used as a molecular probe of the heterogeneity of the interfacial electrostatic potential at the α-Al 2 O 3 (0001)/H 2 O interface. Vibrational sum frequency generation (vSFG) measurements performed in the OH stretching region suggest that thiocyanate species organize interfacial water similarly to halide ions. Changes in the electrostatic potential are then tracked via Stark shifts of the vibrational frequency of the thiocyanate stretch. Our vSFG measurements show that we can simultaneously measure the vSFG response of SCN – ions experiencing charged and neutral surface sites. We assign local potentials of +308 and -154 mV to positively and negatively charged aluminol groups that are present at pH = 4 and pH = 10, respectively. Thiocyanate anions at positively charged surface sites and negatively charged surface sites and those participating in contact ion pairing adopt similar orientations and are oppositely oriented relative to thiocyanate ions near neutral surface sites. All four species followed Langmuir adsorption isotherms. Density functional theory–molecular dynamics (DFT-MD) simulations of SCN – near the neutral α-Al 2 O 3 (0001)/H 2 O interface show that the vSFG response in the C≡N stretch region originates from a SCN–H–O–Al complex, suggesting the surface site specificity of these experiments. To our knowledge this is the first spectroscopic measurement of local potentials associated with a heterogeneously charged surface. The ability to probe the evolution of local charges in situ could provide vital insight into many industrial, electrochemical, and geochemically relevant interfaces.

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

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study

Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.

Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗
DOE OSTI2024

Experimental and Computational Study of Pyrogenic Carbonaceous Matter Facilitated Hydrolysis of 2,4-Dinitroanisole (DNAN)

This study investigated the reaction pathway of 2,4-dinitroanisole (DNAN) on the pyrogenic carbonaceous matter (PCM) to assess the scope and mechanism of PCM-facilitated surface hydrolysis. DNAN degradation was observed at pH 11.5 and 25 °C with a model PCM, graphite, whereas no significant decay occurred without graphite. Experiments were performed at pH 11.5 due to the lack of DNAN decay at pH below 11.0, which was consistent with previous studies. Graphite exhibited a 1.78-fold enhancement toward DNAN decay at 65 °C and pH 11.5 relative to homogeneous solution by lowering the activation energy for DNAN hydrolysis by 54.3 ± 3.9%. This is supported by our results from the computational modeling using Car–Parrinello simulations by ab initio molecular dynamics/molecular mechanics (AIMD/MM) and DFT free energy simulations, which suggest that PCM effectively lowered the reaction barriers by approximately 8 kcal mol –1 compared to a homogeneous solution. Quaternary ammonium (QA)-modified activated carbon performed the best among several PCMs by reducing DNAN half-life from 185 to 2.5 days at pH 11.5 and 25 °C while maintaining its reactivity over 10 consecutive additions of DNAN. We propose that PCM can affect the thermodynamics and kinetics of hydrolysis reactions by confining the reaction species near PCM surfaces, thus making them less accessible to solvent molecules and creating an environment with a weaker dielectric constant that favors nucleophilic substitution reactions. Nitrite formation during DNAN decay confirmed a denitration pathway, whereas demethylation, the preferred pathway in homogeneous solution, produces 2,4-dinitrophenol (DNP). Denitration catalyzed by PCM is advantageous to demethylation because nitrite is less toxic than DNAN and DNP. These findings provide critical insights for reactive adsorbent design that has broad implications for catalyst design and pollutant abatement.

2,4-dinitroanisole (DNAN)↗
DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature

It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. Here, we find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.

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

Multiscale computational understanding and growth of 2D materials: a review

Abstract The successful discovery and isolation of graphene in 2004, and the subsequent synthesis of layered semiconductors and heterostructures beyond graphene have led to the exploding field of two-dimensional (2D) materials that explore their growth, new atomic-scale physics, and potential device applications. This review aims to provide an overview of theoretical, computational, and machine learning methods and tools at multiple length and time scales, and discuss how they can be utilized to assist/guide the design and synthesis of 2D materials beyond graphene. We focus on three methods at different length and time scales as follows: (i) nanoscale atomistic simulations including density functional theory (DFT) calculations and molecular dynamics simulations employing empirical and reactive interatomic potentials; (ii) mesoscale methods such as phase-field method; and (iii) macroscale continuum approaches by coupling thermal and chemical transport equations. We discuss how machine learning can be combined with computation and experiments to understand the correlations between structures and properties of 2D materials, and to guide the discovery of new 2D materials. We will also provide an outlook for the applications of computational approaches to 2D materials synthesis and growth in general.

Momeni, Kasra↗
DOE OSTI2020

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

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

Candidate structure for the H 2 -PRE phase of solid hydrogen

Experimental progress finally reached the metallic solid hydrogen phase, which was predicted by Wigner and Huntington over 80 years ago. However, the different structures in the phase diagram are still being debated due to the difficulty of diffraction experiments for high-pressured hydrogen. The determination of crystal structures under extreme condition is both of the basic condensed matter physics, and in planetary science: the behavior of giant gaseous planets (e.g., Jupiter and Saturn) strongly depends on the properties of inner high-pressured hydrogen. This work describes possible structures appearing under high pressures of 400~600 GPa. We applied a structural search using particle swarm optimization with density functional theory (DFT) to propose several candidate structures. For these structures, we performed fixed-node diffusion Monte Carlo simulations combined with DFT zero-point energy corrections to confirm their relative stability. We found P2 1 /c-8 as a promising candidate structure for the H 2 -PRE phase. P2 1 /c-8 is predicted the most stable at 400 and 500 GPa. P2 1 /c-8 reproduces qualitatively the IR spectrum peaks observed in the H 2 -PRE phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗