SEARCH · Engineering Papers

Results for “SELF-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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Self-Diffusivity Measurement of Eutectic F 7 LiNaK with and without Additives Using Quasi-Elastic Neutron Scattering

The atomic scale relaxation dynamics of eutectic F 7 LiNaK (46.5 LiF–11.5 NaF–42 KF mol %, Li-7 enriched) were measured using quasi-elastic neutron scattering (QENS) over a temperature range of 500–750 °C. Here, the effect of adding 0.988 mol % cerium, 0.499 mol % cesium, and 1.21 mol % zirconium individually to the dynamics of F 7 LiNaK was also investigated. The relaxation process in both pure and doped F 7 LiNaK molten salts was fit with a stretched exponential function and the temperature dependence follows an Arrhenius behavior over a wavevector transfer range of 0.4 Å –1 < Q < 0.9 Å –1 . The measured activation energy for self-diffusion is E a = 0.77 ± 0.02 eV/atom for pure molten F 7 LiNaK. The QENS response with additives added to F 7 LiNaK was also fit with a stretched exponential and the associated Arrhenius behavior was characterized with activation energies of E a = 0.88 ± 0.01 eV/atom for zirconium (1.21 mol %), E a = 1.02 ± 0.02 eV/atom for cerium (0.988 mol %), and E a = 0.71 ± 0.03 eV/atom for cesium (0.499 mol %). The measured diffusivities are compared to those simulated with a neural network force field model by Lee et al. [Lee, S.-C. Comparative Studies of the Structural and Transport Properties of Molten Salt FLiNaK Using the Machine-Learned Neural Network and Reparametrized Classical Forcefields. J. Phys. Chem. B 2021, 125(37), 10562–10570].

FLiNaK↗

DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Nonadditive CO 2 Uptake of Type II Porous Liquids Based on Imine Cages

Type II porous liquids can potentially exploit the fluidity of liquids and sorption properties of porous sorbents, yet CO 2 uptake in porous liquids is still poorly understood. Molecular simulations and experiments are used to examine CO 2 uptake by a prototypical porous liquid composed of porous organic cages (CC13) in 2′-hydroxyacetophenone (2′-HAP). The simulations are in reasonable agreement with experimental measurements of CO 2 solubility and provide unambiguous information on the partitioning of CO 2 within microenvironments in the liquid. Analysis of CO 2 dynamics is performed using these simulations, including assessing the self-diffusivity of CO 2 in both the neat solvent and porous liquid. This offers insights into the kinetics of CO 2 uptake and transport in type II porous liquids based on imine cages. Experiments with type II porous liquids formed by dissolving CC13 in three different size-excluded solvents show nonadditive CO 2 absorption relative to predictions based on ideal volume additivity. This nonadditive absorption behavior is also observed in simulations. Nonadditive CO 2 uptake is also demonstrated in type II porous liquids based on another imine-based porous cage, CC19.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Interfacial Hydrogen-Bond Dynamics in Transition Metal Compounds

Understanding how water behaves when confined within atomic layers of active transition-metal carbides, nitrides, and carbonitrides is essential for uncovering the fundamental principles needed to engineer solid–liquid interfaces at the atomic scale. Yet, how lattice element chemistry and surface termination groups collectively regulate the structure and mobility of such interlayer water remains poorly understood. Here, we present a composition-controlled investigation of interlayer water dynamics in layered transition-metal nitride, carbide, and carbonitride systems using a systematic integration of quasi-elastic neutron scattering (QENS), ab initio molecular dynamics (AIMD) simulations, and density functional theory (DFT) calculations. QENS measurements show that nitride-rich systems host mobile, translationally diffusing water with thermally activated self-diffusion coefficients on the order of 10 –10 m 2 s –1 , whereas mixed C/N lattices confine water to localized, nontranslational motion that is insensitive to temperature. AIMD and DFT reveal that lattice C/N chemistry and surface functional group composition reshape the first hydration layer by modulating the surface electronic structure and termination-dependent hydrogen-bond networks, leading to pronounced differences in water ordering and thermal resilience. On the other hand, fully carbide systems exhibit intermediate behavior, highlighting that water mobility is not primarily controlled by the hydration level alone but by the coupling between lattice composition and surface chemistry. Overall, this study establishes how surface chemistry and lattice composition jointly control interfacial hydrogen bond dynamics, offering a mechanistic framework for designing transition-metal layered materials with tailored interfacial transport properties.

Hydration↗

DOE OSTI2025

Equilibrium Separation of Siloxanes in Metal–Organic Frameworks

We present an in silico assessment of metal–organic frameworks (MOFs) for the equilibrium separation of linear and cyclic siloxanes. Using a combination of configurational bias/continuous fractional component Monte Carlo (CB/CFCMC) simulations and Ideal Adsorbed Solution Theory (IAST), we investigated the adsorption of both equimolar and nonequimolar mixtures of linear and cyclic siloxanes in a selection of synthesizable MOFs with medium to large pore volumes. We showed that configurational entropy effects drive the preferential adsorption of linear siloxanes over cyclic siloxanes. Based on synthesizability metrics we identified ZIF-70 as a promising adsorbent for the separation of linear and cyclic siloxanes. Self-diffusivities of linear and cyclic siloxanes in ZIF-70 calculated from molecular dynamics simulations show that equilibrium can be reached on reasonable time scales. We explored vacuum-temperature swing adsorption (VTSA) as a potential method for the recovery of adsorbed linear siloxanes from ZIF-70, achieving a 75% reduction in adsorbed phase concentration. Additionally, we demonstrated that supercritical CO 2 offers an alternative desorption strategy by displacing adsorbed linear siloxanes in ZIF-70 at pressures above 200 bar driven by entropy effects.

Adsorption↗

DOE OSTI2025

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗

DOE OSTI2025

Determining hexavalent chromium transport properties in alkaline nuclear waste using nuclear magnetic resonance spectroscopy

This study focuses on the transport properties of hexavalent chromium, specifically the chromate anion, to improve predictive models and environmental remediation strategies for Cr(VI) migration. Using 53 Cr Nuclear Magnetic Resonance (NMR) spectroscopy, the research quantifies chromate in multicomponent electrolytes replicating nuclear waste conditions at the Hanford Site in Washington State. The consistency of the 53 Cr NMR signal integral with chromate concentration, despite varying matrix compositions, establishes it as a reliable concentration indicator. The transport properties of chromate in an alkaline solution were assessed using relaxation-based measurements via saturation recovery and Carr-Purcell-Meiboom-Gill experiments, determining spin-lattice and spin-spin relaxation times. These measurements, combined with the Bloembergen-Purcell-Pound equation, helped estimate the rotational correlation time and the 53 Cr self-diffusion coefficient using Stokes-Einstein-Debye and Stokes-Einstein equations. Direct measurements were obtained through pulsed field gradient stimulated echo 53 Cr NMR spectroscopy. Monte Carlo simulations further estimated uncertainty propagation. The results enhance comprehension of chromate transport and highlight prospects for identifying transport properties of NMR-active nuclei, traditionally considered unreachable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Cation selectivities in zwitterion grafted nanopores: effect of zwitterion architecture

Selective separation of monovalent cations is a critical challenge in applications such as water purification and lithium recovery from salt brines. Cross-linked zwitterionic amphiphilic copolymer (ZAC-X) membranes have gained attention for their exceptional anion permselectivity, attributed to self-assembled zwitterion-lined nanodomains that interact preferentially with anions according to their hydrated radii r hyd . However, these membranes show minimal selectivity among monovalent cations, despite significant differences in their hydration structures, motivating studies on the underlying mechanisms of cation 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 zwitterionic (ZI) ligands on cation diffusivities, partitioning, and permeabilities. To this end, we examined two contrasting ZI ligand architectures: Motif A (surface–cation–anion, S–ZI + –ZI − ) and Motif B (surface–anion–cation, S–ZI − –ZI + ). Our results show that in Motif A, the sulfonate groups of the ZI ligands are localized near the pore center radially, leading to strong electrostatic interactions with small bare cations (Mg 2+ and Li + ). This configuration results in high cation partitioning but low cation diffusion, maintaining solution-diffusion tradeoff typical of functionalized membranes. In contrast, Motif B show that sulfonate groups shift radially toward the mid-region of the pore. This shift, especially for small bare cations, introduces steric constraints that weaken their interactions with the sulfonate groups, thereby enhancing hydration and lowering partitioning, while still maintaining their low self-diffusivity. These findings establish zwitterion dipole orientation as a powerful design lever for tuning cation selectivity in membrane systems and offer molecular-level insights for engineering next-generation ion separation materials.

Morishita, Kazuya [Univ. of Texas, Austin, TX (Uni↗

DOE OSTI2026

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗

DOE OSTI2026

Effects of confinement and pressure on the structure and dynamics of carbon dioxide in silica slit pores

An understanding of carbon dioxide fluid properties within geological mesopores is important in applications ranging from carbon sequestration to shale oil recovery. Here, a molecular dynamics study is presented that aims to shed light on these systems by simulation of CO 2 fluid confined in β-cristobalite silica slit pores with different pore widths and pressure conditions. The weakly associating nature of carbon dioxide leads to little difference in structural and interfacial dynamical properties for different pore sizes and pressures. Rather, the behavior is found to be dominated by the entropic effects, namely, how the CO 2 organizes next to the silica surface. The CO 2 self-diffusion coefficient shows the strongest pore size dependence. It is strongly diminished in small pores and does not reach the bulk fluid value even in 6 nm pores. It also decreases with pressure, yielding an activation volume that increases with pore size. These results provide new insight into the behavior of a compressible fluid in nanoscale confinement.

Godahewa, Sahan M. [Univ. of Kansas, Lawrence, KS ↗

DOE OSTI2026

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗

DOE OSTI2025

Ab Initio Insights to Metal Passivation: Diffusion and Defect Formation in Amorphous Zirconia and Alumina

Metal passivation refers to the formation of protective oxide films on metals, which shield them from further corrosion and oxidation, playing a crucial role in maintaining their stability. The mesoscopic Point Defect Model has successfully predicted passivity as a steady state process where oxide growth from oxygen vacancies at the metal/film interface competes with oxide dissolution at the film/environment interface. In this work, informed by the Point Defect Model parameters, we use first-principles calculations to calculate defect formation and atomic diffusion in amorphous materials and correlate these descriptors with the behavior and growth of the oxide film. Focusing on amorphous zirconia and alumina, we demonstrate that defect formation energies exhibit significant variability in amorphous systems. In alumina, vacancies dominate, with cation and anion vacancies occurring at comparable concentrations. Diffusion calculations for stoichiometric amorphous alumina and zirconia, as well as oxygen-deficient zirconia, reveal faster diffusion in the oxygen-deficient case, highlighting the impact of defects on transport. Comparison of calculated self-diffusion coefficients for the dominant defect species with experimentally measured oxide thicknesses shows a clear correlation, suggesting that first-principles-derived diffusivity information can serve as a key descriptor for surface passivation film growth.

Farnell, Mackinzie S. [University of California, B↗

DOE OSTI2025

Microstructure-Based Degradation Modeling in Solid Oxide Cells

The presentation is intended for a graduate class in Georgia Southern University taught by Dr. Hayri Sezer on electrochemistry, SOFC and Li ion battery. It summaries NETL's recent works on microstructure-based modeling of degradation in solid oxide cells. It covers the degradation mechanism of Ni coarsening, Ni migration and the build-up of oxygen partial pressure inside the electrolyte. The phase-field models developed for Ni coarsening and migration are reviewed and the results are compared to experiments in literature to examine different diffusion mechanisms (Ni self-diffusion, Ni(OH)2 diffusion in pores and Ni(OH) surface diffusion) and driving forces (Ni(OH)x concentration and Ni-YSZ contact angle change). It also demonstrates the developed model for oxygen partial pressure and shows that the predicted oxygen partial pressure can explain the Ni oxidation in the hydrogen electrode and crack growth on the electrolyte-oxygen electrode interface under certain operating conditions.

nickel coarsening↗

DOE OSTI2025

Nuclear Magnetic Resonance Dynamics of LiTFSI–Pyrazole Eutectic Solvents

Deep Eutectic Solvents (DESs) have emerged as promising candidates to replace conventional organic solvents in various technological applications due to their low vapor pressure, non-flammability, and ease of preparation at low costs. In particular, Type IV DESs, which are composed of metal salts and hydrogen bond donors, are possible replacements for lithium-ion battery electrolytes. In this study, we investigate the molecular dynamics of solvents of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and pyrazole (PYR) at varying LiTFSI:PYR molar ratios (1:2, 1:3, 1:4, 1:5) using Nuclear Magnetic Resonance Dispersion (NMRD) and Pulsed Field Gradient (PFG) Nuclear Magnetic Resonance (NMR). PFG NMR reveals composition-dependent diffusion trends, while NMRD provides molecular-level insights into the longitudinal relaxation rate (R 1 = 1/T 1 ). Notably, the LiTFSI:PYR (1:2) sample shows distinct behavior across both techniques, exhibiting enhanced relaxation rates and lower self-diffusion for 1 H compared to the other nuclei ( 19 F and 7 Li), suggestive of stronger and more efficient Li + –pyrazole interactions, as confirmed by the modeling of the relaxation profiles. Our study advances understanding of ion dynamics in azole-based eutectic solvents, supporting their potential use in safer battery electrolytes.

FFC-NMR↗

DOE OSTI2015

Molecular dynamics insight to phase transition in n-alkanes with carbon nanofillers

The present work aims to investigate the phase transition, dispersion and diffusion behavior of nanocomposites of carbon nanotube (CNT) and straight chain alkanes. These materials are potential candidates for organic phase change materials(PCMs) and have attracted flurry of research recently. Accurate experimental evaluation of the mass, thermal and transport properties of such composites is both difficult as well as economically taxing. Additionally it is crucial to understand the factors that results in modification or enhancement of their characteristic at atomic or molecular level. Classical molecular dynamics approach has been extended to elucidate the same. Bulk atomistic models have been generated and subjected to rigorous multistage equilibration. To reaffirm the approach, both canonical and constant-temperature, constant- pressure ensembles were employed to simulate the models under consideration. Explicit determination of kinetic, potential, non-bond and total energy assisted in understanding the enhanced thermal and transport property of the nanocomposites from molecular point of view. Crucial parameters including mean square displacement and simulated self diffusion coefficient precisely define the balance of the thermodynamic and hydrodynamic interactions. Radial distribution function also reflected the density variation, strength and mobility of the nanocomposites. It is expected that CNT functionalization could improve the dispersion within n-alkane matrix. This would further ameliorate the mass and thermal properties of the composite. Additionally, the determined density was in good agreement with experimental data. Thus, molecular dynamics can be utilized as a high throughput technique for theoretical investigation of nanocomposites PCMs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Lattice vacancy migration barriers in Fe-Ni alloys, and an indication as to why Ni atoms diffuse slowly: A first-principles study

Lattice vacancy migration barriers in ferromagnetic Fe 𝑥 ⁢Ni 1−𝑥 alloys (0.4 ≤ 𝑥 ≤ 0.6) are accurately quantified within the framework of ab initio electronic structure calculations using the nudged elastic band (NEB) method. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal L⁢1 0 phase—which is under consideration as a material for a rare-earth-free gap magnet for advanced engineering applications—are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni-vacancy interchanges encounter significantly higher energetic barriers than do Fe-vacancy interchanges. We contend that this aspect is a key factor in determining the differences in mobility between Fe and Ni atoms in this ferromagnetic alloy. Moreover, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarized electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarization around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. These results give atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.

density functional theory↗