Results for “diffusion behavior” — Page 3

DOE OSTI2022

Investigation of Molecular Diffusion at Block Copolymer Thin Films Using Maximum Entropy Method-Based Fluorescence Correlation Spectroscopy and Single Molecule Tracking

Fluorescence correlation spectroscopy (FCS) has been widely used to investigate molecular diffusion behavior in various samples. The use of the maximum entropy method (MEM) for FCS data analysis provides a unique means to determine multiple distinct diffusion coefficients without a priori assumption of their number. Comparison of the MEM-based FCS method (MEM-FCS) with another method will reveal its utility and advantage as an analytical tool to investigate diffusion dynamics. Herein, we measured diffusion of fluorescent probes doped into nanostructured thin films using MEM-FCS, and validated the results with single molecule tracking (SMT) data. The efficacy of the MEM code employed was first demonstrated by analyzing simulated FCS data for systems incorporating one and two diffusion modes with broadly distributed diffusion coefficients. The MEM analysis accurately afforded the number of distinct diffusion modes and their mean diffusion coefficients. These results contrasted with those obtained by fitting the simulated data to conventional two-component and anomalous diffusion models, which yielded inaccurate estimates of the diffusion coefficients. Subsequently, the MEM analysis was applied to FCS data acquired from hydrophilic dye molecules incorporated into microphase-separated polystyrene-block-poly(ethylene oxide) (PS-b-PEO) thin films characterized under a water-saturated nitrogen atmosphere. The MEM analysis revealed distinct fast and slow diffusion components attributable to molecules diffusing on the film surface and inside the film, respectively. SMT studies of the same materials yielded trajectories for mobile molecules that appear to follow the curved PEO microdomains. Diffusion coefficients obtained from the SMT data were consistent with those obtained for the slow diffusion component detected by MEM-FCS. Furthermore, these results highlight the utility of MEM-FCS and SMT for gaining complementary information on molecular diffusion processes in heterogeneous material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Effective field theory of fluctuating wall in open systems: from a kink in Josephson junction to general domain wall

We investigate macroscopic behaviors of fluctuating domain walls in nonequilibrium open systems with the help of the effective field theory based on symmetry. Since the domain wall in open systems breaks the translational symmetry, there appears a gapless excitation identified as the Nambu-Goldstone (NG) mode, which shows the non-propagating diffusive behavior in contrast to those in closed systems. After demonstrating the presence of the diffusive NG mode in the (2+1)-dimensional dissipative Josephson junction, we provide a symmetry-based general analysis for open systems breaking the one-dimensional translational symmetry. A general effective Lagrangian is constructed based on the Schwinger-Keldysh formalism, which supports the presence of the gapless diffusion mode in the fluctuation spectrum in the thin wall regime. Besides, we also identify a term peculiar to the open system, which possibly leads to the instability in the thick-wall regime or the nonlinear Kardar-Parisi-Zhang coupling in the thin-wall regime although it is absent in the Josephson junction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

High-temperature and high-pressure NMR investigations of low viscous fluids confined in mesoporous systems

Abstract In this contribution, the relaxation and diffusional behaviors of low viscous fluids, water and methanol confined into mesoporous silica and controlled size pore glass were investigated. The engineered porous systems are relevant to geologically important subsurface energy materials. The engineered porous proxies were characterized by Brunauer–Emmett–Teller (BET) surface analyzer, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (EM) to determine surface area, pore-wall protonation and morphology of these materials, respectively. The confined behavior of the low viscous fluids was studied by varying pore diameter, fluid-to-solid ratio, temperature, and pressure, and then compared to bulk liquid state. Both relaxation and diffusion behaviors for the confined fluids showed increasing deviation from pure bulk fluids as the fluid-to-solid ratio was decreased, and surface-to-volume ratio ( S / V ) was varied. Variable pressure deuteron NMR relaxation of confined D 2 O and confined methanol, deuterated at the hydroxyl or methyl positions, were performed to exploit the sensitivity of the deuteron quadrupole moment to molecular rotation. The methanol results demonstrated greater pressure dependence than those for water only in bulk. The deviations from bulk liquid behavior arise from different reasons such as confinement and the interactions between confined fluid and the nano-pore wall. The results of the present report give insight into the behavior of low viscosity fluid in nano-confined geometries under different state conditions.

Chemistry↗

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↗

DOE OSTI2020

Deuterium diffusion in γ-LiAlO2 pellets irradiated with He+ and D2+ ions

Not Available

Jiang, Weilin↗

DOE OSTI2024

Underlying mechanism of structural transformation between GaSb and GaAs response to intense electronic excitation

Ion irradiation of semiconductors has emerged as a promising approach for fabricating self-organized nanosystems with high atomic precision, despite often being accompanied by undesirable phenomena. Exploring the mechanisms underlying structural transformations is crucial for assessing nanostructure array types under complex irradiation environments. By quantitatively calculating the thermodynamically driven processes and analyzing the impact of intrinsic structural parameters, distinct structural transformations in response to intense electronic excitation are systematically investigated in gallium antimonide (GaSb) and gallium arsenide (GaAs) systems. In high-energy regimes, the nanofibers layer of GaSb exhibits intriguing structural discrepancy, characterized by partial nanofibers with coherent boundaries, interspersed nanopores accompanied by antisite defects and Ga precipitates, distinguishing to a series of discontinuous latent tracks that emerged within cylindrical trajectories in GaAs. Furthermore, significant diffusion behaviors of the nanohillocks are discovered in GaAs, with higher average roughness than GaSb, driven by the gradient stress distribution influenced by the free-surface effects. The deposition energy for melting phase formation, Gibbs free energy, and Ga diffusion coefficients contribute to the distinctive structural features, evidencing relatively stable morphological configurations and higher irradiation resistance in GaAs. Consequently, special optoelectronic properties associated with structural discrepancies facilitate the design and optimization of material functionalities by irradiation technologies.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effect of lithium incorporation on tweaking the electrocatalytic behavior of tantalum-based oxides

To study the effect of lithium incorporation on the electrochemical properties and catalytic behavior of tantalum-based oxides, a series of mesoporous lithium-contained tantalum oxides was developed adopting an inverse micelle-based sol gel synthesis method. Four various phases of amorphous Ta 2 O 5 , polycrystalline LiTaO 3 , Li 5.04 Ta 6.16 O 17.92 , and Li 3 TaO 4 were determined for samples synthesized at Li/Ta ratios of 0.5 (LTO-0.5), 1 (LTO-1), 2 (LTO-2), and 3 (LTO-3) respectively. The surface analysis revealed that LTO-2 sample is featured with lower amount of lithium atoms at the surface rather than the bulk of the sample. The electrochemical analysis of LTO samples showed that LTO-2 possesses lower charge transfer resistance extracted from impedance spectroscopy that consequently leads to a better electrocatalytic performance. The LTO-2 electrode also showed a low open circuit potential (OCP) combined with a relatively higher electroactive surface area, thus capable of running an electrooxidation reaction at a potential lower than all other LTO samples. The more negative OCP value of LTO-2 electrode in AA solution states a more facile electron transfer between the analyte and electrode compared to other synthesized samples. Mott–Schottky plots of LTO electrodes are characterized by positive slopes which indicate n-type semiconductor behavior. The charge carrier density (5.5 × 10 19 cm -3 ) in LTO-2 sample is the lowest among all other samples. The charge carrier density is dependent on the amount of applied lithium in the synthesis confirming various lithium diffusion behaviors in LTO samples. The lower carrier concentration in LTO-2 justifies its higher activity in electrooxidation providing less populated lower unoccupied molecular orbital (LUMO) for the accommodation of electrons extracted from the analyte higher occupied molecular orbital (HOMO). Therefore, LTO-2 showed an improved electrocatalytic performance due to lower charge transfer resistance (higher conductivity) offered by lower amount of lithium atom at its surface and higher oxygen vacancies in the bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗

DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗

DOE OSTI2023

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗

DOE OSTI2022

Influence of Thermal Treatment on Preceramic Polymer Grafted Nanoparticle Network Formation: Implications for Thermal Protection Systems and Aerospace Propulsion Components

Preceramic polymer grafted nanoparticles (PCP GNPs) offer potential advantages in the production of polymer-derived ceramic composites compared to neat preceramic polymers (PCPs), such as reduced thermal shrinkage and increased char yield. In comparison to traditional composites (particles + free polymer), PCPs avoid issues of compatibility, agglomeration, and phase separation during processing and provide routes to control nanoparticle arrangement. Prior literature has established that the conversion efficiency of PCPs and PCP GNPs to ceramics and ceramic composites is linked to the thermal pretreatment of these precursor polymers (e.g., curing). In this paper, we investigate the material transformations that occur when exposing PCP GNPs to a low-temperature thermal treatment (≤250 °C) prior to pyrolysis. The PCP GNPs explored in this study possessed silica nanoparticle cores but were distinct in their PCP coronas, having poly(1,1-dimethylpropylsilane) (allyl-GNP)- or poly(1,1-dimethylbenzylethylsilane) (styryl-GNP)-grafted polymers. After undergoing low-temperature thermal treatment, these PCP GNPs exhibited increased char yield at 800 °C; however, unlike commercial PCPs (e.g., allylhydridopolycarbosilane (SMP-10)), obvious cross-linkable sites are not present in the PCP corona structure. Differential scanning calorimetry, oscillatory rheology, and X-ray photon correlation spectroscopy were utilized to elucidate the thermally induced material changes in the allyl- and styryl-GNPs. Thermally induced changes to the structure and behavior of the PCP GNPs were determined to be linked to the chemical structure of the grafted polycarbosilane chains. Styryl-GNPs experienced a curing event between 180 and 250 °C, effecting the formation of a network, which contributed to char yield improvement (at 800 °C). The nanoscale dynamics of this material also shows a transition from diffuse behavior to ballistic behavior because of the cross-linking event. Conversely, the allyl-GNPs did not exhibit a curing event under heat treatment but did show an increase in thermally induced physical cross-linking. The allyl-GNP physical network is stronger after the first cycle of thermal treatment, and it is considered that this contributed to the improvement in char yield post thermal treatment. The advancements in the understanding of the cross-linking behavior of these hybrid materials are expected to advance the application of these materials to turbine engine, advanced friction, and heat-shielding components.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Reversible to irreversible transitions for ac driven skyrmions on periodic substrates

Abstract Using atomistic simulations, we investigate the dynamical behavior of magnetic skyrmions in dimer and trimer molecular crystal arrangements, as well as bipartite lattices at 3/2 and 5/2 fillings, under ac driving over a square array of anisotropy defects. For low ac amplitudes, at all fillings reversible motion appears in which the skyrmions return to their original positions at the end of each ac drive cycle and the diffusion is zero. We also identify two distinct irreversible regimes. The first is a translating regime in which the skyrmions form channels of flow in opposing directions and translate by one substrate lattice constant per ac drive cycle. The translating state appears in the dimer and trimer arrangements, and produces pronounced peaks in the diffusivity in the direction perpendicular to the external drive. For larger ac amplitudes, we find chaotic irreversible motion in which the skyrmions can randomly exchange places with each other over time, producing long-time diffusive behavior both parallel and perpendicular to the ac driving direction.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Li-Ion Diffusion Correlations in LiAlGeO 4 : Quasielastic Neutron Scattering and Ab Initio Simulation

Here, we investigated the impact of Li stoichiometry and host flexibility on Li + diffusion processes in LiAlGeO 4 at the microscopic level using quasielastic neutron scattering (QENS) and ab initio molecular dynamics (AIMD) simulations. Using sufficiently long AIMD trajectories, we could simulate the observed QENS signal and identify the localized dynamics of Li in crystalline LiAlGeO 4 . Such information is vital to identify the bottleneck of diffusion processes and design materials for battery application. Our AIMD simulations in LiAlGeO 4 reveal that the Li + conductivity can be significantly improved by manipulating the Li stoichiometry and/or host flexibility via amorphization. We determined that excess Li stoichiometry enhances the Coulomb repulsion of neighboring Li sites and softens the host structure to enable faster Li + diffusion along the hexagonal c-axis. In the amorphous structure, random orientations of AlO 4 and GeO 4 polyhedral units create a wide distribution of intersite distances and significantly soften the host structure, greatly enhancing the Li + diffusion. The simulations are used to understand the nature of diffusion, especially the role of the host structure, the possible hopping pathways, and the diffusion behavior in various structures of LiAlGeO 4 .

localized dynamics↗

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 OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗