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

Understanding Solid-Electrolyte Interphase Formation at the Lithium Metal Anode of Lithium-Sulfur Batteries

Lithium-sulfur (Li-S) batteries are one of the most promising chemistries for the next generation of energy storage devices. In a standard Li-S cell, the metallic lithium anode is an essential component of the battery due to its low density, extremely high theoretical specific capacity, and very low negative electrochemical potential. However, several challenges related to the use of Li metal have prevented this battery technology from becoming commercially available. For instance, Li-anodes are highly reactive, which results in the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. This, combined with the non-uniform deposition of Liions during plating at the anode surface, can seriously affect the performance, cycling, and safety of the battery. In addition, it has been suggested that a controlled SEI formation at the metallic Li anode can yield enhanced battery performance. Hence, a comprehensive understanding of how the SEI layer is formed can help elucidate improvements in this battery technology. In this work, the formation of the SEI layer in Li-S batteries is investigated using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. The stability of liquid organic-based electrolyte components, salt concentration, electron-rich environments, and the use of inorganic solid-state electrolyte (SSE) materials are explored in order to provide molecular-level fundamental insights into how the nature and composition of the electrolyte can alter the initial stages of the SEI formation. The products of electrolyte decomposition and electro-/chemical reaction pathways of liquid electrolyte component such as solvent and salts are predicted from AIMD simulations of electrolyte-anode interfaces. Energetics of reactions from DFT ground-state optimizations are also presented to confirm the decomposition mechanisms. Some significant differences are then drawn regarding the use of low and high salt concentrations. Simulations under electron-rich environments show additional multielectron electrochemical reactions of solvent and salt decomposition taking place due to the excess of electrons and the presence of radical anions in the solution. Finally, the initial stages of formation, stability, and main constituents of interfaces between Li-metal and Sbased electrodes with sulfide-based SSEs are characterized in detail.

25 ENERGY STORAGE↗
DOE OSTI2021

Effect of excess Mg to control corrosion in molten MgCl 2 and KCl eutectic salt mixture

Structural alloys may experience corrosion when exposed to molten chloride salts due to selective dissolution of active alloying elements. One way to prevent this is to make the molten salt reducing. For the KCl + MgCl 2 eutectic salt mixture, pure Mg can be added to achieve this. However, Mg can form intermetallic compounds with nickel at high temperatures, which may cause alloy embrittlement. This work shows that an optimum level of excess Mg could be added to the molten salt which will prevent corrosion of alloys like 316 H, while not forming any detectable Ni-Mg intermetallic phases on Ni-rich alloy surfaces.

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

Probing the Mechanism of Cadmium Sulfide Cluster Nucleation and Growth via Sequential Infiltration Synthesis

Few-atom metal chalcogenide clusters may be realized through the sequential infiltration of metal–organic precursors in polymer films. However, the underlying nucleation and growth mechanisms that allow for cluster synthesis with near atomic-scale precision are not fully resolved. The kinetics of the sequential infiltration synthesis (SIS) method that control the nucleation and growth mechanisms of primarily Cd 4 S 4 -core clusters within a poly­(4-vinylpyridine) (P4VP) matrix are probed with in situ UV–visible absorbance spectroscopy. Density functional theory (DFT) calculations that allow simulation of the optical properties of cluster fragments further reveal the thermodynamics that guide cluster growth within the P4VP matrix. Here, we conclude that a reactive capture mechanism for cluster nucleation and growth is dominant, although transient dimethyl cadmium adduction to the polymer backbone may contribute to cluster nucleation under shorter metal–organic purge process conditions. Grazing incidence X-ray diffraction (GI-XRD) and X-ray absorption spectroscopy (XAS) analyses further corroborate the cluster size and atom connectivity throughout the stepwise synthesis.

Jayaweera, Nuwanthaka P. [Argonne National Laborat↗
DOE OSTI2023

A Membrane Contactor Enabling Energy-Efficient CO 2 Capture from Point Sources with Deep Eutectic Solvents

We demonstrate a scalable and energy-efficient hollow fiber membrane contactor (HFMC)-based process using a green solvent for CO 2 capture. This process uses a deep eutectic solvent (DES) in an HFMC to provide close interfacial interactions and contact between the DES and CO 2 . This approach overcomes disadvantages associated with direct absorption in DES and could potentially be applied to a variety of solvent-based CO 2 capture methods. Commercial low-cost polymer hollow fiber membranes (e.g., microporous polypropylene) were evaluated for CO 2 capture with reline, a prototypical DES. Single-gas measurements showed that the DES-based polypropylene HFMC can capture and separate CO 2 while rejecting N2. From a mixed gas containing 50 mol % N 2 and 50 mol % CO ν , the DES-based HFMC separated CO 2 with a purity of 96.9 mol %. The effect of several process parameters including solvent flow rate, pressure, and temperature on the CO 2 separation performance was studied. The flux of the recovered CO 2 was 67.43 mmole/m 2 /h at a feed pressure of 4 bar. In situ Fourier transform infrared (FTIR) measurements combined with density functional theory (DFT)-based molecular dynamics simulations revealed that reline absorbs CO 2 by physical absorption without forming a new chemical compound, and CO 2 separation by reline occurs via the pressure swing mechanism. This research provides fundamental insights about physical solvent-based separation processes and a pathway toward practical deployment.

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

Unraveling Adsorbate-Induced Structural Evolution of Iron Carbide Nanoparticles

Iron carbide (Fe x C y ) nanoparticles (NPs) are promising candidates for replacing platinum group metals in industrial applications, such as high-temperature Fischer–Tropsch synthesis. However, due to their amorphous nature, characterization of the active sites has been challenging experimentally and computationally. Here, using a combined density functional theory (DFT), neural network interatomic potential-assisted global optimization, and ensemble learning study, we evaluate dynamic surface changes associated with syngas (H and CO) interactions. For this purpose, we have developed a general procedure that we use to model an experimentally relevant 270-atom Fe 182 C 88 NP using the neural network-assisted stochastic surface walk global optimization algorithm (SSW-NN). Once generated, the Fe 182 C 88 NP active sites and particle morphology are thoroughly characterized before the effects of syngas adsorbate interactions are explored by using DFT and molecular dynamics simulations. Lastly, we explore correlations between geometric and electronic features of the active sites and the adsorption of H (H ads ), using a regularized random forest machine learning algorithm. In doing so, we identified the Fe–C coordination number and p orbital occupancy as the most important descriptors affecting H ads . Furthermore, using a combined ML and quantum chemistry approach, our work demonstrates a general and efficient procedure for generating and probing complex surface phenomena on binary nanoparticles.

Adsorption↗
DOE OSTI2025

Operando Characterization of Fe in Doped Ni x (Fe 1– x )O y H z Catalysts for Electrochemical Oxygen Evolution

Iron-doped nickel oxyhydroxides, Ni x (Fe 1– x )O y H z , are among the most promising oxygen evolution reaction (OER) electrocatalysts in alkaline environments. Although iron (Fe) significantly enhances the catalytic activity, there is still no clear consensus on whether Fe directly participates in the reaction or merely acts as a promoter. To elucidate the Fe’s role, we performed operando X-ray spectroscopy studies supported by DFT on Ni x (Fe 1– x )O y H z electrocatalysts. We probed the reversible changes in the structure and electronic character of Ni x (Fe 1– x )O y H z as the electrode potential is cycled between the resting (here at 1.10 V RHE ) and operational states (1.66 V RHE ). DFT calculations and XAS simulations on a library of Fe structures in various NiO y H z environments are in favor of a distorted local octahedral Fe(III)O 3 (OH) 3 configuration at the resting state with the NiO y H z scaffold going from α-Ni(OH) 2 to γ-NiOOH as the potential is increased. Under catalytic conditions, EXAFS and HERFD spectra reveal changes in p-d mixing (covalency) relative to the resting state between O/OH ligands and Fe leading to a shift from octahedral to square pyramidal coordination at the Fe site. XES measurements and theoretical simulations further support that the Fe equilibrium structure remains in a formal Fe(III) state under both resting and operational conditions. These spectral changes are attributed to potential dependent structural rearrangements around Fe. The results suggest that ligand dissociation leads to the C 4v symmetry as the most stable intermediate of the Fe during OER. This implies that Fe has a weakly coordinated or easily dissociable ligand that could serve to coordinate the O–O bond formation and, tentatively, play an active role in the Ni x (Fe 1– x )O y H z electrocatalyst.

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

Breaking up CO 2 through pressure-induced redox reaction with rhenium

We present results from density-functional-theory (DFT) based molecular dynamics simulations of fluid mixtures of carbon dioxide (CO 2 ) and rhenium (Re) metal at pressures up to 45 GPa and a temperature of 4000 K. Covalently-bonded carbon chains form as a result of a redox reaction between CO 2 and Re, whereas a pure CO 2 fluid is stable against decomposition at the same conditions. Furthermore, these findings have implications in reducing CO 2 by pressure-induced redox reactions with transition metals, which may play a role in CO 2 destruction and carbon storage as related to possible solutions to the current climate crisis.

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

Multiphase equation of state for Ta 2 O 5

A new equation of state for Ta 2 O 5 is presented. The EOS is constructed using the OpenSesame software and is referred to as SESAME 3530. The EOS uses a combination of density functional theory (DFT) calculations and experimental data. DFT calculations include cold curves and phonons of the solid phases, as well as DFT-based molecular dynamics simulations of the liquid phase. Experimental data includes isobaric, diamond anvil cell, and porous shock Hugoniot data. To fit the data, we create a multiphase EOS consisting of two solid phases and the liquid. Overall agreement with experimental data is shown, and we provide some suggestions for future experiments that could improve our knowledge of the phase diagram.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Understanding Microscopic Mechanisms of LeTID and LID and their Unifying Features by Electron Paramagnetic Resonance

LID degradation involves not only creation of ~ 10^12 cm^-3 recombination centers, but also ~ 10^16 cm^-3 shallow negative-U traps. In Ga-doped Si, LID EPR defects don't appear, but some traps are still created. LeTID: Si DB and H-hyperfine EPR signatures. We postulate that the defect responsible for LeTID is a partially hydrogenated (multivacancy) with a Si dangling bond and H in the vicinity. O involvement is possible yet unclear. We prove that H is related to the structure of the LeTID defect with isotope experiments and its EPR signal is comparable and linear with the Si DB signal upon LeTID degradation. Working on simulating these results with DFT to obtain more detailed defect structure.

Cz Si↗
DOE OSTI2024

Ni and Co Incorporation in Forsterite: A Density Functional Theory Study with Hubbard Correction

Ni and Co are critical elements for the world economy and modern technologies. Mafic and ultramafic deposits represent low-grade yet abundant alternatives to traditional Ni and Co ores. Here, in this work, density functional theory (DFT) with the Hubbard U correction (DFT+U) was used to simulate the incorporation of Ni and Co in forsterite (Mg 2 SiO 4 ), the Mg endmember of olivine, a common mineral in mafic and ultramafic rocks. Hubbard U terms for Ni and Co were parametrized using a series of oxide, hydroxide, carbonate, silicate, and sulfide minerals relevant to extraction and recovery of Ni and Co from mafic and ultramafic deposits. Electronic, energetic, magnetic, and structural properties were considered in the parametrization. For each of Ni and Co, an effective Hubbard correction (U eff ) value that optimized agreement with either experimental data or a hybrid exchange-correlation functional for all of the minerals considered is reported. DFT+U ab initio molecular dynamics (AIMD) simulations of Ni and Co incorporated into the M1 and M2 octahedral sites of forsterite were then performed. Ni and Co substitution in the M1 site was more energetically favorable than substitution in the M2 site, in agreement with published partition coefficients. AIMD trajectories were used to compute extended X-ray absorption fine structure (EXAFS) spectra of Ni in the M1 and M2 sites for direct fitting to a published experimental spectrum of Ni in a natural San Carlos olivine sample. The results of the fit indicated that ordering of Ni in the M1 site was not as strong at the low Ni concentrations relevant to mafic and ultramafic silicate minerals as that at the higher concentrations of the Ni-Mg olivine solid solutions studied to date.

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

Exploring Cr and molten salt interfacial interactions for molten salt applications

Molten salts play an important role in various energy-related applications such as high-temperature heat transfer fluids and reaction media. However, the extreme molten salt environment causes the degradation of materials, raising safety and sustainability challenges. A fundamental understanding of material–molten salt interfacial evolution is needed. This work studies the transformation of metallic Cr in molten 50/50 mol% KCl–MgCl 2 via multi-modal in situ synchrotron X-ray nano-tomography, diffraction and spectroscopy combined with density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. Notably, in addition to the dissolution of Cr in the molten salt to form porous structures, a δ-A15 Cr phase was found to gradually form as a result of the metal–salt interaction. This phase change of Cr is associated with a change in the coordination environment of Cr at the interface. DFT and AIMD simulations provide a basis for understanding the enhanced stability of δ-A15 Cr vs. bcc Cr, by revealing their competitive phase thermodynamics at elevated temperatures and probing the interfacial behavior of the molten salt at relevant facets. This study provides critical insights into the morphological and chemical evolution of metal–molten salt interfaces. Finally, the combination of multimodal synchrotron analysis and atomic simulation also offers an opportunity to explore a broader range of systems critical to energy applications.

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

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

Influence of Cation Size on the Local Atomic Structure and Electronic Properties of Ta Perovskite Oxynitrides

Partial anion substitution in transition metal oxides provides rich opportunities to control and tune physical and chemical properties, for example, combining the merits of oxides and nitrides. In addition, the possibility of resulting anion sublattice order provides a means to target polar and chiral structures based on a wide array of accessible structural archetypes by design. Here, in this work, we investigate the local structures of a family of perovskite tantalum oxynitrides—ATaO 2 N (A = Ba, Sr, and Ca)—using a combination of experimental and theoretical approaches including neutron total scattering, density functional theory (DFT), and ab initio molecular dynamics (AIMD) simulations. We present the first experimental study of chemical short-range order (CSRO) in CaTaO 2 N, confirming local cis N ordering of the anion sub-lattice. Our systematic exploration of a local structure across the A cation size series (from the larger Ba to the smaller Ca) reveals a perovskite motif increasingly distorted with respect to long-range average structures. DFT and AIMD simulations support the observed trends. Overall, structures with cis ordering of the nitrogen anions in each TaO 4 N 2 octahedron are favored over those with trans ordering. With diminishing A cation size, local cis ordering and Ta off-centering play decreasing roles in overall lattice stability, overshadowed by the stabilizing effects of octahedral tilting. The influence of these factors on local dipole formation and frustrated dipole ordering are discussed.

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

Revealing Intermetallic Active Sites of PtNi Nanocatalysts for Reverse Water Gas Shift Reaction

Bimetallic nanoparticles have gained great interest due to their potential applications in catalysis. By combining Lattice Boltzmann Modeling (LBM), Density Functional Theory (DFT) calculations, and kinetic Monte Carlo (KMC) simulation, here we rationalize the superior performance of supported Pt-Ni alloy nanostructures over Pt and Ni alone during the reverse Water Gas Shift (rWGS) reaction, which was observed previously. The LBM simulation including the diffusion and flow of the reacting gas molecules within the supporting nanoporous channels in SBA-15 indicates that the simple combination of contributions from Ni and Pt alone is not enough to result in the observed promotion in rWGS activity and selectivity of the Pt-Ni alloy. Instead, the synergy between Pt and Ni via a unique intermetallic motif is responsible according to the DFT and KMC simulation. In conclusion, compared to the monometallic counterparts, the intermetallic-introduced ensemble effect at the low-coordinated Pt-Ni hybrid site of Pt-Ni alloy nanocatalysts enables the enhanced activation of CO 2 and promoted removal of CO/H 2 O as gas phase product.

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

Many-Body Benchmark of Electronic Charge and Spin Densities for Li 1–x NiO 2

Accurate benchmarks are particularly important for highly correlated oxides as mean-field approximations often fail to describe the subtle balance of charge transfer and magnetism in these materials with an accuracy comparable to experimental needs. Here we present accurate diffusion Monte Carlo (DMC) results of the electronic charge and spin densities for the tunable highly correlated oxide Li 1–x NiO 2 for x = 0, 1/2, and 1. To enable quantitative comparisons, we introduce a robust density-partitioning scheme, extending Voronoi analysis to assign atomic charges from spatially noisy DMC densities. We then benchmark common approximations used in density functional theory (DFT). Comparison against DMC shows that r 2 SCAN delivers the most balanced performance across charge, spin, and radial density descriptors, nearly reproducing DMC results for LiNiO 2 and apical Ni sites in Li 0.5 NiO 2 . Hybrid functionals (PBE0, SCAN0) perform unexpectedly poorly, and PBE + U + V yields inconsistent trends between charge and spin densities. Therefore, the r 2 SCAN functional minimizes errors relative to DMC while capturing the variable valence of the Ni ion and also retaining the computational efficiency of DFT for large-scale simulations of the tunable structural and electronic phases of Li1−xNiO2. Our study highlights the importance of accurate benchmarking of the fundamental quantities involved in DFT to select appropriate DFT approximations in order to advance the predictive modeling of charge-transfer-driven phenomena in correlated electron systems.

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

Modeling separation of lanthanides via heterogeneous ligand binding

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln 3+ ) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln 3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln 3+ coordinated to ligands in water and embedded in metal–organic frameworks (MOFs), and ΔΔE for Ln 3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln 3+ in the periodic table. Static DFT calculations with a full Ln 3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. Finally, we also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

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

Electronic density response of warm dense matter

Matter at extreme temperatures and pressures—commonly known as warm dense matter (WDM)—is ubiquitous throughout our Universe and occurs in astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are probed in x-ray Thomson scattering (XRTS) experiments and are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear response theory and nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques, including path-integral Monte Carlo simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrary complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗