Search NASA⌕ Search

SEARCH · Search NASA

Results for “vacancy diffusion”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Vacancy diffusion barrier spectrum and diffusion correlation in multicomponent alloys

Vacancy diffusion serves a crucial role in many important kinetic behaviors and properties of multicomponent alloys. Essential questions, however, persist regarding how chemical complexity affects diffusion and what unique characteristics, if any, set these alloys apart from traditional metals. Using neural network kinetics model, we study vacancy diffusion in NbMoTa alloy across a broad temperature range (2600 to 800 K). Unlike pure metals, the two key diffusion parameters—diffusion correlation factor f and activation energy ΔG m —are not constant in alloys, but instead substantially decrease with decreasing temperature. This temperature dependence arises from a reduced number of active vacancy jump pathways at lower temperatures, leading to more correlated diffusion. Upon examining vacancy diffusion throughout the entire compositional space of the Nb-Mo-Ta system, we discover that the slowest vacancy diffusion surprisingly occurs in the non-equimolar region, rather than the equimolar concentration where the configurational entropy is highest. The diffusion barrier spectrum, characterizing the diffusion energy landscape, is an intrinsic material characteristic, which controls both f and ΔG m and, thereby, the diffusivity. Lastly, we find that the vacancy diffusion rate drops noticeably in the presence of local chemical order in the NbMoTa system, particularly for MoTa alloys with long-range B2 order.

Defects↗
DOE OSTI2025

Probing the Impact of Vacancy Diffusion on Void Dynamics at the Lithium Metal–Solid Electrolyte Interface

Lithium (Li) metal-based solid-state batteries (SSBs) are considered promising candidates for next-generation energy storage due to their superior energy density and enhanced safety compared to conventional Li-ion systems. However, their practical application is limited by challenges such as void formation at the Li-solid electrolyte (SE) interface, which disrupts ion transport and accelerates interfacial degradation. This work investigates how the coupled effects of electro-dissolution kinetics and surface diffusion at the Li metal surface govern the evolution of interfacial morphology during stripping. This work examines the influence of three distinct surface diffusion modes, which are terrace diffusion, step diffusion, and interlayer diffusion, on maintaining interfacial stability. In addition, how the dominant surface diffusion mechanism can overcome the contact loss due to high reaction kinetics is explored. Furthermore, the roughness of the Li metal anode surface is quantified, and the influence of different diffusion mechanisms on the evolution of the dynamic solid–solid interface is examined. The critical role of temperature in enhancing Li surface diffusivity and expanding the regime of stable contact is highlighted. By identifying distinct regimes of interface stability, this study analyzes how non-uniform electrochemical dynamics dictate void morphology evolution and interfacial contact. These insights offer guiding principles for engineering robust Li–SE interfaces in SSBs.

lithium metal anode↗
DOE OSTI2022

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Larger than uncorrelated vacancy diffusion contributions in chemically disordered crystalline materials

Arising from a variational approach to compute diffusion coefficients, we introduce “superkinetic kinosons”, which are contributions to the diffusion flux of mobile defects and atomic species by single jump mechanisms that can exceed their uncorrelated values. In vacancy-mediated diffusion in crystalline random and high entropy materials, these contributions primarily arise from intermittent jumps of slow moving atoms that remove the vacancy out of correlation traps and enable a quantification of such effects on the diffusion coefficients. They can be significant contributors to diffusion, even when faster moving atoms have infinite percolation networks. Furthermore, their existence and importance provide fundamental insights about underlying aspects of vacancy diffusion such as formation of localized correlation traps and deviation of vacancy diffusion from percolation behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

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

Vacancy-Dependent Diffusion Mechanism in Oxygen-Defective SrFeO 3 Perovskite Materials: First-Principles Density Functional Theory and Experimental Approach

Understanding oxygen diffusion at the atomic scale in SrFeO 3−δ perovskites is crucial for developing oxygen storage materials with optimal performance. Such materials are required to have high stability, corrosion resistance, and acceptable oxygen storage capacity at moderate operating temperatures and pressures. Here, in this study, we used first-principles density functional theory and thermogravimetric analysis to study the vacancy-dependent oxygen diffusion in oxygen-deficient SrFeO 3−δ (δ = 0, 0.065, 0.125, 0.25, 0.5) perovskites. The electronic structures, including the partial- and spin-resolved density of states, for different SrFeO 3−δ phases were calculated and compared with available experimental and theoretical results. By mapping the migration pathways, we investigated diffusion mechanisms and calculated the energy barriers for oxygen diffusion in cubic, orthorhombic, and brownmillerite phases of SrFeO 3−δ perovskites. Using the calculated energy barriers, we deduced the diffusion time scales and diffusion coefficients within SrFeO 3−δ . A diffusion coefficient on the order of 10 –8 m 2 /s was obtained for SrFeO 2.875 . We experimentally investigated the roles of temperature and oxygen partial pressures on the redox kinetics and deduced the kinetics rate and diffusion density, which agreed well with the calculated values for the density of diffusing oxygen vacancy in the lattice. Our results showed that the energy barrier tends to reduce at higher oxygen concentrations. Our results serve as an important guideline for designing oxygen storage materials with optimal redox kinetics.

chemical looping with oxygen uncoupling (CLOU)↗
DOE OSTI2024

Effect of Carbon on Void Nucleation in Iron

The study reports the significance of carbon presence in affecting void nucleation in Fe. Without carbon, void nucleation rates decrease gradually at high temperatures but remain significantly high and almost saturated at low temperatures. With carbon present, even at 1 atomic parts per million, void nucleation rates show a low-temperature cutoff. With higher carbon levels, the nucleation temperature window becomes narrower, the maximum nucleation rate becomes lower, and the temperature of maximum void nucleation shifts to a higher temperature. Fundamentally, this is caused by the change in effective vacancy diffusivity due to the formation of carbon-vacancy complexes. The high sensitivity of void nucleation to carbon comes from the high sensitivity of void nucleation to the vacancy arrival rate in a void. The void nucleation is calculated by first obtaining the effective vacancy diffusivity considering the carbon effect, then calculating the defect concentration and defect flux change considering both carbon effects and pre-existing dislocations, and finally calculating the void nucleation rate based on the recently corrected homogeneous void nucleation theory. The study is important not only in the fundamental understanding of impurity effects in ion/neutron irradiation but also in alloy engineering for judiciously introducing impurities to increase swelling resistance, as well as in the development of simulation and modeling methodologies applicable to other metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Ultrafast Oxygen Conduction in Sillén Oxychlorides

Oxygen ion conductors are crucial for enhancing the efficiency of various clean energy technologies, including fuel cells, solid oxide air batteries, electrolyzers, membranes, sensors, and more. In this study, a structure-similarity analysis of ≈62k oxygen-containing compounds identified the MBi 2 O 4 X (M = rare-earth element, X = halogen element) family as promising candidates for fast oxygen transport. Among these, LaBi 2 O 4 Cl is found as an ultrafast oxygen conductor with an ultralow migration barrier of 0.1 eV based on ab initio studies. Its 2D layered structure, featuring a “triple fluorite” layer, supports diffusion of both oxygen vacancies and interstitials. In addition to vacancy diffusion with a 0.1 eV barrier, ab initio studies show interstitial diffusion exhibits a modest barrier of 0.6–0.8 eV. Frenkel pairs are found to be the dominant defects in intrinsic LaBi 2 O 4 Cl, facilitating significant vacancy-mediated oxygen diffusion at elevated temperatures. With 2.8% oxygen vacancies, LaBi 2 O 4 Cl is predicted to achieve a conductivity of 0.3 S/cm at 25 °C in a single crystal. Experimental synthesis and characterization of polycrystalline LaBi 2 O 4 Cl and Sr-doped LaBi 2 O 4 Cl revealed conductivity exceeding that of YSZ and LSGM below 400 °C, with lower activation energies, achieving a total conductivity of 0.1−0.2 mS/cm at 300 °C. Here, while these results confirm its potential of fast oxygen transport, we suggest further experimental optimization of LaBi 2 O 4 Cl, including aliovalent doping and microstructure refinement, could significantly enhance its performance, facilitating fast oxygen conduction approaching room temperature.

Defects↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2022

A coupled vacancy diffusion-dislocation dynamics model for the climb-glide motion of jogged screw dislocations

Here, we develop a novel model to study the climb/glide motion of jogged screw dislocations within the discrete dislocation dynamics (DDD) framework. We present results for the dependence of the climb velocity on the applied stress and on the jog size and distribution statistics. We show that the model predictions are consistent with experimental data in both γ-TiAl and Zircaloy-4. The ranges of the applied stress and jog spacing that determine the dominance of one of three dislocation mechanisms are identified. These are the jog dragging, dipole dragging, and dipole bypass mechanisms, respectively. The overall dislocation motion in the jog dragging regime is composed of glide of screw segments and climb of jogs, controlling the plastic strain and the creep rate, respectively. Based on current simulations and on a detailed examination of the predicted jog heights compared to experiments, we advance the hypothesis that a combination of jog dragging and dipole bypass mechanisms is necessary to reproduce the high creep rate observed in some experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accelerated molecular dynamics simulations of dislocation climb in nickel

We report the mechanical behavior of materials operating under high temperatures is strongly influenced by creep mechanisms such as dislocation climb, which is controlled by the diffusion of vacancies. However, atomistic simulations of these mechanisms have traditionally been impractical due to the long time scales required. To overcome these time scale challenges, we use Parallel Trajectory Splicing (ParSplice), an accelerated molecular dynamics method, to simulate dislocation climb in nickel. We focus on modeling the activity of a vacancy near a jog on an edge dislocation in order to observe vacancy pipe diffusion and vacancy absorption at the jog. From rigorously constructed trajectories encompassing more than 2000 vacancy absorption events over a simulation time of more than 4μs at 900 K, a comprehensive sampling of available atomistic mechanisms is collated and analyzed further with molecular statics calculations. We estimate average rates for pipe diffusion and vacancy absorption into the jog using data from the dynamic and static calculations, finding very good agreement. Our results strongly suggest that the dominant mechanism for vacancy absorption by jogs is via biased diffusion to the dislocation core followed by fast pipe diffusion to the jog.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular dynamics simulations of fission gas xenon (Xe) diffusion at UO 2 grain-boundaries (Rev.1)

The diffusivity of fission gas xenon (Xe) at UO 2 grain-boundaries is one of the most important parameters in mechanistic modeling of fission gas diffusion in UO 2 based nuclear fuels. In this report, we use molecular dynamics simulations to investigate the Xe diffusivity in UO 2 grain-boundaries, employing the many-body potential developed by Cooper, Rushton and Grimes for UO 2 . Three different types of grain-boundaries are investigated, twist Σ5, tilt Σ5, and a random grain-boundary. Diffusion activation energies in the range of 0.39 – 1.46 eV are obtained for the Xe diffusivity. Comparison to results for the uranium vacancy diffusivity from MD simulations employing the same methodology suggests a weak to moderate attractive Xe-uranium vacancy binding energy depending on the grain-boundary type.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Identifying Hidden Li–Si–O Phases for Lithium‐Ion Batteries via First‐Principle Thermodynamic Calculations

SiO–based materials are promising alloys and conversion‐type anode materials for lithium‐ion batteries and are recently found to be excellent dendrite‐proof layers for lithium‐metal batteries. However, only a small fraction of the Li–Si–O compositional space has been reported, significantly impeding the understanding of the phase transition mechanisms and the rational design of these materials both as anodes and as protection layers for lithium‐metal anodes. Herein, we identify three new thermodynamically stable phases within the Li–Si–O ternary system (Li 2 SiO 5 , Li 4 SiO 6, and Li 4 SiO 8 ) in addition to the existing records via first‐principle calculations. The electronic structure simulation shows that Li 2 SiO 5 and Li 4 SiO 8 phases are metallic in nature, ensuring high electronic conductivity required as electrodes. Moduli calculations demonstrate that the mechanical strength of Li–Si–O phases is much higher than that of lithium metal. The diffusion barriers of interstitial Li range from 0.1 to 0.6 eV and the interstitial Li hopping serves as the dominating diffusion mechanism in the Li–Si–O ternary systems compared with vacancy diffusion. These findings provide a new strategy for future discovery of improved alloying anodes for lithium‐ion batteries and offer important insight towards the understanding of the phase transformation mechanism of alloy‐type protection layers on lithium‐metal anodes.

Qu, Jiale↗
DOE OSTI2021

Persistent and partially mobile oxygen vacancies in Li-rich layered oxides

Increasing the energy density of layered oxide battery electrodes is challenging as accessing high states of delithiation often triggers voltage degradation and oxygen release. Here we utilize transmission-based X-ray absorption spectromicroscopy and ptychography on mechanically cross-sectioned Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ electrodes to quantitatively profile the oxygen deficiency over cycling at the nanoscale. The oxygen deficiency penetrates into the bulk of individual primary particles (~200 nm) and is well-described by oxygen vacancy diffusion. Using an array of characterization techniques, we demonstrate that, surprisingly, bulk oxygen vacancies that persist within the native layered phase are indeed responsible for the observed spectroscopic changes. We additionally show that the arrangement of primary particles within secondary particles (~5 μm) causes considerable heterogeneity in the extent of oxygen release between primary particles. Finally, our work merges an ensemble of length-spanning characterization methods and informs promising approaches to mitigate the deleterious effects of oxygen release in lithium-ion battery electrodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Oxygen vacancy migration and impact on high voltage DC polarization in 0.8BaTiO 3 –0.2BiZn 0.5 Ti 0.5 O 3

Electrical polarization and defect transport are examined in 0.8BaTiO 3 –0.2BiZn 0.5 Ti 0.5 O 3 , an attractive capacitor material for high power electronics. Oxygen vacancies are suggested to be the majority charge carrier at or below 250°C with a grain conduction hopping activation energy of 0.97 eV and 0.92 eV for thermally stimulated depolarization current (TSDC) and impedance spectroscopy measurements, respectively. At higher temperature, thermally generated electronic conduction with an activation energy of 1.6 eV is dominant. Significant oxygen vacancy concentration is indicated (up to ~1%) due to cation vacancy formation (i.e., acceptor defects) from observed Bi (and likely Zn) volatility. Oxygen vacancy diffusivity is estimated to be 10 -12.8 cm 2 /s at 250°C. Low diffusivity and high activation energies are indicative of significant defect interactions. Dipolar oxygen vacancy defects are also indicated, with an activation energy of 0.59 eV from TSDC measurements. In conclusion, the large oxygen vacancy content leads to a short lifetime during high voltage (30 kV/cm), high temperature (250°C) direct current (DC) electrical measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Kinetic pathways of fast lithium transport in solid electrolyte interphases with discrete inorganic components

The transport of lithium ions in the solid electrolyte interphase (SEI) has been previously accepted to proceed in two steps: a fast pore diffusion through the outer, porous organic layer followed by a slow knock-off or vacancy diffusion in the inner, dense inorganic layer. The second step is believed to be the rate-limiting step during fast-charging. In this study, we have intentionally constructed a thicker SEI (SEI-rich) structure on the surface of monoclinic Nb 2 O 5 (H-Nb 2 O 5 ) by adding LiNO 3 into a conventional ethylene carbonate based electrolyte. The electrochemical performance of two electrodes, one SEI-rich and one with few SEI (SEI-lean), was found to be almost the same, including their fast-charging capability and cycling stability, despite the significant difference in their SEI structure. Importantly, analysis using cryogenic scanning/transmission electron microscopy showed the discrete decoration of individual inorganic particles (e.g., Li 2 O) and amorphous species (LiN x O y /organic components) over the surface of H-Nb 2 O 5 . These discrete inorganic particles are in contradiction to the formation of dense inner inorganic layer, which has been commonly postulated. Based on these findings, we propose a new mechanism for Li ion transport through the SEI: one-step pore diffusion, without the second step slow diffusion. This one-step pore diffusion process provides an extremely fast Li ion transport, and effectively removes the kinetic limitation of Li ion transport in the SEI for fast charging. Finally, these results strongly suggest that the influence of SEI structure on the transport kinetics of lithium ions is much less significant than previously accepted. These results offer a new understanding of possible lithium ion transport pathway within SEI and may have implications for the future designs of fast-charging battery materials.

25 ENERGY STORAGE↗