Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “defect 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Kinetic Monte Carlo simulations of aging in δ -Pu

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

Actinides↗
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 OSTI2023

Resonant interactions involving local vibrational modes in crystals

When an impurity with a light mass is inserted into a crystal, it can undergo a high-frequency oscillation referred to as a local vibrational mode (LVM). A Fermi resonance may occur between the LVM and lower-frequency modes of the defect. The LVM may also interact with phonons or the electromagnetic field. Understanding these interactions can help model and control diffusion, defect reactions, and thermal conductivity. LVMs have been probed in semiconductors using pressure and alloying as experimental parameters, resulting in anticrossing between localized and extended vibrational modes. Furthermore, these types of vibrational interactions could play an important role in the stability and thermal properties of organic–inorganic hybrid semiconductors. The coupling between an LVM and electromagnetic wave yields an “LVM polariton,” an excitation that has significant vibrational and electric-field amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Extended defects-enhanced oxygen diffusion in ThO 2

Oxygen self-diffusion is key to understanding stoichiometry and defect structures in oxide nuclear fuels. Experimentally, low activation-barrier oxygen migration was found in ThO 2 , a candidate nuclear fuel, possibly due to short-circuit diffusion mechanisms. Here, we perform extensive molecular dynamics simulations to show that various types of extended defects can enhance oxygen self-diffusion with a much-reduced activation barrier in ThO 2 . In this work, we consider extended defects including 1D (dislocation), 2D (grain boundary), and 3D (void) defects. Due to the distinct characteristics of each type of extended defect, the modulation of oxygen diffusion varies. Further, these results provide a quantitative description of oxygen transport, which is significantly enhanced within a close distance (nanometer scale) from the extended defects. Among all these considered defects, grain boundary, particularly the low-energy $Σ$3 twin boundary, exhibits the strongest effect on increasing oxygen transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Automated calculation and convergence of defect transport tensors

Defect diffusion is a key process in materials science and catalysis, but as migration mechanisms are often too complex to enumerate a priori, calculation of transport tensors typically have no measure of convergence and require significant end-user intervention. These two bottlenecks prevent high-throughput implementations essential to propagate model-form uncertainty from interatomic interactions to predictive simulations. In order to address these issues, we extend a massively parallel accelerated sampling scheme, autonomously controlled by Bayesian estimators of statewide sampling completeness, to build atomistic kinetic Monte Carlo models on a state-space irreducible under exchange and space group symmetries. Focusing on isolated defects, we derive analytic expressions for drift and diffusion coefficients, providing a convergence metric by calculating the Kullback–Leibler divergence across the ensemble of diffusion processes consistent with the sampling uncertainty. The autonomy and efficacy of the method is demonstrated on surface trimers in tungsten and Hexa-interstitials in magnesium oxide, both of which exhibit complex, correlated migration mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coalescence dynamics of prismatic dislocation loops due to vacancy supersaturation

Recovery processes in metals with a high defect concentration are at the heart of a structural alloy's usability in engineering applications. However, such processes take place via nonconservative mechanisms mediated by point defect diffusion and are governed by their interactions with the underlying microstructure. Indeed, one of the main challenges in studying recovery processes is their intrinsic multiscale complexity, which often negates the direct use of experimental techniques. Here, in this Research Letter, we study the mechanism of dislocation loop coalescence in metals by simulating the interaction dynamics between vacancies and prismatic loops using a stochastic framework that seamlessly bridges point defect kinetics with dislocation mechanics. Our simulations reveal two fundamental discoveries: first, that climb processes are governed by local fluctuations that are not captured by continuum elastodiffusion models, and, second, that prismatic loops create internal traps for vacancies, leading to the formation of rings around the perimeter of the loops that affect how coalescence takes place. Our results show an excellent agreement with in situ experimental observations of loop coalescence in Fe under irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimizing ensemble NV − spin properties of fluorescent diamond microparticles by systematic low pressure high temperature annealing

Low pressure high temperature annealing is a means for driving nitrogen and defect diffusion in diamond to reduce internal lattice damage without the need for technically complicated high-pressure cells. Herein, we perform a systematic time (5, 15, and 30 min) and temperature (1200 °C–1800 °C) study of effects of low-pressure high temperature annealing on photoluminescence, spin concentrations, and spin relaxation properties of NV centers in ca. 3 μm synthetic type 1b diamond particles. Annealing in the temperature range of ca. 1400 °C–1700 °C for even 5 min leads to a higher optically detected magnetic resonance contrast as compared to standard annealing at 900 °C for 2 h. Particles annealed at 1700 °C for 5 min exhibit a contrast close to about 13% as compared to about 9% for those annealed at 900 °C for 2 h. A reduction in the zero-field splitting strain parameter from E ≈ 4.5 MHz to ≈ 2.5 MHz and spectral linewidth from Δν ≈ 7 MHz to ≈ 4 MHz are observed even after 5 min annealing at 1700 °C. Improvements in these spectral parameters resulted in a roughly 2-fold reduction in the noise level of temperature monitoring experiment utilizing an ensemble of NV centers in the particles. Annealing in the temperature range of 1600 °C for 15 or 30 min or 1700 °C for 5 min resulted in NV T 1 relaxation times approaching ca. 5 ms typically observed for bulk diamond. Quantitative electron paramagnetic resonance (EPR) allowed for estimations of thermal activation energies of paramagnetic center annihilation. Monitoring the primary defect concentration (P1 and other defects with half integer spins) and utilizing second order kinetic modeling, an activation energy of 3.63 ± 0.28 eV was estimated. Alternatively, using the NV half field EPR signal and first order kinetic modeling, a similar activation energy 3.89 ± 0.29 eV was estimated.

NV centers↗
DOE OSTI2025

Scanning tunneling microscope tip-induced formation of Bi bilayers on Bi 2 Te 3

We report the formation of Bi(111) bilayer islands and crater structures on Bi 2 Te 3 (111) surfaces induced by voltage pulses from a scanning tunneling microscope tip. Pulses above a threshold voltage (+3 V) produce craters ∼0.5μm in diameter, similar to the size of the tip. Redeposited material self-assembles into a network of atomically ordered islands with a lattice constant identical to the underlying Bi2Te3 surface. The island size monotonically decreases over several micrometers from the pulse site, until the pristine Bi 2 Te 3 surface is recovered. We assign these islands to Bi bilayer based on atomic resolution images, analysis of step heights, and tunneling spectroscopy. Here, the dependence of bilayer formation on bias polarity and the evidence for defect diffusion together suggest a mechanism driven by the interplay of field evaporation and tunneling-current-induced Joule heating.

Bi bilayer↗
DOE OSTI2023

Does photoinduced structural change by β-rhombohedral boron occur at room temperature?

β-rhombohedral boron (β-B) is an intriguing solid among the elements. The crystal structure of β-B is widely accepted as a complex structure based on the B 12 clusters, and with a lot of intrinsic structural defects. In the last decade, β-B was proved to be the only elemental semiconductor that exhibits the self-compensation property accompanied by changing structural defects during atom doping. Moreover, it is reported that various physical properties of β-B change discontinuously around 550 K, which may originate from defect diffusion [S. Hoffmann and H. Werheit, Solid State Sciences 14, 1572 (2012)]. We previously reported that a volume expansion of about 1.5% takes place above this temperature (between 650 K and 850 K). On the other hand, this expanded structure is stable in dark conditions for more than 4.5 months. Finally, light irradiation at room temperature reverses the volume back to the initial state, and a possible scenario will be discussed in this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating Diiodomethane-Induced Improvement in Photonically Cured MAPbI 3 Solar Cells

Diiodomethane (CH 2 I 2 ) has been reported to improve the photonically cured perovskite solar cell (PSC) device performance by reducing pinholes. Here we show that even for pinhole-free methylammonium lead iodide (MAPbI 3 ) made by photonic curing, adding CH 2 I 2 promotes significant performance improvement. We elucidate the mechanisms behind the observed improvement. In addition to current-density versus voltage measurements, we perform a wide variety of characterization to compare the crystallinity, grain structure, optical, chemical, and electrical properties of the photonically cured samples with and without CH 2 I 2 in the MAPbI 3 layer to those of thermally annealed devices. Here, the addition of CH 2 I 2 promotes grain growth in the vertical direction, increases the I/Pb ratio in the final film, and removes the I - ionic diffusion and defect signature associated with iodine interstitials. As a result, we achieve a champion efficiency of 15.04% with MAPbI 3 conversion time of 20 ms, comparable to PSCs that have MAPbI 3 layer thermally annealed for 10 min. Understanding the mechanisms behind additive-induced improvements for non-thermal annealing processes is critical to enabling high-speed PSC manufacturing.

14 SOLAR ENERGY↗
DOE OSTI2025

Effects of dose rate on the microstructure and deuterium retention in γ-LiAlO 2 pellets

This report presents experimental findings obtained from November 2024 to September 2025. Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by two competing processes: defect production, driven by dose rate, and defect recovery, controlled by defect diffusion, interaction, and annihilation. At a given dose, the resulting microstructural evolution depends on both the dose rate and irradiation temperature. As a continuation of our FY24 tritium science project, which investigated temperature effects at a fixed dose rate, this study focuses on dose-rate effects at a fixed temperature to provide deeper insights into the irradiated microstructure and compositional changes in γ-LiAlO 2 pellets. The study aims to elucidate the impact of dose rate on microstructure, precipitate morphology, deuterium retention, and lithium volatility in γ-LiAlO 2 pellets under sequential 120 keV He + and 80 keV D 2 + ion irradiation. Three dose rates of 7.3×10 4 , 2.9×10 4 and 6.8×10 5 dpa/s were applied to achieve the same total ion fluence of 2×10 17 (He + +D + )/cm 2 at 500 °C, corresponding to a maximum combined dose of 7.55 dpa at ~255 nm. The irradiated pellets were subsequently characterized using scanning transmission electron microscopy (STEM) and time-of-flight secondary ion mass spectrometry (ToF-SIMS). The microstructural response of γ-LiAlO 2 to ion irradiation was found to be strongly dose-rate dependent. At medium and high dose rates, irradiation produced a surface amorphized layer and an underlying crystalline layer containing LiAl 5 O 8 precipitates, with implanted gases accumulating and forming blisters at the crystalline-amorphous interface. It remains to be investigated whether the amorphized layer was produced by He + ion irradiation prior to D 2 + ion irradiation. In contrast, at low dose rates, the material remained crystalline, with cavities, likely gas-filled, distributed around precipitates, within the γ-LiAlO 2 matrix, and along grain boundaries. While precipitate morphology exhibited anisotropy, their size showed little sensitivity to dose rate in the applied range of this study. This result, however, does not rule out the possibility that further reductions in dose rate could influence precipitate size. High dose-rate irradiation enhanced protonium–deuterium isotopic exchange and lithium depletion in the amorphized region. Collectively, the results show that dose rate governs amorphization, gas redistribution, isotopic exchange, and lithium depletion, providing important insights into the mechanisms underlying structural evolution in γ-LiAlO 2 pellets under reactor-relevant irradiation conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2023

The role of native defects and their diffusivity in CdTe

CdTe is a leading thin-film technology for solar-cell applications. Basic features such as the optimum band gap of 1.5 eV and simple manufacturing process have been driving the research on improving the current record efficiency of 22% for CdTe solar cells toward the theoretical limit of 32%. Low p-type doping and short carrier lifetimes in the CdTe absorber have been suggested as the main limiting factors and attributed to defect compensation and carrier recombination. Understanding the role of point defects and impurities in CdTe is crucial to solving these problems. Despite many years of research, the sources of compensation and the microscopic recombination mechanisms are still unclear. By using hybrid density functional calculations with spin-orbit coupling and large supercells, we investigate the electronic properties, the formation energies, and the diffusivities of the native defects in CdTe. As possible sources of compensation for p-type conductivity, we discuss the structural and electronic properties of Cd interstitial, Te vacancies, and TeCd antisites. The stability of each defect is discussed in terms of both the defect formation energy and the migration barrier. We also address possible mechanisms that lead to non-radiative recombination.

14 SOLAR ENERGY↗
DOE OSTI2021

Molecular dynamics studies of lattice defect effects on tritium diffusion in zirconium

Tritium diffusion in α-Zr containing point defects such as vacancies or self-interstitial atoms (SIAs) is simulated using molecular dynamics. Point defects rapidly aggregate to form extended defects, such as 3D nanoclusters and Frank loops. The geometry of extended defects is affected by the presence of tritium. At low temperature and in the absence of tritium, vacancies aggregate to form stacking fault pyramids. Addition of tritium at these temperatures promotes aggregation of vacancies to form 3D nanoclusters, within which the tritium concentration can be sufficiently high to suggest that these defects may serve as nucleation sites for hydride precipitation. Trapping of tritium in vacancy nanocluster reduces the calculated bulk diffusivity by an amount proportional to the vacancy concentration. At high temperature, vacancy clusters change shape to form planar basal dislocation loops, which bind tritium less strongly, leading to a sharp reduction in the fraction of trapped tritium and a corresponding increase in tritium diffusivity at high temperature. In contrast, SIAs increase tritium diffusion through α-Zr. Analysis of atomic trajectories shows that tritium does not interact directly with SIAs. In conclusion, diffusion enhancement is instead related to expansion of the lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Combining machine-learned and empirical force fields with the parareal algorithm: application to the diffusion of atomistic defects

We numerically investigate an adaptive version of the parareal algorithm in the context of molecular dynamics. This adaptive variant has been originally introduced in [1]. We focus here on test cases of physical interest where the dynamics of the system is modelled by the Langevin equation and is simulated using the molecular dynamics software LAMMPS. In this work, the parareal algorithm uses a family of machine-learning spectral neighbor analysis potentials (SNAP) as fine, reference, potentials and embedded-atom method potentials (EAM) as coarse potentials. We consider a self-interstitial atom in a tungsten lattice and compute the average residence time of the system in metastable states. Our numerical results demonstrate significant computational gains using the adaptive parareal algorithm in comparison to a sequential integration of the Langevin dynamics. We also identify a large regime of numerical parameters for which statistical accuracy is reached without being a consequence of trajectorial accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructure-Dependent Rate Theory Model of Radiation-Induced Segregation in Binary Alloys

Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.

36 MATERIALS SCIENCE↗