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

Oxidation of Silicon Particles Suspended in Mullite Matrix With Dual-Mode Vacancy Diffusion of Oxygen

Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.

silicon oxidation↗

Minimizing the diffusivity difference between vacancies and interstitials in multi-principal element alloys

Interstitial atoms usually diffuse much faster than vacancies, which is often the root cause for the ineffective recombination of point defects in metals under irradiation. Here, via ab initio modeling of single-defect diffusion behavior in the equiatomic NiCoCrFe(Pd) alloy, we demonstrate an alloy design strategy that can reduce the diffusivity difference between the two types of point defects. The two diffusivities become almost equal after substituting the NiCoCrFe base alloy with Pd. The underlying mechanism is that Pd, with a much larger atomic size (hence larger compressibility) than the rest of the constituents, not only heightens the activation energy barrier ( E a ) for interstitial motion by narrowing the diffusion channels but simultaneously also reduces E a for vacancies due to less energy penalty required for bond length change between the initial and the saddle states. Our findings have a broad implication that the dynamics of point defects can be manipulated by taking advantage of the atomic size disparity, to facilitate point-defect annihilation that suppresses void formation and swelling, thereby improving radiation tolerance.

36 MATERIALS SCIENCE↗

Grain boundary self-diffusion and point defect interactions in α -U via molecular dynamics

Though metallic U-Zr fuel has been used in nuclear reactors since the 1960s, many of its fundamental and thermodynamic properties are still unknown. The a-U phase, which has a highly anisotropic crystal structure and physical properties, is present in U-Zr fuel. The character and behavior of a-U grain boundaries will strongly impact fuel thermophysical performance under irradiation. Here, we study the interaction of point defects with grain boundaries, diffusion along grain boundaries, and the predicted diffusional creep behavior of a-U via molecular dynamics. We calculate the segregation energy of vacancies and interstitials to grain boundaries and quantify the biased sink strength of the grain boundaries, and observe that this sink strength is not strongly dependent on the grain boundary orientation. We also find that grain boundary diffusivity is strongly dependent on the grain boundary energy and grain boundary orientation. The presence of point defects within the grain boundary can induce diffusion in grain boundaries with low formation energies and can enhance diffusion in high-energy grain boundaries. We also find that diffusional creep of a-U at prototypical metallic fuel operation conditions is extremely high and could help explain observed metallic fuel swelling behaviors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

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↗

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↗

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↗

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↗

Lithium-ion drifting: Application to the study of point defects in floating-zone silicon

The use of lithium-ion (Li(+)) drifting to study the properties of point defects in p-type Floating-Zone (FZ) silicon crystals is reported. The Li(+) drift technique is used to detect the presence of vacancy-related defects (D defects) in certain p-type FZ silicon crystals. SUPREM-IV modeling suggests that the silicon point defect diffusivities are considerably higher than those commonly accepted, but are in reasonable agreement with values recently proposed. These results demonstrate the utility of Li(+) drifting in the study of silicon point defect properties in p-type FZ crystals. Finally, a straightforward measurement of the Li(+) compensation depth is shown to yield estimates of the vacancy-related defect concentration in p-type FZ crystals.

Silicon↗

Effect of Crystal Defects on Minority Carrier Diffusion Length in 6H SiC Measured Using the Electron Beam Induced Current Method

We report values of minority carrier diffusion length in n-type 6H SiC measured using a planar Electron Beam Induced Current (EBIC) method. Values of hole diffusion length in defect free regions of n-type 6H SiC, with a doping concentration of 1.7El7 1/cu cm, ranged from 1.46 microns to 0.68 microns. We next introduce a novel variation of the planar method used above. This 'planar mapping' technique measured diffusion length along a linescan creating a map of diffusion length versus position. This map is then overlaid onto the EBIC image of the corresponding linescan, allowing direct visualization of the effect of defects on minority carrier diffusion length. Measurements of the above n-type 6H SiC resulted in values of hole diffusion length ranging from 1.2 micron in defect free regions to below 0.1 gm at the center of large defects. In addition, measurements on p-type 6H SiC resulted in electron diffusion lengths ranging from 1.42 micron to 0.8 micron.

Tabib-Azar, Massood↗

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↗

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↗

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↗

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↗

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↗

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↗