Results for “diffuse” — Page 27

DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗

DOE OSTI2025

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗

DOE OSTI2026

Non-classical ion-indiffusion processes in quasi-one-dimensional ionic conductor KTiOPO4

Developments in coercive field engineering and waveguide writing via ion-exchange in the ferroelectric ionic-conductor KTiOPO4 (KTP) have led to breakthroughs within the field of nonlinear optics and quantum photonics. However, the microscopic diffusion dynamics of ion-exchange in KTP are poorly understood and show signs of non-classical diffusion dynamics. The precise control of ion-exchange gradients is critical for developing KTP into a mature, integrated nonlinear optical material but requires a level of control that, in turn, demands a deeper understanding of the ion-exchange process. We employ direct ToF-SIMS mapping to untangle the complex diffusion dynamics of the commonly used ion-exchange dopant ions: Rb, Ba, and K. We map periodic diffusion gradients of two exchanges, one dominated by Rb-ions and the other by Ba-ions. For the first time, we show the interdependent nature of the Ba- and Rb ion in the ion-exchange diffusion gradient. The low incorporation of the Ba ion is mapped and quantified, and dopant–crystal interactions are identified. This work highlights the complex diffusion dynamics of dopants in the KTP structure and lays the foundation for deterministic control of diffusion gradients for the next generation of integrated nonlinear and quantum photonic devices.

Barrett, Laura [ORNL]↗

DOE OSTI2024

Synthesis, Processing, and Use of Isotopically Enriched Epitaxial Oxide Thin Films

Isotopic engineering has emerged as a key approach to study the nucleation, diffusion, phase transitions, and reactions of materials at an atomic level. It aims to uncover mass transport pathways, kinetics, and operational and failure mechanisms of functional materials and devices. Understanding these phenomena leads to deeper insights into important physical processes, such as the transport of ions in energy conversion and storage devices and the role of active sites and supports during heterogeneous catalytic reactions. Likewise, isotopic engineering is being pursued as a means of modifying functionality to enable future technological applications. In this Account, we summarize our recent work employing isotope labeling (e.g., 18 O 2 and 57 Fe) during thin film synthesis and postgrowth processing to reveal growth mechanisms, defect chemistry, and elemental diffusion under working and extreme conditions. Isotope-resolved analysis techniques with nanometer-scale spatial resolution, such as time-of-flight secondary ion mass spectrometry and atom probe tomography, facilitate the accurate quantification of isotopic placement and concentration in our well-defined heterostructures with precisely positioned, isotope-enriched layers. By measuring the nanometer-scale redistribution between natural abundance and isotopically enriched oxygen layers during the deposition of Fe 2 O 3 and Cr 2 O 3 by molecular beam epitaxy, we identified intermixing processes driven by surface adatoms occurring both at the film growth surface and within the first few layers below the surface. Further insights into synthesis mechanisms were gained by studying the tungsten oxide thin films grown by evaporating WO 3 powder in the presence of background 18 O 2 , revealing minimal incorporation of background oxygen during the film formation process. Thermal and radiation-enhanced diffusion in epitaxial Fe and Cr oxides were precisely tracked using 18 O and 57 Fe tracer layers incorporated into model epitaxial oxide thin films. This approach has allowed us to access thermal diffusion behavior at lower temperatures than previously measured, revealing a potential changeover in diffusion mechanism. Understanding radiation-enhanced diffusion in model oxides that represent the surface layers on the structural components of nuclear reactors informs our understanding of their corrosion behavior under irradiation. Isotopic labeling can also provide unique insights into the surface exchange reactions and defect chemistry of electrocatalysts. For instance, tracking the change in 18 O concentration at the surface of an epitaxial LaNiO 3 thin film after the electrocatalytic oxygen evolution reaction revealed the participation of lattice oxygen, confirming a hypothesis that had been proposed previously. Lastly, we highlight a new direction wherein we perform in situ processing studies utilizing isotopic tracers in conjunction with model epitaxial thin films within the atom probe tomography instrument. Additionally, this Account illustrates the great potential of isotopic engineering to enable fundamental mechanistic insights into physical processes and engineer functional properties in epitaxial films, heterostructures, and superlattices.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Toward a Better Understanding of Ni Coarsening in Solid Oxide Cells: NiH on Ni (111) Examined Using a Combined Theoretical Approach

The coarsening of the Ni particles in the hydrogen electrode of solid oxide cells (SOCs) is an important degradation mechanism. Here, in this paper, density-functional theory and kinetic Monte Carlo methods are used to explore our recent hypothesis that the surface diffusion of NiH may cause faster Ni coarsening in electrolysis cell mode under an overpotential. Using both methods, the diffusion constant or diffusivity of NiH on Ni (111) is determined as the product of the surface coverage and single-molecule diffusivity for the first time considering all possible diffusion paths. It is then determined versus overpotential at the triple-phase boundary of the hydrogen electrode assuming a typical operating temperature of the SOC. Under a significant overpotential, the diffusivity of NiH is found to be sufficiently large to support the above hypothesis that NiH may promote Ni coarsening. However, based on the adsorption configurations identified, the dissociation and reformation of NiH on Ni (111) could occur. Thus, more work is needed to develop a model of Ni coarsening considering both molecular and dissociated forms of NiH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Comparison of greenhouse gas emission estimates from six hydropower reservoirs using modeling versus field surveys

As with most aquatic ecosystems, reservoirs play an important role in the global carbon (C) cycle and emit greenhouse gases (GHG) as carbon dioxide (CO 2 ) and methane (CH 4 ). However, GHG emissions from reservoirs are poorly quantified, especially in temperate systems, resulting in high uncertainty. We compared reservoir C emission estimates and uncertainty of diffusive, ebullitive, and degassing pathways in six hydropower reservoirs in the southeastern United States among four data sources: two field-based surveys and two models (including the GHG Reservoir “G-res” Tool). We found that CH 4 diffusion was most similar across data sources (modeled minus observed, bias = - 21 g CO 2-eq m -2 y -1 ) and had low relative uncertainty (coefficient of variation, CV = 0.98). On the other hand, CO 2 diffusion was least consistent across data sources (bias = - 518 g CO 2-eq m -2 y -1 ). Both field surveys indicated strong negative CO 2 diffusion (i.e., CO 2 uptake) at all reservoirs, while G-res estimated positive CO 2 diffusion. By extension, total C emissions showed similar discrepancies, leading to high uncertainty in upscaling and interpreting reservoir source-sink dynamics. Finally, CH 4 ebullition had the highest relative uncertainty (CV = 2.77) due to high variability across sites. We discuss limitations of field surveys and these models, including temperature-based annualization methods, varying definitions of ebullition zones, low sampling resolution, and lack of dynamism. Future field efforts focused on capturing variability in CO 2 diffusion and CH 4 ebullition will be especially valuable in reducing uncertainty and improving models to advance our understanding reservoir GHG emissions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Understanding coarsening of a post-corrosion microstructure in a molten salt by combining phase-field modeling and in situ tomography

Alloys corroding in molten salt have been observed to form bicontinuous, nanoporous microstructures via dealloying, which subsequently undergo coarsening due to facile transport in high-temperature conditions. In this work, we describe a methodology to elucidate the underlying transport mechanisms during coarsening of a bicontinuous microstructure via quantitative comparisons between phase-field simulations and four-dimensional in situ experiments, in this case X-ray nanotomography of the coarsening of a dealloyed 80 wt% Ni-20 wt% Cr microwire in molten KCl-MgCl 2 at 800°C. We conduct phase-field simulations initialized from experimental data to model coarsening via three different transport mechanisms: surface diffusion, solid bulk diffusion, and liquid bulk diffusion. These simulations reproduce key features of the experiment, such as the densification of the outer layer of the dealloyed wire and the reduction in radius over time. We quantitatively compare different microstructural characteristics between the simulations and experiment and extract temporal scaling factors that optimally match the time scales of the simulations to that of the experiment. This allows us to evaluate morphological similarity between the simulations and experiment and relate the experimental coarsening kinetics to fundamental material properties. We find that surface diffusion is most likely to be the dominant coarsening mechanism, and its kinetics imply a surface diffusivity of D S = 8.9 x 10 -20 m 3 /s, which is within the range of reported values for Ni-vacuum interfaces at 800°C. However, the difference between the experiment and the surface diffusion simulation increases substantially at late times, suggesting that other mechanisms, such as the dissolution of residual Cr, may be at play.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗

DOE OSTI2024

Anion Selectivities in Zwitterion Grafted Nanopores: Effect of Zwitterion Architecture

The separation of ions of similar charge is a crucial challenge in many applications, from water treatment to precious metal recovery. Membranes with cross-linked zwitterionic amphiphilic copolymer (ZAC-X) selective layers, which feature self-assembled, zwitterion-lined nanodomains for permeation, offer unique permselectivity between monovalent anions (e.g. $\mathrm{Cl^-}$/$\mathrm{F^-}$). This has motivated studies on the mechanisms of transport and selectivity in this family of materials. In this study, we conducted molecular dynamics simulations of aqueous salt solutions within zwitterion-functionalized nanopores to elucidate the influence of dipole orientation of the ZI ligands on anion diffusivities, partitioning, and permeabilities. Here, our model compares systems with contrasting ZI organization: surface-cation-anion ($\mathrm{S\text{-}ZI^{+}\text{-}ZI^{-}}$, Motif A) and surface-anion-cation ($\mathrm{S\text{-}ZI^{-}\text{-}ZI^{+}}$, Motif B). Our results reveal that Motif A exhibits less pronounced ion pairing due to a spatial separation in the radial profiles of cations and anions. Motif B demonstrates prominent ion pairing for smaller anions owing to their overlap with cation distributions. Further, our potential of mean force profiles reveal that anion partitioning increases with anion size in both ligand motifs, whereas Motif B exhibits significantly higher partitioning selectivity towards larger anions compared to Motif A. Our results for ion diffusivities show that the self-diffusivities of both anions and cations are lower for Motif B compared to Motif A. Such trends in anion partitioning and diffusivities can be explained by differences in the interactions and steric hindrance experienced by the anionic species in Motifs A and B. Finally, our results for anion permselectivity, obtained by combining partitioning and diffusivity, indicate that partitioning trends dominate over diffusivity trends. Consequently, anion permeability increases with anion size, and ligand Motif B yields much higher permselectivity towards larger anions compared to ligand Motif A.

anions↗

DOE OSTI2024

Interpenetrating 3D Electrodes for High-Rate Alkaline Water Splitting

Mass transfer is critical for the reaction kinetics and efficiency of alkaline water splitting (AWS). For AWS to operate at high current densities (hundreds of mA/cm 2 ), the device architecture must ensure a large catalytic surface area, rapid ion diffusion, and minimal solution and charge transfer resistances. Effective electrodes should also facilitate gas bubble detachment and release. 3D-printed electrodes have shown promise, but stacking them increases the ion diffusion length and solution resistance. Here we demonstrate a new device architecture with interpenetrating gyroid electrodes, providing a large ion-accessible surface area and gas diffusion channels. This design significantly reduces the interelectrode distance, lowering ion diffusion length and solution resistance. Simulations show faster ion diffusion and higher current density in the interpenetrating configuration compared with separate electrodes. This improved performance, especially at low temperatures and high current densities, highlights a promising strategy for enhancing AWS and other electrochemical systems limited by slow ion diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Shifts in Carbon Emissions Versus Sequestration From Hydropower Reservoirs in the Southeastern United States

Reservoirs are a significant source of carbon (C) to the atmosphere, but their emission rates vary in space and time. Here we compared C emissions via diffusive and ebullitive pathways at several stations in six large hydropower reservoirs in the southeastern US that were previously sampled in summer 2012. We found that carbon dioxide (CO 2 ) diffusion was the dominant flux pathway during 2012 and 2022, with only three exceptions where methane (CH 4 ) diffusion or CH 4 ebullition dominated. CH 4 diffusion rates were positively associated with water temperature. However, we found no clear predictors of CH 4 ebullition, which had extremely high variability, with rates ranging from 0 to 739 mg C m -2 day -1 . For CO 2 diffusion, the direction of the flux shifted between 2012 and 2022, where all but three stations across all reservoirs emitted CO 2 in summer 2012, but every station sequestered CO 2 in summer 2022. Here, indicators of greater algal production were associated with CO 2 sequestration, including surface chlorophyll- a concentration, surface dissolved oxygen saturation, and pH. Additional sampling campaigns outside the summer season highlighted the importance of seasonal phenology in primary production on the direction of CO 2 diffusive fluxes, which shifted to positive CO 2 fluxes by the end of August as productivity decreased. Our results demonstrate the importance of capturing CO 2 sequestration in field and modeling measurements and understanding the seasonal drivers of these estimates. Measuring C emissions from multiple pathways in reservoirs and understanding their spatiotemporal responses and variability are vital to reducing uncertainties in global upscaling efforts.

13 HYDRO ENERGY↗

DOE OSTI2025

Multiscale, mechanistic modeling of cesium transport in silicon carbide for TRISO fuel performance prediction

Understanding cesium (Cs) transport in TRistructural ISOtropic (TRISO) particle fuel is crucial for predicting fission product release in high-temperature reactors. However, current challenges include significant scatter in diffusivity data and unexplained temperature-dependent diffusion regimes in the silicon carbide layer. This study addresses these challenges by developing a multiscale, mechanistic Cs transport model integrating atomistic simulations and phase field modeling. Our model quantifies temperature and grain size effects on Cs diffusivity, attributing experimentally observed regimes to a transition from bulk-dominated diffusivity at high temperatures to grain boundary-dominated diffusivity at lower temperatures. The model, validated against diffusion measurements and advanced gas reactor (AGR)-1 and AGR-2 post-irradiation fission product release data, enhances the predictive capability of the BISON fuel performance code. This study advances our understanding of Cs release from TRISO particles and its dependence on temperature and silicon carbide grain size, with implications for the safety and efficiency of high-temperature nuclear reactors.

BISON↗

DOE OSTI2025

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗

DOE OSTI2024

Cosmic-Ray Feedback on Bistable Interstellar Medium Turbulence

Abstract While cosmic rays ( E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel Jana Howard1,3, Guang Yang3, Haiyan Zhao1*, Patrick Warren4, Tiankai Yao3, Steven Cavazos4 Elizabeth Sooby Wood4*, Ching-heng Shiau2 1 University of Idaho, Environmental Science, Idaho Falls, Idaho, USA 2 Boise State University, Microscopy and Characterization Suite, Idaho Falls, Idaho, USA 3 Idaho National Laboratory, Idaho Falls, Idaho, USA 4 University of Texas San Antonio, Department of Physics and Astronomy, San Antonio, Texas, USA *Corresponding author: haiyanz@uidaho.edu; elizabeth.soobywood@utsa.edu Tri-structural Isotropic (TRISO) particle fuel is the preferred choice for the newest and next-generation high-temperature nuclear reactors due to its robust construction [1]. The fuel design effectively contains fission products inside the particle at extremely high temperatures [2]. However, the effect of transition metal fission products on fuel performance and structural integrity remains largely unknown, creating uncertainties in predicting fuel performance and potential failure mechanisms [2]. For example, palladium (Pd) has a strong tendency to form intermetallic phases with uranium (U). These new phases tend to be hard and brittle which could lead to fracture of fuel kernel during irradiation [3,4]. Pd is not normally produced in high concentrations in normal fission reactions however, driving the production of Pd into the locality of a Uranium Oxycarbide (UCO) nuclear fuel provides an opportunity to closely observe the diffusion and chemical interactions. This study focuses on detailed transmission electron microscopy characterization of the intermetallic phases formed by the interaction between a UCO kernel and a Pd bar after annealing at 1100 °C for 100 hours. Figure 1 provides an overview of the UCO kernel in contact with the Pd bar. The UCO-Pd sample surface was examined using a Focused Ion Beam (FIB) Quanta 3D FEG FEI in VCD detector mode at 10kV, 83pA to reveal surface features. Several key surface features in the interaction region were observed including: (1) a small area with dendritic microstructure, (2) color contrast near the UCO kernel and Pd boundary, and (3) darker and lighter regions throughout the sample surface. To analyze diffusion behavior, FIB was used to prepare lamellae from five different areas of the UCO-Pd sample, as well as one for the as-received sample. These lamellae were examined using a Scanning Transmission Electron Microscope ThermoFisher Spectra 300 (STEM-Spectra). Energy Dispersive Spectroscopy results confirmed the diffusion between the UCO fuel and Pd bar. Figure 2 shows the identified phases of UC, UO2, and Pd dispersed throughout the interaction region. It was also found that the Pd concentration decreases as the radial distance from the UCO-Pd interface increases. The Pd concentration was 36.93% atomic fraction 11µm from the interaction region then decreased to 7.51% atomic fraction at 257 µm away from the interaction region. Figure 3 shows how concentrations of U and Pd vary across the interaction zone, highlighting the extent of diffusion. This study confirms that Pd interacts with surrounding material to form U-Pd phase and diffusion zones. These diffusion zones vary in composition depending on its radial distance from the point of contact with the Pd bar. These findings contribute to a deeper understanding of the fission product behavior in high temperature nuclear fuels, aiding in the prediction and mitigation of potential fuel degradation mechanisms. A B C Fig. 1. SEM micrographs of UCO fuel kernel in contact with solid Pd bar and the lift out locations. Figure 1A shows UCO fuel kernel in contact with solid Pd bar. Figure 1B shows a higher magnification image of the UCO fuel-Pd interaction zone. Figure 1C shows lift out sites for the lamellae. A B C Fig 2. EDS maps reveal the distribution of U, O, and Pd in location 3. Figure 2A shows the distribution of uranium. Figure 2B shows the distribution of oxygen. Figure 2C shows the distribution of palladium. A B C Fig 3. Palladium and uranium concentration across interaction zone in location 2. Figure 3A shows EDS map of lift out number 2. Figure 3B shows a zoomed in EDS map taken from the interaction zone. Figure 3C shows a line graph of uranium and palladium concentrations across the interaction zone. References: 1. B.E. Wells, N.R. Phillips, K.J. Geelhood. Pacific West Laboratory. (2021). TRISO Fuel: Properties and Failure Modes. https://www.nrc.gov/docs/ML2117/ML21175A152.pdf (Accessed January 16, 2025). 2. American Nuclear Society. TRISO Fuel Development Progresses in INL, ORNL. https://www.gen-4.org/gif/upload/docs/application/pdf/2014-03/nov13nn_fuel_reprint.pdf (Accessed January 16, 2025) 3. Clark R.A., M.A. Con

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel Jana Howard1,3, Guang Yang3, Haiyan Zhao1*, Patrick Warren4, Tiankai Yao3, Steven Cavazos4 Elizabeth Sooby Wood4*, Ching-heng Shiau2 1 University of Idaho, Environmental Science, Idaho Falls, Idaho, USA 2 Boise State University, Microscopy and Characterization Suite, Idaho Falls, Idaho, USA 3 Idaho National Laboratory, Idaho Falls, Idaho, USA 4 University of Texas San Antonio, Department of Physics and Astronomy, San Antonio, Texas, USA *Corresponding author: haiyanz@uidaho.edu; elizabeth.soobywood@utsa.edu Tri-structural Isotropic (TRISO) particle fuel is the preferred choice for the newest and next-generation high-temperature nuclear reactors due to its robust construction [1]. The fuel design effectively contains fission products inside the particle at extremely high temperatures [2]. However, the effect of transition metal fission products on fuel performance and structural integrity remains largely unknown, creating uncertainties in predicting fuel performance and potential failure mechanisms [2]. For example, palladium (Pd) has a strong tendency to form intermetallic phases with uranium (U). These new phases tend to be hard and brittle which could lead to fracture of fuel kernel during irradiation [3,4]. Pd is not normally produced in high concentrations in normal fission reactions however, driving the production of Pd into the locality of a Uranium Oxycarbide (UCO) nuclear fuel provides an opportunity to closely observe the diffusion and chemical interactions. This study focuses on detailed transmission electron microscopy characterization of the intermetallic phases formed by the interaction between a UCO kernel and a Pd bar after annealing at 1100 °C for 100 hours. Figure 1 provides an overview of the UCO kernel in contact with the Pd bar. The UCO-Pd sample surface was examined using a Focused Ion Beam (FIB) Quanta 3D FEG FEI in VCD detector mode at 10kV, 83pA to reveal surface features. Several key surface features in the interaction region were observed including: (1) a small area with dendritic microstructure, (2) color contrast near the UCO kernel and Pd boundary, and (3) darker and lighter regions throughout the sample surface. To analyze diffusion behavior, FIB was used to prepare lamellae from five different areas of the UCO-Pd sample, as well as one for the as-received sample. These lamellae were examined using a Scanning Transmission Electron Microscope ThermoFisher Spectra 300 (STEM-Spectra). Energy Dispersive Spectroscopy results confirmed the diffusion between the UCO fuel and Pd bar. Figure 2 shows the identified phases of UC, UO2, and Pd dispersed throughout the interaction region. It was also found that the Pd concentration decreases as the radial distance from the UCO-Pd interface increases. The Pd concentration was 36.93% atomic fraction 11µm from the interaction region then decreased to 7.51% atomic fraction at 257 µm away from the interaction region. Figure 3 shows how concentrations of U and Pd vary across the interaction zone, highlighting the extent of diffusion. This study confirms that Pd interacts with surrounding material to form U-Pd phase and diffusion zones. These diffusion zones vary in composition depending on its radial distance from the point of contact with the Pd bar. These findings contribute to a deeper understanding of the fission product behavior in high temperature nuclear fuels, aiding in the prediction and mitigation of potential fuel degradation mechanisms. A B C Fig. 1. SEM micrographs of UCO fuel kernel in contact with solid Pd bar and the lift out locations. Figure 1A shows UCO fuel kernel in contact with solid Pd bar. Figure 1B shows a higher magnification image of the UCO fuel-Pd interaction zone. Figure 1C shows lift out sites for the lamellae. A B C Fig 2. EDS maps reveal the distribution of U, O, and Pd in location 3. Figure 2A shows the distribution of uranium. Figure 2B shows the distribution of oxygen. Figure 2C shows the distribution of palladium. A B C Fig 3. Palladium and uranium concentration across interaction zone in location 2. Figure 3A shows EDS map of lift out number 2. Figure 3B shows a zoomed in EDS map taken from the interaction zone. Figure 3C shows a line graph of uranium and palladium concentrations across the interaction zone. References: 1. B.E. Wells, N.R. Phillips, K.J. Geelhood. Pacific West Laboratory. (2021). TRISO Fuel: Properties and Failure Modes. https://www.nrc.gov/docs/ML2117/ML21175A152.pdf (Accessed January 16, 2025). 2. American Nuclear Society. TRISO Fuel Development Progresses in INL, ORNL. https://www.gen-4.org/gif/upload/docs/application/pdf/2014-03/nov13nn_fuel_reprint.pdf (Accessed January 16, 2025) 3. Clark R.A., M.A. Con

36 - MATERIALS SCIENCE↗