Search NASA⌕ Search

SEARCH · Search NASA

Results for “vacancy diffusion kinetics”

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.

100 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

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

Mechanistic Understanding of Methane Combustion over Ni/CeO 2 : A Combined Experimental and Theoretical Approach

Catalytic oxidation of methane (CH 4 ) over nonprecious Ni/CeO 2 catalysts has received a lot of attention due to the large natural gas reserves found in North America and the prohibitive cost of palladium-based catalysts, commonly used for CH 4 oxidation. However, the catalytic mechanism of CH 4 oxidation over Ni/CeO 2 still remains unclear. Moreover, the parameters affecting the reaction rates, the interaction between nickel and CeO 2 , and the reaction intermediates are still not well understood. In this study, kinetic model fitting, CH 4 temperature-programmed reduction-mass spectroscopy (CH 4 TPR-MS), in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and density functional theory (DFT) calculations were combined to elucidate the mechanism of complete oxidation of CH 4 over Ni/CeO 2 . CH 4 TPR-MS showed that the complete oxidation of CH 4 over Ni/CeO 2 requires 55–120 °C lower compared to bare CeO 2 or Ni/quartz sand; complete oxidation of CH 4 took place when the surface oxygen species were abundant, while partial oxidation products (CO, H 2 ) were formed when the oxygen species were depleted. In situ DRIFTS showed that CH 3 , CH 2 , CO, and CO 2 were formed after CH 4 activation over Ni/CeO 2 , while CH 3 O species were not observed. Combining those findings with kinetic model fitting, a redox Mars–van Krevelen (MvK) mechanism showed the best description of the experimental observations. The MvK mechanism involves the reaction of dissociated oxygen species with gas-phase CH 4 while water inhibits the reaction rate by adsorbing on the oxidized sites. Moreover, CH 4 activation leads to the reduction of the active sites and oxygen vacancy formation followed by reoxidation of the active sites by gas-phase O 2 . A CH 4 oxidation reaction pathway over Ni/CeO 2 is proposed by DFT calculations. In summary, the findings shown here suggest that CH 4 oxidation over Ni/CeO 2 follows a redox MvK mechanism and provides guidance for the rational design of non-precious-metal catalysts for CH 4 oxidation reactions.

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

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2021

Maintaining a Flat Li Surface during the Li Stripping Process via Interface Design

Electroplating has been the main focus in mitigating the dendrite growth on the Li-metal electrode; however, the stripping process is equally critical, since the nonsmooth Li surface during stripping will lead to nonuniform local current density, planting the seeds for dendrite growth. In this paper, density functional theory (DFT) and kinetic Monte Carlo (KMC) techniques were combined to investigate the vacancy evolution in Li interfaced with different solid–electrolyte interphase (SEI) materials. It was found that the lithiophilic interface, such as Li/Li 2 O, repels vacancies into the bulk Li, so Li atoms can quickly fill the Li vacancies near the Li/Li 2 O interface and maintain a smooth Li surface. In contrast, the lithiophobic interface, such as Li/LiF, traps Li vacancies toward the interface, and the accumulated Li vacancies form voids and roughen the surface. The predicted critical stripping current density, below which a smooth Li surface will be maintained, is therefore much faster at the lithiophilic interface than that at the lithiophobic interface. It was further revealed that the lithiophilicity at different SEI or coating materials can be ranked as Li/Li 2 O > Li/LiPON > Li/Li 2 CO 3 > Li/LiF based on the calculated interfacial adhesion and accumulation of electron density at the interface. Furthermore, this suggests that interface and coating design at nanoscale can be effective for maintaining a smooth Li surface during the stripping process, solving another challenge to achieving a dendrite-free Li-metal electrode in both liquid and solid electrolytes.

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

Coupling between bulk thermal defects and surface segregation dynamics

Surface segregation is a phenomenon that depends on the delicate interplay between thermodynamic driving forces and kinetic obstacles, for which elevated temperature is often needed to enhance the atom mobility and reach equilibrium. Using the classic system of Cu 3 Au ( 100 ) under nonisothermal conditions, in this study, we show an adatom process underlying transient surface segregation dynamics through the temperature-change-driven creation and annihilation of thermal vacancies in the bulk and the resulting bulk-surface mass exchanges. This is demonstrated by monitoring the surface composition evolution of Cu 3 Au ( 100 ) with temperature changes between 250 and 500 ° C , showing that the increase in temperature decreases monotonically the surface Au concentration as a result of the transfer of more Cu than Au from the bulk to the surface to form Cu-rich clusters of adatoms. Such a bulk thermal defect effect is expected to be universal in inducing the disparity in the bulk-surface mass exchanges of dissimilar atoms in multicomponent materials because of the inherent differences in the vacancy formation energies of the constituent atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ba 1−x Sr x FeO 3−δ as an improved oxygen storage material for chemical looping air separation: a computational and experimental study

Chemical looping air separation (CLAS) is a promising technology to generate oxygen-rich gas streams to enable efficient carbon dioxide capture during fossil fuel combustion or gasification. CLAS relies on the capture and release of oxygen from the atmosphere using the redox properties of an oxygen-selective solid oxide carrier. This study investigates the redox characteristics of Ba 1−x Sr x FeO 3−δ (0.0 ≤ x ≤ 0.417, 0.0 ≤ δ ≤ 0.5) using a combination of density functional theory (DFT) calculations and experimental verification using X-ray diffraction, thermogravimetric analysis, and oxygen-temperature-programmed desorption. The DFT computed energies of the Ba 1−x Sr x FeO 3−δ perovskites reveal a composition-dependent transition from hexagonal to cubic phases as the Sr-concentration or oxygen vacancy concentration increases. Oxygen vacancy formation energies of the cubic perovskites are found to be lower than those of their hexagonal counterparts. A low oxygen diffusion barrier of ∼1 eV combined with the thermodynamic preference of Ba 1−x Sr x FeO 3−δ compositions that form in a cubic phase suggests them as promising candidates for oxygen storage applications. The experimental results corroborate this finding by identifying Ba 0.75 Sr 0.25 FeO 3−δ in the cubic phase as an optimal composition offering low-temperature oxygen storage capacities comparable to that of the state-of-the-art Sr 0.75 Ca 0.25 FeO 3−δ perovskite oxygen storage material at 325 °C and 350 °C.

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

Migration Kinetics of Surface Ions in Oxygen-Deficient Perovskite During Topotactic Transitions

Oxygen diffusivity and surface exchange kinetics underpin the ionic, electronic, and catalytic functionalities of complex multivalent oxides. Towards understanding and controlling the kinetics of oxygen transport in emerging technologies, it is highly desirable to reveal the underlying lattice dynamics and ionic activities related to oxygen variation. In this study, the evolution of oxygen content is identified in real-time during the progress of a topotactic phase transition in La 0.7 Sr 0.3 MnO 3-δ epitaxial thin films, both at the surface and throughout the bulk. Using polarized neutron reflectometry, a quantitative depth profile of the oxygen content gradient is achieved, which, alongside atomic-resolution scanning transmission electron microscopy, uniquely reveals the formation of a novel structural phase near the surface. Surface-sensitive X-ray spectroscopies further confirm a significant change of the electronic structure accompanying the transition. The anisotropic features of this novel phase enable a distinct oxygen diffusion pathway in contrast to conventional observation of oxygen motion at moderate temperatures. The results provide insights furthering the design of solid oxygen ion conductors within the framework of topotactic phase transitions.

36 MATERIALS SCIENCE↗