Engineering Papers⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Interplay of molecular dynamics and radiative decay of a TADF emitter in a glass-forming liquid

Herein we investigate the role of molecular dynamics in the luminescent properties of a prototypical thermally activated delayed fluorescence (TADF) emitter, NAI-DMAC, in solution using a combination of temperature dependent time-resolved photoluminescence and absorption spectroscopies. We use a glass forming liquid, 2-methylfuran, to introduce an abrupt change in the temperature dependent diffusion dynamics of the solvent and examine the influence this has on the emission intensity of NAI-DMAC molecules. Comparison of experiment with first principles molecular dynamics simulations reveals that the emission intensity of NAI-DMAC molecules follows the temperature-dependent self-diffusion dynamics of the solvent. A marked reduction of emission intensity is observed as the temperature decreases toward the glass transition because the rate at which NAI-DMAC molecules can access emissive molecular conformations is greatly reduced. Below the glass transition, the diffusion dynamics of the solvent changes more slowly with temperature, which causes the emission intensity to decrease more slowly as well. The combination of experiment and computation suggests a pathway by which TADF emitters may transiently access a distribution of conformational states and avoid the need for an average conformation that strikes a balance between lower singlet–triplet energy splittings versus higher emission probabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Molecular engineering of fluoroether electrolytes for lithium metal batteries

Fluoroether solvents are promising electrolyte candidates for high-energy-density lithium metal batteries, where high ionic conductivity and oxidative stability are important metrics for design of new systems. Recent experiments have shown that these performance metrics, particularly stability, can be tuned by changing the fraction of ether and fluorine content. However, little is known about how different molecular architectures influence the underlying ion transport mechanisms and conductivity. Here, we use all-atom molecular dynamics simulations to elucidate the ion transport and solvation characteristics of fluoroether chains of varying length, and having different ether segment and fluorine terminal group contents. The design rules that emerge from this effort are that solvent size determines lithium-ion transport kinetics, solvation structure, and solvation energy. In particular, the mechanism for lithium-ion transport is found to shift from ion hopping between solvation sites located in different fluoroether chains in short-chain solvents, to ion–solvent co-diffusion in long-chain solvents, indicating that an optimum exists for molecules of intermediate length, where hopping is possible but solvent diffusion is fast. Consistent with these findings, our experimental measurements reveal a non-monotonic behavior of the effects of solvent size on lithium-ion conductivity, with a maximum occurring for medium-length solvent chains. Finally, a key design principle for achieving high ionic conductivity is that a trade-off is required between relying on shorter fluoroether chains having high self-diffusivity, and relying on longer chains that increase the stability of local solvation shells.

25 ENERGY STORAGE↗
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 OSTI2023

The influence of counterion structure identity on conductivity, dynamical correlations, and ion transport mechanisms in polymerized ionic liquids

We used equilibrium and non-equilibrium atomistic simulations to probe the influence of anion chemistry on the true conductivity, dynamical correlations, and ion transport mechanisms in polymeric ionic liquids. An inverse correlation was found between anion self-diffusivities, ionic mobilities, and the anion size for spherical anions. While some larger asymmetric anions had higher diffusivities than smaller spherical anions, their diffusivities and mobilities did not exhibit a direct correlation to the anion volumes. Here, the conductivity and anion dynamical correlations also followed the same trends as displayed by the diffusivity and mobility of anions. All the systems we examined displayed positively correlated motion among anions, suggesting a contribution that enhances the conductivity beyond the ideal Nernst–Einstein value. Analysis of ion transport mechanisms demonstrated very similar hopping characteristics among the spherical anions despite differences in their sizes.

36 MATERIALS SCIENCE↗
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 OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hindered segmental dynamics in associative protein hydrogels studied by neutron spin-echo spectroscopy

Transient binding between associating macromolecules can cause qualitative changes to chain dynamics, including modes of conformational relaxation and diffusion, through tethering effects imparted by long-range connectivity. Here, in this work, the role of binding on short-time segmental dynamics in associative polymer gels is investigated by neutron spin-echo (NSE) measurements on a class of model artificial coiled-coil proteins with a systematically varied architecture, probing timescales of 0.1–130 ns, and length scales close to the molecular radius of gyration. The results illustrate effects of transient cross-linking on chain dynamics on different timescales, manifested in changes in segmental relaxation behavior with variations in strand length, chain concentration, and sticker distribution (endblock- vs midblock-functionalized). In all gels, a short-time cooperative diffusion mode is seen over all wave vectors, analogous to a semidilute solution, with no transitions seen at any known structural length scale. However, the diffusion coefficients are found to decrease with increasing junction density across all gels, with the strand length and number of stickers per chain in each gel appearing to play a relatively minor role. The slowing of cooperative diffusion with junction density contrasts with classical predictions of a greater restoring force for fluctuation dissipation due to the increased elasticity, suggesting additional effects of the coiled-coil junctions such as an enhancement in local viscosity that slows dynamics. Notably, the relaxation rates for all gels can be rescaled by the interjunction spacing inferred from small-angle neutron scattering, where they collapse onto a master curve suggestive of self-similar dynamics even in networks with different strand lengths and chain architectures. On long timescales (but shorter than the junction exchange time), a slowing of network relaxation is observed, resulting in a nondecaying plateau in the spin-echo amplitude attributed to a freezing of chain dynamics due to tethering. A characteristic length scale corresponding to the extent of dynamic fluctuations is estimated for each gel, which appears to be smaller than the interjunction spacing but similar to the correlation blob size of the overlapping strands. The results indicate an important role of transient binding on molecular-scale dynamics in associative polymer gels, even on timescales shorter than the junction exchange time, in addition to its effects on long-range self-diffusion previously observed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Carbon capture in polymer-based electrolytes

Nanoparticle organic hybrid materials (NOHMs) have been proposed as excellent electrolytes for combined CO 2 capture and electrochemical conversion due to their conductive nature and chemical tunability. However, CO 2 capture behavior and transport properties of these electrolytes after CO 2 capture have not yet been studied. Here, we use a variety of nuclear magnetic resonance (NMR) techniques to explore the carbon speciation and transport properties of branched polyethylenimine (PEI) and PEI-grafted silica nanoparticles (denoted as NOHM-I-PEI) after CO 2 capture. Quantitative 13 C NMR spectra collected at variable temperatures reveal that absorbed CO 2 exists as carbamates (RHNCOO – or RR'NCOO – ) and carbonate/bicarbonate (CO 3 2– /HCO 3 – ). The transport properties of PEI and NOHM-I-PEI studied using 1 H pulsed-field-gradient NMR, combined with molecular dynamics simulations, demonstrate that coulombic interactions between negatively and positively charged chains dominate in PEI, while the self-diffusion in NOHM-I-PEI is dominated by silica nanoparticles. These results provide strategies for selecting adsorbed forms of carbon for electrochemical reduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

Diffusional creep in UO 2 informed by lower length scale simulations

Using molecular dynamics, we predict information at the atomistic scale used to develop a mechanistic UO 2 creep model for use in higher length-scale fuel performance codes. The ultimate objective of the model is to not only to capture the creep rates of UO 2 but to determine the dominant mechanism in the diffusional regime, which is still debated in the literature. It is important to have a model to capture the correct mechanisms for creep in UO 2 as this can be used as the foundation when applying to other fuels, such as doped UO 2 , and when irradiation is accounted. In last years NEAMS milestone (FY22), we developed a prelimnary model, however there were issues, for example, excessively high values of uranium vacancy concentrations at the grain boundary. This year we have addressed the issues with the previous version of the model, added a new term that accounts for the nucleation of dislocations at stress raisers (e.g., triple junctions) within the microstructure and discussed where there was disagreement in the literature about the underpinning physics (uranium self-diffusion at the grain boundary).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Modeling Ni Coarsening and Migration in the Hydrogen Electrode of Solid oxide Cells under Operating Conditions

Ni coarsening and migration is the most important degradation in the hydrogen electrode of solid oxide cells. This presentation reviews our works on Ni coarsening and migration in NETL over the past years. We used simulation techniques including phase-field modeling and density function theory to investigate possible mechanisms of Ni coarsening and migration, including self-diffusion of Ni, diffusion of gaseous Ni(OH)2, surface diffusion of Ni(OH)x and Ni-YSZ wettability change, and examined the effect of operating conditions on these mechanisms. So far, none of the mechanisms can fully explain the experiments. Remaining questions and possible paths forward are summarized and discussed.

Lei, Yinkai↗
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

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Determination of the Diffusion Coefficients of Binary CH4 and C2H6 in a Supercritical CO2 Environment (500–2000 K and 100–1000 atm) by Molecular Dynamics Simulations

The self-diffusion coefficients of carbonaceous fuels in a supercritical CO2 environment provide transport information that can help us understand the Allam Cycle mechanism at a high pressure of 300 atm. The diffusion coefficients of pure CO2 and binary CO2/CH4 and CO2/C2H6 at high temperatures (500 K~2000 K) and high pressures (100 atm~1000 atm) are determined by molecular dynamics simulations in this study. Increasing the temperature leads to an increase in the diffusion coefficient, and increasing the pressure leads to a decrease in the diffusion coefficients for both methane and ethane. The diffusion coefficient of methane at 300 atm is approximately 0.012 cm2/s at 1000 K and 0.032 cm2/s at 1500 K. The diffusion coefficient of ethane at 300 atm is approximately 0.016 cm2/s at 1000 K and 0.045 cm2/s at 1500 K. The understanding of diffusion coefficients potentially leads to the reduction in fuel consumption and minimization of greenhouse gas emissions in the Allam Cycle.

Energy & Fuels↗