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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High spin polarization in the disordered quaternary Heusler alloy FeMnVGa

Here in this paper, we report the successful synthesis of a Fe-based highly spin-polarized quaternary Heusler alloy FeMnVGa and its structural, magnetic, and transport properties probed through different experimental methods and theoretical techniques. Density functional theory (DFT) calculations performed on different types of structures reveal that the structure with Ga at 4a, V at 4b, Mn at 4c, and Fe at 4d (space group F$\bar{4}$3m) possess minimum energy among all the ordered variants. Ab initio simulations in the most stable ordered structure show that the compound is a ferromagnet having a large spin-polarization (89.9%). Neutron diffraction reveals that the compound crystallizes in disordered type-2 structure (space group Fm$\bar{3}$m) in which Ga occupies at 4a, V 4b and Fe/Mn occupy 4c/4d sites with 50:50 proportions. The structural disorder is further confirmed by x-ray diffraction, extended x-ray absorption fine structure, 57 Fe Mössbauer spectrometry results, and DFT calculations. Magnetization studies suggest that the compound orders ferromagnetically below T C ~ 293 K and the saturation magnetization follows the Slater-Pauling rule. Mössbauer spectrometry, along with neutron diffraction, suggest that Mn is the major contributor to the total magnetism in the compound, consistent with the theoretical calculations, which also indicates that spin polarization remains high (81.3%), even in the presence of such large atomic disorder. The robustness of the half-metallic ferromagnetic (HMF) property in the presence of disorder is a quite unique characteristic over other reported HMF in literature and makes this compound quite promising for spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Intermetallic alloy structure–activity descriptors derived from inelastic X-ray scattering

Synchrotron spectroscopy and Density Functional Theory (DFT) are combined to develop a new descriptor for the stability of adsorbed chemical intermediates on metal alloy surfaces. This descriptor probes the separation of occupied and unoccupied d electron density in platinum and is related to shifts in Resonant Inelastic X-ray Scattering (RIXS) signals. Simulated and experimental spectroscopy are directly compared to show that the promoter metal identity controls the orbital shifts in platinum electronic structure. The associated RIXS features are correlated with the differences in the band centers of the occupied and unoccupied d bands, providing chemical intuition for the alloy ligand effect and providing a connection to traditional descriptions of chemisorption. The ready accessibility of this descriptor to both DFT calculations and experimental spectroscopy, and its connection to chemisorption, allow for deeper connections between theory and characterization in the discovery of new catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring the Effects of Node Topology, Connectivity, and Metal Identity on the Binding of Nerve Agents and Their Hydrolysis Products in Metal–Organic Frameworks

Recent studies have shown that metal–organic frameworks (MOFs) built from hexanuclear M(IV) oxide cluster nodes are effective catalysts for nerve agent hydrolysis, where the properties of the active sites on the nodes can strongly influence the reaction energetics. The connectivity and metal identity of these M6 nodes can be easily tuned, offering extensive opportunities for computational screening to predict promising new materials. Thus, we used density functional theory (DFT) to examine the effects of node topology, connectivity, and metal identity on the binding energies of multiple nerve agents and their corresponding hydrolysis products. By computing an optimization metric based on the relative binding strengths of key hydrolysis reaction species (water, agent, and bidentate-bound products), we predicted optimal M6 nodes for hydrolyzing specific nerve agent and simulant molecules, where our results are in qualitative agreement with observed experimental trends. This analysis highlighted the notion that no single metal or node topology is optimal for all possible organophosphates, suggesting that MOFs should be selected based on the agent of interest. Using the large amount of data generated from our DFT calculations, we then derived quantitative structure–activity relationship (QSAR) models to help explain the complex trends observed in the binding energies. Through linear regression, we identified the most important descriptors for describing the binding of nerve agents and their hydrolysis products to M6 nodes. These results suggested that both molecular and node properties, including both structural and chemical features, collectively contribute to the binding energetics. By performing a thorough statistical analysis, we showed that our QSAR models are capable of making quantitatively accurate binding energy predictions for nerve agents and their hydrolysis products in a wide variety of M(IV)-MOFs. In conclusion, the insights gained herein can be used to guide future experiments for the synthesis of MOFs with enhanced catalytic activity for organophosphate hydrolysis.

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

Enhancing the performance of lithium oxygen batteries through combining redox mediating salts with a lithium protecting salt

Li–O 2 batteries have recently emerged to meet nowadays elevated electric energy demands. Redox mediators (RMs) for solution-inducing decomposition of discharge products are one approach to increase energy efficiency and reduce high overpotentials in these batteries. However, multiple obstacles hinder their usage such as redox shuttling, capacity fading, electrolyte degradation, etc. Herein, we present a new chemistry based on a combination of LiNO 3 , TEGDME and an ionic liquid that enables LiI (1 M) to lower the charge potential (3.5V) with a long cycle life of 270 cycles. 0.1 M LiI increases the cyclability up to 500 with a slightly increased charge potential (~4V) for a fixed capacity of 1000 mAh/g. Up to 100 cycles, this battery system retained ~95% Li 2 O 2 capacity with a ~0.8 V charge-discharge polarization gap. The addition of LiNO 3 to the electrolyte provides a protective solid electrolyte interface (SEI) on anode that works in synergy with the LiI RM. Moreover, we found that this electrolyte blend results in domain formation of ionic and neutral species enhancing the discharge and charge processes. Finally, DFT calculations provide a better understanding of the role of the anode SEI layer and the Li 2 O 2 decomposition promoted by the LiI during charge on the cathode.

25 ENERGY STORAGE↗
DOE OSTI2021

Understanding the origin of structure sensitivity in hydrodechlorination of trichloroethylene on a palladium catalyst

Periodic density functional theory (DFT) calculations are employed to understand the origin of structure sensitivity in hydrodechlorination (HDC) of trichloroethylene (TCE) over different facets of a palladium catalyst. The HDC reaction is simulated on the terrace (Pd (111) and Pd (100)) and undercoordinated (Pd (211) and Pd (110)) sites of the Pd catalyst. The most stable binding configuration of TCE on the Pd surfaces is observed to be through the di-σ mode of binding, wherein each carbon atom of the TCE molecule is adsorbed atop of the Pd atom. On comparing TCE adsorption over different facets of Pd, a maximum binding energy of –178 kJ mol –1 is calculated over the Pd (110) surface. TCE, upon adsorption on Pd catalyst, undergoes dechlorination followed by hydrogenation of the hydrocarbon intermediates. The activation energies for C–Cl bond dissociation steps are significantly low when compared to the hydrogenation steps. The chlorine released from dechlorination tends to block the active sites, thereby poisoning the surface with high binding energies (B.E > –160 kJ mol –1 ) on all the surfaces. The trend in chlorine binding energies on Pd facets follows: Pd (110) > Pd (211) > Pd (100) > Pd (111). The removal of surface chlorine is facilitated by its reaction with surface hydrogen to form hydrogen chloride. The activation energy for hydrogen chloride formation is calculated to be 90 kJ mol –1 and 88 kJ mol –1 on Pd (111) and Pd (100) terrace sites, respectively and 109 kJ mol –1 and 118 kJ mol –1 on the step Pd (211) and corrugated Pd (110) facets, respectively. This suggests the ease of removal of Cl as HCl from the terrace sites as compared to the step and corrugated sites. The structure sensitivity in the TCE HDC reaction could possibly arise due to the differences in the energetics of Cl removal on different Pd facets. As a result, this mechanistic understanding could provide a rationale for designing suitable catalysts for the HDC of TCE.

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

Liquid-liquid phase transition of hydrogen and its critical point: Analysis from ab initio simulation and a machine-learned potential

We simulate high-pressure hydrogen in its liquid phase close to molecular dissociation using a machine-learned interatomic potential. The model is trained with density functional theory (DFT) forces and energies, with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. We show that an accurate NequIP model, an E(3)-equivariant neural network potential, accurately reproduces the phase transition present in PBE. Moreover, the computational efficiency of this model allows for substantially longer molecular dynamics trajectories, enabling us to perform a finite-size scaling (FSS) analysis to distinguish between a crossover and a true first-order phase transition. Here, we locate the critical point of this transition, the liquid-liquid phase transition (LLPT), at 1200-1300 K and 155-160 GPa, a temperature lower than most previous estimates and close to the melting transition.

08 HYDROGEN↗
DOE OSTI2020

Prediction of the Al-rich part of the Al-Cu phase diagram using cluster expansion and statistical mechanics

In this work, the thermodynamic properties of α-Al and other phases (GP zones, θ'', θ' and θ) in the Al-rich part of the Al-Cu system have been obtained by means of the cluster expansion formalism in combination with statistical mechanics. This information was used to build the Al-rich part of the Al-Cu phase-diagram taking into account vibrational entropic contributions for θ', as those of the other phases were negligible. The simulation predictions of the phase boundaries between α-Al and either θ'', θ' or θ phases as a function of temperature are in good agreement with experimental data and extend the phase boundaries to a wider temperature range. The DFT calculations reveal the presence of a number of metastable Guinier-Preston-zone type configurations that may coexist with α-Al and θ'' at low temperatures. They also demonstrate that θ' is the stable phase below 550K but it is replaced by θ above this temperature due to the vibrational entropic contribution to the Gibbs energy of θ'. This work shows how the combination of cluster expansion and statistical mechanics can be used to expand our knowledge of the phase diagram of metallic alloys and to provide Gibbs free energies of different phases that can be used as input in mesoscale simulations of precipitation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting Morphological Evolution during Coprecipitation of MnCO 3 Battery Cathode Precursors Using Multiscale Simulations Aided by Targeted Synthesis

The performance of lithium-ion batteries is intimately linked to both the structure and the morphology of the cathode material, which in turn is critically linked to the synthesis conditions. However, few studies focus on understanding synthesis, especially during the coprecipitation of metal oxide precursors, a process that largely determines the final morphology of the material. In this paper, we go beyond the typical equilibrium particle shape analysis conducted in the literature and incorporate kinetic aspects of morphology evolution. We perform these studies using controlled synthesis on a well-defined metal salt system (MnCO 3 ) combined with multiscale simulations and high-resolution microscopy. Results show that with increasing metal concentration, the particles transition from rhombohedral to cubic to spherical shapes. Computational analysis using density functional theory (DFT) reveals that rhombohedral shaped particles evolve under equilibrium conditions. Phase field techniques indicate that at higher metal concentrations, fast growth kinetics of the precipitates result in the transition to cubic and, subsequently, spherical shapes, accompanied by a decrease in particle size. This study, while limited to the one metal salt system, provides an approach to shed light on the synthesis process of mixed transition metal salts, gradient materials, and other cathode materials of interest to the battery community.

25 ENERGY STORAGE↗
DOE OSTI2025

Reactivity of Carbonate Solvent Electrolytes on Lithium Silicon Anodes

Silicon (Si) is promising for lithium-ion battery (LIB) anodes due to their high theoretical capacity and low electrochemical potential. However, significant challenges remain, including severe volumetric expansion during cycling and the electrochemical instability of electrolytes, which leads to the formation of a nonuniform solid electrolyte interphase (SEI). To investigate SEI formation mechanisms, computational molecular dynamics simulations offer valuable insights. In this work, we examine the trajectories and charge transfer behavior of lithium hexafluorophosphate (LiPF 6 ) salt with various solvent compositions using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Among the tested electrolyte systems, LiPF 6 with vinylene carbonate (VC) added to ethyl methyl carbonate (EMC) exhibits the lowest reactivity with the Si anode. In contrast, the effects of fluoroethylene carbonate (FEC) and VC depend on whether the primary solvent is EMC alone or a mixture of ethylene carbonate (EC) and EMC. Moreover, we show that electrolyte reactivity varies with the degree of lithiation of the Si anode (LiSi vs Li 15 Si 4 ) and under different charge states. To decouple electrolyte reactivity from surface effects, we analyze the dissociation and formation energies of individual species from solvated configurations. Overall, these first-principles-based findings provide a strategic foundation for electrolyte design to improve cycling stability and extend calendar life in LIBs using Si anodes.

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

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗
DOE OSTI2020

Extraction of Ln(III) and An(III) from Alkaline Solutions by Sulfonamide Ligands and Analogs - 20216

As part of our effort to develop advanced separations involving combined extraction of cesium, strontium, and actinides from High-Level Waste (HLW) at the Savannah River Site, we are now reporting highly efficient extraction of Sm(III) from aqueous media of various alkalinity by dichloromethane solutions of several ligands, including tri-sulfonamides, o-mono-sulfonamido-phenols and o-phenylenediamine- derived disulfonamides. Distribution experiments and DFT calculations provide insights on the complexation of Samarium vs Americium by these ligands. Further optimization of these ligands for extraction into highly lipophilic solvents and application to HLW simulants are underway. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Prediction of Feasibility of Polaronic OER on (110) Surface of Rutile TiO 2

The polaronic effects at the atomic level hold paramount significance for advancing the efficacy of transition metal oxides in applications pertinent to renewable energy. The lattice–distortion mediated localization of photoexcited carriers in the form of polarons plays a pivotal role in the photocatalysis. This investigation focuses on rutile TiO 2 , an important material extensively explored for solar energy conversion in artificial photosynthesis, specifically targeting the generation of green H 2 through photoelectrochemical (PEC) H 2 O splitting. By employing Hubbard-U corrected and hybrid density functional theory (DFT) methods, we systematically probe the polaronic effects in the catalysis of oxygen evolution reaction (OER) on the (110) surface of rutile TiO 2 . Theoretical understanding of polarons within the surface, coupled with simulations of OER at distinct titanium (Ti) and oxygen (O) active sites, reveals diverse polaron formation energies within the lattice sites with strong preference for bulk and surface bridge (O b ) oxygen sites. Moreover, we provide the evidence for the facilitative role of polarons in OER. Here, we find that hole polarons situated at the equatorial oxygen sites near the Ti–active site, along with bridge site hole polarons distal from the O b active site yield a small reduction in OER overpotential by ~0.06 eV and ~0.12 eV, respectively. However, subsurface, equatorial, and bridge site hole polarons significantly reduce the Ti-active site OER overpotential by ~0.4 eV through the peroxo–type oxygen pathway. We also observe that the presence of hole polarons stabilizes the *OH, *O, and *OOH intermediate species compared to the scenario without hole polarons. Overall, this study provides a detailed mechanistic insight into polaron–mediated OER, offering a promising avenue for improving the catalytic activity of transition metal oxide-based photocatalysts catering to renewable energy requisites.

14 SOLAR ENERGY↗
DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fundamentals of Au(111) Surface Dynamics: Coarsening of Two-Dimensional Au Islands

Au(111) surfaces play a central role in many applications, yet studies of fundamental aspects of their dynamics are limited. Thus, using scanning tunneling microscopy (STM) at 300 K, we analyze the coarsening of first-layer 2D Au islands directly on the Au(111) substrate and also of second-layer 2D Au islands. Specifically, we monitor the decay of Au first-layer islands with areas of about 100–500 nm 2 in the vicinity of larger islands or extended step edges over a period of approximately 40 h, the relevant time scale for this process. Experimentally observed behavior is captured by analytical theory for terrace-diffusion-limited decay incorporating DFT results for the Au terrace diffusion barrier and the adatom formation energy.Lastly, experimental observations of second-layer island decay are also compared with appropriate analytical theory and stochastic simulations, thereby determining the effective Ehrlich–Schwoebel barrier for Au on Au(111).

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

Cation Effects on CO 2 Delivery to Cu Electrode in Reactive Capture of CO 2

The direct electrochemical conversion of captured CO 2 , known as reactive capture of CO 2 (RCC), remains a formidable challenge in heterogeneous catalysis. Given that amines are one of the most widely used capture agents for CO 2 , it would be desirable to electrochemically reduce the resultant adducts, such as carbamate, directly in RCC. However, current understanding suggests that the primary species undergoing reduction in RCC with amines is the CO 2 dissociated from the sorbent. Herein, we employ ab initio molecular dynamics (AIMD) with DFT to analyze how the nature of alkali metal cations in the electrolyte affects carbamate at the Cu surface, thereby assessing the possibility of promoting RCC by cation effects. The simulations show that the carbamate’s orientation with respect to the electrode is governed by the optimal distance between the carbamate and the cation, specifically how this distance aligns with the cation’s hydration spheres. Moreover, the slow-growth AIMD results indicate that the CO 2 dissociation barrier correlates with the orientation of carbamate at the interface. When the carbamate resides beyond the cation’s first hydration sphere, it adopts a flat orientation with respect to the surface that promotes the release of CO 2 from the capture agent. In contrast, when the carbamate disrupts the first hydration sphere and exhibits a strong cation−π interaction, it adopts an upright orientation that is less conducive to CO 2 release. These findings reveal a nontrivial cation effect in RCC, suggesting that it should be possible to optimize RCC via the choice of the electrolyte.

ab initio molecular dynamics↗
DOE OSTI2023

Machine learning interatomic potential for silicon-nitride (Si 3 N 4 ) by active learning

Silicon nitride (Si 3 N 4 ) is an extensively used material in the automotive, aerospace, and semiconductor industries. However, its widespread use is in contrast to the scarce availability of reliable interatomic potentials that can be employed to study various aspects of this material on an atomistic scale, particularly its amorphous phase. In this work, we developed a machine learning interatomic potential, using an efficient active learning technique, combined with the Gaussian approximation potential (GAP) method. Our strategy is based on using an inexpensive empirical potential to generate an initial dataset of atomic configurations, for which energies and forces were recalculated with density functional theory (DFT); thereafter, a GAP was trained on these data and an iterative re-training algorithm was used to improve it by learning on-the-fly. When compared to DFT, our potential yielded a mean absolute error of 8 meV/atom in energy calculations for a variety of liquid and amorphous structures and a speed-up of molecular dynamics simulations by 3–4 orders of magnitude, while achieving a first-rate agreement with experimental results. Our potential is publicly available in an open-access repository.

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

Challenges for density functional theory in simulating metal–metal singlet bonding: A case study of dimerized VO2

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V–V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V–V dimer length. The spin-restricted method tends to overestimate the strength of the V–V bonds, resulting in a small V–V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron–electron repulsion, involved in the metal–insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations.

Chemistry↗