Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atomistic models”

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

Complexions at the iron-magnetite interface

Synthesizing distinct phases and controlling crystalline defects are key concepts in materials design. These approaches are often decoupled, with the former grounded in equilibrium thermodynamics and the latter in nonequilibrium kinetics. By unifying them through defect phase diagrams, we can apply phase equilibrium models to thermodynamically evaluate defects—including dislocations, grain boundaries, and phase boundaries—establishing a theoretical framework linking material imperfections to properties. Using scanning transmission electron microscopy (STEM) with differential phase contrast (DPC) imaging, we achieve the simultaneous imaging of heavy Fe and light O atoms, precisely mapping the atomic structure and chemical composition at the iron-magnetite (Fe/Fe 3 O 4 ) interface. We identify a well-ordered two-layer interface-stabilized phase state (referred to as complexion) at the Fe[001]/Fe 3 O 4 [001] interface. Using density-functional theory (DFT), we explain the observed complexion and map out various interface-stabilized phases as a function of the O chemical potential. The formation of complexions increases interface adhesion by 20% and alters charge transfer between adjacent materials, impacting transport properties. Our findings highlight the potential of tunable defect-stabilized phase states as a degree of freedom in materials design, enabling optimized corrosion protection, catalysis, and redox-driven phase transitions, with applications in materials sustainability, efficient energy conversion, and green steel production.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Probing degradation at solid-state battery interfaces using machine-learning interatomic potential

Solid-state batteries featuring fast ion-conducting solid electrolytes are promising next-generation energy storage technologies, yet challenges remain for practical deployment due to electro-chemo-mechanical instabilities at solid-solid interfaces. These interfaces, which include homogeneous/internal interfaces such as grain boundaries (GBs) and heterogeneous/external interfaces between solid-electrolyte and electrode materials, can impede Li-ion transport, deteriorate performance, and eventually lead to cell failure. Here, in this study, we leverage large-scale molecular simulations, enabled by validated machine-learning interatomic potentials, to directly probe the onset of interfacial degradation at the garnet Li 7 La 3 Zr 2 O 12 (LLZO) solid-electrolyte/LiCoO 2 (LCO) cathode interface. By surveying different interfacial geometries and compositions, it is found that Li-deficient interfaces can lead to severe interfacial disordering with cation mixing and Co interdiffusion from LCO into LLZO. By contrast, Li-sufficient interfaces are less disordered, although elemental segregation with local ordering is observed. As a consequence of Co interdiffusion, Co-rich regions are formed at the GBs of LLZO due to cation segregation and trapping effects. This behavior is independent of the GB tilting axis, degree of disorder at the GBs, and Co concentration, which implies Co clustering at GBs is a general phenomenon in polycrystalline LLZO and can dictate its overall transport and mechanical properties. Our findings elucidate the underlying fundamental mechanisms that give rise to experimentally observed physicochemical properties and provide guidelines for interface design that can mitigate interfacial degradation and improve cycling performance.

25 ENERGY STORAGE↗
DOE OSTI2024

Crystal structure and shape selection in the growth of 3D metallic crystallites on layered materials: Fe on MoS 2

Nucleation and growth of supported 3D metal clusters or crystallites during deposition on MoS 2 , or on other weakly-adhering layered materials, can potentially produce diverse growth shapes, and even crystal structures differing from the bulk metal. For Fe deposition on MoS 2 , SEM and AFM observations reveal three distinct crystallite shapes. By comparison with atomistic structure models incorporating realistic Fe-MoS 2 interface structures, here we conclude that these are: triangular fcc(111) pyramids with sloped {100} side facets; bcc(110) A-frame tents with sloped {100} side facets; and bcc(110) mesas with vertical {100} and {110} side facets. The following picture is proposed for the competitive formation of clusters and crystallites with different structures: (i) small nanoclusters formed at the onset of deposition exhibit facile fluxional dynamics allowing sampling of different crystal structures and shapes; (ii) sufficient fluxionality implies a Boltzmann distribution of sampled structures, and thus coexistence of different structures follows from the demonstrated similar energies for those structures; (iii) growing clusters reach a threshold size above which the characteristic time scale for restructuring exceeds that for cluster growth. Thereafter, clusters are locked-in to a specific crystal structure and shape as revealed by imaging of larger crystallites. Despite a penalty for fcc(111) over bcc(111) pyramids based on bulk energetics, favorable surface and interface energies makes them preferable for smaller sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Grand canonically optimized grain boundary phases in hexagonal close-packed titanium

Grain boundaries (GBs) profoundly influence the properties and performance of materials, emphasizing the importance of understanding the GB structure and phase behavior. As recent computational studies have demonstrated the existence of multiple GB phases associated with varying the atomic density at the interface, we introduce a validated, open-source GRand canonical Interface Predictor (GRIP) tool that automates high-throughput, grand canonical optimization of GB structures. While previous studies of GB phases have almost exclusively focused on cubic systems, we demonstrate the utility of GRIP in an application to hexagonal close-packed titanium. We perform a systematic high-throughput exploration of tilt GBs in titanium and discover previously unreported structures and phase transitions. In low-angle boundaries, we demonstrate a coupling between point defect absorption and the change in the GB dislocation network topology due to GB phase transformations, which has important implications for the accommodation of radiation-induced defects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Feature engineering descriptors, transforms, and machine learning for grain boundaries and variable-sized atom clusters

Abstract Obtaining microscopic structure-property relationships for grain boundaries is challenging due to their complex atomic structures. Recent efforts use machine learning to derive these relationships, but the way the atomic grain boundary structure is represented can have a significant impact on the predictions. Key steps for property prediction common to grain boundaries and other variable-sized atom clustered structures include: (1) describing the atomic structure as a feature matrix, (2) transforming the variable-sized feature matrix to a fixed length common to all structures, and (3) applying a machine learning algorithm to predict properties from the transformed matrices. We examine how these steps and different combinations of engineered features impact the accuracy of grain boundary energy predictions using a database of over 7000 grain boundaries. Additionally, we assess how different engineered features support interpretability, offering insights into the physics of the structure-property relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Substantial oxygen loss and chemical expansion in lithium-rich layered oxides at moderate delithiation

Delithiation of layered oxide electrodes triggers irreversible oxygen loss, one of the primary degradation modes in lithium-ion batteries. However, the delithiation-dependent mechanisms of oxygen loss remain poorly understood. Here we investigate the oxygen non-stoichiometry in Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ electrodes as a function of Li content by using cycling protocols with long open-circuit voltage steps at varying states of charge. Surprisingly, we observe substantial oxygen loss even at moderate delithiation, corresponding to 2.5, 4.0 and 7.6 ml O 2 per gram of Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ after resting at upper capacity cut-offs of 135, 200 and 265 mAh g −1 for 100 h. Our observations suggest an intrinsic oxygen instability consistent with predictions of high oxygen activity at intermediate potentials versus Li/Li + . In addition, we observe a large chemical expansion coefficient with respect to oxygen non-stoichiometry, which is about three times greater than those of classical oxygen-deficient materials such as fluorite and perovskite oxides. Furthermore, our work challenges the conventional wisdom that deep delithiation is a necessary condition for oxygen loss in layered oxide electrodes and highlights the importance of calendar ageing for investigating oxygen stability.

Atomistic models↗
DOE OSTI2024

Atomistic origin of the entropy of melting from inelastic neutron scattering and machine learned molecular dynamics

The latent heat, L, is central to melting, but its atomic origin remains elusive. It is proportional to the entropy of fusion, ΔS fus = L/T m (T m is the melting temperature), which depends on changes of atom configurations, atom vibrations, and thermal electron excitations. Here, we combine inelastic neutron scattering and machine-learned molecular dynamics to separate ΔS fus into these components for Ge, Si, Bi, Sn, Pb, and Li. When the vibrational entropy of melting, ΔS vib , is zero, ΔS fus ≃ 1.2 kB per atom. This result provides a baseline for ΔSconfig and nearly coincides with “Richard’s Rule” of melting. The ΔSfus deviates from this value for most elements, however, and we show that this deviation originates with extra ΔS vib and extra ΔS config . These two components are correlated for positive and negative deviations from Richard’s rule – the extra ΔS config is consistently ~ 80% of ΔS vib . Our results, interpreted with potential energy landscape theory, imply a correlation between the change in the number of basins and the change in the inverse of their curvature for the melting of pure elements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Combined molecular and spin dynamics simulation of BCC iron with vacancy defects

Utilizing an atomistic computational model, which handles both translational and spin degrees of freedom, combined molecular and spin dynamics simulations have been performed to investigate the effect of vacancy defects on spin wave excitations in ferromagnetic iron. Fourier transforms of space- and time-displaced correlation functions yield the dynamic structure factor, providing characteristic frequencies and lifetimes of the spin wave modes. A comparison of the system with a 5% vacancy concentration with pure lattice data shows a decrease in frequency and a decrease in lifetime for all transverse spin wave excitations observed. In addition, the clearly defined transverse spin wave excitations are distorted with the introduction of vacancy defects, and we observe reduced excitation lifetimes due to increased magnon–magnon scattering. We observe further evidence of increased magnon–magnon scattering, as the peaks in the longitudinal spin wave spectrum become less distinct. Finally, similar impacts are observed in the vibrational subsystem, with a decrease in characteristic phonon frequency and flattening of lattice excitation signals due to vacancy defects.

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

Evaluating the diffusion of Kr in UO 2 and ADOPT TM using time-of-flight elastic recoil detection analysis (ToF-erda)

A combination of 300 keV 84 Kr ion implantation and Time-of-Flight Elastic Recoil Detection Analysis is utilized to investigate the diffusion of Kr in UO 2 and ADOPT TM fuels. Composition depth-profiles on the nanometer scale were obtained, both for as-implanted samples and after annealing at 800°C for 1 hour. Observed drifts in the 84 Kr profiles could be associated with short-range diffusion mechanisms. The approach employed here provides the possibility to make direct comparisons with atomistic scale modelling data, and can be of service as a separate effect test in line with the Accelerated Fuel Qualification initiative.

ADOPT UO2↗
DOE OSTI2025

Multi-Scale Modeling of the Evolution of Structure and Properties in Materials for Nuclear Energy Applications [Slides]

Nuclear energy is an important component of an overall strategy to address climate change. Idaho National Laboratory (INL) is the U.S. Department of Energy’s primary facility for research and development in nuclear science and technology for energy generation, supporting the improvement and life extension of the existing reactor fleet and the development and licensing of new reactor designs. Computational modeling is an important component of these activities, particularly in the area of materials for nuclear applications, where experimental data can be very challenging and expensive to acquire, and where data is especially scarce for new reactor designs. INL has used multi-scale modeling – linking atomistic, mesoscale, and engineering scales – to improve the ability to predict the performance of materials for nuclear energy applications. In this talk, I will give an overview of the approach and tools used, and several examples of application, including performance of nuclear fuels, understanding radiation-driven formation of nanoscale void and gas bubble superlattices, and powder densification through electric field assisted sintering.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

The Development of Catalysts for Upgrading of Pyrolysis Vapor for Refinery Feedstocks and Intermediates (CRADA Final Report)

Catalytic fast pyrolysis (CFP) is a versatile technology platform to convert biomass into fungible hydrocarbon transportation fuels and chemical co-products. Key technical barriers to reaching this goal include increasing the product yields and achieving the desired fuel properties for gasoline, diesel, and jet range fuels or blendstocks that would be suitable for introduction into existing refinery unit operations. Overcoming these barriers will require durable catalysts that are effective at upgrading and stabilizing biomass pyrolysis vapors. Towards these goals, this CRADA leveraged NREL experience as a leader in biomass pyrolysis research and Johnson Matthey's (JM) experience as a leader in the production of advanced catalytic materials. The scope spanned CFP catalyst development, characterization, multi-scale reaction testing, and computational modeling. CRADA benefits to DOE, Participant, and U.S. Taxpayer: Assists laboratory in achieving programmatic scope, Uses the laboratory’s core competencies. The purpose of this CRADA was to develop and deploy catalysts for biomass CFP to help achieve cost-competitive biofuels and bio-based products. This was accomplished through a close collaboration between biomass conversion researchers at NREL and catalyst development researchers at JM. Summary of Research Results: Focus Area 1. Foundational research on catalytic conversion and deactivation: Key interactions between pyrolysis vapors and heterogeneous catalysts were probed through catalyst characterization, model compound reaction testing, and atomistic-scale computational modeling. Catalyst development focused on multifunctional materials, which include zeolites, oxides, carbides, and nitrides. Computational modeling identified reaction mechanisms and elucidated surface chemistry to test hypotheses regarding mechanisms of deoxygenation, coupling, cracking, dehydration, coke formation, hydrogen transfer, and aromatic ring reactions. This information was used to design multifunctional catalysts to increase product yields, control product selectivity, and reduce deactivation during CFP and downstream processing steps. The results served to increase fundamental understanding of key catalyst attributes and durability features for the upgrading of biomass pyrolysis vapors. Model compound experiments confirmed the importance of metal-acid bifunctionality for the deoxygenation of lignin-derived phenolic species under hydrodeoxygenation conditions. This insight led to the development of catalysts such as Pt/TiO2 and Mo2C, which were confirmed as high-performing materials during subsequent bench-scale experiments using biomass-derived pyrolysis vapors. This focus area also led to the identification of important catalyst deactivation mechanisms associated with the deposition of inorganic contaminants such as potassium. The molecular-level insight from model compound experiments and computational modeling, shown in Figure 1, informed the development of regeneration procedures that have been shown to be effective for restoration of > 90% of initial catalyst activity. This understanding has subsequently been translated to other catalyst systems, including zeolite materials that can be operated without requirements for co-fed hydrogen.

09 BIOMASS FUELS↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Systematic improvement of redox potential calculation of Fe(III)/Fe(II) complexes using a three-layer micro-solvation model

Electrochemical transformations of metal ions in aqueous media are challenging to model accurately due to the dynamic solvation structure surrounding ions at different charge states. Predictive modeling at the atomistic scale is essential for understanding these solvation architectures but is often computationally prohibitive. In this contribution, we present a simple, fast, and accurate three-layer micro-solvation model to evaluate the redox potential of metal ions in aqueous solutions. Our model, developed and validated for Fe 3+ /Fe 2+ redox potentials, combines the DFT-based geometry optimizations of the octahedral Fe complex with two layers of explicit water molecules to capture solute–solvent interactions and an implicit solvation model to account for bulk solvent effects. This approach yields accurate predictions for Fe 3+ /Fe 2+ redox potentials in water, achieving errors of 0.02 V with ωB97X-V, 0.01 V with ωB97X-D3, 0.04 V with ωB97M-V, and 0.02 V with B3LYP-D3 functionals. We further demonstrate the generality of our model by applying it to additional metal complexes, including the challenging Fe(CN) 6 3−/4− system, where our model successfully achieves close agreement with experimental values, with an error of 0.07 V and an average error of 0.21 V for all five systems. In summary, the presented simple solvation model has broad applicability and potential for enhancing computational efficiency in redox potential predictions across various chemical and industrial processes of metal ions.

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

A Universal Model of Cation Effects in Electrocatalysis

Electrolyte cations are conventionally viewed as inert spectators in electrocatalysis. However, a wealth of observations show that catalytic rates are often highly sensitive to cation identity. Despite their prevalence, these cation effects have resisted a unified mechanistic explanation, with different physical phenomena implicated across reaction chemistries, catalyst compositions, and choice of solvent. In this perspective, we describe a general framework for understanding cation effects in electrocatalysis based on electrostatics. We argue that cations influence reaction rates by modifying the strength of the electric field present at the catalyst surface, which alters the energetics of adsorbed intermediates and transition states according to their dipole moments and polarizabilities. The magnitude of this field depends on how cations arrange at the electrode surface, controlled by their size, shape, solvation, and packing efficiency. Cations that can arrange more densely result in a steeper potential drop at the electrode surface and consequently a stronger electric field. Our model further identifies two criteria for observing cation effects: (1) the operating potential must be negative of the electrode’s potential of zero total charge, ensuring that cations accumulate at the interface, and (2) the energetics of the kinetically relevant elementary step must be field sensitive. This framework reconciles previously inconsistent trends, including why cation effects appear only for some catalysts, why reaction selectivity is sensitive to cation identity, and why activity can increase with cation size on certain metals but decrease on others. Supported by kinetic measurements, spectroscopy, and atomistic simulations, the model provides both conceptual value for building intuition about catalysis at charged interfaces and predictive value for anticipating trends for new reactions, catalysts, and electrolytes. We conclude by highlighting the importance of electric fields across electrochemical, thermochemical, and biological catalysis and propose that considering the electrostatic environment around active sites offers new opportunities for improving activity and selectivity.

cation effects↗
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2026

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗