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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microstructure-dependent rate theory model of defect segregation and phase stability in irradiated polycrystalline LiAlO 2

We report gamma lithium aluminate (LiAlO 2 ) is a breeder material for tritium and is one of key components in a tritium-producing burnable absorber rod (TPBAR). Dissolution and precipitation of second phases such as LiAl 5 O 8 and voids are observed in irradiated LiAlO 2 . Such microstructure changes cause the degradation of thermomechanical properties of LiAlO 2 and affect tritium retention and release kinetics, and hence, the TPBAR performance. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed for investigating the accumulation of species and phase stability in polycrystalline LiAlO 2 structures under irradiation. Three sublattices (i.e. [Li, Al, V] I [O, V o ] II [Li i , Al i , O i , V i ] III ), and concentrations of six diffusive species (i.e. Li; vacancy of Li or Al at [Li, Al, V] I sublattice, O vacancy at [O, V o ] II sublattice, and Li, Al and O interstitials at [Li i , Al i , O i , V i ] III interstitial sublattices; are used to describe spatial and temporal distributions of defects and chemistry. Microstructure-dependent thermodynamic and kinetic properties including the generation, reaction, and chemical potentials of defects and defect mobility are taken into account in the model. The parametric studies demonstrated the capability of the developed RIS model to assess the effect of thermodynamic and kinetic properties of defects on the segregation and depletion of species in polycrystalline structures and to explain the phase stability observed in irradiated LiAlO 2 samples. The developed RIS model will be extended to study the precipitation of LiAl 5 O 8 and voids and tritium retention by integrating the phase-field method.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2023

Fast Surface Oxygen Release Kinetics Accelerate Nanoparticle Exsolution in Perovskite Oxides

Exsolution is a recent advancement for fabricating oxide-supported metal nanoparticle catalysts via phase precipitation out of a host oxide. A fundamental understanding and control of the exsolution kinetics are needed to engineer exsolved nanoparticles to obtain higher catalytic activity toward clean energy and fuel conversion. Since oxygen release via oxygen vacancy formation in the host oxide is behind oxide reduction and metal exsolution, we hypothesize that the kinetics of metal exsolution should depend on the kinetics of oxygen release, in addition to the kinetics of metal cation diffusion. In this report we probe the surface exsolution kinetics both experimentally and theoretically using thin-film perovskite SrTi 0.65 Fe 0.35 O 3 (STF) as a model system. We quantitatively demonstrated that in this system the surface oxygen release governs the metal nanoparticle exsolution kinetics. As a result, by increasing the oxygen release rate in STF, either by reducing the sample thickness or by increasing the surface reactivity, one can effectively accelerate the Fe 0 exsolution kinetics. Fast oxygen release kinetics in STF not only shortened the prereduction time prior to the exsolution onset, but also increased the total quantity of exsolved Fe 0 over time, which agrees well with the predictions from our analytical kinetic modeling. The consistency between the results obtained from in situ experiments and analytical modeling provides a predictive capability for tailoring exsolution, and highlights the importance of engineering host oxide surface oxygen release kinetics in designing exsolved nanocatalysts.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2023

Atomistic models of In and Ga diffusion in Cu(In,Ga)Se 2

To improve the performance of Cu(In,Ga)Se 2 thin-film photovoltaic devices, a robust understanding of the dominant diffusion pathways of the alloy species In and Ga is needed. Here, the most probable defect complexes and mechanisms for In and Ga diffusion are identified with the aid of density functional theory. The binding energies and migration barriers for these complexes are calculated in bulk CuInSe 2 and CuGaSe 2 . Analytic models and kinetic lattice Monte Carlo simulations are employed to predict the diffusivity of In and Ga under variations in composition and temperature. Here, we find that a model based on coulombic interactions between group III antisites and vacancies on the Cu-sublattice produces results that match well with experiment.

14 SOLAR ENERGY↗
DOE OSTI2023

Final-stage densification kinetics of direct current–sinteredZrB2

Final-stage sintering was analyzed for nominally phase pure zirconium diboride synthesized by borothermal reduction of high-purity ZrO2. Analysis was conducted on ZrB2 ceramics with relative densities greater than 90% using the Nabarro–Herring stress–directed vacancy diffusion model. Temperatures of 1900°C or above and an applied uniaxial pressure of 50 MPa were required to fully densify ZrB2 ceramics by direct current sintering. Ram travel data were collected and used to determine the relative density of the specimens during sintering. Specimens sintered between 1900 and 2100°C achieved relative densities greater than 97%, whereas specimens sintered below 1900°C failed to reach the final stage of sintering. The average grain size ranged from 1.0 to 14.7 μm. The activation energy was calculated from the slope of an Arrhenius plot that used the Kalish equation. The activation energy was 162 ± 34 kJ/mol, which is consistent with the activation energy for dislocation movement in ZrB2. The diffusion coefficients for dislocation motion that controls densification were 5.1 × 10−6 cm2/s at 1900°C and 5.1 × 10−5 cm2/s at 2100°C, as calculated from activation energy and average grain sizes. This study provides evidence that the dominant mechanism for final-stage sintering of ZrB2 ceramics is dislocation motion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal diffusivity recovery and defect annealing kinetics of self-ion implanted tungsten probed by insitu transient grating spectroscopy

Tungsten is a promising candidate material for plasma-facing armour components in future fusion reactors. A key concern is irradiation-induced degradation of its normally excellent thermal transport properties. In this comprehensive study, thermal diffusivity degradation in ion-implanted tungsten and its evolution from room temperature (RT) to 1073 K is considered. Five samples were exposed to 20 MeV self-ions at RT to achieve damage levels ranging from 3.2 x 10 -4 to 3.2 displacements per atom (dpa). Transient grating spectroscopy with in-situ heating was then used to study thermal diffusivity evolution as a function of temperature. Using a kinetic theory model, an equivalent point defect density is estimated from the measured thermal diffusivity. The results showed a prominent recovery of thermal diffusivity between 450 K and 650 K, which coincides with the onset of mono-vacancy mobility. After 1073 K annealing samples with initial damage of 3.2 x 10 -3 dpa or less recover close to the pristine value of thermal diffusivity. For doses of 3.2 x 10 -2 dpa or higher, on the other hand, a residual reduction in thermal diffusivity remains even after 1073 K annealing. Transmission electron microscopy reveals that this is associated with extended, irradiation-induced dislocation structures that are retained after annealing. A sensitivity analysis shows that thermal diffusivity provides an efficient tool for the assessing total defect content in tungsten up to 1000 K.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2026

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistically-informed modeling of point defect clustering and evolution in irradiated ThO 2

A cluster dynamics (CD) model has been developed to investigate the nucleation and growth of point defect clusters, i.e., interstitial prismatic loops and nanoscale and sub-nanoscale voids, in ThO 2 during irradiation by energetic particles. The model considers cluster off-stoichiometry due to the asymmetry of point defect generation on the O and Th sublattices under irradiation, as well as the point defect diffusivities and the defect binding energies to clusters. The energies were established using detailed molecular dynamics simulations considering the statistical variability of cluster configuration. A high-order adaptive time-integration has been used to solve the model. The predicted loop density and their average size is in good agreement with reported experimental observations for proton irradiated ThO 2 at 600°C. In conclusion, the model did not predict void evolution due to the sluggish kinetics of cation vacancies, explaining the absence of voids in proton irradiated ThO 2 (and other oxides) at relatively low temperatures.

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

Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – a cluster dynamics and phase field model

Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.

Hu, Shenyang↗
DOE OSTI2021

Bulk Diffusion of Cl through O Vacancies in α-Cr 2 O 3 : A Density Functional Theory Study

Many metals form passive oxides which are broken or weakened by anionic-induced degradation leading to material failure and anionic diffusion is an important step of this degradation process. The diffusion of anionic species through oxides involves a combination of diffusion along grain boundaries, cracks, and channels, as well as bulk diffusion via point defects which is the focus of this study. Using density functional theory, we study bulk diffusion of Cl through O vacancies in α -Cr 2 O 3 as a model system for passive metal oxides. Little is known about Cl diffusion and bonding characteristics, so we benchmark our work through comparison to numerous studies on O diffusion in α -Cr 2 O 3 to analyze similarities and differences between the O and Cl diffusion in the passivation and degradation process of α -Cr 2 O 3 respectively. Unlike O diffusion, the lowest diffusion barrier for Cl is cross-plane diffusion between two (0001) planes through a vacant cation site but the much shorter in-plane diffusion path within the same coordination polyhedron has 35% higher barrier. This work provides the basis for considering the contributions of Cl bulk-diffusion in the overall diffusion kinetics of Cl through α -Cr 2 O 3 .

Oware Sarfo, Kofi↗
DOE OSTI2020

Effect of Confinement of SIA Cluster Diffusion by Impurities on Radiation Defect Accumulation due to 14MeV neutrons in Tungsten

Impurities acting as traps, reduce the self-interstitial atom (SIA) cluster mobility, thereby increasing recombination and in turn enhancing radiation-resistance of polycrystalline materials. However, at elevated temperatures, depending on SIA-trap (impurity) binding energy, SIA clusters may trap and detrap (detach) numerous times. In tungsten, SIA clusters of all sizes glide one-dimensionally (1D) along their Burgers vector direction. Hence, detrapped SIAs would retrace their path, confining their 1D-glide between traps. However, small SIA clusters can change the direction of their 1D-glide by overcoming a rotation barrier. In addition, impurities or solutes are known to reduce this barrier. A lower rotation barrier can result in detrapped SIA clusters that diffuse in a random 1D direction, effectively leading to 3D (net-3D) diffusion. Here we investigate the effect of two cases of SIA diffusion, namely confined-1D diffusion and net-3D diffusion, on the damage accumulation in polycrystalline tungsten under 14 MeV neutron irradiation using the object kinetic Monte Carlo method. Simulations were performed using cascade debris obtained from molecular dynamics simulations with primary knock-on atom (PKA) energies corresponding to 14 MeV neutrons. In the simulations with net-3D diffusion, SIA clusters of all sizes diffuse in a random 1D direction after the clusters detrap. While for the confined 1D diffusion case, SIA clusters larger than size 5 retain their original direction of 1D-glide. In both types of simulations, trapped vacancies are assumed to be permanently immobilized. Unexpectedly, the damage accumulation was lower with confined 1D diffusion than that with net-3D. We present a systematic comparison of the influence of the two cases of SIA diffusion on the radiation damage accumulation as a function of dose rate, detrapping barrier, and impurity concentration at 1025 K.

Accumulation and recovery of lattice damage, tungs↗
DOE OSTI2020

Point defect evolution under irradiation: Finite size effects and spatio-temporal correlations

The evolution of point defect concentrations under irradiation is controlled by their diffusion properties, and by their formation and elimination mechanisms. The latter includes the mutual recombination of vacancies and interstitials, and the elimination of point defects at sinks. Two models are traditionally used to predict the defect concentration evolution, the standard rate theory (SRT) and the atomistic kinetic Monte Carlo (AKMC). We show in this work a large discrepancy in the defect concentrations between both methods when the average number of defects in the AKMC simulation box is close or lower than one. The reason is that AKMC naturally captures strong space and time correlations between vacancies and interstitials generated by the finite size of the periodic simulation box. These correlations strongly affect the recombination rate and the point defect concentrations. SRT fails to predict such correlations, and the corresponding solution deviates from the more accurate solution given by the AKMC simulation under similar conditions. These finite size correlation effects are strong when the elimination of point defects occurs by recombination only, but can still be significant in the presence of sinks. In order to account for the spacio-temporal correlation in a continuum framework, we introduce a Correlated Pair Theory (CPT). This theory fully takes into account the correlations between vacancy and interstitial pairs and predicts point defect concentrations in good agreement with AKMC simulations. Inversely, if the goal is for the AKMC to reproduce defect concentrations in bulk using small simulation boxes, the naturally occurring correlations need to be corrected. We show here that the CPT can be used to modify the elimination rates in the AKMC simulations, so as to yield point defect concentrations in agreement with SRT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2024

Intermittent oxidation kinetics and metal/oxide interfacial undulation

The phenomenon of oxygen-adsorption-induced surface restructuring is widespread across various metal-oxygen systems, yet its impact on initiating bulk oxide formation remains largely unexplored. In this study, through in-situ atomic-resolution electron microscopy observations of surface oxidation of Cu(110) and Cu 85 Au 15 (110), we unveil intermittent oxide-film growth modulated by oxygen-induced surface restructuring. This modulation is evidenced by repeated pinning of Cu 2 O growth front at isolated Cu columns of the c(6×2)-O reconstruction, owing to required long-range diffusion of Cu and O atoms to the Cu 2 O growth front. We reveal that Cu vacancies, generated at the Cu 2 O growth front, are injected into the Cu 2 O/Cu interface, inducing hill-and-valley undulation of the Cu 2 O film. In contrast, atomic vacancies produced during the Cu 85 Au 15 (110) oxidation preferentially migrate into interfaces between Au-rich and Au-poor regions in the bulk, resulting in a flat and adherent Cu 2 O film. These findings demonstrate the critical role of oxygen-induced surface restructuring in modulating oxide-film-growth kinetics and manipulability of the fate of injected vacancies by alloying, thereby offering insights applicable to a broader range of metal-oxygen systems for fine-tuning oxidation kinetics and enhancing oxide/metal interfacial adhesion.

25 ENERGY STORAGE↗
DOE OSTI2024

Phase transformation mechanism in irradiation-induced superlattice formation

Atomic kinetic Monte Carlo simulations were used to model void superlattice formation under irradiation in molybdenum, driven by anisotropic diffusion of self-interstitial atoms. A change in the phase transformation mechanism from nucleation and growth to spinodal decomposition occurred with increasing dose rate, with both mechanisms leading to superlattice formation. Analysis of a rate-theory based analytical model showed that an observed change in the kinetics of vacancy accumulation, the appearance of a region of positive second derivative in the plot of average vacancy concentration versus time, was caused by the onset of spinodal instability. Further, the analytical model showed that for molybdenum and several other metals where void superlattice formation is commonly observed, the phase transformation likely occurs by nucleation and growth. However, nickel may offer the possibility of experimental observation of the transition between phase transformation mechanisms.

36 MATERIALS SCIENCE↗