Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

In‐Operando Spatiotemporal Imaging of Coupled Film‐Substrate Elastodynamics During an Insulator‐to‐Metal Transition

Abstract The drive toward non‐von Neumann device architectures has led to an intense focus on insulator‐to‐metal (IMT) and the converse metal‐to‐insulator (MIT) transitions. Studies of electric field‐driven IMT in the prototypical VO 2 thin‐film channel devices are largely focused on the electrical and elastic responses of the films, but the response of the corresponding TiO 2 substrate is often overlooked, since it is nominally expected to be electrically passive and elastically rigid. Here, in‐operando spatiotemporal imaging of the coupled elastodynamics using X‐ray diffraction microscopy of a VO 2 film channel device on TiO 2 substrate reveals two new surprises. First, the film channel bulges during the IMT, the opposite of the expected shrinking in the film undergoing IMT. Second, a microns thick proximal layer in the substrate also coherently bulges accompanying the IMT in the film, which is completely unexpected. Phase‐field simulations of coupled IMT, oxygen vacancy electronic dynamics, and electronic carrier diffusion incorporating thermal and strain effects suggest that the observed elastodynamics can be explained by the known naturally occurring oxygen vacancies that rapidly ionize (and deionize) in concert with the IMT (MIT). Fast electrical‐triggering of the IMT via ionizing defects and an active “IMT‐like” substrate layer are critical aspects to consider in device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dissecting Disorder: Defect-Driven Structural Complexity in Layered Li3InCl6 Solid Electrolyte

Halide solid electrolytes have emerged as promising candidates for solid-state batteries owing to their high oxidative stability and ionic conductivity. Among them, Li3InCl6 (LIC) has attracted significant attention. However, diffraction patterns of LIC synthesized via different methods exhibit distinct differences particularly at low-angle reflectionsindicative of underlying structural disorder. These variations are attributed to deviations from ideal crystallographic order, especially stacking faults, whose impact on structure and ion transport remains poorly understood. Here, we identify and quantify stacking faults in LIC samples prepared under different synthetic conditions. Using X-ray diffraction and time-of-flight neutron diffraction, we construct and refine stacking fault models that accurately reproduce the experimental diffraction features. LIC samples with higher degrees of stacking faults exhibit only negligible differences in ionic conductivities and activation energies. This indicates that stacking faults have a limited impact on altering the Li+ diffusion pathway along the c-axis, likely due to the high concentration of vacancies in the In layers, while Li+ diffusion remains nearly unchanged in the ab-plane. Our results account for the observed differences in diffraction patterns across samples and provide a quantitative assessment of faulting probabilities and stacking sequences. The insights gained from this study are expected to be broadly applicable to other layered halide solid electrolytes and contribute to a deeper understanding of the role of structural disorder in ion transport.

Liu, Jue [ORNL] (ORCID:000000024453910X)↗
DOE OSTI2024

First principles density functional theory study of tritium species adsorption on Ni(111) surface and diffusion in nickel-sublayer for tritium storage

The nickel-plated zircaloy-4 is used as a tritium ( 3 H) getter in the tritium-producing burnable absorber rods (TPBARs) to capture 3 H produced in the 6 Li-riched annular γ-LiAlO 2 pellet under neutron irradiation. The experimental data and our previous theoretical results showed that the 3 H species produced from the γ-LiAlO 2 pellet were mainly 3 H 2 and 3 H 2 O. These 3 H species diffuse from the surface of the LiAlO 2 pellet across vacuum to the nickel-plated zircaloy-4 getter and then further diffuse into the getter to chemically form metal hydrides. While a number of studies show that oxygen binds strongly as compared to 3 H on the nickel (Ni) layer, the detailed mechanism of 3 H species absorption and diffusion across the Ni plate and Ni/Zr interface are still unclear. By employing density functional theory calculations, here we explored the 3 H 2 and 3 H 2 O species adsorption and dissociation on the Ni(111) surface and diffusion into the Ni sublayer. Our results indicated that the 3 H 2 and 3 H 2 O dissociate on the Ni(111) surface. The NiO x and Ni(O 3 H) x could be formed in the Ni layer due to the higher oxygen (O) diffusion energy barrier and formation of Ni vacancy defects. The oxygen was found to be retained in the Ni layer from diffusing across the Ni–Zr interface. This was revealed by comparing the diffusion barriers for 3 H with O. 3 H was found to have nearly three times smaller diffusion barrier than for O, making 3 H comparatively easier to diffuse through the Ni layer. In conclusion, the obtained results provide guidelines for experimental measurements on 3 H retention behavior in TPBARs and may open further avenues to explore the impurity effects on 3 H diffusion and storage at the Ni/zircaloy interfaces.

Tafen, De Nyago [National Energy Technology Lab. (↗
DOE OSTI2024

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗
DOE OSTI2023

Analysis of position-dependent cavity parameters in irradiated metals to obtain insight on fundamental defect migration phenomena

Motion of point defects is a fundamental process that governs microstructure and properties of materials. Here, we examine the near-surface and grain boundary cavity swelling depth profiles in neutron- and ion-irradiated simple metals (Cu, Ni, and Fe-Cr) and investigate diffusional broadening of implanted Ni ions in Fe-Cr alloys. Vacancy migration energies and radiation-enhanced diffusion were experimentally estimated. Cavity denuded zone widths near planar sinks are shown to be dependent on temperature, damage rate, balance of point defects and sinks, and vacancy migration energies. An enhanced cavity swelling zone adjacent to the void-denuded zone was observed in specimens with low to moderate sink strength and interpreted as evidence of 1D gliding interstitial clusters. Radiation-enhanced diffusional broadening of implanted Ni ions (at 400–550 °C) in Fe and Fe-Cr was up to 250 nm toward the surface, which leads to a much broader near-peak-damage region where cavity swelling is suppressed. Diffusional broadening of implanted ions is calculated to be similarly pronounced for self-ion irradiations, particularly near the peak and higher swelling temperature regimes. Adequately high ion energies (~8–15 MeV) are recommended for ion-irradiation studies to provide a sufficiently broad midrange safe analysis region with minimized surface and injected ion effects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhancing low-temperature sintering in the MgO-LiF system: Mechanistic insights

In the present article, we provide compelling evidence that minor (1 wt%) additions of micron and nanometre-sized LiF particles in MgO decompose leading to free Li diffusing into MgO surfaces enhancing vacancy production during direct current sintering. The addition of nanometre-sized LiF particles leads to a > 500 °C reduction in the sintering temperature and over 99 % theoretical density of final consolidated compacts. Correlating differential scanning calorimetry with in-situ x-ray diffraction, together with Schottky calculations, post sintering electron microscopy and laser induced breakdown spectroscopy, we uncover critical insights into this impressive reduction in sintering temperature. Our quantitative analysis reveals that MgO and LiF interact at low temperatures with the diffusion of Li into the surface of MgO particles due to the intrinsic structural disorder of the LiF and MgO crystallites. Nanometre-sized LiF particles were found to react the most at low temperatures due to their increased structural disorder. Our multimodal characterization points to a Li-promoted densification and sintering mechanism. This impressive reduction in sintering temperature can be harnessed to promote low-temperature fabrication of MgO-based composites for technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The mobility of small, over-pressurized helium bubbles in tungsten at 2000 K

Fusion reactor environments inevitably lead to the formation of high-pressure helium bubbles whose nucleation, growth, and diffusion strongly impact the performance of plasma-facing components. Here this research describes a diffusion mechanism of over-pressurized bubbles via a sequence of Frenkel pair nucleation, self-interstitial migration, and Frenkel pair annihilation. Molecular dynamics was used to simulate the diffusion of small bubbles in tungsten at 2000 K with helium-per-vacancy ratios in the range of 4.5 to 7. The diffusion coefficients are calculated and their dependence on helium content, number of vacancies, and number of attached self-interstitials is characterized. It is found that bubbles are most mobile when the nucleation/annihilation rates of Frenkel pairs are nearly equal and when the bubbles nucleate and annihilate a single self-interstitial. All bubbles experience a peak diffusivity, which can be as high as 10 -11 m 2 /s decreasing with bubble size. The calculated diffusion coefficients provide valuable insight into the mobility of small, high-pressure bubbles, and can be used as input parameters in mesoscale models to improve predictions of plasma-surface interactions. (LA-UR-21-21881)

36 MATERIALS SCIENCE↗
DOE OSTI2023

One dimensional wormhole corrosion in metals

Corrosion is a ubiquitous failure mode of materials. Often, the progression of localized corrosion is accompanied by the evolution of porosity in materials previously reported to be either three-dimensional or two-dimensional. However, using new tools and analysis techniques, we have realized that a more localized form of corrosion, which we call 1D wormhole corrosion, has previously been miscategorized in some situations. Using electron tomography, we show multiple examples of this 1D and percolating morphology. To understand the origin of this mechanism in a Ni-Cr alloy corroded by molten salt, we combined energy-filtered four-dimensional scanning transmission electron microscopy and ab initio density functional theory calculations to develop a vacancy mapping method with nanometer-resolution, identifying a remarkably high vacancy concentration in the diffusion-induced grain boundary migration zone, up to 100 times the equilibrium value at the melting point. Deciphering the origins of 1D corrosion is an important step towards designing structural materials with enhanced corrosion resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Insights into oxygen diffusion in rare earth disilicate environmental barrier coatings

Environmental barrier coatings (EBC) are crucial for the use of SiC-based ceramic matrix composites in high-temperature combustion environments, yet knowledge of oxygen diffusion in these coatings is limited. This study investigates oxygen diffusion dynamics in the β-RE 2 Si 2 O 7 system to minimize oxygen penetration in rare earth disilicates. We analyze defect formation energy under varying oxygen conditions, identifying key diffusion mechanisms. In oxygen-rich environments, the most favorable neutral interstitial oxygen diffuses along the [110] direction. In oxygen-poor conditions, neutral oxygen vacancies rotate around Y and Si atoms, exhibiting a diffusivity of 6.59×10 −22 m 2 /s at 1500 K for β-Yb 2 Si 2 O 7 . Under intermediate oxygen levels, charged interstitial oxygen diffuses via concerted interstitialcy along the [001] direction with a diffusivity of 6.21×10 −17 m 2 /s. Additionally, alloying rare earth Y with Er and Yb increases diffusion barriers, contributing to improved EBC performance in extreme environments. The insights gained provides valuable guidance for designing robust coatings tailored to withstand extreme operational environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Active and passive defects in tetragonal tungsten bronze relaxor ferroelectrics

Tetragonal tungsten bronze (TTB) based oxides constitute a large family of dielectric materials which are known to exhibit complex distortions producing incommensurately modulated superstructures as well as significant local deviations from their average symmetry. The local deviations produce diffuse scattering in diffraction experiments. The structure as well as the charge dynamics of these materials are anticipated to be sensitive to defects, such as cation or oxygen vacancies. In this work, in an effort to understand how the structural and charge dynamical properties respond to these two types of vacancy defects, we have performed measurements of dielectric susceptibilities and single crystal diffraction experiments of two types of TTB materials with both ‘filled’ (Ba 2 NdFeNb 4 O 15 and Ba 2 PrFeNb 4 O 15 ) and ‘unfilled’ (Sr 0.5 Ba 0.5 Nb 2 O 6 ) cation sublattices. We also perform these measurements before and after oxygen annealing, which alters the oxygen vacancy concentrations. Surprisingly, we find that many of the diffuse scattering features that are present in the unfilled structure are also present in the filled structure, suggesting that the random fields and disorder that are characteristic of the unfilled structure are not responsible for many of the local structural features that are reflected in the diffuse scattering. Furthermore, oxygen annealing clearly affected both color and dielectric properties, consistent with a diminishment of the oxygen vacancy concentration, but had little effect on observed diffuse patterns.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2024

Investigation of Low-Energy Lattice Dynamics and Their Role in Superionic Na Diffusion and Ultralow Thermal Conductivity of Na 3 PSe 4 as a Solid-State Electrolyte

The atomic dynamics of Na 3 PSe 4 were investigated using a combination of neutron scattering experiments and ab initio and machine-learned molecular dynamics simulations to probe the interplay of fast ionic diffusion with atomic vibrations (phonons) of the host lattice. Our results reveal the existence of low-energy vibrational modes, simultaneously involving motions of Na + ions and framework polyanion subunits, and show that these modes become strongly overdamped in the superionic regime as they couple with the Na + hopping process. In particular, the Na+ migration energy landscape is strongly impacted by low-energy phonons derived from a soft acoustic branch of the host lattice, which modulates the diameter of the Na + diffusion channel at the bottleneck. We find that an additional factor for the enhanced Na + conductivity in Na 3 PSe 4 is the presence of Na-vacancies, which also affect the low-frequency dynamics and thermal vibration amplitudes, pointing to an interplay between Na + vacancies and host dynamics, jointly enhancing ionic diffusivity. Finally, we investigate the origin of ultralow thermal conductivities in Na 3 PSe 4 and Na 3 PS 4 using Green–Kubo simulations and find that low-energy acoustic phonon modes of the overall crystal framework provide a dominant contribution to the thermal conductivity.

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

The Sluggish Diffusion of Cations in CeO 2 Probed through Molecular Dynamics and Metadynamics Simulations

Cation diffusion in fluorite‐structured CeO 2 , though far slower than anion diffusion, is an important, high‐temperature process because it governs diverse fabrication and degradation phenomena. Herein, cation diffusion is studied by means of classical molecular dynamics and metadynamics simulations. Three different mechanisms are examined: migration involving an isolated cerium vacancy, migration involving a cerium vacancy in a defect associate with an oxygen vacancy, and migration involving a cation divacancy. For each mechanism, defect diffusion coefficients are calculated as a function of temperature, from which the respective activation enthalpy of defect migration is obtained. Through comparisons with experimental cation diffusion data (specifically, of the absolute magnitude of the cation diffusivity as well as its activation enthalpy), it is concluded that cation diffusion takes place predominantly neither by isolated vacancies nor by cation vacancy–oxygen vacancy associates but by cation divacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion Behavior of Oversized Fission Products in bcc Fe Cladding: A First-Principles Study

Fuel-Cladding Chemical Interaction (FCCI) poses significant challenges in nuclear reactors, where fission products from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of oversized fission products, Pr, Nd, Ce, and La, within bcc Fe cladding using density functional theory (DFT), nudged elastic band (NEB) method, and self-consistent mean field (SCMF) theory. Our results reveal significant long-range vacancy binding energies, particularly up to the 6th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The NEB calculations indicate significant high barriers for the dissociation of 1nn vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products was derived in Arrhenius form. The significant trapping effect of vacancies by a very dilute amount of fission products reduces vacancy mobility, leading to an oversaturation of point defects, void nucleation, and swelling. These are critical issues for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. This study provides essential insights for developing advanced cladding materials and design strategies to mitigate FCCI, ultimately enhancing nuclear reactor safety and performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Mechanism of Metal Intercalation under Graphene through Small Vacancy Defects

Metal intercalation under graphene has attracted extensive experimental and theoretical research because of its capability to manipulate the electronic structure and properties of graphene. However, the pathways and mechanisms of intercalation are still not well understood. Here, we systematically investigate the intercalation process of metal atoms through graphene vacancies using first-principles calculations. We show that the energy barrier for metal atom penetration through the vacancies in graphene is small as long as the size of the vacancy is larger than a mono-vacancy. However, metal atoms are strongly bonded to the vacancy so that the detachment energy of a metal atom from the vacancy is extremely high. This inhibits the diffusion of the metal atom into the gallery beneath the surface to complete the intercalation process. On the other hand, our calculation results show that the detachment energy of a metal atom from a metal dimer at small vacancy defects is significantly reduced, making intercalation much easier. Therefore, the key step limiting the intercalation process is the detachment of the metal atoms from vacancy defects. This finding from our study provides useful insight into the defect-assisted intercalation mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗