Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Computational Modeling and Experimental Characterization of Martensitic Transformations in Nicoal for Self-Sensing Materials

Fundamental changes to aero-vehicle management require the utilization of automated health monitoring of vehicle structural components. A novel method is the use of self-sensing materials, which contain embedded sensory particles (SP). SPs are micron-sized pieces of shape-memory alloy that undergo transformation when the local strain reaches a prescribed threshold. The transformation is a result of a spontaneous rearrangement of the atoms in the crystal lattice under intensified stress near damaged locations, generating acoustic waves of a specific spectrum that can be detected by a suitably placed sensor. The sensitivity of the method depends on the strength of the emitted signal and its propagation through the material. To study the transition behavior of the sensory particle inside a metal matrix under load, a simulation approach based on a coupled atomistic-continuum model is used. The simulation results indicate a strong dependence of the particle's pseudoelastic response on its crystallographic orientation with respect to the loading direction and suggest possible ways of optimizing particle sensitivity. The technology of embedded sensory particles will serve as the key element in an autonomous structural health monitoring system that will constantly monitor for damage initiation in service, which will enable quick detection of unforeseen damage initiation in real-time and during onground inspections.

Wallace, T. A.↗
DOE OSTI2025

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗
DOE OSTI2025

A review of displacement cascade simulations using molecular dynamics emphasizing interatomic potentials for TPBAR components

This review explores molecular dynamics simulations for studying radiation damage in Tritium Producing Burnable Absorber Rod (TPBAR) materials, emphasizing the role of interatomic potentials in displacement cascades. Recent machine learning potentials (MLPs), trained on quantum data, enhance prediction accuracy over traditional models like EAM. We highlight temperature, PKA energy, and composition effects on damage evolution in TPBAR components, recommending suitable potentials and discussing advancements for materials in extreme radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-principles investigation of the resistive switching energetics in monolayer MoS 2 : insights into metal diffusion and adsorption

A deeper understanding of resistive switching (RS) in 2D materials is essential for advancing neuromorphic computing. The Dissociation-Diffusion-Adsorption (DDA) model offers a useful framework for probing RS mechanisms in non-volatile memory (NVM) and in-memory computing. We have employed first-principles density functional theory (DFT) to explore dissociation, diffusion, and adsorption phenomena within the DDA model, focusing on the interactions between exemplary metal atoms (Au, Ag, Cu) and monolayer MoS 2 . Nudged elastic band (NEB) calculations evaluated diffusion barriers in pristine and sulfur-vacancy MoS 2 . Charged systems were modeled to assess the impact of applied bias on migration pathways. We also examined metal dissociation from bulk electrodes and adsorption at S vacancies. Ag/MoS 2 shows the lowest dissociation barrier (~0.034 eV), while Au and Cu exhibit similar values (~0.32 eV). These insights highlight Ag as a promising candidate for low-energy RS applications and provide guidance for optimizing switching efficiency in 2D memory devices.

Atomistic models↗
DOE OSTI2026

Ion correlations explain kinetic selectivity in diffusion-limited solid-state synthesis reactions

Establishing viable solid-state synthesis pathways for novel inorganic materials remains a major challenge in materials science. Previous pathway design methods using pairwise reaction approaches have navigated the thermodynamic landscape with first-principles data but lack kinetic information, limiting their effectiveness. This gap leads to suboptimal precursor selection and predictions, especially for reactions forming competing phases with similar formation energies, where ion diffusion is a critical influence. Here we demonstrate an inorganic synthesis framework by incorporating machine learning-derived transport properties through ‘liquid-like’ product layers into a thermodynamic cellular reaction model. In the Ba–Ti–O system, known for its competitive polymorphism, we obtain accurate predictions of phase formation with varying BaO:TiO2 ratios as a function of time and temperature. We find that diffusion–thermodynamics interplay governs phase compositions, with cross-ion transport coefficients critical for predicting diffusion-limited selectivity. This work bridges length scales and timescales by integrating solid-state reaction kinetics with first-principles thermodynamics and spatial reactivity.

Atomistic models↗
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↗